Use of amplified fragment length polymorphism markers to assess genetic diversity of Lolium species from Portugal

Genetic Diversity and Structure of Lolium Species Surveyed on Nuclear Simple Sequence Repeat and Cytoplasmic Markers

To assess the genetic diversity and population structure of Lolium species, we used 32 nuclear simple sequence repeat (SSR) markers and 7 cytoplasmic gene markers to analyze a total of 357 individuals from 162 accessions of 9 Lolium species. This survey revealed a high level of polymorphism, with an average number of alleles per locus of 23.59 and 5.29 and an average PIC-value of 0.83 and 0.54 for nuclear SSR markers and cytoplasmic gene markers, respectively. Analysis of molecular variance (AMOVA) revealed that 16.27 and 16.53% of the total variation was due to differences among species, with the remaining 56.35 and 83.47% due to differences within species and 27.39 and 0% due to differences within individuals in 32 nuclear SSR markers set and 6 chloroplast gene markers set, respectively. The 32 nuclear SSR markers detected three subpopulations among 357 individuals, whereas the 6 chloroplast gene markers revealed three subpopulations among 160 accessions in the STRUCTURE analysis. In the clustering analysis, the three inbred species clustered into a single group, whereas the outbreeding species were clearly divided, especially according to nuclear SSR markers. In addition, almost all Lolium multiflorum populations were clustered into group C4, which could be further divided into three subgroups, whereas Lolium perenne populations primarily clustered into two groups (C2 and C3), with a few lines that instead grouped with L. multiflorum (C4) or Lolium rigidum (C6). Together, these results will useful for the use of Lolium germplasm for improvement and increase the effectiveness of ryegrass breeding.

Genetic-geographic correlation revealed across a broad European ecotypic sample of perennial ryegrass (Lolium perenne) using array-based SNP genotyping

Publically available SNP array increases the marker density for genotyping of forage crop, Lolium perenne. Applied to 90 European ecotypes composed of 716 individuals identifies a significant genetic-geographic correlation. Grassland ecosystems are ubiquitous across temperate and tropical regions, totalling 37% of the terrestrial land cover of the planet, and thus represent a global resource for understanding local adaptations to environment. However, genomic resources for grass species (outside cereals) are relatively poor. The advent of next-generation DNA sequencing and high-density SNP genotyping platforms enables the development of dense marker assays for population genetics analyses and genome-wide association studies. A high-density SNP marker resource (Illumina Infinium assay) for perennial ryegrass (Lolium perenne) was created and validated in a broad ecotype collection of 716 individuals sampled from 90 sites across Europe. Genetic diversity within and between populations was assessed. A strong correlation of geographic origin to genetic structure was found using principal component analysis, with significant correlation to longitude and latitude (P < 0.001). The potential of this array as a resource for studies of germplasm diversity and identifying traits underpinning adaptive variation is highlighted.

Identification of hemiclonal reproduction in three species of Hexagrammos marine reef fishes

Natural hybrids between the boreal species Hexagrammos octogrammus and two temperate species Hexagrammos agrammus and Hexagrammos otakii were observed frequently in southern Hokkaido, Japan. Previous studies revealed that H. octogrammus is a maternal ancestor of both hybrids; the hybrids are all fertile females and they frequently breed with paternal species. Although such rampant hybridization occurs, species boundaries have been maintained in the hybrid zone. Possible explanations for the absence of introgressions, despite the frequent backcrossing, might include clonal reproduction: parthenogenesis, gynogenesis and hybridogenesis. The natural hybrids produced haploid eggs that contained only the H. octogrammus genome (maternal ancestor) with discarded paternal genome and generated F1 -hybrid type offspring by fertilization with the haploid sperm of H. agrammus or H. otakii (paternal ancestor). This reproductive mode was found in an artificial backcross hybrid between the natural hybrid and a male of the paternal ancestor. These findings indicate that the natural hybrids adopt hybridogenesis with high possibility and produce successive generations through hybridogenesis by backcrossing with the paternal ancestor. These hybrids of Hexagrammos represent the first hybridogenetic system found from marine fishes that widely inhabit the North Pacific Ocean. In contrast with other hybridogenetic systems, these Hexagrammos hybrids coexist with all three ancestral species in the hybrid zone. The coexistence mechanism is also discussed.

Ecological barriers to gene flow between riparian and forest species of Ainsliaea (Asteraceae)

Understanding the role of habitat-associated adaptation in reducing gene flow resulting in population differentiation and speciation is a major issue in evolutionary biology. We demonstrate a significant role for habitat divergence in species isolation between two naturally hybridizing riparian and nonriparian plants, Ainsliaea faurieana and A. apiculata (Asteraceae), on Yakushima Island, Japan. By analyzing the fine-scale population structure at six sympatric sites, we found that variations in leaf shape, geography, light conditions, and genotype were strongly correlated across riverbank-forest transitions. No evidence of effective gene flow was found between the two species across the majority of the transition zones, although the NewHybrid clustering analysis confirmed interspecific hybridization. However, a relatively high level of gene flow was observed across one zone with a more diffuse ecotone and intermediate flooding and light conditions, possibly generated by human disturbances. These results suggest that the barriers to gene flow between the riparian and forest species are primarily ecological. Additional common garden experiments indicated that the two species are adaptively differentiated to contrasting flooding and light environments. Overall, our study suggests that adaptations to different habitats can lead to the formation of reproductive isolating barriers and the maintenance of distinct species boundaries.

Hybridization and asymmetric introgression between Rhododendron eriocarpum and R. indicum on Yakushima Island, southwest Japan

We examined natural hybridization between two morphologically and ecologically divergent species on Yakushima Island-the light-purple flowered Rhododendron eriocarpum native to seaside habitats and the red flowered R. indicum native to riverside habitats. By investigation of morphological traits and amplified fragment length polymorphisms (AFLP) we found that hybrid individuals morphologically and genetically intermediate between the two species grow mainly in the seaside areas of the island. The degree of morphological and genetic variation was different among the seaside hybrid populations. Although most pollinator species were exclusive to one of the two flower color groups, the halictine bees of Lasioglossum were observed in both the color types. The crucial interspecific incompatibility after pollination has not been previously described. Geographic distance between the populations is likely to be an important primary factor in bringing about natural hybridization and determining the degree of introgression between R. eriocarpum and R. indicum. The fact that hybrids occur mostly in the seaside area on Yakushima Island indicates that asymmetrical introgression occurs from R. indicum to R. eriocarpum. Alternatively, strong habitat-mediated selection from recurrent floods may prevent the hybrids from colonizing riverside habitats.

Association of candidate genes with flowering time and water-soluble carbohydrate content in Lolium perenne (L.)

We describe a candidate gene approach for associating SNPs with variation in flowering time and water-soluble carbohydrate (WSC) content and other quality traits in the temperate forage grass species Lolium perenne. Three analysis methods were used, which took the significant population structure into account. First, a linear mixed model was used enabling a structured association analysis to be incorporated with the nine populations identified in the structure analysis as random variables. Second, a within-population analysis of variance was performed. Third, a tree-scanning method was used, in which haplotype trees were associated with phenotypes on the basis of inferred haplotypes. Analysis of variance within populations identified several associations between WSC, nitrogen (N), and dry matter digestibility with allelic variants within an alkaline invertase candidate gene LpcAI. These associations were only detected in material harvested in one of the two years. By contrast, consistent associations between the L. perenne homolog (LpHD1) of the rice photoperiod control gene HD1 and flowering time were identified. One SNP, in the immediate upstream region of the LpHD1 coding sequence (C-4443-A), was significant in the linear mixed model. Within-population analysis of variance and tree-scanning analysis confirmed and extended this result to the 2118 polymorphisms in some of the populations. The merits of the tree-scanning method are compared to the single SNP analysis. The potential usefulness of the 4443 SNP in marker-assisted selection is currently being evaluated in test crosses of genotypes from this work with turf-grass varieties.

Extremely high cytoplasmic diversity in natural and breeding populations of Lolium (Poaceae)

Ten chloroplast microsatellite markers were used to characterise chloroplast genetic diversity at allelic and haplotypic level in 104 accessions of Lolium perenne, other Lolium species, Festuca species and x Festulolium cultivars. Furthermore, genetic relationships between the accessions and biogeographic distribution of haplotypes were investigated using a range of Nei's population genetic diversity measures and analysis of molecular variance (AMOVA). An extremely high number (511) of haplotypes was detected in 1575 individuals. Nei's gene diversity values among L. perenne accessions ranged between 0 and 0.333. Much of the L. perenne European ecotype diversity (61%), as calculated using AMOVA, could be attributed to within-population variance and this is likely caused by, and maintained by, high levels of natural and anthropogenic seed dispersal. Plastid gene pools and maternal lineages for L. perenne could be clearly identified. Evidence was found, using AMOVA, to show a likely migration route of L. perenne from Southern regions of Europe northwards.

Optimization of bulked AFLP analysis and its application for exploring diversity of natural and cultivated populations of red clover

Landraces and wild populations of red clover (Trifolium pratense L.) may represent a significant yet poorly characterized genetic resource of temperate grasslands. A bulking strategy with amplified fragment length polymorphism (AFLP) markers was optimized to characterize 120 red clover populations in 6 different groups: Swiss wild clover populations, Mattenklee landraces, Mattenklee cultivars, field clover cultivars, Dutch wild clover populations, and Dutch landraces. Analysis of 2 bulked samples/population consisting of 20 plants each with12 AFLP primer combinations was found optimal for determining genetic diversity and relationships within and among red clover populations and groups. Swiss wild clover populations were clearly separated from all other red clover groups and variability within and among populations was shown to be particularly high in wild clover populations and Mattenklee landraces, emphasising their value as genetic resources for improvement of red clover cultivars, as well as for conservation and restoration of biodiversity. This study shows that the ancestry of red clover landraces is primarily found in introduced cultivars rather than in natural wild clover populations. In addition, the methodological considerations presented here may help improve diversity analyses using bulked samples.

Molecular diversity of local Norwegian meadow fescue (Festuca pratensis Huds.) populations and Nordic cultivars-consequences for management and utilisation

Genetic diversity and relatedness were studied in 30 Norwegian local populations of meadow fescue (Festuca pratensis Huds.) using amplified fragment length polymorphism (AFLP) markers. The populations were also compared with 13 Nordic meadow fescue cultivars in order to analyse the distribution of variation in local populations and cultivars and to elucidate relationships between local populations and cultivars. Analysis of molecular variance (AMOVA: ) analysis showed that most of the variation was present within populations and that little variation was found between local populations and cultivars. Separate AMOVA: analyses of local populations and cultivars revealed a higher level of variation within registered cultivars than within local populations. A cluster analysis based on corrected average pairwise differences between populations showed that the populations could be divided into three clusters, of which one also contained the cultivars. These results were supported by principal coordinates analysis. The results indicate that the Norwegian meadow fescue has a narrow genetic basis and that the local populations in Norway can be divided into three groups following the most probable routes of introduction of the species into Norway. The inland populations are closely related to the cultivars and have most probably been established as a result of migration from sown meadows. The western and southern populations probably originate from human activity-for example, trade-to the coastal western and northern parts of the country and to the central parts of southern Norway.

Genetic diversity in perennial ryegrass and white clover among old Dutch grasslands as compared to cultivars and nature reserves

To support conservation policies for old Dutch grasslands that are still in agricultural use, morphological variation and AFLP-based (amplified fragment length polymorphism-based) genetic diversity was studied in perennial ryegrass and white clover populations and compared with the diversity in reference varieties. In addition, AFLP variation was also studied in grasslands located in nature reserves. From principal component analysis (PCA), it appeared that date of ear emergence in perennial ryegrass and characters related to plant vigour in white clover were the main morphological characters separating the reference varieties from the old Dutch grassland populations, and some of the grassland populations from each other. In both species, intrapopulation variation was lower for the reference varieties. Lower heterogeneity within the reference varieties was also found in the AFLP analysis. All common AFLP's observed in old Dutch grasslands could also be found in the reference varieties and nature reserves. Only a small number of low-frequency alleles found in old Dutch grasslands were absent from the other two groups. However, band frequencies of markers could vary considerably between populations, which may have been caused by selection. Analysis of the AFLP data by PCA distinguished the majority of reference varieties from the old Dutch grasslands, and showed genetic differentiation only between some grasslands. Comparison of old Dutch grasslands with grasslands in nature reserves indicated that basically the same range of genetic variation is covered by the two groups. Our study indicates that the Netherlands harbour a more or less continuous population for major parts of the diversity of perennial ryegrass and white clover. It was concluded that no specific conservation measures are presently needed to maintain genetic diversity of perennial ryegrass and white clover occurring in old Dutch grasslands.

Neighbour joining trees, dominant markers and population genetic structure

Population genetic theory for 'traditional' codominant loci showing low levels of allelic diversity (eg allozymes) has been well characterised and evaluated. In contrast, appropriate methods for the analysis of data from more recently developed marker systems are still being refined. For multilocus dominant markers such as amplified fragment length polymorphisms (AFLPs) and randomly amplified polymorphic DNA (RAPDs), the methods of data analysis can be split into two main categories. In population-based approaches, population allele frequencies are compared to obtain some measure of the partitioning of genetic diversity into within- and between-population components. In contrast, individual-based approaches use individual multilocus genotypes as the unit of analysis. Inferences on population processes such as gene flow are based on inter-relationships among individual samples as visualised on phenetic diagrams such as neighbour joining trees. Using a simulation approach coupled with neighbour joining analyses, we show that while the underlying population genetic structure is an important determinant of tree shape in the analysis of dominant data, the number of loci examined also affects the topology. At low levels of population differentiation (eg FST=0.07), mutually exclusive clustering of individuals into their respective populations can occur when sufficiently large numbers of loci are scored (eg 250 loci, typical of many AFLP studies). In contrast, unresolved star-shaped topologies can be recovered at higher levels of population differentiation (FST= >0.15) when lower numbers of loci are employed (eg 50 loci, typical of many RAPD studies). Thus, the relationship between tree topology and the extent of genetic structuring of populations is contingent upon the number of dominant loci scored. The consequences of these findings for the biological interpretation of individual-based analysis of dominant data sets are discussed.

Genetic diversity of Eurycoma longifolia inferred from single nucleotide polymorphisms

Eurycoma longifolia Jack. is a treelet that grows in the forests of Southeast Asia and is widely used throughout the region because of its reported medicinal properties. Widespread harvesting of wild-grown trees has led to rapid thinning of natural populations, causing a potential decrease in genetic diversity among E. longifolia. Suitable genetic markers would be very useful for propagation and breeding programs to support conservation of this species, although no such markers currently exist. To meet this need, we have applied a genome complexity reduction strategy to identify a series of single nucleotide polymorphisms (SNPs) within the genomes of several E. longifolia accessions. We have found that the occurrence of these SNPs reflects the geographic origins of individual plants and can distinguish different natural populations. This work demonstrates the rapid development of molecular genetic markers in species for which little or no genomic sequence information is available. The SNP markers that we have developed in this study will also be useful for identifying genetic fingerprints that correlate with other properties of E. longifolia, such as high regenerability or the appearance of bioactive metabolites.

Molecular genecology of temperature response in Lolium perenne: 2. association of AFLP markers with ecogeography

Improved winter hardiness is an important breeding objective in the forage grass Lolium perenne. This is a complex trait with several components, including the ability to survive and grow at low temperature, to acclimate to cold, tolerate wind, snow cover and ice encasement. Marker-assisted selection has the potential to increase the efficiency of breeding for improved cold tolerance. Here we describe a genecological approach to identifying molecular markers that are associated with adaptation to low winter temperatures. AFLP was used to assess the genetic diversity in 29 wild populations of ryegrass (Lolium perenne) representing a pan-European temperature cline in terms of their geographical origin. A further 18 populations from a temperature cline in Bulgaria were also analysed. In addition, two varieties and five populations representing parents of mapping families currently in use at IGER were included in the analysis. Principal coordinate (PCoA) and cluster analyses of the molecular marker data showed that the Bulgarian altitude cline populations could be distinguished clearly from the other populations. Two regression analyses were carried out; one to identify AFLP markers that correlated in frequency with low mean January temperature of the geographical origin of the population, and another to identify AFLP markers correlating in frequency with the cold tolerance phenotype of the populations, as determined by LT50 values in freezing tests. In the first analysis six AFLP markers showed significant type II trends with mean January temperature, and in the second analysis 28 bands had a significant univariate relationship with the LT50 value of the accessions. In steps 2 and 3 of the stepwise analysis a further 4 and 5 bands, respectively, improved the fit significantly. The results of the two types of regression analysis are discussed in relation to ecogeography and cold tolerance phenotype of the populations.

Land barriers and open oceans: effects on gene diversity and population structure in Avicennia germinans L. (Avicenniaceae)

Avicennia germinans L. is a widespread mangrove species occupying the west coast of Africa and the Atlantic and Pacific coasts of the Americas from the Bahamas to Brazil and Baja California to Peru. An amplified fragment length polymorphism (AFLP) molecular analysis was carried out to assess genetic architecture within this species and to evaluate the effects of the Atlantic Ocean and the Central American Isthmus (CAI) on population and regional genetic diversity and differentiation. In total, 349 polymorphic AFLP fragments were identified among 144 individuals from 14 populations from the east Atlantic, west Atlantic and east Pacific. Levels of genetic diversity varied considerably among populations, but were generally higher in populations from the east Atlantic. Regional differentiation between the Pacific coast and Atlantic populations was greater than between east and west Atlantic populations, suggesting that the CAI has had an important influence on population genetic structure in this species. The lower level of divergence of east Atlantic from west Atlantic populations suggests some dispersal across the Atlantic Ocean, although migration rates are probably low; Nm from GST equal to 0.41 and accumulation of private and rare alleles in the east Atlantic. Population differentiation did not appear to follow an isolation by distance model and has probably resulted from complex patterns of population bottlenecks, and founder events due to landscape changes during the Pleistocene, particularly in the west Atlantic. The molecular data provide no support for the treatment of east Atlantic populations as a separate species A. africana.

Genetic variation of Trigonobalanus verticillata, a primitive species of Fagaceae, in Malaysia revealed by chloroplast sequences and AFLP markers

The genetic variation of Trigonobalanus verticillata, the most recently described genus of Fagaceae, was studied using chloroplast DNA sequences and AFLP fingerprinting. This species has a restricted distribution that is known to include seven localities in tropical lower montane forests in Malaysia and Indonesia. A total of 75 individuals were collected from Bario, Kinabalu, and Fraser's Hill in Malaysia. The sequences of rbcL, matK, and three non-coding regions (atpB-rbcL spacer, trnL intron, and trnL-trnF spacer) were determined for 19 individuals from these populations. We found a total of 30 nucleotide substitutions and four length variations, which allowed identification of three haplotypes characterizing each population. No substitutions were detected within populations, while the tandem repeats in the trnL -trnF spacer had a variable repeat number of a 20-bp motif only in Kinabalu. The differentiation of the populations inferred from the cpDNA molecular clock calibrated with paleontological data was estimated to be 8.3 MYA between Bario and Kinabalu, and 16.7 MYA between Fraser's Hill and the other populations. In AFLP analysis, four selective primer pairs yielded a total of 431 loci, of which 340 (78.9%) were polymorphic. The results showed relatively high gene diversity (H(S) = 0.153 and H(T) = 0.198) and nucleotide diversity (pi(S) = 0.0132 and pi(T) = 0.0168) both within and among the populations. Although the cpDNA data suggest that little or no gene flow occurred between the populations via seeds, the fixation index estimated from AFLP data (F(ST) = 0.153 and N(ST) = 0.214) implies that some gene flow occurs between populations, possibly through pollen transfer.