Comparative in silico analysis of EST-SSRs in angiosperm and gymnosperm tree genera

Mfind: a tool for DNA barcode analysis in angiosperms and its relationship with microsatellites using a sliding window algorithm

MAIN CONCLUSION

Mfind is a tool to analyze the impact of microsatellite presence on DNA barcode specificity. We found a significant correlation between barcode entropy and microsatellite count in angiosperm. Genetic barcodes and microsatellites are some of the identification methods in taxonomy and biodiversity research. It is important to establish a relationship between microsatellite quantification and genetic information in barcodes. In order to clarify the association between the genetic information in barcodes (expressed as Shannon's Measure of Information, SMI) and microsatellites count, a total of 330,809 DNA barcodes from the BOLD database (Barcode of Life Data System) were analyzed. A parallel sliding-window algorithm was developed to compute the Shannon entropy of the barcodes, and this was compared with the quantification of microsatellites like (AT)n, (AC)n, and (AG)n. The microsatellite search method utilized an algorithm developed in the Java programming language, which systematically examined the genetic barcodes from an angiosperm database. For this purpose, a computational tool named Mfind was developed, and its search methodology is detailed. This comprehensive study revealed a broad overview of microsatellites within barcodes, unveiling an inverse correlation between the sumz of microsatellites count and barcodes information. The utilization of the Mfind tool demonstrated that the presence of microsatellites impacts the barcode information when considering entropy as a metric. This effect might be attributed to the concise length of DNA barcodes and the repetitive nature of microsatellites, resulting in a direct influence on the entropy of the barcodes.

Plastid Phylogenomic Analyses Reveal the Taxonomic Position of Peucedanum franchetii

Peucedanum franchetii is a famous folk medicinal plant in China. However, the taxonomy of the P. franchetii has not been sufficiently resolved. Due to similar morphological features between P. franchetii and Ligusticopsis members, the World Flora Online (WFO) Plant List suggested that this species transformed into the genus Ligusticopsis and merged with Ligusticopsis likiangensis. However, both species are obviously diverse in leaf shape, bracts, and bracteoles. To check the taxonomic position of P. franchetii, we newly sequenced and assembled the plastome of P. franchetii and compared it with nine other plastomes of the genus Ligusticopsis. Ten plastomes were highly conserved and similar in gene order, codon bias, RNA editing sites, IR borders, and SSRs. Nevertheless, 10 mutation hotspot regions (infA, rps8, matK, ndhF, rps15, psbA-trnH, rps2-rpoC2, psbA-trnK, ycf2-trnL, and ccsA-ndhD) were still detected. In addition, both phylogenetic analyses based on plastome data and ITS sequences robustly supported that P. franchetii was not clustered with members of Peucedanum but nested in Ligusticopsis. P. franchetii was sister to L. likiangensis in the ITS topology but clustered with L. capillacea in the plastome tree. These findings implied that P. franchetii should be transferred to genus Ligusticopsis and not merged with L. likiangensis, but as an independent species, which was further verified by morphological evidences. Therefore, transferring P. franchetii under the genus Ligusticopsis as an independent species was reasonable, and a new combination was presented.

The Plastome Sequences of Triticum sphaerococcum (ABD) and Triticum turgidum subsp. durum (AB) Exhibit Evolutionary Changes, Structural Characterization, Comparative Analysis, Phylogenomics and Time Divergence

The mechanism and course of Triticum plastome evolution is currently unknown; thus, it remains unclear how Triticum plastomes evolved during recent polyploidization. Here, we report the complete plastomes of two polyploid wheat species, Triticum sphaerococcum (AABBDD) and Triticum turgidum subsp. durum (AABB), and compare them with 19 available and complete Triticum plastomes to create the first map of genomic structural variation. Both T. sphaerococcum and T. turgidum subsp. durum plastomes were found to have a quadripartite structure, with plastome lengths of 134,531 bp and 134,015 bp, respectively. Furthermore, diploid (AA), tetraploid (AB, AG) and hexaploid (ABD, AGA^m) Triticum species plastomes displayed a conserved gene content and commonly harbored an identical set of annotated unique genes. Overall, there was a positive correlation between the number of repeats and plastome size. In all plastomes, the number of tandem repeats was higher than the number of palindromic and forward repeats. We constructed a Triticum phylogeny based on the complete plastomes and 42 shared genes from 71 plastomes. We estimated the divergence of Hordeum vulgare from wheat around 11.04-11.9 million years ago (mya) using a well-resolved plastome tree. Similarly, Sitopsis species diverged 2.8-2.9 mya before Triticum urartu (AA) and Triticum monococcum (AA). Aegilops speltoides was shown to be the maternal donor of polyploid wheat genomes and diverged ~0.2-0.9 mya. The phylogeny and divergence time estimates presented here can act as a reference framework for future studies of Triticum evolution.

The dynamic history of gymnosperm plastomes: Insights from structural characterization, comparative analysis, phylogenomics, and time divergence

Gymnosperms are among the most endangered groups of plant species; they include ginkgo, pines (Conifers I), cupressophytes (Conifers II), cycads, and gnetophytes. The relationships among the five extant gymnosperm groups remain equivocal. We analyzed 167 available gymnosperm plastomes and investigated their diversity and phylogeny. We found that plastome size, structure, and gene order were highly variable in the five gymnosperm groups, of which Parasitaxus usta (Vieill.) de Laub. and Macrozamia mountperriensis F.M.Bailey had the smallest and largest plastomes, respectively. The inverted repeats (IRs) of the five groups were shown to have evolved through distinctive evolutionary scenarios. The IRs have been lost in all conifers but retained in cycads and gnetophytes. A positive association between simple sequence repeat (SSR) abundance and plastome size was observed, and the SSRs with the most variation were found in Pinaceae. Furthermore, the number of repeats was negatively correlated with IR length; thus, the highest number of repeats was detected in Conifers I and II, in which the IRs had been lost. We constructed a phylogeny based on 29 shared genes from 167 plastomes. With the plastome tree and 13 calibrations, we estimated the tree height between present-day angiosperms and gymnosperms to be ∼380 million years ago (mya). The placement of Gnetales in the tree agreed with the Gnetales-other gymnosperms hypothesis. The divergence between Ginkgo and cycads was estimated as ∼284 mya; the crown age of the cycads was 251 mya. Our time-calibrated plastid-based phylogenomic tree provides a framework for comparative studies of gymnosperm evolution.

De novo transcriptome analysis of white teak (Gmelina arborea Roxb) wood reveals critical genes involved in xylem development and secondary metabolism

BACKGROUND

Gmelina arborea Roxb is a fast-growing tree species of commercial importance for tropical countries due to multiple industrial uses of its wood. Wood is primarily composed of thick secondary cell walls of xylem cells which imparts the strength to the wood. Identification of the genes involved in the secondary cell wall biosynthesis as well as their cognate regulators is crucial to understand how the production of wood occurs and serves as a starting point for developing breeding strategies to produce varieties with improved wood quality, better paper pulping or new potential uses such as biofuel production. In order to gain knowledge on the molecular mechanisms and gene regulation related with wood development in white teak, a de novo sequencing and transcriptome assembly approach was used employing secondary cell wall synthesizing cells from young white teak trees.

RESULTS

For generation of transcriptome, RNA-seq reads were assembled into 110,992 transcripts and 49,364 genes were functionally annotated using plant databases; 5071 GO terms and 25,460 SSR markers were identified within xylem transcripts and 10,256 unigenes were assigned to KEGG database in 130 pathways. Among transcription factor families, C2H2, C3H, bLHLH and MYB were the most represented in xylem. Differential gene expression analysis using leaves as a reference was carried out and a total of 20,954 differentially expressed genes were identified including monolignol biosynthetic pathway genes. The differential expression of selected genes (4CL, COMT, CCoAOMT, CCR and NST1) was validated using qPCR.

CONCLUSIONS

We report the very first de novo transcriptome of xylem-related genes in this tropical timber species of commercial importance and constitutes a valuable extension of the publicly available transcriptomic resource aimed at fostering both basic and breeding studies.

Characterization and Application of EST-SSR Markers Developed From the Transcriptome of Amentotaxus argotaenia (Taxaceae), a Relict Vulnerable Conifer

Amentotaxus argotaenia (Taxaceae) is a vulnerable coniferous species with preference for shade and moist environment. Accurate estimation of genetic variation is crucial for its conservation, especially in the context of global warming. In this study, we acquired a transcriptome from A. argotaenia leaves using Illumina sequencing and de novo assembled 62,896 unigenes, of which 5510 EST-SSRs were detected. Twenty-two polymorphic EST-SSRs were successfully developed and further used to investigate genetic variation, linkage disequilibrium, and bottleneck signatures of A. argotaenia. The results showed that A. argotaenia had moderate genetic variation and high genetic differentiation, which may provide raw material to protect against climatic changes and accelerate local adaptation, respectively. No bottlenecks were found to occur in A. argotaenia. Our study not only showed that these EST markers are very effective in population genetic analysis but also lay a solid foundation for further investigating adaptive evolution and conservation strategies of A. argotaenia.

Characterization of the Rosa roxburghii Tratt transcriptome and analysis of MYB genes

Rosa roxburghii Tratt (Rosaceae) has a fruit that is flavorful, economically valuable, and highly nutritious, providing numerous health benefits. Myeloblastosis (MYB) proteins play key roles in the development and fruit quality of R. roxburghii. However, there is little available genomic and transcriptomic information for R. roxburghii. Here, a normalized cDNA library was constructed from five tissues, including the stem, leaf, flower, young fruit, and mature fruit, using the Illumina HiSeq 3000 platform. De novo assembly was performed, and 470.66 million clean reads were obtained. In total, 63,727 unigenes, with an average GC content of 42.08%, were discovered, 60,406 of which were annotated. In addition, 9,354 unigenes were assigned to Gene Ontology categories, and 20,202 unigenes were assigned to 25 Eukaryotic Ortholog Groups. Additionally, 19,508 unigenes were classified into 140 pathways of the Kyoto Encyclopedia of Genes and Genomes database. Based on the transcriptome, 163 unigenes associated with MYBs were detected. Among these genes, 75 genes were significantly expressed in the various tissues, including 10 R1 MYB, 42 R2R3 MYB, one R1R2R3 MYB, three R4 MYB and 19 atypical MYB-like proteins. The expression levels of the 12 MYB genes randomly selected for quantitative real-time PCR analysis corroborated the RNA sequencing results. A total of 37,545 microsatellites were detected, with an average expressed sequence tag–simple sequence repeat frequency of 0.59 (37,545/63,727). This transcriptome data improves our understanding of the role of MYB in R. roxburghii and will be valuable for identifying genes of interest.

Generation and characterization of expressed sequence tags (ESTs) from coralloid root cDNA library of Cycas debaoensis

A normalized full-length cDNA library was constructed from the coralloid roots of Cycas debaoensis by the DSN (duplex-specific nuclease) normalization method combined with the SMART (Switching Mechanism At 5' end of the RNA Transcript) technique. The titer of the original cDNA library was about 1.5 × 106 cfu·mL-1 and the average insertion size was about 1 kb with a high recombination rate (97%). The 5011 high-quality expressed sequence tags (ESTs) were obtained from 5393 randomly picked cDNA clones. Clustering and assembly of ESTs resulted in 2984 unique sequences, consisting of 618 contigs and 2366 singlets. EST sequence annotation revealed that 2333 and 1901 unigenes were functionally annotated in the NCBI non-redundant database and Swiss-Prot protein database, respectively. Functional analysis demonstrated that 1495 (50.1%) unigenes were associated with 4082 Gene Ontology (GO) terms. A total of 847 unigenes were grouped into 22 Cluster of Orthologous Groups (COG) functional categories. Based on the EST dataset, 22 ESTs that encoded putative receptor-like protein kinase (RLK) genes were screened. Furthermore, a total of 94 simple sequence repeats (SSRs) were discovered, of which 20 loci were successfully amplified in C. debaoensis. This study is the first EST analysis for the coralloid roots of C. debaoensis and provides a valuable genomic resource for novel gene discovery, gene expression and comparative genomics, conservation and management studies as well as applications in C. debaoensis and related cycad species.

Development of EST-SSR markers and their application in an analysis of the genetic diversity of the endangered species Magnolia sinostellata

Magnolia sinostellata is an endemic species of Magnoliaceae that is narrowly distributed in the south of Zhejiang Province, China. To explore the genetic diversity and population structure of this endangered species, this study developed sequence tag-simple sequence repeat (EST-SSR) markers based on transcriptome data of M. sinostellata. In total, 25472 SSRs were identified among 110644 unique assembled sequences with a total of 90.83 Mb and an average frequency of 23.02%. The mononucleotide (33.53%) and dinucleotide (42.08%) motifs appeared to be the most abundant. In total, 150 potential loci were randomly selected to validate the quality of the developed SSR markers; an effective PCR rate of 32.00% and a polymorphism rate of 15.33% were obtained for these loci. After performing sequencing and cloning for validation, 23 pairs of SSR primers were retained and used to characterize the genetic diversity and population structure of M. sinostellata. Overall, 204 alleles were amplified. The results of Shannon's information index (I), heterozygosity (Ho), heterozygosity (He) and Nei's expected heterozygosity (H) indicated rich genetic diversity in M. sinostellata. However, the high inbreeding coefficient and differential coefficient suggest that serious genetic drift occurred within populations, and genetic differentiation is apparent among the populations. Consequently, although M. sinostellata has high genetic diversity among populations, it is still in a serious and dangerous condition. Habitat destruction caused by human activities is the main threat to this species, and enhancing the species abundance by adopting some conservation measures should be favourable for saving the species.

De novo transcriptome assembly, development of EST-SSR markers and population genetic analyses for the desert biomass willow, Salix psammophila

Salix psammophila, a sandy shrub known as desert willow, is regarded as a potential biomass feedstock and plays an important role in maintaining local ecosystems. However, a lack of genomic data and efficient molecular markers limit the study of its population evolution and genetic breeding. In this study, chromosome counts, flow cytometry and SSR analyses indicated that S. psammophila is tetraploid. A total of 6,346 EST-SSRs were detected based on 71,458 de novo assembled unigenes from transcriptome data. Twenty-seven EST-SSR markers were developed to evaluate the genetic diversity and population structure of S. psammophila from eight natural populations in Northern China. High levels of genetic diversity (mean 10.63 alleles per locus; mean HE 0.689) were dectected in S. psammophila. The weak population structure and little genetic differentiation (pairwise FST = 0.006-0.016) were found among Population 1-Population 7 (Pop1-Pop7; Inner Mongolia and Shaanxi), but Pop8 (Ningxia) was clearly separated from Pop1-Pop7 and moderate differentiation (pairwise FST = 0.045-0.055) was detected between them, which may be influenced by local habitat conditions. Molecular variance analyses indicated that most of the genetic variation (94.27%) existed within populations. These results provide valuable genetic informations for natural resource conservation and breeding programme optimisation of S. psammophila.

Intraspecific variations in cardamom (Elettaria cardamomum Maton): assessment of genomic diversity by flow cytometry, cytological studies and ISSR analysis

BACKGROUND

The main goal of the work was to analyse intraspecific variation in Elettaria cardamomum Maton (cardamom) using genome size, cytological studies and molecular marker data. Nuclear DNA content and molecular marker details furnish data on genome size and genetic diversity respectively among the studied accessions and both complement each other for evolutionary and taxonomic studies.

RESULTS

The relative 2C genome size and total number of base pairs of cardamom was determined through flow cytometric analysis using propidium iodide staining. The nuclear DNA content was estimated in various sections of the species representing individuals from wild and cultivar genotypes following Zea mays L. CE-777 (2C = 5.43 pg) as internal reference standard. Chromosome number from growing root tip was examined following standard protocols. Twenty-six ISSR primers that generated polymorphic bands were used for genetic diversity analysis of the thirty accessions of cardamom. Estimated nuclear 2C DNA content ranged from 2.57 to 3.22 pg demonstrating 1.25-fold variation. The mean amount of 2C nuclear DNA of the cardamom was calculated as 2.87 pg which is equivalent of 2806 Mbp as the diploid genome size. The chromosome number was found to be 2n = 48. Among the thirty accessions of cardamom studied using ISSR markers, C53 (feral from Bonacaud) showed a very prominent level of genetic diversity and was lowest for C96 (Avinash-I, a released variety from Indian Institute of Spices Research, Kozhikode).

CONCLUSION

These analyses revealed the existence of genetic variability within the studied cardamom accessions. The plant specimens also differed significantly in their genome size. However, the genetic variability parameters did not show any correlation with genome size.

Cross-Species, Amplifiable EST-SSR Markers for Amentotaxus Species Obtained by Next-Generation Sequencing

Amentotaxus, a genus of Taxaceae, is an ancient lineage with six relic and endangered species. Four Amentotaxus species, namely A. argotaenia, A. formosana, A. yunnanensis, and A. poilanei, are considered a species complex because of their morphological similarities. Small populations of these species are allopatrically distributed in Asian forests. However, only a few codominant markers have been developed and applied to study population genetic structure of these endangered species. In this study, we developed and characterized polymorphic expressed sequence tag-simple sequence repeats (EST-SSRs) from the transcriptome of A. formosana. We identified 4955 putative EST-SSRs from 68,281 unigenes as potential molecular markers. Twenty-six EST-SSRs were selected for estimating polymorphism and transferability among Amentotaxus species, of which 23 EST-SSRs were polymorphic within Amentotaxus species. Among these, the number of alleles ranged from 1-4, the polymorphism information content ranged from 0.000-0.692, and the observed and expected heterozygosity were 0.000-1.000 and 0.080-0.740, respectively. Population genetic structure analyses confirmed that A. argotaenia and A. formosana were separate species and A. yunnanensis and A. poilanei were the same species. These novel EST-SSRs can facilitate further population genetic structure research of Amentotaxus species.

Comparative in silico analysis of SSRs in coding regions of high confidence predicted genes in Norway spruce (Picea abies) and Loblolly pine (Pinus taeda)

Background

Microsatellites or simple sequence repeats (SSRs) are DNA sequences consisting of 1–6 bp tandem repeat motifs present in the genome. SSRs are considered to be one of the most powerful tools in genetic studies. We carried out a comparative study of perfect SSR loci belonging to class I (≥20) and class II (≥12 and <20 bp) types located in coding regions of high confidence genes in Picea abies and Pinus taeda. SSRLocator was used to retrieve SSRs from the full length CDS of predicted genes in both species.

Results

Trimers were the most abundant motifs in class I followed by hexamers in Picea abies, while trimers and hexamers were equally abundant in Pinus taeda class I SSRs. Hexamers were most frequent within class II SSRs followed by trimers, in both species. Although the frequency of genes containing SSRs was slightly higher in Pinus taeda, SSR counts per Mbp for class I was similar in both species (P-value = 0.22); while for class II SSRs, it was significantly higher in Picea abies (P-value = 0.00009). AT-rich motifs were higher in abundance than the GC-rich motifs, within class II SSRs in both the species (P-values = 10⁻⁹ and 0). With reference to class I SSRs, AT-rich and GC-rich motifs were detected with equal frequency in Pinus taeda (P-value = 0.24); while in Picea abies, GC-rich motifs were detected with higher frequency than the AT-rich motifs (P-value = 0.0005).

Conclusions

Our study gives a comparative overview of the genome SSRs composition based on high confidence genes in the two recently sequenced and economically important conifers and, also provides information on functional molecular markers that can be applied in genetic studies in Pinus and Picea species.

Electronic supplementary material

The online version of this article (doi:10.1186/s12863-015-0304-y) contains supplementary material, which is available to authorized users.

Developing expressed sequence tag libraries and the discovery of simple sequence repeat markers for two species of raspberry (Rubus L.)

BACKGROUND

Due to a relatively high level of codominant inheritance and transferability within and among taxonomic groups, simple sequence repeat (SSR) markers are important elements in comparative mapping and delineation of genomic regions associated with traits of economic importance. Expressed sequence tags (ESTs) are a source of SSRs that can be used to develop markers to facilitate plant breeding and for more basic research across genera and higher plant orders.

METHODS

Leaf and meristem tissue from 'Heritage' red raspberry (Rubus idaeus) and 'Bristol' black raspberry (R. occidentalis) were utilized for RNA extraction. After conversion to cDNA and library construction, ESTs were sequenced, quality verified, assembled and scanned for SSRs.  Primers flanking the SSRs were designed and a subset tested for amplification, polymorphism and transferability across species. ESTs containing SSRs were functionally annotated using the GenBank non-redundant (nr) database and further classified using the gene ontology database.

RESULTS

To accelerate development of EST-SSRs in the genus Rubus (Rosaceae), 1149 and 2358 cDNA sequences were generated from red raspberry and black raspberry, respectively. The cDNA sequences were screened using rigorous filtering criteria which resulted in the identification of 121 and 257 SSR loci for red and black raspberry, respectively. Primers were designed from the surrounding sequences resulting in 131 and 288 primer pairs, respectively, as some sequences contained more than one SSR locus. Sequence analysis revealed that the SSR-containing genes span a diversity of functions and share more sequence identity with strawberry genes than with other Rosaceous species.

CONCLUSION

This resource of Rubus-specific, gene-derived markers will facilitate the construction of linkage maps composed of transferable markers for studying and manipulating important traits in this economically important genus.

Mining microsatellite markers from public expressed sequence tags databases for the study of threatened plants

BACKGROUND

Simple Sequence Repeats (SSRs) are widely used in population genetic studies but their classical development is costly and time-consuming. The ever-increasing available DNA datasets generated by high-throughput techniques offer an inexpensive alternative for SSRs discovery. Expressed Sequence Tags (ESTs) have been widely used as SSR source for plants of economic relevance but their application to non-model species is still modest.

METHODS

Here, we explored the use of publicly available ESTs (GenBank at the National Center for Biotechnology Information-NCBI) for SSRs development in non-model plants, focusing on genera listed by the International Union for the Conservation of Nature (IUCN). We also search two model genera with fully annotated genomes for EST-SSRs, Arabidopsis and Oryza, and used them as controls for genome distribution analyses. Overall, we downloaded 16 031 555 sequences for 258 plant genera which were mined for SSRsand their primers with the help of QDD1. Genome distribution analyses in Oryza and Arabidopsis were done by blasting the sequences with SSR against the Oryza sativa and Arabidopsis thaliana reference genomes implemented in the Basal Local Alignment Tool (BLAST) of the NCBI website. Finally, we performed an empirical test to determine the performance of our EST-SSRs in a few individuals from four species of two eudicot genera, Trifolium and Centaurea.

RESULTS

We explored a total of 14 498 726 EST sequences from the dbEST database (NCBI) in 257 plant genera from the IUCN Red List. We identify a very large number (17 102) of ready-to-test EST-SSRs in most plant genera (193) at no cost. Overall, dinucleotide and trinucleotide repeats were the prevalent types but the abundance of the various types of repeat differed between taxonomic groups. Control genomes revealed that trinucleotide repeats were mostly located in coding regions while dinucleotide repeats were largely associated with untranslated regions. Our results from the empirical test revealed considerable amplification success and transferability between congenerics.

CONCLUSIONS

The present work represents the first large-scale study developing SSRs by utilizing publicly accessible EST databases in threatened plants. Here we provide a very large number of ready-to-test EST-SSR (17 102) for 193 genera. The cross-species transferability suggests that the number of possible target species would be large. Since trinucleotide repeats are abundant and mainly linked to exons they might be useful in evolutionary and conservation studies. Altogether, our study highly supports the use of EST databases as an extremely affordable and fast alternative for SSR developing in threatened plants.

Microsatellite resources of Eucalyptus: current status and future perspectives

Eucalyptus is the premier paper pulp, short rotation plantation species grown all over the world. Genetic improvement programs integrating molecular marker tools are in progress in many parts of the globe to increase the productivity. Whole genome sequence and expressed sequence tags (ESTs) of the eucalypts paved way for introduction of molecular genetics and breeding in this genus. Different molecular characterization approaches have been used simultaneously in eucalypts, however, microsatellites or simple sequence repeats (SSRs) with their prolific characteristics could occupy a special niche in Eucalyptus genetic improvement. Further, highly informative SSRs were used for the clonal identity, genetic fidelity and in certification of breeder's rights. Eucalyptus genetic linkage maps generated with microsatellite loci were used successfully to identify quantitative trait loci (QTLs) for various economically important traits. Progressively more numbers of microsatellites are being linked to genes associated with adaptive and functional variations, therefore making their utility broader in genetic applications. Availability of common SSR markers across the species provides an opportunity to validate the expression of QTLs across variable genetic backgrounds and accurately compare the position of QTLs in other species. Recent evidences suggest that the presence of SSRs in micro RNAs of plant species play a role in the quantitative trait expression. Similar studies in eucalypts may provide new insights into the genetic architecture of transcript-level variations and post transcriptional gene regulation. This review on eucalypts microsatellites, highlights the availability and characteristics of genomic and eSSRs and their potential in genetic analysis of natural and breeding populations and also discusses the future prospects in population genetics and marker assisted selection.