Konstantin V. Krutovsky

Association Genetics of Coastal Douglas Fir (Pseudotsuga menziesii var. menziesii, Pinaceae). I. Cold-Hardiness Related Traits

Adaptation to cold is one of the greatest challenges to forest trees. This process is highly synchronized with environmental cues relating to photoperiod and temperature. Here, we use a candidate gene-based approach to search for genetic associations between 384 single-nucleotide polymorphism (SNP) markers from 117 candidate genes and 21 cold-hardiness related traits. A general linear model approach, including population structure estimates as covariates, was implemented for each marker–trait pair. We discovered 30 highly significant genetic associations [false discovery rate (FDR) Q < 0.10] across 12 candidate genes and 10 of the 21 traits. We also detected a set of 7 markers that had elevated levels of differentiation between sampling sites situated across the Cascade crest in northeastern Washington. Marker effects were small (r2 < 0.05) and within the range of those published previously for forest trees. The derived SNP allele, as measured by a comparison to a recently diverged sister species, typically affected the phenotype in a way consistent with cold hardiness. The majority of markers were characterized as having largely nonadditive modes of gene action, especially underdominance in the case of cold-tolerance related phenotypes. We place these results in the context of trade-offs between the abilities to grow longer and to avoid fall cold damage, as well as putative epigenetic effects. These associations provide insight into the genetic components of complex traits in coastal Douglas fir, as well as highlight the need for landscape genetic approaches to the detection of adaptive genetic diversity.

Nucleotide Diversity and Linkage Disequilibrium in Cold-Hardiness- and Wood Quality-Related Candidate Genes in Douglas Fir

Nuclear sequence variation and linkage disequilibrium (LD) were studied in 15 cold-hardiness- and 3 wood quality-related candidate genes in Douglas fir [Pseudotsuga menziesii (Mirb.) Franco]. This set of genes was selected on the basis of its function in other plants and collocation with cold-hardiness-related quantitative trait loci (QTL). The single-nucleotide polymorphism (SNP) discovery panel represented 24 different trees from six regions in Washington and Oregon plus parents of a segregating population used in the QTL study. The frequency of SNPs was one SNP per 46 bp across coding and noncoding regions on average. Haplotype and nucleotide diversities were also moderately high with Hd = 0.827 ± 0.043 and π = 0.00655 ± 0.00082 on average, respectively. The nonsynonymous (replacement) nucleotide substitutions were almost five times less frequent than synonymous ones and substitutions in noncoding regions. LD decayed relatively slowly but steadily within genes. Haploblock analysis was used to define haplotype tag SNPs (htSNPs). These data will help to select SNPs for association mapping, which is already in progress.

Multilocus Patterns of Nucleotide Diversity and Divergence Reveal Positive Selection at Candidate Genes Related to Cold-hardiness in Coastal Douglas-fir (Pseudotsuga menziesii var. menziesii)

Forest trees exhibit remarkable adaptations to their environments. The genetic basis for phenotypic adaptation to climatic gradients has been established through a long history of common garden, provenance, and genecological studies. The identities of genes underlying these traits, however, have remained elusive and thus so have the patterns of adaptive molecular diversity in forest tree genomes. Here, we report an analysis of diversity and divergence for a set of 121 cold-hardiness candidate genes in coastal Douglas fir (Pseudotsuga menziesii var. menziesii). Application of several different tests for neutrality, including those that incorporated demographic models, revealed signatures of selection consistent with selective sweeps at three to eight loci, depending upon the severity of a bottleneck event and the method used to detect selection. Given the high levels of recombination, these candidate genes are likely to be closely linked to the target of selection if not the genes themselves. Putative homologs in Arabidopsis act primarily to stabilize the plasma membrane and protect against denaturation of proteins at freezing temperatures. These results indicate that surveys of nucleotide diversity and divergence, when framed within the context of further association mapping experiments, will come full circle with respect to their utility in the dissection of complex phenotypic traits into their genetic components.

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

BackgroundPrevious 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.ResultsThe 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, He , 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.ConclusionsMany 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.

Estimation of population structure in coastal Douglas-fir (Pseudotsuga menziesii (Mirb.) Franco var. menziesii) using allozyme and microsatellite markers

Characterizing population structure using neutral markers is an important first step in association genetic studies in order to avoid false associations between phenotypes and genotypes that may arise from non-selective demographic factors. Population structure was studied in a wide sample of ∼1,300 coastal Douglas-fir [Pseudotsuga menziesii (Mirb.) Franco var. menziesii] trees from Washington and Oregon. This sample is being used for association mapping between cold hardiness and phenology phenotypes and single-nucleotide polymorphisms in adaptive-trait candidate genes. All trees were genotyped for 25 allozyme and six simple sequence repeat (SSR) markers using individual megagametophytes. Population structure analysis was done separately for allozyme and SSR markers, as well as for both datasets combined. The parameter of genetic differentiation (θ or FST) was standardized to take into account high within-population variation in the SSR loci and to allow comparison with allozyme loci. Genetic distance between populations was positively and significantly correlated with geographic distance, and weak but significant clinal variation was found for a few alleles. Although the STRUCTURE simulation analysis inferred the same number of populations as used in this study and as based on previous analysis of quantitative adaptive trait variation, clustering among populations was not significant. In general, results indicated weak differentiation among populations for both allozyme and SSR loci (θs = 0.006–0.059). The lack of pronounced population subdivision in the studied area should facilitate association mapping in this experimental population, but we recommend taking the STRUCTURE analysis and population assignments for individual trees (Q-matrix) into account in association mapping.

Evolution of exon-intron structure and alternative splicing.

Despite significant advances in high-throughput DNA sequencing, many important species remain understudied at the genome level. In this study we addressed a question of what can be predicted about the genome-wide characteristics of less studied species, based on the genomic data from completely sequenced species. Using NCBI databases we performed a comparative genome-wide analysis of such characteristics as alternative splicing, number of genes, gene products and exons in 36 completely sequenced model species. We created statistical regression models to fit these data and applied them to loblolly pine (Pinus taeda L.), an example of an important species whose genome has not been completely sequenced yet. Using these models, the genome-wide characteristics, such as total number of genes and exons, can be roughly predicted based on parameters estimated from available limited genomic data, e.g. exon length and exon/gene ratio.

Genetic structure and association mapping of adaptive and selective traits in the east Texas loblolly pine (Pinus taeda L.) breeding populations

First-generation selection (FGS) and second-generation selection (SGS) breeding populations of loblolly pine from east Texas were studied to estimate the genetic diversity, population structure, linkage disequilibrium (LD), signatures of selection and association of breeding traits with a genome-wide panel of 4,264 single nucleotide polymorphisms (SNPs). Relatively high levels of observed (Ho = 0.178–0.198) and expected (He = 0.180–0.198) heterozygosities were observed in all populations. The amount of inbreeding was very low with many populations exhibiting a slight excess of heterozygotes. The population structure was weak, but FST indicated more pronounced differentiation in the SGS populations. As expected for outcrossing natural populations, the genome-wide LD was low, but marker density was insufficient to deduce the decay rate. Numerous associations were found between various phenotypic traits and SNPs, but only a few remained significant after false positive correction. Signatures of diversifying and balancing selection were found in markers representing important biological functions. These results present the first step in the application of marker-assisted selection (MAS) to the Western Gulf Forest Tree Improvement Program (WGFTIP) for loblolly pine and will contribute to the knowledgebase necessary for genomic selection technology.

Conserved ortholog sets in forest trees

Putative single-copy genes and conserved ortholog sets (COS) were identified in model plant species thale cress (Arabidopsis thaliana), rice (Oryza sativa ssp. japonica), and poplar [black cottonwood, Populus trichocarpa (Torr. & Gray ex Brayshaw)] and used to find putative COS in four conifers (the Coniferales order). Using expressed sequence tag sequences, unique transcript sets were assembled in loblolly pine (Pinus taeda L.), white spruce [Picea glauca (Moench) Voss], Douglas-fir [Pseudotsuga menziesii (Mirb.) Franco var. menziesii], and sugi [Cryptomeria japonica (Thunberg ex Linnaeus f.) D. Don]. They were compared with COS sets identified in three model plant species using comparative sequence analysis. Almost half of the single-copy genes in herbaceous species (Arabidopsis and rice) had additional copies and homologs in poplar and conifers. The identified tentative COS sets have many applications in evolutionary genomics studies, phylogenetic analysis, and comparative mapping.

Development of Genetic Markers in Eucalyptus Species by Target Enrichment and Exome Sequencing

The advent of next-generation sequencing has facilitated large-scale discovery, validation and assessment of genetic markers for high density genotyping. The present study was undertaken to identify markers in genes supposedly related to wood property traits in three Eucalyptus species. Ninety four genes involved in xylogenesis were selected for hybridization probe based nuclear genomic DNA target enrichment and exome sequencing. Genomic DNA was isolated from the leaf tissues and used for on-array probe hybridization followed by Illumina sequencing. The raw sequence reads were trimmed and high-quality reads were mapped to the E. grandis reference sequence and the presence of single nucleotide variants (SNVs) and insertions/ deletions (InDels) were identified across the three species. The average read coverage was 216X and a total of 2294 SNVs and 479 InDels were discovered in E. camaldulensis, 2383 SNVs and 518 InDels in E. tereticornis, and 1228 SNVs and 409 InDels in E. grandis. Additionally, SNV calling and InDel detection were conducted in pair-wise comparisons of E. tereticornis vs. E. grandis, E. camaldulensis vs. E. tereticornis and E. camaldulensis vs. E. grandis. This study presents an efficient and high throughput method on development of genetic markers for family– based QTL and association analysis in Eucalyptus.

From population genetics to population genomics of forest trees: Integrated population genomics approach

Early works by Altukhov and his associates on pine and spruce laid the foundation for Russian population genetic studies on tree species with the use of molecular genetic markers. In recent years, these species have become especially popular as nontraditional eukaryotic models for population and evolutionary genome-wide research. Tree species with large, cross-pollinating native populations, high genetic and phenotypic variation, growing in diverse environments and affected by environmental changes during hundreds of years of their individual development, are an ideal model for studying the molecular genetic basis of adaptation. The great advance in this field is due to the rapid development of population genomics in the last few years. In the broad sense, population genomics is a novel, fast-developing discipline, combining traditional population genetic approaches with the genome-wide level of analysis. Thousands of genes with known function and sometimes known genome-wide localization can be simultaneously studied in many individuals. This opens new prospects for obtaining statistical estimates for a great number of genes and segregating elements. Mating system, gene exchange, reproductive population size, population disequilibrium, interaction among populations, and many other traditional problems of population genetics can be now studied using data on variation in many genes. Moreover, population genome-wide analysis allows one to distinguish factors that affect individual genes, allelles, or nucleotides (such as, for example, natural selection) from factors affecting the entire genome (e.g., demography). This paper presents a brief review of traditional methods of studying genetic variation in forest tree species and introduces a new, integrated population genomics approach. The main stages of the latter are: (1) selection of genes, which are tentatively involved in variation of adaptive traits, by means of a detailed examination of the regulation and the expression of individual genes and genotypes, with subsequent determination of their complete allelic composition by direct nucleotide sequencing; (2) examination of the phenotypic effects of individual alleles by, e.g., association mapping; and (3) determining the frequencies of the selected alleles in natural population for identification of the adaptive variation pattern in the heterogeneous environment. Through decoding the phenotypic effects of individual alleles and identification of adaptive variation patterns at the population level, population genomics in the future will serve as a very helpful, efficient, and economical tool, essential for developing a correct strategy for conserving and increasing forests and other commercially valuable plant and animal species.