Extensive Variation in Drought-Induced Gene Expression Changes Between Loblolly Pine Genotypes

Drought response is coordinated through expression changes in a large suite of genes. Interspecific variation in this response is common and associated with drought-tolerant and -sensitive genotypes. The extent to which different genetic networks orchestrate the adjustments to water deficit in tolerant and sensitive genotypes has not been fully elucidated, particularly in non-model or woody plants. Differential expression analysis via RNA-seq was evaluated in root tissue exposed to simulated drought conditions in two loblolly pine (Pinus taeda L.) clones with contrasting tolerance to drought. Loblolly pine is the prevalent conifer in southeastern U.S. and a major commercial forestry species worldwide. Significant changes in gene expression levels were found in more than 4,000 transcripts [drought-related transcripts (DRTs)]. Genotype by environment (GxE) interactions were prevalent, suggesting that different cohorts of genes are influenced by drought conditions in the tolerant vs. sensitive genotypes. Functional annotation categories and metabolic pathways associated with DRTs showed higher levels of overlap between clones, with the notable exception of GO categories in upregulated DRTs. Conversely, both differentially expressed transcription factors (TFs) and TF families were largely different between clones. Our results indicate that the response of a drought-tolerant loblolly pine genotype vs. a sensitive genotype to water limitation is remarkably different on a gene-by-gene level, although it involves similar genetic networks. Upregulated transcripts under drought conditions represent the most diverging component between genotypes, which might depend on the activation and repression of substantially different groups of TFs.

Predicting Adaptive Genetic Variation of Loblolly Pine (Pinus taeda L.) Populations Under Projected Future Climates Based on Multivariate Models

Greenhouse gas emission and global warming are likely to cause rapid climate change within the natural range of loblolly pine over the next few decades, thus bringing uncertainty to their adaptation to the environment. Here, we studied adaptive genetic variation of loblolly pine and correlated genetic variation with bioclimatic variables using multivariate modeling methods—Redundancy Analysis, Generalized Dissimilarity Modeling, and Gradient Forests. Studied trees (N = 299) were originally sampled from their native range across eight states on the east side of the Mississippi River. Genetic variation was calculated using a total of 44,317 single-nucleotide polymorphisms acquired by exome target sequencing. The fitted models were used to predict the adaptive genetic variation on a large spatial and temporal scale. We observed east-to-west spatial genetic variation across the range, which presented evidence of isolation by distance. Different key factors drive adaptation of loblolly pine from different geographical regions. Trees residing near the northeastern edge of the range, spanning across Delaware and Maryland and mountainous areas of  Virginia, North Carolina, South Carolina, and northern Georgia, were identified to be most likely impacted by climate change based on the large difference in genetic composition under current and future climate conditions. This study provides new perspectives on adaptive genetic variation of loblolly pine in response to different climate scenarios, and the results can be used to target particular populations while developing adaptive forest management guidelines.

Detecting the genetic basis of local adaptation in loblolly pine (Pinus taeda L.) using whole exome-wide genotyping and an integrative landscape genomics analysis approach

In the Southern United States, the widely distributed loblolly pine contributes greatly to lumber and pulp production, as well as providing many important ecosystem services. Climate change may affect the productivity and range of loblolly pine. Nevertheless, we have insufficient knowledge of the adaptive potential and the genetics underlying the adaptability of loblolly pine. To address this, we tested the association of 2.8 million whole exome-based single nucleotide polymorphisms (SNPs) with climate and geographic variables, including temperature, precipitation, latitude, longitude, and elevation data. Using an integrative landscape genomics approach by combining multiple environmental association and outlier detection analyses, we identified 611 SNPs associated with 56 climate and geographic variables. Longitude, maximum temperature of the warm months and monthly precipitation associated with most SNPs, indicating their importance and complexity in shaping the genetic variation in loblolly pine. Functions of candidate genes related to terpenoid synthesis, pathogen defense, transcription factors, and abiotic stress response. We provided evidence that environment-associated SNPs also composed the genetic structure of adaptive phenotypic traits including height, diameter, metabolite levels, and gene transcript abundance. Our study promotes understanding of the genetic basis of local adaptation in loblolly pine and provides promising tools for selecting genotypes adapted to local environments in a changing climate.

Exploring the genetic basis of gene transcript abundance and metabolite levels in loblolly pine (Pinus taeda L.) using association mapping and network construction

BACKGROUND

Identifying genetic variations that shape important complex traits is fundamental to the genetic improvement of important forest tree species, such as loblolly pine (Pinus taeda L.), which is one of the most commonly planted forest tree species in the southern U.S. Gene transcripts and metabolites are important regulatory intermediates that link genetic variations to higher-order complex traits such as wood development and drought response. A few prior studies have associated intermediate phenotypes including mRNA expression and metabolite levels with a limited number of molecular markers, but the identification of genetic variations that regulate intermediate phenotypes needs further investigation.

RESULTS

We identified 1841 single nucleotide polymorphisms (SNPs) associated with 191 gene expression mRNA phenotypes and 524 SNPs associated with 53 metabolite level phenotypes using 2.8 million exome-derived SNPs. The identified SNPs reside in genes with a wide variety of functions. We further integrated the identified SNPs and the associated expressed genes and metabolites into networks. We described the SNP-SNP interactions that significantly impacted the gene transcript abundance and metabolite level in the networks. Key loci and genes in the wood development and drought response networks were identified and analyzed.

CONCLUSIONS

This work provides new candidate genes for research on the genetic basis of gene expression and metabolism linked to wood development and drought response in loblolly pine and highlights the efficiency of using association-mapping-based networks to discover candidate genes with important roles in complex biological processes.

Assessing the Gene Content of the Megagenome: Sugar Pine (Pinus lambertiana)

Sugar pine (Pinus lambertiana Douglas) is within the subgenus Strobus with an estimated genome size of 31 Gbp. Transcriptomic resources are of particular interest in conifers due to the challenges presented in their megagenomes for gene identification. In this study, we present the first comprehensive survey of the P. lambertiana transcriptome through deep sequencing of a variety of tissue types to generate more than 2.5 billion short reads. Third generation, long reads generated through PacBio Iso-Seq have been included for the first time in conifers to combat the challenges associated with de novo transcriptome assembly. A technology comparison is provided here to contribute to the otherwise scarce comparisons of second and third generation transcriptome sequencing approaches in plant species. In addition, the transcriptome reference was essential for gene model identification and quality assessment in the parallel project responsible for sequencing and assembly of the entire genome. In this study, the transcriptomic data were also used to address questions surrounding lineage-specific Dicer-like proteins in conifers. These proteins play a role in the control of transposable element proliferation and the related genome expansion in conifers.

Exome genotyping, linkage disequilibrium and population structure in loblolly pine (Pinus taeda L.)

Background

Loblolly pine (Pinus taeda L.) is one of the most widely planted and commercially important forest tree species in the USA and worldwide, and is an object of intense genomic research. However, whole genome resequencing in loblolly pine is hampered by its large size and complexity and a lack of a good reference. As a valid and more feasible alternative, entire exome sequencing was hence employed to identify the gene-associated single nucleotide polymorphisms (SNPs) and to genotype the sampled trees.

Results

The exons were captured in the ADEPT2 association mapping population of 375 clonally-propagated loblolly pine trees using NimbleGen oligonucleotide hybridization probes, and then exome-enriched genomic DNA fragments were sequenced using the Illumina HiSeq 2500 platform. Oligonucleotide probes were designed based on 199,723 exons (≈49 Mbp) partitioned from the loblolly pine reference genome (PineRefSeq v. 1.01). The probes covered 90.2 % of the target regions. Capture efficiency was high; on average, 67 % of the sequence reads generated for each tree could be mapped to the capture target regions, and more than 70 % of the captured target bases had at least 10X sequencing depth per tree. A total of 972,720 high quality SNPs were identified after filtering. Among them, 53 % were located in coding regions (CDS), 5 % in 5’ or 3’ untranslated regions (UTRs) and 42 % in non-target and non-coding regions, such as introns and adjacent intergenic regions collaterally captured. We found that linkage disequilibrium (LD) decayed very rapidly, with the correlation coefficient (r²) between pairs of SNPs linked within single scaffolds decaying to half maximum (r² = 0.22) within 55 bp, to r² = 0.1 within 192 bp, and to r² = 0.05 within 451 bp. Population structure analysis using unlinked SNPs demonstrated the presence of two main distinct clusters representing western and eastern parts of the loblolly pine range included in our sample of trees.

Conclusions

The obtained results demonstrated the efficiency of exome capture for genotyping species such as loblolly pine with a large and complex genome. The highly diverse genetic variation reported in this study will be a valuable resource for future genetic and genomic research in loblolly pine.

Electronic supplementary material

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

Decoding the massive genome of loblolly pine using haploid DNA and novel assembly strategies

Background

The size and complexity of conifer genomes has, until now, prevented full genome sequencing and assembly. The large research community and economic importance of loblolly pine, Pinus taeda L., made it an early candidate for reference sequence determination.

Results

We develop a novel strategy to sequence the genome of loblolly pine that combines unique aspects of pine reproductive biology and genome assembly methodology. We use a whole genome shotgun approach relying primarily on next generation sequence generated from a single haploid seed megagametophyte from a loblolly pine tree, 20-1010, that has been used in industrial forest tree breeding. The resulting sequence and assembly was used to generate a draft genome spanning 23.2 Gbp and containing 20.1 Gbp with an N50 scaffold size of 66.9 kbp, making it a significant improvement over available conifer genomes. The long scaffold lengths allow the annotation of 50,172 gene models with intron lengths averaging over 2.7 kbp and sometimes exceeding 100 kbp in length. Analysis of orthologous gene sets identifies gene families that may be unique to conifers. We further characterize and expand the existing repeat library based on the de novo analysis of the repetitive content, estimated to encompass 82% of the genome.

Conclusions

In addition to its value as a resource for researchers and breeders, the loblolly pine genome sequence and assembly reported here demonstrates a novel approach to sequencing the large and complex genomes of this important group of plants that can now be widely applied.

Association of loblolly pine xylem development gene expression with single-nucleotide polymorphisms

Variation in the expression of genes with putative roles in wood development was associated with single-nucleotide polymorphisms (SNPs) using a population of loblolly pine (Pinus taeda L.) that included individuals from much of the native range. Association studies were performed using 3938 SNPs and expression data obtained using quantitative real-time polymerase chain reaction (PCR) (qRT-PCR) for 106 xylem development genes in 400 clonally replicated loblolly pine individuals. A general linear model (GLM) approach, which takes the underlying population structure into consideration, was used to discover significant associations. After adjustment for multiple testing using a false discovery rate correction, 88 statistically significant associations (Q<0.05) were observed for 80 SNPs with the expression data of 33 xylem development genes. Thirty SNPs caused nonsynonymous mutations, 18 resulted in synonymous mutations, 11 were in 3' untranslated regions (UTRs), 1 was in a 5' UTR and 20 were in introns. Using AraNet, we found that Arabidopsis genes with high similarity to the loblolly pine genes involved in 21 of the 88 statistically significant associations are connected in functional gene networks. Comparisons of gene expression values revealed that in most cases the average expression in plants homozygous for the rare SNP allele was lower than that of plants that were heterozygous or homozygous for the abundant allele. Although there are association studies of SNPs and expression profiles for humans, Arabidopsis and white spruce, to the best of our knowledge, this is the first example of such an association genetic study in pines. Functional validation of these associations will lead to a deeper understanding of the molecular basis of phenotypic differences in wood development among individuals in conifer populations.

The Pinus taeda genome is characterized by diverse and highly diverged repetitive sequences

BACKGROUND

In today's age of genomic discovery, no attempt has been made to comprehensively sequence a gymnosperm genome. The largest genus in the coniferous family Pinaceae is Pinus, whose 110-120 species have extremely large genomes (c. 20-40 Gb, 2N = 24). The size and complexity of these genomes have prompted much speculation as to the feasibility of completing a conifer genome sequence. Conifer genomes are reputed to be highly repetitive, but there is little information available on the nature and identity of repetitive units in gymnosperms. The pines have extensive genetic resources, with approximately 329000 ESTs from eleven species and genetic maps in eight species, including a dense genetic map of the twelve linkage groups in Pinus taeda.

RESULTS

We present here the Sanger sequence and annotation of ten P. taeda BAC clones and Genome Analyzer II whole genome shotgun (WGS) sequences representing 7.5% of the genome. Computational annotation of ten BACs predicts three putative protein-coding genes and at least fifteen likely pseudogenes in nearly one megabase of sequence. We found three conifer-specific LTR retroelements in the BACs, and tentatively identified at least 15 others based on evidence from the distantly related angiosperms. Alignment of WGS sequences to the BACs indicates that 80% of BAC sequences have similar copies (> or = 75% nucleotide identity) elsewhere in the genome, but only 23% have identical copies (99% identity). The three most common repetitive elements in the genome were identified and, when combined, represent less than 5% of the genome.

CONCLUSIONS

This study indicates that the majority of repeats in the P. taeda genome are 'novel' and will therefore require additional BAC or genomic sequencing for accurate characterization. The pine genome contains a very large number of diverged and probably defunct repetitive elements. This study also provides new evidence that sequencing a pine genome using a WGS approach is a feasible goal.

Seasonal variation in gene expression for loblolly pines (Pinus taeda) from different geographical regions

In developing xylem, gene expression levels vary in different genotypes, at different stages of development, throughout a growing season, and in response to stresses. Commercially important characteristics such as wood-specific gravity are known to differ with seed source. For example, when grown on a common site, the specific gravity of Arkansas loblolly pine (Pinus taeda L.) trees is greater than that of Louisiana loblolly pine, and Texas loblolly pines have a greater specific gravity than loblolly pines from the Atlantic coast. A microarray analysis was performed to examine variation in gene expression among trees from different geographical sources when grown on a common site, and seasonal variation in gene expression in each seed source. We used microarrays containing 2171 expressed sequence tags (ESTs) with putative functions of interest, selected from several loblolly pine xylem partial cDNA libraries and a shoot tip library. Genes with significant variation in expression for each factor were identified. Many genes preferentially expressed in latewood compared with earlywood were for proteins involved in cell wall biosynthesis. Variation in gene expression among trees from the two seed sources in each growing season suggests that there may be more differences between South Arkansas trees and South Louisiana trees in latewood than in earlywood. Variation in gene expression among trees from different regions may reflect adaptation to different environments.

An arabinogalactan protein associated with secondary cell wall formation in differentiating xylem of loblolly pine

Arabinogalactan proteins (AGPs) are abundant plant proteoglycans implicated in plant growth and development. Here, we report the genetic characterization, partial purification and immunolocalization of a classical AGP (PtaAGP6, accession number AF101785) in loblolly pine (Pinus taeda L.). A PtaAGP6 full-length cDNA clone was expressed in bacteria. PtaAGP6 resembles tomato LeAGP-1 and Arabidopsis AtAGP17-19 in that they all possess a subdomain composed of basic amino acids. The accessibility of this domain in the glycoprotein makes it possible to label the PtaAGP6 epitopes on the cell surface or in the cell wall with polyclonal antibodies raised against this subdomain. The antibodies recognize the peptide of the basic subdomain and bind to the intact protein molecule. A soluble protein-containing fraction was purified from the differentiating xylem of pine trees by using beta-glucosyl Yariv reagent (beta-glcY) and was recognized by antibodies against the basic subdomain. Immunolocalization studies showed that the PtaAGP6 epitopes are restricted to a file of cells that just precede secondary cell wall thickening, suggesting roles in xylem differentiation and wood formation. The location of apparent labeling of the PtaAGP6 epitopes is separated from the location of lignin deposition. Multiple single nucleotide polymorphisms (SNPs) were detected in EST variants. Denaturing HPLC analysis of PCR products suggests that PtaAGP6 is encoded by a single gene. Mobility variation in denaturing gel electrophoresis was used to map PtaAGP6 SNPs to a site on linkage group 5.

Purification and cloning of an arabinogalactan-protein from xylem of loblolly pine

An arabinogalactan-protein (AGP) was purified from differentiating xylem of loblolly pine (Pinus taeda L.) and the N-terminal sequence used to identify a cDNA clone. The protein, PtaAGP3, was not coded for by any previously identified AGP-like genes. Moreover, PtaAGP3 was abundantly and preferentially expressed in differentiating xylem. The encoded protein contains four domains, a signal peptide, a cleaved hydrophilic region, a region rich in serine, alanine, and proline/hydroxyproline, and a hydrophobic C-terminus. It is postulated to contain a GPI (glycosylphosphatidylinositol) anchor site. If the protein is cleaved at the putative GPI anchor site, as has been observed in other classical AGPs, all but the Ser-Ala-Pro/Hyp-rich domain may be missing from the mature protein. Xylem-specific AGPs are hypothesized to be involved in xylem development.

Two pine endo-beta-1,4-glucanases are associated with rapidly growing reproductive structures

Two cDNA clones encoding endo-beta-1,4-glucanases (EGases) were isolated from a radiata pine (Pinus radiata) cDNA library prepared from immature female strobili. The cDNAs PrCel1 (inus adiata cellulase ) and PrCel2 encode proteins 509 and 515 amino acids in length, respectively, including putative signal peptides. Both proteins contain domains conserved in plant and bacterial EGases. The proteins PRCEL1 and PRCEL2 showed strong similarity to each other (76% amino acid identity), and higher similarity to TPP18 (73 and 67%, respectively), an EGase cloned from tomato (Lycopersicon esculentum) pistils, than to any other reported EGases. Northern-blot analyses indicated that both genes displayed a similar pattern of expression. The only significant difference was in the level of expression. In situ hybridizations were used to demonstrate that, within differentiating pine reproductive structures, PrCel1 expression was greatest in microsporangia in pollen strobili and near the developing ovule in the seed strobili. Expression was also found in vegetative tissues, especially in regions experiencing cell elongation, such as the elongating region of root tips. Both proteins have an ability to degrade carboxymethylcellulose in vitro. Genomic-blot analysis indicated the presence of a family of EGase genes in the radiata pine genome, and that PrCel1 and PrCel2 are transcribed from distinct one-copy genes.

Xylem-specific gene expression in loblolly pine

Two genes preferentially expressed in differentiating xylem of loblolly pine (Pinus taeda L.) were cloned from cDNA and genomic libraries and designated PtX3H6 and PtX14A9. Transcripts of PtX3H6 and PtX14A9 are very abundant in differentiating xylem, less abundant in needles, and very low or non-detectable in embryos and megagametophytes. PtX3H6 contains a putative signal peptide, a threonine-rich region, a proline-rich region, and a hydrophobic tail. Repeats of Pro-Pro-Pro-Val-X-X are similar to repeats found in proline-rich cell wall proteins. The amino acid compositions of PtX3H6 and PtX14A9 are similar to those of arabinogalactan proteins (AGPs). PtX14A9 contains an 8 amino acid sequence similar to amino terminal sequences of ryegrass, carrot and rose AGPs. Upstream sequences have been determined from genomic clones encoding PtX3H6 and PtX14A9. A 7 bp sequence found in the 5' flanking regions of both genes has previously been shown to be involved in the vascular-specific expression of GRP 1.8, a glycine-rich protein found in bean. The sequence is also present upstream of another glycine-rich protein from bean, GRP 1.0, and may be partially responsible for the xylem-specific expression of pTx3H6 and PtX14A9.

Agrobacterium-mediated DNA transfer in sugar pine

DNA transfer using Agrobacterium tumefaciens has been demonstrated in sugar pine, Pinus lambertiana Dougl. Shoots derived from cytokinin-treated cotyledons formed galls after inoculation with A. tumefaciens strains containing the plasmid pTiBo542. A selectable marker, neomycin phosphotransferase II, conferring resistance to kanamycin, was transferred into sugar pine using a binary armed vector system. Callus proliferated from the galls grew without hormones and in some cases, kanamycin-resistant callus could be cultured. Southern blots provided evidence of physical transfer of T-DNA and the nptII gene. Expression of the nptII gene under control of the nos promoter was demonstrated by neomycin phosphotransferase assays. Several aspects of DNA transfer were similar to those previously observed in angiosperms transformed by A. tumefaciens. This is the first evidence for DNA transfer by Agrobacterium in this species and the first physical evidence for transfer in any pine. These results bring us closer to genetic engineering in this commercially important genus of forest trees.