Next Article in Journal
Comparative Transcriptomics Reveals Clues for Differences in Pathogenicity between Hysterothylacium aduncum, Anisakis simplex sensu stricto and Anisakis pegreffii
Next Article in Special Issue
Climate Change and Green Sea Turtle Sex Ratio—Preventing Possible Extinction
Previous Article in Journal
Genome-Wide Analysis and Function Prediction of Long Noncoding RNAs in Sheep Pituitary Gland Associated with Sexual Maturation
Previous Article in Special Issue
Asymmetric Isolation and the Evolution of Behaviors Influencing Dispersal: Rheotaxis of Guppies above Waterfalls
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Genes
Volume 11
Issue 3
10.3390/genes11030322
genes-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Identification of Genic SSRs Provide a Perspective for Studying Environmental Adaptation in the Endemic Shrub Tetraena mongolica

by
Zhenhua Dang
1,†,
Lei Huang
1,†,
Yuanyuan Jia
1,
Peter J. Lockhart
2,
Yang Fong
2 and
Yunyun Tian
3,*
1
Inner Mongolia Key Laboratory of Grassland Ecology & Ministry of Education Key Laboratory of Ecology and Resource Use of the Mongolian Plateau, School of Ecology and Environment, Inner Mongolia University, Hohhot 010021, China
2
School of Fundamental Sciences, College of Sciences, Massey University, Palmerston North 4442, New Zealand
3
Ministry of Education Key Laboratory of Herbage & Endemic Crop Biotechnology, School of Life Sciences, Inner Mongolia University, Hohhot 010021, China
*
Author to whom correspondence should be addressed.
†
These authors contributed equally.
Genes 2020, 11(3), 322; https://0-doi-org.brum.beds.ac.uk/10.3390/genes11030322
Submission received: 1 February 2020 / Revised: 10 March 2020 / Accepted: 16 March 2020 / Published: 18 March 2020
(This article belongs to the Special Issue Rapid Evolution)
Browse Figures
Versions Notes

Abstract

:
Tetraena mongolica is a xerophytic shrub endemic to desert regions in Inner Mongolia. This species has evolved distinct survival strategies that allow it to adapt to hyper-drought and heterogeneous habitats. Simple sequence repeats (SSRs) may provide a molecular basis in plants for fast adaptation to environmental change. Thus, identifying SSRs and their possible effects on gene behavior has the potential to provide valuable information for studies of adaptation. In this study, we sequenced six individual transcriptomes of T. mongolica from heterogeneous habitats, focused on SSRs located in genes, and identified 811 polymorphic SSRs. Of the identified SSRs, 172, 470, and 76 were located in 5′ UTRs, CDSs, and 3′ UTRs in 591 transcripts; and AG/CT, AAC/GTT, and AT/AT were the most abundant repeats in each gene region. Functional annotation showed that many of the identified polymorphic SSRs were in genes that were enriched in several GO terms and KEGG pathways, suggesting the functional significance of these genes in the environmental adaptation process. The identification of polymorphic genic SSRs in our study lays a foundation for future studies investigating the contribution of SSRs to regulation of genes in natural populations of T. mongolica and their importance for adaptive evolution of this species.

1. Introduction

Tetraena mongolica (2n = 2x = 28) [1], a monotypic species of the Zygophyllaceae family, is a xerophytic shrub that originated from the Tethys Ocean, and is now endemic to the Eastern Alxa–Western Ordos desert region in Inner Mongolia, China [2,3]. The area is characterized by hyper-drought conditions, low annual temperatures, and high soil salinity [4]. To adapt to the semi-desert environment, T. mongolica has evolved distinct survival strategies, such as an extensive root system, thick cuticles on its stems, and succulent leaves. T. mongolica has important ecological roles as a windbreak and in the fixation of local ecosystem soils [3]. Moreover, high levels of triacylglycerol in its stems makes it combustible even in the fresh state, so it has been used as firewood and is named locally as “oil firewood” [5,6]. T. mongolica is a typical species that grows in this area and is a good candidate for studies of plant origin, evolution, and adversity adaptation in such desert regions. Around these topics, a series of fundamental research investigations have been implemented, the results of which are of importance for comprehensively understanding the properties of T. mongolica. This pioneering work includes the studies of systematic, geographic, and taxonomic aspects [7,8,9]; morphology and anatomy [10,11,12]; physiological ecology [13,14,15,16]; chemometrics and biochemistry [17,18]; population genetics [5]; conservation biology [19,20,21]; abiotic stress response [3]; and medicinal phytochemistry [22,23]. In recent decades, the plant has suffered from serious habitat destruction in size and has even disappeared, leading to the general decline of the species and making the survival of T. mongolica much harder than before [24]. Correspondingly, one consideration for the ongoing survival of the species that has been paid much attention recently is an apparent decrease in genetic diversity due to loss of demographics, inbreeding depression, stochastic events, and other factors [5,24]. However, contrary to the assumption of diversity loss, it has been observed that the species has an intermediate or high level of genetic diversity as detected from analyses of inter-simple sequence repeats and nuclear microsatellites, respectively [5,24]. Therefore, further exploration of these hereditary characteristics and a potential link to adaptation strategies of the species is needed.
Genic simple sequence repeats (SSRs) (also called genic microsatellites or expressed sequence tag (EST)-SSRs), a relatively newly established concept, are simple sequence repeats located in genes, including protein-coding genes (CDSs) and their untranslated regions (UTRs) [25]. With the rapid accumulation of genetic data from expressed sequence tags (ESTs), a vast number of genic microsatellites have been identified and studied in many plants, especially for non-model species, which have a limited genomic background. Several outstanding characteristics of genic SSRs have led to these markers being widely used for the study of germplasm characterization [26], genetic mapping and breeding studies [27], and functional genome analysis [28]. The merits of these molecular resource include their high quality and relative ease-of-use for comparisons between genome counterparts, the high frequency and non-random distribution within transcribed regions of the genomes, and the transferability of markers to related species or genera [25]. It has been revealed that variations of these SSRs may change gene activity or regulate gene expression, thus affecting biological and cellular processes, such as protein structure, sensing and signaling, and gene transcription [29]. For example, CGG/CCG repeats have been associated with DNA methylation, indicating that this repeat may be the trigger for gene silencing [30]. ATTCT/AGAAT repeats can readily form unpaired regions in supercoiled plasmids, which have roles in the initiation of replication [31]. The GAA/TTC repeats form triplexes and related structures known as “sticky DNA”, which trap RNA polymerase on the template, thus blocking transcription elongation [32,33]. Therefore, genic SSRs have been proposed as evolutionary tuning knobs that may provide a molecular basis for fast adaptation to environmental changes in organisms [34]. Identification and in-depth study of these genomic elements, especially those in natural populations that inhabit ecologically heterogeneous and stressful environments, can provide new insights into the adaptive evolutionary mechanisms of plants [35,36].
Next-generation sequencing technologies have revolutionized studies in molecular ecology [37]. These technologies and relevant bioinformatics tools have helped in the identification of many genetic variations and provided insights into adaptive evolution, population genomics, and conservation genomics [37,38]. In this study, we sequenced, assembled, and compared the transcriptomes of six individual T. mongolica plants collected from heterogeneous habitats, and identified, validated, and analyzed the characteristics of polymorphic genic SSRs in this species. Our results will facilitate further studies into the functions of the identified SSRs, adaptive evolution mechanisms of T. mongolica, and the protection and exploration of the germplasm resources of this species.

2. Materials and Methods

2.1. Plant Materials

Leaves of T. mongolica were collected from six natural populations in June 2015 from the Eastern Alxa–Western Ordos area in Inner Mongolia, China. Twenty individuals from each population were selected at intervals of 30 m to avoid sampling progeny from the same maternal parent. Detailed geographical information about the sampling location is given in Table 1. After collection, the leaf samples were snap frozen in liquid nitrogen, then kept in a −80 °C refrigerator until they were transported to the laboratory.

2.2. RNA Preparation, cDNA Library Construction, and Sequencing

T. mongolica leaves (one sample per population) were used for RNA isolation. Total RNA was extracted using a Plant Plus RNA Kit (DP437, Tiangen, Beijing, China), according to the manufacturer’s instructions. The extracted RNA was treated with deoxyribonuclease I (TaKaRa Bio Inc., Otsu, Shiga, Japan) for 30 min at 37 °C to remove residual DNA. RNA quality was monitored using an Agilent 2100 Bioanalyzer with a minimum RNA integrity number value of 8.0. Six cDNA libraries (D1–D6) were constructed and sequenced on a paired-end flow cell device using an Illumina HiSeqTM 4000 sequencer at Biomarker Biotechnology Co., Ltd (Beijing, China). The sequencing data were called and quality controlled through the Illumina data processing pipeline.

2.3. De novo Assembly and Sequence Clustering

The adaptor reads, low-quality reads, and reads in which the proportion of “N” was greater than 5% were filtered to obtain the clean reads. De novo transcriptome assembly was performed for each sample using Trinity software (v2.4.0). TGICL v2.1 (http://sourceforge.net/projects/tgicl/files/tgicl%20v2.1/) was used to cluster the homologous transcripts and form an unextendible unigene. Finally, to eliminate redundancy and splicing among the six assembled transcriptomes, we further merged the unigenes to obtain non-redundant unigenes (all-unigenes).

2.4. Identification of Potential Polymorphic Genic SSRs

CandiSSR was used to identify potential polymorphic genic SSRs among the six unigenes datasets of T. mongolica [39]. The parameters used in the pipeline were: flanking sequence length 100 bp, blast e-value cutoff 1.0 × e−10, blast identity cutoff 95%, and blast coverage cutoff 95% [39,40]. The SSR search criteria were designed for di-, tri-, tetra-, penta-, and hexa-nucleotide repeats, with a minimum repeat number of six, five, five, four, and four, respectively, while SSRs with mono-nucleotide repeats and motif bigger than six nucleotides were not examined in this study. Primers for the identified SSRs were designed automatically in the pipeline using Primer3 [41].

2.5. Location Prediction and Functional Annotation of the Identified SSRs

We first determined the protein-coding sequences (CDSs) of the SSR-containing sequences. The sequences were aligned to protein databases in the priority order of NR, Swiss-Prot, KEGG, and COG using BLASTX (http://blast.ncbi.nlm.nih.gov/Blast.cgi), with a significance threshold e-value <1.0 × 10−5. The highest ranked proteins in blast results were used to predict the CDSs of the SSR-containing sequences. The results were then used to identify the distributions of the SSRs. For unigenes with full-length CDSs, the locations of the SSRs were analyzed based on the relative positions of the SSRs with respect to the start (ATG) and stop (TAA, TAG, TGA) codons. For unigenes with partial CDSs, we manually performed nucleotide blast analysis to identify full-length homologous genes and then predicted the locations of the SSRs according to their positions in the query sequences. To obtain functional information about the genic-SSR-containing unigenes, gene ontology (GO) functional categories were identified using the Blast2GO software v2.5.0 (http://www.blast2go.com/b2ghome) and pathways were assigned by sequence searches against the KEGG database using the BLASTX algorithm, with an e-value threshold of 1.0 × 10−5.

2.6. SSR Polymorphic Validation

Genomic DNA was extracted from T. mongolica leaves (eight individuals per population) using a Plant Genomic DNA kit (DP305, TianGen, Beijing, China). DNA quality was checked on 1% TAE agarose gels and NanoDrop 2000c (Gene Company Limit, Beijing, China). PCR amplifications were performed on the ABI2720 Thermal cycler in 25 μl reaction mixtures that included 1 μl of template DNA (50 ng/μL), 12.5 μL of Premix Taq (TaKaRa Biotechnology Co., Dalian, Liaoning Province, China), 0.5 μl (10 pM) of forward primer, 0.5 μl (10 pM) of reverse primer, and ddH2O. PCR amplifications were carried out as follows: 5 min at 94 °C, followed by 30 cycles of 30 s at 94 °C, 30 s at the primer-specific annealing temperature (Table 2), 30 s at 72 °C, and the final extension step at 72 °C for 10 min. For the successfully amplified primer pairs, the 5′ end of each forward primer was tagged with one of three fluorescent dyes (6-carboxy-fluorescine, hexachloro-6-carboxy-fluorescine, and 6-carboxy-X-rhodamine) and used for amplifications with the same protocol. The labeled PCR products were analyzed on an ABI 3730xl DNA analyzer with a GeneScan 500 LIZ Size Standard (Applied Biosystems, Beijing, China). Allele sizes were called using GeneMarker v2.6.0 (SoftGenetics, State College, Pennsylvania, USA). Several genetic parameters, including mean number of alleles per locus (NA) and observed and expected heterozygosity (HO and HE), were estimated across all the identified loci using GenAlEx software (v6.502) [42]. Polymorphism information content (PIC) of each locus was measured by PowerMarker program (v3.0) [43].

2.7. Data Availability

The raw sequencing data have been deposited in the Sequence Read Archive under accession number PRJNA597711.

3. Results

3.1. Sequencing Outputs and de novo Assembly

A total of 51.48 gigabases of clean data were obtained from the six cDNA libraries. As shown in Table 2, the clean reads in all samples had an average CG content of 44.24%, and the quality of these data were quite high, with an average Q20 percentage nearly of 97.15%. After de novo assembly, 73,977 to 92,301 unigenes in the six datasets were generated, and the mean length and N50 were 806 bp and 1536 bp, respectively. Finally, the assembled unigenes of the six transcriptomes were further clustered into a total of 119,603 all-unigene sequences, which had an average length of 1098 bp and average N50 value of 1843 bp.

3.2. Identification of Polymorphic SSRs

The CandiSSR analysis detected 811 polymorphic SSRs that were present in almost every dataset. The polymorphic levels of each SSR in these datasets were random (Table S1). The identified SSRs represented 89 motif types in 646 of the all-unigene sequences. The most common SSRs were tri-nucleotide repeats (522, 64.36%), followed by di-nucleotide repeats (261, 32.18%); quad-, penta-, and hexa-nucleotide repeats were far less frequent. Among the tri-nucleotide repeats, AAC/GTT (81, 15.52%) was the most abundant, followed by AGC/GCT (52, 9.96%), ACT/AGT (48, 9.20%), and AAG/CTT (43, 8.24%). Among the di-nucleotide repeats, AT/AT (119, 45.60%) was the most abundant, followed by AG/CT (99, 37.93%) and AC/GT (43, 16.48%).

3.3. Location Prediction and Frequency Analysis of the Polymorphic SSRs

Among the polymorphic SSR-containing sequences, 620 were aligned with known proteins and 763 SSRs were embedded in these transcripts. Of these SSRs, 718 were located in the coding and untranslated regions of 591 unigenes. Among them, 172, 470, and 76 SSRs in 150, 372, and 69 unigenes were distributed in 5′ UTRs, coding regions, and 3′ UTRs, respectively (Table S2). Most of the tri-nucleotide repeats (90.12%) were located in CDSs, while most of the di-nucleotide repeats (91.33%) were located in UTRs. In the CDSs, AAC/GTT, AGC/CTG, AAG/CTT, ATC/ATG, AAG/CTT, ACC/GGT, and ATC/GAT were relatively abundant, together accounting for 62.55% of all the SSRs located in CDSs. AG/CT (62, 36.05%) was the most abundant repeat in 5′ UTRs and AT/AT (42, 55.26%) was the most abundant repeat in 3′ UTRs (Figure 1).

3.4. Functional Annotation of the Polymorphic SSR-Containing Sequences

Based on the sequence similarity, the 424 unigenes containing 525 SSRs were mapped to 3561 terms in the GO database by blast. The terms under the three main GO categories were divided among 45 sub-categories. For 5′ UTRs that contained SSRs, the most abundant terms in each main categories were “metabolic process”, “cell”, and “binding”, respectively. For CDSs that contained SSRs, “cellular process” was the most abundant term, followed by “metabolic process”, “cell”, “cell part”, and “single-organism process”. For the 3′ UTRs that contained SSRs, the top three terms were “metabolic process”, “cellular process”, and “single-organism process” (Figure 2, Table S3).
Pathway analysis assigned 228 genic-SSR-containing sequences to 77 KEGG pathways (Figure 3, Table S4). For the 5′ UTRs, CDSs, and 3′ UTRs that contained SSRs, “metabolic pathways” was the most common and abundant term, followed by “plant hormone signal transduction”, “endocytosis”, and “biosynthesis of secondary metabolites” in each gene region.

3.5. Validation of the Polymorphic SSRs

Primer pairs were successfully designed for 795 identified SSRs using their flanking sequences (Table S5). To validate the effectiveness of these primers, 20 of the identified polymorphic SSRs were selected randomly and validated by PCR. Seventeen SSRs were successfully amplified with polymorphic loci. The loci size polymorphism levels were further assessed by capillary electrophoresis (Figure 4). A total of 76 alleles were detected with 2 to 7 (average 4.471) alleles obtained per locus (Table 3). The expected heterozygosity (HE), ranged from 0.274 to 0.774 and the observed heterozygosity (HO) from 0.167 to 0.729, with mean values of 0.587 and 0.418, respectively. The polymorphism information content (PIC) value ranged from 0.240 to 0.728, with an average value of 0.531 (Table 3).

4. Discussion

4.1. Identifying and General Profiling of Genic SSRs in T. mongolica Transcriptomes

SSRs have been recognized as an efficient tool for linking phenotypic and genotypic variations. They have been used widely in plants, for linkage map development, quantitative trait loci mapping, marker-assisted selection, parentage analysis, cultivar fingerprinting, genetic diversity studies, gene flow, and evolutionary studies [25,37]. Although microsatellite markers are versatile, developing polymorphic and robust SSR markers has been the main limitation for their use, especially for species that lack genomic background information. In recent years, next-generation sequencing technologies and a variety of computer programs (e.g., MISA [45], SSRFinder [46], SSRIT [47]) have promoted the identification of SSRs in model and non-model organisms [37]. However, the downstream development and identification of usable and researchable polymorphic SSRs is still a manual and time-consuming process [39]. Recently, an integrated procedure called CandiSSR was developed and used to obtain polymorphic SSRs in silico using genomic or transcriptomic sequence data [40]. In the present study, using the bioinformatics pipeline, we identified 811 potential polymorphic SSRs in six sequenced transcriptomes of T. mongolica. By comparing our finding with some other SSR development projects performed using CandiSSR, we found a similar polymorphic genic SSRs generation rate among these different studies. For example, 497 polymorphic genic SSRs were identified in the comparison of the transcriptomes representing the northern and southern population of Parrotia subaequalis [40]; and 1663 polymorphic EST-SSRs were identified by comparing the transcriptomes of tea plant and 19 Camellia species [48]. The identified SSRs in these studies were relatively small. One reason was probably attributed to the search criterial used in CandiSSR (SSRs with mono-nucleotide repeats and motif bigger than six nucleotides were excluded). Another reason was that the pipeline was designed to identify polymorphic SSRs directly based on the transcriptome or genome sequences comparison, thus the SSRs number would be significantly lower when compared to those SSR mining projects that were performed using only one dataset. Experimental validation showed that 17 of 20 randomly selected SSRs had 2 to 7 (average 4.471) polymorphic alleles per locus. The average HE and HO were 0.587 and 0.418, respectively (Table 3 and Figure 4). Although these genic SSRs may have low levels of polymorphism compared with the levels in genome SSRs (mean number of alleles per locus was 26.0, and mean HE and HO were 0.868 and 0.840, respectively [24]) in the T. mongolica genome, the genic SSRs have some intrinsic advantages (e.g., high transferability and functional importance) that make them a valuable resource that can be used for genetic and evolutionary studies of the species or its phylogenetically closed species and genus, such as species Larrea tridentata and genus Viviania marifolia [8,9]. In addition, Zhi et al. (2018) indicated that the highly polymorphic nuclear microsatellites might be over-estimating the genetic diversity of the present demography of T. mongolica [24], thus we considered that a rational genetic variation reserve could further be evaluated in T. mongolica by combing genetic data obtained from either the nuclear or genic SSRs. However, analyses of gene regions rather than non-coding genome regions evenly revealed a high level of genetic diversity in this endangered species [5,24], suggesting that the genome of the species is highly heterozygous to some extent and saturated with plenty of genetic variations that are potentially important for environmental adaptation and ongoing survival of the species. For the other three unamplified SSRs, we have speculated several possible reasons that might explain failure of the validations. Firstly, the PCR conditions might not have been suitable enough for amplifying the three SSRs. Secondly, the specificity of the primers was low. Thirdly, the primer binding sites may cross the boundary between introns and exons of the SSR-containing genes. Overall, whilst some markers could not be validated for technical reasons, we nevertheless were able to conclude that a large number of the identified SSRs in the present study are probably polymorphic and suitable for further study in T. mongolica transcripts.
Strong evidence has shown that SSRs are non-randomly distributed across CDSs and UTRs [29]. In many species, exons contain many more tri- and hexa-nucleotide SSRs than di- and tetra-nucleotide SSRs [49,50]. This distribution bias may be explained by the low number of non-trimeric SSRs found in coding regions, possibly caused by frameshift mutations [51], that might be related to features of adaptation or responses to environment stress, as was previously found in Escherichia coli K12 [35]. In T. mongolica, the distribution analysis of transcribed sequences showed that > 90.12% of the identified tri-nucleotide SSRs were located in coding regions, whereas 91.33% of the di-nucleotide SSRs were distributed in UTRs (Figure 1). The SSRs in transcribed regions also showed bias to a specific nucleotide composition [29]. It was reported that the UTRs, in particular the 5′ UTRs, were strongly biased toward AG/CT in Arabidopsis thaliana [49,52]. In this study, almost all the identified AG/CT repeats were distributed in 5′ UTRs and most of the AT/AT repeats were located in 3′ UTRs. The SSRs in coding regions also exhibited a strong bias toward specific triplet repeats, such as AAG, AAC, ATC, AGC, AGG, and ACG, and these triplet repeats are known to be the most common repeats in the exons of dicotyledonous plants [53]. As shown in Figure 1, our results are consistent with these observations, so the identified SSRs in T. mongolica may have similar functions to their counterparts.

4.2. Probable Functions of Genic SSRs in the Environmental Adaptation of T. mongolica

Habitat heterogeneity and environmental pressures cause various stress responses that involve changes in gene transcription and translation in multiple cellular processes and metabolic pathways, and plants have evolved effective and efficient strategies to adapt to environmental changes for their survival [54,55]. Although narrowly distributed in the Ordos desert region, the eco-physiology characteristics of T. mongolica populations in their heterogeneous habitats (e.g., upland, hill, tableland, high plain) are different. Previous studies have shown that T. mongolica populations have different water parameters (water potential in saturated point, water potential in turgor loss point, available water capacity), endogenous hormones (indole-3-acetic acid, abscisic acid, gibberellin, zeatin riboside), antioxidant activities (superoxide dismutase, peroxidases, catalases, glutathione reductase), and photosynthesis, which suggests that the populations have evolved distinct adaptative mechanisms under heterogeneous growth conditions [14,15,16]. In recent decades, anthropogenic disturbances such as mineral exploitation, road building, overgrazing, and farming have led to habitat isolation, thus further accelerating the differentiation and survival adaptation of T. mongolica [19,20]. Genic SSRs have been proposed as fast-evolving genomic elements that might be of importance for adaptive evolution [34,36]. Therefore, it is reasonable to assume that the genes that contain polymorphic SSRs are important drivers of a plant’s response to environmental pressure and fluctuations. In the present study, we identified 646 transcripts that contained 811 polymorphic SSRs that may play significant roles in the local adaptation processes of T. mongolica. Of the SSR containing sequences, the majority (620) returned at least one blast hit with the known protein database at the e-value > 1.0 × 10-5, while the remaining 26 unigenes did not match known genes, suggesting that these unannotated unigenes might represent unique genetic resources of T. mongolica. Substantial data indicate that SSR expansion or contraction may affect gene expression and activity, resulting in concomitant physiological consequences and observable phenotypic changes [29,55]. Variations in 5′ UTRs have been found to regulate gene expression by affecting transcription and translation. In coding regions, variations can lead to a gain or loss of gene function via frameshift mutations or expanded toxic mRNA, and in 3′ UTRs variations can cause transcription slippage and produce expanded mRNA, which can disrupt splicing and disrupt other cellular functions [25,29]. Variations in CT repeats in the 5′ UTR of the starch-synthesis genes in waxy rice (Oryza sativa) were correlated with the physicochemical properties of starch [56]. Changes in the number of TGA repeats in the cellulase synthase gene (PtoCesA10) in Populus tomentosa resulted in an insertion mutation in the coding region and the addition of aspartate to the polypeptide chain of the enzyme, which is uniquely associated with the lignin content [57]. Myotonic dystrophy type 1 is caused by CTG repeat expansions in the 3′ UTR of DMPK [58,59]. In our analysis, 150, 372, and 69 unigenes were found to contain 172, 470, and 76 SSRs in 5′ UTRs, coding regions, and 3′ UTRs, respectively (Table S2). The GO functional annotations showed that the identified SSR-containing sequences were enrichened in “cellular process”, “metabolic process”, “cell”, “cell part”, and “single-organism process” terms under the three main GO categories, suggesting such groups of T. mongolica genes might first be regulated by SSRs mutations under the heterogeneous environments (Figure 2, Table S3). The KEGG analysis showed that among the SSR-containing 5′ UTRs and CDSs, 35 were distributed in “plant hormone signal transduction” (Figure 3, Table S4), implying that the regulation of endogenous hormones may play significant roles in the local adaptation of T. mongolica. This prediction is partially supported by the eco-physiology observations described above [14,15,16]. In future studies, we plan to further investigate how changes in the SSR repeats in natural T. mongolica populations impact gene functions and clarify the probable functions of these SSRs, which would provide hints for the study of environmental adaptation of the species and information to answer the question of whether genic SSRs play significant roles in adaptive evolution.

Supplementary Materials

The following are available online at https://0-www-mdpi-com.brum.beds.ac.uk/2073-4425/11/3/322/s1: Table S1: Summary of the identified polymorphic SSRs by CandiSSR. Table S2: Summary of the location information of the identified polymorphic SSRs. Table S3: GO functional annotation of the polymorphic SSR containing sequences. Table S4: KEGG functional classification of the SSR containing sequences. Table S5: Primers for the polymorphic SSRs identified by CandiSSR.

Author Contributions

Project design, Z.D. and Y.T; sampling, Z.D., L.H., and Y.J.; wet lab work, L.H. and Y.J.; data submission, Y.T.; bioinformatic analyses, data analysis, writing, and review, Z.D., L.H, Y.T., P.J.L., and Y.F.; All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the National Natural Science Foundation of China (31860078).

Acknowledgments

We thank Margaret Biswas, from Liwen Bianji, Edanz Group China (www.liwenbianji.cn/ac), for editing the English text of a draft of this manuscript.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Guo, X.H.; Ci, Z.L.; Sun, J.; Yang, J. Chromosome variations during the tissue culture of Tetraena mongolia Maxim. J. Inn. Mong. Agric. Univ. (Nat. Sci. Edn.) 2001, 22, 55–59. [Google Scholar]
  2. Zhang, Y.J.; Yang, C. Comparative analysis of genetic diversity in the endangered shrub Tetraena mongolica and its related congener Zygophyllum xanthoxylon. Acta Phytoecol. Sin. 2000, 24, 425–429. [Google Scholar] [CrossRef]
  3. Chen, N.M.; Feng, J.C.; Song, B.; Tang, S.; He, J.Q.; Zhou, Y.J.; Shi, S.; Xu, X.J. De novo transcriptome sequencing and identification of genes related to salt and PEG stress in Tetraena mongolica Maxim. Trees 2019, 33, 1639–1656. [Google Scholar] [CrossRef]
  4. Dang, Z.H.; Zheng, L.L.; Wang, J.; Gao, Z.; Wu, S.B.; Qi, Z.; Wang, Y.C. Transcriptomic profiling of the salt-stress response in the wild recretohalophyte Reaumuria trigyna. BMC Genom. 2013, 14, 29. [Google Scholar] [CrossRef] [Green Version]
  5. Ge, X.J.; Yu, Y.; Zhao, N.X.; Chen, H.S.; Qi, W.Q. Genetic variation in the endangered Inner Mongolia endemic shrub Tetraena mongolica Maxim (Zygophyllaceae). Biol. Conserv. 2003, 111, 427–434. [Google Scholar] [CrossRef]
  6. Wang, G.L.; Lin, Q.Q.; Xu, Y.N. Tetraena mongolica Maxim can accumulate large amounts of triacylglycerol in phloem cells and xylem parenchyma of stems. Phytochemistry 2007, 68, 2112–2117. [Google Scholar] [CrossRef] [PubMed]
  7. Zhang, Y.J.; Yang, C.; Cheng, W.L. A comparison of ITS sequences of nuclear ribosomal DNA between endemic species Tetraena mongolica Maxim in Ordos Platean and three species of zygophyllaceae in the same habitat. Acta Sci. Natl. Univ. Neimong. 2003, 34, 420–424. [Google Scholar] [CrossRef]
  8. Ma, X.; Chang, J.Y.; Li, Z.H.; Zhai, W.; Yu, X.X.; Feng, Y.L. The complete chloroplast genome of Tetraena mongolica (Zygophyllaceae), an endangered shrub endemic to China. Mitochondrial DNA B 2019, 4, 1030–1031. [Google Scholar] [CrossRef] [Green Version]
  9. Duan, Y.Z.; Du, Z.Y.; Wang, H.T. Chloroplast genome characteristics of endangered relict plant Tetraena mongolica in the arid region of northwest China. Bull. Bot. Res. 2019, 39, 653–663. [Google Scholar] [CrossRef]
  10. Xu, Q.; Jiang, C.Q.; Liu, S.R.; Guo, Q.S. Study on Pollination ecology of endangered plant Tetraena mongolica population. Forest Res. 2003, 16, 391–397. [Google Scholar] [CrossRef]
  11. Jiang, S.; Ren, X.H.; Gu, S.; Li, J.; Xu, K. Wood anatomy of Tetraena mongolica Maxim (Zygophyllaceae). J. Trop. Subtrop. Bot. 2008, 16, 466–471. [Google Scholar] [CrossRef]
  12. Xu, Q.; Guo, Q.S.; Liu, S.R.; Jiang, C.Q.; Hao, Y.G. A study on the relationship between fruiting characteristics and reproductive age, habitat of endangered species Tetraena mongolica. Sci. Silv. Sin. 2003, 39, 26–32. [Google Scholar] [CrossRef]
  13. Li, X.; Wang, Y.C.; Zheng, R. Water parameters of xeric shrubs in west erdos region (I). Chin. J. Desert Res. 2005, 25, 581–586. [Google Scholar] [CrossRef]
  14. Shi, S.L.; Wang, Y.C.; Zhou, J.H. The change of water parameters in Tetraena mongolica with season and habitat. Acta Ecol. Sin. 2008, 28, 6079–6089. [Google Scholar] [CrossRef]
  15. Shi, S.L.; Wang, Y.C.; Zhou, J.H. The seasonal change of endogenous phytohormone and differentiation of populations in Tetraena mongolica in different habitats. Acta Ecol. Sin. 2009, 29, 2252–2262. [Google Scholar] [CrossRef]
  16. Shi, S.L.; Wang, Y.C.; Li, X.; Zhou, H.B. Variety of antioxidant system of Tetraena mongolica in different growth stages and habitats. Chin. J. Desert Res. 2012, 32, 771–779. [Google Scholar]
  17. Shi, G.R.; Ding, L.L.; Liu, Q.; TANG, S.A.; Duan, H.Q. Chemical constituents contained in Tetraena mongolica. China J. Chin. Mater. Med. 2012, 37, 1579–1580. [Google Scholar] [CrossRef]
  18. Zhi, Y.B.; Yang, C.; Wang, Z.L.; Yao, Y.P.; Gao, T.Y.; Hao, X.; Liu, J.P. Characteristic analysis of amino acid contents of Tetraena mongolica Maxim. Acta Sci. Natl. Univ. Neimongol 2005, 36, 306–312. [Google Scholar] [CrossRef]
  19. Zhi, Y.B.; Yang, C.; Wang, Z.S.; An, S.Q.; Wang, Z.L.; Li, H.L.; Su, Z.A.; Wang, Q. The endangered characteristics and mechanism of the endemic relict shrub Tetraena mongolica Maxim. Acta Ecol. Sin. 2008, 28, 767–776. [Google Scholar] [CrossRef]
  20. Ge, X.J.; Hwang, C.C.; Liu, Z.H.; Huang, C.C.; Huang, W.H.; Huang, K.H.; Wang, W.K.; Chiang, T.Y. Conservation genetics and phylogeography of endangered and endemic shrub Tetraena mongolica (Zygophyllaceae) in Inner Mongolia, China. BMC Genet. 2011, 12, 1. [Google Scholar] [CrossRef] [Green Version]
  21. Zhen, J.H.; Liu, G.H. Population structure characteristics of Tetraena mongolica in different habitats. Acta Ecol. Sin. 2008, 28, 1829–1841. [Google Scholar] [CrossRef]
  22. Wu, Z.G.; Wei, W.; Cheng, K.; Zheng, L.L.; Ma, C.M.; Wang, Y.C. Insecticidal activity of triterpenoids and volatile oil from the stems of Tetraena mongolica. Pestic. Biochem. Phys. 2020. [Google Scholar] [CrossRef]
  23. Wu, Z.G.; Wei, W.; Xu, H.Y.; Zheng, L.L.; Ma, C.M.; Wang, Y.C. Constituents from the leaves of Tetraena mongolica and their protective activity in HEK 293t cells damaged by CdCl2. J. Nat. Prod. 2019, 82, 2707–2712. [Google Scholar] [CrossRef]
  24. Zhi, Y.B.; Sun, Z.L.; Sun, P.; Zhao, K.; Guo, Y.N.; Zhang, D.J.; Zhang, B.W. How much genetic variation is stored in the endangered and fragmented shrub Tetraena mongolica Maxim? PeerJ 2018, 6, e5645. [Google Scholar] [CrossRef] [Green Version]
  25. Varshney, R.K.; Graner, A.; Sorrells, M.E. Genic microsatellite markers in plants: Features and applications. Trends Biotechnol. 2005, 23, 48–55. [Google Scholar] [CrossRef]
  26. Nachimuthu, V.V.; Muthurajan, R.; Duraialaguraja, S.; Sivakami, R.; Pandian, B.A.; Ponniah, G.; Gunasekaran, K.; Swaminathan, M.; Suji, K.K.; Sabariappan, R. Analysis of population structure and genetic diversity in rice germplasm using SSR markers: An initiative towards association mapping of agronomic traits in Oryza Sativa. Rice 2015, 8, 30. [Google Scholar] [CrossRef] [Green Version]
  27. Pan, L.; Huang, T.; Yang, Z.F.; Tang, L.; Cheng, Y.J.; Wang, J.P.; Ma, X.; Zhang, X.Q. EST-SSR marker characterization based on RNA-sequencing of Lolium multiflorum and cross transferability to related species. Mol. Breed. 2018, 38, 80. [Google Scholar] [CrossRef]
  28. Qu, J.B.; Huang, C.Y.; Zhang, J.X. Genome-wide functional analysis of SSR for an edible mushroom Pleurotus ostreatus. Gene 2016, 575, 524–530. [Google Scholar] [CrossRef] [Green Version]
  29. Li, Y.C.; Korol, A.B.; Fahima, T.; Nevo, E. Microsatellites within genes: Structure, function, and evolution. Mol. Biol. Evol. 2004, 21, 991–1007. [Google Scholar] [CrossRef]
  30. Usdin, K. The biological effects of simple tandem repeats: Lessons from the repeat expansion diseases. Genome Res. 2008, 18, 1011–1019. [Google Scholar] [CrossRef] [Green Version]
  31. Potaman, V.N.; Bissler, J.J.; Hashem, V.I.; Oussatcheva, E.A.; Lu, L.; Shlyakhtenko, L.S.; Lyubchenko, Y.L.; Matsuura, T.; Ashizawa, T.; Leffak, M.; et al. Unpaired structures in SCA10 (ATTCT)n·(AGAAT)n repeats. J. Mol. Biol. 2003, 326, 1095–1111. [Google Scholar] [CrossRef]
  32. Vetcher, A.A.; Napierala, M.; Wells, R.D. Sticky DNA: Effect of the polypurine·polypyrimidine sequence. J. Biol. Chem. 2002, 227, 39228–39234. [Google Scholar] [CrossRef] [Green Version]
  33. Potaman, V.N.; Oussatcheva, E.A.; Lyubchenko, Y.L.; Shlyakhtenko, L.S.; Bidichandani, S.I.; Ashizawa, T.; Sinden, R.R. Length-dependent structure formation in Friedreich ataxia (GAA)n·(TTC)n repeats at neutral pH. Nucleic Acids Res. 2004, 32, 1224–1231. [Google Scholar] [CrossRef] [Green Version]
  34. Kashi, Y.; King, D.G. Simple sequence repeats as advantageous mutators in evolution. Trends Genet. 2006, 22, 253–259. [Google Scholar] [CrossRef] [Green Version]
  35. Rocha, E.P.C.; Matic, I.; Taddei, F. Over-representation of repeats in stress response genes: A strategy to increase versatility under stressful conditions? Nucleic Acids Res. 2002, 30, 1886–1894. [Google Scholar] [CrossRef] [Green Version]
  36. Li, Y.C.; Korol, A.B.; Fahima, T.; Beiles, A.; Nevo, E. Microsatellites: Genomic distribution, putative functions and mutational mechanisms: A review. Mol. Ecol. 2002, 11, 2453–2465. [Google Scholar] [CrossRef]
  37. Ekblom, R.; Galindo, J. Applications of next generation sequencing in molecular ecology of non-model organisms. Heredity 2011, 107, 1–15. [Google Scholar] [CrossRef] [Green Version]
  38. Gayral, P.; Weinert, L.; Chiari, Y.; Tsagkogeorga, G.; Ballenghien, M.; Galtier, N. Next-generation sequencing of transcriptomes: A guide to RNA isolation in nonmodel animals. Mol. Ecol. Resour. 2011, 11, 650–661. [Google Scholar] [CrossRef]
  39. Xia, E.H.; Yao, Q.Y.; Zhang, H.B.; Jiang, J.J.; Zhang, L.P.; Gao, L.Z. CandiSSR: An efficient pipeline used for identifying candidate polymorphic SSRs based on multiple assembled sequences. Front. Plant Sci. 2016, 6, 1171. [Google Scholar] [CrossRef] [Green Version]
  40. Zhang, Y.Y.; Zhang, M.Y.; Hu, Y.M.; Zhuang, X.; Xu, W.Q.; Li, P.F.; Wang, Z.S. Mining and characterization of novel EST-SSR markers of Parrotia subaequalis (Hamamelidaceae) from the first Illumina-based transcriptome datasets. PLoS ONE 2019, 14, e0215874. [Google Scholar] [CrossRef]
  41. Untergasser, A.; Cutcutache, I.; Koressaar, T.; Ye, J.; Faircloth, B.C.; Remm, M.; Rozen, S.G. Primer3-new capabilities and interfaces. Nucleic Acids Res. 2012, 40, e115. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  42. Peakall, R.; Smouse, P.E. GenAlEx 6.5: Genetic analysis in Excel. Population genetic software for teaching and research-an update. Bioinformatics 2012, 28, 2537–2539. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  43. Liu, K.J.; Muse, S.V. PowerMarker: An integrated analysis environment for genetic marker analysis. Bioinformatics 2005, 21, 2128–2129. [Google Scholar] [CrossRef] [Green Version]
  44. Hall, T.A. BioEdit: A user-friendly biological sequence alignment editor and analysis program for Windows 95/98/NT. Nucl. Acids. Symp. Ser. 1999, 41, 95–98. [Google Scholar] [CrossRef]
  45. Thiel, T.; Michalek, W.; Varshney, R.K.; Graner, A. Exploiting EST databases for the development and characterization of gene-derived SSR-markers in barley (Hordeum vulgare L.). Theor. Appl. Genet. 2003, 106, 411–422. [Google Scholar] [CrossRef]
  46. Gao, L.F.; Tang, J.F.; Li, H.W.; Jia, J.Z. Analysis of microsatellites in major crops assessed by computational and experimental approaches. Mol. Breed. 2003, 12, 245–261. [Google Scholar] [CrossRef]
  47. Temnykh, S.; DeClerck, G.; Lukashova, A.; Lipovich, L.; Cartinhour, S.; McCouch, S. Computational and experimental analysis of microsatellites in rice (Oryza sativa L.): Frequency, length variation, transposon associations, and genetic marker potential. Genome Res. 2001, 11, 1441–1452. [Google Scholar] [CrossRef] [Green Version]
  48. Xia, E.H.; Li, F.D.; Tong, W.; Li, P.H.; Wu, Q.; Zhao, H.J.; Ge, R.P.; Li, R.P.; Li, Y.Y.; Zhang, Z.Z.; et al. Tea plant information archive: A comprehensive genomics and bioinformatics platform for tea plant. Plant Biotechnol. J. 2019, 17, 1938–1953. [Google Scholar] [CrossRef]
  49. Morgante, M.; Hanafey, M.; Powell, W. Microsatellites are preferentially associated with nonrepetitive DNA in plant genomes. Nat. Genet. 2002, 30, 194–200. [Google Scholar] [CrossRef]
  50. Cordeiro, G.M.; Casu, R.; McIntyre, C.L.; Manners, J.M.; Henry, R.J. Microsatellite markers from sugarcane (Saccharum spp.) ESTs cross transferable to erianthus and sorghum. Plant Sci. 2001, 160, 1115–1123. [Google Scholar] [CrossRef]
  51. Metzgar, D.; Bytof, J.; Wills, C. Selection against frameshift mutations limits microsatellite expansion in coding DNA. Genome Res. 2000, 10, 72–80. [Google Scholar] [CrossRef]
  52. Lawson, M.J.; Zhang, L.Q. Distinct patterns of SSR distribution in the Arabidopsis thaliana and rice genomes. Genome Biol. 2006, 7, R14. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  53. Sonah, H.; Deshmukh, R.K.; Sharma, A.; Singh, V.P.; Gupta, D.K.; Gacche, R.N.; Rana, J.C.; Singh, N.K.; Sharma, T.R. Genome-wide distribution and organization of microsatellites in plants: An insight into marker development in Brachypodium. PLoS ONE 2011, 6, e21298. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  54. Yotoko, K.S.C.; Dornelas, M.C.; Togni, P.D.; Fonseca, T.C.; Salzano, F.M.; Bonatto, S.L.; Freitas, L.B. Does variation in genome sizes reflect adaptive or neutral processes? New clues from Passiflora. PLoS ONE 2011, 6, e18212. [Google Scholar] [CrossRef] [Green Version]
  55. Gao, C.H.; Ren, X.D.; Mason, A.S.; Li, J.N.; Wang, W.; Xiao, M.L.; Fu, D.H. Revisiting an important component of plant genomes: Microsatellites. Funct. Plant Biol. 2013, 40, 645–661. [Google Scholar] [CrossRef]
  56. Bao, J.S.; Corke, H.; Sun, M. Microsatellites in starch-synthesizing genes in relation to starch physicochemical properties in waxy rice (Oryza sativa L.). Theor. Appl. Genet. 2002, 105, 898–905. [Google Scholar] [CrossRef] [PubMed]
  57. Du, Q.Z.; Pan, W.; Xu, B.H.; Li, B.L.; Zhang, D.Q. Polymorphic simple sequence repeat (SSR) loci within cellulose synthase (PtoCesA) genes are associated with growth and wood properties in Populus tomentosa. New Phytol. 2013, 197, 763–776. [Google Scholar] [CrossRef]
  58. Ranum, L.P.W.; Cooper, T.A. RNA-mediated neuromuscular disorders. Annu. Rev. Neurosci. 2006, 29, 259–277. [Google Scholar] [CrossRef]
  59. Ranum, L.P.W.; Day, J.W. Myotonic dystrophy: Clinical and molecular parallels between myotonic dystrophy type 1 and type 2. Curr. Neurol. Neurosci. Rep. 2002, 2, 465–470. [Google Scholar] [CrossRef]
Genes 11 00322 g001 550
Figure 1. Distribution analysis of the identified genic simple sequence repeats (SSRs). The y-axis indicates the numbers of identified genic SSRs; the x-axis indicates the distribution of the polymorphic genic SSRs and motif sequence types.
Figure 1. Distribution analysis of the identified genic simple sequence repeats (SSRs). The y-axis indicates the numbers of identified genic SSRs; the x-axis indicates the distribution of the polymorphic genic SSRs and motif sequence types.
Genes 11 00322 g001
Genes 11 00322 g002 550
Figure 2. Gene ontology (GO) functional annotation of the polymorphic SSR-containing sequences. The three lines of bubbles represent SSRs that were located in the 5′ untranslated regions (UTRs), protein-coding sequences, and 3′ UTRs, respectively. GO terms that contained unigenes more than or equal to ten in one of the gene regions are shown in the figure. The relative number of unigenes assigned to each term is indicated by the size of each bubble. The grey hollow circle indicates no unigene was assigned to the relevant GO term.
Figure 2. Gene ontology (GO) functional annotation of the polymorphic SSR-containing sequences. The three lines of bubbles represent SSRs that were located in the 5′ untranslated regions (UTRs), protein-coding sequences, and 3′ UTRs, respectively. GO terms that contained unigenes more than or equal to ten in one of the gene regions are shown in the figure. The relative number of unigenes assigned to each term is indicated by the size of each bubble. The grey hollow circle indicates no unigene was assigned to the relevant GO term.
Genes 11 00322 g002
Genes 11 00322 g003 550
Figure 3. KEGG functional classification of the SSR-containing sequences. The y-axis indicates the numbers of SSR-containing sequences enriched in KEGG pathways; the x-axis indicates the top ten enriched pathways assigned to the 5′ UTR, CDSs, and 3′ UTR SSR-containing sequences.
Figure 3. KEGG functional classification of the SSR-containing sequences. The y-axis indicates the numbers of SSR-containing sequences enriched in KEGG pathways; the x-axis indicates the top ten enriched pathways assigned to the 5′ UTR, CDSs, and 3′ UTR SSR-containing sequences.
Genes 11 00322 g003
Genes 11 00322 g004 550
Figure 4. Representative polymorphic SSRs identified by CandiSSR and capillary electrophoresis. (A) Five unigenes assembled in the D1–D5 cDNA libraries that correspond to a polymorphic SSR identified by CandiSSR [39]. Multiple sequence alignment was performed using the Bioedit software (v7.0.9) [44]. (B) The capillary electrophoresis results for the polymorphic SSRs in (A).
Figure 4. Representative polymorphic SSRs identified by CandiSSR and capillary electrophoresis. (A) Five unigenes assembled in the D1–D5 cDNA libraries that correspond to a polymorphic SSR identified by CandiSSR [39]. Multiple sequence alignment was performed using the Bioedit software (v7.0.9) [44]. (B) The capillary electrophoresis results for the polymorphic SSRs in (A).
Genes 11 00322 g004
Table
Table 1. Location information for six population of T. mongolica.
Table 1. Location information for six population of T. mongolica.
PopulationLongitude (E)Latitude (N)Altitude (m)HabitatsSoil Water Content (%) *
D1106°53′43″39°21′57″1212.5Foothills5.48
D2106°53′52″39°22′30″1185.6Tableland5.24
D3106°53′31″39°29′52″1216.5Foothills4.93
D4107°05′45″40°14′58″1150.5Tableland3.89
D5106°52′07″40°08′03″1036.9Sandy Land2.43
D6106°55′07″40°08′02″1049.5Piedmont Plain2.43
*Cited from [16].
Table
Table 2. Summary of sequencing and assembly results.
Table 2. Summary of sequencing and assembly results.
CR (No.)CN (nt)Q20 (%)GC (%)Ug (No.)ML (bp)N50 (bp)
D153,284,2547,992,638,10097.22%43.76%80,4097911499
D255,473,3868,321,007,90097.19%43.99%80,8298241579
D364,363,3729,654,505,80096.88%45.15%77,6417861516
D452,017,9547,802,693,10096.94%44.27%84,6738511600
D563,352,4309,502,864,50097.36%43.72%92,3017941534
D654,704,6888,205,703,20097.31%44.54%73,9777881489
All 119,60310981843
CR, CN, GC, and Ug represent clean read, clean nucleotide, GC content, and unigenes, respectively. Q20 represents the clean reads that had Phred-like quality scores at the Q20 level (an error probability of 1%). ML represents the mean length of assembled sequences and N50 indicates that 50% of the assembled bases were incorporated into sequences with a length of N50 or longer.
Table
Table 3. Characteristics of 17 validated microsatellites for T. mongolica.
Table 3. Characteristics of 17 validated microsatellites for T. mongolica.
Gene IDPS (5′–3′)RMASTa (°C)PFNAHOHEPIC
CL3279.Ct4F: GTAGTACTACTGCTGCATCGTATCCT
R: CAACCCTATCTTCATCATCATCG		TGC101–11354protoporphyrinogen oxidase, Vitis vinifera40.4360.5100.458
CL4993.Ct1F: ACTCCTCTCATCCATCCATTAAG
R: GGAGTTTAACGCTGTCATTGTG		TC102–11455-60.5000.6950.644
CL6111.Ct3F: TGGAGTCTGAAGGCAGTGAG
R: ACTTGAACTTCTTGATTCCACC		TAG118–12455-30.1670.5900.508
CL8609.Ct3F: GCATTAGAGGAGCGAATCGAAG
R: GCCTCGCTTCTCATTTCTCAAC		GA164–17459-60.4170.6170.584
CL8025.Ct4F: CATCGCCGCCTTTCATAGAC
R: GACGCTTAGAATTGGAAGATGATG		TC175–18155cyclin-dependent kinase G-2-like, Citrus sinensis30.3540.4720.422
Ug20261F: GGGGAAAGATGCTGTTATGGAG
R: TAGCATCCGAGCCACTACCAC		AGG186–19559-40.5000.6250.552
CL6305.Ct2F: CGCTTGCTTTAACGACGAACC
R: TGTGGTGGGTCGGATGATGTT		GCA176–18855serine/threonine-protein kinase RIO1-like isoform X1, Citrus sinensis50.2790.6860.621
CL7264.Ct1F: GTTGTGGCGGCGTAGTTTATG
R: AACTCGCAAACCAAGAGCATAAC		TG193–20558-50.3330.7740.728
CL9244.Ct2F: CTGAGATTTGTTGGTGGGTTTG
R: CCAGTATCTCCGAACCACCTCT		AGG373–38256Glutaredoxin 4 isoform 1, Theobroma cacao30.4380.5510.457
CL8609.Ct2F: GGAGCTGAATTAGAGCATTAGAGG
R: GAAATCTCTCTTGTTCAATCCACC		GA202–21255-60.4580.6630.626
Ug13288F: AGCATTACATTATCCCTTCCTCAC
R: CAGAGACGGTGTCGTATTGGA		TAA240–25855peptide chain release factor 1-like, Glycine max50.5330.7100.649
Ug19883F: GAGTTATGAATGACGCTACACGAG
R: GCCTGCTTTGCGTTTCTTC		TGC351–36055-30.3190.2740.240
Ug13288F: CATCGCCGCCTTTCATAGAC
R: GACGCTTAGAATTGGAAGATGATG		TAA242–26355peptide chain release factor 1-like, Glycine max70.5750.7580.713
Ug31697F: CAACAGAAAGCACCAACCCAG
R: GCATCCACCCTGTTCAGCAT		CTC241–25360-50.7290.7700.722
Ug4409F: CATCGGCCTCTGCTCATACAC
R: CGCTTCAGGCTCTCATATTCAG		TCA272–27555-20.2920.3990.317
CL12118.Ct3F: CAGAGAGAATAATAGCAGCCATAG
R: CCCAAGCATCCAACAATAAC		AG289–29955Ethylene-responsive transcription factor, Morus notabilis50.3130.3490.331
Ug31697F: GGAGGTGATGGAGAAGGTGAGA
R: CAACCCATCACAATCTCACATCA		GA302–31259-40.4580.5270.459
Mean 4.4710.4180.5870.531
Ug, Cl, Ct, PS, RM, AS, Ta and PF represent unigene, cluster, contig, primer sequence, repeat motif, allele size, annealing temperature and putative function respectively; “-” represents no blast hits with known proteins deposited in the public databases.

Share and Cite

MDPI and ACS Style

Dang, Z.; Huang, L.; Jia, Y.; Lockhart, P.J.; Fong, Y.; Tian, Y. Identification of Genic SSRs Provide a Perspective for Studying Environmental Adaptation in the Endemic Shrub Tetraena mongolica. Genes 2020, 11, 322. https://0-doi-org.brum.beds.ac.uk/10.3390/genes11030322

AMA Style

Dang Z, Huang L, Jia Y, Lockhart PJ, Fong Y, Tian Y. Identification of Genic SSRs Provide a Perspective for Studying Environmental Adaptation in the Endemic Shrub Tetraena mongolica. Genes. 2020; 11(3):322. https://0-doi-org.brum.beds.ac.uk/10.3390/genes11030322

Chicago/Turabian Style

Dang, Zhenhua, Lei Huang, Yuanyuan Jia, Peter J. Lockhart, Yang Fong, and Yunyun Tian. 2020. "Identification of Genic SSRs Provide a Perspective for Studying Environmental Adaptation in the Endemic Shrub Tetraena mongolica" Genes 11, no. 3: 322. https://0-doi-org.brum.beds.ac.uk/10.3390/genes11030322

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop