In silico mining of EST-SSRs in <Emphasis Type="Italic">Arachis hypogaea</Emphasis> L. and their utilization for genetic structure and diversity analysis in cultivars/breeding lines in Odisha,India

A total of 26,685 unutilized public domain expressed sequence tags (ESTs) of Arachis hypogaea L. were analyzed to give a total of 4442 EST-SSRs, in which 517 ESTs contained more than one simple sequence repeat (SSR). Of these EST-SSRs, 2542 were mononucleotide repeats (MNRs), 803 were dinucleotide repeats (DNRs), 1043 were trinucleotide repeats (TNRs), 40 were tetranucleotide repeats (TtNRs), six were pentanucleotide repeats (PNRs) and eight were hexanucleotide repeats (HNRs). Out of these 4442 EST-SSRs, only 1160 were found to be successful in non-redundant primer design; 1060 were simple SSRs, while the remaining 100 were compound forms. Among all the motifs, MNRs were abundant, followed by TNRs and DNRs. The AAG/CTT motif was the most abundant (~33 %) TNR, while AG/CT was the most abundant DNR. For redundancy and novelty, a stringent criterion deploying three different strategies was used and a total of 782 novel EST-SSRs were added to the public domain of peanut. These novel EST-SSR markers will be useful for qualitative and quantitative trait mapping, marker-assisted selection and genetic diversity studies in cultivated peanut as well as related Arachis species. A subset of 30 novel EST-SSRs was further randomly selected for validation and genotyping studies with eight well-known cultivars and 32 advanced breeding lines (ADBX lines, ADBY lines and ADBZ lines) from Odisha state, India. The number of polymorphic markers among accessions of A. hypogaea was low; however, a set of informative EST-SSR markers detected considerable levels of genetic variability in peanut cultivars and uncharacterized breeding lines collected from Odisha. The 30 newly developed EST-SSRs from Arachis spp. showed ~97 % amplification in Cicer arientinum and 93 % in pigeon pea. Thus, the EST-SSRs developed in this study will be a very useful asset for genetic analysis, comparative genome mapping, population genetic structure and phylogenetic inferences among wild and allied species of Arachis.     

Genetic diversity and population structure of common walnut (<Emphasis Type="Italic">Juglans regia</Emphasis>) in China based on EST-SSRs and the nuclear gene phenylalanine ammonia-lyase (<Emphasis Type="Italic">PAL</Emphasis>)

Common walnut (Juglans regia L.) is an economically important temperate tree species valued for both its nut and wood. We investigated the genetic diversity and population structure of J. regia germplasm from 13 locations in China using 10 markers derived from expressed sequences (EST-SSR) and sequence polymorphisms within the phenylalanine ammonia-lyase (PAL) gene. Analysis of the population genetic structure based on EST-SSRs showed distinct populations in northern versus southern China that were not reflected in the spatial distribution of PAL haplotypes. High levels of population differentiation were probably the result of reproductive isolation and in southern China, hybridization with Juglans sigillata. The results indicate the possible presence of distinct evolutionary lineages of J. regia in the Qinling and Daba Mountains of China and in Yunnan province that may require ecological management if they are to be retained as in situ resources.     

Transferability of the EST-SSRs developed on Nules clementine (<Emphasis Type="Italic">Citrus clementina</Emphasis> Hort ex Tan) to other <Emphasis Type="Italic">Citrus</Emphasis> species and their effectiveness for genetic mapping

Background

During the last decade, numerous microsatellite markers were developed for genotyping and to identify closely related plant genotypes. In citrus, previously developed microsatellite markers were arisen from genomic libraries and more often located in non coding DNA sequences. To optimize the use of these EST-SSRs as genetic markers in genome mapping programs and citrus systematic analysis, we have investigated their polymorphism related to the type (di or trinucleotide) or their position in the coding sequences.

Results

Among 11000 unigenes from a Clementine EST library, we have found at least one microsatellite sequence (repeated units size ranged from 2 to 6 nucleotides) in 1500 unigenes (13.6%). More than 95% of these SSRs were di or trinucleotides. If trinucleotide microsatellites were encountered trough all part of EST sequences, dinucleotide microsatellites were preferentially (50%) concentrated in the 5' 100th nucleotides. We assessed the polymorphism of 41 EST-SSR, by PCR amplification droved with flanking primers among ten Citrus species plus 3 from other genera. More than 90% of EST-SSR markers were polymorphic. Furthermore, dinucleotide microsatellite markers were more polymorphic than trinucleotide ones, probably related to their distribution that was more often located in the 5' UnTranslated Region (UTR). We obtained a good agreement of diversity relationships between the citrus species and relatives assessed with EST-SSR markers with the established taxonomy and phylogeny. To end, the heterozygosity of each genotype and all dual combinations were studied to evaluate the percentage of mappable markers. Higher values (> 45%) were observed for putative Citrus inter-specific hybrids (lime lemon, or sour orange) than for Citrus basic true species (mandarin, pummelo and citron) (<30%). Most favorable combinations for genome mapping were observed in those involving interspecific hybrid genotypes. Those gave higher levels of mappable markers (>70%) with a significant proportion suitable for synteny analysis.

Conclusion

Fourty one new EST-SSR markers were produced and were available for citrus genetic studies. Whatever the position of the SSR in the ESTs the EST-SSR markers we developed are powerful to investigate genetic diversity and genome mapping in citrus.

     

Mining and characterization of EST-SSR markers for <Emphasis Type="Italic">Zingiber officinale</Emphasis> Roscoe with transferability to other species of Zingiberaceae

Zingiber officinale is a model spice herb, well known for its medicinal value. It is primarily a vegetatively propagated commercial crop. However, considerable diversity in its morphology, fiber content and chemoprofiles has been reported. The present study explores the utility of EST-derived markers in studying genetic diversity in different accessions of Z. officinale and their cross transferability within the Zingiberaceae family. A total of 38,115 ESTs sequences were assembled to generate 7850 contigs and 10,762 singletons. SSRs were searched in the unigenes and 515 SSR-containing ESTs were identified with a frequency of 1 SSR per 25.21 kb of the genome. These ESTs were also annotated using BLAST2GO. Primers were designed for 349 EST-SSRs and 25 primer pairs were randomly picked for EST SSR study. Out of these, 16 primer pairs could be optimized for amplification in different accessions of Z. officinale as well as other species belonging to Zingiberaceae. GES454, GES466, GES480 and GES486 markers were found to exhibit 100% cross-transferability among different members of Zingiberaceae.     

A novel set of 223 EST-SSR markers in <Emphasis Type="Italic">Casuarina</Emphasis> L. ex Adans.: polymorphisms,cross-species transferability,and utility for commercial clone genotyping

Simple sequence repeat (SSR) markers are very useful for genetic applications in plants, but SSR resource for the important tree genus Casuarina L. ex Adans. is still limited. In this study, we report a novel set of 223 SSR markers in Casuarina developed from expressed sequence tag (EST) resource of GenBank. The 223 EST-SSR markers were polymorphic among 10 unrelated individuals of C. equisetifolia L. Johnson, with the number of alleles per locus (N_a), observed heterozygosity (H_o), expected heterozygosity (H_e), and polymorphic information content (PIC) averaging at 5.5, 0.72, 0.86, and 0.63, respectively. The rates of cross-species transferability ranged from 96.9% (C. glauca Sieber ex Sprengel) through 97.8% (C. cunninghamiana Miquel) to 99.1% (C. junghuhniana Miquel). Fifty-five C. equisetifolia clones widely planted in China were successfully genotyped with a subset of 20 EST-SSRs. These newly developed markers will have a great potential for genetic and breeding applications in Casuarina species and related taxa.     

Development of microsatellite markers and their use in genetic diversity and population structure analysis in <Emphasis Type="Italic">Casuarina</Emphasis>

Casuarina is a widely cultivated plantation tree species in coastal India, primarily due to its fast growth, high productivity and suitable for pulp and paper production. However, genetic studies of Casuarina have been hindered by lack of genomic resources and genetic markers. Knowledge of the genetic diversity and population structure of Casuarina germplasms will provide the basis for utilizing and improving resource in the breeding program. Keeping this in view, in the present study, we have identified a total of 11,503 simple sequence repeat (SSR) makers from 86,415 expressed sequence tags (ESTs) of Casuarina equisetifolia and C. junghuhniana after redundancy elimination. Dinucleotide repeats were the most abundant accounting for 72.5 % of all microsatellites, followed by trimer (23.4 %), hexamer (1.7 %), tetramer (1.5 %), and very few pentamer (0.6 %) repeats. Of these, 50 markers were used to estimate genetic diversity and population structure among 96 accessions of C. cunninghamiana and C. junghuhniana. EST-SSR markers revealed high level of polymorphism, detecting a total of 829 alleles with an average of 17 alleles per locus. Polymorphic information content (PIC) values ranged from 0.32 to 0.93, with an average of 0.78 per locus. The average observed (H _o ) and expected heterozygosity (H _e ) obtained was high and fairly similar in C. cunninghamiana and C. junghuhniana, thereby suggesting highly heterogeneous nature of Casuarina. Population structure using a Bayesian model-based clustering approach identified clear delineation between C. cunninghamiana and C junghuhniana. Further, these markers were also evaluated in four species of Casuarina confirming high rate of cross-species transferability. The results of this study can provide valuable insights for genetic and genomic research in Casuarina.     

Low genetic differentiation among altitudes in wild <Emphasis Type="Italic">Camellia oleifera</Emphasis>, a subtropical evergreen hexaploid plant

Camellia oleifera is a subtropical evergreen plant. Cultivated C. oleifera is the most important woody oil crop in China. Wild C. oleifera is an essential genetic resource for breeding. The patterns of genetic differentiation among altitudes/latitudes in wild C. oleifera are still unknown. Camellia oleifera may be predominantly hexaploid. The characteristics of polyploidy may lead to considerable biases in estimates of genetic diversity and differentiation. Our study used C. oleifera as a case study for analysing genetic diversity, structure and differentiation in polyploid plants using simple sequence repeats (SSRs). Wild C. oleifera samples were collected at different altitudes on the Jinggang and Lu mountains of China. The ploidy levels were determined with flow cytometry analysis. Eight highly polymorphic SSRs were used to genotype the samples. Genetic diversity and structure were analysed. Various estimates of genetic differentiation were compared. The flow cytometry results indicated that wild C. oleifera samples were all hexaploid at various altitudes of the Jinggang and Lu mountains. High levels of genetic diversity were found on both the Jinggang and Lu mountains. Genetic structure analyses indicated clear genetic differentiation between the Jinggang and Lu mountains and lower genetic differentiation among altitudes within each mountain. Classical genetic differentiation estimates of F_st failed to discriminate genetic differentiation between and within mountains. The Rho statistic showed a moderate level of genetic differentiation between mountains and lower levels of genetic differentiation within each mountain. Our study demonstrates that Rho is the statistic of choice for estimating genetic differentiation in polyploids.     

Genetic diversity of cucumber estimated by morpho-physiological and EST-SSR markers

In the present study, genetic variation among 40 cucumber genotypes was analyzed by means of morpho-physiological traits and 21 EST-SSR markers. Diversity was observed for morpho-physiological characters like days to 50% female flowering (37–46.9, number of fruits/plant (1.33–5.80), average fruit weight (41–333), vine length (36–364), relative water content (58.5–92.7), electrolyte leakage (15.9–37.1), photosynthetic efficiency (0.40–0.75) and chlorophyll concentration index (11.1–28.6). The pair wise Jaccard similarity coefficient ranged from 0.00 to 0.27 for quantitative traits and 0.24 to 0.96 for EST-SSR markers indicating that the accessions represent genetically diverse populations. With twenty-one EST-SSR markers, polymorphism revealed among 40 cucumber genotypes, number of alleles varied 2–6 with an average 3.05. Polymorphism information content varied from 0.002 to 0.989 (mean = 0.308). The number of effective allele (Ne), expected heterozygosity (He) and unbiased expected heterozygosity (uHe) of these EST-SSRs were 1.079–1.753, 0.074–0.428 and 0.074–0.434, respectively. Same 21 EST-SSR markers transferability checked in four other Cucumis species: snapmelon (Cucumis melo var. momordica), muskmelon (Cucumis melo L.), pickling melon (Cucumis melo var. conomon) and wild muskmelon (Cucumis melo var. agrestis) with frequency of 61.9, 95.2, 76.2, and 76.2%, respectively. Present study provides useful information on variability, which can assist geneticists with desirable traits for cucumber germplasm utilization. Observed physiological parameters may assists in selection of genotype for abiotic stress tolerance also, EST-SSR markers may be useful for genetic studies in related species.     

Genetic diversity and population structure of a drought-tolerant species of <Emphasis Type="Italic">Eucalyptus</Emphasis>, using microsatellite markers

Given the impact of climate change on the availability of water resources, it becomes necessary the use of plant species well suited to planting on dryland sites. Eucalyptus cladocalyx, a native tree of South Australia, is capable of growing under relatively dry environments and saline soils. Two hundred twenty simple sequence repeat (microsatellites) markers, from a consensus linkage map of Eucalyptus, were selected to examine genetic diversity and population structure in a collection of E. cladocalyx introduced to southern Atacama Desert, Chile. A total of 130 microsatellites were successfully amplified, some of which are associated with quantitative traits of interest in Eucalyptus. Genetic analysis revealed a total of 457 alleles, ranging from 2 to 8 alleles per locus. A moderate level of genetic diversity (He = 0.492) and differentiation (FST = 0.086) was found among the populations. Mount Remarkable and Marble Range showed the highest and lowest level of genetic diversity, respectively. The Bayesian clustering analysis revealed three homogeneous genetic groups confirming that the individuals of E. cladocalyx from natural forest are highly and significantly structured. These results provide a novel information for the development of breeding strategies in E. cladocalyx by using marker-assisted selection in regions with low rainfall patterns.     

Genome-wide identification of simple sequence repeats and development of polymorphic SSR markers for genetic studies in tea plant (<Emphasis Type="Italic">Camellia sinensis</Emphasis>)

The tea plant (Camellia sinensis (L.) O. Kuntze) is one of the most popular non-alcoholic beverage crops worldwide. The availability of complete genome sequences for the Camellia sinensis var. ‘Shuchazao’ has provided the opportunity to identify all types of simple sequence repeat (SSR) markers by genome-wide scan. In this study, a total of 667,980 SSRs were identified in the ~?3.08 Gb genome, with an overall density of 216.88 SSRs/Mb. Dinucleotide repeats were predominant among microsatellites (72.25%), followed by trinucleotide repeats (15.35%), while the remaining SSRs accounted for less than 13%. The motif AG/CT (49.96%) and AT/TA (40.14%) were the most and the second most abundant among all identified SSR motifs, respectively; meanwhile, AAT/ATT (41.29%) and AAAT/ATTT (67.47%) were the most common among trinucleotides and tetranucleotides, respectively. A total of 300 primer pairs were designed to screen six tea cultivars for polymorphisms of SSR markers using the five selected repeat types of microsatellite sequences. The resulting 96 SSR markers that yielded polymorphic and unambiguous bands were further deployed on 47 tea cultivars for genetic diversity assessment, demonstrating high polymorphism of these SSR markers. Remarkably, the dendrogram revealed that the phylogenetic relationships among these tea cultivars are highly consistent with their genetic backgrounds or places of origin. The identified genome-wide SSRs and newly developed SSR markers will provide a powerful means for genetic researches in tea plant, including genetic diversity and evolutionary origin analysis, fingerprinting, QTL mapping, and marker-assisted selection for breeding.     

Development and characterization of 24 chloroplast microsatellite markers for two species of <Emphasis Type="Italic">Eranthis</Emphasis> (Ranunculaceae)

Chloroplast microsatellites for two Korean endemic species, Eranthis byunsanensis and E. pungdoensis (Ranunculaceae), were isolated to address the questions of their distributional patterns and evolutionary relationships, using next-generation sequencing. Twenty-four polymorphic chloroplast microsatellite markers for these two species were developed, and then characterized in 65 individuals (55 individuals of E. byunsanensis and 10 individuals of E. pungdoensis). The number of alleles per locus ranged from 2 to 9; the average number of alleles across all the loci scored 4.792. The unbiased diversity per locus ranged from 0.089 to 0.880; the unbiased diversity averaged over all the loci was 0.646. The developed markers were successfully amplified for three congeneric species, E. stellata, E. pinnatifida, and E. longistipitata. The markers developed in this study can provide a valuable and important tool for understanding genetic variations, population structures, evolutionary histories and phylogeography of E. byunsanensis, E. pungdoensis, and related species.     

Microsatellite markers for the Ussuri white-toothed shrew (Soricidae: <Emphasis Type="Italic">Crocidura lasiura</Emphasis>) developed by Ion Torrent sequencing and their application to the shrew populations in disturbed forests

Understanding the impacts of forest management practices on genetic diversity is essential for effective animal management and conservation. We characterized novel microsatellite loci in the Ussuri white-toothed shrew (Crocidura lasiura Dobson, 1890) to test the impacts of anthropogenic thinning of forest trees on the shrew populations and their genetic diversity. Using Ion Torrent sequencing technology, we characterized 611 potential microsatellite markers with complete di- to tetra-nucleotide motifs, identifying nine polymorphic loci. The observed and expected heterozygosities across the nine loci were 0.526 and 0.527, respectively. Mean allelic diversity was 5.2 alleles per locus, with the mean polymorphism information content at 0.498. In comparison among shrew populations, which inhabited in the forests thinned in 2004 (CLA; n = 10), 2008 (CLB; n = 9) and 2014 (CLC; n = 3), the observed heterozygosities are similar among the three populations (0.525 at CLA, 0.532 at CLB and 0.519 at CLC), whereas the expected heterozygosities were much lower in population of CLC (0.377) than that of CLA (0.509) and CLB (0.533). The small sample size at CLC limited effective comparison and evaluation of the impact of forest thinning on genetic diversity in this shrew population. Future application of the species-specific microsatellite markers described here to a larger sample size would be valuable in estimating the ecological parameters of shrew populations associated with existing forest management practices.