Bidirectional introgression between loblolly pine (Pinus taeda L.) and shortleaf pine (P. echinata Mill.) has increased since the 1950s

Hybridization between loblolly pine and shortleaf pine has been shown in sympatric and allopatric portions of the species native ranges. In this study, we used 25 microsatellite markers to determine the hybrid status of 165 and 151 loblolly pine and shortleaf pine trees sampled across their ranges, respectively, and to estimate population differentiation within these species. Estimated differentiation (?? _PT) for these current-day samples was significantly higher in both species??0.115 in loblolly pine and 0.146 in shortleaf pine??than for trees planted from seed collected from the same locations in the 1950s. These increases are likely due to anthropogenic causes such as habitat fragmentation. In addition, the proportion of hybrids rose dramatically in both species: 27.3% hybrids in loblolly pine populations and 46.7% hybrids in shortleaf pine populations compared to rates of 4.5% and 3.3%, respectively, in the 1950s populations. Our results suggest that shortleaf pine and remnant naturally regenerated loblolly pine are at risk.     

Natural hybridization within seed sources of shortleaf pine (Pinus echinata Mill.) and loblolly pine (Pinus taeda L.)

Shortleaf and loblolly pine trees (n = 93 and 102, respectively) from 22 seed sources of the Southwide Southern Pine Seed Source Study plantings or equivalent origin were evaluated for amplified fragment length polymorphism (AFLP) variation. These sampled trees represent shortleaf pine and loblolly pine, as they existed across their native geographic ranges before intensive forest management. Using 17 primer pairs, a total of 96 AFLPs between shortleaf pine and loblolly pine were produced and scored on the sample trees and two control-pollinated F1 interspecies hybrids and their parents. In addition, the well known isocitrate dehydrogenase (IDH) isozyme marker was scored for all trees. IDH detected two putative hybrids among the loblolly pine samples and two among the shortleaf pine samples, while either 13 or 12 putative hybrids were detected using all AFLP markers and IDH and either NewHybrids or Structure software, respectively. Results of this study show that later generation hybrids can be reliably identified using AFLP markers and confirmed that IDH is not a definitive marker for detecting hybrids; that is, at least in some seed sources, the alternative species’ IDH allele resides in the source species. Based on all the markers, hybridization frequency varied geographically, ranging from 30% in an Arkansas seed source to 0% in several other seed sources. The hybridization level was higher in populations west of the Mississippi River than in populations east of the river; the shortleaf pine hybridization rates were 16.3% and 2.4% and the loblolly pine rates were 4.5% and 3.3%, west and east of the river, respectively. The results suggest that hybridization between these two species is significant but varies by seed source and species, and the potential for the unintended creation of hybrids should be considered in forest management decisions regarding both natural and artificial regeneration.     

Genetic diversity within and among populations of shortleaf pine (Pinus echinata Mill.) and loblolly pine (Pinus taeda L.)

Shortleaf pine (n = 93) and loblolly pine (n = 112) trees representing 22 seed sources or 16 physiographic populations were sampled from Southwide Southern Pine Seed Source Study plantings located in Oklahoma, Arkansas, and Mississippi. The sampled trees were grown from shortleaf pine and loblolly pine seeds formed in 1951 and 1952, prior to the start of intensive forest management across their native ranges. Amplification fragment length polymorphism (AFLP) markers were developed and used to study genetic diversity and its structure in these pine species. After screening 48 primer pairs, 17 and 21 pairs were selected that produced 794 and 647 AFLPs in shortleaf pine and loblolly pine, respectively. High-AFLP-based genetic diversity exists within shortleaf pine and loblolly pine, and most (84.73% in shortleaf pine; 87.69% in loblolly pine) of this diversity is maintained within physiographic populations. The high value of unbiased measures of genetic identity and low value of genetic distance for all pairwise comparisons indicates that the populations have similar genetic structures. For shortleaf pine, there was no significant correlation between geographic distance and genetic distance (r = 0.28), while for loblolly pine there was a weak but significant correlation (r = 0.51).     

Pollen pool heterogeneity in shortleaf pine, Pinus echinata Mill

Pollen is the dominant vector of gamete exchange for most temperate tree species. Because pollen movement influences the creation, maintenance and erosion of genetic structure in adult populations, it is important to understand what factors influence the process of pollen movement. Isolation by distance in pollen donor populations can create highly structured pollen polls by increased sampling of local fathers. Extrinsic factors, such as the intervening vegetative structure and local pollen donor densities, can also influence the genetic composition of local pollen pools. Using paternally inherited chloroplast microsatellite markers, we examined the structure and diversity of pollen pools in Pinus echinata Mill. in southern Missouri, USA. Our analysis is based on a multivariate AMOVA analysis of stands ( approximately 1 ha; six per region) nested within regions (approximately 800 ha; four each). Significant multilocus structure of the pollen pool within regions (phiSR = 0.095), but not among regions (phiRT = 0.010), indicates that pollen movement is relatively restricted. Furthermore, the significant correlation between pairwise genetic and physical distances (Mantel correlation; rho = 0.32) provided support for the isolation by distance hypothesis. Our results indicated that availability of pollen donors did not affect diversity of the pollen pool, measured by the number of unique multilocus genotypes at each stand. However, pollen pool diversity was negatively associated with vegetative structure, measured as total forest tree density. Our findings indicated that on-going pollen movement within continuous forest is relatively restricted as a result of both isolation by distance and vegetative structure.     

Association genetics of the loblolly pine (Pinus taeda, Pinaceae) metabolome

The metabolome of a plant comprises all small molecule metabolites, which are produced during cellular processes. The genetic basis for metabolites in nonmodel plants is unknown, despite frequently observed correlations between metabolite concentrations and stress responses. A quantitative genetic analysis of metabolites in a nonmodel plant species is thus warranted. Here, we use standard association genetic methods to correlate 3563 single nucleotide polymorphisms (SNPs) to concentrations of 292 metabolites measured in a single loblolly pine (Pinus taeda) association population. A total of 28 single locus associations were detected, representing 24 and 20 unique SNPs and metabolites, respectively. Multilocus Bayesian mixed linear models identified 2998 additional associations for a total of 1617 unique SNPs associated to 255 metabolites. These SNPs explained sizeable fractions of metabolite heritabilities when considered jointly (56.6% on average) and had lower minor allele frequencies and magnitudes of population structure as compared with random SNPs. Modest sets of SNPs (n = 1-23) explained sizeable portions of genetic effects for many metabolites, thus highlighting the importance of multi-SNP models to association mapping, and exhibited patterns of polymorphism consistent with being linked to targets of natural selection. The implications for association mapping in forest trees are discussed.     

Segregation and linkage relationships of isoenzymes in shortleaf pine (Pinus echinata Mill.)

 Segregation and linkage relationships were analyzed between 28 isoenzyme loci in ten natural stands representing much of the natural range of Pinus echinata Mill. (shortleaf pine). A total of 203 possible two-locus combinations were tested. Three linkage groups were revealed in this study at a linkLOD of 4.0. The first linkage group (A) consisted of Pgi and Adh-1; Gdh, Idh, Skdh-2, G6pd-2 and Aco were mapped to the second linkage group (B); the third group (C) had 2 loci: Mdh-2 and Mdh-3. A moderate linkage between Mnr-2 and Dia-2 and weak linkages between Mnr-1 and Dia-1, and Got-2 and 6pgd-2 were also detected. The significance of these results in shortleaf pine is discussed and compared to linkage maps previously reported in other conifers, including pines. Received: 7 April 1997 / Accepted: 25 April 1997     

Development of DNA methylation-based epigenetic age predictors in loblolly pine (Pinus taeda)

Biological ageing is connected to life history variation across ecological scales and informs a basic understanding of age-related declines in organismal function. Altered DNA methylation dynamics are a conserved aspect of biological ageing and have recently been modelled to predict chronological age among vertebrate species. In addition to their utility in estimating individual age, differences between chronological and predicted ages arise due to acceleration or deceleration of epigenetic ageing, and these discrepancies are linked to disease risk and multiple life history traits. Although evidence suggests that patterns of DNA methylation can describe ageing in plants, predictions with epigenetic clocks have yet to be performed. Here, we resolve the DNA methylome across CpG, CHG, and CHH-methylation contexts in the loblolly pine tree (Pinus taeda) and construct epigenetic clocks capable of predicting ages in this species within 6% of its maximum lifespan. Although patterns of CHH-methylation showed little association with age, both CpG and CHG-methylation contexts were strongly associated with ageing, largely becoming hypomethylated with age. Among age-associated loci were those in close proximity to malate dehydrogenase, NADH dehydrogenase, and 18S and 26S ribosomal RNA genes. This study reports one of the first epigenetic clocks in plants and demonstrates the universality of age-associated DNA methylation dynamics which can inform conservation and management practices, as well as our ecological and evolutionary understanding of biological ageing in plants.     

Cold stratification and growth of radicles of loblolly pine (Pinus taeda) embryos

The mechanism of germination enhancement by cold stratification was examined in seeds of loblolly pine ( Pinus taeda L. ), Removal of the seed coat permitted elongation of radicles from unstratified embryos, but both rates of germination and radicle elongation were increased by stratification. Radicles of both stratified and unstratified embryos excised from the megagamethophyte elongated only when in contact with solid incubation media supplemented with sucrose. Stratification of embryos either in the presence or absence of the megagametophyte resulted in similar enhancement of radicle elongation. Elongation rates of radicles were increased after stratification independent of sucrose concentration, and changes in sucrose content in the megagamethophyte during stratification or incubations subsequent to stratification were insufficient to regulate radicle growth. Our results support the hypothesis that the embryos of pine seeds perceive the low temperature stimulus directly and this stimulus results in a growth potential increase in the embryonic axes. We propose that this growth potential increase enables the embryos to overcome the mechanical restraint of the seeds coats and to germinate.     

Monoterpene synthases of loblolly pine (Pinus taeda) produce pinene isomers and enantiomers

The turpentine fraction of conifer oleoresin is a complex mixture of monoterpene olefins and plays important roles in defense and in the mediation of chemical communication between conifer hosts and insect predators. The stereochemistry of the turpentine monoterpenes is critical in these interactions, influencing host recognition, toxicity, and potency of derived pheromones, and the stereochemical composition of these compounds lends insight into their biogenetic origin, with implications for the numbers and types of enzymes responsible and their corresponding genes. Analysis of the oleoresin from several tissues of loblolly pine (Pinus taeda) showed the derived turpentine to consist mainly of (+)-(3R:5R)-alpha-pinene and (-)-(3S:5S)-beta-pinene. Cell-free extracts from xylem tissue yielded three monoterpene synthases which together account for the monoterpene isomer and enantiomer content of the turpentine of this tissue. The major products of these enzymes, produced from the universal precursor of monoterpenes, geranyl diphosphate, were shown to be (+)-alpha-pinene, (-)-alpha-pinene, and (-)-beta-pinene, respectively. In most properties (molecular mass of approximately 60 kDa, K(m) for geranyl diphosphate of 3 microM, requirement for monovalent and divalent cations), these enzymes resemble other monoterpene synthases from conifer species.     

Comparative mapping in Pinus: sugar pine (Pinus lambertiana Dougl.) and loblolly pine (Pinus taeda L.)

The majority of genomic research in conifers has been conducted in the Pinus subgenus Pinus mostly due to the high economic importance of the species within this taxon. Genetic maps have been constructed for several of these pines and comparative mapping analyses have consistently revealed notable synteny. In contrast, little genomic research has been conducted on the Pinus subgenus Strobus, even though these pines have strong ecological relevance. We report a consensus genetic linkage map for sugar pine (Pinus lambertiana Dougl.) constructed with 399 single nucleotide polymorphisms markers derived from annotated genes. The map is 1,231 cM in length and organized into 19 linkage groups. Two of the mapping populations were derived from trees that were segregating for the major gene of resistance (Cr1) to Cronartium ribicola, the fungal pathogen responsible for white pine blister rust. The third mapping population was derived from a full-sib cross segregating for partial resistance to white pine blister rust. In addition, we report the first comparative mapping study between subgenera Strobus and Pinus. Sixty mapped markers were found in common between sugar pine and the loblolly pine reference map with 56 of them (93%) showing collinearity. All 19 linkage groups of the sugar pine consensus map coaligned to the 12 linkage groups of the loblolly pine reference map. The syntenic relationship observed between these two clades of pines provides a foundation for advancing genomic research and genetic resources in subgenus Strobus.     

Induction of somatic embryogenesis in loblolly pine (Pinus taeda L.)

Summary Several factors that may affect induction of somatic embryogenesis in loblolly pine (Pinus taeda L.) were investigated in 1994 and 1995. Megagametophytes containing immature zygotic embryos were excised from seeds as explants. Potassium chloride, silver nitrate, myo-inositol, coconut water, or polyamine was added to the control media (U.S. patent no. 5,036,007) to determine the effects of each single ingredient or their combinations on the initiation of embryogenic tissue. Supplements of myo-inositol at 22.2 mM resulted in increases in frequencies of cell mass extrusion and proliferation compared with the control media in consecutive years. Addition of silver nitrate showed the potential to promote initiation of embryogenic culture. The combination of 10 mM potassium with 29.4 μM silver nitrate achieved the highest frequencies in both extrusion and proliferation of embryogenic tissue. The combination of silver nitrate at 29.4 μM with addition of myo-inositol at 11.1 or 22.2 mM achieved a higher conversion rate from extrusion to proliferation. Polyamine did not significantly affect the induction of somatic embryogenesis, but coconut water was inhibitory. Published with approval of the Director of Arkansas Agricultural Experimental Station.     

Intergenotypic interactions among families of loblolly pine (Pinus taeda L.)

Summary The effects of competition on the growth of families of loblolly pine (Pinus taeda, L.) seedlings were investigated. The experimental design made it possible to evaluate the effects of crowding on growth and to determine the types and magnitudes of intergenotypic interactions among pairs of families. The results showed that intergenotypic interactions were both highly variable and pronounced in their effect on early growth. Evidence was also found for precompetition cooperating interactions occurring among seedlings surrounded by neighbors of the same family.     

Family 34 glycosyltransferase (GT34) genes and proteins in Pinus radiata (radiata pine) and Pinus taeda (loblolly pine)

Using a functional genomics approach, four candidate genes (PtGT34A, PtGT34B, PtGT34C and PtGT34D) were identified in Pinus taeda. These genes encode CAZy family GT34 glycosyltransferases that are involved in the synthesis of cell‐wall xyloglucans and heteromannans. The full‐length coding sequences of three orthologs (PrGT34A, B and C) were isolated from a xylem‐specific cDNA library from the closely related Pinus radiata. PrGT34B is the ortholog of XXT1 and XXT2, the two main xyloglucan (1→6)‐α‐xylosyltransferases in Arabidopsis thaliana. PrGT34C is the ortholog of XXT5 in A. thaliana, which is also involved in the xylosylation of xyloglucans. PrGT34A is an ortholog of a galactosyltransferase from fenugreek (Trigonella foenum‐graecum) that is involved in galactomannan synthesis. Truncated coding sequences of the genes were cloned into plasmid vectors and expressed in a Sf9 insect cell‐culture system. The heterologous proteins were purified, and in vitro assays showed that, when incubated with UDP‐xylose and cellotetraose, cellopentaose or cellohexaose, PrGT34B showed xylosyltransferase activity, and, when incubated with UDP‐galactose and the same cello‐oligosaccharides, PrGT34B showed some galactosyltransferase activity. The ratio of xylosyltransferase to galactosyltransferase activity was 434:1. Hydrolysis of the galactosyltransferase reaction products using galactosidases showed the linkages formed were α‐linkages. Analysis of the products of PrGT34B by MALDI‐TOF MS showed that up to three xylosyl residues were transferred from UDP‐xylose to cellohexaose. The heterologous proteins PrGT34A and PrGT34C showed no detectable enzymatic activity.     

Anchored reference loci in loblolly pine (Pinus taeda L.) for integrating pine genomics

Anchored reference loci provide a framework for comparative mapping. They are landmarks to denote conserved chromosomal segments, allowing the synthesis of genetic maps from multiple sources. We evaluated 90 expressed sequence tag polymorphisms (ESTPs) from loblolly pine (Pinus taeda L.) for this function. Primer sets were assayed for amplification and polymorphism in six pedigrees, representing two subgenera of Pinus and a distant member of the Pinaceae, Douglas-fir (Pseudotsuga menziesii [Mirb.] Franco). On average, 89% of primer sets amplified in four species of subgenus Pinus, 49% in one species of subgenus Strobus, and 22% in Douglas-fir. Polymorphisms were detected for 37-61% of the ESTPs within each pedigree. Comparative mapping in loblolly and slash pine (P. elliottii Englm.) revealed that ESTPs mapped to the same location. Disrupted synteny or significant disruptions in colinearity were not detected. Thirty-five ESTPs met criteria established for anchor loci. The majority of those that did not meet these criteria were excluded when map location was known in only a single species. Anchor loci provide a unifying tool for the community, facilitating the creation of a "generic" pine map and serving as a foundation for studies on genome organization and evolution.     

High-efficiency Agrobacterium-mediated transformation of Norway spruce (Picea abies) and loblolly pine (Pinus taeda)

Agrobacterium-mediated gene transfer is the method of choice for many plant biotechnology laboratories; however, large- scale use of this organism in conifer transformation has been limited by difficult propagation of explant material, selection efficiencies and low transformation frequency. We have analyzed co-cultivation conditions and different disarmed strains of Agrobacterium to improve transformation. Additional copies of virulence genes were added to three common disarmed strains. These extra virulence genes included either a constitutively active virG or extra copies of virG and virB, both from pTiBo542. In experiments with Norway spruce, we increased transformation efficiencies 1000-fold from initial experiments where little or no transient expression was detected. Over 100 transformed lines expressing the marker gene -glucuronidase (GUS) were generated from rapidly dividing embryogenic suspension-cultured cells co- cultivated with Agrobacterium. GUS activity was used to monitor transient expression and to further test lines selected on kanamycin-containing medium. In loblolly pine, transient expression increased 10-fold utilizing modified Agrobacterium strains. Agrobacterium-mediated gene transfer is a useful technique for large-scale generation of transgenic Norway spruce and may prove useful for other conifer species.     

Anatomy of seedling tap roots of loblolly pine (Pinus taeda L.)

The development of tap root anatomical features was investigated in seedlings of loblolly pine (Pinus taeda L.) under both pot and pouch growth regimes. The roots possessed the three anatomical zones previously observed in jack pine (Pinus banksiana Lamb) and Eucalyptus pilularis Sm. - white, condensed tannin (CT), and cork - suggesting that this developmental sequence is preserved over species and growth conditions. Xylem development was centripetal and similar to that found earlier in P. sylvestris. Tracheids with lignified, secondary walls were detected distal to the point of endodermal Casparian band deposition. However, tests for ability to conduct fluid indicated that the protoxylem was capable of transport only proximal to the Casparian bands. Detailed examination of suberin lamella deposition in the endodermis demonstrated that passage cells were present through the white and CT zones. Progressive, centripetal cortical death in the CT zone did not include the endodermis, which remained alive until the cork layer formed, at which point the endodermis was crushed. Therefore, passage cells remain as functional portals for nutrient and water uptake in the CT zone even though the central cortex is dead. Tracer tests indicated that the endodermis provides an apoplastic barrier to tracer diffusion into the stele and that this function was taken over by the young cork layers. Results of this study point to a strong role for the endodermis in the regulation of nutrient and water uptake until the maturation of the first cork layer.     

Water status and growth of loblolly pine (Pinus taeda L.) callus

Water status of Pinus taeda L. callus supported on Murashige and Skoog (MS) liquid medium was characterized over an 8 week period using thermocouple psychrometry. Medium with 30 gl⁻¹ sucrose was used to produce a high water potential (Ψ_w) of −0.4 MPa (H), and the same medium was used to create a moderate Ψ_w of −0.7 MPa (M) by the addition of 10% polyethylene glycol (PEG, w/v, MW=8000). Calli were produced from cotyledon explants on H medium for 2 weeks and then transferred to either M or H medium. Callus absorption of PEG accounted for 40% of the callus dry weight and less than 7% of the callus fresh weight. Callus dry weight (without the PEG fraction) on M medium was 40% of that observed on H medium. Fresh weight on M medium was only 15% of that observed on H medium. The Ψ_w of both H and M media remained constant throughout the culture period. On H medium, callus Ψ_w and osmotic potential (Ψ_s) both increased 0.05 MPa/week with the callus Ψ_w approaching that of the external medium. On M medium, callus Ψ_w and Ψ_s both decreased more than 0.1 MPa/week with the callus Ψ_w decreasing greatly below that of the external medium. The latter was attributed to a rapidly produced osmotic shock induced upon callus transfer and/or PEG which caused less callus hydration and resulted in reduced growth. Callus turgor potential (Ψ_p) was estimated to be +0.02 to +0.09 MPa and turgor was maintained as callus Ψ_w increased or decreased. After 8 weeks, cell volumes from callus on M medium were 50 to 60% less than on H medium, suggesting that reduced cell volumes were related to turgor maintenance.