RNA-Seq Uncovers SNPs and Alternative Splicing Events in Asian Lotus (Nelumbo nucifera)

Resequencing of 296 cultivated and wild lotus accessions unravels its evolution and breeding history

Lotus (family: Nelumbonaceae) are perennial aquatic plants that represent one of the most ancient basal dicots. In the present study, we resequenced 296 lotus accessions from various geographical locations and germplasms to explore their genomic diversity and population structure. This germplasm set consisted of four accessions of American wild lotus and 292 accessions of Asian lotus, which were divided into four subgroups: wild, rhizome, flower and seed. Total single nucleotide polymorphisms (SNPs) suggested that the wild lotus had the highest variant number (7 191 010). Population structure and genome diversity analysis indicated that the American wild lotus demonstrated a distant genetic relationship with the Asian lotus. Furthermore, the seed and rhizome lotus groups had not originated from a single source but rather had a more complex multisource origin. Besides that, the seed lotus showed higher genetic diversity, which might have been due to the gene flow from the flower lotus to seed lotus by artificial crossing, and the rhizome lotus showed a much lower genetic diversity than the other groups. The present study provides SNP markers for lotus genomic diversity analysis, which will be useful for guiding lotus breeding.

Transcription Analysis of Recombinant Trichoderma reesei HJ-48 to Compare the Molecular Basis for Fermentation of Glucose and Xylose

Profiling the transcriptome changes involved in xylose metabolism by the fungus Trichoderma reesei allows for the identification of potential targets for ethanol production processing. In the present study, the transcriptome of T. reesei HJ-48 grown on xylose versus glucose was analyzed using nextgeneration sequencing technology. During xylose fermentation, numerous genes related to central metabolic pathways, including xylose reductase (XR) and xylitol dehydrogenase (XDH), were expressed at higher levels in T. reesei HJ-48. Notably, growth on xylose did not fully repress the genes encoding enzymes of the tricarboxylic acid and respiratory pathways. In addition, increased expression of several sugar transporters was observed during xylose fermentation. This study provides a valuable dataset for further investigation of xylose fermentation and provides a deeper insight into the various genes involved in this process.

Integrated Omics Analyses Identify Key Pathways Involved in Petiole Rigidity Formation in Sacred Lotus

Sacred lotus (Nelumbo nucifera Gaertn.) is a relic aquatic plant with two types of leaves, which have distinct rigidity of petioles. Here we assess the difference from anatomic structure to the expression of genes and proteins in two petioles types, and identify key pathways involved in petiole rigidity formation in sacred lotus. Anatomically, great variation between the petioles of floating and vertical leaves were observed. The number of collenchyma cells and thickness of xylem vessel cell wall was higher in the initial vertical leaves’ petiole (IVP) compared to the initial floating leaves’ petiole (IFP). Among quantified transcripts and proteins, 1021 and 401 transcripts presented 2-fold expression increment (named DEGs, genes differentially expressed between IFP and IVP) in IFP and IVP, 421 and 483 proteins exhibited 1.5-fold expression increment (named DEPs, proteins differentially expressed between IFP and IVP) in IFP and IVP, respectively. Gene function and pathway enrichment analysis displayed that DEGs and DEPs were significantly enriched in cell wall biosynthesis and lignin biosynthesis. In consistent with genes and proteins expressions in lignin biosynthesis, the contents of lignin monomers precursors were significantly different in IFP and IVP. These results enable us to understand lotus petioles rigidity formation better and provide valuable candidate genes information on further investigation.

The complexity of alternative splicing and landscape of tissue-specific expression in lotus (Nelumbo nucifera) unveiled by Illumina- and single-molecule real-time-based RNA-sequencing

Alternative splicing (AS) plays a critical role in regulating different physiological and developmental processes in eukaryotes, by dramatically increasing the diversity of the transcriptome and the proteome. However, the saturation and complexity of AS remain unclear in lotus due to its limitation of rare obtainment of full-length multiple-splice isoforms. In this study, we apply a hybrid assembly strategy by combining single-molecule real-time sequencing and Illumina RNA-seq to get a comprehensive insight into the lotus transcriptomic landscape. We identified 211,802 high-quality full-length non-chimeric reads, with 192,690 non-redundant isoforms, and updated the lotus reference gene model. Moreover, our analysis identified a total of 104,288 AS events from 16,543 genes, with alternative 3ʹ splice-site being the predominant model, following by intron retention. By exploring tissue datasets, 370 tissue-specific AS events were identified among 12 tissues. Both the tissue-specific genes and isoforms might play important roles in tissue or organ development, and are suitable for ‘ABCE’ model partly in floral tissues. A large number of AS events and isoform variants identified in our study enhance the understanding of transcriptional diversity in lotus, and provide valuable resource for further functional genomic studies.

The Latest Studies on Lotus (Nelumbo nucifera)-an Emerging Horticultural Model Plant

Lotus (Nelumbo nucifera) is a perennial aquatic basal eudicot belonging to a small family Nelumbonaceace, which contains only one genus with two species. It is an important horticultural plant, with its uses ranging from ornamental, nutritional to medicinal values, and has been widely used, especially in Southeast Asia. Recently, the lotus obtained a lot of attention from the scientific community. An increasing number of research papers focusing on it have been published, which have shed light on the mysteries of this species. Here, we comprehensively reviewed the latest advancement of studies on the lotus, including phylogeny, genomics and the molecular mechanisms underlying its unique properties, its economic important traits, and so on. Meanwhile, current limitations in the research of the lotus were addressed, and the potential prospective were proposed as well. We believe that the lotus will be an important model plant in horticulture with the generation of germplasm suitable for laboratory operation and the establishment of a regeneration and transformation system.

Integrated Analysis of LncRNA-mRNA Coexpression in the Extracellular Matrix of Developing Deciduous Teeth in Miniature Pigs

Miniature pigs, a valuable alternative model for understanding human tooth development, have deciduous teeth from all four tooth families that are replaced once by permanent molars. The extracellular matrix (ECM) supports cells and maintains the integrity of tooth germs during tooth development. However, details on the role of the ECM in tooth development are poorly understood. Here, we performed long noncoding RNA (lncRNA) and messenger RNA (mRNA) expression profiles in the ECM components of deciduous tooth germs by RNA sequencing in miniature pigs. From the early cap to the late bell stages, we identified 4,562 and 3,238 differentially expressed genes (DEGs) from E40 to E50 and E50 to E60, respectively. In addition, a total of 1,464 differentially expressed lncRNAs from E40 to E50 and 969 differentially expressed lncRNAs from E50 to E60 were obtained. Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis showed that DEGs were enriched significantly for multiple signaling pathways, especially for the ECM pathway. We then outlined the detailed dynamic gene expression profiling of ECM components during deciduous molar development. Comparison of the cap and bell stages revealed that the structure and functions of the ECM dynamically changed. The ECM-related genes, including THBS1, COL4A5, COL4A6, COL1A1, CHAD, TNR, GP1BA, and ITGA3, were significantly changed, and some were shown to enrich during the bell stage development. Finally, we outlined the coexpression of lncRNAs and ECM properties during tooth development. We showed that the interplay of key lncRNAs could change ECM processes and influence the ECM establishment of tooth patterns to accomplish full tooth formation. These results might provide information to elucidate the regulation network of the lncRNA and ECM in tooth development.

Alternative splicing and translation play important roles in hypoxic germination in rice

Post-transcriptional mechanisms (PTMs), including alternative splicing (AS) and alternative translation initiation (ATI), may explain the diversity of proteins involved in plant development and stress responses. Transcriptional regulation is important during the hypoxic germination of rice seeds, but the potential roles of PTMs in this process have not been characterized. We used a combination of proteomics and RNA sequencing to discover how AS and ATI contribute to plant responses to hypoxia. In total, 10 253 intron-containing genes were identified. Of these, ~1741 differentially expressed AS (DAS) events from 811 genes were identified in hypoxia-treated seeds compared with controls. Over 95% of these were not present in the list of differentially expressed genes. In particular, regulatory pathways such as the spliceosome, ribosome, endoplasmic reticulum protein processing and export, proteasome, phagosome, oxidative phosphorylation, and mRNA surveillance showed substantial AS changes under hypoxia, suggesting that AS responses are largely independent of transcriptional regulation. Considerable AS changes were identified, including the preferential usage of some non-conventional splice sites and enrichment of splicing factors in the DAS data sets. Taken together, these results not only demonstrate that AS and ATI function during hypoxic germination but they have also allowed the identification of numerous novel proteins/peptides produced via ATI.

Transcriptomic analysis identifies the key genes involved in stamen petaloid in lotus (Nelumbo nucifera)

Background

Flower morphology, a phenomenon regulated by a complex network, is one of the vital ornamental features in Nelumbo nucifera. Stamen petaloid is very prevalent in lotus flowers. However, the mechanism underlying this phenomenon is still obscure.

Results

Here, the comparative transcriptomic analysis was performed among petal, stamen petaloid and stamen through RNA-seq. Using pairwise comparison analysis, a large number of genes involved in hormonal signal transduction pathways and transcription factors, especially the MADS-box genes, were identified as candidate genes for stamen petaloid in lotus.

Conclusions

Taken together, these results provide an insight into the molecular networks underlying lotus floral organ development and stamen petaloid.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-4950-0) contains supplementary material, which is available to authorized users.

Gene-based SNP identification and validation in soybean using next-generation transcriptome sequencing

Gene-based molecular markers are increasingly used in crop breeding programs for marker-assisted selection. However, identification of genetic variants associated with important agronomic traits has remained a difficult task in soybean. RNA-Seq provides an efficient way, other than assessing global expression variations of coding genes, to discover gene-based SNPs at the whole genome level. In this study, RNA isolated from four soybean accessions each with three replications was subjected to high-throughput sequencing and a range of 44.2-65.9 million paired-end reads were generated for each library. A total of 75,209 SNPs were identified among different genotypes after combination of replications, 89.1% of which were located in expressed regions and 27.0% resulted in amino acid changes. GO enrichment analysis revealed that most significant enriched genes with nonsynonymous SNPs were involved in ribonucleotide binding or catalytic activity. Of 22 SNPs subjected to PCR amplification and Sanger sequencing, all of them were validated. To test the utility of identified SNPs, these validated SNPs were also assessed by genotyping a relative large population with 393 wild and cultivated soybean accessions. These SNPs identified by RNA-Seq provide a useful resource for genetic and genomic studies of soybean. Moreover, the collection of nonsynonymous SNPs annotated with their predicted functional effects also provides a valuable asset for further discovery of genes, identification of gene variants, and development of functional markers.

Comparative genomics and transcriptome analysis of Lactobacillus rhamnosus ATCC 11443 and the mutant strain SCT-10-10-60 with enhanced L-lactic acid production capacity

Mechanisms for high L-lactic acid production remain unclear in many bacteria. Lactobacillus rhamnosus SCT-10-10-60 was previously obtained from L. rhamnosus ATCC 11443 via mutagenesis and showed improved L-lactic acid production. In this study, the genomes of strains SCT-10-10-60 and ATCC 11443 were sequenced. Both genomes are a circular chromosome, 2.99 Mb in length with a GC content of approximately 46.8%. Eight split genes were identified in strain SCT-10-10-60, including two LytR family transcriptional regulators, two Rex redox-sensing transcriptional repressors, and four ABC transporters. In total, 60 significantly up-regulated genes (log₂fold-change ≥ 2) and 39 significantly down-regulated genes (log₂fold-change ≤ - 2) were identified by a transcriptome comparison between strains SCT-10-10-60 and ATCC 11443. KEGG pathway enrichment analysis revealed that "pyruvate metabolism" was significantly different (P < 0.05) between the two strains. The split genes and the differentially expressed genes involved in the "pyruvate metabolism" pathway are probably responsible for the increased L-lactic acid production by SCT-10-10-60. The genome and transcriptome sequencing information and comparison of SCT-10-10-60 with ATCC 11443 provide insights into the anabolism of L-lactic acid and a reference for improving L-lactic acid production using genetic engineering.

Proteogenomic analysis reveals alternative splicing and translation as part of the abscisic acid response in Arabidopsis seedlings

In eukaryotes, mechanisms such as alternative splicing (AS) and alternative translation initiation (ATI) contribute to organismal protein diversity. Specifically, splicing factors play crucial roles in responses to environment and development cues; however, the underlying mechanisms are not well investigated in plants. Here, we report the parallel employment of short-read RNA sequencing, single molecule long-read sequencing and proteomic identification to unravel AS isoforms and previously unannotated proteins in response to abscisic acid (ABA) treatment. Combining the data from the two sequencing methods, approximately 83.4% of intron-containing genes were alternatively spliced. Two AS types, which are referred to as alternative first exon (AFE) and alternative last exon (ALE), were more abundant than intron retention (IR); however, by contrast to AS events detected under normal conditions, differentially expressed AS isoforms were more likely to be translated. ABA extensively affects the AS pattern, indicated by the increasing number of non-conventional splicing sites. This work also identified thousands of unannotated peptides and proteins by ATI based on mass spectrometry and a virtual peptide library deduced from both strands of coding regions within the Arabidopsis genome. The results enhance our understanding of AS and alternative translation mechanisms under normal conditions, and in response to ABA treatment.

Molecular Mechanisms of Developmental Toxicity Induced by Graphene Oxide at Predicted Environmental Concentrations

Developmental toxicity is a critical issue in nanotoxicity. However, very little is known about the effects of graphene oxide (GO, a widely used carbon material) at predicted environmental concentrations on biological development or the specific molecular mechanisms. The present study established that the development of zebrafish embryos exposed to trace concentrations (1-100 μg/L) of GO was impaired because of DNA modification, protein carbonylation and excessive generation of reactive oxygen species (ROS), especially the superoxide radical. Noticeably, there was a nonmonotonic response of zebrafish developmental toxicity to GO at μg/L to mg/L levels. Transcriptomics analysis revealed that disturbing collagen- and matrix metalloproteinase (MMP)-related genes affected the skeletal and cardiac development of zebrafish. Moreover, metabolomics analysis showed that the inhibition of amino acid metabolism and the ratios of unsaturated fatty acids (UFAs) to saturated fatty acids (SFAs) contributed to the above developmental toxicity. The present work verifies the developmental toxicity of GO at trace concentrations and illustrates for the first time the specific molecular mechanisms thereof. Because of the potential developmental toxicity of GO at trace concentrations, government administrators and nanomaterial producers should consider its potential risks prior to the widespread environmental exposure to GO.

ACE: an efficient and sensitive tool to detect insecticide resistance-associated mutations in insect acetylcholinesterase from RNA-Seq data

Background

Insecticide resistance is a substantial problem in controlling agricultural and medical pests. Detecting target site mutations is crucial to manage insecticide resistance. Though PCR-based methods have been widely used in this field, they are time-consuming and inefficient, and typically have a high false positive rate. Acetylcholinesterases (Ace) is the neural target of the widely used organophosphate (OP) and carbamate insecticides. However, there is not any software available to detect insecticide resistance associated mutations in RNA-Seq data at present.

Results

A computational pipeline ACE was developed to detect resistance mutations of ace in insect RNA-Seq data. Known ace resistance mutations were collected and used as a reference. We constructed a Web server for ACE, and the standalone software in both Linux and Windows versions is available for download. ACE was used to analyse 971 RNA-Seq data from 136 studies in 7 insect pests. The mutation frequency of each RNA-Seq dataset was calculated. The results indicated that the resistance frequency was 30%–44% in an eastern Ugandan Anopheles population, thus suggesting this resistance-conferring mutation has reached high frequency in these mosquitoes in Uganda. Analyses of RNA-Seq data from the diamondback moth Plutella xylostella indicated that the G227A mutation was positively related with resistance levels to organophosphate or carbamate insecticides. The wasp Nasonia vitripennis had a low frequency of resistant reads (<5%), but the agricultural pests Chilo suppressalis and Bemisia tabaci had a high resistance frequency. All ace reads in the 30 B. tabaci RNA-Seq data were resistant reads, suggesting that insecticide resistance has spread to very high frequency in B. tabaci.

Conclusions

To the best of our knowledge, the ACE pipeline is the first tool to detect resistance mutations from RNA-Seq data, and it facilitates the full utilization of large-scale genetic data obtained by using next-generation sequencing.

Electronic supplementary material

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

Small RNA and degradome profiling reveals miRNA regulation in the seed germination of ancient eudicot Nelumbo nucifera

BACKGROUND

MicroRNAs (miRNAs) play important roles in plant growth and development. MiRNAs and their targets have been widely studied in model plants, but limited knowledge is available concerning this small RNA population and their targets in sacred lotus (Nelumbo nucifera Gaertn.).

RESULTS

In this study, a total of 145 known miRNAs belonging to 47 families and 78 novel miRNAs were identified during seed germination using high-throughput small RNA sequencing. Furthermore, some miRNA families which have not yet been reported in monocot or eudicot species were detected in N. nucifera, indicating that these miRNAs was divergence from monocots and core eudicots during evolution. Using degradome sequencing, 2580 targets were detected for all the miRNAs. GO (Gene Ontology) and KEGG pathway analyses showed that many target genes enriched in "regulation of transcription" and involved in "carbohydrate", "amino acid and energy metabolism". Nine miRNAs and three corresponding targets of them were further validated by quantitative RT-PCR.

CONCLUSIONS

The results present here suggested that many miRNAs were involved in the regulation of seed germination of sacred lotus, providing a foundation for future studies of sacred lotus seed longevity. Comparative analysis of miRNAs from different plants also provided insight into the evolutionary gains and losses of miRNAs in plants.

Metabolomic and Proteomic Profiles Reveal the Dynamics of Primary Metabolism during Seed Development of Lotus (Nelumbo nucifera)

Sacred lotus (Nelumbo nucifera) belongs to the Nelumbonaceae family. Its seeds are widely consumed in Asian countries as snacks or even medicine. Besides the market value, lotus seed also plays a crucial role in the lotus life cycle. Consequently, it is essential to gain a comprehensive understanding of the development of lotus seed. During its development, lotus seed undergoes cell division, expansion, reserve accumulation, desiccation, and maturation phases. We observed morphological and biochemical changes from 10 to 25 days after pollination (DAP) which corresponded to the reserve synthesis and accumulation phase. The volume of the seed expanded until 20 DAP with the color of the seed coat changing from yellow-green to dark green and gradually fading again. Starch and protein rapidly accumulated from 15 to 20 DAP. To further reveal metabolic adaptation, primary metabolites and proteins profiles were obtained using mass spectrometry based platforms. Metabolites and enzymes involved in sugar metabolism, glycolysis, TCA cycle and amino acid metabolism showed sequential dynamics enabling the clear separation of the different metabolic states during lotus seed development. The integration of the data revealed a highly significant metabolic switch at 15 DAP going through a transition of metabolically highly active tissue to the preparation of storage tissue. The results provide a reference data set for the evaluation of primary metabolism during lotus seed development.

Development and characterization of polymorphic microRNA-based microsatellite markers in Nelumbo nucifera (Nelumbonaceae)

PREMISE OF THE STUDY

Polymorphic microRNA (miRNA)-based microsatellite markers were developed to investigate the genetic diversity and population structure of Nelumbo nucifera (Nelumbonaceae).

METHODS AND RESULTS

A total of 485 miRNA-based microsatellites were found from the genomic DNA sequences of N. nucifera. After several rounds of screening, 21 primer pairs flanking di-, tri-, or pentanucleotide repeats were identified that revealed high levels of genetic diversity in four populations with two to five alleles per locus. The observed and expected heterozygosity per locus ranged from 0.000 to 1.000 and from 0.000 to 0.803, respectively.

CONCLUSIONS

The polymorphic microsatellite markers will be useful for studying the genetic diversity and population structure of N. nucifera.

Genome-Wide Identification of SSR and SNP Markers Based on Whole-Genome Re-Sequencing of a Thailand Wild Sacred Lotus (Nelumbo nucifera)

Genomic resources such as single nucleotide polymorphism (SNPs), insertions and deletions (InDels) and SSRs (simple sequence repeats) are essential for crop improvement and better utilization in genetic breeding. However, the resources for the sacred lotus (Nelumbo nucifera Gaertn.) are still limited. In the present study, to dissect large-scale genomic molecular marker resources for sacred lotus, we re-sequenced a Thailand sacred lotus cultivar ‘Chiang Mai wild lotus’ and compared with the reported lotus genome ‘Middle lake wild lotus’. A total of 3,180,059 SNPs, 328, 251 InDels and 14,191 SVs were found between the two genomes. The functional impact analyses of these SNPs indicated that they may be involved in metabolic processes, binding, catalytic activity, etc. Mining the genome sequences for SSRs showed that 191,657 SSRs were identified with a frequency of one SSR per 4.23 kb and 103,656 SSR primer pairs were designed. Furthermore, 14, 502 EST-SSRs were also indentified using the available RNA-seq data in the NCBI. A subset of 150 SSRs (genomic and EST-SSRs) was randomly selected for validation and genetic diversity analysis. The genotypes could be easily distinguished using these SSR markers and the ‘Chiang Mai wild lotus’ was obviously differentiated from the other Chinese accessions. This study provides considerable amounts of genomic resources and markers for the quantitative trait locus (QTL) identification and molecular selection of the species, which could have a potential role in various applications in sacred lotus breeding.

Cracking the Code of Human Diseases Using Next-Generation Sequencing: Applications, Challenges, and Perspectives

Next-generation sequencing (NGS) technologies have greatly impacted on every field of molecular research mainly because they reduce costs and increase throughput of DNA sequencing. These features, together with the technology's flexibility, have opened the way to a variety of applications including the study of the molecular basis of human diseases. Several analytical approaches have been developed to selectively enrich regions of interest from the whole genome in order to identify germinal and/or somatic sequence variants and to study DNA methylation. These approaches are now widely used in research, and they are already being used in routine molecular diagnostics. However, some issues are still controversial, namely, standardization of methods, data analysis and storage, and ethical aspects. Besides providing an overview of the NGS-based approaches most frequently used to study the molecular basis of human diseases at DNA level, we discuss the principal challenges and applications of NGS in the field of human genomics.