Triplet-repeat microsatellites shared among hard and soft pines

High-quality chromosome-level genome assembly and multi-omics analysis of rosemary (Salvia rosmarinus) reveals new insights into the environmental and genome adaptation

High-quality genome of rosemary (Salvia rosmarinus) represents a valuable resource and tool for understanding genome evolution and environmental adaptation as well as its genetic improvement. However, the existing rosemary genome did not provide insights into the relationship between antioxidant components and environmental adaptability. In this study, by employing Nanopore sequencing and Hi-C technologies, a total of 1.17 Gb (97.96%) genome sequences were mapped to 12 chromosomes with 46 121 protein-coding genes and 1265 non-coding RNA genes. Comparative genome analysis reveals that rosemary had a closely genetic relationship with Salvia splendens and Salvia miltiorrhiza, and it diverged from them approximately 33.7 million years ago (MYA), and one whole-genome duplication occurred around 28.3 MYA in rosemary genome. Among all identified rosemary genes, 1918 gene families were expanded, 35 of which are involved in the biosynthesis of antioxidant components. These expanded gene families enhance the ability of rosemary adaptation to adverse environments. Multi-omics (integrated transcriptome and metabolome) analysis showed the tissue-specific distribution of antioxidant components related to environmental adaptation. During the drought, heat and salt stress treatments, 36 genes in the biosynthesis pathways of carnosic acid, rosmarinic acid and flavonoids were up-regulated, illustrating the important role of these antioxidant components in responding to abiotic stresses by adjusting ROS homeostasis. Moreover, cooperating with the photosynthesis, substance and energy metabolism, protein and ion balance, the collaborative system maintained cell stability and improved the ability of rosemary against harsh environment. This study provides a genomic data platform for gene discovery and precision breeding in rosemary. Our results also provide new insights into the adaptive evolution of rosemary and the contribution of antioxidant components in resistance to harsh environments.

An annotated genetic map of loblolly pine based on microsatellite and cDNA markers

Background

Previous loblolly pine (Pinus taeda L.) genetic linkage maps have been based on a variety of DNA polymorphisms, such as AFLPs, RAPDs, RFLPs, and ESTPs, but only a few SSRs (simple sequence repeats), also known as simple tandem repeats or microsatellites, have been mapped in P. taeda. The objective of this study was to integrate a large set of SSR markers from a variety of sources and published cDNA markers into a composite P. taeda genetic map constructed from two reference mapping pedigrees. A dense genetic map that incorporates SSR loci will benefit complete pine genome sequencing, pine population genetics studies, and pine breeding programs. Careful marker annotation using a variety of references further enhances the utility of the integrated SSR map.

Results

The updated P. taeda genetic map, with an estimated genome coverage of 1,515 cM_(Kosambi) across 12 linkage groups, incorporated 170 new SSR markers and 290 previously reported SSR, RFLP, and ESTP markers. The average marker interval was 3.1 cM. Of 233 mapped SSR loci, 84 were from cDNA-derived sequences (EST-SSRs) and 149 were from non-transcribed genomic sequences (genomic-SSRs). Of all 311 mapped cDNA-derived markers, 77% were associated with NCBI Pta UniGene clusters, 67% with RefSeq proteins, and 62% with functional Gene Ontology (GO) terms. Duplicate (i.e., redundant accessory) and paralogous markers were tentatively identified by evaluating marker sequences by their UniGene cluster IDs, clone IDs, and relative map positions. The average gene diversity, H_e, among polymorphic SSR loci, including those that were not mapped, was 0.43 for 94 EST-SSRs and 0.72 for 83 genomic-SSRs. The genetic map can be viewed and queried at http://www.conifergdb.org/pinemap.

Conclusions

Many polymorphic and genetically mapped SSR markers are now available for use in P. taeda population genetics, studies of adaptive traits, and various germplasm management applications. Annotating mapped genes with UniGene clusters and GO terms allowed assessment of redundant and paralogous EST markers and further improved the quality and utility of the genetic map for P. taeda.

Transferability of Simple Sequence Repeat (SSR) Markers Developed in Litchi chinensis to Blighia sapida (Sapindaceae)

Ackee (Blighia sapida, Sapindaceae) is a multipurpose fruit tree species of high economic importance, native to the Guinean forests of West Africa, and belongs to the same family as that of lychee (Litchi chinensis). In this study, a set of 12 primer pairs for simple sequence repeats (SSRs) previously developed for lychee has been evaluated for polymorphism in 16 ackee trees from different populations. Seven primer pairs have been found to be transferable, and four have revealed polymorphisms. However, the average number of alleles per locus has dropped from 4.9 for lychee to 3.7 for ackee. Characterization of the four polymorphic markers in 279 individuals belonging to14 different ackee populations from Benin has revealed that the numbers of alleles per locus range from two to 14 with a mean number of 5.8. The observed and expected heterozygosities range between 0.020 to 0.359 and 0.020 to 0.396, respectively.

Validating a QTL region characterized by multiple haplotypes

Validation of a quantitative trait locus (QTL) for outcrossing perennial plants is rarely reported due to complexity of plausible genetic models and reliance on field designs already available. Here, a particular marker-QTL haplotype exerted a large, positive effect on height for Pinus taeda and its origin could be traced to a founder, GP(3), in a three-generation QTL pedigree. To validate this QTL effect, we used an extended GP(3)-based pedigree. In the validation cross, each of the 46 offspring was clonally propagated from developing seeds using somatic embryogenesis technology. Subsequent analyses were conducted separately for seedlings and for other somatic emblings. For seedlings, the original QTL effect could not be fully validated. For somatic emblings, a strong negative QTL effect was detected in the validation cross; some evidence from another cross-supported the original positive QTL effect. From this part of the analysis, three distinct marker-QTL haplotypes at a single locus could be inferred. Validating QTL haplotypes in readily available field tests was feasible despite the genetic model complexity inherent to outcrossing long-lived perennials.

Evolution of a triplet repeat in a conifer

The opportunity to trace the evolution of a triplet repeat is rare, especially for seed-plant lineages with a well-defined fossil record. Microsatellite PtTX2133 sequences from 18 species in 2 conifer genera were used to calibrate the birth of a CAGn repeat, from its protomicrosatellite origins to its repeat expansion. Birth occurred in the hard-pine genome ~ 136 million years ago, or 14 million generations ago, then expanded as a polymorphic triplet repeat 136-100 million years before a major North American vicariance event. Calibration of the triplet-repeat birth and expansion is supported by the shared allelic lineages among Old and New World hard pines and the shared alleles solely among North American diploxylon or hard pines. Five CAGn repeat units appeared to be the expansion threshold for Old and New World diploxylon pines. Haploxylon pine species worldwide did not undergo birth and repeat expansion, remaining monomorphic, with a single imperfect 198-bp allele. A sister genus, Picea, had only a region of cryptic simplicity, preceding a proto-microsatellite region. The polymorphic triplet repeat in hard pines is older than some long-lived microsatellites reported for reptiles, yet younger than those reported for insects. Some cautionary points are raised about phylogenetic applications for this long-lived microsatellite.

Fossil Calibration of Molecular Divergence Infers a Moderate Mutation Rate and Recent Radiations for Pinus

Silent mutation rate estimates for Pinus vary 50-fold, ranging from angiosperm-like to among the slowest reported for plants. These differences either reflect extraordinary genomic processes or inconsistent fossil calibration, and they have important consequences for population and biogeographical inferences. Here we estimate mutation rates from 4 Pinus species that represent the major lineages using 11 nuclear and 4 chloroplast loci. Calibration was tested at the divergence of Pinus subgenera with the oldest leaf fossil from subg. Strobus (Eocene; 45 MYA) or a recently published subg. Strobus wood fossil (Cretaceous; 85 MYA). These calibrations place the origin of Pinus 190-102 MYA and give absolute silent rate estimates of 0.70-1.31x10(-9) and 0.22-0.42x10(-9).site-1.year-1 for the nuclear and chloroplast genomes, respectively. These rates are approximately 4- to 20-fold slower than angiosperms, but unlike many previous estimates, they are more consistent with the high per-generation deleterious mutation rates observed in pines. Chronograms from nuclear and chloroplast genomes show that the divergence of subgenera accounts for about half of the time since Pinus diverged from Picea, with subsequent radiations occurring more recently. By extending the sampling to encompass the phylogenetic diversity of Pinus, we predict that most extant subsections diverged during the Miocene. Moreover, subsect. Australes, Ponderosae, and Contortae, containing over 50 extant species, radiated within a 5 Myr time span starting as recently as 18 MYA. An Eocene divergence of pine subgenera (using leaf fossils) does not conflict with fossil-based estimates of the Pinus-Picea split, but a Cretaceous divergence using wood fossils accommodates Oligocene fossils that may represent modern subsections. Because homoplasy and polarity of character states have not been tested for fossil pine assignments, the choice of fossil and calibration node represents a significant source of uncertainty. Based on several lines of evidence (including agreement with ages inferred using calibrations outside of Pinus), we conclude that the 85 MYA calibration at the divergence of pine subgenera provides a reasonable lower bound and that further refinements in age and mutation rate estimates will require a synthetic examination of pine fossil history.

High degree of conservation of nuclear microsatellite loci in the genus Clusia

In plants, microsatellites and their flanking DNA are rarely conserved across a whole genus, let alone other genera in the same family. Therefore, the possibility of using microsatellite primers developed for a species across a large number of plant species in the same genus is often limited. Remarkably, dinucleotide nuclear microsatellites developed for Clusia minor and for Clusia nemorosa amplified homologous microsatellites in species across the whole genus Clusia. In this present study, we report on the DNA sequence variation across the genus of 3 microsatellite loci with varying levels of variation. Compared over the species, there was a correlation between the lengths of the microsatellite loci. Interrupts occurred multiple times and did not seem to lead to the death of the microsatellite. These highly conserved markers will be useful for studying the variable reproductive systems in the genus Clusia.

Pine microsatellite markers allow roots and ectomycorrhizas to be linked to individual trees

Linking roots and ectomycorrhizas (EcM) to individual host trees in the field is required to test whether individual trees support different ectomycorrhizal communities. Here we describe a method that identifies the source of EcM roots by PCR of polymorphic pine nuclear microsatellite loci using fluorescently labelled primers and high-throughput fragment analysis. ITS-PCR can also be performed on the same EcM DNA extract for fungal identification. The method was tested on five neighbouring Scots pine (Pinus sylvestris var scotica) trees in native woodland. Successful host tree identification from DNA extracts of EcM root tips was achieved for 93% of all root fragments recovered from soil cores. It was estimated that each individual mature pine sampled was colonised by between 15 and 19 EcM fungi. The most abundant fungal species were found on all five trees, and within the constraints of the sampling scheme, no differences between trees in EcM fungal community structure or composition were detected.

An ancient bottleneck in the Lost Pines of central Texas

The retreating edge hypothesis for species responding to climate change predicts severe bottlenecks and eventual extinction. The disjunct Lost Pines population at the westernmost edge of the widespread Pinus taeda range is well suited for testing this prediction. The occurrence of one or more genetic bottlenecks in the Lost Pines population was tested using 34 nuclear microsatellite markers and a control sample from the larger, more continuous east Texas P. taeda forests. The Lost Pines population has undergone drastic contractions in effective population size between 3000 and 30 000 years bp. These results were supported by: (i) detection of transient heterozygosity excess, (ii) a mode-shift indicator of allele frequencies, and (iii) a ratio of allele number to allele size range. No bottleneck was detected for the east Texas control using any of the three methods. The distribution of allele frequencies was skewed for the Lost Pines population compared to the control, indicating a loss of rare alleles. However, allelic diversity was similar between the Lost Pines population and its east Texas control; the mean allele number per locus was 5.29 and 5.38, respectively. It is proposed that the Lost Pines population was the western refugium for P. taeda during Pleistocene glaciation and that East Texas P. taeda forests descended from the bottlenecked Lost Pines population.

Cross-species transferability and mapping of genomic and cDNA SSRs in pines

Two unigene datasets of Pinus taeda and Pinus pinaster were screened to detect di-, tri- and tetranucleotide repeated motifs using the SSRIT script. A total of 419 simple sequence repeats (SSRs) were identified, from which only 12.8% overlapped between the two sets. The position of the SSRs within their coding sequences were predicted using FrameD. Trinucleotides appeared to be the most abundant repeated motif (63 and 51% in P. taeda and P. pinaster, respectively) and tended to be found within translated regions (76% in both species), whereas dinucleotide repeats were preferentially found within the 5'- and 3'-untranslated regions (75 and 65%, respectively). Fifty-three primer pairs amplifying a single PCR fragment in the source species (mainly P. taeda), were tested for amplification in six other pine species. The amplification rate with other pine species was high and corresponded with the phylogenetic distance between species, varying from 64.6% in P. canariensis to 94.2% in P. radiata. Genomic SSRs were found to be less transferable; 58 of the 107 primer pairs (i.e. 54%) derived from P. radiata amplified a single fragment in P. pinaster. Nine cDNA-SSRs were located to their chromosomes in two P. pinaster linkage maps. The level of polymorphism of these cDNA-SSRs was compared to that of previously and newly developed genomic-SSRs. Overall, genomic SSRs tend to perform better in terms of heterozygosity and number of alleles. This study suggests that useful SSR markers can be developed from pine ESTs.

Extracting DNA from submerged pine wood

A DNA extraction protocol for submerged pine logs was developed with the following properties: (i) high molecular weight DNA, (ii) PCR amplification of chloroplast and nuclear sequences, and (iii) high sequence homology to voucher pine specimens. The DNA extraction protocol was modified from a cetyltrimehtylammonium bromide (CTAB) protocol by adding stringent electrophoretic purification, proteinase K, RNAse, polyvinyl pyrrolidone (PVP), and Gene Releaser®. Chloroplast rbcL (ribulose-1,5-bisphosphate carboxylase) could be amplified. Nuclear ribosomal sequences had >95% homology to Pinus taeda and Pinus palustris. Microsatellite polymorphism for PtTX2082 matched 2 of 14 known P. taeda alleles. Our results show DNA analysis for submerged conifer wood is feasible.Key words: conifers, wood, polymerase chain reaction, sequencing.

Single-copy, species-transferable microsatellite markers developed from loblolly pine ESTs

Microsatellites, or simple sequence repeats (SSRs), are usually regarded as the "markers of choice" in population genetics research because they exhibit high variability. The development cost of these markers is usually high. In addition, microsatellite primers developed for one species often do not cross-amplify in related species, requiring separate development for each species. However, microsatellites found in expressed sequence tags (ESTs) might better cross-amplify as they reside in or near conserved coding DNA. In this study, we identified 14 Pinus taeda (loblolly pine) EST-SSRs from public EST databases and tested for their cross-species transferability to P. contorta ssp. latifolia, P. ponderosa, and P. sylvestris. As part of our development of a P. contorta microsatellite set, we also compared their transferability to that of 99 traditional microsatellite markers developed in P. taeda and tested on P. contorta ssp. latifolia. Compared to traditional microsatellites, EST-SSRs had higher transfer rates across pine species; however, the level of polymorphism of microsatellites derived from ESTs was lower. Sequence analyses revealed that the frequencies of insertions/deletions and base substitutions were lower in EST-SSRs than in other types of microsatellites, confirming that EST-SSRs are more conserved than traditional SSRs. Our results also provide a battery of 23 polymorphic, robust microsatellite primer pairs for lodgepole pine.

Cross-amplification and sequence variation of microsatellite loci in Eurasian hard pines

Microsatellite transfer across coniferous species is a valued methodology because de novo development for each species is costly and there are many species with only a limited commodity value. Cross-species amplification of orthologous microsatellite regions provides valuable information on mutational and evolutionary processes affecting these loci. We tested 19 nuclear microsatellite markers from Pinus taeda L. (subsection Australes) and three from P. sylvestris L. (subsection Pinus) on seven Eurasian hard pine species ( P. uncinata Ram., P. sylvestris L., P. nigra Arn., P. pinaster Ait., P. halepensis Mill., P. pinea L. and P. canariensis Sm.). Transfer rates to species in subsection Pinus (36-59%) were slightly higher than those to subsections Pineae and Pinaster (32-45%). Half of the trans-specific microsatellites were found to be polymorphic over evolutionary times of approximately 100 million years (ten million generations). Sequencing of three trans-specific microsatellites showed conserved repeat and flanking regions. Both a decrease in the number of perfect repeats in the non-focal species and a polarity for mutation, the latter defined as a higher substitution rate in the flanking sequence regions close to the repeat motifs, were observed in the trans-specific microsatellites. The transfer of microsatellites among hard pine species proved to be useful for obtaining highly polymorphic markers in a wide range of species, thereby providing new tools for population and quantitative genetic studies.

Conserved simple sequence repeats for the Limnanthaceae (Brassicales)

The Limnanthaceae (Order Brassicales) is a family of 18 taxa of Limnanthes (meadowfoam) native to California, Oregon, and British Columbia. Cultivated meadowfoam ( L. alba Benth.), a recently domesticated plant, has been the focus of research and development as an industrial oilseed for three decades. The goal of the present research was to develop several hundred simple sequence repeat (SSR) markers for genetic mapping, molecular breeding, and genomics research in wild and cultivated meadowfoam taxa. We developed 389 SSR markers for cultivated meadowfoam by isolating and sequencing 1,596 clones from L. alba genomic DNA libraries enriched for AG(n) or AC(n) repeats, identifying one or more unique SSRs in 696 clone sequences, and designing and testing primers for 624 unique SSRs. The SSR markers were screened for cross- taxa utility and polymorphisms among ten of 17 taxa in the Limnanthaceae; 373 of these markers were polymorphic and 106 amplified loci from every taxon. Cross-taxa amplification percentages ranged from 37.3% in L. douglasii ssp. rosea (145/389) to 85.6% in L. montana (333/389). The SSR markers amplified 4,160 unique bands from 14 genotypes sampled from ten taxa (10.7 unique bands per SSR marker), of which 972 were genotype-specific. Mean and maximum haplotype heterozygosities were 0.71 and 0.90, respectively, among six L. alba genotypes and 0.63 and 0.93, respectively, among 14 genotypes (ten taxa). The SSR markers supply a critical mass of high-throughput DNA markers for biological and agricultural research across the Limnanthaceae and open the way to the development of a genetic linkage map for meadowfoam ( x = 5).

Comparative mapping between quercus and castanea using simple-sequence repeats (SSRs)

Simple sequence repeat (SSR) markers from Quercus and Castanea were used for comparative mapping between Quercus robur (L.) and Castanea sativa (Mill.). We tested the transferability of SSRs developed in Quercus to Castanea and vice-versa. In total, 47% (25) of the Quercus SSRs and 63% (19) of the Castanea SSRs showed a strong amplification product in the non-source species. From these 44 putative comparative anchor tags, 19 (15 from Quercus and 4 from Castanea) were integrated in two previously established genetic linkage maps for the two genera. SSR loci were sequenced to confirm the orthology of the markers. The combined information from both genetic mapping and sequence analysis were used to determine the homeology between seven linkage groups, aligned on the basis of pairs or triplets of common markers, while two additional groups were matched using a single microsatellite marker. Orthologous loci identified between Q. robur and C. sativa will be useful as anchor loci for comparative mapping studies within the Fagaceae family.

No clustering for linkage map based on low-copy and undermethylated microsatellites

Clustering has been reported for conifer genetic maps based on hypomethylated or low-copy molecular markers, resulting in uneven marker distribution. To test this, a framework genetic map was constructed from three types of microsatellites: low-copy, undermethylated, and genomic. These Pinus taeda L. microsatellites were mapped using a three-generation pedigree with 118 progeny. The microsatellites were highly informative; of the 32 markers in intercross configuration, 29 were segregating for three or four alleles in the progeny. The sex-averaged map placed 51 of the 95 markers in 15 linkage groups at LOD > 4.0. No clustering or uneven distribution across the genome was observed. The three types of P. taeda microsatellites were randomly dispersed within each linkage group. The 51 microsatellites covered a map distance of 795 cM, an average distance of 21.8 cM between markers, roughly half of the estimated total map length. The minimum and maximum distances between any two bins was 4.4 and 45.3 cM, respectively. These microsatellites provided anchor points for framework mapping for polymorphism in P. taeda and other closely related hard pines.Key words: simple sequence repeats, genome duplication, conifers, Pinus taeda L. hypomethylation, low-copy kinetic component.

Cross-species amplification from crop soybean Glycine max provides informative microsatellite markers for the study of inbreeding wild relatives

The development of microsatellite markers through transfer of primers from related species (cross-species amplification) remains a little-explored alternative to the de novo method in plants. In this study of 100 microsatellite loci from Glycine max, we examined two aspects of primer transfer. First, we tested if source locus properties can predict primer transfer and polymorphism in Glycine cyrtoloba and Glycine clandestina. We transferred 23 primers to G. cyrtoloba and 42 to G. clandestina, with 19 loci polymorphic within G. clandestina. However, we could not predict transfer or polymorphism from the source locus properties. Second, we evaluated the subset of 11 polymorphic loci for study in G. clandestina populations representing two local morphotypes. All loci were informative within populations (population mean He +/- SE = 0.58 +/- 0.04). We directly sequenced 28 alleles at 4 representative loci. The allelic patterns and sequencing results established that 8 of 11 loci were typical microsatellites, confirming the utility of primer transfer as an alternative to de novo development. Additionally, we found that morphotypic differentiation between populations was paralleled by changes in polymorphism level at six loci and size homoplasy at one locus. We interpret these patterns as being a product of selfing in G. clandestina. Our results demonstrate the value of allele sequence knowledge for the most effective use of microsatellites.

Isolation, characterization, and cross-species utility of microsatellites in yellow cedar (Chamaecyparis nootkatensis)

Chamaecyparis nootkatensis is an ecologically and economically important conifer of the north Pacific coastal forests. To aid in studies of clonal structure and genetic differentiation of this and related species, we isolated and characterized microsatellites from C. nootkatensis. A microsatellite-enriched library yielded 75 repeat-containing sequences for which primer pairs were designed. Only five showed reliable amplification and polymorphism, with an average of 13.7 alleles/locus and a mean expected heterozygosity of 0.592. In progeny tests with four families, few null alleles were directly detected and loci segregated according to Mendelian expectations. However, in one primer pair, high heterozygote deficiency was observed, suggesting the presence of a null allele. The ability of primer pairs to cross amplify was tested on 18 species of the Cupressaceae sensu lato; three primer pairs yielded polymorphic loci in Cupressus and Juniperus species, but not in other Chamaecyparis species. This also supports recent findings of a closer affinity of C. nootkatensis with Cupressus over other Chamaecyparis species.