Deep sequencing of Lotus corniculatus L. reveals key enzymes and potential transcription factors related to the flavonoid biosynthesis pathway

Illuminating the Plant Rhabdovirus Landscape through Metatranscriptomics Data

Rhabdoviruses infect a large number of plant species and cause significant crop diseases. They have a negative-sense, single-stranded unsegmented or bisegmented RNA genome. The number of plant-associated rhabdovirid sequences has grown in the last few years in concert with the extensive use of high-throughput sequencing platforms. Here, we report the discovery of 27 novel rhabdovirus genomes associated with 25 different host plant species and one insect, which were hidden in public databases. These viral sequences were identified through homology searches in more than 3000 plant and insect transcriptomes from the National Center for Biotechnology Information (NCBI) Sequence Read Archive (SRA) using known plant rhabdovirus sequences as the query. The identification, assembly and curation of raw SRA reads resulted in sixteen viral genome sequences with full-length coding regions and ten partial genomes. Highlights of the obtained sequences include viruses with unique and novel genome organizations among known plant rhabdoviruses. Phylogenetic analysis showed that thirteen of the novel viruses were related to cytorhabdoviruses, one to alphanucleorhabdoviruses, five to betanucleorhabdoviruses, one to dichorhaviruses and seven to varicosaviruses. These findings resulted in the most complete phylogeny of plant rhabdoviruses to date and shed new light on the phylogenetic relationships and evolutionary landscape of this group of plant viruses. Furthermore, this study provided additional evidence for the complexity and diversity of plant rhabdovirus genomes and demonstrated that analyzing SRA public data provides an invaluable tool to accelerate virus discovery, gain evolutionary insights and refine virus taxonomy.

2020 taxonomic update for phylum Negarnaviricota (Riboviria: Orthornavirae), including the large orders Bunyavirales and Mononegavirales

In March 2020, following the annual International Committee on Taxonomy of Viruses (ICTV) ratification vote on newly proposed taxa, the phylum Negarnaviricota was amended and emended. At the genus rank, 20 new genera were added, two were deleted, one was moved, and three were renamed. At the species rank, 160 species were added, four were deleted, ten were moved and renamed, and 30 species were renamed. This article presents the updated taxonomy of Negarnaviricota as now accepted by the ICTV.

The Research Progress of Chalcone Isomerase (CHI) in Plants

Chalcone isomerase (CHI) is the second rate-limiting and the first reported enzyme involved in the biosynthetic pathway of flavonoids. It catalyzes the intramolecular cyclization reaction, converting the bicyclic chalcone into tricyclic (2S)-flavanone. In this paper, we obtained and analyzed 916 DNA sequences, 1310 mRNA sequences, and 2403 amino acid sequences of CHI registered in NCBI by Jan 2018. The full length of CHI DNA sequences ranges from 218 to 3758 bp, CHI mRNA sequences ranges from 265 to 1436 bp, and CHI amino acid sequences ranges from 35 to 465 amino acid residues. Forty representative species were selected from each family to construct the maximum likelihood tree and analyze the evolutionary relationship. According to the medicinal and agricultural use, 13 specific species were selected, and their physicochemical properties were analyzed. The molecular weight of CHI ranges from 23 to 26 kD, and the isoelectric point of CHI ranges from 4.93 to 5.85. All the half-life periods of CHI are 30 h in mammalian reticulocytes in vitro, 20 h in yeast, and 10 h in E. coli in vivo, theoretically. The consistency of the 13 CHI amino acid sequences is 63.55%. According to the similarity between each sequence, we selected four CHI sequences of Paeonia suffruticosa, Paeonia lactiflora, Taxus wallichiana, and Tradescantia hirsutiflora for secondary structure, three-dimensional protein models, conserved domains, transmembrane structure, and signal peptide prediction analysis. It was found that CHI sequences of Paeonia suffruticosa and Paeonia lactiflora owned a higher similarity; they both share the template 4doi.1.A. The four CHI all have no signal peptides, and they exert their activities in cytoplasm. Then, PubMed, Web of Science, Science Direct, and Research Gate were used as information sources through the search terms 'chalcone isomerase', 'biosynthesis', 'expression', and their combinations to get the latest and comprehensive information of CHI, mainly from the year 2010 to 2018. More than 300 papers were searched and 116 papers were reviewed in the present work. We summarized the classification of CHI, catalytic reaction mechanism of CHI, and progress of genetic engineering regarding CHI clone, expression, and exogenous stimulator regulation. This paper will lay a foundation for further studies of CHI and other functional genes involved in flavonoids biosynthetic pathway.

The R2R3-MYB TT2b and the bHLH TT8 genes are the major regulators of proanthocyanidin biosynthesis in the leaves of Lotus species

By exploiting interspecific hybrids and their progeny, we identified key regulatory and transporter genes intimately related to proanthocyanidin biosynthesis in leaves of Lotus spp. Proanthocyanidins (PAs), known as condensed tannins, are polymeric flavonoids enriching forage legumes of key nutritional value to prevent bloating in ruminant animals. Unfortunately, major forage legumes such as alfalfa and clovers lack PAs in edible tissues. Therefore, engineering the PA trait in herbage of forage legumes is paramount to improve both ecological and economical sustainability of cattle production system. Progresses on the understanding of genetic determinants controlling PA biosynthesis and accumulation have been mainly made studying mutants of Arabidopsis, Medicago truncatula and Lotus japonicus, model species unable to synthesize PAs in the leaves. Here, we exploited interspecific hybrids between Lotus corniculatus, with high levels of PAs in the leaves, and Lotus tenuis, with no PAs in these organs, and relative F₂ progeny, to identify among candidate PA regulators and transporters the genes mainly affecting this trait. We found that the levels of leaf PAs significantly correlate with the expression of MATE1, the putative transporter of glycosylated PA monomers, and, among the candidate regulatory genes, with the expression of the MYB genes TT2a, TT2b and MYB14 and the bHLH gene TT8. The expression levels of TT2b and TT8 also correlated with those of all key structural genes of the PA pathways investigated, MATE1 included. Our study unveils a different involvement of the three Lotus TT2 paralogs to the PA trait and highlights differences in the regulation of this trait in our Lotus genotypes with respect to model species. This information opens new avenues for breeding bloat safe forage legumes.

Comparative transcriptome investigation of global gene expression changes caused by miR156 overexpression in Medicago sativa

Background

Medicago sativa (alfalfa) is a low-input forage and potential bioenergy crop, and improving its yield and quality has always been a focus of the alfalfa breeding industry. Transgenic alfalfa plants overexpressing a precursor of alfalfa microRNA156 (MsmiR156) were recently generated by our group. These plants (miR156OE) showed enhanced biomass yield, reduced internodal length, increased shoot branching and trichome density, and a delay in flowering time. Transcripts of three SQUAMOSA-PROMOTER BINDING PROTEIN-LIKE (SPL) genes (MsSPL6, MsSPL12, and MsSPL13) were found to be targeted for cleavage by MsmiR156 in alfalfa.

Results

To further illustrate the molecular mechanisms underlying the effects of miR156 in alfalfa, two miR156OE genotypes (A11a and A17) were subjected to Next Generation RNA Sequencing with Illumina HiSeq. More than 1.11 billion clean reads were obtained from our available sequenced samples. A total of 160,472 transcripts were generated using Trinity de novo assembly and 4,985 significantly differentially expressed genes were detected in miR156OE plants A11a and A17 using the Medicago truncatula genome as reference. A total of 17 genes (including upregulated, downregulated, and unchanged) were selected for quantitative real-time PCR (qRT-PCR) validation, which showed that gene expression levels were largely consistent between qRT-PCR and RNA-Seq data. In addition to the established SPL genes MsSPL6, MsSPL12 and MsSPL13, four new SPLs; MsSPL2, MsSPL3, MsSPL4 and MsSPL9 were also down-regulated significantly in both miR156OE plants. These seven SPL genes belong to genes phylogeny clades VI, IV, VIII, V and VII, which have been reported to be targeted by miR156 in Arabidopsis thaliana. The gene ontology terms characterized electron transporter, starch synthase activity, sucrose transport, sucrose-phosphate synthase activity, chitin binding, sexual reproduction, flavonoid biosynthesis and lignin catabolism correlate well to the phenotypes of miR156OE alfalfa plants.

Conclusions

This is the first report of changes in global gene expression in response to miR156 overexpression in alfalfa. The discovered miR156-targeted SPL genes belonging to different clades indicate miR156 plays fundamental and multifunctional roles in regulating alfalfa plant development.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-016-3014-6) contains supplementary material, which is available to authorized users.

De novo transcriptome assembly of Ipomoea nil using Illumina sequencing for gene discovery and SSR marker identification

Ipomoea nil is widely used as an ornamental plant due to its abundance of flower color, but the limited transcriptome and genomic data hinder research on it. Using illumina platform, transcriptome profiling of I. nil was performed through high-throughput sequencing, which was proven to be a rapid and cost-effective means to characterize gene content. Our goal is to use the resulting information to facilitate the relevant research on flowering and flower color formation in I. nil. In total, 268 million unique illumina RNA-Seq reads were produced and used in the transcriptome assembly. These reads were assembled into 220,117 contigs, of which 137,307 contigs were annotated using the GO and KEGG database. Based on the result of functional annotations, a total of 89,781 contigs were assigned 455,335 GO term annotations. Meanwhile, 17,418 contigs were identified with pathway annotation and they were functionally assigned to 144 KEGG pathways. Our transcriptome revealed at least 55 contigs as probably flowering-related genes in I. nil, and we also identified 25 contigs that encode key enzymes in the phenylpropanoid biosynthesis pathway. Based on the analysis relating to gene expression profiles, in the phenylpropanoid biosynthesis pathway of I. nil, the repression of lignin biosynthesis might lead to the redirection of the metabolic flux into anthocyanin biosynthesis. This may be the most likely reason that I. nil has high anthocyanins content, especially in its flowers. Additionally, 15,537 simple sequence repeats (SSRs) were detected using the MISA software, and these SSRs will undoubtedly benefit future breeding work. Moreover, the information uncovered in this study will also serve as a valuable resource for understanding the flowering and flower color formation mechanisms in I. nil.

An insight into the genes involved in secoiridoid biosynthesis in Gentiana macrophylla by RNA-seq

The dried root of Gentiana macrophylla is a well-known traditional Chinese herbal medicine for treating jaundice, hepatitis, and stomachic and choleretic ailments. However, natural sources are now in short supply. A lack of information about its genetic background has been a great hindrance to producing its active constituents via genetic engineering. We performed RNA-seq to obtain 42,918 unigenes (average length = 667 bp) in its transcriptome. Of these, 32,141 (74.89 %) were annotated and 2,339 unigenes were assigned to secondary-metabolite pathways. In all, 114 putative unigenes involved in secoiridoid biosynthesis were identified in our transcriptome library. A Blast X search against the Arabidopsis gene regulatory information server showed that 4,413 unigenes are homologous to transcription factor genes from Arabidopsis. Organ-specific genes and candidate gene expression profiles were also investigated with digital gene expression technology. Quantitative PCR was used to verify the expression patterns of several novel transcripts involved in secoiridoid biosynthesis. Our results not only enrich the gene resource but will also benefit research into the molecular genetics and functional genomics of this species.