Transcriptomic Identification of Drought-Related Genes and SSR Markers in Sudan Grass Based on RNA-Seq

Identification and characterization of functionally relevant SSR markers in natural Dalbergia odorifera populations

BACKGROUND

Dalbergia odorifera is a rare and precious rosewood specie, which is valued for its amber tones, abstract figural patterns, and impermeability to water and insects. However, the information on genetic diversity and marker-assisted selection breeding of D. odorifera is still limited. Simple sequence repeat (SSR) markers are an ideal tool for genetic diversity analysis and marker-assisted molecular breeding for complex traits.

RESULTS

Here, we have developed SSR markers within candidate genes and used them to explore the genetic diversity among D. odorifera germplasm resources. A total of 635 SSR loci were identified. The proportions of mono-, di- and tri-nucleotide repeat motifs were 52.28%, 22.99% and 21.42%, respectively. From these, a total of 114 SSR primers were synthesized, of which 24 SSR markers displayed polymorphism (polymorphic information content (PIC) > 0.25). Subsequently, these polymorphic markers were used for the genetic diversity analysis of 106 D. odorifera individuals from 11 natural populations. According to the genetic diversity analysis of D. odorifera natural populations, the average observed heterozygosity (Ho) was 0.500, the average expected heterozygosity (He) was 0.524, and the average Shannon's information index (I) was 0.946. These indicated that the natural populations had moderate genetic diversity. AMOVA analysis showed that 5% of the total variation was within the individuals of a population, whereas 95% of the variation was among the individuals of the populations, indicating a high degree of genetic variation between populations. On the basis of their genetic structures, these populations could be divided into four groups.

CONCLUSIONS

Our study provides important experimental resources for genetic studies and assists in the program of molecular breeding of D. odorifera wood formation.

Molecular Mechanism Underlying the Sorghum sudanense (Piper) Stapf. Response to Osmotic Stress Determined via Single-Molecule Real-Time Sequencing and Next-Generation Sequencing

Drought, as a widespread environmental factor in nature, has become one of the most critical factors restricting the yield of forage grass. Sudangrass (Sorghum sudanense (Piper) Stapf.), as a tall and large grass, has a large biomass and is widely used as forage and biofuel. However, its growth and development are limited by drought stress. To obtain novel insight into the molecular mechanisms underlying the drought response and excavate drought tolerance genes in sudangrass, the first full-length transcriptome database of sudangrass under drought stress at different time points was constructed by combining single-molecule real-time sequencing (SMRT) and next-generation transcriptome sequencing (NGS). A total of 32.3 Gb of raw data was obtained, including 20,199 full-length transcripts with an average length of 1628 bp after assembly and correction. In total, 11,921 and 8559 up- and down-regulated differentially expressed genes were identified between the control group and plants subjected to drought stress. Additionally, 951 transcription factors belonging to 50 families and 358 alternative splicing events were found. A KEGG analysis of 158 core genes exhibiting continuous changes over time revealed that 'galactose metabolism' is a hub pathway and raffinose synthase 2 and β-fructofuranosidase are key genes in the response to drought stress. This study revealed the molecular mechanism underlying drought tolerance in sudangrass. Furthermore, the genes identified in this study provide valuable resources for further research into the response to drought stress.

A Comparison of Differential Gene Expression in Response to the Onset of Water Stress Between Three Hybrid Brachiaria Genotypes

Brachiaria (Trin.) Griseb. (syn. Urochloa P. Beauv.) is a C4 grass genus belonging to the Panicoideae. Native to Africa, these grasses are now widely grown as forages in tropical areas worldwide and are the subject of intensive breeding, particularly in South America. Tolerance to abiotic stresses such as aluminum and drought are major breeding objectives. In this study, we present the transcriptomic profiling of leaves and roots of three Brachiaria interspecific hybrid genotypes with the onset of water stress, Br12/3659-17 (gt-17), Br12/2360-9 (gt-9), and Br12/3868-18 (gt-18), previously characterized as having good, intermediate and poor tolerance to drought, respectively, in germplasm evaluation programs. RNA was extracted from leaf and root tissue of plants at estimated growing medium water contents (EWC) of 35, 15, and 5%. Differentially expressed genes (DEGs) were compared between different EWCs, 35/15, 15/5, and 35/5 using DESeq2. Overall, the proportions of DEGs enriched in all three genotypes varied in a genotype-dependent manner in relation to EWC comparison, with intermediate and sensitive gt-9 and gt-18 being more similar to each other than to drought tolerant gt-17. More specifically, GO terms relating to carbohydrate and cell wall metabolism in the leaves were enriched by up-regulated DEGs in gt-9 and gt-18, but by down-regulated DEGs in gt-17. Across all genotypes, analysis of DEG enzyme activities indicated an excess of down-regulated putative apoplastic peroxidases in the roots as water stress increased. This suggests that changes in root cell-wall architecture may be an important component of the response to water stress in Brachiaria.

Abundance and diversity of epiphytic microbiota on forage crops and their fermentation characteristic during the ensiling of sterile sudan grass

This study aimed to evaluate the effects of exogenous epiphytic microbiota inoculation on the fermentation quality and microbial community of sudan grass silage. Gamma irradiated sudan grass was ensiled with distilled water (STR), epiphytic microbiota of sudan grass (SUDm), forage sorghum (FSm), napier grass (NAPm) and whole crop corn (WCCm). The FSm inoculated silage have significantly lower lactic acid (LA) concentration and higher pH during early ensiling, while LA concentration gradually and significantly increased with the progression of ensiling and have lower pH in relation to other treatments for terminal silage. Inoculation of NAPm resulted in lower LA and higher acetic acid (AA) concentrations, higher pH, ammonia-N and dry matter losses for terminal silage, followed by SUDm silage. Inoculations of WCCm significantly increased LA production and pH decline during early ensiling and have higher LA and pH then NAPm and SUDm silages during final ensiling. The early fermentation of SUDm silage was dominated by genus of Pediococcus. The genera of Lactobacillus were predominant in WCCm and NAPm silages during 3 days of ensiling, while Weissella dominated initial microbial community of FS silage. The terminal silage of NAPm was dominated by Enterobacter and Rosenbergiella, while Enterobacter and Lactobacillus dominated terminal SUDm silage. The final silage of FSm was dominated by Lactobacillus, Weissella and Pediococcus, while Lactobacillus and Acetobacter dominated terminal WCCm silages. The results demonstrated that among the four forages the epiphytic microbiota from forage sorghum positively influenced the microbial community and fermentability of sudan grass silage.

Comparative analysis of drought-responsive transcriptomes of sugarcane genotypes with differential tolerance to drought

Water stress causes considerable yield losses in sugarcane. To investigate differentially expressed genes under water stress, two sugarcane genotypes were subjected to three water-deficit levels (mild, moderate, and severe) and subsequent recovery and leaf transcriptome was generated using Illumina NextSeq sequencing. Among the differentially expressed genes, the tolerant genotype Co 06022 generated 2970 unigenes (p ≤ 0.05, functionally known, non-redundant DEGs) at 2-day stress, and there was a progressive decrease in the expressed genes as the stress period increased with 2109 unigenes at 6-day stress and 2307 unigenes at 10-day stress. There was considerable reduction at recovery with 1334 unigenes expressed at 10 days after recovery. However, in the susceptible genotype Co 8021, the number of unigenes expressed at 2 days was lower (2025) than the tolerant genotype and a further reduction was seen at 6-day stress (1552). During recovery, more differentially expressed genes were observed in the susceptible cultivar indicating that the cultivar has to activate more functions/processes to recover from the damage caused by stress. Comparison of DEGs between all stages of stress and recovery in both genotypes revealed that, the commonly up- and down-regulated genes across different stages were approximately double in the tolerant genotype. The most enriched gene ontology classes were heme binding, peroxidase activity and metal ion binding in the biological process and response to oxidative stress, hydrogen peroxide catabolic process and response to stress in the molecular function category. The cellular component was enriched with DEGs involved in extracellular region followed by integral component of membrane. The KEGG pathway analysis revealed important metabolic activities and functionally important genes involved in mitigating water-deficit stress in both the varieties. In addition, several unannotated genes in important pathways were detected and together may provide novel insights into water-deficit tolerance mechanisms in sugarcane. The reliability of the observed expression patterns was confirmed by qRT-PCR. The results of this study will help to identify useful genes for improving drought tolerance in sugarcane.

Transcriptome sequencing of the apricot (Prunus armeniaca L.) and identification of differentially expressed genes involved in drought stress

Apricot (Prunus armeniaca L.) is an important fruit crop that is widely planted throughout the world. But drought affects both yield and quality of apricot. In order to study the effects of long-term drought on the molecular and physiological mechanisms of apricot, we used transcriptome sequencing and measured physiological indices. First, 322 million high-quality clean reads were obtained, and 74,892 unigenes were generated for the transcriptome. Among the assembled unigenes, 18,671 simple sequence repeats (SSRs) and 5581 differentially expressed genes (DEGs) were identified. The Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis of the DEGs revealed that starch and sucrose metabolism, plant-pathogen interaction and plant hormone signal transduction pathways are enriched. Additionally, we used quantitative real-time PCR (qRT-PCR) to confirm the RNA-seq results with 11 drought-related DEGs. Second, through the physiological analysis of apricot leaves under constant drought stress, and the results show the internal microstructure of apricot leaves changed to withstand drought stress. At the same time, plants exposed to long-term drought stress showed higher degree of membrane damage, which reduced photosynthesis in the damaged leaves. Our findings enrich the genome resources for apricot and refine our understanding of the molecular and physiological mechanisms of drought response in this fruit crop, providing insights into drought adaptation of the apricot.

De Novo Assembly and Discovery of Genes That Involved in Drought Tolerance in the Common Vetch

The common vetch (Vicia sativa) is often used as feed for livestock because of its high nutritional value. However, drought stress reduces forage production through plant damage. Here, we studied the transcriptional profiles of common vetch exposed to drought in order to understand the molecular mechanisms of drought tolerance in this species. The genome of the common vetch has not been sequenced, therefore we used Illumina sequencing to generate de novo transcriptomes. Nearly 500 million clean reads were used to generate 174,636 transcripts, including 122,299 unigenes. In addition, 5313 transcription factors were identified and these transcription factors were classified into 79 different gene families. We also identified 11,181 SSR loci from di- to hexa-nucleotides whose repeat number was greater than five. On the basis of differentially expressed genes, Gene Ontology analysis identified many drought-relevant categories, including “oxidation-reduction process”, “lipid metabolic process” and “oxidoreductase activity”. In addition to these, Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis identified pathways, such as “Plant hormone signal transduction”, “Glycolysis/Gluconeogenesis” and “Phenylpropanoid biosynthesis”, as differentially expressed in the plants exposed to drought. The expression results in this study will be useful for further extending our knowledge on the drought tolerance of common vetch.

High-throughput SSR marker development and its application in a centipedegrass (Eremochloa ophiuroides (Munro) Hack.) genetic diversity analysis

Centipedegrass (Eremochloa ophiuroides (Munro) Hack.) is a perennial, warm-season C4 grass species that shows great potential for use as a low-maintenance turfgrass species in tropical and subtropical regions. However, limited genetic and genomic information is available for this species, which has impeded systematic studies on the enhancement of its turf quality and resistance against biotic and abiotic stress. In this study, Illumina HiSeq high-throughput sequencing technology was performed to generate centipedegrass transcriptome sequences. A total of 352,513 assembled sequences were used to search for simple sequence repeat (SSR) loci, and 64,470 SSR loci were detected in 47,638 SSR containing sequences. The tri-nucleotides were the most frequent repeat motif, followed by di-nucleotides, tetra-nucleotides hexnucleotides, and pentanucleotides. A total of 48,061 primer pairs were successfully designed in the flanking sequences of the SSRs, and 100 sets of primers were randomly selected for the initial validation in four centipedegrass accessions. In total, 56 (56.0%) of the 100 primer pairs tested successfully amplified alleles from all four centipedegrass accessions, while 50 were identified as polymorphic markers and were then used to assess the level of genetic diversity among 43 centipedegrass core collections. The genetic diversity analysis exhibited that the number of alleles (Na) per locus ranged from 3 to 13, and the observed heterozygosity (Ho) ranged from 0.17 to 0.83. The polymorphism information content (PIC) value of the markers ranged from 0.15 to 0.78, and the genetic distances (coefficient Nei72) between the accessions varied from 0.07 to 0.48. The UPGMA-based dendrogram clustered all 43 core collections into two main groups and six subgroups, which further validated the effectiveness of these newly developed SSR markers. Hence, these newly developed SSR markers will be valuable and potentially useful for future genetic and genomic studies of E. ophiuroides.

Dehydration induced transcriptomic responses in two Tibetan hulless barley (Hordeum vulgare var. nudum) accessions distinguished by drought tolerance

Background

The harsh environment on the Qinghai-Tibetan Plateau gives Tibetan hulless barley (Hordeum vulgare var. nudum) great ability to resist adversities such as drought, salinity, and low temperature, and makes it a good subject for the analysis of drought tolerance mechanism. To elucidate the specific gene networks and pathways that contribute to its drought tolerance, and for identifying new candidate genes for breeding purposes, we performed a transcriptomic analysis using two accessions of Tibetan hulless barley, namely Z772 (drought-tolerant) and Z013 (drought-sensitive).

Results

There were more up-regulated genes of Z772 than Z013 under both mild (5439-VS-2604) and severe (7203-VS-3359) dehydration treatments. Under mild dehydration stress, the pathways exclusively enriched in drought-tolerance genotype Z772 included Protein processing in endoplasmic reticulum, tricarboxylic acid (TCA) cycle, Wax biosynthesis, and Spliceosome. Under severe dehydration stress, the pathways that were mainly enriched in Z772 included Carbon fixation in photosynthetic organisms, Pyruvate metabolism, Porphyrin and chlorophyll metabolism. The main differentially expressed genes (DEGs) in response to dehydration stress and genes whose expression was different between tolerant and sensitive genotypes were presented in this study, respectively. The candidate genes for drought tolerance were selected based on their expression patterns.

Conclusions

The RNA-Seq data obtained in this study provided an initial overview on global gene expression patterns and networks that related to dehydration shock in Tibetan hulless barley. Furthermore, these data provided pathways and a targeted set of candidate genes that might be essential for deep analyzing the molecular mechanisms of plant tolerance to drought stress.

Electronic supplementary material

The online version of this article (10.1186/s12864-017-4152-1) contains supplementary material, which is available to authorized users.