Simple sequence repeat DNA markers in alfalfa and perennial and annual Medicago species

Eco-Geographical, Morphological and Molecular Characterization of a Collection of the Perennial Endemic Species Medicago tunetana (Murb.) A.W. Hill (Fabaceae) from Tunisia

In order to characterize and conserve the endemic pastoral species Medicago tunetana, many prospecting missions were carried out in mountainous regions of the Tunisian ridge. Twenty-seven eco-geographical and morphological traits were studied for six M. tunetana accessions and followed by molecular analysis using seven Simple Sequence Repeat (SSR). Only five markers were polymorphic and reproductible in the six M. tunetana populations. A total of 54 alleles were observed with an average of 10.8 bands/primer/genotype. Mean Polymorphism Information Content (PIC), Nei gene diversity (h) Shannon’s information index (I) indicated the high level of polymorphism. The generated dendrogram with hierarchical UPGMA cluster analysis grouped accessions into two main groups with various degree of subclustring. All the studied accessions shared 57% of genetic similarity. Analysis of variance showed high significant difference between morphological traits among M. tunetana populations where MT3 from Kesra showed different morphological patterns regarding leaf, pod and seeds traits. Canonical correspondence analysis (CCA) showed two principal groups of M. tunetana populations based on potassium, total and active lime contents in soil. Our results suggest that SSR markers developed in M. truncatula could be a valuable tool to detect polymorphism in M. tunetana. Furthermore, the studied morphological markers showed a large genetic diversity among M. tunetana populations. This approach may be applicable for the analysis of intra specific variability in M. tunetana accessions. Our study could help in the implementation of an effective and integrated conservation programs of perennial endemic Medicago.

High-density linkage map construction and mapping QTL for yield and yield components in autotetraploid alfalfa using RAD-seq

BACKGROUND

Alfalfa (Medicago sativa L.) is an important forage crop grown worldwide. Alfalfa is called the "queen of forage crops" due to its high forage yield and nutritional characteristics. The aim of this study was to undertake quantitative trait loci (QTL) mapping of yield and yield-related traits in an F₁ population of two alfalfa varieties that differ in their yield and yield-related traits.

RESULTS

We constructed a high-density linkage map using single nucleotide polymorphism (SNP) markers generated by restriction-site associated DNA sequencing (RAD-seq). The linkage map contains 4346 SNP and 119 simple sequence repeat (SSR) markers, with 32 linkage groups for each parent. The average marker distances were 3.00 and 1.32 cM, with coverages of 3455 cM and 4381 cM for paternal and maternal linkage maps, respectively. Using these maps and phenotypic data, we identified a total of 21 QTL for yield and yield components, including five for yield, five for plant height, five for branch number, and six for shoot diameter. Among them, six QTL were co-located for more than one trait. Five QTL explained more than 10% of the phenotypic variation.

CONCLUSIONS

We used RAD-seq to construct a linkage map for alfalfa that greatly enhanced marker density compared to previous studies. This high-density linkage map of alfalfa is a useful reference for mapping yield-related traits. Identified yield-related loci could be used to validate their usefulness in developing markers for maker-assisted selection in breeding populations to improve yield potential in alfalfa.

Development and characterization of simple sequence repeat (SSR) markers based on RNA-sequencing of Medicago sativa and in silico mapping onto the M. truncatula genome

Sufficient codominant genetic markers are needed for various genetic investigations in alfalfa since the species is an outcrossing autotetraploid. With the newly developed next generation sequencing technology, a large amount of transcribed sequences of alfalfa have been generated and are available for identifying SSR markers by data mining. A total of 54,278 alfalfa non-redundant unigenes were assembled through the Illumina HiSeqTM 2000 sequencing technology. Based on 3,903 unigene sequences, 4,493 SSRs were identified. Tri-nucleotide repeats (56.71%) were the most abundant motif class while AG/CT (21.7%), AGG/CCT (19.8%), AAC/GTT (10.3%), ATC/ATG (8.8%), and ACC/GGT (6.3%) were the subsequent top five nucleotide repeat motifs. Eight hundred and thirty- seven EST-SSR primer pairs were successfully designed. Of these, 527 (63%) primer pairs yielded clear and scored PCR products and 372 (70.6%) exhibited polymorphisms. High transferability was observed for ssp falcata at 99.2% (523) and 71.7% (378) in M. truncatula. In addition, 313 of 527 SSR marker sequences were in silico mapped onto the eight M. truncatula chromosomes. Thirty-six polymorphic SSR primer pairs were used in the genetic relatedness analysis of 30 Chinese alfalfa cultivated accessions generating a total of 199 scored alleles. The mean observed heterozygosity and polymorphic information content were 0.767 and 0.635, respectively. The codominant markers not only enriched the current resources of molecular markers in alfalfa, but also would facilitate targeted investigations in marker-trait association, QTL mapping, and genetic diversity analysis in alfalfa.

Complex patterns of autopolyploid evolution in alfalfa and allies (Medicago sativa; Leguminosae)

PREMISE OF THE STUDY

Although there is growing evidence that autopolyploidy is a widespread and important evolutionary phenomenon, it has received less attention than allopolyploidy. Medicago sativa comprises several diploid and autopolyploid taxa, including autotetraploid cultivated alfalfa, and affords an opportunity to elucidate the evolutionary history of a morphologically and genetically complex autopolyploid system.

METHODS

Phylogenies and haplotype networks were constructed from two chloroplast noncoding regions (rpl20-rps12 and trnS-trnG spacers) across seven diploid and polyploid infraspecific taxa of M. sativa and five additional closely related Medicago species, and genetic differentiation was estimated.

KEY RESULTS

The two most prominent M. sativa autopolyploids have contrasting evolutionary histories. Chloroplast data support a simple autopolyploid origin of subsp. sativa (alfalfa) from diploid subsp. caerulea, from which it is distinguishable in several quantitative characters. In contrast, morphologically identical diploid and autopolyploid cytotypes of subsp. falcata were found to possess very different chloroplast haplotypes, suggesting past introgression from M. prostrata into the polyploid. Despite the presence of hybrids between tetraploid subspecies falcata and sativa, there was little evidence of introgression of chloroplast genomes from either subspecies into the other.

CONCLUSIONS

Autopolyploid evolution in M. sativa is complicated and has followed very different paths in different subspecific taxa. The potential exists for gene flow in virtually all combinations of subspecies both within and between ploidies, yet despite the existence of hybrids, morphologically and genetically distinctive subspecies persist.

Using RNA-Seq for gene identification, polymorphism detection and transcript profiling in two alfalfa genotypes with divergent cell wall composition in stems

BACKGROUND

Alfalfa, [Medicago sativa (L.) sativa], a widely-grown perennial forage has potential for development as a cellulosic ethanol feedstock. However, the genomics of alfalfa, a non-model species, is still in its infancy. The recent advent of RNA-Seq, a massively parallel sequencing method for transcriptome analysis, provides an opportunity to expand the identification of alfalfa genes and polymorphisms, and conduct in-depth transcript profiling.

RESULTS

Cell walls in stems of alfalfa genotype 708 have higher cellulose and lower lignin concentrations compared to cell walls in stems of genotype 773. Using the Illumina GA-II platform, a total of 198,861,304 expression sequence tags (ESTs, 76 bp in length) were generated from cDNA libraries derived from elongating stem (ES) and post-elongation stem (PES) internodes of 708 and 773. In addition, 341,984 ESTs were generated from ES and PES internodes of genotype 773 using the GS FLX Titanium platform. The first alfalfa (Medicago sativa) gene index (MSGI 1.0) was assembled using the Sanger ESTs available from GenBank, the GS FLX Titanium EST sequences, and the de novo assembled Illumina sequences. MSGI 1.0 contains 124,025 unique sequences including 22,729 tentative consensus sequences (TCs), 22,315 singletons and 78,981 pseudo-singletons. We identified a total of 1,294 simple sequence repeats (SSR) among the sequences in MSGI 1.0. In addition, a total of 10,826 single nucleotide polymorphisms (SNPs) were predicted between the two genotypes. Out of 55 SNPs randomly selected for experimental validation, 47 (85%) were polymorphic between the two genotypes. We also identified numerous allelic variations within each genotype. Digital gene expression analysis identified numerous candidate genes that may play a role in stem development as well as candidate genes that may contribute to the differences in cell wall composition in stems of the two genotypes.

CONCLUSIONS

Our results demonstrate that RNA-Seq can be successfully used for gene identification, polymorphism detection and transcript profiling in alfalfa, a non-model, allogamous, autotetraploid species. The alfalfa gene index assembled in this study, and the SNPs, SSRs and candidate genes identified can be used to improve alfalfa as a forage crop and cellulosic feedstock.

Using RNA-Seq for gene identification, polymorphism detection and transcript profiling in two alfalfa genotypes with divergent cell wall composition in stems

Background

Alfalfa, [Medicago sativa (L.) sativa], a widely-grown perennial forage has potential for development as a cellulosic ethanol feedstock. However, the genomics of alfalfa, a non-model species, is still in its infancy. The recent advent of RNA-Seq, a massively parallel sequencing method for transcriptome analysis, provides an opportunity to expand the identification of alfalfa genes and polymorphisms, and conduct in-depth transcript profiling.

Results

Cell walls in stems of alfalfa genotype 708 have higher cellulose and lower lignin concentrations compared to cell walls in stems of genotype 773. Using the Illumina GA-II platform, a total of 198,861,304 expression sequence tags (ESTs, 76 bp in length) were generated from cDNA libraries derived from elongating stem (ES) and post-elongation stem (PES) internodes of 708 and 773. In addition, 341,984 ESTs were generated from ES and PES internodes of genotype 773 using the GS FLX Titanium platform. The first alfalfa (Medicago sativa) gene index (MSGI 1.0) was assembled using the Sanger ESTs available from GenBank, the GS FLX Titanium EST sequences, and the de novo assembled Illumina sequences. MSGI 1.0 contains 124,025 unique sequences including 22,729 tentative consensus sequences (TCs), 22,315 singletons and 78,981 pseudo-singletons. We identified a total of 1,294 simple sequence repeats (SSR) among the sequences in MSGI 1.0. In addition, a total of 10,826 single nucleotide polymorphisms (SNPs) were predicted between the two genotypes. Out of 55 SNPs randomly selected for experimental validation, 47 (85%) were polymorphic between the two genotypes. We also identified numerous allelic variations within each genotype. Digital gene expression analysis identified numerous candidate genes that may play a role in stem development as well as candidate genes that may contribute to the differences in cell wall composition in stems of the two genotypes.

Conclusions

Our results demonstrate that RNA-Seq can be successfully used for gene identification, polymorphism detection and transcript profiling in alfalfa, a non-model, allogamous, autotetraploid species. The alfalfa gene index assembled in this study, and the SNPs, SSRs and candidate genes identified can be used to improve alfalfa as a forage crop and cellulosic feedstock.

The use of neutral and non-neutral SSRs to analyse the genetic structure of a Tunisian collection of Medicago truncatula lines and to reveal associations with eco-environmental variables

In this study, we investigated the genetic diversity of a collection of 136 Medicago truncatula lines from 10 Tunisian natural populations collected in well-defined locations and in various ecological conditions of soil, salinity and water availability. The genetic diversity was evaluated using a set of 18 microsatellites (SSRs), representing the 8 chromosomes of M. truncatula. A neutrality test showed that 7 SSRs were non-neutral with evidence of balancing selection. The 11 neutral SSRs revealed a geographical pooling with the Tunisian Dorsale axis restricting migration of alleles. The 7 non-neutral alleles demonstrate a correlation with rainfall, altitude and salinity environmental variables suggesting that these SSRs are linked to genes involved in water use efficiency, resistance to salinity or adaptation to altitude, and that there is local adaptation of M. truncatula to these variables. This demonstrates that the choice of so-called neutral markers should be carefully evaluated in population genetic studies. This study illustrates the genetic diversity occurring in natural Tunisian populations of M. truncatula and describes the first collection of this species dedicated to natural variation involved in adaptation to the environment.

Genetic diversity of some Mediterranean populations of the cultivated alfalfa (Medicago sativa L.) using SSR markers

This species study was to investigate the differentiation level among 26 populations in which 12 are locals originating from the Tunisian South and 14 introduced from Italy, Austerely, France and Morocco with two SSR markers. These highly polymorphic and co dominant markers, together with recent population genetic statistic extended to autotetraploids, offer tools to analyse genetic diversity in alfalfa. The number of alleles per locus varied between 8 and 9. The genetic similarity between these various populations is estimated by the index of Rogers and Tanimoto. Genetic diversity is analysed by two statistical procedures: Hierarchical classification and Correspondence Factorial Analysis (CFA). Four large groups were obtained.

Development and characterization of EST-SSR markers for the crown rust pathogen of ryegrass (Puccinia coronata f.sp. lolii)

The causative organism of crown rust in ryegrasses (Puccinia coronata f.sp. lolii) is an obligate biotroph that causes significant economic losses within the temperate grazing industries of dairy, meat, and wool production. This study reports on the development, transferability, and utility of gene-associated simple sequence repeat (SSR) molecular markers for crown rust. Analysis of 1,100 expressed sequence tag (EST) sequences from a urediniospore-derived cDNA library detected 55 SSR loci. The majority of EST-SSR arrays contained perfect trinucleotide repeats with consistently low repeat numbers, and the motifs (ACC)n and (CAT)n were most commonly represented. DNA extraction from single pustules, in conjunction with multiple displacement amplification, provided the basis for PCR-based screening to evaluate genetic marker performance. An example of the identification of intraspecific genetic diversity was obtained from the analysis of 16 P. coronata isolates originating from the United Kingdom, Australia, New Zealand, and Japan. A subset of 12 robust EST-SSR markers was informative for determination of pathogen diversity within and between these localities. It was also demonstrated that crown rust EST-SSR markers were capable of cross-amplification in closely related fungal taxa (Puccinia spp.) and filamentous fungi within the Ascomycota.

Genetic diversity among alfalfa (Medicago sativa) cultivars coming from a breeding program, using SSR markers

Alfalfa (Medicago sativa) is an autotetraploid, allogamous and heterozygous species whose cultivars are synthetic populations. The breeders apply selection pressure for some agronomic traits within a breeding pool to increase the frequency of favorable individuals. The objective of this study was to investigate the differentiation level among seven cultivars originating from one breeding program, and between these cultivars and the breeding pool, with eight SSR markers. These highly polymorphic and codominant markers, together with recent population genetic statistics extended to autotetraploids, offer tools to analyse genetic diversity in alfalfa. The number of alleles per locus varied between 3 and 24. All loci were at a panmictic equilibrium in the cultivars, except one, probably because of null alleles. With seven SSR loci, each cultivar was at panmictic equilibrium. The mean gene diversity was high, ranging from 0.665 to 0.717 in the cultivars. The parameter F(ST) indicated a low but significant diversity among cultivars. Among 21 pairs of cultivars, 15 were significantly different. The breeding pool also had a high diversity, and was significantly different from each cultivar except the most recent one. Considering the characteristics of the breeding program and the mode of cultivar elaboration, we found that they were unable to generate a large variety differentiation. Estimation of population genetics parameters at SSR loci can be applied for assessing the differences between cultivars or populations, either for variety distinction or the management of genetic resources.

A genetic linkage map in autotetraploid lucerne adapted to northern Australia, and use of the map to identify DNA markers linked to resistance to Phytophthora medicaginis

Phytophthora root rot, caused by Phytophthora medicaginis, is a major limitation to lucerne (Medicago sativa L.) production in Australia and North America. Quantitative trait loci (QTLs) involved in resistance to P. medicaginis were identified in a lucerne backcross population of 120 individuals. A genetic linkage map was constructed for tetraploid lucerne using 50 RAPD (randomly amplified polymorphic DNA), 104 AFLP (amplified fragment length polymorphism) markers, and one SSR (simple sequence repeat or microsatellite) marker, which originated from the resistant parent (W116); 13 markers remain unlinked. The linkage map contains 18 linkage groups covering 2136.5 cM, with an average distance of 15.0 cM between markers. Four of the linkage groups contained only either 2 or 3 markers. Using duplex markers and repulsion phase linkages the map condensed to 7 homology groups and 2 unassigned linkage groups. Three regions located on linkage groups 2, 14, and 18, were identified as associated with root reaction and the QTLs explained 6–15% of the phenotypic variation. The research also indicates that different resistance QTLs are involved in conferring resistance in different organs. Two QTLs were identified as associated with disease resistance expressed after inoculation of detached leaves. The marker, W11-2 on group 18, identified as associated with root reaction, contributed 7% of the phenotypic variation in leaf response in our population. This marker appears to be linked to a QTL encoding a resistance factor contributing to both root and leaf reaction. One other QTL, not identified as associated with root reaction, was positioned on group 1 and contributed to 6% of the variation. This genetic linkage map provides an entry point for future molecular-based improvement of lucerne in Australia, and markers linked to the QTLs we have reported should be useful for marker-assisted selection for partial resistance to P. medicaginis in lucerne.

Medicago truncatula EST-SSRs reveal cross-species genetic markers for Medicago spp

Expressed sequence tags (ESTs) are important resources for gene discovery and molecular marker development. From over 147,000 ESTs of Medicago truncatula, we have identified 4,384 ESTs containing perfect simple sequence repeats (EST-SSR) of di-, tri-, tetra- or pentanucleotides. Six hundred sixteen primer pairs (PPs) were designed and screened over a panel of eight genotypes representing six Medicago spp. and subspecies. Nearly, 74% (455) of the PPs produced characteristic SSR bands of expected size length in at least one Medicago species. Four hundred six (89%) of these 455 PPs produced SSR bands in all eight genotypes tested. Only 17 PPs were M. truncatula -specific. High levels of polymorphism (>70%) were detected for these markers in alfalfa, M. truncatula, and other annual medics. About 48% of the reported markers are part of gene transcripts linked to putative functions. Our results indicate that the SSR markers developed from M. truncatula ESTs are valuable genetic markers for the Medicago genus. These markers will be useful in establishing the genomic relationships of M. truncatula to important forage legume crops such as alfalfa and other annual medics.

High variability and disomic segregation of microsatellites in the octoploid Fragaria virginiana Mill. (Rosaceae)

The objectives of the present study were to develop microsatellite markers for the wild strawberry, Fragaria virginiana, to evaluate segregation patterns of microsatellite alleles in this octoploid species, and assess genetic variability at microsatellite loci in a wild population. A genomic library was screened for microsatellite repeats and several PCR primers were designed and tested. We also tested the use of heterologous primers and found that F. virginiana primers amplified products in cultivated strawberry, Fragaria x ananassa Duch. and Fragaria chiloensis. Similarly, microsatellite loci developed from cultivated strawberry also successfully amplified F. virginiana loci. We investigated four microsatellite loci in detail, three developed from F. virginiana and one from cultivated strawberry. A survey of 100 individuals from a population of F. virginiana in Pennsylvania demonstrated high heterozygosities (H(e) or gene diversity ranged from 0.80 to 0.88 per locus) and allelic diversity (12-17 alleles per locus), but individual plants had no more than two alleles per locus. Segregation patterns in parents and progeny of two controlled crosses at these four loci were consistent with disomic Mendelian inheritance. Together these findings suggest that the genome of F. virginiana is "highly diploidized" and at least a subset of microsatellite loci can be treated as codominant, diploid markers. Significant heterozygote deficiencies were found at three of the four loci for hermaphroditic individuals but for only one locus among females in this gynodioecious species.

Genomic organization of rDNA loci in natural populations of Medicago truncatula Gaertn

Medicago truncatula Gaertn. is an annual self-pollinating species characterized by a diploid complement 2n = 16 and low DNA content. It responds very well to transformation methods so it is used as a model species for Leguminosae. In contrast with the advanced studies in molecular biology, cytogenetic research has remained limited even though it is an extremely valuable approach to the analysis of the genome structure. In the present study we examined the chromosomal distribution of rDNA sequences in five natural populations of M. truncatula, explored the genomic diversity of this species and found markers for chromosome identification. FISH experiments revealed three distribution patterns of rDNA sequences, distinguished by one, two and three loci of 5S genes; 18S-5.8S-25S genes were always localized at a single locus. The results add information to the genome structure of M. truncatula, revealing a pattern of distribution of rDNA genes unobserved previously, which consists of 5S genes clustered at a single locus. The physical mapping of rDNA sequences is a first contribution towards the construction of a detailed molecular karyotype of M. truncatula.

AFLP-based assessment of genetic diversity among nine alfalfa germplasms using bulk DNA templates

Improving commercial utilization of perennial Medicago collections requires developing approaches that can rapidly and accurately characterize genetic diversity among large numbers of populations. This study evaluated the potential of using amplified fragment length polymorphism (AFLP) DNA markers, in combination with DNA bulking over multiple genotypes, as a strategy for high-throughput characterization of genetic distances (D) among alfalfa (Medicago sativa L.) accessions. Bulked DNA templates from 30 genotypes within each of nine well-recognized germplasms (African, Chilean, Flemish, Indian, Ladak, Medicago sativa subsp. falcata, Medicago sativa subsp. varia, Peruvian, and Turkistan) were evaluated using 34 primer combinations. A total of 3754 fragments were identified, of which 1541 were polymorphic. The number of polymorphic fragments detected per primer combination ranged from 20 to 85. Pairwise D estimates among the nine germplasms ranged from 0.52 to 1.46 with M. sativa subsp. falcata being the most genetically dissimilar. Unweighted pair-group method arithmetic average (UPGMA) analysis of the marker data produced two main clusters, (i) M. sativa subsp. sativa and M. sativa subsp. varia, and (ii) M. sativa subsp. falcata. Cluster-analysis results and D estimates among the Chilean, Peruvian, Flemish, and M. sativa subsp. varia germplasms supported the hypothesis that Peruvian was more similar to original Spanish introductions into Central and South America than Chilean. Hierarchical arrangement of the nine germplasms was supported by their respective geographic, subspecific, and intersubspecific hybrid origins. Subsets of as few as seven highly informative primer pairs were identified that produced comparable D estimates and similar heirarchical arrangements compared with the complete dataset. The results indicate that use of primer-pair subsets for AFLP analysis of bulk DNA templates could serve as a high-throughput system for accurately characterizing genetic diversity among large numbers of alfalfa populations.

Inheritance and mapping of 2n-egg production in diploid alfalfa

The production of eggs with the sporophytic chromosome number (2n eggs) in diploid alfalfa (Medicago spp.) is mainly associated with the absence of cytokinesis after restitutional meiosis. The formation of 2n eggs through diplosporic apomeiosis has also been documented in a diploid mutant of M. sativa subsp. falcata (L.) Arcang. (2n = 2x = 16), named PG-F9. Molecular tagging of 2n-egg formation appears to be an essential step towards marker-assisted breeding and map-based cloning strategies aimed at investigating and manipulating reproductive mutants of the M. sativa complex. We made controlled crosses between PG-F9 and three wild type plants of M. sativa subsp. coerulea (Less.) Schm. (2n = 2x = 16) and then hand-pollinated the F1 progenies with tetraploid plants of M. sativa subsp. sativa L. (2n = 4x = 32). As a triploid embryo block prevents the formation of 3x progenies in alfalfa because of endosperm imbalance, and owing to the negligible selfing rate, seed set in 2x-4x crosses was used to discriminate the genetic capacity for 2n-egg production. F1 plants that exhibited null or very low seed sets were classified as normal egg producers and plants with high seed sets as 2n-egg producers. A bulked segregant analysis (BSA) with RAPD (random amplified polymorphic DNA), ISSR (inter-simple sequence repeat), and AFLP (amplified fragment length polymorphism) markers was employed to identify a genetic linkage group related to the 2n-egg trait using one of the three F1 progenies. This approach enabled us to detect a paternal ISSR marker of 610 bp, generated by primer (CA)8-GC, located 9.8 cM from a putative gene (termed Tne1, two-n-eggs) that in its recessive form determines 2n eggs and a 30% recombination genomic window surrounding the target locus. Eight additional RAPD and AFLP markers, seven of maternal, and one of paternal origin, significantly co-segregated with the trait under investigation. The minimum number of quantitative trait loci (QTLs) controlling seed set in 2x-4x crosses was estimated by ANOVA and regression analysis. Four maternal and three paternal independent molecular markers significantly affected the trait. A paternal RAPD marker allele, mapped in the same linkage group of Tne1, explained 43% of the variation for seed set in 2x-4x crosses indicating the presence of a major QTL. A map of the PG-F9 chromosome regions carrying the minor genes that determine the expression level of 2n eggs was constructed using selected RAPD and AFLP markers. Two of these genes were linked to previously mapped RFLP loci belonging to groups 1 and 8. Molecular and genetic evidence support the involvement of at least five genes.