Gradients in cellular metabolism regulate development of tonotopy along the cochlea. In vertebrates with elongated auditory organs, mechanosensory hair cells (HCs) are organised such that complex sounds are broken down into their component frequencies along the basal-to-apical long (tonotopic) axis. To generate the frequency-specific characteristics required at the appropriate positions, nascent HCs must determine their tonotopic properties during development. This relies on complex signalling between the developing HC and its local niche along the axis within the auditory organ. Here we apply NAD(P)H fluorescence lifetime imaging (FLIM) and live imaging of mitochondria to reveal metabolic gradients along the tonotopic axis of the developing cochlea. We further show that re-shaping these gradients during development, by inhibiting cytosolic glucose metabolism, alters normal Bmp7 and Chordin like-1 signalling leading to flattening of tonotopic patterning along the axis. Our work supports a causal link between morphogen signalling and metabolic reprogramming in specifying tonotopic morphologies in auditory HCs.
Precartilage condensation during limb skeletogenesis occurs by tissue phase separation controlled by a bistable cell-state switch with suppressed oscillatory dynamics The tetrapod limb skeleton is initiated in unpatterned limb bud mesenchyme by the formation of precartilage condensations. Here, based on time-lapse videographic analysis of a forming condensation in a high-density culture of chicken limb bud mesenchyme, we observe a phase transition to a more fluidized state for cells within spatial compacted foci (protocondensations that will progress to condensations), as reflected in their spatial confinement, cell-substratum interaction and speed of motion. Previous work showed that galectin-8 and galectin-1A, two proteins of the galactoside-binding galectin family, are the earliest determinants of this process in the chicken limb bud, and that their interactions in forming skeletogenic patterns of condensations can be interpreted mathematically through a reaction-diffusion-adhesion framework. Based on this framework, we use an ordinary differential equation-based approach to analyze the core switching modality of the galectin reaction network and characterize the states of the network independent of the diffusive and adhesive arms of the patterning mechanism. We identify two steady states where the concentrations of both galectins are respectively, negligible, and very high. An explicit Lyapunov function shows that there are no periodic solutions. For sigmoidal galectin production terms, the model exhibits a bistable switch that arises from a monostable state via saddle-node bifurcation. Our model therefore predicts that the galectin network exists in low and high expression states separated in space or time without any intermediate states. This provides a causal basis for the observed outside vs. inside transition observed in the in vitro video data. We performed a quantitative analysis of the distribution of galectin-1A in cultures of condensing chick limb mesenchymal cells and found that the interior of the protocondensations had concentrations of this protein (compared to the immediate exterior) over and above that expected from its higher cell density, consistent with the models predictions. The galectin-based patterning network is thus suggested, on theoretical grounds, to incorporate a core switch independent of any spatial or temporal dynamics, that drives the chondrogenic cell state transition in limb skeletogenesis.
Fine-tuning the Spike: Role of the nature and topology of the glycan shield in the structure and dynamics of SARS-CoV-2 S The dense glycan shield is an essential feature of the SARS-CoV-2 spike (S) architecture, key to immune evasion and to the activation of the prefusion conformation. Recent studies indicate that the occupancy and structures of the SARS-CoV-2 S glycans depend not only on the nature of the host cell, but also on the structural stability of the trimer; a point that raises important questions about the relative competence of different glycoforms. Moreover, the functional role of the glycan shield in the SARS-CoV-2 pathogenesis suggests that the evolution of the sites of glycosylation is potentially intertwined with the evolution of the protein sequence to affect optimal activity. Our results from multi-microsecond molecular dynamics simulations indicate that the type of glycosylation at N234, N165 and N343 greatly affects the stability of the receptor binding domain (RBD) open conformation, and thus its exposure and accessibility. Furthermore, our results suggest that the loss of glycosylation at N370, a newly acquired modification in the SARS-CoV-2 S glycan shields topology, may have contributed to increase the SARS-CoV-2 infectivity as we find that N-glycosylation at N370 stabilizes the closed RBD conformation by binding a specific cleft on the RBD surface. We discuss how the absence of the N370 glycan in the SARS-CoV-2 S frees the RBD glycan binding cleft, which becomes available to bind cell-surface glycans, potentially increases host cell surface localization. Diverse stem-chondrichthyan oral structures and evidence for an independently acquired acanthodid dentition The teeth of sharks famously form a series of transversely-organised files with conveyor-belt replacement that are borne directly on the jaw cartilages, in contrast to the dermal plate-borne dentition of bony fishes that undergoes site-specific replacement. A major obstacle in understanding how this system evolved is the poorly understood relationships of the earliest chondrichthyans and the profusion of morphologically and terminologically diverse bones, cartilages, splints and whorls that they possess. Here we use tomographic methods to investigate mandibular structures in several early branching  acanthodian-grade stem-chondrichthyans. We show that the dentigerous jaw bones of disparate genera of ischnacanthids are united by a common construction, being growing bones with non-shedding dentition. Mandibular splints, which support the ventro-lateral edge of the Meckels cartilage in some taxa, are formed from dermal bone and may be an acanthodid synapomorphy. We demonstrate that the teeth of Acanthodopsis are borne directly on the mandibular cartilage and that this taxon is deeply nested within an edentulous radiation, representing an unexpected independent origin of teeth. Many or even all of the range of unusual oral structures may be apomorphic, but they should nonetheless be considered when building hypotheses of tooth and jaw evolution, both in chondrichthyans and more broadly. A new method for the sampling and preservation of placental specimens in low-resource settings for the identification of P. falciparum and analysis of nucleic acids Collection, preservation, and shipment of histological specimens in low-resource settings is challenging. We present a novel method that achieved excellent preservation of placental specimens from rural Mali by using formalin fixation, ethanol dehydration, and long-term storage in a solar-powered freezer. Sample preservation success was 92%, permitting evaluation of current and past malaria infection, anemia, placental maturity, and inflammation. Using RNAscope(R) hybridization we were able to visualize cell-specific gene expression patterns in the formalin-fixed paraffin-embedded (FFPE) specimens. Additionally, our method entailed mirrored sampling from the two cut faces of a cotyledon, one for the FFPE workflows and the other for storage in RNAlater and RNA-seq. Juvenile Gouldian Finches (Erythrura gouldiae) form sibling sub-groups during social integration The formation of social relationships in complex groups is critical in shaping patterns of social organization and behavioral development. In many birds, young individuals remain dependent on their parents for extended periods but must abruptly transition to navigating interactions in the wider group after independence. While lack of social relationships during this period is detrimental in the development of later social skills, little is known about the social relationships juveniles form after independence in many bird species. In this study we describe patterns of social interactions in juvenile Gouldian Finches (Erythrura gouldiae) during transitions from family groups into flocks of unfamiliar individuals. A group of 20 juveniles from 4 families were introduced into two flocks. After introductions juveniles showed a gradient of approach rates with most approaches directed towards siblings, followed by juvenile peers, adult females, and lastly adult males. Significant preferences for siblings resulted in the emergence of sibling subgroups within the larger social network. This active self-assortment of siblings suggests that sibling sub-groups are an important bridge linking social connections within the family and the wider group. Such findings suggests that sibling relationships have a critical role in the socialization after independence, as well as structuring the social organization of Gouldian finch flocks. Sablefish (Anoplopoma fimbria Pallas, 1814) plasma biochemistry and hematology reference intervals including blood cell morphology Plasma biochemistry and hematology reference intervals are integral health assessment tools in all medical fields, including aquatic animal health. As sablefish (Anoplopoma fimbria) are becoming aquaculturally and economically more important, this manuscript provides essential reference intervals (RI) for their plasma biochemistry and hematology along with reference photomicrographs of blood cells in healthy, fasted sablefish. Blood cell morphology can differ between fish species. In addition, blood cell counts and blood chemistry can vary between fish species, demographics, water conditions, seasons, diets, and culture systems, which precludes the use of RIs from other fish species. For this study, blood was collected for plasma biochemistry and hematology analysis between June 20 and July 18, 2019, from healthy, yearling sablefish, hatched and reared in captivity on a commercial diet. Overnight fast of 16-18 hours did not sufficiently reduce lipids in the blood, which led to visible lipemia and frequent rupture of blood cells during analysis. Therefore, sablefish should be fasted for 24 to 36 hours before blood is collected to reduce hematology artifacts or possible reagent interference in plasma biochemistry analysis. Lymphocytes were the most dominant leukocytes (98%), while eosinophils were rare, and basophils were not detected in sablefish. Neutrophils were very large cells with Dohle bodies. In mammals and avian species, Dohle bodies are usually signs of toxic change from inflammation, but no such association was found in these fish. In conclusion, lipemia can interfere with sablefish blood analysis, and available removal methods should be evaluated as fasting for up to 36 h might not always be feasible. Also, more studies are required to establish RI for different developmental stages and rearing conditions. Diverse classes of constraints enable broader applicability of a linear programming-based dynamic metabolic modeling framework Current metabolic modeling tools suffer from a variety of limitations, from scalability to simplifying assumptions, that preclude their use in many applications. We recently created a modeling framework, LK-DFBA, that addresses a key gap: capturing metabolite dynamics and regulation while retaining a potentially scalable linear programming structure. Key to this frameworks success are the linear kinetics and regulatory constraints imposed on the system. However, while the linearity of these constraints reduces computational complexity, it may not accurately capture the behavior of many biochemical systems. Here, we developed three new classes of LK-DFBA constraints to better model interactions between metabolites and the reactions they regulate. We tested these new approaches on several synthetic and biological systems, and also performed the first-ever comparison of LK-DFBA predictions to experimental data. We found that no single constraint approach was optimal across all systems examined, and systems with the same topological structure but different parameters were often best modeled by different types of constraints. However, we did find that the optimal constraint approach was generally robust to local perturbations of the system, indicating that just a single wild-type dataset could allow identification of the ideal constraint for a given system. These results suggest that the availability of multiple constraint approaches will allow LK-DFBA to model a wider range of metabolic systems. Comprehensive-immunopeptidomics analysis reveals presentation of potential neoantigens carrying cancer driver mutations Understanding the properties of human leukocyte antigen (HLA) peptides (HLAp) is essential for cancer precision medicine, while the direct identification of HLAp from a tiny piece of clinical tissues by mass spectrometry (MS) is still confronted with technical challenges. Here, to overcome that hindrances, high-field asymmetric waveform ion mobility spectrometry (FAIMS) is introduced to conduct differential ion mobility (DIM)-MS by seamless gas-phase fractionation optimal for scarce samples. The global-immunopeptidomics analysis enables in-depth identification of 44,785 HLAp, including 2 neoantigens with KRAS-G12V and CPPED1-R228Q from stage IV colorectal cancer (CRC) tissues obtained from 17 subjects. Comparison of tissue-based personal immunopeptidome disclosed cancer-specific processing of HLAp. Since the direct identification of neoantigens from CRC tissues suggests more potential neoantigens have yet to been identified, the targeted-immunopeptidomics by screening-oriented parallel reaction monitoring (PRM) has been established and identified two more neoantigens with oncogenic KRAS. The targeted-immunopeptidomics analysis clarifies the presentation of aimed neoantigens, including neoantigens that are unpredictable by computational algorithms. The comprehensive-immunopeptidomics approach combining the FAIMS-assisted DIM-MS and the screening-oriented PRM contributes the direct and effective identification of HLAp and emerging neoantigens from clinical tissue samples, as an advanced strategy to identify the authentic and recurrent neoantigens for cancer immunotherapy. Common coupling map advances GPCR-G protein selectivity Two-thirds of human hormones and one-third of clinical drugs act on membrane receptors that couple to G proteins to achieve appropriate functional responses. While G protein transducers from literature are annotated in the Guide to Pharmacology database, two recent large-scale datasets now expand the receptor-G protein  couplome. However, these three datasets differ in scope and reported G protein couplings giving different coverage and conclusions on GPCR-G protein signaling. Here, we report a common coupling map uncovering novel couplings supported by both large-scale studies, the selectivity/promiscuity of GPCRs and G proteins, and how the co-coupling and co-expression of G proteins compare to the families from phylogenetic relationships. The coupling map and insights on GPCR-G protein selectivity will catalyze advances in receptor research and cellular signaling towards the exploitation of G protein signaling bias in design of safer drugs. Universal protein misfolding intermediates can bypass the proteostasis network and remain soluble and less-functional Misfolded protein conformations with decreased functionality can bypass the proteostasis machinery and remain soluble in vivo. This is an unexpected phenomenon as several cellular quality control mechanisms have evolved to rid cells of misfolded proteins. Three questions, then, are: how is it structurally possible for long-lived, soluble, misfolded proteins to bypass the proteostasis machinery and processes? How widespread are these soluble, misfolded states across the proteome? And how long do they persist for? Here, we address these questions using coarse-grain molecular dynamics simulations of the synthesis, termination, and post-translational dynamics of a representative set of cytosolic E. coli proteins. We predict that half of all proteins exhibit subpopulations of misfolded conformations that are likely to bypass molecular chaperones, avoid aggregation, and not be rapidly degraded. These misfolded states may persist for months or longer for some proteins. Structurally characterizing these misfolded states, we observe they have a large amount of native structure, but also contain localized misfolded regions from non-native changes in entanglement, in which a protein segment threads through a loop formed by another portion of the protein that is not found in the native state. The surface properties of these misfolded states are native like, suggesting they may bypass the proteostasis machinery and its regulatory processes to remain soluble, while their entanglements make these states long-lived kinetic traps, as disentanglement requires unfolding of already folded portions of the protein. In terms of function, we predict that one-third of proteins have subpopulations that misfold into less-functional states that have structurally perturbed functional sites yet remain soluble. Data from limited-proteolysis mass spectrometry experiments, which interrogate the misfolded conformations populated by proteins upon unfolding and refolding, are consistent with the structural changes seen in the entangled states of glycerol-3-phosphate dehydrogenase upon misfolding. These results provide an explanation for how proteins can misfold into soluble conformations with reduced functionality that can bypass cellular quality controls, and indicate, unexpectedly, this may be a wide-spread phenomenon in proteomes. Such entanglements are observed in many native structures, suggesting the non-native entanglements we observe are plausible. More broadly, these near-native entangled structures suggest a hypothesis for how synonymous mutations can modulate downstream protein structure and function, with these mutations partitioning nascent proteins between these kinetically trapped states. Naming Human Diseases: Ethical Principles of Curating Exclusive Substitute for Inopportune Nosology BackgroundIn the medical sphere, understanding naming conventions strengthen the integrity of naming human diseases remains nominal rather than substantial yet. Since the current nosology-based standard for human diseases could not offer a one-size-fits-all corrective mechanism, many idiomatic but flawed names frequently appear in scientific literature and news outlets at the cost of sociocultural impacts.

ObjectiveWe attempt to examine the ethical oversights of current naming practices and propose heuristic rationales and approaches to determine a pithy name instead of an inopportune nosology.

MethodsFirst, we examined the compiled global online news volumes and emotional tones on some inopportune nosology like German measles, Middle Eastern Respiratory Syndrome, Spanish flu, Hong Kong flu, and Huntingtons disease in the wake of COVID-19. Second, we prototypically scrutinize the lexical dynamics and pathological differentials of German measles and common synonyms by leveraging the capacity of the Google Books Ngram Corpus. Third, we demonstrated the empirical approaches to curate an exclusive substitute for an anachronistic nosology German measles based on deep learning models and post-hoc explanations.

ResultsThe infodemiological study shows that the public informed the offensive names with extremely negative tones in textual and visual narratives. The findings of the historiographical study indicate that many synonyms of German measles did not survive, while German measles became an anachronistic usage, and rubella has taken the dominant place since 1994. The PubMedBERT model could identify rubella as a potential substitution for German measles with the highest semantic similarity. The results of the semantic drift experiments further indicate that rubella tends to survive during the ebb and flow of semantic drift.

ConclusionsOur findings indicate that the nosological evolution of anachronistic names could result in sociocultural impacts without a corrective mechanism. To mitigate such impacts, we introduce some ethical principles for formulating an improved naming scheme. Based on deep learning models and post-hoc explanations, our illustrated experiments could provide hallmark references to the remedial mechanism of naming practices and pertinent credit allocations. Positive Selection and Heat-Response Transcriptomes Reveal Adaptive Features of the Brassicaceae Desert Model, Anastatica hierochuntica O_LIPlant adaptation to a desert environment and its endemic heat stress is poorly understood at the molecular level. The naturally heat-tolerant Brassicaceae species Anastatica hierochuntica is an ideal extremophyte model to identify genetic adaptations that have evolved to allow plants to tolerate heat stress and thrive in deserts.
C_LIO_LIWe generated an A. hierochuntica reference transcriptome and pinpointed extremophyte adaptations by comparing Arabidopsis thaliana and A. hierochuntica transcriptome responses to heat and identifying positively selected genes in A. hierochuntica.
C_LIO_LIThe two species exhibit similar transcriptome adjustment in response to heat and the A. hierochuntica transcriptome does not exist in a constitutive heat "stress-ready" state. Furthermore, the A. hierochuntica global transcriptome as well as heat-responsive orthologs, display a lower basal and higher heat-induced expression than in A. thaliana. Genes positively selected in multiple extremophytes are associated with stomatal opening, nutrient acquisition, and UV-B induced DNA repair while those unique to A. hierochuntica are consistent with its photoperiod-insensitive, early-flowering phenotype.
C_LIO_LIWe suggest that evolution of a flexible transcriptome confers the ability to quickly react to extreme diurnal temperature fluctuations characteristic of a desert environment while positive selection of genes involved in stress tolerance and early flowering could facilitate an opportunistic desert lifestyle.
C_LI Pathways of savannization in a mesic African savanna-forest mosaic following an extreme fire O_LIFires in savannas limit tree cover, thereby promoting flammable grass accumulation and fuelling further frequent fires. Meanwhile, forests and thickets form dense canopies that reduce C4-grass fuel loads and creating a humid microclimate, thereby excluding fires under typical climatic conditions.
C_LIO_LIHowever, extreme fires occasionally burn into these closed-canopy systems. Although these rare fires cause substantial tree mortality and can make repeat fires more likely, the long-term consequences of an extreme fire for closed canopy vegetation structure and potential to convert to savanna (hereafter "savannization") remain largely unknown.
C_LIO_LIHere, we analysed whether an extreme fire could, alone, alter species composition, vegetation structure, and fire regimes of closed-canopy ecosystems in an intact savanna-forest-thicket mosaic, or whether successive fires after an initial extreme fire were necessary to trigger a biome transition between from forest to savanna.
C_LIO_LIWe found that forests that only burned once recovered, whereas those that burned again following an initial extreme fire transitioned from closed-canopy forests towards open, grassy savannas.
C_LIO_LIWhile thickets had less tree mortality in fires than forests, repeat fires nonetheless precipitated a transition towards savannas.
C_LIO_LIColonization of the savanna tree community lagged behind the grass community, but also began to transition.
C_LI

SynthesisOur results suggest that rare extreme fires, followed by repeated burning can indeed result in savannization in places where savanna and forest represent alternative stable states. Spatial and temporal autocorrelation weave complexity in brain networks High-throughput experimental methods in neuroscience have led to an explosion of techniques for measuring complex interactions and multi-dimensional patterns. However, whether sophisticated measures of emergent phenomena can be traced back to simpler low-dimensional statistics is largely unknown. To explore this question, we examine resting state fMRI (rs-fMRI) data using complex topology measures from network neuroscience. We show that spatial and temporal autocorrelation are reliable statistics which explain numerous measures of network topology. Surrogate timeseries with subject-matched spatial and temporal autocorrelation capture nearly all reliable individual and regional variation in these topology measures. Network topology changes during aging are driven by spatial autocorrelation, and multiple serotonergic drugs causally induce the same topographic change in temporal autocorrelation. This reductionistic interpretation of widely-used complexity measures may help link them to neurobiology. The BioImage Archive - building a home for life-sciences microscopy data Despite the huge impact of data resources in genomics and structural biology, until now there has been no central archive for biological data for all imaging modalities. The BioImage Archive is a new data resource at the European Bioinformatics Institute (EMBL-EBI) designed to fill this gap. In its initial development BioImage Archive accepts bioimaging data associated with publications, in any format, from any imaging modality from the molecular to the organism scale, excluding medical imaging. The BioImage Archive will ensure reproducibility of published studies that derive results from image data and reduce duplication of effort. Most importantly, the BioImage Archive will help scientists to generate new insights through reuse of existing data to answer new biological questions, and provision of training, testing and benchmarking data for development of tools for image analysis. The archive is available at https://www.ebi.ac.uk/bioimage-archive/.

HighlightsO_LIThe BioImage Archive is a new archival data resource at the European Bioinformatics Institute (EMBL-EBI).
C_LIO_LIThe BioImage Archive aims to accept all biological imaging data associated with peer-reviewed publications using approaches that probe biological structure, mechanism and dynamics, as well as other important datasets that can serve as reference examples for particular biological or technical domains.
C_LIO_LIThe BioImage Archive aims to encourage the use of valuable imaging data, to improve reproducibility of published results that rely on image data, and to facilitate extraction of novel biological insights from existing data and development of new image analysis methods.
C_LIO_LIThe BioImage Archive forms the foundation for an ecosystem of related databases, supporting those resources with storage infrastructure and indexing across databases.
C_LIO_LIAcross this ecosystem, the BioImage Archive already stores and provides access to over 1.5 petabytes of image data from many different imaging modalities and biological domains.
C_LIO_LIFuture development of the BioImage Archive will support the fast-emerging next generation file formats (NGFFs) for bioimaging data, providing access mechanisms tailored toward modern visualisation and data exploration tools, as well as unlocking the power of modern AI-based image-analysis approaches.
C_LI Single cell genomics reveals plastid-lacking Picozoa are close relatives of red algae The endosymbiotic origin of plastids from cyanobacteria gave eukaryotes photosynthetic capabilities and launched the diversification of countless forms of algae. These primary plastids are found in members of the eukaryotic supergroup Archaeplastida. All known archaeplastids still retain some form of primary plastids, which are widely assumed to have a single origin. Here, we used single-cell genomics from natural samples combined with phylogenomics to infer the evolutionary origin of the phylum Picozoa, a globally distributed but seemingly rare group of marine microbial heterotrophic eukaryotes. Strikingly, the analysis of 43 single-cell genomes shows that Picozoa belong to Archaeplastida, specifically related to red algae and the phagotrophic rhodelphids. These picozoan genomes support the hypothesis that Picozoa lack a plastid, and further reveal no evidence of an early cryptic endosymbiosis with cyanobacteria. These findings change our understanding of plastid evolution as they either represent the first complete plastid loss in a free-living taxon, or indicate that red algae and rhodelphids obtained their plastids independently of other archaeplastids. Cell-type specific circadian bioluminescence rhythms in Dbp reporter mice Circadian rhythms are endogenously generated physiological and molecular rhythms with a cycle length of about 24 h. Bioluminescent reporters have been exceptionally useful for studying circadian rhythms in numerous species. Here, we report development of a reporter mouse generated by modification of a widely expressed and highly rhythmic gene encoding D-site albumin promoter binding protein (Dbp). In this line of mice, firefly luciferase is expressed from the Dbp locus in a Cre-recombinase-dependent manner, allowing assessment of bioluminescence rhythms in specific cellular populations. A mouse line in which luciferase expression was Cre-independent was also generated. The Dbp reporter alleles do not alter Dbp gene expression rhythms in liver or circadian locomotor activity rhythms. In vivo and ex vivo studies show the utility of the reporter alleles for monitoring rhythmicity. Our studies reveal cell-type specific characteristics of rhythms among neuronal populations within the suprachiasmatic nuclei ex vivo. In vivo studies show Dbp-driven bioluminescence rhythms in the liver of Albumin-Cre;DbpKI/+ "liver reporter" mice. After a shift of the lighting schedule, locomotor activity achieved the proper phase relationship with the new lighting cycle more rapidly than hepatic bioluminescence did. As previously shown, restricting food access to the daytime altered the phase of hepatic rhythmicity. Our model allowed assessment of the rate of recovery from misalignment once animals were provided with food ad libitum. These studies confirm the previously demonstrated circadian misalignment following environmental perturbations and reveal the utility of this model for minimally invasive, longitudinal monitoring of rhythmicity from specific mouse tissues. Microbial loop of a Proterozoic ocean analogue Meromictic Lake Cadagno, an ancient ocean analogue, is known for its permanent stratification and persistent anoxygenic microbial bloom within the chemocline. Although the anaerobic microbial ecology of the lake has been extensively studied for at least 25 years, a comprehensive picture of the microbial food web linking the bacterial layer to phytoplankton and viruses, with explicit measures of primary and secondary production, is still missing. This study sought to better understand the abundances and productivity of microbes in the context of nutrient biogeochemical cycling across the stratified zones of Lake Cadagno. Photosynthetic pigments and chloroplast 16S rRNA gene phylogenies suggested the presence of eukaryotic phytoplankton through the water column. Evidence supported high abundances of Ankyra judayi, a high-alpine adapted chlorophyte, in the oxic mixolimnion where oxygenic-primary production peaked. Through the no- and low-oxygen chemocline and monimolimnion, chlorophytes related to Closteriopsis acicularis, a known genus of meromictic lakes, and Parachlorella kessleri were observed. Chromatium, anoxygenic phototrophic sulfur bacteria, dominated the chemocline along with Lentimicrobium, a genus of known fermenters whose abundance was newly reported in Lake Cadagno. Secondary production peaked in the chemocline suggesting primary producers depend on hetrotrophs for nutrient remineralization. As previously observed, sulfur reducing bacteria (SRBs), especially Desulfocapsa and Desulfobulbus, were present in the chemocline and anoxic monimolimnion. Virus-to-microbe ratios (VMR) peaked in the zone of phytoplankton, yet were at a minimum at the peak of Chromatium. These dynamic trends suggest viruses may play a role in the modulation of oxygenic and anoxygenic photo- and chemosynthesis in Lake Cadagno and other permanently stratified systems.

ImportanceAs a window to the past, the study offers insights into the role of microbial guilds of Proterozoic ocean chemoclines in the production and recycling of organic matter of sulfur- and ammonia-containing ancient oceans. The new observations described here suggest that eukaryotic algae were persistent in the low oxygen upper-chemocline in association with purple and green sulfur bacteria in the lower half of the chemocline. Further, this study provides the first insights into Lake Cadagno viral ecology. High viral abundances suggested viruses may be essential components of the chemocline where their activity may result in the release and recycling of organic matter. The framework developed in this study through the integration of diverse geochemical and biological data types lays the foundation for future studies to quantitatively resolve the processes performed by discrete populations comprising the microbial loop in this early anoxic ocean analogue. Omicron BA.1 and BA.2 Variants Increase the Interactions of SARS-CoV-2 Spike Glycoprotein with ACE2 SARS-CoV-2 infection is initiated by binding of the receptor-binding domain (RBD) of its spike glycoprotein to the peptidase domain (PD) of angiotensin-converting enzyme 2 (ACE2) receptors in host cells. Recently detected Omicron variant of SARS-CoV-2 (B.1.1.529) is heavily mutated on RBD. Currently, the most common Omicron variants are the original BA.1 Omicron strain and the BA.2 variant, which became more prevalent since it first appeared. To investigate how these mutations affect RBD-PD interactions, we performed all-atom molecular dynamics simulations of the BA.1 and BA.2 RBD-PD in the presence of full-length glycans, explicit water and ions. Simulations revealed that RBDs of BA.1 and BA.2 variants exhibit a more dispersed interaction network and make an increased number of salt bridges and hydrophobic interactions with PD compared to wild-type RBD. Although BA.1 and BA.2 differ in two residues at the RBD-ACE2 interface, no major difference in RBD-PD interactions and binding strengths were observed between these variants. Using the conformations sampled in each trajectory, the Molecular Mechanics Poisson-Boltzmann Surface Area (MMPBSA) method estimated ~34% and ~51% stronger binding free energies for BA.1 and BA.2 RBD, respectively, than wild-type RBD, which may result in higher binding efficiency of the Omicron variant to infect host cells. FACT regulates pluripotency through distal regulation of gene expression in murine embryonic stem cells The FACT complex is a conserved histone chaperone with essential roles in transcription and histone deposition. FACT is essential in pluripotent and cancer cells, but otherwise dispensable for most mammalian cell types. FACT deletion or inhibition can block induction of pluripotent stem cells, yet the mechanism through which FACT regulates cell fate decisions remains unclear. To determine this mechanism, we used inducible depletion of FACT subunit SPT16 in murine embryonic stem cells paired with genomic factor localization, nascent transcription, and chromatin accessibility analyses. Over a timecourse of SPT16 depletion, nucleosomes invade loci bound by master pluripotency factors and gene-distal DNaseI hypersensitive sites. Simultaneously, transcription of Pou5f1 (OCT4), Sox2, Nanog, and enhancer RNAs produced at the genes' associated enhancers are downregulated, suggesting that FACT regulates expression of the pluripotency factors themselves. We find that FACT maintains cellular pluripotency through a precise nucleosome-based regulatory mechanism for appropriate expression of both coding and non-coding transcripts associated with pluripotency. The obstacles and potential clues of prime editing applications in tomato, a dicot plant Precision genome editing is highly desired for crop improvement. The recently emerged CRISPR/Cas technology offers great potential applications in precision plant genome engineering. A prime editing (PE) approach combining a reverse transcriptase (RT) with a Cas9 nickase and a "priming" extended guide RNA has shown a high frequency for precise genome modification in mammalian cells and several plant species. Nevertheless, the applications of the PE approach in dicot plants are still limited and inefficient. We designed and tested prime editors for precision editing of a synthetic sequence in a transient assay and for desirable alleles of 10 loci in tomato by stable transformation. Our data obtained by targeted deep sequencing also revealed only low PE efficiencies in both the tobacco and tomato systems. Further assessment of the activities of the PE components uncovered potential reasons for the inefficiency of the PE complexes. Modifying the pegRNA sequences by shortening or introducing mismatches to the primer binding sequences (PBS) in order to reduce their melting temperatures (Tm) did not enhance the PE efficiency at the SlBMP21, SlALC and SlALS1 loci. Our data show challenges of PE approach in tomato, indicating a further improvement of the PE system for the successful applications such as use of improved expression systems. Our work provides an important clue for the successful application of the PE approach in crop improvement. Estimation of differential cell cycle kinetics in higher plant root meristem with cellular fate and positional resolution. Plant root development is a complex spatial-temporal process that originates in the root apical meristem (RAM). To shape the organs structure signaling between the different cells and cell files must be highly coordinated. Thereby, diverging kinetics of chromatin remodeling and cell growth in these files need to be integrated and balanced by differential cell growth and local differences in cell proliferation frequency. Understanding the local differences in cell cycle duration in the RAM and its correlation with chromatin organization is crucial to build a holistic view on the different regulatory processes and requires a quantitative estimation of the chromatin geometry and underlying mitotic cell cycle phases timing at every cell file and every position. Unfortunately, so far precise methods for such analysis are missing.

This study presents a robust and straightforward pipeline to determine in parallel the duration of cell cycles key stages in all cell layers of a plants root and their nuclei organization. The methods combine marker-free techniques based on the detection of the nucleus, deep analysis of the chromatin phase transition, incorporation of 5-ethynyl-2'-deoxyuridine (EdU), and mitosis with a deep-resolution plant phenotyping platform to analyze all key cell cycle events kinetics.

In the Arabidopsis thaliana L. RAM S-phase duration was found to be as short as 20-30 minutes in all cell files. The subsequent G2-phase duration however depends on the cell type/position and varies from 3.5 hours in the pericycle to more than 4.5 hours in the epidermis. Overall, S+G2+M duration in Arabidopsis under our condition is 4 hours in the pericycle and up to 5.5 hours in the epidermis.

Endocycle duration was determined as the time required to achieve 100% EdU index in the transition zone and estimated to be in the range of 3-4 hours.

Besides Arabidopsis, we show that the presented technique is applicable also to root tips of other dicot and monocot plants (tobacco (Nicotiana tabacum L.), tomato (Lycopersicon esculentum L.) and wheat (Triticum aestivum L.). Cis-Regulatory Hubs: a new 3D Model of Complex Disease Genetics with an Application to Schizophrenia The 3-dimensional (3D) conformation of the chromatin creates complex networks of noncoding regulatory regions (distal elements) and promoters impacting gene regulation. Despite the importance of the role of noncoding regions in complex diseases, little is known about their interplay within regulatory hubs and implication in multigenic diseases like schizophrenia. Here we show that cis-regulatory hubs (CRHs) in neurons highlight functional interactions between distal elements and promoters, providing a model to explain epigenetic mechanisms involved in complex diseases. CRHs represent a new 3D model, where distal elements interact to create a complex network of active genes. In a disease context, CRHs highlighted strong enrichments in schizophrenia-associated genes, schizophrenia-associated SNPs and schizophrenia heritability compared to equivalent structures. Finally, CRHs exhibit larger proportions of genes differentially expressed in schizophrenia compared to promoter-distal element pairs or TADs. CRHs thus capture causal regulatory processes improving the understanding of complex disease etiology such as schizophrenia. These multiple lines of genetic and statistical evidence support CRHs as 3D models to study dysregulation of gene expression in complex diseases more generally.

Summary BlurbGenes and their regulatory elements are organized within 3D networks in neurons which model functional structures and explain schizophrenia genetic etiology. Non-Newtonian Blood Rheology Impacts Left Atrial Stasis in Patient-Specific Simulations The lack of mechanically effective contraction of the left atrium (LA) during atrial fibrillation (AF) disturbs blood flow, increasing the risk of thrombosis and ischemic stroke. Thrombosis is most likely in the left atrial appendage (LAA), a small narrow sac where blood is prone to stagnate. Slow flow promotes the formation of erythrocyte aggregates in the LAA, also known as rouleaux, causing viscosity gradients that are usually disregarded in patient-specific simulations. To evaluate these non-Newtonian effects, we built atrial models derived from 4D computed tomography scans of patients and carried out computational fluid dynamics simulations using the Carreau-Yasuda constitutive relation. We examined six patients, three of whom had AF and LAA thrombosis or a history of transient ischemic attacks (TIAs). We modeled the effects of hematocrit and rouleaux formation kinetics by varying the parameterization of the Carreau-Yasuda relation and modulating non-Newtonian viscosity changes based on residence time. Comparing non-Newtonian and Newtonian simulations indicates that slow flow in the LAA increases blood viscosity, altering secondary swirling flows and intensifying blood stasis. While some of these effects can be subtle when examined using instantaneous metrics like shear rate or kinetic energy, they are manifested in the blood residence time, which accumulates over multiple heartbeats. Our data also reveal that LAA blood stasis worsens when hematocrit increases, offering a potential new mechanism for the clinically reported correlation between hematocrit and stroke incidence. In summary, we submit that hematocrit-dependent non-Newtonian blood rheology should be considered in calculating patient-specific blood stasis indices by computational fluid dynamics. The first sheep graph-based pan-genome reveals the spectrum of structural variations and their effects on tail phenotypes Structural variations (SVs) are a major contributor to genetic diversity and phenotypic variations, but their prevalence and functions in domestic animals are largely unexplored. Here, we assembled 26 haplotype-resolved genome assemblies from 13 genetically diverse sheep using PacBio HiFi sequencing. We constructed a graph-based ovine pan-genome and discovered 142,422 biallelic insertions and deletions, 7,028 divergent alleles and 13,419 multiallelic variations. We then used a graph-based approach to genotype the biallelic SVs in 684 individuals from 45 domestic breeds and two wild species. Integration with RNA-seq data allows to identify candidate expression-associated SVs. We demonstrate a direct link of SVs and phenotypes by localizing the putative causative insertion in HOXB13 gene responsible for the long-tail trait and identifying multiple large SVs associated with the fat-tail. Beyond generating a benchmark resource for ovine structural variants, our study highlights that animal genetic research will greatly benefit from using a pan-genome graph rather than a single reference genome.