Development of SSR Markers Based on Transcriptome Sequencing and Association Analysis with Drought Tolerance in Perennial Grass Miscanthus from China

Association Analysis of Tiller-Related Traits with EST-SSR Markers in Psathyrostachys juncea

Psathyrostachys juncea is a long-lived perennial Gramineae grass with dense basal tillers and soft leaves. It is used widely in cold and dry areas of Eurasia and North America to establish grazing pasture and is even used as an ideal plant for revegetation and ecological restoration. Plant architecture, especially tillering traits, is critical for bunch grasses in breeding programs, and these traits in plants are mostly quantitative traits. In this study, the genetic diversity, population structure, and linkage disequilibrium of 480 individual lines were analyzed using 127 pairs of the EST-SSR marker, and a significant association between ten plant-architecture-related traits of P. juncea and molecular markers was found. The results of the genetic diversity analysis showed that the number of observed alleles was 1.957, the number of effective alleles was 1.682, Shannon's information index was 0.554, observed heterozygosity was 0.353, expected heterozygosity was 0.379, and the polymorphism information content was 0.300. A total of 480 individual lines were clustered into five groups based on population genetic structure, principal coordinate analysis, and unweighted pair group method with arithmetic mean analysis (UPGMA). The linkage disequilibrium coefficient (r2) was between 0.00 and 0.68, with an average of 0.04, which indicated a relatively low level of linkage disequilibrium among loci. The results of the association analysis revealed 55 significant marker-trait associations (MTA). Moreover, nine SSR markers were associated with multiple traits. This study provides tools with promising applications in the molecular selection and breeding of P. juncea germplasm.

Transcriptional Regulation of Quinoa Seed Quality: Identification of Novel Candidate Genetic Markers for Increased Protein Content

Quinoa (Chenopodium quinoa Willd.) is a crop that has great potential for increased cultivation in diverse climate regions. The seed protein quality obtained from this crop is high concerning the requirements to meet human nutritional needs, but the seed protein content is relatively low if compared to crops such as grain legumes. Increased seed protein content is desirable for increasing the economic viability of this crop in order for it to be used as a protein crop. In this study, we characterized three genotypes of quinoa with different levels of seed protein content. By performing RNA sequencing of developing seeds, we determined the genotype differences in gene expression and identified genetic polymorphisms that could be associated with increased protein content. Storage nutrient analyses of seeds of three quinoa genotypes (Titicaca, Pasankalla, and Regalona) from different ecoregions grown under controlled climate conditions showed that Pasankalla had the highest protein content (20%) and the lowest starch content (46%). Our seed transcriptome analyses revealed highly differentially expressed transcripts (DETs) in Pasankalla as compared to the other genotypes. These DETs encoded functions in sugar transport, starch and protein synthesis, genes regulating embryo size, and seed transcription factors. We selected 60 genes that encode functions in the central carbon metabolism and transcription factors as potential targets for the development of high-precision markers. Genetic polymorphisms, such as single nucleotide polymorphisms (SNPs) and base insertions and deletions (InDels), were found in 19 of the 60 selected genes, which can be further evaluated for the development of genetic markers for high seed protein content in quinoa. Increased cultivation of quinoa can contribute to a more diversified agriculture and support the plant protein diet shift. The identification of quinoa genotypes with contrasting seed quality can help establish a model system that can be used for the identification of precise breeding targets to improve the seed quality of quinoa. The data presented in this study based on nutrient and transcriptome analyses contribute to an enhanced understanding of the genetic regulation of seed quality traits in quinoa and suggest high-precision candidate markers for such traits.

Genome-Wide Investigation of the NAC Transcription Factor Family in Miscanthus sinensis and Expression Analysis Under Various Abiotic Stress

The NAC transcription factor family is deemed to be a large plant-specific gene family that plays important roles in plant development and stress response. Miscanthus sinensis is commonly planted in vast marginal land as forage, ornamental grass, or bioenergy crop which demand a relatively high resistance to abiotic stresses. The recent release of a draft chromosome-scale assembly genome of M. sinensis provided a basic platform for the genome-wide investigation of NAC proteins. In this study, a total of 261 M. sinensis NAC genes were identified and a complete overview of the gene family was presented, including gene structure, conserved motif compositions, chromosomal distribution, and gene duplications. Results showed that gene length, molecular weights (MW), and theoretical isoelectric points (pI) of NAC family were varied, while gene structure and motifs were relatively conserved. Chromosomal mapping analysis found that the M. sinensis NAC genes were unevenly distributed on 19 M. sinensis chromosomes, and the interchromosomal evolutionary analysis showed that nine pairs of tandem duplicates genes and 121 segmental duplications were identified, suggesting that gene duplication, especially segmental duplication, is possibly associated with the amplification of M. sinensis NAC gene family. The expression patterns of 14 genes from M. sinensis SNAC subgroup were analyzed under high salinity, PEG, and heavy metals, and multiple NAC genes could be induced by the treatment. These results will provide a very useful reference for follow-up study of the functional characteristics of NAC genes in the mechanism of stress-responsive and potential roles in the development of M. sinensis.

Transcriptome characterization of candidate genes related to chromium uptake, transport and accumulation in Miscanthus sinensis

Miscanthus sinensis is a C4 perennial grass species that is widely used as forage, ornamental grass and bioenergy crop due to its broad adaption and great biological traits. Recent studies indicated that M. sinensis could also grow in marginal lands which were contaminated with heavy metals, and exhibited important ecological restoration potential. In this study, transcriptome characterization of candidate genes related to chromium (Cr) uptake, transport and accumulation in M. sinensis were employed to investigate the molecular mechanism of plant tolerance to heavy metal stress. The result showed that following treatment of 200 mg/L of Cr, plant roots could accumulate most Cr and localize mainly in cell walls and soluble fractions, whereas Cr in stems and leaves was primarily in soluble fractions. A total of 83,645 differentially expressed genes (DEGs) were obtained after the treatment. Many genes involved in heavy metal transport, metal ion chelation and photosynthesis were found to be Cr-induced DEGs. Co-expression and weighted correlation network analysis revealed that Glutathion metabolism and ABC transporters pathways play an important role in Cr tolerance of M. sinensis. A hypothesis schematic diagram for the Cr uptake, transport and accumulation of M. sinensis cells were suggested, which could provide a molecular and genetic basis for future candidate genes validation and breeding of such crops.

The development of SSR markers based on RNA-sequencing and its validation between and within Carex L. species

BACKGROUND

Carex L. is one of the largest genera in the Cyperaceae family and an important vascular plant in the ecosystem. However, the genetic background of Carex is complex and the classification is not clear. In order to investigate the gene function annotation of Carex, RNA-sequencing analysis was performed. Simple sequence repeats (SSRs) were generated based on the Illumina data and then were utilized to investigate the genetic characteristics of the 79 Carex germplasms.

RESULTS

In this study, 36,403 unigenes with a total length of 41,724,615 bp were obtained and annotated based on GO, KOG, KEGG, NR databases. The results provide a theoretical basis for gene function exploration. Out of 8776 SSRs, 96 pairs of primers were randomly selected. One hundred eighty polymorphic bands were amplified with a polymorphism rate of 100% based on 42 pairs of primers with higher polymorphism levels. The average band number was 4.3 per primer, the average distance value was 0.548, and the polymorphic information content was ranged from 0.133 to 0.494. The number of observed alleles (Na), effective alleles (Ne), Nei's (1973) gene diversity (H), and the Shannon information index (I) were 2.000, 1.376, 0.243, and 0.391, respectively. NJ clustering divided into three groups and the accessions from New Zealand showed a similar genetic attribute and clustered into one group. UPGMA and PCoA analysis also revealed the same result. The analysis of molecular variance (AMOVA) revealed a superior genetic diversity within accessions than between accessions based on geographic origin cluster and NJ cluster. What's more, the fingerprints of 79 Carex species are established in this study. Different combinations of primer pairs can be used to identify multiple Carex at one time, which overcomes the difficulties of traditional identification methods.

CONCLUSIONS

The transcriptomic analysis shed new light on the function categories from the annotated genes and will facilitate future gene functional studies. The genetic characteristics analysis indicated that gene flow was extensive among 79 Carex species. These markers can be used to investigate the evolutionary history of Carex and related species, as well as to serve as a guide in future breeding projects.

MicroRNA-Mediated Responses to Chromium Stress Provide Insight Into Tolerance Characteristics of Miscanthus sinensis

Chromium (Cr) is a heavy metal in nature, which poses a potential risk to toxicity to both animals and plants when releasing into the environment. However, the regulation of microRNA (miRNA)-mediated response to heavy metal Cr has not been studied in Miscanthus sinensis. In this study, based on high-throughput miRNA sequencing, a total of 104 conserved miRNAs and 158 nonconserved miRNAs were identified. Among them, there were 45 differentially expressed miRNAs in roots and 13 differentially expressed miRNAs in leaves. The hierarchical clustering analysis showed that these miRNAs were preferentially expressed in a certain tissue. There were 833 differentially expressed target genes of 45 miRNAs in roots and 280 differentially expressed target genes of 13 miRNA in leaves. After expression trend analysis, five significantly enriched modules were obtained in roots, and three significantly enriched trend blocks in leaves. Based on the candidate gene annotation and gene ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) function analysis, miR167a, novel_miR15, and novel_miR22 and their targets were potentially involved in Cr transportation and chelation. Besides, miR156a, miR164, miR396d, and novel_miR155 were identified as participating in the physiological and biochemical metabolisms and the detoxification of Cr of plants. The results demonstrated the critical role of miRNA-mediated responses to Cr treatment in M. sinensis, which involves ion uptake, transport, accumulation, and tolerance characteristics.

Molecular Markers Improve Abiotic Stress Tolerance in Crops: A Review

Plants endure many abiotic stresses, such as temperature (heat or frost), drought, and salt. Such factors are primary and frequent stressors that reduce agriculture crop yields. Often alterations in nutrient management and constituents, along with variations in biosynthetic capacity, ultimately reduce or halt plant growth. Genetically, stress is an environmental condition that interferes with complete genetic expression. A vast range of molecular genomic markers is available for the analysis of agricultural crops. These markers are classified into various groups based on how the markers are used: RAPD (Random amplified polymorphic DNA) markers serve to identify and screen hybrids based on salinity and drought stress tolerance, while simple sequence repeat (SSR) markers are excellent for the assessment of stress tolerance. Such markers also play an important role in the QTL (Quantitative trait loci) mapping of stress-related genes. Dehydrins for drought and saltol for salinity stresses are primitive genes which regulate responses to these conditions. Further, a focus on traits using single-gene single nucleotide polymorphisms (SNP) markers supports genetic mapping and the sequencing of stress-related traits in inbred lines. DNA markers facilitate marker-assisted breeding to enhance abiotic stress tolerance using advanced techniques and marker modification.

Development of EST-SSR markers and association mapping with floral traits in Syringa oblata

BACKGROUND

Lilac (Syringa oblata) is an important woody plant with high ornamental value. However, very limited genetic marker resources are currently available, and little is known about the genetic architecture of important ornamental traits for S. oblata, which is hindering its genetic studies. Therefore, it is of great significance to develop effective molecular markers and understand the genetic architecture of complex floral traits for the genetic research of S. oblata.

RESULTS

In this study, a total of 10,988 SSRs were obtained from 9864 unigene sequences with an average of one SSR per 8.13 kb, of which di-nucleotide repeats were the dominant type (32.86%, 3611). A set of 2042 primer pairs were validated, out of which 932 (45.7%) exhibited successful amplifications, and 248 (12.1%) were polymorphic in eight S. oblata individuals. In addition, 30 polymorphic EST-SSR markers were further used to assess the genetic diversity and the population structure of 192 cultivated S. oblata individuals. Two hundred thirty-four alleles were detected, and the PIC values ranged from 0.23 to 0.88 with an average of 0.51, indicating a high level of genetic diversity within this cultivated population. The analysis of population structure showed two major subgroups in the association population. Finally, 20 significant associations were identified involving 17 markers with nine floral traits using the mixed linear model. Moreover, marker SO104, SO695 and SO790 had significant relationship with more than one trait.

CONCLUSION

The results showed newly developed markers were valuable resource and provided powerful tools for genetic breeding of lilac. Beyond that, our study could serve an efficient foundation for further facilitate genetic improvement of floral traits for lilac.

Development of SSR markers based on transcriptome data and association mapping analysis for fruit shell thickness associated traits in oil palm (Elaeis guineensis Jacq.)

Present study mainly aimed to ascertain the distribution characteristics of gene-based microsatellite loci and to develop polymorphic SSR markers from the already available transcriptome data of Elaeis guineensis Jacq, an important oil crop. From this study, we identified the sum of 5791 SSRs across 51,425 unigenes from the transcripts of oil palm. We were able to evaluate 331primer pairs and characterized 183 polymorphic gene-based SSR markers. We identified a total of 506 alleles from the 183 polymorphic SSR loci, with an average of 2.77 alleles per locus. The characterized gene-based SSR markers from the transcriptome data of oil palm exhibited moderate levels of polymorphism with a significant level of heterozygosity ranges from 0.096 to 0.594 (mean = 0.336 ± 0.11). Among the identified SSR markers, sixty polymorphic markers were used to analyze genotypes of 55 oil palm accessions selected from three different provinces of China. Association mapping analysis provided the information of four markers that are associated with fruit shell thickness trait of oil palm. Among the four markers identified from association analysis, one SSR marker obtained from Unigene17150 is strictly associated with the oil palm fruit shell thickness trait.

Identification of candidate reference genes for quantitative RT-PCR in Miscanthus sinensis subjected to various abiotic stresses

Quantitative real-time PCR (qRT-PCR) has been widely used for studying gene expression at the transcript level. Its accuracy usually relies on the reference genes that are utilized for data normalization. Miscanthus sinensis, a perennial C4 grass with high biomass and strong resistance to adversities, is often utilized as a high value energy crop. However, no reliable reference genes have been investigated for normalizing gene expression for this species. In this study, 12 candidate reference genes were selected to identify their stability under five different abiotic stress treatments (drought, salt, cadmium, chromium and arsenic) by using geNorm, NormFinder, BestKeeper and RefFinder softwares. The results showed that 18S rRNA and Unigene33312 were the best reference genes under drought treatments. Unigene33312 and Unigene33024 were found to be the most stably expressed genes under salt stress and Cd stress. Moreover, Unigene33024 and PP2A were the most suitable reference genes under Cr stress and Unigene33024 and Sb09g019750 were deemed more suitable reference genes under As stress. In total, considering all the samples, Unigene33024 and PP2A were the most stable genes while ACTIN and Unigene26576 were the least stable reference genes for internal control. The expression patterns of two target genes (Cu/Zn SOD and CAT) were used to further verify those selected reference genes under different conditions. The results showed that the most and the least stable reference genes had clearly different expression patterns. This work comprehensively estimated the stability of reference genes in M. sinensis which may give insight to the reference genes selection in other tissues as well as other related varieties. These suggested reference genes would assist in further putative gene expression validation in M. sinensis.

Characteristics of Microsatellites Mined from Transcriptome Data and the Development of Novel Markers in Paeonia lactiflora

The insufficient number of available simple sequence repeats (SSRs) inhibits genetic research on and molecular breeding of Paeonia lactiflora, a flowering crop with great economic value. The objective of this study was to develop SSRs for P. lactiflora with Illumina RNA sequencing and assess the role of SSRs in gene regulation. The results showed that dinucleotides with AG/CT repeats were the most abundant type of repeat motif in P. lactiflora and were preferentially distributed in untranslated regions. Significant differences in SSR size were observed among motif types and locations. A large number of unigenes containing SSRs participated in catalytic activity, metabolic processes and cellular processes, and 28.16% of all transcription factors and 21.74% of hub genes for inflorescence stem straightness were found to contain SSRs. Successful amplification was achieved with 89.05% of 960 pairs of SSR primers, 55.83% of which were polymorphic, and most of the 46 tested primers had a high level of transferability to the genus Paeonia. Principal component and cluster dendrogram analyses produced results consistent with known genealogical relationships. This study provides a set of SSRs with abundant information for future accession identification, marker-trait association and molecular assisted breeding in P. lactiflora.

Transcriptome analysis of terpenoid biosynthetic genes and simple sequence repeat marker screening in Eucommia ulmoides

Trans-polyisoprene rubber is produced in the tissues of leaves, bark, and fruit of Eucommia ulmoides and is considered an important energy source. Transcript profiles of two tissues from E. ulmoides cv. Qinzhong No. 3, leaf and fruit, were analysed using the Illumina HiSeq 2000 system. In total, 104 million clean reads were obtained and assembled into 58,863 unigenes. Through gene functional classification, 28,091 unigenes (47.72%) were annotated and 65 unigenes have been hypothesized to encode proteins involved in terpenoid biosynthesis. In addition, 10,041 unigenes were detected as differentially expressed unigenes, and 29 of them were putatively related to terpenoid biosynthesis. The synthesis of trans-polyisoprene rubbers in E. ulmoides was hypothesised to be dominated by the mevalonate pathway. Farnesyl diphosphate synthase 2 (FPPS2) was considered a key component in the biosynthesis of trans-polyprenyl diphosphate. Rubber elongation factor 3 (REF3) might be involved in stabilising the membrane of rubber particles in E. ulmoides. To date, 351 simple sequence repeats (SSRs) were validated as polymorphisms from eight E. ulmoides plants (two parent plants and six F₁ individuals), and these could act as molecular markers for genetic map density increase and breeding improvement of E. ulmoides.

RNAseq analysis reveals drought-responsive molecular pathways with candidate genes and putative molecular markers in root tissue of wheat

Drought is one of the major impediments in wheat productivity. Traditional breeding and marker assisted QTL introgression had limited success. Available wheat genomic and RNA-seq data can decipher novel drought tolerance mechanisms with putative candidate gene and marker discovery. Drought is first sensed by root tissue but limited information is available about how roots respond to drought stress. In this view, two contrasting genotypes, namely, NI5439 41 (drought tolerant) and WL711 (drought susceptible) were used to generate ~78.2 GB data for the responses of wheat roots to drought. A total of 45139 DEGs, 13820 TF, 288 miRNAs, 640 pathways and 435829 putative markers were obtained. Study reveals use of such data in QTL to QTN refinement by analysis on two model drought-responsive QTLs on chromosome 3B in wheat roots possessing 18 differentially regulated genes with 190 sequence variants (173 SNPs and 17 InDels). Gene regulatory networks showed 69 hub-genes integrating ABA dependent and independent pathways controlling sensing of drought, root growth, uptake regulation, purine metabolism, thiamine metabolism and antibiotics pathways, stomatal closure and senescence. Eleven SSR markers were validated in a panel of 18 diverse wheat varieties. For effective future use of findings, web genomic resources were developed. We report RNA-Seq approach on wheat roots describing the drought response mechanisms under field drought conditions along with genomic resources, warranted in endeavour of wheat productivity.

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.