Age and Sex Effects on Advanced White Matter Microstructure Measures in 15,628 Older Adults: A UK Biobank Study A comprehensive characterization of the brains white matter is critical for improving our understanding of healthy and diseased aging. Here we used diffusion-weighted magnetic resonance imaging (dMRI) to estimate age and sex effects on white matter microstructure in a cross-sectional sample of 15,628 adults aged 45-80 years old (47.6% male, 52.4% female). Microstructure was assessed using the following four models: a conventional single-shell model, diffusion tensor imaging (DTI); a more advanced single-shell model, the tensor distribution function (TDF); an advanced multi-shell model, neurite orientation dispersion and density imaging (NODDI); and another advanced multi-shell model, mean apparent propagator MRI (MAPMRI). Age was modeled using a data-driven statistical approach, and normative centile curves were created to provide sex-stratified white matter reference charts. Participant age and sex substantially impacted many aspects of white matter microstructure across the brain, with the advanced dMRI models TDF and NODDI detecting such effects the most sensitively. These findings and the normative reference curves provide an important foundation for the study of healthy and diseased brain aging.
Bayesian inference of circular variables in ring attractor networks Working memories are thought to be held in attractor networks in the brain. Because working memories are often based on uncertain information, memories should ideally come with a representation of this uncertainty for strategic use in behavior. However, the attractor states that hold these memories in attractor networks commonly do not represent such uncertainty. Focusing here on ring attractor networks for encoding head direction, we show that these networks in fact feature all the motifs required to represent uncertainty in head direction estimates. Specifically, they could do so by transiently modulating their overall activity by uncertainty, in line with a circular Kalman filter that performs near-optimal statistical circular estimation. More generally, we show that ring attractors can perform near-optimal Bayesian computation if they can flexibly deviate from their attractor states. Finally, we show that the basic network motifs sufficient for such statistical inference are already known to be present in the brain. Overall, our work demonstrates that ring attractors can in principle implement a dynamic Bayesian inference algorithm in a biologically plausible manner.
Understanding the mechanism of monoADP ribosylation in OsSRT1 and its linkage to the DNA repair system under stress conditions The role of sirtuins in plants are slowly unraveling. There are only reports of H3K9Ac deacetylation by OsSRT1. Here our studies shade light on its dual enzyme capability with preference for mono ADP ribosylation over deacetylation. OsSRT1 can specifically transfer the single ADP ribose group on its substrates in an enzymatic manner. This mono ADPr effect is not well known in plants, more so for deacetylases. The products of this reaction (NAM and ADP ribose) have immense negative effect on this enzyme suggesting a tighter regulation. Resveratrol, a natural plant polyphenol proves to be a strong activator of this enzyme at 150 M concentration. Under different abiotic stress conditions, we could link this ADP ribosylase activity to the DNA repair pathway by activating the enzyme PARP1. Metal stress in plants also influences these enzyme activities.

HighlightsO_LIOsSRT1 can transfer a single moiety of ADP-ribose on itself as well as other nuclear proteins like histones H3 and H2A.
C_LIO_LINAM, ADP-ribose and certain metal ions negatively regulate this ADP-ribose transfer.
C_LIO_LIADPr of OsPARP1 and OsPARP2 links OsSRT1 to DNA damage repair pathways.
C_LIO_LIOsSRT1 positively regulates the activity of OsPARP1 by ADP ribosylating it.
C_LIO_LIOn plants exposure to H2O2 (oxidative stress) and Arsenic toxicity, there is a link between the increased activity of the players of DNA damage repair system and overexpression of OsSRT1.
C_LI 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 dynamical systems treatment of transcriptomic trajectories in hematopoiesis Inspired by Waddingtons illustration of an epigenetic landscape, cell-fate transitions have been envisioned as bifurcating dynamical systems, wherein the dynamics of an exogenous signal couples to a cells enormously complex signaling and transcriptional machinery, eliciting a qualitative transition in the collective state of a cell - its fate. It remains unclear, however, whether the dynamical systems framework can go beyond a word-based caricature of the system and provide sharp quantitative insights that further our understanding of differentiation. Single-cell RNA sequencing (scRNA-seq), which measures the distributions of possible transcriptional states in large populations of differentiating cells, provides an alternate view, in which development is marked by the individual concentration variations of a myriad of genes. Here, starting from formal mathematical derivations, we challenge these transcriptomic trajectories to a rigorous statistical evaluation of whether they display signatures consistent with bifurcations. After pinpointing bifurcations along transcriptomic trajectories of the neutrophil branch of hematopoeitic differentiation we are able to further leverage the primitive features of a linear instability to identify the single-direction in gene expression space along which the bifurcation unfolds and identify possible gene contributors. This scheme identifies transcription factors long viewed to play a crucial role in the process of neutrophil differentiation in addition to identifying a host of other novel genetic players. Most broadly speaking, we provide evidence that, though very high-dimensional, a bifurcating dynamical systems formalism might be appropriate for the process of cellular differentiation and that it can be leveraged to provide insights. Ambitiously, our work attempts to take a step beyond data-analysis and towards the construction of falsifiable mathematical models that describe the dynamics of the entire transcriptome. SPIRO - the automated Petri plate imaging platform designed by biologists, for biologists The imaging of plant seedlings, fungal mycelia and bacterial colonies grown on Petri plates is common in phenotyping assays, and is typically done manually despite the procedures being time-consuming and laborious. The main reason for this is the still limited availability of existing automated phenotyping tools and facilities. Additionally, constructing a custom-made automated solution is a daunting task for most research groups specializing in biology.

Here, we describe SPIRO, the Smart Plate Imaging Robot, an automated platform that acquires time-lapse photos of up to four vertically oriented Petri plates in a single experiment. SPIRO was designed for biologists by biologists; thus, its assembly does not require experience in engineering or programming and its operation is sufficiently intuitive to be carried out without training. SPIRO has a small footprint optimal for fitting into standard incubators for plants and microbes, and is equipped with an LED light source for imaging in the dark, thus allowing acquisition of photos under optimal growth conditions. SPIROs web-based user interface allows setting up experiments and downloading data remotely, without interfering with samples growth. SPIROs excellent image quality is optimal for automated image processing, which we demonstrate with two semi-automated assays for analysis of commonly used phenotypic traits: seed germination and root growth.

Moreover, the robot can be easily customized for a specific use, as all information about SPIRO, including the models for 3D-printed structural components, control software, and scripts for image analysis, are released under permissive open-source licenses. Can double PPO mutations exist in the same allele and are such mutants functional? BackgroundResistance to PPO-inhibiting herbicides is primarily endowed by target site mutations at the PPO2 gene that compromise binding of the herbicide to the catalytic domain. In Amaranthus ssp. PPO2, the most prevalent target mutations are deletion of the G210 codon and the R128G and G339A substitutions. These mutations strongly affect the dynamic of the PPO2 binding pocket resulting in reduced affinity with the ligand. Here we investigated the likelihood of co-occurrence of the most widespread target site mutations in the same ppo2 allele.

ResultsPlants carrying R128G+/+ {Delta}G210+/-, where + indicates presence of the mutation, were crossed with each other. The ppo2 of the offspring was subjected to pyrosequencing and E. coli- based Sanger sequencing to determine mutation frequencies and allele co-occurrence. The data show that R128G {Delta}G210 can occur only in one allele; the second allele carries only one mutation. Double mutation in both alleles is less likely because of significant loss of enzyme activity. The segregation of offspring populations derived from a cross between heterozygous plants carrying {Delta}G210 G399A also showed no co-occurrence in the same allele. The offspring exhibited the expected mutation distribution patterns with few exceptions.

ConclusionsHomozygous double-mutants are not physiologically viable. Double-mutant plants can only exist in a heterozygous state. Alternatively, if two mutations are detected in one plant, each mutation would occur in a separate allele.

NomenclaturePalmer amaranth, Amaranthus palmeri S. Wats.; protoporphyrinogen oxidase, PPO; tall waterhemp, Amaranthus tuberculatus (Moq.) J.D.Sauer Activation of the type 3 secretion system of enteropathogenic E. coli leads to remodeling of its membrane composition The cell envelope of gram-negative bacteria is a complex structure, essential for bacterial survival and for resistance to many antibiotics. Channels that cross the bacterial envelope and the host cell membrane form secretion systems that are activated upon attachment to host, enabling bacteria to inject effector molecules into the host cell, required for bacterial-host interaction. The type III secretion system (T3SS) is critical for the virulence of several pathogenic bacteria, including enteropathogenic E. coli (EPEC). The EPEC T3SS activation is associated with repression of carbon storage regulator (CsrA), resulting in gene expression remodeling, which is known to affect EPEC central carbon metabolism and contributes to the adaptation to a cell-adherent lifestyle in a poorly understood manner. We reasoned that the changes in bacterial envelope upon attachment to host and the activation of a secretion system may involve a modification of the lipid composition of bacterial envelope. Accordingly, we performed a lipidomics analysis on mutant strains that simulate T3SS activation. We saw a shift in glycerophospholipid metabolism towards the formation of lysophospholipids, attributed to corresponding upregulation of the phospholipase pldA and the acyltransferase ygiH upon T3SS activation in EPEC. We also detected a shift from menaquinones and ubiquinones to undecaprenyl lipids, concomitant to abnormal synthesis of O-antigen. The remodeling of lipid metabolism is mediated by CsrA and associated with increased bacteria cell size and Zeta potential, and a corresponding alteration in EPEC permeability to vancomycin, increasing the sensitivity of T3SS-activated strains and of adherent wild type EPEC to the antibiotic.

ImportanceThe characterization of EPEC membrane lipid metabolism upon attachment to host is an important step towards a better understanding the shift of EPEC, a notable human pathogen, from a planktonic to adherent life style. It may also apply to other pathogenic bacteria that use this secretion system. We predict that upon attachment to host cells the lipid remodeling upon T3SS activation contributes to bacterial fitness and promotes host colonization, and show that it is associated with increased cell permeability and higher sensitivity to vancomycin.

To the best of our knowledge, this is the first demonstration of a bacterial lipid remodeling due to activation of a secretion system. Reconstitution of minimal motility system based on Spiroplasma swimming by expressing two bacterial actins in synthetic minimal bacterium Motility is one of the most important features of life, but its evolutionary origin remains unknown. In this study, we focused on Spiroplasma, commensal, or parasitic bacteria. They swim by switching the helicity of a ribbon-like cytoskeleton that comprises six proteins, each of which evolved from a nucleosidase and bacterial actin called MreB. We expressed these proteins in a synthetic, non-motile minimal bacterium, JCVI-syn3.0B, whose reduced genome was computer-designed and chemically synthesized. The synthetic bacterium exhibited swimming motility with features characteristic of Spiroplasma swimming. Moreover, some combinations of the two proteins produced a helical cell shape and swimming, suggesting that the swimming originated from the differentiation and coupling of bacterial actins, and we also obtained a minimal system for motility of the synthetic bacterium.

One-Sentence SummaryThe minimal system comprised two bacterial actins that provided cell helicity and swimming to the synthetic minimal bacterium. Depletion of Sun1/2 induces heterochromatin accrual in mesenchymal stem cells during adipogenesis Critical to the mechano-regulation of mesenchymal stem cells (MSC), Linker of the Nucleoskeleton and Cytoskeleton (LINC) complex transduces cytoskeletal forces to the nuclei. The LINC complex contains outer nuclear membrane Nesprin proteins that associate with the cytoskeleton and their inner nuclear membrane couplers, Sun proteins. In addition to coupling Nesprin-associated cytoskeletal elements to inner nuclear membrane, Sun proteins also function in regulating nuclear mechanics and chromatin tethering to inner nuclear membrane. This suggests that release of LINC-mediated cytoskeletal connections from cell nuclei may have different effects on chromatin organization and MSC differentiation than those due to ablation of intra-nuclear Sun proteins. To test this hypothesis we compared Sun1/2 depletion with expression of the dominant-negative KASH domain (dnKASH) that inhibits Nesprin-Sun association. In cells cultured under adipogenic conditions, disconnecting LINC from the cytoskeleton through dnKASH expression, increased adipogenic gene expression and fat droplet formation; heterochromatin H3K27 and H3K9 methylation was unaltered. In contrast, Sun1/2 depletion inhibited adipogenic gene expression and fat droplet formation; as well the anti-adipogenic effect of Sun1/2 depletion was accompanied by increased H3K9me3, which was enriched on Adipoq, silencing this fat locus. We conclude that releasing the nucleus from cytoskeletal constraints via dnKASH accelerates adipogenesis while depletion of Sun1/2 increases heterochromatin accrual on adipogenic genes in a fashion independent of LINC complex function. Therefore, while these two approaches both disable LINC functions, their divergent effects on the epigenetic landscape indicate they cannot be used interchangeably to study mechanical regulation of cell differentiation. Niche differences, not fitness differences, explain coexistence across ecological groups Understanding the drivers of species coexistence is an important objective in ecology. Yet, the multitude of methods to study coexistence hampers cross-community comparisons. Here, we standardized niche and fitness differences (i.e how species limit themselves compared to others and their competitive ability, respectively) across 1018 species pairs to investigate species coexistence across ecological groups and methodological settings (experimental setup, natural co-occurrence, population model used, and growth method). We find that, first, coexistence is driven by large niche differences, not by small fitness differences. Second, species group into clear clusters of coexisting and non-coexisting species along the niche axis. Finally, these clusters are not driven by ecological or methodological settings. This suggests differences between coexisting and non-coexisting communities transcending those measured in our empirical systems. Overall, our results show that species coexistence is mainly influenced by mechanisms acting on niche differences. Citation needed?Wikipedia bibliometrics during the first waveof the COVID-19 pandemic. BackgroundWith the COVID-19 pandemics outbreak, millions flocked to Wikipedia for updated information. Amid growing concerns regarding an "infodemic", ensuring the quality of information is a crucial vector of public health. Investigating if and how Wikipedia remained up to date and in line with science is key to formulating strategies to counter misinformation. Using citation analyses, we asked: which sources informed Wikipedias COVID-19-related articles before and during the pandemics first wave (January-May 2020).

ResultsWe found that coronavirus-related articles referenced trusted media sources and high-quality academic research. Moreover, despite a surge in COVID-19 preprints, Wikipedia had a clear preference for open-access studies published in respected journals and made little use of preprints. Building a timeline of English COVID-19 articles from 2001-2020 revealed a nuanced trade-off between quality and timeliness. It further showed how preexisting articles on key topics related to the virus created a framework for integrating new knowledge. Supported by a rigid sourcing policy, this "scientific infrastructure" facilitated contextualization and regulated the influx of new information. Lastly, we constructed a network of DOI-Wikipedia articles, which showed the shifting landscape of pandemic-related knowledge on Wikipedia and how academic citations create a web of shared knowledge supporting topics like COVID-19 vaccine development.

ConclusionsUnderstanding how scientific research interacts with the digital knowledge-sphere during the pandemic provides insight into how Wikipedia can facilitate access to science. It also reveals how, aided by what we term its "citizen encyclopedists", it successfully fended off COVID-19 disinformation and how this unique model may be deployed in other contexts. A novel two-generation approach to the population dynamics of gregarious parasitoids Parasitoids, as important natural enemies, help maintain balance in natural ecosystems. Their population dynamics is generally predicted from the number of individuals. Here, using gregarious parasitoids as models, we show that this traditional approach omits one important parameter: the mothers manipulation of offspring fertility due to the clutch size-body size- fertility correlation. As a result of this correlation, when females deliberately adjust the number of offspring laid in a host, they determine not only the number but also the body sizes and reproductive potentials of those offspring. For the first time in the model species Anaphes flavipes, we determined the parasitoids offspring fertility from clutch size. Using this, we experimentally clarified the advantage of specific clutch size combinations and we show that identical fertility in the F1 generation can lead to distinctly different fertility values in the F2 generation. Even with the same number of hosts, lower fertility in the F1 generation can cause higher fertility in the F2 generation. Based on these results, we propose a novel two-generation approach which includes the clutch size-body size-fertility correlation. Our novel approach provides a new perspective for determining the individual fitness levels of gregarious parasitoids with new options for the modelling of parasitoid population dynamics. Synthetic reconstruction of the hunchback promoter specifies the role of Bicoid, Zelda and Hunchback in the dynamics of its transcription For over 40 years, the Bicoid-hunchback (Bcd-hb) system in the fruit fly embryo has been used as the model to study how positional information in morphogen concentration gradients is robustly translated into step-like responses. A body of quantitative comparisons between theory and experiment have since questioned the initial paradigm that the sharp hb transcription pattern emerges solely from diffusive biochemical interactions between the Bicoid transcription factor and the gene promoter region. Several alternative mechanisms have been proposed, such as additional sources of positional information, positive feedback from Hb proteins or out-of-equilibrium transcription activation.

By using the MS2-MCP RNA-tagging system and analysing in real time, the transcription dynamics of synthetic reporters for Bicoid and/or its two partners Zelda and Hunchback, we show that all the hb expression pattern features and temporal dynamics are compatible with an equilibrium model with a short decay length Bicoid activity gradient as a sole source of positional information. Meanwhile, Bicoids partners speed-up the process by different means: Zelda lowers the Bicoid concentration threshold required for transcriptional activation while Hunchback reduces burstiness and increases the polymerase firing rate. Lower limb in sprinters with larger relative mass but not larger normalized moment of inertia Inhomogeneous muscle hypertrophy in response to daily motor executions is a morphological adaptation in humans. Inhomogeneous morphological adaptation would change the mass distribution, which in turn may affect the mechanical difficulty in moving from an inertia perspective; however, the morphological adaptation from the inertia perspective has not been examined. Here we show no corresponding differences in the normalized mass and normalized moment of inertia between the sprinters and the untrained non-sprinters. We analyzed the fat- and water-separated magnetic resonance images from the lower limbs of 11 male sprinters (100 m best time of 10.44-10.83 s) and 12 untrained non-sprinters. The lower-limb relative mass was significantly larger in sprinters (18.7{+/-}0.7% body mass) than in non-sprinters (17.6{+/-}0.6% body mass), while the normalized moment of inertia of the lower limb around the hip in the anatomical position was similar between them (0.044{+/-}0.002 vs. 0.042{+/-}0.002 [a. u.]). The thigh relative mass in sprinters (12.9{+/-}0.4% body mass) was significantly larger than that in non-sprinters (11.9{+/-}0.4% body mass), whereas the shank and foot relative masses were similar. We revealed that the adaptation to sprinting in humans does not relatively increase the mechanical difficulty to swing the lower limb in terms of inertia, even though the lower-limb mass increases, reflecting muscular hypertrophy. This is suggested to be one of the significant advantages of human morphology for adaptability to explosive motor tasks.

NEW & NOTEWORTHYWe compared the lower-limb inertia properties between male sprinters and untrained non-sprinters by using the fat- and water-separated magnetic resonance images. Sprinters had lower limbs with a larger relative mass but not a larger normalized moment of inertia than non-sprinters. We suggest one of the practical advantages of human morphological characteristics for their adaptation to explosive motor tasks, which provides a novel perspective to the literature for understanding human morphological adaptability. Mouse oocytes do not contain a Balbiani body Oocytes spend the majority of their lifetime in a primordial, dormant state. Unlike many somatic cell types and mature oocytes, the cellular and molecular biology of primordial oocytes is largely unexplored. Yet, studying their cellular biology is necessary to understand the mechanisms through which oocytes maintain cellular fitness for decades, and why they eventually fail with age.

A hallmark of primordial oocytes in many species is the Balbiani body, a non-membrane bound compartment that contains the majority of mitochondria in the oocyte cytoplasm. The Balbiani body has been proposed to be essential for maintaining mitochondria in a healthy state during long-lasting dormancy, however, the architecture and function of the mammalian Balbiani body remains unknown.

Here, we develop enabling methods for live-imaging based comparative characterization of Xenopus, mouse and human primordial oocytes. We show that primordial oocytes in all three vertebrate species contain active mitochondria, Golgi apparatus and lysosomes. We further demonstrate that human and Xenopus oocytes have a Balbiani body characterized by a dense accumulation of mitochondria in their cytoplasm. However, despite previous reports, we did not find a Balbiani body in mouse oocytes. Instead, we demonstrate what was previously used as a marker for the Balbiani body in mouse primordial oocytes is in fact a ring-shaped Golgi apparatus that is not functionally associated with oocyte dormancy. Our work provides the first insights into the organisation of the cytoplasm in mammalian primordial oocytes, and clarifies relative advantages and limitations of choosing different model organisms for studying oocyte dormancy. Imbalanced segregation of recombinant haplotypes in hybrid populations reveals inter- and intrachromosomal Dobzhansky-Muller incompatibilities Dobzhansky-Muller incompatibilities (DMIs) are a major component of reproductive isolation between species. DMIs imply negative epistasis, and are exposed when two diverged populations hybridize. Mapping the locations of DMIs has largely relied on classical genetic mapping. Approaches to date are hampered by low power and the challenge of identifying DMI loci on the same chromosome, because strong initial linkage of parental haplotypes weakens statistical tests. Here, we propose new statistics to infer negative epistasis from haplotype frequencies in hybrid populations. When two divergent populations hybridize, the variance in heterozygosity at two loci decreases faster with time at DMI loci than at random pairs of loci. When two populations hybridize at near-even admixture proportions, the deviation of the observed variance from its expectation becomes negative for the DMI pair. This negative deviation enables us to detect intermediate to strong negative epistasis both within and between chromosomes. In practice, the detection window in hybrid populations depends on the demographic scenario, the recombination rate, and the strength of epistasis. When the initial proportion of the two parental populations is uneven, only strong DMIs can be detected with our method unless migration prevents parental haplotypes from being lost. We use the new statistics to infer candidate DMIs from three hybrid populations of swordtail fish. We identify numerous new DMI candidates, some of which are inferred to interact with several loci within and between chromosomes. Moreover, we discuss our results in the context of an expected enrichment in intrachromosomal over interchromosomal DMIs.

Author SummaryGenetic incompatibility in the form of (Bateson-)Dobzhansky-Muller incompatibilities (DMIs) is an important component of reproductive isolation between species. However, the evolutionary role of DMIs during the process of speciation is contentious. DMIs occur when two or more genetic variants interact to reduce their carriers fitness. Once recombination combines incompatible variants in hybrids, selection acts to remove these variants from the population. One step towards addressing the evolutionary role of DMIs is to quantify the prevalence of DMIs in incipient and hybridizing species. Here, we present statistics that are sensitive to the resulting recombinant imbalance and that can indicate the location of DMIs in hybrid genomes in various demographic scenarios. We use simulations to show that the time window during which a DMI is detectable depends on its genomic location, the severity of the DMI, and the populations demography. Importantly, our statistic distinguishes genetic associations arising due to physical linkage from those arising due to gene interactions, which allows for the inference of both inter- and intrachromosomal DMIs. Applying our statistics to three hybrid populations of swordtail fish, we confirm previously known DMIs and identify new candidate incompatibilities. LinearTurboFold: Linear-Time Global Prediction of Conserved Structures for RNA Homologs with Applications to SARS-CoV-2 The constant emergence of COVID-19 variants reduces the effectiveness of existing vaccines and test kits. Therefore, it is critical to identify conserved structures in SARS-CoV-2 genomes as potential targets for variant-proof diagnostics and therapeutics. However, the algorithms to predict these conserved structures, which simultaneously fold and align multiple RNA homologs, scale at best cubically with sequence length, and are thus infeasible for coronaviruses, which possess the longest genomes ([~]30,000 nt) among RNA viruses. As a result, existing efforts on modeling SARS-CoV-2 structures resort to single sequence folding as well as local folding methods with short window sizes, which inevitably neglect long-range interactions that are crucial in RNA functions. Here we present LinearTurboFold, an efficient algorithm for folding RNA homologs that scales linearly with sequence length, enabling unprecedented global structural analysis on SARS-CoV-2. Surprisingly, on a group of SARS-CoV-2 and SARS-related genomes, LinearTurbo-Folds purely in silico prediction not only is close to experimentally-guided models for local structures, but also goes far beyond them by capturing the end-to-end pairs between 5 and 3 UTRs ([~]29,800 nt apart) that match perfectly with a purely experimental work. Furthermore, LinearTurboFold identifies novel conserved structures and conserved accessible regions as potential targets for designing efficient and mutation-insensitive small-molecule drugs, antisense oligonucleotides, siRNAs, CRISPR-Cas13 guide RNAs and RT-PCR primers. LinearTurboFold is a general technique that can also be applied to other RNA viruses and full-length genome studies, and will be a useful tool in fighting the current and future pandemics.

Significance StatementConserved RNA structures are critical for designing diagnostic and therapeutic tools for many diseases including COVID-19. However, existing algorithms are much too slow to model the global structures of full-length RNA viral genomes. We present LinearTurboFold, a linear-time algorithm that is orders of magnitude faster, making it the first method to simultaneously fold and align whole genomes of SARS-CoV-2 variants, the longest known RNA virus ([~]30 kilobases). Our work enables unprecedented global structural analysis and captures long-range interactions that are out of reach for existing algorithms but crucial for RNA functions. LinearTurboFold is a general technique for full-length genome studies and can help fight the current and future pandemics. Inter-tissue convergence of gene expression during ageing suggests age-related loss of tissue and cellular identity Developmental trajectories of gene expression may reverse in their direction during ageing, a phenomenon previously linked to cellular identity loss. Our analysis of cerebral cortex, lung, liver and muscle transcriptomes of 16 mice, covering development and ageing intervals, revealed widespread but tissue-specific ageing-associated expression reversals. Cumulatively, these reversals create a unique phenomenon: mammalian tissue transcriptomes diverge from each other during postnatal development, but during ageing, they tend to converge towards similar expression levels, a process we term Divergence followed by Convergence, or DiCo. We found that DiCo was most prevalent among tissue-specific genes and associated with loss of tissue identity, which is confirmed using data from independent mouse and human datasets. Further, using publicly available single-cell transcriptome data, we showed that DiCo could be driven both by alterations in tissue cell type composition and also by cell-autonomous expression changes within particular cell types. Cell-to-cell variations in subcellular localisation of V2R are independent of its expression level G protein coupled receptors (GPCRs) translate the actions of hormones into intracellular signalling events. Mutations in GPCRs can prevent their correct expression and trafficking to the cell surface and cause disease. We use single cell measurements in HEK293 cells to show that the balance between endoplasmic reticulum (ER) and cell surface localisation of the Vasopressin 2 receptor (V2R) varies significantly from cell to cell. We find that mutations in the V2R affect the proportion of cells able to send this GPCR to the cell surface but do not prevent all cells in the population from correctly trafficking the mutant receptors. These findings reveal that the ability of cells to correctly traffic V2R to the cell surface depends not only on the expressed V2R mutant but also on the individual cell environment.

Significance statementMissense mutations in the Vasopressin 2 Receptor (V2R) cause Nephrogenic Diabetes Insipidus. Some of these mutations prevent correct expression and trafficking of V2R to the cell surface resulting in a loss-of-function. We show -using single cell measurements-that the balance between endoplasmic reticulum and cell surface localisation of the V2R varies significantly from cell to cell, independent from its expression level. Mutations affect the proportion of cells able to send V2R to the cell surface but do not prevent all cells in the population from correctly trafficking the mutant receptors. Hence, the ability of cells to correctly traffic V2R to the cell surface depends not only on expressed V2R mutant but also on the cell environment. A phenomics approach for in vitro antiviral drug discovery BackgroundThe current COVID-19 pandemic has highlighted the need for new and fast methods to identify novel or repurposed therapeutic drugs. Here we present a method for untargeted phenotypic drug screening of virus-infected cells, combining Cell Painting with antibody-based detection of viral infection in a single assay. We designed an image analysis pipeline for segmentation and classification of virus-infected and non-infected cells, followed by extraction of morphological properties.

ResultsWe show that the methodology can successfully capture virus-induced phenotypic signatures of MRC-5 human lung fibroblasts infected with Human coronavirus 229E (CoV-229E). Moreover, we demonstrate that our method can be used in phenotypic drug screening using a panel of nine host- and virus-targeting antivirals. Treatment with effective antiviral compounds reversed the morphological profile of the host cells towards a non-infected state.

ConclusionsThe method can be used in drug discovery for morphological profiling of novel antiviral compounds on both infected and non-infected cells. T Cell Receptor Beta Germline Variability is Revealed by Inference From Repertoire Data BackgroundT and B cell receptor (TCR, BCR) repertoires constitute the foundation of adaptive immunity. Adaptive immune receptor repertoire sequencing (AIRR-seq) is a common approach to study immune system dynamics. Understanding the genetic factors influencing the composition and dynamics of these repertoires is of major scientific and clinical importance. The chromosomal loci encoding for the variable regions of TCRs and BCRs are challenging to decipher due to repetitive elements and undocumented structural variants.

MethodsTo confront this challenge, AIRR-seq-based methods have recently been developed for B cells, enabling genotype and haplotype inference and discovery of undocumented alleles. However, this approach relies on complete coverage of the receptors variable regions, whereas most T cell studies sequence a small fraction of that region. Here, we adapted a B cell pipeline for undocumented alleles, genotype, and haplotype inference for full and partial TCR sequences. The pipeline also deals with gene assignment ambiguities, which is especially important in the analysis of data-sets of partial sequences.

ResultsWe identified 39 undocumented polymorphisms in T cell receptor Beta V (TRBV) and 31 undocumented 5 UTR sequences. A subset of these inferences was also observed using independent genomic approaches. We found that a single nucleotide polymorphism differentiating between the two documented T cell receptor Beta D2 (TRBD2) alleles is strongly associated with dramatic changes in the expressed repertoire.

ConclusionsWe reveal a rich picture of germline variability, and demonstrate how a single nucleotide polymorphism dramatically affects the composition of the whole repertoire. Our findings provide a basis for annotation of TCR repertoires for future basic and clinical studies. Mate sharing in male stump-tailed macaques as a possible case of coalition-like behavior to modify the group-wise fitness distribution In multi-male multi-female groups of animals, male reproductive success is often skewed toward top-ranking males. Coalition formation by middle- to low-rankers can be seen as a collaborative effort to modify the distribution of reproductive success within the group, so that they gain more than they would do otherwise. It has been pointed out, on the other hand, that a coalition of top-ranking males could hardly be profitable in the sense that they would gain little additional benefit from making collaborative effort. Here we report our novel observation of facultative sharing of mating opportunities among males in a wild population of stump-tailed macaques (Macaca arctoides) as a possible case of coalition-like behavior in which dominant males jointly guard females from mating with subordinate males and actively share mating opportunities within the allies. First, we report our novel observation of facultative sharing of mating opportunities in male stump-tailed macaques, where two or more males remain in close proximity to and copulate with a female in turn without contesting or sneaking. Second, considering the kind of coalition formation in which dominant males collaboratively exclude subordinates from mating competition and thereby strengthen the reproductive skew that already exists, we specify, by means of mathematical modeling, the condition for this kind of coalition formation to be selectively favored. Finally, we derive predictions about the occurrence of the coalition-like behavior depending on ecological and demographic factors, and test them empirically using data from the five groups of stump-tailed macaques in our study population. Recurrent exon-deleting activating mutations in AHR act as drivers of urinary tract cancer Bladder cancer has a high recurrence rate and low survival of advanced stage patients. Few genetic drivers of bladder cancer have thus far been identified. We performed in-depth structural variant analysis on whole-genome sequencing data of 206 metastasized urinary tract cancers. In ~10% of the patients, we identified recurrent in-frame deletions of exons 8 and 9 in the aryl hydrocarbon receptor gene (AHR{Delta}e8-9), which codes for a ligand-activated transcription factor. Pan-cancer analyses show that AHR{Delta}e8-9 is highly specific to urinary tract cancer and mutually exclusive with other bladder cancer drivers. In AHR{Delta}e8-9 the ligand-binding domain is disrupted and we show that this results in ligand-independent AHR-pathway activation. In bladder organoids, AHR{Delta}e8-9 induces a transformed phenotype that is characterized by upregulation of AHR target genes, downregulation of differentiation markers and upregulation of genes associated with stemness and urothelial cancer. Furthermore, AHR{Delta}e8-9 expression results in anchorage independent growth of bladder organoids, indicating tumorigenic potential. DNA-binding deficient AHR{Delta}e8-9 fails to induce transformation, suggesting a role for AHR target genes in the acquisition of the oncogenic phenotype. In conclusion, we show that AHR{Delta}e8-9 is a novel driver of urinary tract cancer and that the AHR pathway could be an interesting therapeutic target. Second Heart Field-derived Cells Contribute to Angiotensin II-mediated Ascending Aortopathies BackgroundThe ascending aorta is a common location for aneurysm and dissection. This aortic region is populated by a mosaic of medial and adventitial cells that are embryonically derived from either the second heart field (SHF) or the cardiac neural crest. SHF-derived cells populate areas that coincide with the spatial specificity of thoracic aortopathies. The purpose of this study was to determine whether and how SHF-derived cells contribute to ascending aortopathies.

MethodsAscending aortic pathologies were examined in patients with sporadic thoracic aortopathies and angiotensin II (AngII)-infused mice. Ascending aortas without overt pathology from AngII-infused mice were subjected to mass spectrometry assisted proteomics, and molecular features of SHF-derived cells were determined by single cell transcriptomic analyses. Genetic deletion of either low-density lipoprotein receptor-related protein 1 (Lrp1) or transforming growth factor-{beta} receptor 2 (Tgfbr2) in SHF- derived cells was conducted to examine the impact of SHF-derived cells on vascular integrity.

ResultsPathologies in human ascending aortic aneurysmal tissues were predominant in outer medial layers and adventitia. This gradient was mimicked in mouse aortas following AngII infusion that was coincident with the distribution of SHF-derived cells. Proteomics indicated that brief AngII infusion, prior to overt pathology, evoked downregulation of SMC proteins and differential expression of extracellular matrix proteins, including several LRP1 ligands. LRP1 deletion in SHF-derived cells augmented AngII-induced ascending aortic aneurysm and rupture. Single cell transcriptomic analysis revealed that brief AngII infusion decreased Lrp1 and Tgfbr2 mRNA abundance in SHF-derived cells and induced a unique fibroblast population with low abundance of Tgfbr2 mRNA. SHF-specific Tgfbr2 deletion led to embryonic lethality at E12.5 with dilatation of the outflow tract and retroperitoneal hemorrhage. Integration of proteomic and single cell transcriptomics results identified plasminogen activator inhibitor 1 (PAI1) as the most increased protein in SHF-derived SMCs and fibroblasts during AngII infusion. Immunostaining revealed a transmural gradient of PAI1 in both ascending aortas of AngII-infused mice and human ascending aneurysmal aortas that mimicked the gradient of medial and adventitial pathologies.

ConclusionSHF-derived cells exert a critical role in maintaining vascular integrity through LRP1 and TGF-{beta} signaling associated with increases of aortic PAI1.



O_FIG O_LINKSMALLFIG WIDTH=200 HEIGHT=89 SRC="FIGDIR/small/930917v6_ufig1.gif" ALT="Figure 1">
View larger version (34K):
org.highwire.dtl.DTLVardef@185e3d9org.highwire.dtl.DTLVardef@131d702org.highwire.dtl.DTLVardef@11afa42org.highwire.dtl.DTLVardef@242786_HPS_FORMAT_FIGEXP  M_FIG C_FIG Clinical PerspectiveO_ST_ABSWhat is new?C_ST_ABS- SHF-derived SMCs and fibroblasts associate with AngII-induced aortic pathologies.
- AngII induces a distinct fibroblast sub-cluster that is less abundant for mRNAs related to major extracellular components and TGF{beta} ligands and receptors, but more abundant for proliferative genes.
- TGFBR2 deletion in SHF-derived cells are embryonic lethal with significant dilatation of the outflow tract in mice.
- SHF-specific deletion of LRP1 leads to aortic pathologies in mice, supporting the importance of SHF-derived cells in maintaining ascending aortic wall integrity.


What are the clinical implications?- Heterogeneity of the embryonic origins of SMCs and fibroblasts contributes to complex mechanisms of vasculopathy formation, which should be considered when investigating the pathogenesis of thoracic aortopathies. Predictability of drug sensitivity in multi-species environments Microbes associate in nature forming complex communities, but they are often studied in purified form. Here I show that differences in antibiotic sensitivity between co-existing bacterial species can modify the diffusion of antibiotic molecules, and alter the drug efficacy measured in conventional in vitro assays. I developed a model, validated experimentally using an artificial two-species community, showing antibiotic molecules re-distribute evenly across all species if they have similar drug sensitivity. The result is lower drug efficacy with respect to conventional in vitro assays. However, when species have different sensitivities, the re-distribution of antibiotic is uneven and most sensitive species accumulate more drug. The result is increased drug efficacy against these species. Despite this discrepancy with conventional laboratory assays, my data suggests the changes in drug efficacy are predictable and provides insight into how communities can affect the diffusion of other molecules in the environment.