Undermethylated DNA as a source of microsatellites from a conifer genome

Molecular characterisation of Pinus sylvestris (L.) in Ireland at the western limit of the species distribution

BACKGROUND

Scots pine (Pinus sylvestris L.) underwent significant population declines across much of northwest Europe during the mid-to-late Holocene and was thought to have become extirpated in Ireland from about 400 AD. However, most extant populations are plantations reintroduced from Scotland. Others are naturalised therefrom and one in Western Ireland is a putative relict. In this paper, Scots pine in Ireland are genetically described for the first time.

RESULTS

Using two mitochondrial (mtDNA) loci, eight chloroplast (cpSSR) and 18 nuclear (nSSR) loci, the genetic composition and diversity of 19 Irish Scots pine populations is described and compared to other European populations. All trees sampled in Ireland were fixed for mitotype a, which is the most common across northwest Europe. By contrast, cpSSR (HCP = 0.967) and nSSR (He = 0.540) variation was high, and comparable with estimates for other regions across the species range. Differentiation at both sets of loci were similarly low (cpSSR FST = 0.019; nSSR FST = 0.018), but populations from continental Europe were significantly differentiated from all Irish populations based on nSSR variation.

CONCLUSIONS

All Irish Scots pine are likely part of a common Irish-Scottish gene pool which diverged from continental Scots pine following post-glacial recolonisation. A high genetic diversity and an absence of evidence of inbreeding suggests the regional decline of Scots pine did not critically reduce allelic variation. The post-glacial relationship between Irish and Scottish pine is discussed, and a suggestion from recent palaeoecological work that reintroduced Scots pine be managed as a native species is now further supported by genetic data.

Measuring DNA modifications with the comet assay: a compendium of protocols

The comet assay is a versatile method to detect nuclear DNA damage in individual eukaryotic cells, from yeast to human. The types of damage detected encompass DNA strand breaks and alkali-labile sites (e.g., apurinic/apyrimidinic sites), alkylated and oxidized nucleobases, DNA-DNA crosslinks, UV-induced cyclobutane pyrimidine dimers and some chemically induced DNA adducts. Depending on the specimen type, there are important modifications to the comet assay protocol to avoid the formation of additional DNA damage during the processing of samples and to ensure sufficient sensitivity to detect differences in damage levels between sample groups. Various applications of the comet assay have been validated by research groups in academia, industry and regulatory agencies, and its strengths are highlighted by the adoption of the comet assay as an in vivo test for genotoxicity in animal organs by the Organisation for Economic Co-operation and Development. The present document includes a series of consensus protocols that describe the application of the comet assay to a wide variety of cell types, species and types of DNA damage, thereby demonstrating its versatility.

The development and validation of EpiComet-Chip, a modified high-throughput comet assay for the assessment of DNA methylation status

DNA damage and alterations in global DNA methylation status are associated with multiple human diseases and are frequently correlated with clinically relevant information. Therefore, assessing DNA damage and epigenetic modifications, including DNA methylation, is critical for predicting human exposure risk of pharmacological and biological agents. We previously developed a higher-throughput platform for the single cell gel electrophoresis (comet) assay, CometChip, to assess DNA damage and genotoxic potential. Here, we utilized the methylation-dependent endonuclease, McrBC, to develop a modified alkaline comet assay, "EpiComet," which allows single platform evaluation of genotoxicity and global DNA methylation [5-methylcytosine (5-mC)] status of single-cell populations under user-defined conditions. Further, we leveraged the CometChip platform to create an EpiComet-Chip system capable of performing quantification across simultaneous exposure protocols to enable unprecedented speed and simplicity. This system detected global methylation alterations in response to exposures which included chemotherapeutic and environmental agents. Using EpiComet-Chip on 63 matched samples, we correctly identified single-sample hypermethylation (≥1.5-fold) at 87% (20/23), hypomethylation (≥1.25-fold) at 100% (9/9), with a 4% (2/54) false-negative rate (FNR), and 10% (4/40) false-positive rate (FPR). Using a more stringent threshold to define hypermethylation (≥1.75-fold) allowed us to correctly identify 94% of hypermethylation (17/18), but increased our FPR to 16% (7/45). The successful application of this novel technology will aid hazard identification and risk characterization of FDA-regulated products, while providing utility for investigating epigenetic modes of action of agents in target organs, as the assay is amenable to cultured cells or nucleated cells from any tissue. Environ. Mol. Mutagen. 58:508-521, 2017. © 2017 Wiley Periodicals, Inc.

Lineage-Specific Methyltransferases Define the Methylome of the Globally Disseminated Escherichia coli ST131 Clone

Escherichia coli sequence type 131 (ST131) is a clone of uropathogenic E. coli that has emerged rapidly and disseminated globally in both clinical and community settings. Members of the ST131 lineage from across the globe have been comprehensively characterized in terms of antibiotic resistance, virulence potential, and pathogenicity, but to date nothing is known about the methylome of these important human pathogens. Here we used single-molecule real-time (SMRT) PacBio sequencing to determine the methylome of E. coli EC958, the most-well-characterized completely sequenced ST131 strain. Our analysis of 52,081 methylated adenines in the genome of EC958 discovered three (m6)A methylation motifs that have not been described previously. Subsequent SMRT sequencing of isogenic knockout mutants identified the two type I methyltransferases (MTases) and one type IIG MTase responsible for (m6)A methylation of novel recognition sites. Although both type I sites were rare, the type IIG sites accounted for more than 12% of all methylated adenines in EC958. Analysis of the distribution of MTase genes across 95 ST131 genomes revealed their prevalence is highly conserved within the ST131 lineage, with most variation due to the presence or absence of mobile genetic elements on which individual MTase genes are located.

IMPORTANCE

DNA modification plays a crucial role in bacterial regulation. Despite several examples demonstrating the role of methyltransferase (MTase) enzymes in bacterial virulence, investigation of this phenomenon on a whole-genome scale has remained elusive until now. Here we used single-molecule real-time (SMRT) sequencing to determine the first complete methylome of a strain from the multidrug-resistant E. coli sequence type 131 (ST131) lineage. By interrogating the methylome computationally and with further SMRT sequencing of isogenic mutants representing previously uncharacterized MTase genes, we defined the target sequences of three novel ST131-specific MTases and determined the genomic distribution of all MTase target sequences. Using a large collection of 95 previously sequenced ST131 genomes, we identified mobile genetic elements as a major factor driving diversity in DNA methylation patterns. Overall, our analysis highlights the potential for DNA methylation to dramatically influence gene regulation at the transcriptional level within a well-defined E. coli clone.

Highly informative single-copy nuclear microsatellite DNA markers developed using an AFLP-SSR approach in black spruce (Picea mariana) and red spruce (P. rubens)

Background

Microsatellites or simple sequence repeats (SSRs) are highly informative molecular markers for various biological studies in plants. In spruce (Picea) and other conifers, the development of single-copy polymorphic genomic microsatellite markers is quite difficult, owing primarily to the large genome size and predominance of repetitive DNA sequences throughout the genome. We have developed highly informative single-locus genomic microsatellite markers in black spruce (Picea mariana) and red spruce (Picea rubens) using a simple but efficient method based on a combination of AFLP and microsatellite technologies.

Principal Findings

A microsatellite-enriched library was constructed from genomic AFLP DNA fragments of black spruce. Sequencing of the 108 putative SSR-containing clones provided 94 unique sequences with microsatellites. Twenty-two of the designed 34 primer pairs yielded scorable amplicons, with single-locus patterns. Fourteen of these microsatellite markers were characterized in 30 black spruce and 30 red spruce individuals drawn from many populations. The number of alleles at a polymorphic locus ranged from 2 to 18, with a mean of 9.3 in black spruce, and from 3 to 15, with a mean of 6.2 alleles in red spruce. The polymorphic information content or expected heterozygosity ranged from 0.340 to 0.909 (mean = 0.67) in black spruce and from 0.161 to 0.851 (mean = 0.62) in red spruce. Ten SSR markers showing inter-parental polymorphism inherited in a single-locus Mendelian mode, with two cases of distorted segregation. Primer pairs for almost all polymorphic SSR loci resolved microsatellites of comparable size in Picea glauca, P. engelmannii, P. sitchensis, and P. abies.

Significance

The AFLP-based microsatellite-enriched library appears to be a rapid, cost-effective approach for isolating and developing single-locus informative genomic microsatellite markers in black spruce. The markers developed should be useful in black spruce, red spruce and other Picea species for various genetics, genomics, breeding, forensics, conservation studies and applications.

Comparative mapping among subsectionAustrales(genusPinus, family Pinaceae)

Comparative mapping in conifers has not yet been used to test for small-scale genomic disruptions such as inversions, duplications, and deletions occurring between closely related taxa. Using comparative mapping to probe this smaller scale of inquiry may provide clues about speciation in a phylogenetically problematic taxon, the diploxylon pine subsection Australes (genus Pinus , family Pinaceae). Genetic maps were constructed for two allopatric species of Australes, P. elliottii var. elliottii and P. caribaea var. hondurensis , using microsatellites and an F1hybrid. A third map was generated directly from the meiotic products of an adult F1hybrid, eliminating the need for an F2generation. Numerous small-scale disruptions were detected in addition to synteny and collinearity, and these included (1) map shrinkage, (2) a paracentric inversion, (3) transmission ratio distortion, and (4) mild selection against a parental haplotype. Such cryptic signatures of genomic divergence between closely related interfertile species are useful in elucidating this problematic evolutionary history.

Genome-wide analysis of DNA methylation patterns

Cytosine methylation is the most common covalent modification of DNA in eukaryotes. DNA methylation has an important role in many aspects of biology, including development and disease. Methylation can be detected using bisulfite conversion, methylation-sensitive restriction enzymes, methyl-binding proteins and anti-methylcytosine antibodies. Combining these techniques with DNA microarrays and high-throughput sequencing has made the mapping of DNA methylation feasible on a genome-wide scale. Here we discuss recent developments and future directions for identifying and mapping methylation, in an effort to help colleagues to identify the approaches that best serve their research interests.

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.

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.

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.

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.

Gene enrichment in plant genomic shotgun libraries

The Arabidopsis genome (about 130 Mbp) has been completely sequenced; whereas a draft sequence of the rice genome (about 430 Mbp) is now available and the sequencing of this genome will be completed in the near future. The much larger genomes of several important crop species, such as wheat (about 16,000 Mbp) or maize (about 2500 Mbp), may not be fully sequenced with current technology. Instead, sequencing-analysis strategies are being developed to obtain sequencing and mapping information selectively for the genic fraction (gene space) of complex plant genomes.