A novel genome‐wide microsatellite resource for species of Eucalyptus with linkage‐to‐physical correspondence on the reference genome sequence

Keystone species in their native ranges, eucalypts, are ecologically and genetically very diverse, growing naturally along extensive latitudinal and altitudinal ranges and variable environments. Besides their ecological importance, eucalypts are also the most widely planted trees for sustainable forestry in the world. We report the development of a novel collection of 535 microsatellites for species of Eucalyptus, 494 designed from ESTs and 41 from genomic libraries. A selected subset of 223 was evaluated for individual identification, parentage testing, and ancestral information content in the two most extensively studied species, Eucalyptus grandis and Eucalyptus globulus. Microsatellites showed high transferability and overlapping allele size range, suggesting they have arisen still in their common ancestor and confirming the extensive genome conservation between these two species. A consensus linkage map with 437 microsatellites, the most comprehensive microsatellite‐only genetic map for Eucalyptus, was built by assembling segregation data from three mapping populations and anchored to the Eucalyptus genome. An overall colinearity between recombination‐based and physical positioning of 84% of the mapped microsatellites was observed, with some ordering discrepancies and sporadic locus duplications, consistent with the recently described whole genome duplication events in Eucalyptus. The linkage map covered 95.2% of the 605.8‐Mbp assembled genome sequence, placing one microsatellite every 1.55 Mbp on average, and an overall estimate of physical to recombination distance of 618 kbp/cM. The genetic parameters estimates together with linkage and physical position data for this large set of microsatellites should assist marker choice for genome‐wide population genetics and comparative mapping in Eucalyptus.     

Positioning of the major locus for <Emphasis Type="Italic">Puccinia psidii</Emphasis> rust resistance (<Emphasis Type="Italic">Ppr1</Emphasis>) on the <Emphasis Type="Italic">Eucalyptus</Emphasis> reference map and its validation across unrelated pedigrees

In this report the major locus for Puccinia psidii rust resistance, Ppr1, was positioned on the reference genetic map for Eucalyptus. Additionally, its position was validated by association genetics in a related and two unrelated pedigrees involving different Eucalyptus grandis resistant trees crossed to individuals of two other species, Eucalyptus tereticornis and Eucalyptus camaldulensis. These results are consistent with the hypothesis that Ppr1 controls a large proportion of the variation for rust resistance, strengthening its role as a major locus in Eucalyptus and providing its unequivocal genomic position on linkage group 3. A localized map with 19 microsatellite loci was built around Ppr1. Multiallelic profiles were observed at several mapped microsatellites suggesting recent tandem duplications in the genomic landscape surrounding Ppr1. Markers EMBRA125 and EMBRA1071 flank Ppr1 at 9.5% and 7% recombination, respectively, and were found to be in linkage equilibrium in a E. grandis breeding population, consistent with the expectations in outcrossed Eucalyptus. Their potential use for MAS will specifically be directed to identifying resistant offspring of P. psidii resistant parent trees that are heterozygous at Ppr1. In these circumstances, a significant amount of LD is expected to occur between specific alleles at flanking microsatellites and the resistance allele at Ppr1. Moreover, the positional information of Ppr1 paves the way for prospective undertakings in this genomic region with the upcoming availability of a draft genome for E. grandis.     

Genome-wide linkage analysis of age at onset of alcohol dependence: a comparison between microsatellites and single-nucleotide polymorphisms

Background

Using the dataset provided for Genetic Analysis Workshop 14 by the Collaborative Study on the Genetics of Alcoholism, we performed genome-wide linkage analysis of age at onset of alcoholism to compare the utility of microsatellites and single-nucleotide polymorphisms (SNPs) in genetic linkage study.

Methods

A multipoint nonparametric variance component linkage analysis method was applied to the survival distribution function obtained from semiparametric proportional hazards model of the age at onset phenotype of alcoholism. Three separate linkage analyses were carried out using 315 microsatellites, 2,467 and 9,467 SNPs, spanning the 22 autosomal chromosomes.

Results

Heritability of age at onset was estimated to be approximately 12% (p < 0.001). We observed weak correlation, both in trend and strength, of genome-wide linkage signals between microsatellites and SNPs. Results from SNPs revealed more and stronger linkage signals across the genome compared with those from microsatellites. The only suggestive evidence of linkage from microsatellites was on chromosome 1 (LOD of 1.43). Differences in map densities between the two sets of SNPs used in this study did not appear to confer an advantage in terms of strength of linkage signals.

Conclusion

Our study provided support for better performance of dense SNP maps compared with the sparse mirosatellite maps currently available for linkage analysis of quantitative traits. This better performance could be attributable to precise definition and high map resolutions achievable with dense SNP maps, thus resulting in increased power to detect possible loci affecting given trait or disease.

     

A microsatellite linkage map of Drosophila mojavensis

Background

Drosophila mojavensishas been a model system for genetic studies of ecological adaptation and speciation. However, despite its use for over half a century, no linkage map has been produced for this species or its close relatives.

Results

We have developed and mapped 90 microsatellites in D. mojavensis, and we present a detailed recombinational linkage map of 34 of these microsatellites. A slight excess of repetitive sequence was observed on the X-chromosome relative to the autosomes, and the linkage groups have a greater recombinational length than the homologous D. melanogaster chromosome arms. We also confirmed the conservation of Muller's elements in 23 sequences between D. melanogaster and D. mojavensis.

Conclusions

The microsatellite primer sequences and localizations are presented here and made available to the public. This map will facilitate future quantitative trait locus mapping studies of phenotypes involved in adaptation or reproductive isolation using this species.     

Development, characterization and mapping of microsatellite markers in Eucalyptus grandis and E. urophylla

 We report on the development, genetic characterization and linkage mapping of a battery of SSR (simple sequence repeat) loci in Eucalyptus grandis and E. urophylla. This study reveals the abundance of SSRs in Eucalyptus, the very high information content of these markers for mapping and individual identification, and demonstrates the feasibility of constructing a comprehensive microsatellite-based linkage map for Eucalyptus. Primer sequence for a set of 20 highly informative EMBRA (Eucalyptus microsatellites from Brazil) loci are made available together with their map position and estimates of the expected heterozygosity and allele size range in these two species. Using genomic library enrichment and anchored-PCR screening prior to sequencing, the efficiency of SSR marker locus development was 63% from sequencing data to operationally useful SSR loci. Absolute transportability between the two species and very high levels of allelic variability and expected heterozygosity (H) were seen at all SSR loci surveyed. The number of alleles per locus ranged from 9 to 26 with an average of 16.3±4.8. The average H of 15 loci was 0.86±0.04, 0.83±0.08 and 0.89±0.04, respectively, for E. urophylla, E. grandis and the combined two-species estimate. In the mapping analysis 16 out of 20 marker loci segregated in a fully informative configuration, allowing the determination of synteny of six homologous linkage groups between the two species. The availability of transportable, multiallelic, PCR-based co-dominant SSR loci represents a dramatic improvement in our ability to carry out detailed population genetic analysis and to search, understand, and manipulate allelic variation at QTLs (quantitative trait loci) in species of Eucalyptus. Received: 16 March 1998 / Accepted: 22 March 1998     

A set of chloroplast microsatellite primers for Eucalyptus (Myrtaceae)

We present a set of 35 chloroplast microsatellite primers for Eucalyptus. Ten of the microsatellites displayed intraspecific polymorphism, and identified nine haplotypes among 16 Eucalyptus globulus individuals. Primer conservation was high, with a polymerase chain reaction (PCR) success rate of 98% when tested on four other Eucalyptus species and seven additional myrtaceous genera. Chloroplast microsatellites have applications in phylogeographic studies, fingerprinting and progeny analysis.     

A SSR-based composite genetic linkage map for the cultivated peanut (<Emphasis Type="Italic">Arachis hypogaea</Emphasis> L.) genome

Background  

The construction of genetic linkage maps for cultivated peanut (Arachis hypogaea L.) has and continues to be an important research goal to facilitate quantitative trait locus (QTL) analysis and gene tagging for use in a marker-assisted selection in breeding. Even though a few maps have been developed, they were constructed using diploid or interspecific tetraploid populations. The most recently published intra-specific map was constructed from the cross of cultivated peanuts, in which only 135 simple sequence repeat (SSR) markers were sparsely populated in 22 linkage groups. The more detailed linkage map with sufficient markers is necessary to be feasible for QTL identification and marker-assisted selection. The objective of this study was to construct a genetic linkage map of cultivated peanut using simple sequence repeat (SSR) markers derived primarily from peanut genomic sequences, expressed sequence tags (ESTs), and by "data mining" sequences released in GenBank.     

A radiation hybrid map of the RN region in pigs demonstrates conserved gene order compared with the human and mouse genomes

We recently constructed a 7000-rad porcine whole-genome radiation hybrid (RH) panel with the primary objective of integrating linkage maps of microsatellites with evolutionary conserved genes into one ordered map. In order to evaluate the resolution of this RH panel, we have now constructed a radiation hybrid map of the Chromosome (Chr) 15q2.3-q2.6 region containing the RN gene. This gene has large effects on glycogen content in muscle and meat quality. Ten microsatellites covering a region of 55 centiMorgans and eight genes (AE3, FN1, IGFBP5, INHA, IRS1, PAX3, TNP1, and VIL1) were placed on the Sscr15 RH map. All the genes, except IRS1, were mapped on the RH map between microsatellites located in 15q2.5. The relative order of AE3 and INHA was inverted on the porcine physical map in comparison with the mouse linkage map. The order of other genes already mapped in the mouse (FN1, IGFBP5, TNP1, VIL1, INHA/AE3, and PAX3) was identical in pigs. We found no clear difference between the gene order on pig Chr 15 and human Chr 2q. Received: 4 November 1998 / Accepted: 8 February 1999     

Sex-specific population structure,natural selection,and linkage disequilibrium in a wild bird population as revealed by genome-wide microsatellite analyses

Background  

Sexual dimorphism in ecologically important traits is widespread, yet the differences in the genomic architecture between the two sexes are largely unexplored. We employed a genome-wide multilocus approach to examine the sexual differences in population subdivision, natural selection and linkage disequilibrium (LD) in a wild Siberian jay (Perisoreus infaustus) population, using genotypes at a total of 107 autosomal and Z-chromosomal microsatellites.     

Genomic Characterization of DArT Markers Based on High-Density Linkage Analysis and Physical Mapping to the Eucalyptus Genome

Diversity Arrays Technology (DArT) provides a robust, high throughput, cost-effective method to query thousands of sequence polymorphisms in a single assay. Despite the extensive use of this genotyping platform for numerous plant species, little is known regarding the sequence attributes and genome-wide distribution of DArT markers. We investigated the genomic properties of the 7,680 DArT marker probes of a Eucalyptus array, by sequencing them, constructing a high density linkage map and carrying out detailed physical mapping analyses to the Eucalyptus grandis reference genome. A consensus linkage map with 2,274 DArT markers anchored to 210 microsatellites and a framework map, with improved support for ordering, displayed extensive collinearity with the genome sequence. Only 1.4 Mbp of the 75 Mbp of still unplaced scaffold sequence was captured by 45 linkage mapped but physically unaligned markers to the 11 main Eucalyptus pseudochromosomes, providing compelling evidence for the quality and completeness of the current Eucalyptus genome assembly. A highly significant correspondence was found between the locations of DArT markers and predicted gene models, while most of the 89 DArT probes unaligned to the genome correspond to sequences likely absent in E. grandis, consistent with the pan-genomic feature of this multi-Eucalyptus species DArT array. These comprehensive linkage-to-physical mapping analyses provide novel data regarding the genomic attributes of DArT markers in plant genomes in general and for Eucalyptus in particular. DArT markers preferentially target the gene space and display a largely homogeneous distribution across the genome, thereby providing superb coverage for mapping and genome-wide applications in breeding and diversity studies. Data reported on these ubiquitous properties of DArT markers will be particularly valuable to researchers working on less-studied crop species who already count on DArT genotyping arrays but for which no reference genome is yet available to allow such detailed characterization.     

High-density genetic linkage maps with over 2,400 sequence-anchored DArT markers for genetic dissection in an F2 pseudo-backcross of Eucalyptus grandis × E. urophylla

Traits that differentiate cross-fertile plant species can be dissected by genetic linkage analysis in interspecific hybrids. Such studies have been greatly facilitated in Eucalyptus tree species by the recent development of Diversity Arrays Technology (DArT) markers. DArT is an affordable, high-throughput marker technology for the construction of high-density genetic linkage maps. Eucalyptus grandis and Eucalyptus urophylla are commonly used to produce fast-growing, disease tolerant hybrids for clonal eucalypt plantations in tropical and subtropical regions. We analysed 7,680 DArT markers in an F2 pseudo-backcross mapping pedigree based on an F1 hybrid clone of E. grandis and E. urophylla. A total of 2,440 markers (31.7%) were polymorphic and could be placed in linkage maps of the F1 hybrid and two pure-species backcross parents. An integrated genetic linkage map was constructed for the pedigree resulting in 11 linkage groups (n = 11) with 2,290 high-confidence (LOD ≥ 3.0) markers and a total map length of 1,107.6 cM. DNA sequence analysis of the mapped DArT marker fragments revealed that 43% were located in protein coding regions and 90% could be placed in the recently completed draft genome assembly of E. grandis. Together with the anchored genomic sequence information, this linkage map will allow detailed genetic dissection of quantitative traits and hybrid fitness characters segregating in the F2 progeny and will facilitate the development of markers for molecular breeding in Eucalyptus.     

AFLP genetic maps of Eucalyptus globulus and E. tereticornis

 Amplified fragment length polymorphism (AFLP) analysis is a rapid and efficient technique for detecting large numbers of DNA markers in eucalypts. We have used AFLP markers in a two-way pseudo-testcross strategy to generate genetic maps of two clones of different Eucalyptus species (E. tereticornis and E. globulus). Of 606 polymorphic fragments scored, 487 segregated in a 1 : 1 ratio, corresponding to DNA polymorphisms heterozygous in one parent and null in the other. In the maternal E. tereticornis map, 268 markers were ordered in 14 linkage groups (919 cM); the paternal E. globulus map had 200 markers in 16 linkage groups (967 cM). Results from PGRI software were compared with MAPMAKER. The average density of markers was approximately 1 per 3.9 cM. Framework markers were ordered with an average confidence level of 90%, covering 80–100% of the estimated Eucalyptus genome size. In order to investigate the homologies between the E. tereticornis and the E. globulus genetic linkage maps, we included 19 markers segregating 3 : 1 in the analysis. Some homeologous linkage groups were recognized. The linkage data developed in these maps will be used to detect loci controlling commercially important traits. Received: 17 July 1997 / Accepted: 13 October 1997     

A new genomic resource dedicated to wood formation in Eucalyptus

Background  

Renowned for their fast growth, valuable wood properties and wide adaptability, Eucalyptus species are amongst the most planted hardwoods in the world, yet they are still at the early stages of domestication because conventional breeding is slow and costly. Thus, there is huge potential for marker-assisted breeding programs to improve traits such as wood properties. To this end, the sequencing, analysis and annotation of a large collection of expressed sequences tags (ESTs) from genes involved in wood formation in Eucalyptus would provide a valuable resource.     

Construction of a high density SNP linkage map of kelp (Saccharina japonica) by sequencing Taq I site associated DNA and mapping of a sex determining locus

Background

Kelp (Saccharina japonica) has been intensively cultured in China for almost a century. Its genetic improvement is comparable with that of rice. However, the development of its molecular tools is extremely limited, thus its genes, genetics and genomics. Kelp performs an alternative life cycle during which sporophyte generation alternates with gametophyte generation. The gametophytes of kelp can be cloned and crossed. Due to these characteristics, kelp may serve as a reference for the biological and genetic studies of Volvox, mosses and ferns.

Results

We constructed a high density single nucleotide polymorphism (SNP) linkage map for kelp by restriction site associated DNA (RAD) sequencing. In total, 4,994 SNP-containing physical (tag-defined) RAD loci were mapped on 31 linkage groups. The map expanded a total genetic distance of 1,782.75 cM, covering 98.66% of the expected (1,806.94 cM). The length of RAD tags (85 bp) was extended to 400–500 bp with Miseq method, offering us an easiness of developing SNP chips and shifting SNP genotyping to a high throughput track. The number of linkage groups was in accordance with the documented with cytological methods. In addition, we identified a set of microsatellites (99 in total) from the extended RAD tags. A gametophyte sex determining locus was mapped on linkage group 2 in a window about 9.0 cM in width, which was 2.66 cM up to marker_40567 and 6.42 cM down to marker_23595.

Conclusions

A high density SNP linkage map was constructed for kelp, an intensively cultured brown alga in China. The RAD tags were also extended so that a SNP chip could be developed. In addition, a set of microsatellites were identified among mapped loci, and a gametophyte sex determining locus was mapped. This map will facilitate the genetic studies of kelp including for example the evaluation of germplasm and the decipherment of the genetic bases of economic traits.

Electronic supplementary material

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

Gene-based SSR markers for common bean (Phaseolus vulgaris L.) derived from root and leaf tissue ESTs: an integration of the BMc series

Background  

Sequencing of cDNA libraries for the development of expressed sequence tags (ESTs) as well as for the discovery of simple sequence repeats (SSRs) has been a common method of developing microsatellites or SSR-based markers. In this research, our objective was to further sequence and develop common bean microsatellites from leaf and root cDNA libraries derived from the Andean gene pool accession G19833 and the Mesoamerican gene pool accession DOR364, mapping parents of a commonly used reference map. The root libraries were made from high and low phosphorus treated plants.     

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.     

A genetic linkage map for Eucalyptus globulus with candidate loci for wood,fibre, and floral traits

A genetic linkage map containing potential candidate loci for wood, fibre and floral traits has been constructed for Eucalyptus globulus (Labill.) based on the segregation of 249 codominant loci in an outbred F₁ population of 148 individuals. The map contains 204 RFLP loci, including 31 cambium-specific expressed sequence tags (ESTs) and 14 known function genes, and 40 microsatellite and five isozyme loci. Independent male and female maps were constructed, and the 98 loci (39%) that segregated in both parents were used to combine the parental maps into an integrated map. The 249 loci mapped to 11 major linkage groups (n=11 in eucalypts) and a 12th small linkage group containing three loci that segregated in the male parent only. Total map distance is 1375 cM with an average interval of 6 cM. Forty one of the mapped loci identify known proteins (five isozymes) or sequences with known function (14 genes and 22 ESTs). The mapped genes include enzymes involved in lignin and cell-wall polysaccharide biosynthesis, and floral-development genes. This map will be used to locate quantitative trait loci for wood, fibre, and other traits in Eucalyptus. Received: 30 August 2000 / Accepted: 23 March 2001     

A microsatellite-based analysis for the detection of selection on BTA1 and BTA20 in northern Eurasian cattle (Bos taurus) populations

Background

Microsatellites surrounding functionally important candidate genes or quantitative trait loci have received attention as proxy measures of polymorphism level at the candidate loci themselves. In cattle, selection for economically important traits is a long-term strategy and it has been reported that microsatellites are linked to these important loci.

Methods

We have investigated the variation of seven microsatellites on BTA1 (Bos taurus autosome 1) and 16 on BTA20, using bovine populations of typical production types and horn status in northern Eurasia. Genetic variability of these loci and linkage disequilibrium among these loci were compared with those of 28 microsatellites on other bovine chromosomes. Four different tests were applied to detect molecular signatures of selection.

Results

No marked difference in locus variability was found between microsatellites on BTA1, BTA20 and the other chromosomes in terms of different diversity indices. Average D'' values of pairwise syntenic markers (0.32 and 0.28 across BTA 1 and BTA20 respectively) were significantly (P < 0.05) higher than for non-syntenic markers (0.15). The Ewens-Watterson test, the Beaumont and Nichol''s modified frequentist test and the Bayesian F_ST-test indicated elevated or decreased genetic differentiation, at SOD1 and AGLA17 markers respectively, deviating significantly (P < 0.05) from neutral expectations. Furthermore, lnRV, lnRH and lnRθ'' statistics were used for the pairwise population comparison tests and were significantly less variable in one population relative to the other, providing additional evidence of selection signatures for two of the 51 loci. Moreover, the three Finnish native populations showed evidence of subpopulation divergence at SOD1 and AGLA17. Our data also indicate significant intergenic linkage disequilibrium around the candidate loci and suggest that hitchhiking selection has played a role in shaping the pattern of observed linkage disequilibrium.

Conclusion

Hitchhiking due to tight linkage with alleles at candidate genes, e.g. the POLL gene, is a possible explanation for this pattern. The potential impact of selective breeding by man on cattle populations is discussed in the context of selection effects. Our results also suggest that a practical approach to detect loci under selection is to simultaneously apply multiple neutrality tests based on different assumptions and estimations.     

An enhanced microsatellite map of diploid Fragaria

A total of 45 microsatellites (SSRs) were developed for mapping in Fragaria. They included 31 newly isolated codominant genomic SSRs from F. nubicola and a further 14 SSRs, derived from an expressed sequence tagged library (EST-SSRs) of the cultivated strawberry, F. × ananassa. These, and an additional 64 previously characterised but unmapped SSRs and EST-SSRs, were scored in the diploid Fragaria interspecific F₂ mapping population (FV×FN) derived from a cross between F. vesca 815 and F. nubicola 601. The cosegregation data of these 109 SSRs, and of 73 previously mapped molecular markers, were used to elaborate an enhanced linkage map. The map is composed of 182 molecular markers (175 microsatellites, six gene specific markers and one sequence-characterised amplified region) and spans 424 cM over seven linkage groups. The average marker spacing is 2.3 cM/marker and the map now contains just eight gaps longer than 10 cM. The transferability of the new SSR markers to the cultivated strawberry was demonstrated using eight cultivars. Because of the transferable nature of these markers, the map produced will provide a useful reference framework for the development of linkage maps of the cultivated strawberry and for the development of other key resources for Fragaria such as a physical map. In addition, the map now provides a framework upon which to place transferable markers, such as genes of known function, for comparative mapping purposes within Rosaceae.