title: Genetic Conservation of SARS-CoV-2 RNA Replication Complex in Globally Circulating Isolates and Recently Emerged Variants from Humans and Minks Suggests Minimal Pre-Existing Resistance to Remdesivir

abstract: Remdesivir (RDV) exhibits potent antiviral activity against SARS-CoV-2 and is currently the only drug approved for the treatment of COVID-19. However, little is currently known about the potential for pre-existing resistance to RDV and the possibility of SARS-CoV-2 genetic diversification that might impact RDV efficacy as the virus continue to spread globally. In this study, >90,000 SARS-CoV-2 sequences from globally circulating clinical isolates, including sequences from recently emerged United Kingdom and South Africa variants, and >300 from mink isolates were analyzed for genetic diversity in the RNA replication complex (nsp7, nsp8, nsp10, nsp12, nsp13, and nsp14) with a focus on the RNA-dependent RNA polymerase (nsp12), the molecular target of RDV. Overall, low genetic variation was observed with only 12 amino acid substitutions present in the entire RNA replication complex in [≥]0.5% of analyzed sequences with the highest overall frequency (82.2%) observed for nsp12 P323L that consistently increased over time. Low sequence variation in the RNA replication complex was also observed among the mink isolates. Importantly, the coronavirus Nsp12 mutations previously selected in vitro in the presence of RDV were identified in only 2 isolates (0.002%) within all the analyzed sequences. In addition, among the sequence variants observed in [≥]0.5% clinical isolates, including P323L, none were located near the established polymerase active site or sites critical for the RDV mechanism of inhibition. In summary, the low diversity and high genetic stability of the RNA replication complex observed over time and in the recently emerged SARS-CoV-2 variants suggests a minimal global risk of pre-existing SARS-CoV-2 resistance to RDV.

category: genomics


title: Vertical inheritance governs biosynthetic gene cluster evolution and chemical diversification

abstract:
While specialized metabolites are thought to mediate ecological interactions, the evolutionary processes driving their diversification, particularly among closely related lineages, remain poorly understood. Here, we examine the evolutionary dynamics governing the distribution of natural product biosynthetic gene clusters (BGCs) using 118 strains within the marine actinomycete genus Salinispora. While previous evidence indicated that horizontal gene transfer (HGT) largely contributed to BGC diversity, we find that a majority of BGCs in Salinispora genomes are conserved through processes of vertical descent. In particular, vertical inheritance maintained BGCs over evolutionary timescales (millions of years) allowing for BGC diversification among Salinispora species. By coupling the genomic analyses with high-resolution tandem mass spectrometry (LC-MS/MS), we identified that BGC evolution in Salinispora proceeds largely through gene gain/loss events and constrained recombination that contributes to interspecies diversity at the gene, pathway, and metabolite levels. Consequently, the evolutionary processes driving BGC diversification had direct consequences for compound production and contributed to chemical diversification, as exemplified in our case study of the medically relevant proteosome inhibitors, the salinosporamides. Together, our results support the concept that specialized metabolites, and their cognate BGCs, represent functional traits associated with niche differentiation among Salinispora species.

GRAPHICAL ABSTRACT

O_FIG O_LINKSMALLFIG WIDTH=200 HEIGHT=107 SRC="FIGDIR/small/423547v2_ufig1.gif" ALT="Figure 1">
View larger version (21K):
org.highwire.dtl.DTLVardef@aef190org.highwire.dtl.DTLVardef@17536baorg.highwire.dtl.DTLVardef@5bf336org.highwire.dtl.DTLVardef@bc7d2f_HPS_FORMAT_FIGEXP  M_FIG C_FIG O_TEXTBOXSIGNIFICANCENatural products are traditionally exploited for their pharmaceutical potential; yet what is often overlooked is that the evolution of the biosynthetic gene clusters (BGCs) encoding these small molecules likely affects the diversification of the produced compounds and implicitly has an impact on the compounds activities and ecological functions. And while the prevailing dogma in natural product research attributes frequent and widespread horizontal gene transfer (HGT) as an integral driver of BGC evolution, we find that the majority of BGC diversity derives from processes of vertical descent, with HGT events being rare. This understanding can facilitate informed biosynthetic predictions to identify novel natural products, in addition to uncovering how these specialized metabolites contribute to the environmental distribution of microbes.

C_TEXTBOX

category: microbiology


title: Long term survival of Dehalococcoides mccartyi strains in mixed cultures under electron acceptor and ammonium limitation

abstract:
Few strains of Dehalococcoides mccartyi harbour and express the vinyl chloride reductase (VcrA) that catalyzes the dechlorination of vinyl chloride (VC), a carcinogenic soil and groundwater contaminant. The vcrA operon is found on a Genomic Island (GI) and therefore believed to participate in horizontal gene transfer. To try to induce horizontal gene transfer of the vcrA-GI, we blended two enrichment cultures in medium without ammonium while providing VC. We hypothesized that these conditions would select for a mutant strain of D. mccartyi that could both fix nitrogen and respire VC. However, after more than 4 years of incubation, we found no evidence for horizontal gene transfer of the vcrA-GI. Rather, we observed VC-dechlorinating activity attributed to the trichloroethene reductase TceA. Sequencing and protein modelling revealed a mutation in the predicted active site of TceA which may have influenced substrate specificity. We also identified two nitrogen-fixing D. mccartyi strains in the KB-1 culture. The presence of multiple strains of D. mccartyi with distinct phenotypes is a feature of natural environments and certain enrichment cultures (such as KB-1) and may enhance bioaugmentation success. The fact that multiple distinct strains persist in the culture for decades and that we could not induce horizontal gene transfer of the vcrA-GI suggests that it is not as mobile as predicted, or that mobility is restricted in ways yet to be discovered to specific sub-clades of Dehalococcoides.

TOC Art

O_FIG O_LINKSMALLFIG WIDTH=200 HEIGHT=106 SRC="FIGDIR/small/423565v2_ufig1.gif" ALT="Figure 1">
View larger version (36K):
org.highwire.dtl.DTLVardef@80423forg.highwire.dtl.DTLVardef@7d6484org.highwire.dtl.DTLVardef@5223b2org.highwire.dtl.DTLVardef@a58aab_HPS_FORMAT_FIGEXP  M_FIG C_FIG


title: Naa12 compensates for Naa10 in mice in the amino-terminal acetylation pathway

abstract: Amino-terminal acetylation is catalyzed by a set of N-terminal acetyltransferases (NATs). The NatA complex (including X-linked Naa10 and Naa15) is the major acetyltransferase, with 40-50% of all mammalian proteins being potential substrates. However, the overall role of amino-terminal acetylation on a whole-organism level is poorly understood, particularly in mammals. Male mice lacking Naa10 show no globally apparent in vivo amino-terminal acetylation impairment and do not exhibit complete embryonic lethality. Rather Naa10 nulls display increased neonatal lethality, and the majority of surviving undersized mutants exhibit a combination of hydrocephaly, cardiac defects, homeotic anterior transformation, piebaldism and urogenital anomalies. Naa12 is a previously unannotated Naa10-like paralogue with NAT activity that genetically compensates for Naa10. Mice deficient for Naa12 have no apparent phenotype, whereas mice deficient for Naa10 and Naa12 display embryonic lethality. The discovery of Naa12 adds to the currently known machinery involved in amino-terminal acetylation in mice.

category: developmental biology


title: The evolution of social parasitism in Formica ants revealed by a global phylogeny

abstract: Studying the behavioral and life history transitions from a cooperative, eusocial life history to exploitative social parasitism allows for deciphering the conditions under which changes in behavior and social organization lead to diversification. The Holarctic ant genus Formica is ideally suited for studying the evolution of social parasitism because half of its 176 species are confirmed or suspected social parasites, which includes all three major classes of social parasitism known in ants. However, the life-history transitions associated with the evolution of social parasitism in this genus are largely unexplored. To test competing hypotheses regarding the origins and evolution of social parasitism, we reconstructed the first global phylogeny of Formica ants and representative formicine outgroups. The genus Formica originated in the Old World during the Oligocene ([~]30 Ma ago) and dispersed multiple times to the New World. Within Formica, the capacity for dependent colony foundation and temporary social parasitism arose once from a facultatively polygynous, independently colony founding ancestor. Within this parasitic clade, dulotic social parasitism evolved once from a facultatively temporary parasitic ancestor that likely practiced colony budding frequently. Permanent social parasitism evolved twice from temporary social parasitic ancestors that rarely practiced colony budding, demonstrating that obligate social parasitism can originate from different facultative parasitic backgrounds in socially polymorphic organisms. In contrast to inquiline ant species in other genera, the high social parasite diversity in Formica likely originated via allopatric speciation, highlighting the diversity of convergent evolutionary trajectories resulting in nearly identical parasitic life history syndromes.

category: evolutionary biology


title: Application of the Anatomical Fiducials Framework to a Clinical Dataset of Patients with Parkinsons Disease

abstract: Establishing spatial correspondence between subject and template images is necessary in neuroimaging research and clinical applications such as brain mapping and stereotactic neurosurgery. Our anatomical fiducials (AFIDs) framework has recently been validated to serve as a quantitative measure of image registration based on salient anatomical features. In this study, we sought to apply the AFIDs protocol to the clinic, focusing on structural magnetic resonance images obtained from patients with Parkinsons Disease (PD). We confirmed AFIDs could be placed to millimetric accuracy in the PD dataset with results comparable to those in normal control subjects. We evaluated subject-to-template registration using this framework by aligning the clinical scans to standard template space using a robust open preprocessing workflow. We found that registration errors measured using AFIDs were higher than previously reported, suggesting the need for optimization of image processing pipelines for clinical grade datasets. Finally, we examined the utility of using point-to-point distances between AFIDs as a morphometric biomarker of PD, finding evidence of reduced distances between AFIDs that circumscribe regions known to be affected in PD including the substantia nigra. Overall, we provide evidence that AFIDs can be successfully applied in a clinical setting and utilized to provide localized and quantitative measures of registration error. AFIDs provide clinicians and researchers with a common, open framework for quality control and validation of spatial correspondence and the location of anatomical structures, facilitating aggregation of imaging datasets and comparisons between various neurological conditions.

category: neuroscience


title: Genome evolution in bacteria isolated from million-year-old subseafloor sediments

abstract:
Beneath the seafloor, microbial life subsists in isolation from the surface world under persistent energy limitation. The nature and extent of genomic evolution in subseafloor microbes has been unknown. Here we show that the genomes of Thalassospira bacterial populations cultured from million-year-old subseafloor sediments evolve by point mutation, with a relatively low rate of homologous recombination and a high frequency of pseudogenes. Ratios of synonymous to non-synonymous mutation rates correlate with the accumulation of pseudogenes, consistent with a dominant role for genetic drift in the subseafloor strains, but not in type strains of Thalassospira isolated from the surface world. Our findings demonstrate that the long term physical isolation of these bacteria, in the absence of recombination, has resulted in clonal populations that evolve consistent with  Mullers Ratchet, whereby reduced access to novel genetic material from neighbors has resulted in fixation of new mutations that accumulate in genomes over millions of years.

Significance statementThe nature and extent of genomic evolution in subseafloor microbial populations subsisting for millions of years below the seafloor is unknown. Subseafloor populations have ultra-slow metabolic rates that are hypothesized to restrict reproduction and, consequently, the spread of new traits. Our findings demonstrate that genomes of cultivated bacterial strains from the genus Thalassospira isolated from million-year-old abyssal sediment exhibit greatly reduced levels of homologous recombination, elevated numbers of pseudogenes, and genome-wide evidence of relaxed purifying selection. These substitutions and pseudogenes are fixed into the population, suggesting the genome evolution of these bacteria has been dominated by genetic drift, whereby under long-term physical isolation in small population sizes, and in the absence of homologous recombination, newly acquired mutations accumulate in the genomes of clonal populations over millions of years.


title: Cell patterning by secretion-induced plasma membrane flows

abstract:
Cells self-organize using reaction-diffusion and fluid-flow principles. Whether bulk membrane flows contribute to cell patterning has not been established. Here, using mathematical modelling, optogenetics and synthetic probes, we show that polarized exocytosis causes lateral membrane flows away from regions of membrane insertion. Plasma membrane-associated proteins with sufficiently low diffusion and/or detachment rates couple to the flows and deplete from areas of exocytosis. In rod-shaped fission yeast cells, zones of Cdc42 GTPase activity driving polarized exocytosis are limited by GTPase activating proteins (GAPs). We show that membrane flows pattern the GAP Rga4 distribution and coupling of a synthetic GAP to membrane flows is sufficient to establish the rod shape. Thus, membrane flows induced by Cdc42-dependent exocytosis form a negative feedback restricting the zone of Cdc42 activity.

One Sentence SummaryExocytosis causes bulk membrane flows that drag associated proteins and form a negative feedback restricting the exocytic site.

category: cell biology


title: Alterations in excitatory and inhibitory synaptic development within the mesolimbic dopamine pathway in a mouse model of prenatal drug exposure

abstract: The rise in rates of opioid abuse in recent years has led to an increase in the incidence of neonatal abstinence syndrome (NAS). Despite having a greater understanding of NAS and its symptoms, there still remains a lack of information surrounding the long-term effects of prenatal exposure to drugs of abuse on neurological development. One potential outcome of prenatal drug exposure that has been increasingly explored is disruption in normal synaptogenesis within the central nervous system. Both opioids and gabapentin, an antiepileptic drug commonly co-abused by opioid abuse disorder patients, have been shown to interfere with the normal functioning of astrocytes, non-neuronal glial cells known to serve many functions, including regulation of synaptic development. The goal of this study was to investigate the effects of prenatal drug exposure on synaptogenesis within brain regions associated with the mesolimbic dopamine pathway, the primary reward pathway within the brain associated with drug abuse and addiction, in a pregnant mouse model. Immunohistochemistry (IHC) and confocal fluorescence microscopy imaging studies on the brains of postnatal day 21 (P21) mouse pups revealed a significant increase in the mean number of excitatory synapses within the anterior cingulate cortex (ACC), nucleus accumbens (NAc), and prefrontal cortex (PFC) in mice that were prenatally exposed to either the opioid drug buprenorphine or gabapentin. These studies also revealed a significant decrease in the mean number of inhibitory synapses within the NAc and PFC of mice treated with buprenorphine. This observed net increase in excitatory signaling capability within the developing mesolimbic dopamine pathway suggests that exposure to drugs of abuse in utero can trigger maladaptive neuronal connectivity that persists beyond the earliest stages of life.


title: Predicting speech from a cortical hierarchy of event-based timescales

abstract: How can anticipatory neural processes structure the temporal unfolding of context in our natural environment? We here provide evidence for a neural coding scheme that sparsely updates contextual representations at the boundary of events and gives rise to a hierarchical, multi-layered organization of predictive language comprehension. Training artificial neural networks to predict the next word in a story at five stacked timescales and then using model-based functional MRI, we observe a sparse, event-based "surprisal hierarchy". The hierarchy evolved along a temporo-parietal pathway, with model-based surprisal at longest timescales represented in inferior parietal regions. Along this hierarchy, surprisal at any given timescale gated bottom-up and top-down connectivity to neighbouring timescales. In contrast, surprisal derived from a continuously updated context influenced temporo-parietal activity only at short timescales. Representing context in the form of increasingly coarse events constitutes a network architecture for making predictions that is both computationally efficient and semantically rich.


title: Memory precision of object-location binding is unimpaired in APOE ε4-carriers with spatial navigation deficits

abstract: Research suggests that tests of memory fidelity, feature binding and spatial navigation are promising for early detection of subtle behavioural changes related to Alzheimers disease (AD). In the absence of longitudinal data, one way of testing the early detection potential of cognitive tasks is through the comparison of individuals at different genetic risk for AD. Most studies have done so using samples aged 70 years or older. Here, we tested whether memory fidelity of long-term object-location binding may be a sensitive marker even among cognitively healthy individuals in their mid-60s by comparing participants at low and higher risk based on presence of the {varepsilon}4-allele of the apolipoprotein gene (n=26 {varepsilon}3{varepsilon}3, n=20 {varepsilon}3{varepsilon}4 carriers). We used a continuous report paradigm in a visual memory task that required participants to recreate the spatial position of objects in a scene. We employed mixture modelling to estimate the two distinct memory processes that underpin the trial-by-trial variation in localisation errors: retrieval success which indexes the proportion of trials where participants recalled any information about an objects position and the precision with which participants retrieved this information. Prior work has shown that these memory paradigms that separate retrieval success from precision are capable of detecting subtle differences in mnemonic fidelity even when retrieval success could not. Nonetheless, a Bayesian analysis found good evidence that {varepsilon}3{varepsilon}4 carriers did not remember fewer object locations (F(1, 42)=.450, p=.506, BF01 =3.02), nor was their precision for the spatial position of objects reduced compared to {varepsilon}3{varepsilon}3 carriers (F(1, 42)=.12, p=.726, BF01 =3.19). Because the participants in the sample presented here were a subset of a study on APOE effects on spatial navigation in the Sea Hero Quest game (Coughlan et al., 2019. PNAS, 116(9)), we obtained these data to contrast APOE effects on the two tasks within the same sample (n=33). Despite the smaller sample size, wayfinding deficits among {varepsilon}3{varepsilon}4 could be replicated (F(1, 33) =5.60, p=.024, BF10 =3.44). Object-location memory metrics and spatial navigation scores were not correlated (all r<.25, p>.1, 0<BF10 <3). These findings show spared object-location binding in the presence of a detrimental APOE {varepsilon}4 effect on spatial navigation. This suggests that the sensitivity of memory fidelity and binding tasks may not extend to individuals with one {varepsilon}4-allele in their early to mid-60s. The results provide further support to prior proposals that spatial navigation may be a sensitive marker for the earliest AD-dependent cognitive changes, even before episodic memory.


title: Polyamine biosynthesis and eIF5A hypusination are modulated by the DNA tumor virus KSHV and promote KSHV viral infection

abstract: Polyamines are critical metabolites involved in various cellular processes and often dysregulated in cancers. Kaposis sarcoma associated Herpesvirus (KSHV) is a defined oncogenic virus belonging to the sub-family of human gamma-herpesviruses. KSHV infection leads to the profound alteration of host metabolic landscape to favor the development of KSHV-associated malignancies. In our studies, we identified that polyamine biosynthesis and eIF5A hypusination are dynamically regulated by KSHV infection likely through the modulation of key enzymes of these pathways, such as ODC1, and that in return these metabolic pathways are required for both KSHV lytic switch from latency and de novo infection. The further analysis unraveled that translation of critical KSHV latent and lytic proteins (LANA, RTA) depends on eIF5A hypusination. We also demonstrated that KSHV infection can be efficiently and specifically suppressed by using inhibitors targeting either polyamine biosynthesis or eIF5A hypusination. Above all, our results illustrated that the dynamic and profound interaction of a DNA tumor virus (KSHV) with host polyamine biosynthesis and eIF5A hypusination metabolic pathways promote viral propagation and oncogenesis, which serve as new therapeutic targets to treat KSHV-associated malignancies.


title: Linker histone H1.8 inhibits chromatin-binding of condensins and DNA topoisomerase II to tune chromosome length and individualization

abstract: DNA loop extrusion by condensins and decatenation by DNA topoisomerase II (topo II) are thought to drive mitotic chromosome compaction and individualization. Here, we reveal that the linker histone H1.8 antagonizes condensins and topo II to shape mitotic chromosome organization. In vitro chromatin reconstitution experiments demonstrate that H1.8 inhibits binding of condensins and topo II to nucleosome arrays. Accordingly, H1.8 depletion in Xenopus egg extracts increased condensins and topo II levels on mitotic chromatin. Chromosome morphology and Hi-C analyses suggest that H1.8 depletion makes chromosomes thinner and longer through shortening the average loop size and reducing the DNA amount in each layer of mitotic loops. Furthermore, excess loading of condensins and topo II to chromosomes by H1.8 depletion causes hyper-chromosome individualization and dispersion. We propose that condensins and topo II are essential for chromosome individualization, but their functions are tuned by the linker histone to keep chromosomes together until anaphase.


title: Disabling de novo DNA methylation in embryonic stem cells allows an illegitimate fate trajectory

abstract:
Genome remethylation is essential for mammalian development but specific reasons are unclear. Here we examined embryonic stem (ES) cell fate in the absence of de novo DNA methyltransferases. We observed that ES cells deficient for both Dnmt3a and Dnmt3b are rapidly eliminated from chimaeras. On further investigation we found that in vivo and in vitro the formative pluripotency transition is derailed towards production of trophoblast. This aberrant trajectory is associated with failure to suppress activation of Ascl2. Ascl2 encodes a bHLH transcription factor expressed in placenta. Misexpression of Ascl2 in ES cells provokes transdifferentiation to trophoblast-like cells. Conversely, Ascl2 deletion rescues formative transition of Dnmt3a/b mutants and improves contribution to chimaeric epiblast. Thus, de novo DNA methylation safeguards against ectopic activation of Ascl2. However, Dnmt3a/b-deficient cells remain defective in ongoing embryogenesis. We surmise that multiple developmental transitions may be secured by DNA methylation silencing potentially disruptive genes.

SIGNIFICANCE STATEMENTMammalian DNA is widely modified by methylation of cytosine residues. This modification is added to DNA during early development. If methylation is prevented, the embryo dies by mid-gestation with multiple abnormalities. In this study we found that stem cells lacking the DNA methylation enzymes do not differentiate efficiently into cell types of the embryo and are diverted into producing placental cells. This switch in cell fate is driven by a transcription factor, Ascl2, which should only be produced in placenta. In the absence of DNA methylation, the Ascl2 gene is mis-expressed. Removing Ascl2 redirects embryonic fate but not full differentiation potential, suggesting that methylation acts at multiple developmental transitions to restrict activation of disruptive genes.


title: Quantifying shifts in natural selection on codon usagebetween protein regions: A population genetics approach

abstract:
BackgroundCodon usage bias (CUB), the non-uniform usage of synonymous codons, occurs across all domains of life. Adaptive CUB is hypothesized to result from selection for efficient ribosome elongation, accurate translation, mRNA secondary structure, and/or protein folding. Given the critical link between protein folding and protein function, numerous studies have analyzed the relationship between codon usage and protein structure. The results from these studies have often been contradictory, likely reflecting the differing methods used for measuring codon usage and the failure to appropriately control for confounding factors, such as differences in amino acid usage between protein structures and changes in the frequency of different structures with gene expression.

ResultsHere we take an explicit population genetics approach to quantify codon-specific shifts in natural selection related to protein structure. As a result, our approach avoids the problem of needing to choose a specific metric of codon usage, naturally takes amino acid usage of a region into account, and explicitly takes variation in gene expression into account. Our results reveal a weak relationship between codon usage and protein structure, indicating that differences in selection between structures are subtle and/or intermittent. While the magnitude of differences in selection are slight, our results indicate the relationship between codon usage and protein structure is more complex than previously believed.

ConclusionsWhile our results clearly indicate that selection on codon usage shifts between secondary structure, the overall impact of these shifts on codon usage are small. Nevertheless, our work demonstrates the statistical power and benefits of studying selective shifts on codon usage or other genomic features from an explicitly evolutionary approach. Limitations of this approach are discussed.


title: Evolutionary coupling range varies widely among enzymes depending on selection pressure

abstract:
Recent studies proposed that enzyme active sites induce evolutionary constraints at long distances. The physical origin of such long-range evolutionary coupling is unknown. Here, I use a recent biophysical model of evolution to study the relationship between physical and evolutionary couplings on a diverse data set of monomeric enzymes. I show that evolutionary coupling is not universally long-range. Rather, range varies widely among enzymes, from 2[A] to 20[A]. Furthermore, the evolutionary coupling range of an enzyme does not inform on the underlying physical coupling, which is short-range for all enzymes. Rather, evolutionary coupling range is determined by functional selection pressure.

SIGNIFICANCEUntil recently, only residues near enzyme active sites were thought to be evolutionarily constrained. However, recent studies proposed that active sites induce long-range evolutionary constraints. This seems to conflict with the common finding that physical couplings in proteins are short-range. This raises the question of how short-range physical couplings may cause long-range evolutionary couplings. Here, I show that the function that maps physical coupling into evolutionary coupling depends on functional selection pressure. Under weak selection, both couplings are similarly short-range; under strong selection, short-range physical coupling is non-linearly turned into long-range evolutionary coupling. Thus, due to a huge variation of selection pressure, evolutionary coupling range varies widely among enzymes, from very short (2 [A]) to very long (20 [A]).


title: QTL Mapping in Outbred Tetraploid (and Diploid) Diallel Populations

abstract: Over the last decade, multiparental populations have become a mainstay of genetics research in diploid species. Our goal was to extend this paradigm to autotetraploids by developing software for quantitative trait locus (QTL) mapping in connected F1 populations derived from a set of shared parents. For QTL discovery, phenotypes are regressed on the dosage of parental haplotypes to estimate additive effects. Statistical properties of the model were explored by simulating half-diallel diploid and tetraploid populations with different population sizes and numbers of parents. Across scenarios, the number of progeny per parental haplotype (pph) largely determined the statistical power for QTL detection and accuracy of the estimated haplotype effects. Multi-allelic QTL with heritability 0.2 were detected with 90% probability at 25 pph and genome-wide significance level 0.05, and the additive haplotype effects were estimated with over 90% accuracy. Following QTL discovery, the software enables a comparison of models with multiple QTL and non-additive effects. To illustrate, we analyzed potato tuber shape in a half-diallel population with 3 tetraploid parents. A well-known QTL on chromosome 10 was detected, for which the inclusion of digenic dominance lowered the Deviance Information Criterion (DIC) by 17 points compared to the additive model. The final model also contained a minor QTL on chromosome 1, but higher order dominance and epistatic effects were excluded based on the DIC. In terms of practical impacts, the software is already being used to select offspring based on the effect and dosage of particular haplotypes in breeding programs.

category: genetics


title: Sensitive detection of DNA contamination in tumor samples via microhaplotypes

abstract: Low levels of sample contamination with other human DNAs can have disastrous effects on the accurate identification of somatic variation in tumor samples. Detection of sample contamination in DNA is often based on low frequency variants that indicate if more than a single source of DNA is present. This strategy works with standard DNA samples but can be problematic in solid tumor FFPE samples because there are often huge variations in allele frequency (AF) due to copy number changes arising from gains and losses across the genome. The variable AFs make detection of contamination challenging. To avoid this, we counted microhaplotypes to assess sample contamination. Microhaplotypes are sets of variants on the same sequencing read that can be unambiguously phased. Instead of measuring AF, the number of microhaplotypes is determined. Contamination detection becomes based on fundamental genomic properties, linkage disequilibrium (LD) and the diploid nature of human DNA, rather than variant frequencies. We optimized microhaplotype panel content and selected 164 SNV sets located in regions already being sequenced within a cancer panel. Thus, contamination detection uses existing sequence data. LD data from the 1000 Genomes Project is used to make the panel ancestry agnostic, providing the same sensitivity for contamination detection with samples from individuals of African, East Asian, and European ancestry. Detection of 1% contamination with no matching normal sample is possible. The methods described here can also be extended to other DNA mixtures such as forensic and non-invasive prenatal testing samples where DNA mixes can be similarly detected. The microhaplotype method allows sensitive detection of DNA contamination in FFPE tumor and other samples when deep coverage with Illumina or other high accuracy NGS is used.


title: Co-translational folding of the first transmembrane domain of ABC-transporter CFTR is supported by assembly with the first cytosolic domain

abstract: ABC-transporters transport a wealth of molecules across membranes and consist of transmembrane and cytosolic domains. Their activity cycle involves a tightly regulated and concerted domain choreography. Regulation is driven by the cytosolic domains and function by the transmembrane domains. Folding of these polytopic multidomain proteins to their functional state is a challenge for cells, which is mitigated by co-translational and sequential events. We here reveal the first stages of co-translational domain folding and assembly of CFTR, the ABC-transporter defective in the most abundant rare inherited disease cystic fibrosis. We have combined biosynthetic radiolabeling with protease-susceptibility assays and domain-specific antibodies. The most N-terminal domain, TMD1 (transmembrane domain 1), folds both its hydrophobic and soluble helices during translation: the transmembrane helices pack tightly and the cytosolic N- and C-termini assemble with the first cytosolic helical loop ICL1, leaving only ICL2 exposed. This N-C-ICL1 assembly is strengthened by two independent events: i) assembly of ICL1 with the N-terminal subdomain of the next domain, cytosolic NBD1 (nucleotide-binding domain 1); and ii) in the presence of corrector drug VX-809, which rescues cell-surface expression of a range of disease-causing CFTR mutants. Both lead to increased shielding of the CFTR N-terminus, and their additivity implies different modes of action. Early assembly of NBD1 and TMD1 is essential for CFTR folding and positions both domains for the required assembly with TMD2. Altogether, we have gained insights into this first, nucleating, VX-809-enhanced domain-assembly event during and immediately after CFTR translation, involving structures conserved in type-I ABC exporters.

category: biochemistry


title: Implications of Spike-glycoprotein processing at S1/S2 by Furin, at S2' by Furin and/or TMPRSS2 and shedding of ACE2: cell-to-cell fusion, cell entry and infectivity of SARS-CoV-2

abstract:
The spike (S)-protein of SARS-CoV-2 binds ACE2 and requires proteolytic "priming" at PRRAR685{downarrow} into S1 and S2 (cleavage at S1/S2), and "fusion-activation" at a S2 site for viral entry. In vitro, Furin cleaved peptides mimicking the S1/S2 cleavage site more efficiently than at the putative S2, whereas TMPRSS2 inefficiently cleaved both sites. In HeLa cells Furin-like enzymes mainly cleaved at S1/S2 during intracellular protein trafficking, and S2 processing by Furin at KPSKR815{downarrow} was strongly enhanced by ACE2, but not for the optimized S2 KRRKR815{downarrow} mutant (S2), whereas individual/double KR815AA mutants were retained in the endoplasmic reticulum. Pharmacological Furin-inhibitors (Boston Pharmaceuticals, BOS-inhibitors) effectively blocked endogenous S-protein processing in HeLa cells. Furthermore, we show using pseudotyped viruses that while entry by a "pH-dependent" endocytosis pathway in HEK293 cells did not require Furin processing at S1/S2, a "pH-independent" viral entry in lung-derived Calu-3 cells was sensitive to inhibitors of Furin (BOS) and TMPRSS2 (Camostat). Consistently, these inhibitors potently reduce infectious viral titer and cytopathic effects, an outcome enhanced when both compounds were combined. Quantitative analyses of cell-to-cell fusion and spike processing revealed the key importance of the Furin sites for syncytia formation. Our assays showed that TMPRSS2 enhances fusion and proteolysis at S2 in the absence of cleavage at S1/S2, an effect that is linked to ACE2 shedding by TMPRSS2. Overall, our results indicate that Furin and TMPRSS2 play synergistic roles in generating fusion-competent S-protein, and in promoting viral entry, supporting the combination of Furin and TMPRSS2 inhibitors as potent antivirals against SARS-CoV-2.

IMPORTANCESARS-CoV-2 is the etiological agent of COVID-19 that resulted in >5 million deaths. The spike protein (S) of the virus directs infection of the lungs and other tissues by binding the angiotensin-converting enzyme 2 (ACE2) receptor. For effective infection, the S-protein is cleaved at two sites: S1/S2 and S2. Cleavage at S1/S2, induces a conformational change favoring the recognition of ACE2. The S2 cleavage is critical for cell-to-cell fusion and virus entry into host cells. Our study contributes to a better understanding of the dynamics of interaction between Furin and TMPRSS2 during SARS-CoV-2 entry and suggests that the combination of a non-toxic Furin inhibitor with a TMPRSS2 inhibitor could significantly reduce viral entry in lung cells, as evidenced by an average synergistic [~]95% reduction of viral infection. This represents a powerful novel antiviral approach to reduce viral spread in individuals infected by SARS-CoV-2 or future related coronaviruses.


title: SARS-CoV-2 spike-glycoprotein processing at S1/S2 and S2' and shedding of the ACE2 viral receptor: roles of Furin and TMPRSS2 and implications for viral infectivity and cell-to-cell fusion


title: Endosomal escape of delivered mRNA from endosomal recycling tubules visualized at the nanoscale

abstract: Delivery of exogenous mRNA using lipid nanoparticles (LNP) is a promising strategy for therapeutics. However, a bottleneck remains the poor understanding of the parameters that correlate with endosomal escape vs. cytotoxicity. To address this problem, we compared the endosomal distribution of six LNP-mRNA formulations of diverse chemical composition and efficacy, similar to those employed in mRNA-based vaccines, in primary human adipocytes, fibroblasts and HeLa cells. Surprisingly, we found that total uptake is not a sufficient predictor of delivery and different LNP vary considerably in endosomal distributions. Prolonged uptake impaired endosomal acidification, a sign of cytotoxicity, and caused mRNA to accumulate in compartments defective in cargo transport and unproductive for delivery. In contrast, early endocytic/recycling compartments have the highest probability for mRNA escape. By super-resolution microscopy we could resolve single LNP-mRNA within sub-endosomal compartments and capture events of mRNA escape from endosomal recycling tubules. Our results change the view of the mechanisms of endosomal escape and define quantitative parameters to guide the development of mRNA formulations towards higher efficacy and lower cytotoxicity.

category: bioengineering


title: The ratio of exhausted to resident infiltrating lymphocytes is prognostic for colorectal cancer patient outcome

abstract: Immunotherapy success in colorectal cancer (CRC) is mainly limited to patients whose tumours exhibit high microsatellite instability (MSI). However, there is variability in treatment outcomes within this group, which is in part driven by the frequency and characteristics of tumour infiltrating immune cells. Indeed, the presence of specific infiltrating immune cell subsets has been shown to correlate with immunotherapy responses and is in many cases prognostic of treatment outcome. Tumour-infiltrating lymphocytes (TILs) can undergo distinct differentiation programs, acquiring features of tissue-residency or exhaustion, a process during which T cells upregulate inhibitory receptors such as PD-1 and loose functionality. While residency and exhaustion programs of CD8 T cells are relatively well-studied, these programs have only recently been appreciated in CD4 T cells and remain largely unknown in tumour-infiltrating natural killer (NK) cells. In this study, we use single cell RNA-seq data to identify signatures of residency and exhaustion in CRC infiltrating lymphocytes, including CD8, CD4 and NK cells. We then test these signatures in independent single cell data from tumour and normal tissue infiltrating immune cells. Further, we use versions of these signatures designed for bulk RNA-seq data to explore tumour intrinsic mutations associated with residency and exhaustion from TCGA data. Finally, using two independent transcriptomic data sets from patients with colon adenocarcinoma, we show that combinations of these signatures, in particular combinations of NK activity signatures, together with tumour-associated signatures, such as TGF-{beta} signalling, are associated with distinct survival outcomes in colorectal cancer patients.

category: cancer biology


title: The bacterial virulence factors rhamnolipids and their (R)-3-hydroxyalkanoate precursors activate Arabidopsis innate immunity through two independent mechanisms

abstract:
Plant innate immunity is activated upon perception of invasion pattern molecules by plant cell-surface immune receptors. Several bacteria of the genera Pseudomonas and Burkholderia produce rhamnolipids (RLs) from L- rhamnose and (R)-3-hydroxyalkanoate precursors (HAAs). RL and HAA secretion is required to modulate bacterial surface motility, biofilm development, and thus successful colonization of hosts. Here, we show that the lipidic secretome from the opportunistic pathogen Pseudomonas aeruginosa mostly comprising RLs and HAAs stimulates Arabidopsis immunity. We demonstrate that HAAs are sensed by the bulb-type lectin receptor kinase LIPOOLIGOSACCHARIDE-SPECIFIC REDUCED ELICITATION/S-DOMAIN-1-29 (LORE/SD1-29) that also mediates medium-chain 3-hydroxy fatty acid (mc-3-OH-FA) perception in the plant Arabidopsis thaliana. HAA sensing induces canonical immune signaling and local resistance to plant pathogenic Pseudomonas infection. By contrast, RLs trigger an atypical immune response and resistance to Pseudomonas infection independent of LORE. Thus, the glycosyl moieties of RLs, albeit abolishing sensing by LORE, do not impair their ability to trigger plant defense. In addition, our results show that RL-triggered immune response is affected by the sphingolipid composition of the plasma membrane. In conclusion, RLs and their precursors released by bacteria can both be perceived by plants but through distinct mechanisms.

SignificanceActivation of plant innate immunity relies on the perception of microorganisms through self and nonself elicitors. Rhamnolipids and their precursor HAAs are exoproducts produced by beneficial and pathogenic bacteria. They are involved in bacterial surface dissemination and biofilm development. As these compounds are released in the extracellular milieu, they have the potential to be perceived by the plant immune system. Our work shows that both compounds independently activate plant immunity. We demonstrate that HAAs are perceived by the receptor protein kinase LORE. By contrast, rhamnolipids are not senses by LORE but activate a non-canonical immune response affected by the sphingolipid composition of the plant plasma membrane. Thus, plants are able to sense bacterial molecules as well as their direct precursors to trigger a distinct immune response.

category: plant biology


title: Arabidopsis root responses to salinity depend on pectin modification and cell wall sensing

abstract: Soil salinity is an increasing worldwide problem for agriculture, affecting plant growth and yield. To understand the molecular mechanisms activated in response to salt in plants, we investigated the Catharanthus roseus Receptor like Kinase 1 Like (CrRLK1L) family, which contains well described sensors previously shown to be involved in maintaining and sensing the structural integrity of the cell walls. We found that herk1the1-4 double mutants, lacking the function of the Arabidopsis thaliana Receptor like Kinase HERKULES1 combined with a gain-of-function allele of THESEUS1, phenocopied the phenotypes previously reported in plants lacking FERONIA (FER) function. We report that both fer-4 and herk1the1-4 mutants respond strongly to salt application, resulting in a more intense activation of early and late stress responses. We also show that salt triggers de-methyl esterification of loosely bound pectins, responsible for the activation of several salt response signaling pathways. Addition of calcium chloride or chemically inhibiting pectin methyl esterase (PME) activity reduced activation of the early signaling protein Mitogen Activated Protein Kinase 6 (MPK6) as well as amplitude of salt-induced marker gene induction. MPK6 is required for the full induction of the salt-induced gene expression markers we tested. The sodium specific root halotropism response on the other hand, appears independent of MPK6 or calcium application, and is only mildly influenced by the cell wall sensors FER/HERK1/THE1-4 or alteration of PME activity. We hypothesize a model where salt-triggered modification of pectin requires the functionality of FER alone or the HERK1/THE1 combination to attenuate salt responses. Collectively, our results show the complexity of salt stress responses and salt sensing mechanisms and their connection to cell wall modifications, responsible for several salt response pathways and ultimately plant resilience to salinity.


title: Artificial light at night leads to circadian disruption in a songbird: integrated evidence from behavioural, genomic and metabolomic data

abstract: Globally increasing levels of artificial light at night (ALAN) are associated with shifts in circadian rhythms of behaviour in many wild species. However, it is still unclear whether changes in behavioural timing are underlined by parallel shifts in the molecular clock, and whether such internal shifts may differ between different tissues and physiological pathways, which could highlight circadian disruption. We tackled these questions in a comprehensive study that integrated behavioural, gene expression and metabolomic analyses. We exposed captive male great tits (Parus major) to three ALAN intensities or to dark nights, recorded their activity rhythms and obtained mid-day and midnight samples of brain, liver, spleen and blood. ALAN advanced wake-up time, and this shift was paralleled by an advance in the expression of the clock gene BMAL1 in all tissues, suggesting close links of brain and peripheral clock gene expression with activity rhythms. However, several metabolic and immune genes were desynchronised the shifted BMAL1 expression, suggesting circadian disruption of behaviour and physiology. This result was reinforced by untargeted metabolomic profiling, which showed that only 9.7% of the 755 analysed metabolites followed the behavioural shift. We suggest circadian as a key mediator of the health impacts of ALAN on wild animals.

category: zoology


title: Cerebral μ-opioid and CB1-receptor systems have distinct roles in human feeding behavior

abstract: Eating behavior varies greatly between healthy individuals, but the neurobiological basis of these trait-like differences in feeding remains unknown. Central {micro}-opioid receptors (MOR) and cannabinoid CB1-receptors (CB1R) regulate energy balance via multiple neural pathways, promoting food intake and reward. Because obesity and eating disorders have been associated with alterations in the brains opioid and endocannabinoid signaling, the variation in MOR and CB1R system function could potentially underlie distinct eating behavior phenotypes. In this retrospective positron emission tomography (PET) study, we analyzed [11C]carfentanil PET scans of MORs from 92 healthy subjects (70 males and 22 females), and [18F]FMPEP-d2 scans of CB1Rs from 35 subjects (all males, all also included in the [11C]carfentanil sample). Eating styles were measured with the Dutch Eating Behavior Questionnaire (DEBQ). We found that lower cerebral MOR availability was associated with increased external eating - individuals with low MORs reported being more likely to eat in response to environments palatable food cues. CB1R availability was associated with multiple eating behavior traits. We conclude that although MORs and CB1Rs overlap anatomically and functionally in the brain, they have distinct roles in mediating individual feeding patterns.


title: Age-specific survivorship and fecundity shapegenetic diversity in marine fishes

abstract:
Genetic diversity varies among species due to a range of eco-evolutionary processes that are not fully understood. The neutral theory predicts that the amount of variation in the genome sequence between different individuals of the same species should increase with its effective population size (Ne). In real populations, multiple factors that modulate the variance in reproductive success among individuals cause Ne to differ from the total number of individuals (N). Among these, age-specific mortality and fecundity rates are known to have a direct impact on the [Formula] ratio. However, the extent to which vital rates account for differences in genetic diversity among species remains unknown. Here, we addressed this question by comparing genome-wide genetic diversity across 16 marine fish species with similar geographic distributions but contrasted lifespan and age-specific survivorship and fecundity curves. We sequenced the whole genome of 300 individuals to high coverage and assessed their genome-wide heterozygosity with a reference-free approach. Genetic diversity varied from 0.2 to 1.4% among species, and showed a negative correlation with adult lifespan, with a large negative effect (slope = - 0.089 per additional year of lifespan) that was further increased when brooding species providing intense parental care were removed from the dataset (slope = -0.129 per additional year of lifespan). Using published vital rates for each species, we showed that the [Formula] ratio resulting simply from life tables parameters can predict the observed differences in genetic diversity among species. Using simulations, we further found that the extent of reduction in [Formula] with increasing adult lifespan is particularly strong under Type III survivorship curves (high juvenile and low adult mortality) and increasing fecundity with age, a typical characteristic of marine fishes. Our study highlights the importance of vital rates as key determinants of species genetic diversity levels in nature.

Author SummaryUnderstanding how and why genetic diversity varies across species has important implications for evolutionary and conservation biology. Although genomics has vastly improved our ability to document intraspecific DNA sequence variation at the genome level, the range and determinants of genetic diversity remain partially understood. At a broad taxonomic scale in eukaryotes, the main determinants of diversity are reproductive strategies distributed along a trade-off between the quantity and the size of offspring, which likely affect the long-term effective population size. Long-lived species also tend to show lower genetic diversity, a result which has however not been reported by comparative studies of genetic diversity at lower taxonomic scales. Here, we compared genetic diversity across 16 European marine fish species showing marked differences in longevity. Adult lifespan was the best predictor of genetic diversity, with genome-wide average heterozygosity ranging from 0.2% in the black anglerfish (L. budegassa) to 1.4% in the European pilchard (S. pilchardus). Using life tables summarizing age-specific mortality and fecundity rates for each species, we showed that the variance in lifetime reproductive success resulting from age structure, iteroparity and overlapping generations can predict the range of observed differences in genetic diversity among marine fish species. We then used computer simulations to explore how combinations of vital rates characterizing different life histories affect the relationship between adult lifespan and genetic diversity. We found that marine fishes that display high juvenile but low adult mortality, and increasing fecundity with age, are typically expected to show reduced genetic diversity with increased adult lifespan. However, the impact of adult lifespan vanished using bird and mammal-like vital rates. Our study shows that variance in lifetime reproductive success can have a major impact on species genetic diversity and explains why this effect varies widely across taxonomic groups.


title: Standardizing a protocol for studying pH-dependent transient conformations using computer simulations

abstract: 1.Understanding pH-dependent protein stability is important in biological - transport, storage, and delivery, in vivo conditions such as aggregation potential in neurodegenerative disease, and in studying the folding/unfolding of proteins. Using computer simulations, we can replace complex experimental determination and provide an atomistic-level interpretation of the cause and effect of pH on protein stability. Here, we standardize a method that provides a framework through which we examined pH-dependent transient conformations during unfolding simulations of proteins. Constant pH simulations utilized in the prediction of pKa values of charged groups of the peptide. The calculated pKa values employed to fix the appropriate protonation state of the amino acid to simulate the effect of pH on the system. Trajectories from multiple high-temperature MD simulations of the protein sample the conformational space during unfolding for a given pH state. The ensemble of conformations is analyzed from its free energy landscape to identify transient and stable conformations both at a given pH and between different pH. As a test system RN80, a protein fragment analog of the C-peptide from bovine pancreatic ribonuclease-A used to measure the accuracy of the predictions from simulations. Experimental measures of the helix content determined as a function of pH display a bell-shaped curve, i.e. RN80 alpha-helix formation is maximum at pH5 with a subsequent loss in helicity at higher and lower pH. The main forces stabilizing the alpha-helix are a salt-bridge formed between Glu-2 and Arg-10 and cation-pi-interaction between Tyr-8 and His-12. Our protocol includes constant pH calculations, optimal high-temperature simulations, and Free Energy landscape analysis exhibited the agreement with the experimental observations.

category: bioinformatics


title: Post-transcriptional gene regulation by the RNA binding protein IGF2BP3 is critical for MLL-AF4 mediated leukemogenesis

abstract: Despite recent advances in therapeutic approaches, patients with MLL-rearranged leukemia still have poor outcomes and a high risk of relapse. Here, we found that MLL-AF4, the most common MLL fusion protein in patients, transcriptionally induces IGF2BP3 and that IGF2BP3 strongly amplifies MLL-Af4 mediated leukemogenesis. Deletion of Igf2bp3 significantly increases the survival of mice with MLL-Af4 driven leukemia and greatly attenuates disease, with a minimal impact on baseline hematopoiesis. At the cellular level, MLL-Af4 leukemia-initiating cells require Igf2bp3 for their function in leukemogenesis. eCLIP and transcriptome analysis of MLL-Af4 transformed stem and progenitor cells and MLL-Af4 bulk leukemia cells reveals a complex IGF2BP3-regulated post-transcriptional operon governing leukemia cell survival and proliferation. Regulated mRNA targets include important leukemogenic genes such as those in the Hoxa locus and numerous genes within the Ras signaling pathway. Together, our findings show that IGF2BP3 is an essential positive regulator of MLL-AF4 mediated leukemogenesis and represents an attractive therapeutic target in this disease.


title: Inferring brain-wide interactions using data-constrained recurrent neural network models

abstract: Behavior arises from the coordinated activity of numerous anatomically and functionally distinct brain regions. Modern experimental tools allow unprecedented access to large neural populations spanning many interacting regions brain-wide. Yet, understanding such large-scale datasets necessitates both scalable computational models to extract meaningful features of inter-region communication and principled theories to interpret those features. Here, we introduce Current-Based Decomposition (CURBD), an approach for inferring brain-wide interactions using data-constrained recurrent neural network models that directly reproduce experimentally-obtained neural data. CURBD leverages the functional interactions inferred by such models to reveal directional currents between multiple brain regions. We first show that CURBD accurately isolates inter-region currents in simulated networks with known dynamics. We then apply CURBD to multi-region neural recordings obtained from mice during running, macaques during Pavlovian conditioning, and humans during memory retrieval to demonstrate the widespread applicability of CURBD to untangle brain-wide interactions underlying behavior from a variety of neural datasets.


title: Density and temperature controlled fluid extraction in a bacterial biofilm is determined by poly-γ-glutamic acid production

abstract: A hallmark of microbial biofilms is the self-production of extracellular matrix that encases the cells resident within the community. The matrix provides protection from the environment, while spatial heterogeneity of expression influences the structural morphology and colony spreading dynamics. Bacillus subtilis is a model bacterial system used to uncover the regulatory pathways and key building blocks required for biofilm growth and development. Previous reports have suggested that poly-{gamma}-glutamic acid (PGA) production is suppressed during biofilm formation and does not play a major role in biofilm morphology of the undomesticated isolate NCIB 3610. In this work we report on the observation of multiple travelling fronts that develop during the early stage of B. subtilis colony biofilm formation. We find the emergence of a highly motile population of bacteria that is facilitated by the extraction of fluid from the underlying agar substrate. Motility develops behind a moving front of fluid that propagates from the boundary of the biofilm towards the interior. The extent of proliferation is strongly modulated by the presence of extracellular polysaccharides (EPS). We trace the origin of this moving front of fluid to the production of PGA. We find that PGA production is correlated with higher temperatures, resulting in a mature biofilm morphology that is distinct from the biofilm architecture typically associated with B. subtilis. Our results suggest that B. subtilis NCIB 3610 produces distinct biofilm matrices in response to environmental conditions.


title: Truncated denitrifiers dominate the denitrification pathway in tundra soil metagenomes

abstract:
BackgroundIn contrast to earlier assumptions, there is now mounting evidence for the role of tundra soils as important sources of the greenhouse gas nitrous oxide (N2O). However, the microorganisms involved in the cycling of N2O in this system remain largely uncharacterized. Since tundra soils are variable sources and sinks of N2O, we aimed at investigating differences in community structure across different soil ecosystems in the tundra.

ResultsWe analysed 1.4 Tb of metagenomic data from soils in northern Finland covering a range of ecosystems from dry upland soils to water-logged fens and obtained 796 manually binned and curated metagenome-assembled genomes (MAGs). We then searched for MAGs harbouring genes involved in denitrification, an important process driving N2O emissions. Communities of potential denitrifiers were dominated by microorganisms with truncated denitrification pathways (i.e., lacking one or more denitrification genes) and differed across soil ecosystems. Upland soils showed a strong N2O sink potential and were dominated by members of the Alphaproteobacteria such as Bradyrhizobium and Reyranella. Fens, which had in general net-zero N2O fluxes, had a high abundance of poorly characterized taxa affiliated with the Chloroflexota lineage Ellin6529 and the Acidobacteriota subdivision Gp23.

ConclusionsBy coupling an in-depth characterization of microbial communities with in situ measurements of N2O fluxes, our results suggest that the observed spatial patterns of N2O fluxes in the tundra are related to differences in the composition of denitrifier communities.


title: In-depth characterization of denitrifier communities across different soil ecosystems in the tundra


title: Evidence of a dysregulated Vitamin D pathway in SARS-CoV-2 infected patient's lung cells

abstract: Although a defective vitamin D pathway has been widely suspected to be associated in SARS-CoV-2 pathobiology, the status of the vitamin D pathway and vitamin D-modulated genes in lung cells of patients infected with SARS-CoV-2 remains unknown. To understand the significance of the vitamin D pathway in SARS-CoV-2 pathobiology, computational approaches were applied to transcriptomic datasets from bronchoalveolar lavage fluid (BALF) cells of such patients or healthy individuals. Levels of vitamin D receptor, retinoid X receptor, and CYP27A1 in BALF cells of patients infected with SARS-CoV-2 were found to be reduced. Additionally, 107 differentially expressed, predominantly downregulated genes modulated by vitamin D were identified in transcriptomic datasets from patients cells. Further analysis of differentially expressed genes provided eight novel genes with a conserved motif with vitamin D-responsive elements, implying the role of both direct and indirect mechanisms of gene expression by the dysregulated vitamin D pathway in SARS-CoV-2-infected cells. Network analysis of differentially expressed vitamin D-modulated genes identified pathways in the immune system, NF-KB/cytokine signaling, and cell cycle regulation as top predicted pathways that might be affected in the cells of such patients. In brief, the results provided computational evidence to implicate a dysregulated vitamin D pathway in the pathobiology of SARS-CoV-2 infection.


title: From infection to immunity: understanding the response to SARS-CoV2 through in-silico modeling

abstract:
BackgroundImmune system conditions of the patient is a key factor in COVID-19 infection survival. A growing number of studies have focused on immunological determinants to develop better biomarkers for therapies.

AimThe dynamics of the insurgence of immunity is at the core of the both SARS-CoV-2 vaccine development and therapies. This paper addresses a fundamental question in the management of the infection: can we describe the insurgence (and the span) of immunity in COVID-19? The in-silico model developed here answers this question at individual (personalized) and population levels.

We simulate the immune response to SARS-CoV-2 and analyze the impact of infecting viral load, affinity to the ACE2 receptor and age in the artificially infected population on the course of the disease.

MethodsWe use a stochastic agent-based immune simulation platform to construct a virtual cohort of infected individuals with age-dependent varying degree of immune competence. We use a parameter setting to reproduce known inter-patient variability and general epidemiological statistics.

ResultsWe reproduce in-silico a number of clinical observations and we identify critical factors in the statistical evolution of the infection. In particular we evidence the importance of the humoral response over the cytotoxic response and find that the antibody titers measured after day 25 from the infection is a prognostic factor for determining the clinical outcome of the infection.

Our modeling framework uses COVID-19 infection to demonstrate the actionable effectiveness of simulating the immune response at individual and population levels. The model developed is able to explain and interpret observed patterns of infection and makes verifiable temporal predictions.

Within the limitations imposed by the simulated environment, this work proposes in a quantitative way that the great variability observed in the patient outcomes in real life can be the mere result of subtle variability in the infecting viral load and immune competence in the population.

In this work we i) show the power of model predictions, ii) identify the clinical end points that could be more suitable for computational modeling of COVID-19 immune response, iii) define the resolution and amount of data required to empower this class of models for translational medicine purposes and, iv) we exemplify how computational modeling of immune response provides an important view to discuss hypothesis and design new experiments, in particular paving the way to further investigations about the duration of vaccine-elicited immunity especially in the view of the blundering effect of immunesenescence.


title: Enhanced reproductive thermotolerance is associated with increased accumulation of flavonols in pollen of the tomato high-pigment 2 mutant.

abstract: Climate change has created an environment where heat stress conditions are becoming more frequent as temperatures continue to rise in crop production areas around the world. This situation leads to decreased crop production due to plant sensitivity to heat stress. Reproductive success is critically dependent on plants ability to produce functional pollen grains, which are the most thermo-sensitive tissue. Flavonols are plant secondary metabolites known for their potent antioxidative activity, essential for male fertility in several species including tomato, and implicated in heat stress tolerance. Since flavonols are highly abundant in fruits of the tomato high pigment-2 (hp2) mutant, we tested the level of flavonols in pollen of this mutant, under the hypothesis that increased accumulation of flavonols would render pollen more tolerant to heat stress. Indeed, pollen from three alleles of the hp2 mutant were found to have flavonols levels increased by 40% on average compared with wild-type under moderate chronic heat stress conditions. This mutant produced on average 7.8-fold higher levels of viable pollen and displayed better germination competence under heat stress conditions. The percentage of fully seeded fruits and the number of seeds per fruit were maintained in the mutant under heat stress conditions while decreased in wild-type plants. Our results strongly suggest that increased pollen flavonols content enhances pollen thermotolerance and reproductive success under heat stress conditions. Thus, the high flavonols trait may help frame the model for improving crop resilience to heat stress.


title: Test-retest reliability of neural entrainment in the human auditory system

abstract: Auditory stimuli are often rhythmic in nature. Brain activity synchronizes with auditory rhythms via neural entrainment, and entrainment seems to be beneficial for auditory perception. However, it is not clear to what extent neural entrainment in the auditory system is reliable over time - a necessary prerequisite for targeted intervention. The current study aimed to establish the reliability of neural entrainment over time and to predict individual differences in auditory perception from associated neural activity. Across two different sessions, human listeners detected silent gaps presented at different phase locations of a 2-Hz frequency modulated (FM) noise while EEG activity was recorded. As expected, neural activity was entrained by the 2-Hz FM noise. Moreover, gap detection was sinusoidally modulated by the phase of the 2-Hz FM into which the gap fell. Critically, both the strength of neural entrainment as well as the modulation of performance by the stimulus rhythm were highly reliable over sessions. Moreover, gap detection was predictable from pre-gap neural 2-Hz phase and alpha amplitude. Our results demonstrate that neural entrainment in the auditory system and the resulting behavioral modulation are reliable over time, and that both entrained delta and non-entrained alpha oscillatory activity contribute to near-threshold stimulus perception. The latter suggests that improving auditory perception might require simultaneously targeting entrained brain rhythms as well as the alpha rhythm.


title: The SARS-CoV-2 Spike protein disrupts human cardiac pericytes function through CD147-receptor-mediated signalling: a potential non-infective mechanism of COVID-19 microvascular disease

abstract:
Severe coronavirus disease 2019 (COVID-19) manifests as a life-threatening microvascular syndrome. The severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) uses the Spike (S) protein to engage with its receptors and infect host cells. To date, it is still not known whether heart vascular pericytes (PCs) are infected by SARS-CoV-2, and if the S protein alone provokes PC dysfunction. Here, we aimed to investigate the effects of the S protein on primary human cardiac PC signalling and function. Results show, for the first time, that cardiac PCs are not permissive to SARS-CoV-2 infection in vitro, whilst a recombinant S protein alone elicits functional alterations in PCs. This was documented as: (1) increased migration, (2) reduced ability to support endothelial cell (EC) network formation on Matrigel, (3) secretion of pro-inflammatory molecules typically involved in the cytokine storm, and (4) production of pro-apoptotic factors responsible for EC death. Next, adopting a blocking strategy against the S protein receptors angiotensin-converting enzyme 2 (ACE2) and CD147, we discovered that the S protein stimulates the phosphorylation/activation of the extracellular signal-regulated kinase 1/2 (ERK1/2) through the CD147 receptor, but not ACE2, in PCs. The neutralisation of CD147, either using a blocking antibody or mRNA silencing, reduced ERK1/2 activation and rescued PC function in the presence of the S protein. In conclusion, our findings suggest that circulating S protein prompts vascular PC dysfunction, potentially contributing to establishing microvascular injury in organs distant from the site of infection. This mechanism may have clinical and therapeutic implications.

Clinical perspectiveO_LISevere COVID-19 manifests as a microvascular syndrome, but whether SARS-CoV-2 infects and damages heart vascular pericytes (PCs) remains unknown.
C_LIO_LIWe provide evidence that cardiac PCs are not infected by SARS-CoV-2. Importantly, we show that the recombinant S protein alone elicits cellular signalling through the CD147 receptor in cardiac PCs, thereby inducing cell dysfunction and microvascular disruption in vitro.
C_LIO_LIThis study suggests that soluble S protein can potentially propagate damage to organs distant from sites of infection, promoting microvascular injury. Blocking the CD147 receptor in patients may help protect the vasculature not only from infection, but also from the collateral damage caused by the S protein.
C_LI


title: Soil microbial communities are not altered by garlic mustard in recently invaded central Illinois forests

abstract: The invasive forest plant garlic mustard (Alliaria petiolata) has been shown to alter soil microbial communities in the northeastern part of its invaded range in the United States, and this disruption of soil communities may contribute to its invasion success. However, garlic mustard allelochemistry can vary with invasion age, and it is not clear whether garlic mustards impacts on soil microbes are consistent over its invaded range. Here, we compare the composition and diversity of soil fungal, bacterial, and archaeal communities among garlic mustard present, absent, and removed treatments in replicated blocks across five forests in the midwestern United States with relatively young garlic mustard invasions (approximately 17-26 years old, with consistent management). We collected samples in May and August, corresponding to garlic mustard active and senescent life history stages. While soil fungal and bacterial/ archaeal communities (based on ITS2 region and 16S rRNA gene DNA sequencing, respectively) differed significantly between different blocks/ forests and over time, we found no significant effect of garlic mustard treatment on soil microbial community composition or the relative abundance of mycorrhizal, saprotrophic, or pathogenic fungal guilds. The lack of garlic mustard impacts on the soil microbial community in recently invaded central Illinois forests suggests that these well-documented impacts in the northeastern United States and in older invasions cannot necessarily be generalized across all environmental contexts.

category: ecology


title: The molecular framework of heterophylly in Callitriche palustris L. differs from that in other amphibious plants

abstract: Heterophylly refers to the development of different leaf forms in a single plant depending on the environmental conditions. It is often observed in amphibious aquatic plants that can grow under aerial and submerged conditions. Although heterophylly is well recognized in aquatic plants, the associated developmental mechanisms and the molecular basis remain unclear. In this study, we analyzed heterophyllous leaf formation in an aquatic plant, Callitriche palustris, to clarify the underlying developmental and molecular mechanisms. Morphological analyses revealed extensive cell elongation and the rearrangement of cortical microtubules in the elongated submerged leaves of C. palustris. Our observations also suggested that gibberellin, ethylene, and abscisic acid regulate the formation of submerged leaves. However, the perturbation of one or more of the hormones was insufficient to induce the formation of submerged leaves under aerial conditions. Finally, we analyzed gene expression changes during aerial and submerged leaf development and narrowed down the candidate genes controlling heterophylly via transcriptomic comparisons, including a comparison with a closely related terrestrial species. We revealed that the molecular mechanism regulating heterophylly in C. palustris is associated with hormonal changes and diverse transcription factor gene expression profiles, suggesting differences from the corresponding mechanisms in previously investigated amphibious plants.


title: Automation of Leaf Counting in Maize and Sorghum Using Deep Learning

abstract: Leaf number and leaf emergence rate are phenotypes of interest to plant breeders, plant geneticists, and crop modelers. Counting the extant leaves of an individual plant is straightforward even for an untrained individual, but manually tracking changes in leaf numbers for hundreds of individuals across multiple time points is logistically challenging. This study generated a dataset including over 150,000 maize and sorghum images for leaf counting projects. A subset of 17,783 images also includes annotations of the positions of individual leaf tips. With these annotated images, we evaluate two deep learning-based approaches for automated leaf counting: the first based on counting-by-regression from whole image analysis and a second based on counting-by-detection. Both approaches can achieve RMSE (root of mean square error) smaller than one leaf, only moderately inferior to the RMSE between human annotators of between 0.57 and 0.73 leaves. The counting-by-regression approach based on CNNs (convolutional neural networks) exhibited lower accuracy and increased bias for plants with extreme leaf numbers which are underrepresented in this dataset. The counting-by-detection approach based on Faster R-CNN object detection models achieve near human performance for plants where all leaf tips are visible. The annotated image data and model performance metrics generated as part of this study provide large scale resources for the comparison and improvement of algorithms for leaf counting from image data in grain crops.


title: Genomic analysis provides novel insights into diversification and taxonomy of Allorhizobium vitis (i.e. Agrobacterium vitis)

abstract:
BackgroundAllorhizobium vitis (formerly named Agrobacterium vitis or Agrobacterium biovar 3) is the primary causative agent of crown gall disease of grapevine worldwide. We obtained and analyzed whole-genome sequences of diverse All. vitis strains to get insights into their diversification and taxonomy.

ResultsPairwise genome comparisons and phylogenomic analysis of various All. vitis strains clearly indicated that All. vitis is not a single species, but represents a species complex composed of several genomic species. Thus, we emended the description of All. vitis, which now refers to a restricted group of strains within the All. vitis species complex (i.e. All. vitis sensu stricto) and proposed a description of a novel species, All. ampelinum sp. nov. The type strain of All. vitis sensu stricto remains the current type strain of All. vitis, K309T. The type strain of All. ampelinum sp. nov. is S4T. We also identified sets of gene clusters specific to the All. vitis species complex, All. vitis sensu stricto and All. ampelinum, respectively, for which we predicted the biological function and infer the role in ecological diversification of these clades, including some we could experimentally validate. All. vitis species complex-specific genes confer tolerance to different stresses, including exposure to aromatic compounds. Similarly, All. vitis sensu stricto-specific genes confer the ability to degrade 4-hydroxyphenylacetate and a putative compound related to gentisic acid. All. ampelinum-specific genes have putative functions related to polyamine metabolism and nickel assimilation. Congruently with the genome-based classification, All. vitis sensu stricto and All. ampelinum were clearly delineated by MALDI-TOF MS analysis. Moreover, our genome-based analysis indicated that Allorhizobium is clearly separated from other genera of the family Rhizobiaceae.

ConclusionsComparative genomics and phylogenomic analysis provided novel insights into the diversification and taxonomy of Allorhizobium vitis species complex, supporting our redefinition of All. vitis sensu stricto and description of All. ampelinum. Our pan-genome analyses suggest that these species have differentiated ecologies, each relying on specialized nutrient consumption or toxic compound degradation to adapt to their respective niche.


title: A global map of the impact of deletion of Post-Translational Modification sites in genetic diseases

abstract:
BackgroundThere are >200 protein post-translational modification (PTMs) types described in eukaryotes, having diverse species conservation levels, proteome coverage, number of high-throughput experiments and functional roles. From a clinical perspective, a number of diseases have been associated to deregulated PTM sites and missense rare variants are globally enriched in PTMs. We hypothesize that some genetic diseases may be caused by the deregulation of particular functions produced by the removal of a specific PTM type by genomic variants.

ResultsWe collected >320,000 human PTMs of 59 types and cross them with >4M missense DNA variants annotated with pathogenic predictions and disease associations. We report >1.74M PTM-variant concurrences in >16,500 proteins that an enrichment analysis distributed in 217 pairwise significant associations between 18 PTM types and 150 genetic diseases. Around 23% of these associations are already described in the literature, 34% have partial evidences based on single variants, related diseases or regulatory evidences, and 43% are novel. Removal of acetylation presents the highest effect, still low studied PTM types like S-glutathionylation or S-nitrosylation show relevance. A network of PTM types and phenotypes associations is also discussed. Using pathogenicity predictions we identified potential PTM sites to produce particular diseases if genomic variants remove them.

ConclusionsOur results show an important impact of PTM removal producing genetic diseases and phenotypes that is PTM type specific. We describe for the first time a general scenario of PTM types and genetic diseases direct associations, many of them novel, that provides new capacities to understand and diagnose these disorders.

category: systems biology


title: Progressing adaptation of SARS-CoV-2 to humans

abstract: The second and subsequent waves of coronavirus disease 2019 (COVID-19) have caused problems worldwide 1. Here, using an objective analytical method 2, we present the changes that occurred in the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), the causative virus of COVID-19, over time. The virus has mutated in three major directions, resulting in three groups to date. Analysis of the basic structure of the group of viruses was completed by April and shared across all continents. However, the virus continued to mutate independently in each country after the borders were closed. In particular, the virus mutated before the occurrence of the second and subsequent peaks. It seems that the mutations conferred higher infectivity to the virus, because of which the virus overcame previously effective protection and caused second waves of the disease. Currently, each country may possess such a unique, stronger variant. Some of them slowly entered other countries and caused epidemics. These viruses could also serve as sources of further mutations by exchanging parts of the genome, which could create variants with superior infectivity.


title: Cetylpyridinium chloride-containing mouthwashes reduce the infectivity of SARS-CoV-2 variants in vitro

abstract: Oral mouthwashes decrease the infectivity of several respiratory viruses including SARS-CoV-2. However, the precise agents with antiviral activity present in these oral rinses and their exact mechanism of action remain unknown. Here we show that Cetylpyridinium chloride (CPC), a quaternary ammonium compound present in many oral mouthwashes, reduces SARS-CoV-2 infectivity by inhibiting the viral fusion step with target cells after disrupting the integrity of the viral envelope. We also found that CPC-containing mouth rinses decreased more than a thousand times the infectivity of SARS-CoV-2 in vitro, while the corresponding vehicles had no effect. This activity was effective for different SARS-CoV-2 variants, including the B.1.1.7 variant, predominant in UK, also in the presence of sterilized saliva. CPC-containing mouth rinses could therefore represent a cost-effective measure to reduce SARS-CoV-2 infectivity in saliva, aiding to reduce viral transmission from infected individuals regardless of the variants they are infected with.


title: A De Novo regulation design shows an effectiveness in altering plant secondary metabolism

abstract:
IntroductionTranscription factors (TFs) and cis-regulatory elements (CREs) control gene transcripts involved in various biological processes. We hypothesize that TFs and CREs can be effective molecular tools for De Novo regulation designs to engineer plants.

ObjectivesWe selected two Arabidopsis TF types and two tobacco CRE types to design a De Novo regulation and evaluated its effectiveness in plant engineering.

MethodsG-box and MYB recognition elements (MREs) were identified in four Nicotiana tabacum JAZs (NtJAZs) promoters. MRE-like and G-box like elements were identified in one nicotine pathway gene promoter. TF screening led to select Arabidopsis Production of Anthocyanin Pigment 1 (PAP1/MYB) and Transparent Testa 8 (TT8/bHLH). Two NtJAZ and two nicotine pathway gene promoters were cloned from commercial Narrow Leaf Madole (NL) and KY171 (KY) tobacco cultivars. Electrophoretic mobility shift assay (EMSA), cross-linked chromatin immunoprecipitation (ChIP), and dual luciferase assays were performed to test the promoter binding and activation by PAP1 (P), TT8 (T), PAP1/TT8 together, and the PAP1/TT8/Transparent Testa Glabra 1 (TTG1) complex. A DNA cassette was designed and then synthesized for stacking and expressing PAP1 and TT8 together. Three years of field trials were performed by following industrial and GMO protocols. Gene expression and metabolic profiling were completed to characterize plant secondary metabolism.

ResultsPAP1, TT8, PAP1/TT8, and the PAP1/TT8/TTG1 complex bound to and activated NtJAZ promoters but did not bind to nicotine pathway gene promoters. The engineered red P+T plants significantly upregulated four NtJAZs but downregulated the tobacco alkaloid biosynthesis. Field trials showed significant reduction of five tobacco alkaloids and four carcinogenic tobacco specific nitrosamines in most or all cured leaves of engineered P+T and PAP1 genotypes.

ConclusionG-boxes, MREs, and two TF types are appropriate molecular tools for a De Novo regulation design to create a novel distant-pathway cross regulation for altering plant secondary metabolism.


title: RedOak: a reference-free and alignment-freestructure for indexing a collection of similargenomes

abstract:
BackgroundAs the cost of DNA sequencing decreases, high-throughput sequencing technologies become increasingly accessible to many laboratories. Consequently, new issues emerge that require new algorithms, including tools for indexing and compressing hundred to thousands of complete genomes.

ResultsThis paper presents RedOak, a reference-free and alignment-free software package that allows for the indexing of a large collection of similar genomes. RedOak can also be applied to reads from unassembled genomes, and it provides a nucleotide sequence query function. This software is based on a k-mer approach and has been developed to be heavily parallelized and distributed on several nodes of a cluster. The source code of our RedOak algorithm is available at https://gite.lirmm.fr/doccy/RedOak.

ConclusionsRedOak may be really useful for biologists and bioinformaticians expecting to extract information from large sequence datasets.


title: Robust data storage in DNA by de Bruijn graph-based decoding

abstract:
DNA data storage is a rapidly developing technology with great potential due to its high density, long-term durability, and low maintenance cost. The major technical challenges include various errors, such as the strand breaks, rearrangements, and indels that frequently arise during DNA synthesis, amplification, sequencing, and preservation. Through a de novo assembly strategy, we developed an algorithm based on the de Bruijn graph and greedy path search (DBGPS) to address these issues. DBGPS shows distinct advantages in handling DNA breaks, rearrangements, and indels. The robustness of DBGPS is demonstrated by accelerated aging, multiple independent data retrievals, deep error-prone PCR, and large data scale simulations. Remarkably, 6.8 MB of data can be retrieved accurately from a seriously corrupted sample that has been treated at 70 {degrees}C for 70 days. With DBGPS, we were able to achieve a logical density of 1.30 bits/cycle and a physical density of 295 PB/g.

One-Sentence SummaryA de Bruijn graph-based de novo assembly algorithm for DNA data storage enables fast and robust data readouts even with DNA samples that have been severely corrupted.


category: synthetic biology


title: Robust retrieval of data stored in DNA by de Bruijn graph-based de novo strand assembly


title: Global community effect: large-scale cooperation yields collective survival of differentiating embryonic stem cells

abstract: Cells can help each other replicate by communicating with diffusible molecules. In cell cultures, molecules may diffuse within a cell colony or between adjacent or distant colonies. Determining which cell helps which cells replication is challenging. We developed a systematic approach, integrating modeling and experiments, for determining the length-scales of cell-cell communication (from microns to centimeters). With this approach, we discovered that differentiating murine ES cells, scattered across centimeters on a dish, communicate over millimeters to form one macroscopic entity that survives if and only if its centimeter-scale population-density is above a threshold value. Single-cell-level measurements, transcriptomics, and modeling revealed that this "macroscopic quorum sensing" arises from differentiating ES cells secreting and sensing survival-promoting FGF4 that diffuses over millimeters and activates YAP1-induced survival mechanisms. Through the same mechanism, a lone macroscopic, but not microscopic, colony survives differentiation. Our work rigorously establishes that in vitro ES-cell differentiation relies on macroscopic cooperation.


title: Centimeter-scale quorum sensing dictates collective survival of differentiating embryonic stem cells


title: Deep and accurate detection of m6A RNA modifications using miCLIP2 and m6Aboost machine learning

abstract:
N6-methyladenosine (m6A) is the most abundant internal RNA modification in eukaryotic mRNAs and influences many aspects of RNA processing. miCLIP (m6A individual-nucleotide resolution UV crosslinking and immunoprecipitation) is an antibody-based approach to map m6A sites with single-nucleotide resolution. However, due to broad antibody reactivity, reliable identification of m6A sites from miCLIP data remains challenging. Here, we present miCLIP2 in combination with machine learning to significantly improve m6A detection. The optimised miCLIP2 results in high-complexity libraries from less input material. Importantly, we established a robust computational pipeline to tackle the inherent issue of false positives in antibody-based m6A detection. The analyses are calibrated with Mettl3 knockout cells to learn the characteristics of m6A deposition, including m6A sites outside of DRACH motifs. To make our results universally applicable, we trained a machine learning model, m6Aboost, based on the experimental and RNA sequence features. Importantly, m6Aboost allows prediction of genuine m6A sites in miCLIP2 data without filtering for DRACH motifs or the need for Mettl3 depletion. Using m6Aboost, we identify thousands of high-confidence m6A sites in different murine and human cell lines, which provide a rich resource for future analysis. Collectively, our combined experimental and computational methodology greatly improves m6A identification.

HighlightsO_LImiCLIP2 produces complex libraries to map m6A RNA modifications
C_LIO_LIMettl3 KO miCLIP2 allows to identify Mettl3-dependent RNA modification sites
C_LIO_LIMachine learning predicts genuine m6A sites from human and mouse miCLIP2 data without Mettl3 KO
C_LIO_LIm6A modifications occur outside of DRACH motifs and associate with alternative splicing
C_LI


title: The ENIGMA Toolbox: Cross-disorder integration and multiscale neural contextualization of multisite neuroimaging datasets

abstract:
Among  big data initiatives, the ENIGMA (Enhancing NeuroImaging Genetics through Meta-Analysis) Consortium--a worldwide alliance of over 2,000 scientists diversified into over 50 Working Groups--has yielded some of the largest studies of the healthy and diseased brain. Integration of multisite datasets to assess transdiagnostic similarities and differences and to contextualize findings with respect to neural organization, however, have been limited. Here, we introduce the ENIGMA Toolbox, a Python/Matlab ecosystem for (i) accessing 100+ ENIGMA datasets, facilitating cross-disorder analysis, (ii) visualizing data on brain surfaces, and (iii) contextualizing findings at the microscale (postmortem cytoarchitecture and gene expression) and macroscale (structural and functional connectomes). Our Toolbox equips scientists with tutorials to explore molecular, histological, and network correlates of noninvasive neuroimaging markers of brain disorders. Moreover, our Toolbox bridges the gap between standardized data processing protocols and analytic workflows and facilitates cross-consortia initiatives. The Toolbox is documented and openly available at http://enigma-toolbox.readthedocs.io.

O_FIG O_LINKSMALLFIG WIDTH=200 HEIGHT=63 SRC="FIGDIR/small/423838v3_ufig1.gif" ALT="Figure 1">
View larger version (16K):
org.highwire.dtl.DTLVardef@183a3f3org.highwire.dtl.DTLVardef@180d368org.highwire.dtl.DTLVardef@b752a3org.highwire.dtl.DTLVardef@1928c53_HPS_FORMAT_FIGEXP  M_FIG C_FIG


title: Unbiased population heterozygosity estimates from genome-wide sequence data

abstract:
O_LIHeterozygosity is a metric of genetic variability frequently used to inform the management of threatened taxa. Estimating observed and expected heterozygosities from genome-wide sequence data has become increasingly common, and these estimates are often derived directly from genotypes at single nucleotide polymorphism (SNP) markers. While many SNP markers can provide precise estimates of genetic processes, the results of  downstream analysis with these markers may depend heavily on  upstream filtering decisions.
C_LIO_LIHere we explore the downstream consequences of sample size, rare allele filtering, missing data thresholds and known population structure on estimates of observed and expected heterozygosity using two reduced-representation sequencing datasets, one from the mosquito Aedes aegypti (ddRADseq) and the other from a threatened grasshopper, Keyacris scurra (DArTseq).
C_LIO_LIWe show that estimates based on polymorphic markers only (i.e. SNP heterozygosity) are always biased by global sample size (N), with smaller N producing larger estimates. By contrast, results are unbiased by sample size when calculations consider monomorphic as well as polymorphic sequence information (i.e. genome-wide or autosomal heterozygosity). SNP heterozygosity is also biased when differentiated populations are analysed together, while autosomal heterozygosity remains unbiased. We also show that when nucleotide sites with missing genotypes are included, observed and expected heterozygosity estimates diverge in proportion to the amount of missing data permitted at each site.
C_LIO_LIWe make three recommendations for estimating genome-wide heterozygosity: (i) autosomal heterozygosity should be reported instead of (or in addition to) SNP heterozygosity; (ii) sites with any missing data should be omitted; (iii) populations should be analysed in independent runs. This should facilitate comparisons within and across studies and between observed and expected measures of heterozygosity.
C_LI


title: Polaratio: A magnitude-contingent monotonic correlation metric and its improvements to scRNA-seq clustering

abstract:
MotivationSingle-cell RNA sequencing (scRNA-seq) technologies and analysis tools have allowed researchers to achieve remarkably detailed understandings of the roles and relationships between cells and genes. However, conventional distance metrics, such as Euclidean, Pearson, and Spearman distances, fail to simultaneously take into account the high dimensionality, monotonicity, and magnitude of gene expression data. To address several shortcomings in these commonly used metrics, we present a magnitude-contingent monotonic correlation metric called Polaratio which is designed to enhance the quality of scRNA-seq data analysis.

ResultsWe integrate three interpretable clustering algorithms - Single-Cell Consensus Clustering (SC3), Hierarchical Clustering (HC), and K-Medoids (KM) - through a consensus cell clustering procedure, which we evaluate on various biological datasets to benchmark Polaratio against several well-known metrics. Our results demonstrate Polaratios ability to improve the accuracy of cell clustering on 5 out of 7 publicly available datasets.

Availabilityhttps://github.com/dubai03nsr/Polaratio

Contactpcicalese@uh.edu


title: A unified evolutionary origin for the ubiquitous protein transporters SecY and YidC

abstract: Cells use transporters to move protein across membranes, but the origins of the most ancient transporters are unknown. Here, we analyse the ubiquitous protein-conducting channel SecY. Features conserved by its two duplicated halves suggest that their common ancestor was an antiparallel homodimeric channel. Structural searches with SecYs halves detect exceptional similarity with the only other ubiquitous protein transporter, YidC. Their shared fold comprises a three-helix bundle interrupted by a helical hairpin. In YidC this hairpin is cytoplasmic and facilitates substrate delivery, whereas in SecY it is transmembrane and forms the substrate-binding lateral gate helices. In both, the three-helix bundle forms a protein-conducting hydrophilic groove, delimited by a conserved hydrophobic residue. We propose that SecY originated as a homodimeric YidC homolog. Many YidC homologs now use this interface to heterodimerise with a conserved partner. Unification of the two ubiquitous protein transporters would reconstruct a key step in the evolution of cells.

category: molecular biology