Dynamic regulation of genome-wide pre-mRNA splicing and stress tolerance by the Sm-like protein LSm5 in Arabidopsis

Strain-specific impacts of probiotics are a significant driver of gut microbiome development in very preterm infants

The development of the gut microbiome from birth plays important roles in short- and long-term health, but factors influencing preterm gut microbiome development are poorly understood. In the present study, we use metagenomic sequencing to analyse 1,431 longitudinal stool samples from 123 very preterm infants (<32 weeks' gestation) who did not develop intestinal disease or sepsis over a study period of 10 years. During the study period, one cohort had no probiotic exposure whereas two cohorts were given different probiotic products: Infloran (Bifidobacterium bifidum and Lactobacillus acidophilus) or Labinic (B. bifidum, B. longum subsp. infantis and L. acidophilus). Mothers' own milk, breast milk fortifier, antibiotics and probiotics were significantly associated with the gut microbiome, with probiotics being the most significant factor. Probiotics drove microbiome transition into different preterm gut community types (PGCTs), each enriched in a different Bifidobacterium sp. and significantly associated with increased postnatal age. Functional analyses identified stool metabolites associated with PGCTs and, in preterm-derived organoids, sterile faecal supernatants impacted intestinal, organoid monolayer, gene expression in a PGCT-specific manner. The present study identifies specific influencers of gut microbiome development in very preterm infants, some of which overlap with those impacting term infants. The results highlight the importance of strain-specific differences in probiotic products and their impact on host interactions in the preterm gut.

Tumor microenvironment gene expression profiles associated to complete pathological response and disease progression in resectable NSCLC patients treated with neoadjuvant chemoimmunotherapy

BACKGROUND

Neoadjuvant chemoimmunotherapy for non-small cell lung cancer (NSCLC) has improved pathological responses and survival rates compared with chemotherapy alone, leading to Food and Drug Administration (FDA) approval of nivolumab plus chemotherapy for resectable stage IB-IIIA NSCLC (AJCC 7th edition) without ALK or EGFR alterations. Unfortunately, a considerable percentage of tumors do not completely respond to therapy, which has been associated with early disease progression. So far, it is impossible to predict these events due to lack of knowledge. In this study, we characterized the gene expression profile of tumor samples to identify new biomarkers and mechanisms behind tumor responses to neoadjuvant chemoimmunotherapy and disease recurrence after surgery.

METHODS

Tumor bulk RNA sequencing was performed in 16 pretreatment and 36 post-treatment tissue samples from 41 patients with resectable stage IIIA NSCLC treated with neoadjuvant chemoimmunotherapy from NADIM trial. A panel targeting 395 genes related to immunological processes was used. Tumors were classified as complete pathological response (CPR) and non-CPR, based on the total absence of viable tumor cells in tumor bed and lymph nodes tested at surgery. Differential-expressed genes between groups and pathway enrichment analysis were assessed using DESeq2 and gene set enrichment analysis. CIBERSORTx was used to estimate the proportions of immune cell subtypes.

RESULTS

CPR tumors had a stronger pre-established immune infiltrate at baseline than non-CPR, characterized by higher levels of IFNG, GZMB, NKG7, and M1 macrophages, all with a significant area under the receiver operating characteristic curve (ROC) >0.9 for CPR prediction. A greater effect of neoadjuvant therapy was also seen in CPR tumors with a reduction of tumor markers and IFNγ signaling after treatment. Additionally, the higher expression of several genes, including AKT1, BST2, OAS3, or CD8B; or higher dendritic cells and neutrophils proportions in post-treatment non-CPR samples, were associated with relapse after surgery. Also, high pretreatment PD-L1 and tumor mutational burden levels influenced the post-treatment immune landscape with the downregulation of proliferation markers and type I interferon signaling molecules in surgery samples.

CONCLUSIONS

Our results reinforce the differences between CPR and non-CPR responses, describing possible response and relapse immune mechanisms, opening the possibility of therapy personalization of immunotherapy-based regimens in the neoadjuvant setting of NSCLC.

Immune and Genetic Signatures of Breast Carcinomas Triggering Anti-Yo–Associated Paraneoplastic Cerebellar Degeneration

Background and Objectives

Paraneoplastic cerebellar degeneration (PCD) with anti-Yo antibodies is a cancer-related autoimmune disease directed against neural antigens expressed by tumor cells. A putative trigger of the immune tolerance breakdown is genetic alteration of Yo antigens. We aimed to identify the tumors' genetic and immune specificities involved in Yo-PCD pathogenesis.

Methods

Using clinicopathologic data, immunofluorescence (IF) imaging, and whole-transcriptome analysis, 22 breast cancers (BCs) associated with Yo-PCD were characterized in terms of oncologic characteristics, genetic alteration of Yo antigens, differential gene expression profiles, and morphofunctional specificities of their in situ antitumor immunity by comparing them with matched control BCs.

Results

Yo-PCD BCs were invasive carcinoma of no special type, which early metastasized to lymph nodes. They overexpressed human epidermal growth factor receptor 2 (HER2) but were hormone receptor negative. All Yo-PCD BCs carried at least 1 genetic alteration (variation or gain in copy number) on CDR2L, encoding the main Yo antigen that was found aberrantly overexpressed in Yo-PCD BCs. Analysis of the differentially expressed genes found 615 upregulated and 54 downregulated genes in Yo-PCD BCs compared with HER2-driven control BCs without PCD. Ontology enrichment analysis found significantly upregulated adaptive immune response pathways in Yo-PCD BCs. IF imaging confirmed an intense immune infiltration with an overwhelming predominance of immunoglobulin G–plasma cells.

Discussion

These data confirm the role of genetic alterations of Yo antigens in triggering the immune tolerance breakdown but also outline a specific biomolecular profile in Yo-PCD BCs, suggesting a cancer-specific pathogenesis.

Differential microRNA expression analyses across two brain regions in Alzheimer's disease

Dysregulation of microRNAs (miRNAs) is involved in the pathogenesis of neurodegenerative diseases, including Alzheimer's disease (AD). Hitherto, sample sizes from differential miRNA expression studies in AD are exceedingly small aggravating any biological inference. To overcome this limitation, we investigated six candidate miRNAs in a large collection of brain samples. Brain tissue was derived from superior temporal gyrus (STG) and entorhinal cortex (EC) from 99 AD patients and 91 controls. MiRNA expression was examined by qPCR (STG) or small RNA sequencing (EC). Brain region-dependent differential miRNA expression was investigated in a transgenic AD mouse model using qPCR and FISH. Total RNA sequencing was used to assess differential expression of miRNA target genes. MiR-129-5p, miR-132-5p, and miR-138-5p were significantly downregulated in AD vs. controls both in STG and EC, while miR-125b-5p and miR-501-3p showed no evidence for differential expression in this dataset. In addition, miR-195-5p was significantly upregulated in EC but not STG in AD patients. The brain region-specific pattern of miR-195-5p expression was corroborated in vivo in transgenic AD mice. Total RNA sequencing identified several novel and functionally interesting target genes of these miRNAs involved in synaptic transmission (GABRB1), the immune-system response (HCFC2) or AD-associated differential methylation (SLC16A3). Using two different methods (qPCR and small RNA-seq) in two separate brain regions in 190 individuals we more than doubled the available sample size for most miRNAs tested. Differential gene expression analyses confirm the likely involvement of miR-129-5p, miR-132-5p, miR-138-5p, and miR-195-5p in AD pathogenesis and highlight several novel potentially relevant target mRNAs.

Transcriptome-guided annotation and functional classification of long non-coding RNAs in Arabidopsis thaliana

Long non-coding RNAs (lncRNAs) are a prominent class of eukaryotic regulatory genes. Despite the numerous available transcriptomic datasets, the annotation of plant lncRNAs remains based on dated annotations that have been historically carried over. We present a substantially improved annotation of Arabidopsis thaliana lncRNAs, generated by integrating 224 transcriptomes in multiple tissues, conditions, and developmental stages. We annotate 6764 lncRNA genes, including 3772 that are novel. We characterize their tissue expression patterns and find 1425 lncRNAs are co-expressed with coding genes, with enriched functional categories such as chloroplast organization, photosynthesis, RNA regulation, transcription, and root development. This improved transcription-guided annotation constitutes a valuable resource for studying lncRNAs and the biological processes they may regulate.

The gill transcriptome of threatened European freshwater mussels

Genomic tools applied to non-model organisms are critical to design successful conservation strategies of particularly threatened groups. Freshwater mussels of the Unionida order are among the most vulnerable taxa and yet almost no genetic resources are available. Here, we present the gill transcriptomes of five European freshwater mussels with high conservation concern: Margaritifera margaritifera, Unio crassus, Unio pictorum, Unio mancus and Unio delphinus. The final assemblies, with N50 values ranging from 1069–1895 bp and total BUSCO scores above 90% (Eukaryote and Metazoan databases), were structurally and functionally annotated, and made available. The transcriptomes here produced represent a valuable resource for future studies on these species’ biology and ultimately guide their conservation.

Measurement(s)	transcriptomics
Technology Type(s)	Illumina sequencing
Sample Characteristic - Organism	Margaritifera margaritifera • Unio crassus • Unio delphinus • Unio mancus • Unio pictorum
Sample Characteristic - Location	Europe

Duplex sequencing identifies genomic features that determine susceptibility to benzo(a)pyrene-induced in vivo mutations

Exposure to environmental mutagens increases the risk of cancer and genetic disorders. We used Duplex Sequencing (DS), a high-accuracy error-corrected sequencing technology, to analyze mutation induction across twenty 2.4 kb intergenic and genic targets in the bone marrow of MutaMouse males exposed to benzo(a)pyrene (BaP), a widespread environmental pollutant. DS revealed a linear dose-related induction of mutations across all targets with low intra-group variability. Heterochromatic and intergenic regions exhibited the highest mutation frequencies (MF). C:G > A:T transversions at CCA, CCC and GCC trinucleotides were enriched in BaP-exposed mice consistent with the known etiology of BaP mutagenesis. However, GC-content had no effect on mutation susceptibility. A positive correlation was observed between DS and the “gold-standard” transgenic rodent gene mutation assay. Overall, we demonstrate that DS is a promising approach to study in vivo mutagenesis and yields critical insight into the genomic features governing mutation susceptibility, spectrum, and variability across the genome.

Tissue-specific transcriptome profiles identify functional differences key to understanding whole plant response to life in variable salinity

Plants endure environmental stressors via adaptation and phenotypic plasticity. Studying these mechanisms in seagrasses is extremely relevant as they are important primary producers and functionally significant carbon sinks. These mechanisms are not well understood at the tissue level in seagrasses. Using RNA-seq, we generated transcriptome sequences from tissue of leaf, basal leaf meristem and root organs of Posidonia australis, establishing baseline in situ transcriptomic profiles for tissues across a salinity gradient. Samples were collected from four P. australis meadows growing in Shark Bay, Western Australia. Analysis of gene expression showed significant differences between tissue types, with more variation among leaves than meristem or roots. Gene ontology enrichment analysis showed the differences were largely due to the role of photosynthesis, plant growth and nutrient absorption in leaf, meristem and root organs, respectively. Differential gene expression of leaf and meristem showed upregulation of salinity regulation processes in higher salinity meadows. Our study highlights the importance of considering leaf meristem tissue when evaluating whole-plant responses to environmental change.

This article has an associated First Person interview with the first author of the paper.

Summary: Differences in seagrass leaf, meristem and root transcriptomes across variable salinities are due to tissue-specific processes. Leaf meristem contained the broadest process range, indicating preferential use for inferring plant-wide activity.

Metabolome and transcriptome profiling reveal regulatory network and mechanism of flavonoid biosynthesis during color formation of Dioscorea cirrhosa L

Dioscorea cirrhosa is a plant that is used as a dye as well as in medicine. Many metabolites with pharmacological activity exist in the tubers of D. cirrhosa. However, little is known about the mechanism regulating biosynthesis in these metabolites. In this study, transcriptome and metabolome profiling were performed in four color tubers. A total of 531 metabolites, including 62 flavonoids, were identified. Epicatechin and proanthocyanin B2 were the key metabolites that exhibited high content levels in the four tubers. These metabolites were divided into nine classes with distinct change patterns. A total of 22,865 differentially expressed genes (DEGs) were identified by transcriptome analysis. Among these DEGs, we identified 67 candidate genes related to the flavonoid biosynthesis pathway and three genes that played pivotal roles in proanthocyanin (PA) synthesis. A weighted gene co-expression network analysis (WGCNA) revealed that the two modules, "MEblue" and "MEblack," were two key gene sets strongly associated with phenylpropanoid and flavonoid biosynthesis. We also found that the plant hormone signal transduction biological process exhibited activity in the late stage of tuber color formation. Additionally, we identified 37 hub transcript factors related to flavonoid biosynthesis, of which 24 were found to be highly associated with flavonoid pathway genes. In addition to the MYB-bHLH-WD40 (MBW) genes, we found that the plant hormone gene families exhibited high expression levels. This study provides a reference for understanding the synthesis of D. cirrhosa tuber metabolites at the molecular level and provides a foundation for the further development of D. cirrhosa related plant pigments as well as its further use in the pharmaceutical industry.

Constitutive expression of a pea apyrase, psNTP9, increases seed yield in field-grown soybean

To address the demand for food by a rapidly growing human population, agricultural scientists have carried out both plant breeding and genetic engineering research. Previously, we reported that the constitutive expression of a pea apyrase (Nucleoside triphosphate, diphosphohydrolase) gene, psNTP9, under the control of the CaMV35S promoter, resulted in soybean plants with an expanded root system architecture, enhanced drought resistance and increased seed yield when they are grown in greenhouses under controlled conditions. Here, we report that psNTP9-expressing soybean lines also show significantly enhanced seed yields when grown in multiple different field conditions at multiple field sites, including when the gene is introgressed into elite germplasm. The transgenic lines have higher leaf chlorophyll and soluble protein contents and decreased stomatal density and cuticle permeability, traits that increase water use efficiency and likely contribute to the increased seed yields of field-grown plants. These altered properties are explained, in part, by genome-wide gene expression changes induced by the transgene.

Relevance and Regulation of Alternative Splicing in Plant Heat Stress Response: Current Understanding and Future Directions

Alternative splicing (AS) is a major mechanism for gene expression in eukaryotes, increasing proteome diversity but also regulating transcriptome abundance. High temperatures have a strong impact on the splicing profile of many genes and therefore AS is considered as an integral part of heat stress response. While many studies have established a detailed description of the diversity of the RNAome under heat stress in different plant species and stress regimes, little is known on the underlying mechanisms that control this temperature-sensitive process. AS is mainly regulated by the activity of splicing regulators. Changes in the abundance of these proteins through transcription and AS, post-translational modifications and interactions with exonic and intronic cis-elements and core elements of the spliceosomes modulate the outcome of pre-mRNA splicing. As a major part of pre-mRNAs are spliced co-transcriptionally, the chromatin environment along with the RNA polymerase II elongation play a major role in the regulation of pre-mRNA splicing under heat stress conditions. Despite its importance, our understanding on the regulation of heat stress sensitive AS in plants is scarce. In this review, we summarize the current status of knowledge on the regulation of AS in plants under heat stress conditions. We discuss possible implications of different pathways based on results from non-plant systems to provide a perspective for researchers who aim to elucidate the molecular basis of AS under high temperatures.

Local and systemic transcriptomic responses from acute exercise induced muscle damage of the human knee extensors

Skeletal muscle is adaptable to a direct stimulus of exercise-induced muscle damage (EIMD). Local muscle gene networks and systemic circulatory factors respond to EIMD within days, mediating anti-inflammation and cellular proliferation. Here we show in humans that local EIMD of one muscle group is associated with a systemic response of gene networks that regulate muscle structure and cellular development in non-local homologous muscle not directly altered by EIMD. In the non-dominant knee-extensors of seven males, EIMD was induced through voluntary contractions against an electric motor that lengthened muscles. Neuromuscular assessments, vastus lateralis muscle biopsies and blood draws occurred at two days prior, and one and two days post the EIMD intervention. From the muscle and blood plasma samples, RNA-seq measured transcriptome changes of differential expression using bioinformatic analyses.Relative to the time of the EIMD intervention, local muscle that was mechanically damaged had 475 genes differentially expressed, as compared to 33 genes in the non-local homologous muscle. Gene and network analysis showed that activity of the local muscle was related to structural maintenance, repair, and energetic processes, whereas gene and network activity of the non-local muscle (that was not directly modified by the EIMD) were related to muscle cell development, stress response, and structural maintenance. Altered expression of two novel miRNAs related to the EIMD response supported that systemic factors were active. Together, these results indicate that the expression of genes and gene networks that control muscle contractile structure can be modified in response to non-local EIMD in humans.

Use of a graph neural network to the weighted gene co-expression network analysis of Korean native cattle

In the general framework of the weighted gene co-expression network analysis (WGCNA), a hierarchical clustering algorithm is commonly used to module definition. However, hierarchical clustering depends strongly on the topological overlap measure. In other words, this algorithm may assign two genes with low topological overlap to different modules even though their expression patterns are similar. Here, a novel gene module clustering algorithm for WGCNA is proposed. We develop a gene module clustering network (gmcNet), which simultaneously addresses single-level expression and topological overlap measure. The proposed gmcNet includes a “co-expression pattern recognizer” (CEPR) and “module classifier”. The CEPR incorporates expression features of single genes into the topological features of co-expressed ones. Given this CEPR-embedded feature, the module classifier computes module assignment probabilities. We validated gmcNet performance using 4,976 genes from 20 native Korean cattle. We observed that the CEPR generates more robust features than single-level expression or topological overlap measure. Given the CEPR-embedded feature, gmcNet achieved the best performance in terms of modularity (0.261) and the differentially expressed signal (27.739) compared with other clustering methods tested. Furthermore, gmcNet detected some interesting biological functionalities for carcass weight, backfat thickness, intramuscular fat, and beef tenderness of Korean native cattle. Therefore, gmcNet is a useful framework for WGCNA module clustering.

3-oxo-C12:2-HSL, quorum sensing molecule from human intestinal microbiota, inhibits pro-inflammatory pathways in immune cells via bitter taste receptors

In the gut ecosystem, microorganisms regulate group behaviour and interplay with the host via a molecular system called quorum sensing (QS). The QS molecule 3-oxo-C12:2-HSL, first identified in human gut microbiota, exerts anti-inflammatory effects and could play a role in inflammatory bowel diseases where dysbiosis has been described. Our aim was to identify which signalling pathways are involved in this effect. We observed that 3-oxo-C12:2-HSL decreases expression of pro-inflammatory cytokines such as Interleukine-1β (− 35%) and Tumor Necrosis Factor-α (TNFα) (− 40%) by stimulated immune RAW264.7 cells and decreased TNF secretion by stimulated PBMC in a dose-dependent manner, between 25 to 100 µM. Transcriptomic analysis of RAW264.7 cells exposed to 3-oxo-C12:2-HSL, in a pro-inflammatory context, highlighted JAK-STAT, NF-κB and TFN signalling pathways and we confirmed that 3-oxo-C12:2-HSL inhibited JAK1 and STAT1 phosphorylation. We also showed through a screening assay that 3-oxo-C12:2-HSL interacted with several human bitter taste receptors. Its anti-inflammatory effect involved TAS2R38 as shown by pharmacologic inhibition and led to an increase in intracellular calcium levels. We thus unravelled the involvement of several cellular pathways in the anti-inflammatory effects exerted by the QS molecule 3-oxo-C12:2-HSL.

Microbial Risk Score for Capturing Microbial Characteristics, Integrating Multi-omics Data, and Predicting Disease Risk

Background:

With the rapid accumulation of microbiome-wide association studies, a great amount of microbiome data are available to study the microbiome’s role in human disease and advance the microbiome’s potential use for disease prediction. However, the unique features of microbiome data hinder its utility for disease prediction.

Methods:

Motivated from the polygenic risk score framework, we propose a microbial risk score (MRS) framework to aggregate the complicated microbial profile into a summarized risk score that can be used to measure and predict disease susceptibility. Specifically, the MRS algorithm involves two steps: 1) identifying a sub-community consisting of the signature microbial taxa associated with disease, and 2) integrating the identified microbial taxa into a continuous score. The first step is carried out using the existing sophisticated microbial association tests and pruning and thresholding method in the discovery samples. The second step constructs a community-based MRS by calculating alpha diversity on the identified sub-community in the validation samples. Moreover, we propose a multi-omics data integration method by jointly modeling the proposed MRS and other risk scores constructed from other omics data in disease prediction.

Results:

Through three comprehensive real data analyses using the NYU Langone Health COVID-19 cohort, the gut microbiome health index (GMHI) multi-study cohort, and a large type 1 diabetes cohort separately, we exhibit and evaluate the utility of the proposed MRS framework for disease prediction and multi-omics data integration. In addition, the disease-specific MRSs for colorectal adenoma, colorectal cancer, Crohn’s disease, and rheumatoid arthritis based on the relative abundances of 5, 6, 12, and 6 microbial taxa respectively are created and validated using the GMHI multi-study cohort. Especially, Crohn’s disease MRS achieves AUCs of 0.88 ([0.85–0.91]) and 0.86 ([0.78–0.95]) in the discovery and validation cohorts, respectively.

Conclusions:

The proposed MRS framework sheds light on the utility of the microbiome data for disease prediction and multi-omics integration, and provides great potential in understanding the microbiome’s role in disease diagnosis and prognosis.

The Rice Serine/Arginine Splicing Factor RS33 Regulates Pre-mRNA Splicing during Abiotic Stress Responses

Abiotic stresses profoundly affect plant growth and development and limit crop productivity. Pre-mRNA splicing is a major form of gene regulation that helps plants cope with various stresses. Serine/arginine (SR)-rich splicing factors play a key role in pre-mRNA splicing to regulate different biological processes under stress conditions. Alternative splicing (AS) of SR transcripts and other transcripts of stress-responsive genes generates multiple splice isoforms that contribute to protein diversity, modulate gene expression, and affect plant stress tolerance. Here, we investigated the function of the plant-specific SR protein RS33 in regulating pre-mRNA splicing and abiotic stress responses in rice. The loss-of-function mutant rs33 showed increased sensitivity to salt and low-temperature stresses. Genome-wide analyses of gene expression and splicing in wild-type and rs33 seedlings subjected to these stresses identified multiple splice isoforms of stress-responsive genes whose AS are regulated by RS33. The number of RS33-regulated genes was much higher under low-temperature stress than under salt stress. Our results suggest that the plant-specific splicing factor RS33 plays a crucial role during plant responses to abiotic stresses.

Plant hormone regulation of abiotic stress responses

Plant hormones are signalling compounds that regulate crucial aspects of growth, development and environmental stress responses. Abiotic stresses, such as drought, salinity, heat, cold and flooding, have profound effects on plant growth and survival. Adaptation and tolerance to such stresses require sophisticated sensing, signalling and stress response mechanisms. In this Review, we discuss recent advances in understanding how diverse plant hormones control abiotic stress responses in plants and highlight points of hormonal crosstalk during abiotic stress signalling. Control mechanisms and stress responses mediated by plant hormones including abscisic acid, auxin, brassinosteroids, cytokinins, ethylene and gibberellins are discussed. We discuss new insights into osmotic stress sensing and signalling mechanisms, hormonal control of gene regulation and plant development during stress, hormone-regulated submergence tolerance and stomatal movements. We further explore how innovative imaging approaches are providing insights into single-cell and tissue hormone dynamics. Understanding stress tolerance mechanisms opens new opportunities for agricultural applications.

Comparative whole transcriptome analysis of Parkinson's disease focusing on the efficacy of zonisamide

OBJECTIVE

Interindividual variations in responsiveness to zonisamide in patients with Parkinson's disease (PD) have been observed in clinical settings. To decipher the molecular mechanisms determining the efficacy of zonisamide, we conducted whole transcriptome sequencing analysis of patients with PD.

METHODS

We selected 23 super-responders (SRs) and 25 non-responders (NRs) to zonisamide from patients with PD who had participated in a previous clinical trial for the approval of zonisamide for the treatment of 'wearing-off'. Whole transcriptome analysis of peripheral blood was conducted on samples taken before and 12 weeks after zonisamide treatment. We performed differential gene expression analysis to compare between the SRs and NRs at each time point.

RESULTS

Differentially expressed genes in the pre-treatment samples were significantly enriched for glutamatergic synapses and insulin-like growth factor binding (P_adj=7.8 × 10^-3 and 0.029, respectively). The gene sets associated with these functions changed more dynamically by treatment in SRs than NRs (p=7.2 × 10^-3 and 8.2 × 10^-3, respectively).

CONCLUSIONS

Our results suggest that the efficacy of zonisamide in PD patients is associated with glutamate-related synaptic modulation and p53-mediated dopaminergic neural loss. Their transcriptomic differences could be captured before treatment, which would lead to the realisation of future personalised treatment.

An obesogenic feedforward loop involving PPARγ, acyl-CoA binding protein and GABAA receptor

Acyl-coenzyme-A-binding protein (ACBP), also known as a diazepam-binding inhibitor (DBI), is a potent stimulator of appetite and lipogenesis. Bioinformatic analyses combined with systematic screens revealed that peroxisome proliferator-activated receptor gamma (PPARγ) is the transcription factor that best explains the ACBP/DBI upregulation in metabolically active organs including the liver and adipose tissue. The PPARγ agonist rosiglitazone-induced ACBP/DBI upregulation, as well as weight gain, that could be prevented by knockout of Acbp/Dbi in mice. Moreover, liver-specific knockdown of Pparg prevented the high-fat diet (HFD)-induced upregulation of circulating ACBP/DBI levels and reduced body weight gain. Conversely, knockout of Acbp/Dbi prevented the HFD-induced upregulation of PPARγ. Notably, a single amino acid substitution (F77I) in the γ2 subunit of gamma-aminobutyric acid A receptor (GABA_AR), which abolishes ACBP/DBI binding to this receptor, prevented the HFD-induced weight gain, as well as the HFD-induced upregulation of ACBP/DBI, GABA_AR γ2, and PPARγ. Based on these results, we postulate the existence of an obesogenic feedforward loop relying on ACBP/DBI, GABA_AR, and PPARγ. Interruption of this vicious cycle, at any level, indistinguishably mitigates HFD-induced weight gain, hepatosteatosis, and hyperglycemia.

Bovine host genome acts on rumen microbiome function linked to methane emissions

Our study provides substantial evidence that the host genome affects the comprehensive function of the microbiome in the rumen of bovines. Of 1,107/225/1,141 rumen microbial genera/metagenome assembled uncultured genomes (RUGs)/genes identified from whole metagenomics sequencing, 194/14/337 had significant host genomic effects (heritabilities ranging from 0.13 to 0.61), revealing that substantial variation of the microbiome is under host genomic control. We found 29/22/115 microbial genera/RUGs/genes host-genomically correlated (|0.59| to |0.93|) with emissions of the potent greenhouse gas methane (CH4), highlighting the strength of a common host genomic control of specific microbial processes and CH4. Only one of these microbial genes was directly involved in methanogenesis (cofG), whereas others were involved in providing substrates for archaea (e.g. bcd and pccB), important microbial interspecies communication mechanisms (ABC.PE.P), host-microbiome interaction (TSTA3) and genetic information processes (RP-L35). In our population, selection based on abundances of the 30 most informative microbial genes provided a mitigation potential of 17% of mean CH4 emissions per generation, which is higher than for selection based on measured CH4 using respiration chambers (13%), indicating the high potential of microbiome-driven breeding to cumulatively reduce CH4 emissions and mitigate climate change.

Sucrose addition directionally enhances bacterial community convergence and network stability of the shrimp culture system

Sucrose is an effective carbon source for creating more reliable and environmentally friendly conditions for shrimp growth by regulating bacteria in biofloc-based culture systems. However, the influence of sucrose addition on the interaction, co-occurrence networks, and assembly mechanisms of bacterial communities in biofloc-based culture systems remains largely unknown. Here, we comprehensively investigated the effects of sucrose addition on bacterial communities in three habitats (water, bioflocs, and gut). The bacterial community structures and compositions of these three habitats became more similar in groups with sucrose addition, compared with those in controls. More than 50% gut bacterial communities were mainly derived from water and biofloc communities in the sucrose addition groups, but only about 33% bacterial communities migrated from water and biofloc to the gut in the control culture system. Sucrose addition accordantly enriched core taxa belonging to the phylum Actinobacteria and the families Rhodobacteraceae and Flavobacteriaceae in water, biofloc, and gut habitats. These core taxa were important for maintaining bacterial network stability in the sucrose addition culture systems and some were identified as keystone taxa for improving shrimp growth. Furthermore, after sucrose addition, gut bacterial community assembly from water and biofloc was dominated by the heterogeneous select with the ratios of 55–91% and 67–83%, respectively, indicating that sucrose addition can directionally shape the bacterial assembly of the shrimp culture system. These results provide a basis for selectively regulating certain beneficial taxa to improve shrimp growth in culture systems.

Ubiquitin ligase STUB1 destabilizes IFNγ-receptor complex to suppress tumor IFNγ signaling

The cytokine IFNγ differentially impacts on tumors upon immune checkpoint blockade (ICB). Despite our understanding of downstream signaling events, less is known about regulation of its receptor (IFNγ-R1). With an unbiased genome-wide CRISPR/Cas9 screen for critical regulators of IFNγ-R1 cell surface abundance, we identify STUB1 as an E3 ubiquitin ligase for IFNγ-R1 in complex with its signal-relaying kinase JAK1. STUB1 mediates ubiquitination-dependent proteasomal degradation of IFNγ-R1/JAK1 complex through IFNγ-R1K285 and JAK1K249. Conversely, STUB1 inactivation amplifies IFNγ signaling, sensitizing tumor cells to cytotoxic T cells in vitro. This is corroborated by an anticorrelation between STUB1 expression and IFNγ response in ICB-treated patients. Consistent with the context-dependent effects of IFNγ in vivo, anti-PD-1 response is increased in heterogenous tumors comprising both wildtype and STUB1-deficient cells, but not full STUB1 knockout tumors. These results uncover STUB1 as a critical regulator of IFNγ-R1, and highlight the context-dependency of STUB1-regulated IFNγ signaling for ICB outcome.

Cell environment shapes TDP-43 function with implications in neuronal and muscle disease

TDP-43 (TAR DNA-binding protein 43) aggregation and redistribution are recognised as a hallmark of amyotrophic lateral sclerosis and frontotemporal dementia. As TDP-43 inclusions have recently been described in the muscle of inclusion body myositis patients, this highlights the need to understand the role of TDP-43 beyond the central nervous system. Using RNA-seq, we directly compare TDP-43-mediated RNA processing in muscle (C2C12) and neuronal (NSC34) mouse cells. TDP-43 displays a cell-type-characteristic behaviour targeting unique transcripts in each cell-type, which is due to characteristic expression of RNA-binding proteins, that influence TDP-43's performance and define cell-type specific splicing. Among splicing events commonly dysregulated in both cell lines, we identify some that are TDP-43-dependent also in human cells. Inclusion levels of these alternative exons are altered in tissues of patients suffering from FTLD and IBM. We therefore propose that TDP-43 dysfunction contributes to disease development either in a common or a tissue-specific manner.

The PIWI/piRNA response is relaxed in a rodent that lacks mobilizing transposable elements

Transposable elements (TEs) are genomic parasites that can propagate throughout host genomes. Mammalian genomes are typically dominated by LINE retrotransposons and their associated SINEs, and germline mobilization is a challenge to genome integrity. There are defenses against TE proliferation and the PIWI/piRNA defense is among the most well understood. However, the PIWI/piRNA system has been investigated largely in animals with actively mobilizing TEs and it is unclear how the PIWI/piRNA system functions in the absence of mobilizing TEs. The 13-lined ground squirrel provides the opportunity to examine PIWI/piRNA and TE dynamics within the context of minimal, and possibly nonexistent, TE accumulation. To do so, we compared the PIWI/piRNA dynamics in squirrels to observations from the rabbit and mouse. Despite a lack of young insertions in squirrels, TEs were still actively transcribed at higher levels compared to mouse and rabbit. All three Piwi genes were not expressed, prior to P8 in squirrel testis, and there was little TE expression change with the onset of Piwi expression. We also demonstrated there was not a major expression change in the young squirrel LINE families in the transition from juvenile to adult testis in contrast to young mouse and rabbit LINE families. These observations lead us to conclude that PIWI suppression, was weaker for squirrel LINEs and SINEs and did not strongly reduce their transcription. We speculate that, although the PIWI/piRNA system is adaptable to novel TE threats, transcripts from TEs that are no longer threatening receive less attention from PIWI proteins.

Transcriptional responses of Trichodesmium to natural inverse gradients of Fe and P availability

The filamentous diazotrophic cyanobacterium Trichodesmium is responsible for a significant fraction of marine di-nitrogen (N₂) fixation. Growth and distribution of Trichodesmium and other diazotrophs in the vast oligotrophic subtropical gyres is influenced by iron (Fe) and phosphorus (P) availability, while reciprocally influencing the biogeochemistry of these nutrients. Here we use observations across natural inverse gradients in Fe and P in the North Atlantic subtropical gyre (NASG) to demonstrate how Trichodesmium acclimates in situ to resource availability. Transcriptomic analysis identified progressive upregulation of known iron-stress biomarker genes with decreasing Fe availability, and progressive upregulation of genes involved in the acquisition of diverse P sources with decreasing P availability, while genes involved in N₂ fixation were upregulated at the intersection under moderate Fe and P availability. Enhanced N₂ fixation within the Fe and P co-stressed transition region was also associated with a distinct, consistent metabolic profile, including the expression of alternative photosynthetic pathways that potentially facilitate ATP generation required for N₂ fixation with reduced net oxygen production. The observed response of Trichodesmium to availability of both Fe and P supports suggestions that these biogeochemically significant organisms employ unique molecular, and thus physiological responses as adaptations to specifically exploit the Fe and P co-limited niche they construct.

Transcriptome landscape of early inflorescence developmental stages identifies key flowering time regulators in chickpea

Transcriptome landscape during early inflorescence developmental stages identified candidate flowering time regulators including Early Flowering 3a. Further genomics approaches validated the role of this gene in flowering time regulation. The early stages of inflorescence development in plants are as crucial as the later floral developmental stages. Several traits, such as inflorescence architecture and flower developmental timings, are determined during those early stages. In chickpea, diverse forms of inflorescence architectures regarding meristem determinacy and the number of flowers per node are observed within the germplasm. Transcriptome analysis in four desi chickpea accessions with such unique inflorescence characteristics identifies the underlying shared regulatory events leading to inflorescence development. The vegetative to reproductive stage transition brings about major changes in the transcriptome landscape. The inflorescence development progression associated genes identified through co-expression network analysis includes both protein-coding genes and long non-coding RNAs (lncRNAs). Few lncRNAs identified in our study positively regulate flowering-related mRNA stability by acting competitively with miRNAs. Bulk segregrant analysis and association mapping narrowed down an InDel marker regulating flowering time in chickpea. Deletion of 11 bp in first exon of a negative flowering time regulator, Early Flowering 3a gene, leads to early flowering phenotype in chickpea. Understanding the key players involved in vegetative to reproductive stage transition and floral meristem development will be useful in manipulating flowering time and inflorescence architecture in chickpea and other legumes.

ABA Mediates Plant Development and Abiotic Stress via Alternative Splicing

Alternative splicing (AS) exists in eukaryotes to increase the complexity and adaptability of systems under biophysiological conditions by increasing transcriptional and protein diversity. As a classic hormone, abscisic acid (ABA) can effectively control plant growth, improve stress resistance, and promote dormancy. At the transcriptional level, ABA helps plants respond to the outside world by regulating transcription factors through signal transduction pathways to regulate gene expression. However, at the post-transcriptional level, the mechanism by which ABA can regulate plant biological processes by mediating alternative splicing is not well understood. Therefore, this paper briefly introduces the mechanism of ABA-induced alternative splicing and the role of ABA mediating AS in plant response to the environment and its own growth.

Knowledge graph-based recommendation framework identifies drivers of resistance in EGFR mutant non-small cell lung cancer

Resistance to EGFR inhibitors (EGFRi) presents a major obstacle in treating non-small cell lung cancer (NSCLC). One of the most exciting new ways to find potential resistance markers involves running functional genetic screens, such as CRISPR, followed by manual triage of significantly enriched genes. This triage process to identify ‘high value’ hits resulting from the CRISPR screen involves manual curation that requires specialized knowledge and can take even experts several months to comprehensively complete. To find key drivers of resistance faster we build a recommendation system on top of a heterogeneous biomedical knowledge graph integrating pre-clinical, clinical, and literature evidence. The recommender system ranks genes based on trade-offs between diverse types of evidence linking them to potential mechanisms of EGFRi resistance. This unbiased approach identifies 57 resistance markers from >3,000 genes, reducing hit identification time from months to minutes. In addition to reproducing known resistance markers, our method identifies previously unexplored resistance mechanisms that we prospectively validate.

Resistance to EGFR inhibitors presents a major obstacle in treating non-small cell lung cancer. Here, the authors develop a recommender system ranking genes based on trade-offs between diverse types of evidence linking them to potential mechanisms of EGFRi resistance.

Newcastle disease virus genotype VII gene expression in experimentally infected birds

Newcastle disease virus genotype VII (NDV-GVII) is a highly contagious pathogen responsible for pandemics that have caused devastating economic losses in the poultry industry. Several features in the transcription of NDV mRNA, including differentially expressed genes across the viral genome, are shared with that for other single, non-segmented, negative-strand viruses. Previous studies measuring viral gene expression using northern blotting indicated that the NDV transcription produced non-equimolar levels of viral mRNAs. However, deep high-throughput sequencing of virus-infected tissues can provide a better insight into the patterns of viral transcription. In this report, the transcription pattern of virulent NDV-GVII was analysed using RNA-seq and qRT-PCR. This study revealed the transcriptional profiling of these highly pathogenic NDV-GVII genes: NP:P:M:F:HN:L, in which there was a slight attenuation at the NP:P and HN:L gene boundaries. Our result also provides a fully comprehensive qPCR protocol for measuring viral transcript abundance that may be more convenient for laboratories where accessing RNA-seq is not feasible.

Oncolytic Orf virus licenses NK cells via cDC1 to activate innate and adaptive antitumor mechanisms and extends survival in a murine model of late-stage ovarian cancer

BACKGROUND

Novel therapies are needed to improve outcomes for women diagnosed with ovarian cancer. Oncolytic viruses are multifunctional immunotherapeutic biologics that preferentially infect cancer cells and stimulate inflammation with the potential to generate antitumor immunity. Herein we describe Parapoxvirus ovis (Orf virus (OrfV)), an oncolytic poxvirus, as a viral immunotherapy for ovarian cancer.

METHODS

The immunotherapeutic potential of OrfV was tested in the ID8 orthotopic mouse model of end-stage epithelial ovarian carcinoma. Immune cell profiling, impact on secondary lesion development and survival were evaluated in OrfV-treated mice as well as in Batf3 knockout, mice depleted of specific immune cell subsets and in mice where the primary tumor was removed. Finally, we interrogated gene expression datasets from primary human ovarian tumors from the International Cancer Genome Consortium database to determine whether the interplay we observed between natural killer (NK) cells, classical type 1 dendritic cells (cDC1s) and T cells exists and influences outcomes in human ovarian cancer.

RESULTS

OrfV was an effective monotherapy in a murine model of advanced-stage epithelial ovarian cancer. OrfV intervention relied on NK cells, which when depleted abrogated antitumor CD8⁺ T-cell responses. OrfV therapy was shown to require cDC1s in experiments with BATF3 knockout mice, which do not have mature cDC1s. Furthermore, cDC1s governed antitumor NK and T-cell responses to mediate antitumor efficacy following OrfV. Primary tumor removal, a common treatment option in human patients, was effectively combined with OrfV for optimal therapeutic outcome. Analysis of human RNA sequencing datasets revealed that cDC1s correlate with NK cells in human ovarian cancer and that intratumoral NK cells correlate positively with survival.

CONCLUSIONS

The data herein support the translational potential of OrfV as an NK stimulating immunotherapeutic for the treatment of advanced-stage ovarian cancer.

Apoplastic class III peroxidases PRX62 and PRX69 promote Arabidopsis root hair growth at low temperature

Root Hairs (RHs) growth is influenced by endogenous and by external environmental signals that coordinately regulate its final cell size. We have recently determined that RH growth was unexpectedly boosted when Arabidopsis thaliana seedlings are cultivated at low temperatures. It was proposed that RH growth plasticity in response to low temperature was linked to a reduced nutrient availability in the media. Here, we explore the molecular basis of this RH growth response by using a Genome Wide Association Study (GWAS) approach using Arabidopsis thaliana natural accessions. We identify the poorly characterized PEROXIDASE 62 (PRX62) and a related protein PRX69 as key proteins under moderate low temperature stress. Strikingly, a cell wall protein extensin (EXT) reporter reveals the effect of peroxidase activity on EXT cell wall association at 10 °C in the RH apical zone. Collectively, our results indicate that PRX62, and to a lesser extent PRX69, are key apoplastic PRXs that modulate ROS-homeostasis and cell wall EXT-insolubilization linked to RH elongation at low temperature.

Computational identification of host genomic biomarkers highlighting their functions, pathways and regulators that influence SARS-CoV-2 infections and drug repurposing

The pandemic threat of COVID-19 has severely destroyed human life as well as the economy around the world. Although, the vaccination has reduced the outspread, but people are still suffering due to the unstable RNA sequence patterns of SARS-CoV-2 which demands supplementary drugs. To explore novel drug target proteins, in this study, a transcriptomics RNA-Seq data generated from SARS-CoV-2 infection and control samples were analyzed. We identified 109 differentially expressed genes (DEGs) that were utilized to identify 10 hub-genes/proteins (TLR2, USP53, GUCY1A2, SNRPD2, NEDD9, IGF2, CXCL2, KLF6, PAG1 and ZFP36) by the protein–protein interaction (PPI) network analysis. The GO functional and KEGG pathway enrichment analyses of hub-DEGs revealed some important functions and signaling pathways that are significantly associated with SARS-CoV-2 infections. The interaction network analysis identified 5 TFs proteins and 6 miRNAs as the key regulators of hub-DEGs. Considering 10 hub-proteins and 5 key TFs-proteins as drug target receptors, we performed their docking analysis with the SARS-CoV-2 3CL protease-guided top listed 90 FDA approved drugs. We found Torin-2, Rapamycin, Radotinib, Ivermectin, Thiostrepton, Tacrolimus and Daclatasvir as the top ranked seven candidate drugs. We investigated their resistance performance against the already published COVID-19 causing top-ranked 11 independent and 8 protonated receptor proteins by molecular docking analysis and found their strong binding affinities, which indicates that the proposed drugs are effective against the state-of-the-arts alternatives independent receptor proteins also. Finally, we investigated the stability of top three drugs (Torin-2, Rapamycin and Radotinib) by using 100 ns MD-based MM-PBSA simulations with the two top-ranked proposed receptors (TLR2, USP53) and independent receptors (IRF7, STAT1), and observed their stable performance. Therefore, the proposed drugs might play a vital role for the treatment against different variants of SARS-CoV-2 infections.

A novel single-cell RNA-sequencing approach and its applicability connecting genotype to phenotype in ageing disease

Single cell multi-omics analysis has the potential to yield a comprehensive understanding of the cellular events that underlie the basis of human diseases. The cardinal feature to access this information is the technology used for single-cell isolation, barcoding, and sequencing. Most currently used single-cell RNA-sequencing platforms have limitations in several areas including cell selection, documentation and library chemistry. In this study, we describe a novel high-throughput, full-length, single-cell RNA-sequencing approach that combines the CellenONE isolation and sorting system with the ICELL8 processing instrument. This method offers substantial improvements in single cell selection, documentation and capturing rate. Moreover, it allows the use of flexible chemistry for library preparations and the analysis of living or fixed cells, whole cells independent of sizing and morphology, as well as of nuclei. We applied this method to dermal fibroblasts derived from six patients with different segmental progeria syndromes and defined phenotype associated pathway signatures with variant associated expression modifiers. These results validate the applicability of our method to highlight genotype-expression relationships for molecular phenotyping of individual cells derived from human patients.

Discovery of antimicrobials by massively parallelized growth assays (Mex)

The number of newly approved antimicrobial compounds has been steadily decreasing over the past 50 years emphasizing the need for novel antimicrobial substances. Here we present Me^x, a method for the high-throughput discovery of novel antimicrobials, that relies on E. coli self-screening to determine the bioactivity of more than ten thousand naturally occurring peptides. Analysis of thousands of E. coli growth curves using next-generation sequencing enables the identification of more than 1000 previously unknown antimicrobial peptides. Additionally, by incorporating the kinetics of growth inhibition, a first indication of the mode of action is obtained, which has implications for the ultimate usefulness of the peptides in question. The most promising peptides of the screen are chemically synthesized and their activity is determined in standardized susceptibility assays. Ten out of 15 investigated peptides efficiently eradicate bacteria at a minimal inhibitory concentration in the lower µM or upper nM range. This work represents a step-change in the high-throughput discovery of functionally diverse antimicrobials.

Deciphering the link between Diabetes mellitus and SARS-CoV-2 infection through differential targeting of microRNAs in the human pancreas

PURPOSE

Coronavirus Disease 2019 (COVID-19) severity and Diabetes mellitus affect each other bidirectionally. However, the cause of severe acute respiratory syndrome coronavirus 2 (SARS-CoV2) infection on the incidence of diabetes is unclear. In the SARS-CoV-2-infected cells, host microRNAs (miRNAs) may target the native gene transcripts as well as the viral genomic and subgenomic RNAs. Here, we investigated the role of miRNAs in linking Diabetes to SARS-CoV-2 infection in the human pancreas.

METHODS

Differential gene expression and disease enrichment analyses were performed on an RNA-Seq dataset of human embryonic stem cell-derived (hESC) mock-infected and SARS-CoV-2-infected pancreatic organoids to obtain the dysregulated Diabetes-associated genes. The miRNA target prediction for the Diabetes-associated gene transcripts and the SARS-CoV-2 RNAs has been made to determine the common miRNAs targeting them. Minimum Free Energy (MFE) analysis was done to identify the miRNAs, preferably targeting SARS-CoV-2 RNAs over the Diabetes-associated gene transcripts.

RESULTS

The gene expression and disease enrichment analyses of the RNA-Seq data have revealed five biomarker genes, i.e., CP, SOCS3, AGT, PSMB8 and CFB that are associated with Diabetes and get significantly upregulated in the pancreas following SARS-CoV-2-infection. Four miRNAs, i.e., hsa-miR-298, hsa-miR-3925-5p, hsa-miR-4691-3p and hsa-miR-5196-5p, showed preferential targeting of the SARS-CoV-2 genome over the cell's Diabetes-associated messenger RNAs (mRNAs) in the human pancreas.

CONCLUSION

Our study proposes that the differential targeting of the Diabetes-associated host genes by the miRNAs may lead to diabetic complications or new-onset Diabetes that can worsen the condition of COVID-19 patients.

Bacillus velezensis stimulates resident rhizosphere Pseudomonas stutzeri for plant health through metabolic interactions

Trophic interactions play a central role in driving microbial community assembly and function. In gut or soil ecosystems, successful inoculants are always facilitated by efficient colonization; however, the metabolite exchanges between inoculants and resident bacteria are rarely studied, particularly in the rhizosphere. Here, we used bioinformatic, genetic, transcriptomic, and metabonomic analyses to uncover syntrophic cooperation between inoculant (Bacillus velezensis SQR9) and plant-beneficial indigenous Pseudomonas stutzeri in the cucumber rhizosphere. We found that the synergistic interaction of these two species is highly environmental dependent, the emergence of syntrophic cooperation was only evident in a static nutrient-rich niche, such as pellicle biofilm in addition to the rhizosphere. Our results identified branched-chain amino acids (BCAAs) biosynthesis pathways are involved in syntrophic cooperation. Genome-scale metabolic modeling and metabolic profiling also demonstrated metabolic facilitation among the bacterial strains. In addition, biofilm matrix components from Bacillus were essential for the interaction. Importantly, the two-species consortium promoted plant growth and helped plants alleviate salt stress. In summary, we propose a mechanism in which synergic interactions between a biocontrol bacterium and a partner species promote plant health.

The novel activity of Argonautes in intron splicing: A transcriptome-wide survey in plants

The importance of the evolutionarily conserved Argonaute (AGO) proteins has been well recognized for their involvement in the RNA interference pathways. Recent discoveries in animals demonstrated that AGOs also participate in alternative splicing (AS). Motivated by the question whether the AGO proteins are also functional in RNA splicing in plants, we searched for the introns excised through an AGO-dependent manner in Arabidopsis (Arabidopsis thaliana). RNA sequencing (RNA-seq) data analysis uncovered hundreds of the introns up- or down-regulated in the ago1 and ago4 mutants, respectively. For different genes, AGOs might play either a positive or a negative role in intron excision, which was further validated by reverse transcription-polymerase chain reaction (RT-PCR). Some introns were specifically regulated by one of the AGO proteins, while some were regulated by both AGOs. Besides, a large portion of the AGO-dependent introns were organ-specifically regulated. RNA immunoprecipitation combined with high-throughput sequencing (RIP-seq) revealed that both AGOs preferentially bound to the intronic regions, supporting their high intron binding affinities. Immunoprecipitation followed by mass spectrometry (IP-MS) was performed to identify the proteins potentially interacting with the two AGOs. Six novel interactors (two interacting with AGO1 and four with both AGOs) involved in mRNA binding were uncovered, which might facilitate AGO-intron recognition. Analysis of the RNA-seq data from the rice (Oryza sativa) ago18 mutants revealed that hundreds of the introns were expressed in an AGO18-dependent manner. In summary, our results point to the novel role of the plant AGOs in intron splicing, paving a way for further studies on the mechanisms underlying AGO-mediated RNA splicing.

Analysis of Genomic Alternative Splicing Patterns in Rat under Heat Stress Based on RNA-Seq Data

Heat stress is one of the most severe challenges faced in livestock production in summer. Alternative splicing as an important post-transcriptional regulation is rarely studied in heat-stressed animals. Here, we performed and analyzed RNA-sequencing assays on the liver of Sprague-Dawley rats in control (22 °C, n = 5) and heat stress (4 °C for 120 min, H120; n = 5) groups, resulting in the identification of 636 differentially expressed genes. Identification analysis of the alternative splicing events revealed that heat stress-induced alternative splicing events increased by 20.18%. Compared with other types of alternative splicing events, the alternative start increased the most (43.40%) after heat stress. Twenty-eight genes were differentially alternatively spliced (DAS) between the control and H120 groups, among which Acly, Hnrnpd and mir3064 were also differentially expressed. For DAS genes, Srebf1, Shc1, Srsf5 and Ensa were associated with insulin, while Cast, Srebf1, Tmem33, Tor1aip2, Slc39a7 and Sqstm1 were enriched in the composition of the endoplasmic reticulum. In summary, our study conducts a comprehensive profile of alternative splicing in heat-stressed rats, indicating that alternative splicing is one of the molecular mechanisms of heat stress response in mammals and providing reference data for research on heat tolerance in mammalian livestock.

Differential gene expression and network analysis in head and neck squamous cell carcinoma

Head and neck squamous cell carcinoma (HNSCC) is a prevalent malignancy with a poor prognosis, whose biomarkers have not been studied in great detail. We have collected genomic data of HNSCC patients from The Cancer Genome Atlas (TCGA) and analyzed them to get deeper insights into the gene expression pattern. Initially, 793 differentially expressed genes (DEGs) were categorized, and their enrichment analysis was performed. Later, a protein-protein interaction network for the DEGs was constructed using the STRING plugin in Cytoscape to study their interactions. A set of 10 hub genes was selected based on Maximal Clique Centrality score, and later their survival analysis was studied. The elucidated set of 10 genes, i.e., PRAME, MAGEC2, MAGEA12, LHX1, MAGEA3, CSAG1, MAGEA6, LCE6A, LCE2D, LCE2C, referred to as potential candidates to be explored as HNSCC biomarkers. The Kaplan-Meier overall survival of the selected genes suggested that the alterations in the candidate genes were linked to the decreased survival of the HNSCC patients. Altogether, the results of this study signify that the genomic alterations and differential expression of the selected genes can be explored in therapeutic interpolations of HNSCC, exploiting early diagnosis and target-propelled therapy.

Beware to ignore the rare: how imputing zero-values can improve the quality of 16S rRNA gene studies results

BACKGROUND

16S rRNA-gene sequencing is a valuable approach to characterize the taxonomic content of the whole bacterial population inhabiting a metabolic and spatial niche, providing an important opportunity to study bacteria and their role in many health and environmental mechanisms. The analysis of data produced by amplicon sequencing, however, brings very specific methodological issues that need to be properly addressed to obtain reliable biological conclusions. Among these, 16S count data tend to be very sparse, with many null values reflecting species that are present but got unobserved due to the multiplexing constraints. However, current data workflows do not consider a step in which the information about unobserved species is recovered.

RESULTS

In this work, we evaluate for the first time the effects of introducing in the 16S data workflow a new preprocessing step, zero-imputation, to recover this lost information. Due to the lack of published zero-imputation methods specifically designed for 16S count data, we considered a set of zero-imputation strategies available for other frameworks, and benchmarked them using in silico 16S count data reflecting different experimental designs. Additionally, we assessed the effect of combining zero-imputation and normalization, i.e. the only preprocessing step in current 16S workflow. Overall, we benchmarked 35 16S preprocessing pipelines assessing their ability to handle data sparsity, identify species presence/absence, recovery sample proportional abundance distributions, and improve typical downstream analyses such as computation of alpha and beta diversity indices and differential abundance analysis.

CONCLUSIONS

The results clearly show that 16S data analysis greatly benefits from a properly-performed zero-imputation step, despite the choice of the right zero-imputation method having a pivotal role. In addition, we identify a set of best-performing pipelines that could be a valuable indication for data analysts.

Constitutive activation of canonical Wnt signaling disrupts choroid plexus epithelial fate

The choroid plexus secretes cerebrospinal fluid and is critical for the development and function of the brain. In the telencephalon, the choroid plexus epithelium arises from the Wnt- expressing cortical hem. Canonical Wnt signaling pathway molecules such as nuclear β-CATENIN are expressed in the mouse and human embryonic choroid plexus epithelium indicating that this pathway is active. Point mutations in human β-CATENIN are known to result in the constitutive activation of canonical Wnt signaling. In a mouse model that recapitulates this perturbation, we report a loss of choroid plexus epithelial identity and an apparent transformation of this tissue to a neuronal identity. Aspects of this phenomenon are recapitulated in human embryonic stem cell derived organoids. The choroid plexus is also disrupted when β-Catenin is conditionally inactivated. Together, our results indicate that canonical Wnt signaling is required in a precise and regulated manner for normal choroid plexus development in the mammalian brain.

Capture-C: a modular and flexible approach for high-resolution chromosome conformation capture

Chromosome conformation capture (3C) methods measure the spatial proximity between DNA elements in the cell nucleus. Many methods have been developed to sample 3C material, including the Capture-C family of protocols. Capture-C methods use oligonucleotides to enrich for interactions of interest from sequencing-ready 3C libraries. This approach is modular and has been adapted and optimized to work for sampling of disperse DNA elements (NuTi Capture-C), including from low cell inputs (LI Capture-C), as well as to generate Hi-C like maps for specific regions of interest (Tiled-C) and to interrogate multiway interactions (Tri-C). We present the design, experimental protocol and analysis pipeline for NuTi Capture-C in addition to the variations for generation of LI Capture-C, Tiled-C and Tri-C data. The entire procedure can be performed in 3 weeks and requires standard molecular biology skills and equipment, access to a next-generation sequencing platform, and basic bioinformatic skills. Implemented with other sequencing technologies, these methods can be used to identify regulatory interactions and to compare the structural organization of the genome in different cell types and genetic models.

A masculinizing supergene underlies an exaggerated male reproductive morph in a spider

In many species, individuals can develop into strikingly different morphs, which are determined by a simple Mendelian locus. How selection shapes loci that control complex phenotypic differences remains poorly understood. In the spider Oedothorax gibbosus, males develop either into a 'hunched' morph with conspicuous head structures or as a fast-developing 'flat' morph with a female-like appearance. We show that the hunched-determining allele contains a unique genomic fragment of approximately 3 megabases that is absent in the flat-determining allele. This fragment comprises dozens of genes that duplicated from genes found at the same as well as different chromosomes. All functional duplicates, including a duplicate of the key sexual differentiation regulatory gene doublesex, show male-specific expression, which illustrates their integrated role as a masculinizing supergene. Our findings demonstrate how extensive indel polymorphisms and duplications of regulatory genes may contribute to the evolution of co-adapted gene clusters, sex-limited reproductive morphs and the enigmatic evolution of exaggerated sexual traits in general.

Cohort profile: InTraUterine sampling in early pregnancy (ITU), a prospective pregnancy cohort study in Finland: study design and baseline characteristics

PURPOSE

The InTraUterine sampling in early pregnancy (ITU) is a prospective pregnancy cohort study. The overarching aim of ITU is to unravel genomic, epigenomic, transcriptomic, endocrine, inflammatory and metabolic maternal-placental-fetal mechanisms involved in the programming of health and disease after exposure to prenatal environmental adversity, such as maternal malnutrition, cardiometabolic disorders, infections, medical interventions, mental disorders and psychosocial stress. This paper describes the study protocol, design and baseline characteristics of the cohort.

PARTICIPANTS

We included 944 pregnant Finnish women, their partners and children born alive between April 2012 and December 2017. The women were recruited through the national, voluntary trisomy 21 screening between 9⁺⁰ and 21⁺⁶ gestational weeks. Of the participating women, 543 were screen positive and underwent fetal chromosomal testing. Test result of these women suggested no fetal chromosomal abnormality. Further, we recruited 401 women who were screen negative and who did not undergo fetal chromosomal testing.

FINDINGS TO DATE

We have collected chorionic villi and amniotic fluid from the screen-positive women; blood, urine, buccal swabs and diurnal salivary samples from all women; blood and buccal swabs from all partners; and placenta, cord blood and buccal swabs from all newborns for analyses of the genome, epigenome, transcriptome, and endocrine, inflammatory and metabolic markers. These data are coupled with comprehensive phenotypes, including questions on demographic characteristics, health and well-being of the women and their partners during pregnancy and of the women and their children at the child's age of 1.7 and 3 years. Data also come from patient records and nationwide registers covering health, lifestyle and medication data.

FUTURE PLANS

Multiple layers of ITU data allow integrative data analyses, which translate to biomarker identification and allow risk stratification and understanding of the biological mechanisms involved in prenatal programming of health and disease.

Pathogenesis-related protein 1 (PR-1) genes in soybean: Genome-wide identification, structural analysis and expression profiling under multiple biotic and abiotic stresses

Plant pathogenesis-related (PR) proteins are a large group of proteins, classified in 17 families, that are induced by pathological conditions. Here, we characterized the soybean PR-1 (GmPR-1) gene repertoire at the sequence, structural and expression levels. We found 24 GmPR-1 genes, clustered in two phylogenetic groups. GmPR-1 genes are under strong purifying selection, particularly those that emerged by tandem duplications. GmPR-1 promoter regions are abundant in cis-regulatory elements associated with major stress-related transcription factor families, namely WRKY, ERF, HD-Zip, C2H2, NAC, and GATA. We observed that 23 GmPR-1 genes are induced by stress conditions or exclusively expressed upon stress. We explored 1972 transcriptome samples, including 26 stress conditions, revealing that most GmPR-1 genes are differentially expressed in a plethora of biotic and abiotic stresses. Our findings highlight stress-responsive GmPR-1 genes with potential biotechnological applications, such as the development of transgenic lines with increased resistance to biotic and abiotic stresses.

Salt Stress Modulates the Landscape of Transcriptome and Alternative Splicing in Date Palm (Phoenix dactylifera L.)

Date palm regards as a valuable genomic resource for exploring the tolerance genes due to its ability to survive under the sever condition. Although a large number of differentiated genes were identified in date palm responding to salt stress, the genome-wide study of alternative splicing (AS) landscape under salt stress conditions remains unknown. In the current study, we identified the stress-related genes through transcriptomic analysis to characterize their function under salt. A total of 17,169 genes were differentially expressed under salt stress conditions. Gene expression analysis confirmed that the salt overly sensitive (SOS) pathway genes, such as PdSOS2;1, PdSOS2;2, PdSOS4, PdSOS5, and PdCIPK11 were involved in the regulation of salt response in date palm, which is consistent with the physiological analysis that high salinity affected the Na+/K+ homeostasis and amino acid profile of date palm resulted in the inhibition of plant growth. Interestingly, the pathway of "spliceosome" was enriched in the category of upregulation, indicating their potential role of AS in date palm response to salt stress. Expectedly, many differentially alternative splicing (DAS) events were found under salt stress conditions, and some splicing factors, such as PdRS40, PdRSZ21, PdSR45a, and PdU2Af genes were abnormally spliced under salt, suggesting that AS-related proteins might participated in regulating the salt stress pathway. Moreover, the number of differentially DAS-specific genes was gradually decreased, while the number of differentially expressed gene (DEG)-specific genes was increased with prolonged salt stress treatment, suggesting that AS and gene expression could be distinctively regulated in response to salt stress. Therefore, our study highlighted the pivotal role of AS in the regulation of salt stress and provided novel insights for enhancing the resistance to salt in date palm.

Identification of an integrated stress and growth response signaling switch that directs vertebrate intestinal regeneration

BACKGROUND

Snakes exhibit extreme intestinal regeneration following months-long fasts that involves unparalleled increases in metabolism, function, and tissue growth, but the specific molecular control of this process is unknown. Understanding the mechanisms that coordinate these regenerative phenotypes provides valuable opportunities to understand critical pathways that may control vertebrate regeneration and novel perspectives on vertebrate regenerative capacities.

RESULTS

Here, we integrate a comprehensive set of phenotypic, transcriptomic, proteomic, and phosphoproteomic data from boa constrictors to identify the mechanisms that orchestrate shifts in metabolism, nutrient uptake, and cellular stress to direct phases of the regenerative response. We identify specific temporal patterns of metabolic, stress response, and growth pathway activation that direct regeneration and provide evidence for multiple key central regulatory molecules kinases that integrate these signals, including major conserved pathways like mTOR signaling and the unfolded protein response.

CONCLUSION

Collectively, our results identify a novel switch-like role of stress responses in intestinal regeneration that forms a primary regulatory hub facilitating organ regeneration and could point to potential pathways to understand regenerative capacity in vertebrates.

Inferring the genetic responses to acute drought stress across an ecological gradient

BACKGROUND

How do xerophytic species thrive in environments that experience extreme annual drought? Although critical to the survival of many species, the genetic responses to drought stress in many non-model organisms has yet to be explored. We investigated this question in Mentzelia section Bartonia (Loasaceae), which occurs throughout western North America, including arid lands. To better understand the genetic responses to drought stress among species that occur in different habitats, the gene expression levels of three species from Mentzelia were compared across a precipitation gradient. Two de novo reference transcriptomes were generated and annotated. Leaf and root tissues were collected from control and drought shocked plants and compared to one another for differential expression. A target-gene approach was also implemented to better understand how drought-related genes from model and crop species function in non-model systems.

RESULTS

When comparing the drought-shock treatment plants to their respective control plants, we identified 165 differentially expressed clusters across all three species. Differentially expressed genes including those associated with water movement, photosynthesis, and delayed senescence. The transcriptome profiling approach was coupled with a target genes approach that measured expression of 90 genes associated with drought tolerance in model organisms. Comparing differentially expressed genes with a ≥ 2 log-fold value between species and tissue types showed significant differences in drought response. In pairwise comparisons, species that occurred in drier environments differentially expressed greater genes in leaves when drought shocked than those from wetter environments, but expression in the roots mostly produced opposite results.

CONCLUSIONS

Arid-adapted species mount greater genetic responses compared to the mesophytic species, which has likely evolved in response to consistent annual drought exposure across generations. Drought responses also depended on organ type. Xerophytes, for example, mounted a larger response in leaves to downregulate photosynthesis and senescence, while mobilizing carbon and regulating water in the roots. The complexity of drought responses in Mentzelia suggest that whole organism responses need to be considered when studying drought and, in particular, the physiological mechanisms in which plants regulate water, carbon, cell death, metabolism, and secondary metabolites.