Chuyu Ye

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

Endogenous target mimecry of an miRNA affects nicotine biosynthesis. 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.

Echinochloa Chloroplast Genomes: Insights into the Evolution and Taxonomic Identification of Two Weedy Species

The genus Echinochloa (Poaceae) includes numerous problematic weeds that cause the reduction of crop yield worldwide. To date, DNA sequence information is still limited in the genus Echinochloa. In this study, we completed the entire chloroplast genomes of two Echinochloa species (Echinochloa oryzicola and Echinochloa crus-galli) based on high-throughput sequencing data from their fresh green leaves. The two Echinochloa chloroplast genomes are 139,891 and 139,800 base pairs in length, respectively, and contain 131 protein-coding genes, 79 indels and 466 substitutions helpful for discrimination of the two species. The divergence between the genus Echinochloa and Panicum occurred about 21.6 million years ago, whereas the divergence between E. oryzicola and E. crus-galli chloroplast genes occurred about 3.3 million years ago. The two reported Echinochloa chloroplast genome sequences contribute to better understanding of the diversification of this genus.

PlantcircBase: A Database for Plant Circular RNAs

Circular RNAs (circRNAs) are a kind of non-coding RNAs; they were first found over 35 years ago but have only become a research hotspot in recent years. CircRNAs are covalently closed loops derived from precursor mRNAs by non-sequential back-splicing. Based on their genomic location, circRNAs can be classified into exonic, intronic, intergenic, and exon-intronic circRNAs (Chen, 2016). By taking advantages of recent advances in bioinformatics algorithms, a large number of circRNAs have been detected among diverse organisms, including plants such as Arabidopsis thaliana and Oryza sativa (Memczak et al., 2013; Ye et al., 2015; Chen, 2016).

Genome re-sequencing suggested a weedy rice origin from domesticated indica-japonica hybridization: a case study from southern China.

Main conclusionWhole-genome re-sequencing of weedy rice from southern China reveals that weedy rice can originate from hybridization of domesticatedindicaandjaponicarice.AbstractWeedy rice (Oryza sativa f. spontanea Rosh.), which harbors phenotypes of both wild and domesticated rice, has become one of the most notorious weeds in rice fields worldwide. While its formation is poorly understood, massive amounts of rice genomic data may provide new insights into this issue. In this study, we determined genomes of three weedy rice samples from the lower Yangtze region, China, and investigated their phylogenetics, population structure and chromosomal admixture patterns. The phylogenetic tree and principle component analysis based on 46,005 SNPs with 126 other Oryza accessions suggested that the three weedy rice accessions were intermediate between japonica and indica rice. An ancestry inference study further demonstrated that weedy rice had two dominant genomic components (temperate japonica and indica). This strongly suggests that weedy rice originated from indica-japonica hybridization. Furthermore, 22,443 novel fixed single nucleotide polymorphisms were detected in the weedy genomes and could have been generated after indica-japonica hybridization for environmental adaptation.

Genome-wide identification of non-coding RNAs interacted with microRNAs in soybean

A wide range of RNA species interacting with microRNAs (miRNAs) form a complex gene regulation network and play vital roles in diverse biological processes. In this study, we performed a genome-wide identification of endogenous target mimics (eTMs) for miRNAs and phased-siRNA-producing loci (PHAS) in soybean with a focus on those involved in lipid metabolism. The results showed that a large number of eTMs and PHAS genes could be found in soybean. Additionally, we found that lipid metabolism related genes were potentially regulated by 28 miRNAs, and nine of them were potentially further regulated by a number of eTMs with expression evidence. Thirty-three miRNAs were found to trigger production of phasiRNAs from 49 PHAS genes, which were able to target lipid metabolism related genes. Degradome data supported miRNA- and/or phasiRNA-mediated cleavage of genes involved in lipid metabolism. Most eTMs for miRNAs involved in lipid metabolism and phasiRNAs targeting lipid metabolism related genes showed a tissue-specific expression pattern. Our bioinformatical evidences suggested that lipid metabolism in soybean is potentially regulated by a complex non-coding network, including miRNAs, eTMs, and phasiRNAs, and the results extended our knowledge on functions of non-coding RNAs.

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.

Genome Re-Sequencing of Semi-Wild Soybean Reveals a Complex Soja Population Structure and Deep Introgression

Semi-wild soybean is a unique type of soybean that retains both wild and domesticated characteristics, which provides an important intermediate type for understanding the evolution of the subgenus Soja population in the Glycine genus. In this study, a semi-wild soybean line (Maliaodou) and a wild line (Lanxi 1) collected from the lower Yangtze regions were deeply sequenced while nine other semi-wild lines were sequenced to a 3-fold genome coverage. Sequence analysis revealed that (1) no independent phylogenetic branch covering all 10 semi-wild lines was observed in the Soja phylogenetic tree; (2) besides two distinct subpopulations of wild and cultivated soybean in the Soja population structure, all semi-wild lines were mixed with some wild lines into a subpopulation rather than an independent one or an intermediate transition type of soybean domestication; (3) high heterozygous rates (0.19–0.49) were observed in several semi-wild lines; and (4) over 100 putative selective regions were identified by selective sweep analysis, including those related to the development of seed size. Our results suggested a hybridization origin for the semi-wild soybean, which makes a complex Soja population structure.

PcircRNA_finder: a software for circRNA prediction in plants

Motivation: Recent studies reveal an important role of non-coding circular RNA (circRNA) in the control of cellular processes. Because of differences in the organization of plant and mammal genomes, the sensitivity and accuracy of circRNA prediction programs using algorithms developed for animals and humans perform poorly for plants. Results: A circRNA prediction software for plants (termed PcircRNA_finder) was developed that is more sensitive in detecting circRNAs than other frequently used programs (such as find_circ and CIRCexplorer), Based on analysis of simulated and real rRNA-/RNAase R RNA-Seq data from Arabidopsis thaliana and rice PcircRNA_finder provides a more comprehensive sensitive, precise prediction method for plants circRNAs. Availability and Implementation: http://ibi.zju.edu.cn/bioinplant/tools/manual.htm. Contact: fanlj@zju.edu.cn Supplementary information: Supplementary data are available at Bioinformatics online.

Full-length sequence assembly reveals circular RNAs with diverse non-GT/AG splicing signals in rice

ABSTRACT Circular RNAs (circRNAs) have been identified in diverse eukaryotic species and are characterized by RNA backsplicing events. Current available methods for circRNA identification are able to determine the start and end locations of circRNAs in the genome but not their full-length sequences. In this study, we developed a method to assemble the full-length sequences of circRNAs using the backsplicing RNA-Seq reads and their corresponding paired-end reads. By applying the method to an rRNA-depleted/RNase R-treated RNA-Seq dataset, we for the first time identified full-length sequences of nearly 3,000 circRNAs in rice. We further showed that alternative circularization of circRNA is a common feature in rice and, surprisingly, found that the junction sites of a large number of rice circRNAs are flanked by diverse non-GT/AG splicing signals while most human exonic circRNAs are flanked by canonical GT/AG splicing signals. Our study provides a method for genome-wide identification of full-length circRNAs and expands our understanding of splicing signals of circRNAs.

A diverse set of miRNAs responsive to begomovirus-associated betasatellite in Nicotiana benthamiana

BackgroundRoles of microRNAs (miRNAs) and short interfering RNAs (siRNAs) in biotic stress responses, e.g., viral infection, have been demonstrated in plants by many studies. Tomato yellow leaf curl China virus (TYLCCNV) is a monopartite begomovirus that can systemically infect Solanaceae plants, and induces leaf curling, yellowing and enation symptoms when co-inoculated with a betasatellite (TYLCCNB). The released genome sequence of Nicotiana benthamiana provides an opportunity to identify miRNAs and siRNAs responsive to begomovirus-associated betasatellite in N. benthamiana.ResultsmiRNAs were identified in three small RNA libraries generated using RNA isolated from N. benthamiana plants systemically infected with TYLCCNV (Y10A) alone, co-infected with Y10A and its betasatellite TYLCCNB (Y10β) or a TYLCCNB mutant (Y10mβ) that contains a mutated βC1, the sole betasatellite-encoded protein. A total of 196 conserved miRNAs from 38 families and 197 novel miRNAs from 160 families were identified. Northern blot analysis confirmed that expression of species-specific miRNAs was much lower than that of conserved miRNAs. Several conserved and novel miRNAs were found to be responsive to co-infection of Y10A and Y10β but not to co-infection of Y10A and Y10mβ, suggesting that these miRNAs might play a role unique to interaction between Y10β and N. benthamiana. Additionally, we identified miRNAs that can trigger the production of phased secondary siRNAs (phasiRNAs).ConclusionsIdentification of miRNAs with differential expression profiles in N. benthamiana co-infected with Y10A and Y10β and co-infected with Y10A and Y10mβ indicates that these miRNAs are betasatellite-responsive. Our result also suggested a potential role of miRNA-mediated production of phasiRNAs in interaction between begomovirus and N. benthamiana.