Proteogenomic analysis of chemo-refractory high-grade serous ovarian cancer.

To improve the understanding of chemo-refractory high-grade serous ovarian cancers (HGSOCs), we characterized the proteogenomic landscape of 242 (refractory and sensitive) HGSOCs, representing one discovery and two validation cohorts across two biospecimen types (formalin-fixed paraffin-embedded and frozen). We identified a 64-protein signature that predicts with high specificity a subset of HGSOCs refractory to initial platinum-based therapy and is validated in two independent patient cohorts. We detected significant association between lack of Ch17 loss of heterozygosity (LOH) and chemo-refractoriness. Based on pathway protein expression, we identified 5 clusters of HGSOC, which validated across two independent patient cohorts and patient-derived xenograft (PDX) models. These clusters may represent different mechanisms of refractoriness and implicate putative therapeutic vulnerabilities.

Prevalence of persistent SARS-CoV-2 in a large community surveillance study.

Persistent SARS-CoV-2 infections may act as viral reservoirs that could seed future outbreaks1-5, give rise to highly divergent lineages6-8 and contribute to cases with post-acute COVID-19 sequelae (long COVID)9,10. However, the population prevalence of persistent infections, their viral load kinetics and evolutionary dynamics over the course of infections remain largely unknown. Here, using viral sequence data collected as part of a national infection survey, we identified 381 individuals with SARS-CoV-2 RNA at high titre persisting for at least 30 days, of which 54 had viral RNA persisting at least 60 days. We refer to these as 'persistent infections' as available evidence suggests that they represent ongoing viral replication, although the persistence of non-replicating RNA cannot be ruled out in all. Individuals with persistent infection had more than 50% higher odds of self-reporting long COVID than individuals with non-persistent infection. We estimate that 0.1-0.5% of infections may become persistent with typically rebounding high viral loads and last for at least 60 days. In some individuals, we identified many viral amino acid substitutions, indicating periods of strong positive selection, whereas others had no consensus change in the sequences for prolonged periods, consistent with weak selection. Substitutions included mutations that are lineage defining for SARS-CoV-2 variants, at target sites for monoclonal antibodies and/or are commonly found in immunocompromised people11-14. This work has profound implications for understanding and characterizing SARS-CoV-2 infection, epidemiology and evolution.

ATGL is a biosynthetic enzyme for fatty acid esters of hydroxy fatty acids.

Branched fatty acid (FA) esters of hydroxy FAs (HFAs; FAHFAs) are recently discovered lipids that are conserved from yeast to mammals1,2. A subfamily, palmitic acid esters of hydroxy stearic acids (PAHSAs), are anti-inflammatory and anti-diabetic1,3. Humans and mice with insulin resistance have lower PAHSA levels in subcutaneous adipose tissue and serum1. PAHSA administration improves glucose tolerance and insulin sensitivity and reduces inflammation in obesity, diabetes and immune-mediated diseases1,4-7. The enzyme(s) responsible for FAHFA biosynthesis in vivo remains unknown. Here we identified adipose triglyceride lipase (ATGL, also known as patatin-like phospholipase domain containing 2 (PNPLA2)) as a candidate biosynthetic enzyme for FAHFAs using chemical biology and proteomics. We discovered that recombinant ATGL uses a transacylation reaction that esterifies an HFA with a FA from triglyceride (TG) or diglyceride to produce FAHFAs. Overexpression of wild-type, but not catalytically dead, ATGL increases FAHFA biosynthesis. Chemical inhibition of ATGL or genetic deletion of Atgl inhibits FAHFA biosynthesis and reduces the levels of FAHFA and FAHFA-TG. Levels of endogenous and nascent FAHFAs and FAHFA-TGs are 80-90 per cent lower in adipose tissue of mice in which Atgl is knocked out specifically in the adipose tissue. Increasing TG levels by upregulating diacylglycerol acyltransferase (DGAT) activity promotes FAHFA biosynthesis, and decreasing DGAT activity inhibits it, reinforcing TGs as FAHFA precursors. ATGL biosynthetic transacylase activity is present in human adipose tissue underscoring its potential clinical relevance. In summary, we discovered the first, to our knowledge, biosynthetic enzyme that catalyses the formation of the FAHFA ester bond in mammals. Whereas ATGL lipase activity is well known, our data establish a paradigm shift demonstrating that ATGL transacylase activity is biologically important.

LncRNA FLAIL affects alternative splicing and represses flowering in Arabidopsis.

How the noncoding genome affects cellular functions is a key biological question. A particular challenge is to distinguish the effects of noncoding DNA elements from long noncoding RNAs (lncRNAs) that coincide at the same loci. Here, we identified the flowering-associated intergenic lncRNA (FLAIL) in Arabidopsis through early flowering flail mutants. Expression of FLAIL RNA from a different chromosomal location in combination with strand-specific RNA knockdown characterized FLAIL as a trans-acting RNA molecule. FLAIL directly binds to differentially expressed target genes that control flowering via RNA-DNA interactions through conserved sequence motifs. FLAIL interacts with protein and RNA components of the spliceosome to affect target mRNA expression through co-transcriptional alternative splicing (AS) and linked chromatin regulation. In the absence of FLAIL, splicing defects at the direct FLAIL target flowering gene LACCASE 8 (LAC8) correlated with reduced mRNA expression. Double mutant analyses support a model where FLAIL-mediated splicing of LAC8 promotes its mRNA expression and represses flowering. Our study suggests lncRNAs as accessory components of the spliceosome that regulate AS and gene expression to impact organismal development.

Beneficial metabolic effects of PAHSAs depend on the gut microbiota in diet-induced obese mice but not in chow-fed mice.

Dietary lipids play an essential role in regulating the function of the gut microbiota and gastrointestinal tract, and these luminal interactions contribute to mediating host metabolism. Palmitic Acid Hydroxy Stearic Acids (PAHSAs) are a family of lipids with antidiabetic and anti-inflammatory properties, but whether the gut microbiota contributes to their beneficial effects on host metabolism is unknown. Here, we report that treating chow-fed female and male germ-free (GF) mice with PAHSAs improves glucose tolerance, but these effects are lost upon high fat diet (HFD) feeding. However, transfer of feces from PAHSA-treated, but not vehicle-treated, chow-fed conventional mice increases insulin sensitivity in HFD-fed GF mice. Thus, the gut microbiota is necessary for, and can transmit, the insulin-sensitizing effects of PAHSAs in HFD-fed GF male mice. Analyses of the cecal metagenome and lipidome of PAHSA-treated mice identified multiple lipid species that associate with the gut commensal Bacteroides thetaiotaomicron (Bt) and with insulin sensitivity resulting from PAHSA treatment. Supplementing live, and to some degree, heat-killed Bt to HFD-fed female mice prevented weight gain, reduced adiposity, improved glucose tolerance, fortified the colonic mucus barrier and reduced systemic inflammation compared to HFD-fed controls. These effects were not observed in HFD-fed male mice. Furthermore, ovariectomy partially reversed the beneficial Bt effects on host metabolism, indicating a role for sex hormones in mediating the Bt probiotic effects. Altogether, these studies highlight the fact that PAHSAs can modulate the gut microbiota and that the microbiota is necessary for the beneficial metabolic effects of PAHSAs in HFD-fed mice.

Composition and function of stress granules and P-bodies in plants.

Stress Granules (SGs) and Processing-bodies (P-bodies) are biomolecular condensates formed in the cell with the highly conserved purpose of maintaining balance between storage, translation, and degradation of mRNA. This balance is particularly important when cells are exposed to different environmental conditions and adjustments have to be made in order for plants to respond to and tolerate stressful conditions. While P-bodies are constitutively present in the cell, SG formation is a stress-induced event. Typically thought of as protein-RNA aggregates, SGs and P-bodies are formed by a process called liquid-liquid phase separation (LLPS), and both their function and composition are very dynamic. Both foci are known to contain proteins involved in translation, protein folding, and ATPase activity, alluding to their roles in regulating mRNA and protein expression levels. From an RNA perspective, SGs and P-bodies primarily consist of mRNAs, though long non-coding RNAs (lncRNAs) have also been observed, and more focus is now being placed on the specific RNAs associated with these aggregates. Recently, metabolites such as nucleotides and amino acids have been reported in purified plant SGs with implications for the energetic dynamics of these condensates. Thus, even though the field of plant SGs and P-bodies is relatively nascent, significant progress has been made in understanding their composition and biological role in stress responses. In this review, we discuss the most recent discoveries centered around SG and P-body function and composition in plants.

Linking discoveries, mechanisms, and technologies to develop a clearer perspective on plant long noncoding RNAs.

Long noncoding RNAs (lncRNAs) are a large and diverse class of genes in eukaryotic genomes that contribute to a variety of regulatory processes. Functionally characterized lncRNAs play critical roles in plants, ranging from regulating flowering to controlling lateral root formation. However, findings from the past decade have revealed that thousands of lncRNAs are present in plant transcriptomes, and characterization has lagged far behind identification. In this setting, distinguishing function from noise is challenging. However, the plant community has been at the forefront of discovery in lncRNA biology, providing many functional and mechanistic insights that have increased our understanding of this gene class. In this review, we examine the key discoveries and insights made in plant lncRNA biology over the past two and a half decades. We describe how discoveries made in the pregenomics era have informed efforts to identify and functionally characterize lncRNAs in the subsequent decades. We provide an overview of the functional archetypes into which characterized plant lncRNAs fit and speculate on new avenues of research that may uncover yet more archetypes. Finally, this review discusses the challenges facing the field and some exciting new molecular and computational approaches that may help inform lncRNA comparative and functional analyses.

Glucose hypometabolism prompts RAN translation and exacerbates C9orf72-related ALS/FTD phenotypes.

The most prevalent genetic cause of both amyotrophic lateral sclerosis and frontotemporal dementia is a (GGGGCC)n nucleotide repeat expansion (NRE) occurring in the first intron of the C9orf72 gene (C9). Brain glucose hypometabolism is consistently observed in C9-NRE carriers, even at pre-symptomatic stages, but its role in disease pathogenesis is unknown. Here, we show alterations in glucose metabolic pathways and ATP levels in the brains of asymptomatic C9-BAC mice. We find that, through activation of the GCN2 kinase, glucose hypometabolism drives the production of dipeptide repeat proteins (DPRs), impairs the survival of C9 patient-derived neurons, and triggers motor dysfunction in C9-BAC mice. We also show that one of the arginine-rich DPRs (PR) could directly contribute to glucose metabolism and metabolic stress by inhibiting glucose uptake in neurons. Our findings provide a potential mechanistic link between energy imbalances and C9-ALS/FTD pathogenesis and suggest a feedforward loop model with potential opportunities for therapeutic intervention.

Global chemical effects of the microbiome include new bile-acid conjugations.

A mosaic of cross-phylum chemical interactions occurs between all metazoans and their microbiomes. A number of molecular families that are known to be produced by the microbiome have a marked effect on the balance between health and disease1-9. Considering the diversity of the human microbiome (which numbers over 40,000 operational taxonomic units10), the effect of the microbiome on the chemistry of an entire animal remains underexplored. Here we use mass spectrometry informatics and data visualization approaches11-13 to provide an assessment of the effects of the microbiome on the chemistry of an entire mammal by comparing metabolomics data from germ-free and specific-pathogen-free mice. We found that the microbiota affects the chemistry of all organs. This included the amino acid conjugations of host bile acids that were used to produce phenylalanocholic acid, tyrosocholic acid and leucocholic acid, which have not previously been characterized despite extensive research on bile-acid chemistry14. These bile-acid conjugates were also found in humans, and were enriched in patients with inflammatory bowel disease or cystic fibrosis. These compounds agonized the farnesoid X receptor in vitro, and mice gavaged with the compounds showed reduced expression of bile-acid synthesis genes in vivo. Further studies are required to confirm whether these compounds have a physiological role in the host, and whether they contribute to gut diseases that are associated with microbiome dysbiosis.

Viral burden is associated with age, vaccination, and viral variant in a population-representative study of SARS-CoV-2 that accounts for time-since-infection-related sampling bias.

In this study, we evaluated the impact of viral variant, in addition to other variables, on within-host viral burden, by analysing cycle threshold (Ct) values derived from nose and throat swabs, collected as part of the UK COVID-19 Infection Survey. Because viral burden distributions determined from community survey data can be biased due to the impact of variant epidemiology on the time-since-infection of samples, we developed a method to explicitly adjust observed Ct value distributions to account for the expected bias. By analysing the adjusted Ct values using partial least squares regression, we found that among unvaccinated individuals with no known prior exposure, viral burden was 44% lower among Alpha variant infections, compared to those with the predecessor strain, B.1.177. Vaccination reduced viral burden by 67%, and among vaccinated individuals, viral burden was 286% higher among Delta variant, compared to Alpha variant, infections. In addition, viral burden increased by 17% for every 10-year age increment of the infected individual. In summary, within-host viral burden increases with age, is reduced by vaccination, and is influenced by the interplay of vaccination status and viral variant.

Correction: Viral burden is associated with age, vaccination, and viral variant in a population-representative study of SARS-CoV-2 that accounts for time-since-infection-related sampling bias.

[This corrects the article DOI: 10.1371/journal.ppat.1011461.].

Development of a mobile, high-throughput, and low-cost image-based plant growth phenotyping system.

Nondestructive plant phenotyping forms a key technique for unraveling molecular processes underlying plant development and response to the environment. While the emergence of high-throughput phenotyping facilities can further our understanding of plant development and stress responses, their high costs greatly hinder scientific progress. To democratize high-throughput plant phenotyping, we developed sets of low-cost image- and weight-based devices to monitor plant shoot growth and evapotranspiration. We paired these devices to a suite of computational pipelines for integrated and straightforward data analysis. The developed tools were validated for their suitability for large genetic screens by evaluating a cowpea (Vigna unguiculata) diversity panel for responses to drought stress. The observed natural variation was used as an input for a genome-wide association study, from which we identified nine genetic loci that might contribute to cowpea drought resilience during early vegetative development. The homologs of the candidate genes were identified in Arabidopsis (Arabidopsis thaliana) and subsequently evaluated for their involvement in drought stress by using available T-DNA insertion mutant lines. These results demonstrate the varied applicability of this low-cost phenotyping system. In the future, we foresee these setups facilitating the identification of genetic components of growth, plant architecture, and stress tolerance across a wide variety of plant species.

Identification and functional annotation of long intergenic non-coding RNAs in Brassicaceae.

Long intergenic noncoding RNAs (lincRNAs) are a large yet enigmatic class of eukaryotic transcripts that can have critical biological functions. The wealth of RNA-sequencing (RNA-seq) data available for plants provides the opportunity to implement a harmonized identification and annotation effort for lincRNAs that enables cross-species functional and genomic comparisons as well as prioritization of functional candidates. In this study, we processed >24 Tera base pairs of RNA-seq data from >16,000 experiments to identify ∼130,000 lincRNAs in four Brassicaceae: Arabidopsis thaliana, Camelina sativa, Brassica rapa, and Eutrema salsugineum. We used nanopore RNA-seq, transcriptome-wide structural information, peptide data, and epigenomic data to characterize these lincRNAs and identify conserved motifs. We then used comparative genomic and transcriptomic approaches to highlight lincRNAs in our data set with sequence or transcriptional conservation. Finally, we used guilt-by-association analyses to assign putative functions to lincRNAs within our data set. We tested this approach on a subset of lincRNAs associated with germination and seed development, observing germination defects for Arabidopsis lines harboring T-DNA insertions at these loci. LincRNAs with Brassicaceae-conserved putative miRNA binding motifs, small open reading frames, or abiotic-stress modulated expression are a few of the annotations that will guide functional analyses into this cryptic portion of the transcriptome.

CyVerse: Cyberinfrastructure for open science.

CyVerse, the largest publicly-funded open-source research cyberinfrastructure for life sciences, has played a crucial role in advancing data-driven research since the 2010s. As the technology landscape evolved with the emergence of cloud computing platforms, machine learning and artificial intelligence (AI) applications, CyVerse has enabled access by providing interfaces, Software as a Service (SaaS), and cloud-native Infrastructure as Code (IaC) to leverage new technologies. CyVerse services enable researchers to integrate institutional and private computational resources, custom software, perform analyses, and publish data in accordance with open science principles. Over the past 13 years, CyVerse has registered more than 124,000 verified accounts from 160 countries and was used for over 1,600 peer-reviewed publications. Since 2011, 45,000 students and researchers have been trained to use CyVerse. The platform has been replicated and deployed in three countries outside the US, with additional private deployments on commercial clouds for US government agencies and multinational corporations. In this manuscript, we present a strategic blueprint for creating and managing SaaS cyberinfrastructure and IaC as free and open-source software.

Telomerase RNA in Hymenoptera (Insecta) switched to plant/ciliate-like biogenesis.

In contrast to the catalytic subunit of telomerase, its RNA subunit (TR) is highly divergent in size, sequence and biogenesis pathways across eukaryotes. Current views on TR evolution assume a common origin of TRs transcribed with RNA polymerase II in Opisthokonta (the supergroup including Animalia and Fungi) and Trypanosomida on one hand, and TRs transcribed with RNA polymerase III under the control of type 3 promoter, found in TSAR and Archaeplastida supergroups (including e.g. ciliates and Viridiplantae taxa, respectively). Here, we focus on unknown TRs in one of the largest Animalia order - Hymenoptera (Arthropoda) with more than 300 available representative genomes. Using a combination of bioinformatic and experimental approaches, we identify their TRs. In contrast to the presumed type of TRs (H/ACA box snoRNAs transcribed with RNA Polymerase II) corresponding to their phylogenetic position, we find here short TRs of the snRNA type, likely transcribed with RNA polymerase III under the control of the type 3 promoter. The newly described insect TRs thus question the hitherto assumed monophyletic origin of TRs across Animalia and point to an evolutionary switch in TR type and biogenesis that was associated with the divergence of Arthropods.

Detection of the Arabidopsis Proteome and Its Post-translational Modifications and the Nature of the Unobserved (Dark) Proteome in PeptideAtlas.

This study describes a new release of the Arabidopsis thaliana PeptideAtlas proteomics resource (build 2023-10) providing protein sequence coverage, matched mass spectrometry (MS) spectra, selected post-translational modifications (PTMs), and metadata. 70 million MS/MS spectra were matched to the Araport11 annotation, identifying ∼0.6 million unique peptides and 18,267 proteins at the highest confidence level and 3396 lower confidence proteins, together representing 78.6% of the predicted proteome. Additional identified proteins not predicted in Araport11 should be considered for the next Arabidopsis genome annotation. This release identified 5198 phosphorylated proteins, 668 ubiquitinated proteins, 3050 N-terminally acetylated proteins, and 864 lysine-acetylated proteins and mapped their PTM sites. MS support was lacking for 21.4% (5896 proteins) of the predicted Araport11 proteome: the "dark" proteome. This dark proteome is highly enriched for E3 ligases, transcription factors, and for certain (e.g., CLE, IDA, PSY) but not other (e.g., THIONIN, CAP) signaling peptides families. A machine learning model trained on RNA expression data and protein properties predicts the probability that proteins will be detected. The model aids in discovery of proteins with short half-life (e.g., SIG1,3 and ERF-VII TFs) and for developing strategies to identify the missing proteins. PeptideAtlas is linked to TAIR, tracks in JBrowse, and several other community proteomics resources.

Cellular and behavioral effects of altered NaV1.2 sodium channel ion permeability in Scn2aK1422E mice.

Genetic variants in SCN2A, encoding the NaV1.2 voltage-gated sodium channel, are associated with a range of neurodevelopmental disorders with overlapping phenotypes. Some variants fit into a framework wherein gain-of-function missense variants that increase neuronal excitability lead to developmental and epileptic encephalopathy, while loss-of-function variants that reduce neuronal excitability lead to intellectual disability and/or autism spectrum disorder (ASD) with or without co-morbid seizures. One unique case less easily classified using this framework is the de novo missense variant SCN2A-p.K1422E, associated with infant-onset developmental delay, infantile spasms and features of ASD. Prior structure-function studies demonstrated that K1422E substitution alters ion selectivity of NaV1.2, conferring Ca2+ permeability, lowering overall conductance and conferring resistance to tetrodotoxin (TTX). Based on heterologous expression of K1422E, we developed a compartmental neuron model incorporating variant channels that predicted reductions in peak action potential (AP) speed. We generated Scn2aK1422E mice and characterized effects on neurons and neurological/neurobehavioral phenotypes. Cultured cortical neurons from heterozygous Scn2aK1422E/+ mice exhibited lower current density with a TTX-resistant component and reversal potential consistent with mixed ion permeation. Recordings from Scn2aK1442E/+ cortical slices demonstrated impaired AP initiation and larger Ca2+ transients at the axon initial segment during the rising phase of the AP, suggesting complex effects on channel function. Scn2aK1422E/+ mice exhibited rare spontaneous seizures, interictal electroencephalogram abnormalities, altered induced seizure thresholds, reduced anxiety-like behavior and alterations in olfactory-guided social behavior. Overall, Scn2aK1422E/+ mice present with phenotypes similar yet distinct from other Scn2a models, consistent with complex effects of K1422E on NaV1.2 channel function.

Disrupting direct inputs from the dorsal subiculum to the granular retrosplenial cortex impairs flexible spatial memory in the rat.

In a changing environment, animals must process spatial signals in a flexible manner. The rat hippocampal formation projects directly upon the retrosplenial cortex, with most inputs arising from the dorsal subiculum and terminating in the granular retrosplenial cortex (area 29). The present study examined whether these same projections are required for spatial working memory and what happens when available spatial cues are altered. Consequently, injections of iDREADDs were made into the dorsal subiculum of rats. In a separate control group, GFP-expressing adeno-associated virus was injected into the dorsal subiculum. Both groups received intracerebral infusions within the retrosplenial cortex of clozapine, which in the iDREADDs rats should selectively disrupt the subiculum to retrosplenial projections. When tested on reinforced T-maze alternation, disruption of the subiculum to retrosplenial projections had no evident effect on the performance of those alternation trials when all spatial-cue types remained present and unchanged. However, the same iDREADDs manipulation impaired performance on all three alternation conditions when there was a conflict or selective removal of spatial cues. These findings reveal how the direct projections from the dorsal subiculum to the retrosplenial cortex support the flexible integration of different spatial cue types, helping the animal to adopt the spatial strategy that best meets current environmental demands.

Comparative analyses suggest a link between mRNA splicing, stability, and RNA covalent modifications in flowering plants.

BACKGROUND: In recent years, covalent modifications on RNA nucleotides have emerged as pivotal moieties influencing the structure, function, and regulatory processes of RNA Polymerase II transcripts such as mRNAs and lncRNAs. However, our understanding of their biological roles and whether these roles are conserved across eukaryotes remains limited. RESULTS: In this study, we leveraged standard polyadenylation-enriched RNA-sequencing data to identify and characterize RNA modifications that introduce base-pairing errors into cDNA reads. Our investigation incorporated data from three Poaceae (Zea mays, Sorghum bicolor, and Setaria italica), as well as publicly available data from a range of stress and genetic contexts in Sorghum and Arabidopsis thaliana. We uncovered a strong enrichment of RNA covalent modifications (RCMs) deposited on a conserved core set of nuclear mRNAs involved in photosynthesis and translation across these species. However, the cohort of modified transcripts changed based on environmental context and developmental program, a pattern that was also conserved across flowering plants. We determined that RCMs can partly explain accession-level differences in drought tolerance in Sorghum, with stress-associated genes receiving a higher level of RCMs in a drought tolerant accession. To address function, we determined that RCMs are significantly enriched near exon junctions within coding regions, suggesting an association with splicing. Intriguingly, we found that these base-pair disrupting RCMs are associated with stable mRNAs, are highly correlated with protein abundance, and thus likely associated with facilitating translation. CONCLUSIONS: Our data point to a conserved role for RCMs in mRNA stability and translation across the flowering plant lineage.