EMLI-ICC: an ensemble machine learning-based integration algorithm for metastasis prediction and risk stratification in intrahepatic cholangiocarcinoma.

Robust strategies to identify patients at high risk for tumor metastasis, such as those frequently observed in intrahepatic cholangiocarcinoma (ICC), remain limited. While gene/protein expression profiling holds great potential as an approach to cancer diagnosis and prognosis, previously developed protocols using multiple diagnostic signatures for expression-based metastasis prediction have not been widely applied successfully because batch effects and different data types greatly decreased the predictive performance of gene/protein expression profile-based signatures in interlaboratory and data type dependent validation. To address this problem and assist in more precise diagnosis, we performed a genome-wide integrative proteome and transcriptome analysis and developed an ensemble machine learning-based integration algorithm for metastasis prediction (EMLI-Metastasis) and risk stratification (EMLI-Prognosis) in ICC. Based on massive proteome (216) and transcriptome (244) data sets, 132 feature (biomarker) genes were selected and used to train the EMLI-Metastasis algorithm. To accurately detect the metastasis of ICC patients, we developed a weighted ensemble machine learning method based on k-Top Scoring Pairs (k-TSP) method. This approach generates a metastasis classifier for each bootstrap aggregating training data set. Ten binary expression rank-based classifiers were generated for detection of metastasis separately. To further improve the accuracy of the method, the 10 binary metastasis classifiers were combined by weighted voting based on the score from the prediction results of each classifier. The prediction accuracy of the EMLI-Metastasis algorithm achieved 97.1% and 85.0% in proteome and transcriptome datasets, respectively. Among the 132 feature genes, 21 gene-pair signatures were developed to establish a metastasis-related prognosis risk-stratification model in ICC (EMLI-Prognosis). Based on EMLI-Prognosis algorithm, patients in the high-risk group had significantly dismal overall survival relative to the low-risk group in the clinical cohort (P-value < 0.05). Taken together, the EMLI-ICC algorithm provides a powerful and robust means for accurate metastasis prediction and risk stratification across proteome and transcriptome data types that is superior to currently used clinicopathological features in patients with ICC. Our developed algorithm could have profound implications not just in improved clinical care in cancer metastasis risk prediction, but also more broadly in machine-learning-based multi-cohort diagnosis method development. To make the EMLI-ICC algorithm easily accessible for clinical application, we established a web-based server for metastasis risk prediction (http://ibi.zju.edu.cn/EMLI/).

Are miRNAs applicable for balancing crop growth and defense trade-off?

Securing agricultural supplies for the increasing population without negative impacts on environment demands new crop varieties with higher yields, better quality, and stronger stress resilience. But breeding such super crop varieties is restrained by growth-defense (G-D) trade-off. MicroRNAs (miRNAs) are versatile regulators of plant growth and immune responses, with several being demonstrated to simultaneously regulate crop growth and defense against biotic stresses and to balance G-D trade-off. Increasing evidence also links miRNAs to the metabolism and signaling of phytohormones, another type of master regulator of plant growth and defense. Here, we synthesize the reported functions of miRNAs in crop growth, development, and responses to bio-stressors, summarize the regulatory scenarios of miRNAs based on their relationship with target(s), and discuss how miRNAs, particularly those involved in crosstalk with phytohormones, can be applied in balancing G-D trade-off in crops. We also propose several open questions to be addressed for adopting miRNAs in balancing crop G-D trade-off.

MicroRNAs constitute an additional layer in plant response to simultaneous bio- and abiotic stresses as exemplified by UV-B radiation and flg22-treatment on Arabidopsis thaliana.

Plants are confronted with various environmental stresses and develop sophisticated adaptive mechanisms. Our previous work demonstrated that the crosstalk of flg22 and ultraviolet (UV)-B-induced signalling cascades reprograms the expression of flavonol pathway genes (FPGs), benefiting plant defence responses. Although several transcription factors have been identified to be involved in this crosstalk, the underlying mechanism is largely unclear. Here, we analyzed microRNAs (miRNAs) and identified 126, 129 and 113 miRNAs with altered abundances compared to untreated control in flg22-, UV-B- and flg22/UV-B-treated seedlings, respectively. Two distinct modules were identified: The first consists of 10 miRNAs repressed by UV-B but up-regulated by flg22, and the second with five miRNAs repressed by flg22 but up-regulated by UV-B. In Arabidopsis, the knockdown of miR858a, a representative of module I, increased the abundance of CHS (a marker gene for FPGs), whereas its overexpression reduced CHS. Conversely, knockout of miR164b from module II decreased CHS and its overexpression increased CHS transcript levels. These data suggest a decisive role of miRNAs in the crosstalk. In the next, we described the interaction between miR858a and its target MYB111 (a positive regulator of FPGs) from module I in detail. We showed that MYB111 was profoundly post-transcriptionally regulated by miR858a during the crosstalk, whose expression was specifically but antagonistically controlled by UVR8- and FLS2-mediated signallings. Moreover, transcriptional monitoring using the GUS reporter gene demonstrates that miRNA-mediated posttranscriptional regulation is the main driving force in reprogramming the expression of FPGs and regulates plant adaptation to multiple concurrent environmental stresses.

The Conservation of Long Intergenic Non-Coding RNAs and Their Response to Verticillium dahliae Infection in Cotton.

Long intergenic non-coding RNAs (lincRNAs) have been demonstrated to be vital regulators of diverse biological processes in both animals and plants. While many lincRNAs have been identified in cotton, we still know little about the repositories and conservativeness of lincRNAs in different cotton species or about their role in responding to biotic stresses. Here, by using publicly available RNA-seq datasets from diverse sources, including experiments of Verticillium dahliae (Vd) infection, we identified 24,425 and 17,713 lincRNAs, respectively, in Gossypium hirsutum (Ghr) and G. barbadense (Gba), the two cultivated allotetraploid cotton species, and 6933 and 5911 lincRNAs, respectively, in G. arboreum (Gar) and G. raimondii (Gra), the two extant diploid progenitors of the allotetraploid cotton. While closely related subgenomes, such as Ghr_At and Gba_At, tend to have more conserved lincRNAs, most lincRNAs are species-specific. The majority of the synthetic and transcribed lincRNAs (78.2%) have a one-to-one orthologous relationship between different (sub)genomes, although a few of them (0.7%) are retained in all (sub)genomes of the four species. The Vd responsiveness of lincRNAs seems to be positively associated with their conservation level. The major functionalities of the Vd-responsive lincRNAs seem to be largely conserved amongst Gra, Ghr, and Gba. Many Vd-responsive Ghr-lincRNAs overlap with Vd-responsive QTL, and several lincRNAs were predicted to be endogenous target mimicries of miR482/2118, with a pair being highly conserved between Ghr and Gba. On top of the confirmation of the feature characteristics of the lincRNAs previously reported in cotton and other species, our study provided new insights into the conservativeness and divergence of lincRNAs during cotton evolution and into the relationship between the conservativeness and Vd responsiveness of lincRNAs. The study also identified candidate lincRNAs with a potential role in disease response for functional characterization.

Expansion of MIR482/2118 by a class-II transposable element in cotton.

Some plant microRNA (miRNA) families contain multiple members generating identical or highly similar mature miRNA variants. Mechanisms underlying the expansion of miRNA families remain elusive, although tandem and/or segmental duplications have been proposed. In this study of two tetraploid cottons, Gossypium hirsutum and Gossypium barbadense, and their extant diploid progenitors, Gossypium arboreum and Gossypium raimondii, we investigated the gain and loss of members of the miR482/2118 superfamily, which modulates the expression of nucleotide-binding site leucine-rich repeat (NBS-LRR) disease resistance genes. We found significant expansion of MIR482/2118d in G. barbadense, G. hirsutum and G. raimondii, but not in G. arboreum. Several newly expanded MIR482/2118d loci have mutated to produce different miR482/2118 variants with altered target-gene specificity. Based on detailed analysis of sequences flanking these MIR482/2118 loci, we found that this expansion of MIR482/2118d and its derivatives resulted from an initial capture of an MIR482/2118d by a class-II DNA transposable element (TE) in G. raimondii prior to the tetraploidization event, followed by transposition to new genomic locations in G. barbadense, G. hirsutum and G. raimondii. The 'GosTE' involved in the capture and proliferation of MIR482/2118d and its derivatives belongs to the PIF/Harbinger superfamily, generating a 3-bp target site duplication upon insertion at new locations. All orthologous MIR482/2118 loci in the two diploids were retained in the two tetraploids, but mutation(s) in miR482/2118 were observed across all four species as well as in different cultivars of both G. barbadense and G. hirsutum, suggesting a dynamic co-evolution of miR482/2118 and its NBS-LRR targets. Our results provide fresh insights into the mechanisms contributing to MIRNA proliferation and enrich our knowledge on TEs.

miR159 Represses a Constitutive Pathogen Defense Response in Tobacco.

MicroR159 (miR159) regulation of GAMYB expression is highly conserved in terrestrial plants; however, its functional role remains poorly understood. In Arabidopsis (Arabidopsis thaliana), although GAMYB-like genes are constitutively transcribed during vegetative growth, their effects are suppressed by strong and constitutive silencing by miR159. GAMYB expression occurs only if miR159 function is inhibited, which results in detrimental pleiotropic defects, questioning the purpose of the miR159-GAMYB pathway. Here, miR159 function was inhibited in tobacco (Nicotiana tabacum) and rice (Oryza sativa) using miRNA MIM159 technology. Similar to observations in Arabidopsis, inhibition of miR159 in tobacco and rice resulted in pleiotropic defects including stunted growth, implying functional conservation of the miR159-GAMYB pathway among angiosperms. In MIM159 tobacco, transcriptome profiling revealed that genes associated with defense and programmed cell death were strongly activated, including a suite of 22 PATHOGENESIS-RELATED PROTEIN (PR) genes that were 100- to 1,000-fold upregulated. Constitutive expression of a miR159-resistant GAMYB transgene in tobacco resulted in phenotypes similar to that of MIM159 tobacco and activated PR gene expression, verifying the dependence of the above-mentioned changes on GAMYB expression. Consistent with the broad defense response, MIM159 tobacco appeared immune to Phytophthora infection. These findings suggest that the tobacco miR159-GAMYB pathway functions in the biotic defense response, which becomes activated upon miR159 inhibition. However, PR gene expression was not upregulated in Arabidopsis or rice when miR159 was inhibited, suggesting that miR159-GAMYB pathway functional differences exist between species, or factors in addition to miR159 inhibition are required in Arabidopsis and rice to activate this broad defense response.

SRSF6-regulated alternative splicing that promotes tumour progression offers a therapy target for colorectal cancer.

OBJECTIVE: To investigate the molecular function of splicing factor SRSF6 in colorectal cancer (CRC) progression and discover candidate chemicals for cancer therapy through targeting SRSF6. DESIGN: We performed comprehensive analysis for the expression of SRSF6 in 311 CRC samples, The Cancer Genome Atlas and Gene Expression Omnibus (GEO) database. Functional analysis of SRSF6 in CRC was performed in vitro and in vivo. SRSF6-regulated alternative splicing (AS) and its binding motif were identified by next-generation RNA-sequencing and RNA immunoprecipitation sequencing (RIP-seq), which was validated by gel shift and minigene reporter assay. ZO-1 exon23 AS was investigated to mediate the function of SRSF6 in vitro and in vivo. Based on the analysis of domain-specific role, SRSF6-targeted inhibitor was discovered de novoby virtual screening in 4855 FDA-approved drugs and its antitumour effects were evaluated in vitroand in vivo. RESULTS: SRSF6 was frequently upregulated in CRC samples and associated with poor prognosis, which promoted proliferation and metastasis in vitro and in vivo. We identified SRSF6-regulated AS targets and discovered the SRSF6 binding motif. Particularly, SRSF6 regulates ZO-1 aberrant splicing to function as an oncogene by binding directly to its motif in the exon23. Based on the result that SRSF6 RRM2 domain plays key roles in regulating AS and biological function, indacaterol, a ß2-adrenergic receptor agonist approved for chronic obstructive pulmonary disease treatment, is identified as the inhibitor of SRSF6 to suppress CRC tumourigenicity. CONCLUSIONS: SRSF6 functions the important roles in mediating CRC progression through regulating AS, and indacaterol is repositioned as an antitumour drug through targeting SRSF6. ACCESSION NUMBERS: The accession numbers for sequencing data are SRP111763 and SRP111797.

A transcriptomic profile of topping responsive non-coding RNAs in tobacco roots (Nicotiana tabacum).

BACKGROUND: Non-coding RNAs (ncRNAs), including microRNAs (miRNAs), long ncRNAs (lncRNAs) and circular RNAs (circRNAs), accomplish remarkable variety of biological functions. However, the composition of ncRNAs and their interactions with coding RNAs in modulating and controlling of cellular process in plants is largely unknown. Using a diverse group of high-throughput sequencing strategies, the mRNA, miRNA, lncRNA and circRNA compositions of tobacco (Nicotiana tabacum) roots determined and their alteration and potential biological functions in response to topping treatment analyzed. RESULTS: A total of 688 miRNAs, 7423 non-redundant lncRNAs and 12,414 circRNAs were identified, among which, some selected differentially expressed RNAs were verified by quantitative real-time PCR. Using the differentially expressed RNAs, a co-expression network was established that included all four types of RNAs. The number of circRNAs identified were higher than that of miRNAs and lncRNAs, but only two circRNAs were present in the co-expression network. LncRNAs appear to be the most active ncRNAs based on their numbers presented in the co-expression network, but none of them seems to be an eTM (endogenous Target Mimicry) of miRNAs. Integrated with analyses of sequence interaction, several mRNA-circRNA-miRNA interaction networks with a potential role in the regulation of nicotine biosynthesis were uncovered, including a QS-circQS-miR6024 interaction network. In this network miR6024 was significantly down-regulated, while the expression levels of its two targets, circQS and its host gene QS, were sharply increased following the topping treatment. CONCLUSIONS: These results illustrated the transcriptomic profiles of tobacco roots, the organ responsible for nicotine biosynthesis. mRNAs always play the most important roles, while ncRNAs are also expressed extensively for topping treatment response, especially circRNAs are the most activated in the ncRNA pool. These studies also provided insights on the coordinated regulation module of coding and non-coding RNAs in a single plant biological sample. The findings reported here indicate that ncRNAs appear to form interaction complex for the regulation of stress response forming regulation networks with transcripts involved in nicotine biosynthesis in tobacco.

Genome-wide selection footprints and deleterious variations in young Asian allotetraploid rapeseed.

Brassica napus (AACC, 2n = 38) is an important oilseed crop grown worldwide. However, little is known about the population evolution of this species, the genomic difference between its major genetic groups, such as European and Asian rapeseed, and the impacts of historical large-scale introgression events on this young tetraploid. In this study, we reported the de novo assembly of the genome sequences of an Asian rapeseed (B. napus), Ningyou 7, and its four progenitors and compared these genomes with other available genomic data from diverse European and Asian cultivars. Our results showed that Asian rapeseed originally derived from European rapeseed but subsequently significantly diverged, with rapid genome differentiation after hybridization and intensive local selective breeding. The first historical introgression of B. rapa dramatically broadened the allelic pool but decreased the deleterious variations of Asian rapeseed. The second historical introgression of the double-low traits of European rapeseed (canola) has reshaped Asian rapeseed into two groups (double-low and double-high), accompanied by an increase in genetic load in the double-low group. This study demonstrates distinctive genomic footprints and deleterious SNP (single nucleotide polymorphism) variants for local adaptation by recent intra- and interspecies introgression events and provides novel insights for understanding the rapid genome evolution of a young allopolyploid crop.

Analysis of transcriptional and epigenetic changes in hybrid vigor of allopolyploid Brassica napus uncovers key roles for small RNAs.

Heterosis is a fundamental biological phenomenon characterized by the superior performance of a hybrid compared with its parents. The underlying molecular basis for heterosis, particularly for allopolyploids, remains elusive. In this study we analyzed the transcriptomes of Brassica napus parental lines and their F1 hybrids at three stages of early flower development. Phenotypically, the F1 hybrids show remarkable heterosis in silique number and grain yield. Transcriptome analysis revealed that various phytohormone (auxin and salicylic acid) response genes are significantly altered in the F1 hybrids relative to the parental lines. We also found evidence for decreased expression divergence of the homoeologous gene pairs in the allopolyploid F1 hybrids and suggest that high-parental expression-level dominance plays an important role in heterosis. Small RNA and methylation studies aimed at examining the epigenetic effect of the changes in gene expression level in the F1 hybrids showed that the majority of the small interfering RNA (siRNA) clusters had a higher expression level in the F1 hybrids than in the parents, and that there was an increase in genome-wide DNA methylation in the F1 hybrid. Transposable elements associated with siRNA clusters had a higher level of methylation and a lower expression level in the F1 hybrid, implying that the non-additively expressed siRNA clusters resulted in lower activity of the transposable elements through DNA methylation in the hybrid. Our data provide insights into the role that changes in gene expression pattern and epigenetic mechanisms contribute to heterosis during early flower development in allopolyploid B. napus.

Genome-wide identification of oil biosynthesis-related long non-coding RNAs in allopolyploid Brassica napus.

BACKGROUND: Long noncoding RNAs (lncRNAs) are transcripts longer than 200 bp that do not encode proteins but nonetheless have been shown to play important roles in various biological processes in plants. Brassica napus is an important seed oil crop worldwide and the target of many genetic improvement activities. To understand better the function of lncRNAs in regulating plant metabolic activities, we carried out a genome-wide lncRNA identification of lncRNAs in Brassica napus with a focus on lncRNAs involved in lipid metabolism. Twenty ribosomal RNA depleted strand specific RNA-seq (ssRNA-seq) datasets were generatred using RNAs isolated from B. napus seeds at four developmental stages. For comparison we also included 30 publically available RNA-seq datasets generated from poly(A) enriched mRNAs isolated from from various Brassica napus tissues in our analysis. RESULTS: A total of 8905 lncRNA loci were identified, including 7100 long intergenic noncoding RNA (lincRNA) loci and 1805 loci generating long noncoding natural antisense transcript (lncNAT). Many lncRNAs were identified only in the ssRNA-seq and poly(A) RNA-seq dataset, suggesting that B. napus has a large lncRNA repertoire and it is necessary to use libraries prepared from different tissues and developmental stages as well as different library preparation approaches to capture the whole spectrum of lncRNAs. Analysis of coexpression networks revealed that among the regulatory modules are networks containing lncRNAs and protein-coding genes related to oil biosynthesis indicating a possible role of lncRNAs in the control of lipid metabolism. One such example is that several lncRNAs are potential regulators of BnaC08g11970D that encodes oleosin1, a protein found in oil bodies and involved in seed lipid accumulation. We also observed that the expression levels of B. napus lncRNAs is positively correlated with their conservation levels. CONCLUSIONS: We demonstrated that the B. napus genome has a large number of lncRNA and that these lncRNAs are expressed broadly across many developmental times and in different tissue types. We also provide evidence indicating that specific lncRNAs appear to be important regulators of lipid biosynthesis forming regulatory networks with transcripts involved in lipid biosynthesis. We also provide evidence that these lncRNAs are conserved in other species of the Brassicaceae family.

Comparative transcriptome profiling of two Brassica napus cultivars under chromium toxicity and its alleviation by reduced glutathione.

BACKGROUND: Chromium (Cr) being multifarious industrial used element, is considered a potential environmental threat. Cr found to be a prospective water and soil pollutant, and thus it is a current area of concern. Oilseed rape (Brassica napus L.) is well known as a major source of edible oil around the globe. Due to its higher growth, larger biomass and capability to uptake toxic materials B. napus is considered a potential candidate plant against unfavorable conditions. To date, no study has been done that described the Cr and GSH mechanism at RNA-Seq level. RESULTS: Both digital gene expression (DGE) and transcriptome profile analysis (TPA) approaches had opened new insights to uncover the several number of genes related to Cr stress and GSH alleviating mechanism in two leading cultivars (ZS 758 and Zheda 622) of B. napus plants. Data showed that Cr inhibited KEGG pathways i.e. stilbenoid, diarlyheptanoid and gingerol biosynthesis; limonene and pentose degradation and glutathione metabolism in ZS 758; and ribosome and glucosinolate biosynthesis in Zheda-622. On the other hand, vitamin B6, tryptophan, sulfur, nitrogen and fructose and manose metabolisms were induced in ZS 758, and zeatin biosynthesis, linoleic acid metabolism, arginine and proline metabolism, and alanine, asparate and glutamate metabolism pathways in Zheda 622. Cr increased the TFs that were related to hydralase activity, antioxidant activity, catalytic activity phosphatase and pyrophosphatase activity in ZS 758, and vitamin binding and oxidoreductase activity in Zheda 622. Cr also up-regulated the promising proteins related to intracellular membrane bounded organelles, nitrile hyrdatase activity, cytoskeleton protein binding and stress response. It also uncovered, a novel Cr-responsive protein (CL2535.Contig1_All) that was statistically increased as compared to control and GSH treated plants. Exogenously applied GSH successfully not only recovered the changes in metabolic pathways but also induced cysteine and methionine metabolism in ZS 758 and ubiquinone and other terpenoid-quinone biosynthesis pathways in Zheda 622. Furthermore, GSH increased the level of TFs i.e. the gene expression of antioxidant and catalytic activities, iron ion binding and hydrolase activity as compared with Cr. Moreover, results pointed out a novel GSH responsive protein (CL827.Contig3_All) whose expression was found to be significantly increased when compared than Cr stress. Results further delineated that GSH induced TFs such as glutathione disulphide oxidoreducatse and aminoacyl-tRNA ligase activity, and beta glucosidase activity in ZS 758. Similarly in Zheda 622, GSH induced the TFs for instance DNA binding and protein dimerization activity. GSH also highlighted the proteins that were involved in transportation, photosynthesis process, RNA polymerase activity, and against the metal toxicity. These results indicated that cultivar ZS 758 had better metabolism and showed higher tolerance against Cr toxicity. CONCLUSION: The responses of ZS 758 and Zheda 622 differed considerably at both physiological and transcriptional level. Moreover, RNA-Seq method explored the hazardous behavior of Cr as well as GSH up-regulating mechanism by activating plant metabolism, stress responsive genes, TFs and protein encyclopedia.

Regulation of Nicotine Biosynthesis by an Endogenous Target Mimicry of MicroRNA in Tobacco.

The interaction between noncoding endogenous target mimicry (eTM) and its corresponding microRNA (miRNA) is a newly discovered regulatory mechanism and plays pivotal roles in various biological processes in plants. Tobacco (Nicotiana tabacum) is a model plant for studying secondary metabolite alkaloids, of which nicotine accounts for approximately 90%. In this work, we identified four unique tobacco-specific miRNAs that were predicted to target key genes of the nicotine biosynthesis and catabolism pathways and an eTM, novel tobacco miRNA (nta)-eTMX27, for nta-miRX27 that targets QUINOLINATE PHOSPHORIBOSYLTRANSFERASE2 (QPT2) encoding a quinolinate phosphoribosyltransferase. The expression level of nta-miRX27 was significantly down-regulated, while that of QPT2 and nta-eTMX27 was significantly up-regulated after topping, and consequently, nicotine content increased in the topping-treated plants. The topping-induced down-regulation of nta-miRX27 and up-regulation of QPT2 were only observed in plants with a functional nta-eTMX27 but not in transgenic plants containing an RNA interference construct targeting nta-eTMX27. Our results demonstrated that enhanced nicotine biosynthesis in the topping-treated tobacco plants is achieved by nta-eTMX27-mediated inhibition of the expression and functions of nta-miRX27. To our knowledge, this is the first report about regulation of secondary metabolite biosynthesis by an miRNA-eTM regulatory module in plants.

Identification, evolution, and expression partitioning of miRNAs in allopolyploid Brassica napus.

The recently published genome of Brassica napus offers for the first time the opportunity to gain insights into the genomic organization and the evolution of miRNAs in oilseed rape. In this study, 12 small RNA libraries from two B. napus cultivars (Tapidor and Ningyou7) and their four double-haploid lines were sequenced, employing the newly sequenced B. napus genome, together with genomes of its progenitors Brassica rapa and Brassica oleracea. A total of 645 miRNAs including 280 conserved and 365 novel miRNAs were identified. Comparative analysis revealed a high level of genomic conservation of MIRNAs (75.9%) between the subgenomes of B. napus and its two progenitors' genomes, and MIRNA lost/gain events (133) occurred in B. napus after its speciation. Furthermore, significant partitioning of miRNA expressions between the two subgenomes in B. napus was detected. The data of degradome sequencing, miRNA-mediated cleavage, and expression analyses support specific interactions between miRNAs and their targets in the modulation of diverse physiological processes in roots and leaves, as well as in biosynthesis of, for example, glucosinolates and lipids in oilseed rape. These data provide a first genome-wide view on the origin, evolution, and genomic organization of B. napus MIRNAs.

High-throughput single-microbe RNA sequencing reveals adaptive state heterogeneity and host-phage activity associations in human gut microbiome.

Microbial communities such as those residing in the human gut are highly diverse and complex, and many with important implications in health and diseases. The effects and functions of these microbial communities are determined not only by their species compositions and diversities but also by the dynamic intra- and inter-cellular states at the transcriptional level. Powerful and scalable technologies capable of acquiring single-microbe-resolution RNA sequencing information in order to achieve comprehensive understanding of complex microbial communities together with their hosts is therefore utterly needed. Here we report the development and utilization of a droplet-based smRNA-seq (single-microbe RNA sequencing) method capable of identifying large species varieties in human samples, which we name smRandom-seq2. Together with a triple-module computational pipeline designed for the bacteria and bacteriophage sequencing data by smRandom-seq2 in four human gut samples, we established a single-cell level bacterial transcriptional landscape of human gut microbiome, which included 29,742 single microbes and 329 unique species. Distinct adaptive responses states among species in Prevotella and Roseburia genus and intrinsic adaptive strategy heterogeneity in Phascolarctobacterium succinatutens were uncovered. Additionally, we identified hundreds of novel host-phage transcriptional activity associations in the human gut microbiome. Our results indicated the smRandom-seq2 is a high-throughput and high-resolution smRNA-seq technique that is highly adaptable to complex microbial communities in real-word situations and promises new perspectives in the understanding of human microbiomes.

A syntelog-based pan-genome provides insights into rice domestication and de-domestication.

BACKGROUND: Asian rice is one of the world's most widely cultivated crops. Large-scale resequencing analyses have been undertaken to explore the domestication and de-domestication genomic history of Asian rice, but the evolution of rice is still under debate. RESULTS: Here, we construct a syntelog-based rice pan-genome by integrating and merging 74 high-accuracy genomes based on long-read sequencing, encompassing all ecotypes and taxa of Oryza sativa and Oryza rufipogon. Analyses of syntelog groups illustrate subspecies divergence in gene presence-and-absence and haplotype composition and identify massive genomic regions putatively introgressed from ancient Geng/japonica to ancient Xian/indica or its wild ancestor, including almost all well-known domestication genes and a 4.5-Mbp centromere-spanning block, supporting a single domestication event in main rice subspecies. Genomic comparisons between weedy and cultivated rice highlight the contribution from wild introgression to the emergence of de-domestication syndromes in weedy rice. CONCLUSIONS: This work highlights the significance of inter-taxa introgression in shaping diversification and divergence in rice evolution and provides an exploratory attempt by utilizing the advantages of pan-genomes in evolutionary studies.

Dynamic associations between the respiratory tract and gut antibiotic resistome of patients with COVID-19 and its prediction power for disease severity.

The antibiotic resistome is the collection of all antibiotic resistance genes (ARGs) present in an individual. Whether an individual's susceptibility to infection and the eventual severity of coronavirus disease 2019 (COVID-19) is influenced by their respiratory tract antibiotic resistome is unknown. Additionally, whether a relationship exists between the respiratory tract and gut ARGs composition has not been fully explored. We recruited 66 patients with COVID-19 at three disease stages (admission, progression, and recovery) and conducted a metagenome sequencing analysis of 143 sputum and 97 fecal samples obtained from them. Respiratory tract, gut metagenomes, and peripheral blood mononuclear cell (PBMC) transcriptomes are analyzed to compare the gut and respiratory tract ARGs of intensive care unit (ICU) and non-ICU (nICU) patients and determine relationships between ARGs and immune response. Among the respiratory tract ARGs, we found that Aminoglycoside, Multidrug, and Vancomycin are increased in ICU patients compared with nICU patients. In the gut, we found that Multidrug, Vancomycin, and Fosmidomycin were increased in ICU patients. We discovered that the relative abundances of Multidrug were significantly correlated with clinical indices, and there was a significantly positive correlation between ARGs and microbiota in the respiratory tract and gut. We found that immune-related pathways in PBMC were enhanced, and they were correlated with Multidrug, Vancomycin, and Tetracycline ARGs. Based on the ARG types, we built a respiratory tract-gut ARG combined random-forest classifier to distinguish ICU COVID-19 patients from nICU patients with an AUC of 0.969. Cumulatively, our findings provide some of the first insights into the dynamic alterations of respiratory tract and gut antibiotic resistome in the progression of COVID-19 and disease severity. They also provide a better understanding of how this disease affects different cohorts of patients. As such, these findings should contribute to better diagnosis and treatment scenarios.

The complete chloroplast genome of weedy rye Secale cereale subsp. segetale.

Weedy rye (Secale cereale subsp. segetale Zhukov 1928) is a problematic weed species in wheat field. However, it can potentially provide valuable genetics resources to increase the genetic variations and introduce desirable genes for rye and wheat breeding. Here, we assembled the complete chloroplast genome of S. cereale subsp. segetale. The chloroplast genome is 137,051 bp in length, containing a large single copy region (81,090 bp), a small single copy region (12,795 bp) and two separated inverted repeat regions (21,583 bp). A total of 131 unique genes were annotated, consisting of 82 protein-coding genes, 41 tRNA genes, and 8 rRNA genes. The phylogenetic analysis showed that Secale cereale subsp. segetale (weedy rye) and S. cereale subsp. cereale (rye) clustered together as sisters to other Secale species.

The complete chloroplast genome of Nicotiana plumbaginifolia.

Nicotiana plumbaginifolia Vivianiis 1802 is an annual herb, native to Mexico and South America. It is one of the most widely distributed tobacco species. As a wild tobacco, N. plumbaginifolia has provided several economically important disease-resistance genes to cultivated tobacco. We assembled the complete chloroplast genome of N. plumbaginifolia. The chloroplast genome is 155,945 bp in length, which includes a large single copy region (86,621 bp), a small single copy region (18,528 bp) and two separated inverted repeat regions (25,398 bp). A total of 117 unique genes were annotated, consisting of 84 protein-coding genes, 29 tRNA genes and 4 rRNA genes. Based on chloroplast genomes of 17 Nicotiana species, phylogenetic analyses indicated that N. plumbaginifolia was closely related to N. suaveolens and N. amplexicaulis.

Distinctive signatures of pathogenic and antibiotic resistant potentials in the hadal microbiome.

BACKGROUND: Hadal zone of the deep-sea trenches accommodates microbial life under extreme energy limitations and environmental conditions, such as low temperature, high pressure, and low organic matter down to 11,000 m below sea level. However, microbial pathogenicity, resistance, and adaptation therein remain unknown. Here we used culture-independent metagenomic approaches to explore the virulence and antibiotic resistance in the hadal microbiota of the Mariana Trench. RESULTS: The results indicate that the 10,898 m Challenger Deep bottom sediment harbored prosperous microbiota with contrasting signatures of virulence factors and antibiotic resistance, compared with the neighboring but shallower 6038 m steep wall site and the more nearshore 5856 m Pacific basin site. Virulence genes including several famous large translocating virulence genes (e.g., botulinum neurotoxins, tetanus neurotoxin, and Clostridium difficile toxins) were uniquely detected in the trench bottom. However, the shallower and more nearshore site sediment had a higher abundance and richer diversity of known antibiotic resistance genes (ARGs), especially for those clinically relevant ones (e.g., fosX, sul1, and TEM-family extended-spectrum beta-lactamases), revealing resistance selection under anthropogenic stresses. Further analysis of mobilome (i.e., the collection of mobile genetic elements, MGEs) suggests horizontal gene transfer mediated by phage and integrase as the major mechanism for the evolution of Mariana Trench sediment bacteria. Notably, contig-level co-occurring and taxonomic analysis shows emerging evidence for substantial co-selection of virulence genes and ARGs in taxonomically diverse bacteria in the hadal sediment, especially for the Challenger Deep bottom where mobilized ARGs and virulence genes are favorably enriched in largely unexplored bacteria. CONCLUSIONS: This study reports the landscape of virulence factors, antibiotic resistome, and mobilome in the sediment and seawater microbiota residing hadal environment of the deepest ocean bottom on earth. Our work unravels the contrasting and unique features of virulence genes, ARGs, and MGEs in the Mariana Trench bottom, providing new insights into the eco-environmental and biological processes underlying microbial pathogenicity, resistance, and adaptative evolution in the hadal environment.