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.

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.

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.

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.

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/).

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.

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.

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.

Genomic insights into the evolution of Echinochloa species as weed and orphan crop.

As one of the great survivors of the plant kingdom, barnyard grasses (Echinochloa spp.) are the most noxious and common weeds in paddy ecosystems. Meanwhile, at least two Echinochloa species have been domesticated and cultivated as millets. In order to better understand the genomic forces driving the evolution of Echinochloa species toward weed and crop characteristics, we assemble genomes of three Echinochloa species (allohexaploid E. crus-galli and E. colona, and allotetraploid E. oryzicola) and re-sequence 737 accessions of barnyard grasses and millets from 16 rice-producing countries. Phylogenomic and comparative genomic analyses reveal the complex and reticulate evolution in the speciation of Echinochloa polyploids and provide evidence of constrained disease-related gene copy numbers in Echinochloa. A population-level investigation uncovers deep population differentiation for local adaptation, multiple target-site herbicide resistance mutations of barnyard grasses, and limited domestication of barnyard millets. Our results provide genomic insights into the dual roles of Echinochloa species as weeds and crops as well as essential resources for studying plant polyploidization, adaptation, precision weed control and millet improvements.

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.

Population genomic analysis reveals domestication of cultivated rye from weedy rye.

Rye (Secale cereale) is an important crop with multiple uses and a valuable genetic resource for wheat breeding. However, due to its complex genome and outcrossing nature, the origin of cultivated rye remains elusive. The geneticist N.I. Vavilov proposed that cultivated rye had been domesticated from weedy rye, rather than directly from wild species like other crops. Unraveling the domestication history of rye will extend our understanding of crop evolution and upend our inherent understanding of agricultural weeds. To this end, in this study we generated the 8.5 Tb of whole-genome resequencing data from 116 worldwide accessions of wild, weedy, and cultivated rye, and demonstrated that cultivated rye was domesticated directly from weedy relatives with a similar but enhanced genomic selection by humans. We found that a repertoire of genes that experienced artificial selection is associated with important agronomic traits, including shattering, grain yield, and disease resistance. Furthermore, we identified a composite introgression in cultivated rye from the wild perennial Secale strictum and detected a 2-Mb introgressed fragment containing a candidate ammonium transporter gene with potential effect on the grain yield and plant growth of rye. Taken together, our findings unravel the domestication history of cultivated rye, suggest that interspecific introgression serves as one of the likely causes of obscure species taxonomy of the genus Secale, and provide an important resource for future rye and wheat breeding.

Subtyping analysis reveals new variants and accelerated evolution of Clostridioides difficile toxin B.

Clostridioides difficile toxins (TcdA and TcdB) are major exotoxins responsible for C. difficile infection (CDI) associated diseases. The previously reported TcdB variants showed distinct biological features, immunoactivities, and potential pathogenicity in disease progression. Here, we performed global comparisons of amino acid sequences of both TcdA and TcdB from 3,269 C. difficile genomes and clustered them according to the evolutionary relatedness. We found that TcdB was much diverse and could be divided into eight subtypes, of which four were first described. Further analysis indicates that the tcdB gene undergoes accelerated evolution to maximize diversity. By tracing TcdB subtypes back to their original isolates, we found that the distribution of TcdB subtypes was not completely aligned with the phylogeny of C. difficile. These findings suggest that the tcdB genes not only frequently mutate, but also continuously transfer and exchange among C. difficile strains.

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.

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.

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.