QTL associations for density and diameter in Pinus radiata and the potential for marker-aided selection

High density linkage maps, genetic architecture, and genomic prediction of growth and wood properties in Pinus radiata

Background

The growing availability of genomic resources in radiata pine paves the way for significant advances in fundamental and applied genomic research. We constructed robust high-density linkage maps based on exome-capture genotyping in two F₁ populations, and used these populations to perform quantitative trait locus (QTL) scans, genomic prediction and quantitative analyses of genetic architecture for key traits targeted by tree improvement programmes.

Results

Our mapping approach used probabilistic error correction of the marker data, followed by an iterative approach based on stringent parameters. This approach proved highly effective in producing high-density maps with robust marker orders and realistic map lengths (1285–4674 markers per map, with sizes ranging from c. 1643–2292 cM, and mean marker intervals of 0.7–2.1 cM). Colinearity was high between parental linkage maps, although there was evidence for a large chromosomal rearrangement (affecting ~ 90 cM) in one of the parental maps. In total, 28 QTL were detected for growth (stem diameter) and wood properties (wood density and fibre properties measured by Silviscan) in the QTL discovery population, with 1–3 QTL of small to moderate effect size detected per trait in each parental map. Four of these QTL were validated in a second, unrelated F₁ population. Results from genomic prediction and analyses of genetic architecture were consistent with those from QTL scans, with wood properties generally having moderate to high genomic heritabilities and predictive abilities, as well as somewhat less complex genetic architectures, compared to growth traits.

Conclusions

Despite the economic importance of radiata pine as a plantation forest tree, robust high-density linkage maps constructed from reproducible, sequence-anchored markers have not been published to date. The maps produced in this study will be a valuable resource for several applications, including the selection of marker panels for genomic prediction and anchoring a recently completed de novo whole genome assembly. We also provide the first map-based evidence for a large genomic rearrangement in radiata pine. Finally, results from our QTL scans, genomic prediction, and genetic architecture analyses are informative about the genomic basis of variation in important phenotypic traits.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-022-08950-6.

A high-density genetic map of Schima superba based on its chromosomal characteristics

BACKGROUND

Schima superba (Theaceae) is a popular woody tree in China. The obscure chromosomal characters of this species are a limitation in the development of high-density genetic linkage maps, which are valuable resources for molecular breeding and functional genomics.

RESULTS

We determined the chromosome number and the karyotype of S. superba as 2n = 36 = 36 m, which is consistent with the tribe Schimeae (n = 18). A high-density genetic map was constructed using genotyping by sequencing (GBS). A F1 full-sib with 116 individuals and their parents (LC31 × JO32) were sequenced on the Illumina HiSeq™ platform. Overall, 343.3 Gb of raw data containing 1,191,933,474 paired-end reads were generated. Based on this, 99,966 polymorphic SNP markers were developed from the parents, and 2209 markers were mapped onto the integrated genetic linkage map after data filtering and SNP genotyping. The map spanned 2076.24 cM and was distributed among 18 linkage groups. The average marker interval was 0.94 cM. A total of 168 quantitative trait loci (QTLs) for 14 growth traits were identified.

CONCLUSIONS

The chromosome number and karyotype of S. superba was 2n = 36 = 36 m and a linkage map with 2209 SNP markers was constructed to identify QTLs for growth traits. Our study provides a basis for molecular-assisted breeding and genomic studies, which will contribute towards the future research and genetic improvement of S. superba.

Linkage and Association Mapping for Two Major Traits Used in the Maritime Pine Breeding Program: Height Growth and Stem Straightness

Background

Increasing our understanding of the genetic architecture of complex traits, through analyses of genotype-phenotype associations and of the genes/polymorphisms accounting for trait variation, is crucial, to improve the integration of molecular markers into forest tree breeding. In this study, two full-sib families and one breeding population of maritime pine were used to identify quantitative trait loci (QTLs) for height growth and stem straightness, through linkage analysis (LA) and linkage disequilibrium (LD) mapping approaches.

Results

The populations used for LA consisted of two unrelated three-generation full-sib families (n = 197 and n = 477). These populations were assessed for height growth or stem straightness and genotyped for 248 and 217 markers, respectively. The population used for LD mapping consisted of 661 founders of the first and second generations of the breeding program. This population was phenotyped for the same traits and genotyped for 2,498 single-nucleotide polymorphism (SNP) markers corresponding to 1,652 gene loci. The gene-based reference genetic map of maritime pine was used to localize and compare the QTLs detected by the two approaches, for both traits. LA identified three QTLs for stem straightness and two QTLs for height growth. The LD study yielded seven significant associations (P ≤ 0.001): four for stem straightness and three for height growth. No colocalisation was found between QTLs identified by LA and SNPs detected by LD mapping for the same trait.

Conclusions

This study provides the first comparison of LA and LD mapping approaches in maritime pine, highlighting the complementary nature of these two approaches for deciphering the genetic architecture of two mandatory traits of the breeding program.

Comparative Genomics Analyses Reveal Extensive Chromosome Colinearity and Novel Quantitative Trait Loci in Eucalyptus

Dense genetic maps, along with quantitative trait loci (QTLs) detected on such maps, are powerful tools for genomics and molecular breeding studies. In the important woody genus Eucalyptus, the recent release of E. grandis genome sequence allows for sequence-based genomic comparison and searching for positional candidate genes within QTL regions. Here, dense genetic maps were constructed for E. urophylla and E. tereticornis using genomic simple sequence repeats (SSR), expressed sequence tag (EST) derived SSR, EST-derived cleaved amplified polymorphic sequence (EST-CAPS), and diversity arrays technology (DArT) markers. The E. urophylla and E. tereticornis maps comprised 700 and 585 markers across 11 linkage groups, totaling at 1,208.2 and 1,241.4 cM in length, respectively. Extensive synteny and colinearity were observed as compared to three earlier DArT-based eucalypt maps (two maps with E. grandis × E. urophylla and one map of E. globulus) and with the E. grandis genome sequence. Fifty-three QTLs for growth (10-56 months of age) and wood density (56 months) were identified in 22 discrete regions on both maps, in which only one colocalizaiton was found between growth and wood density. Novel QTLs were revealed as compared with those previously detected on DArT-based maps for similar ages in Eucalyptus. Eleven to 585 positional candidate genes were obained for a 56-month-old QTL through aligning QTL confidence interval with the E. grandis genome. These results will assist in comparative genomics studies, targeted gene characterization, and marker-assisted selection in Eucalyptus and the related taxa.

QTL mapping of growth-related traits in a full-sib family of rubber tree (Hevea brasiliensis) evaluated in a sub-tropical climate

The rubber tree (Hevea spp.), cultivated in equatorial and tropical countries, is the primary plant used in natural rubber production. Due to genetic and physiological constraints, inbred lines of this species are not available. Therefore, alternative approaches are required for the characterization of this species, such as the genetic mapping of full-sib crosses derived from outbred parents. In the present study, an integrated genetic map was obtained for a full-sib cross family with simple sequence repeats (SSRs) and expressed sequence tag (EST-SSR) markers, which can display different segregation patterns. To study the genetic architecture of the traits related to growth in two different conditions (winter and summer), quantitative trait loci (QTL) mapping was also performed using the integrated map. Traits evaluated were height and girth growth, and the statistical model was based in an extension of composite interval mapping. The obtained molecular genetic map has 284 markers distributed among 23 linkage groups with a total length of 2688.8 cM. A total of 18 QTLs for growth traits during the summer and winter seasons were detected. A comparison between the different seasons was also conducted. For height, QTLs detected during the summer season were different from the ones detected during winter season. This type of difference was also observed for girth. Integrated maps are important for genetics studies in outbred species because they represent more accurately the polymorphisms observed in the genitors. QTL mapping revealed several interesting findings, such as a dominance effect and unique segregation patterns that each QTL could exhibit, which were independent of the flanking markers. The QTLs identified in this study, especially those related to phenotypic variation associated with winter could help studies of marker-assisted selection that are particularly important when the objective of a breeding program is to obtain phenotypes that are adapted to sub-optimal regions.

Towards decoding the conifer giga-genome

Several new initiatives have been launched recently to sequence conifer genomes including pines, spruces and Douglas-fir. Owing to the very large genome sizes ranging from 18 to 35 gigabases, sequencing even a single conifer genome had been considered unattainable until the recent throughput increases and cost reductions afforded by next generation sequencers. The purpose of this review is to describe the context for these new initiatives. A knowledge foundation has been acquired in several conifers of commercial and ecological interest through large-scale cDNA analyses, construction of genetic maps and gene mapping studies aiming to link phenotype and genotype. Exploratory sequencing in pines and spruces have pointed out some of the unique properties of these giga-genomes and suggested strategies that may be needed to extract value from their sequencing. The hope is that recent and pending developments in sequencing technology will contribute to rapidly filling the knowledge vacuum surrounding their structure, contents and evolution. Researchers are also making plans to use comparative analyses that will help to turn the data into a valuable resource for enhancing and protecting the world's conifer forests.

Annotated genetic linkage maps of Pinus pinaster Ait. from a Central Spain population using microsatellite and gene based markers

BACKGROUND

Pinus pinaster Ait. is a major resin producing species in Spain. Genetic linkage mapping can facilitate marker-assisted selection (MAS) through the identification of Quantitative Trait Loci and selection of allelic variants of interest in breeding populations. In this study, we report annotated genetic linkage maps for two individuals (C14 and C15) belonging to a breeding program aiming to increase resin production. We use different types of DNA markers, including last-generation molecular markers.

RESULTS

We obtained 13 and 14 linkage groups for C14 and C15 maps, respectively. A total of 211 and 215 markers were positioned on each map and estimated genome length was between 1,870 and 2,166 cM respectively, which represents near 65% of genome coverage. Comparative mapping with previously developed genetic linkage maps for P. pinaster based on about 60 common markers enabled aligning linkage groups to this reference map. The comparison of our annotated linkage maps and linkage maps reporting QTL information revealed 11 annotated SNPs in candidate genes that co-localized with previously reported QTLs for wood properties and water use efficiency.

CONCLUSIONS

This study provides genetic linkage maps from a Spanish population that shows high levels of genetic divergence with French populations from which segregating progenies have been previously mapped. These genetic maps will be of interest to construct a reliable consensus linkage map for the species. The importance of developing functional genetic linkage maps is highlighted, especially when working with breeding populations for its future application in MAS for traits of interest.

Identification of putative candidate genes for juvenile wood density in Pinus radiata

Wood formation is a complex developmental process driven by the annual activity of the vascular cambium. Conifers usually produce juvenile wood at young ages followed by mature wood for the rest of their lifetime. Juvenile wood exhibits poorer wood quality (i.e., lower density) compared with mature wood and can account for up to 50% of short-rotation harvested logs, thus representing a major challenge for commercial forestry globally. Wood density is an important quality trait for many timber-related products. Understanding the molecular mechanisms involved in the regulation of juvenile wood density is critical for the improvement of juvenile wood quality via marker-aided selection. A previous study has identified several candidate genes affecting mature wood density in Picea sitchensis (Bong.) Carr.; however, genes associated with juvenile wood density in conifers remain poorly characterized. Here, cDNA microarrays containing 3320 xylem unigenes were used to investigate genes differentially transcribed in juvenile wood with high (HD) and low density (LD) in Pinus radiata D.Don. In total, 814 xylem unigenes with differential transcription were identified in at least one of two microarray experiments and 73 genes (45 for HD, 28 for LD) were identified in both experiments, thus representing putative candidate genes for juvenile wood density. Interestingly, cellulose synthases (PrCesA3, PrCesA11) and sucrose synthase (SuSy), which are involved in secondary cell wall formation, had stronger transcription in juvenile wood with HD, while genes functioning in primary wall formation (pectin synthesis, cell expansion and other modifications) were more transcribed in LD wood. Cell wall genes encoding monolignol biosynthesis enzymes, arabinogalactan proteins, actins and tubulins were differentially transcribed in either HD or LD juvenile wood; however, the latter had exclusively greater transcription of genes involved in monolignol polymerization (laccase and peroxidase). The identified candidate genes also included many non-cell-wall genes (transcription factors, environmental-responsive genes, hormone signalling, etc.) and genes with unknown functions, suggesting complex gene pathways in the regulation of juvenile wood density. Interestingly, 19 out of 73 candidate genes for wood density were among the 108 candidate genes previously identified for microfibril angle, and 16 genes appeared to influence both traits in a synergistic manner for wood stiffness.

Biotechnologies for the management of genetic resources for food and agriculture

In recent years, the land area under agriculture has declined as also has the rate of growth in agricultural productivity while the demand for food continues to escalate. The world population now stands at 7 billion and is expected to reach 9 billion in 2045. A broad range of agricultural genetic diversity needs to be available and utilized in order to feed this growing population. Climate change is an added threat to biodiversity that will significantly impact genetic resources for food and agriculture (GRFA) and food production. There is no simple, all-encompassing solution to the challenges of increasing productivity while conserving genetic diversity. Sustainable management of GRFA requires a multipronged approach, and as outlined in the paper, biotechnologies can provide powerful tools for the management of GRFA. These tools vary in complexity from those that are relatively simple to those that are more sophisticated. Further, advances in biotechnologies are occurring at a rapid pace and provide novel opportunities for more effective and efficient management of GRFA. Biotechnology applications must be integrated with ongoing conventional breeding and development programs in order to succeed. Additionally, the generation, adaptation, and adoption of biotechnologies require a consistent level of financial and human resources and appropriate policies need to be in place. These issues were also recognized by Member States at the FAO international technical conference on Agricultural Biotechnologies for Developing Countries (ABDC-10), which took place in March 2010 in Mexico. At the end of the conference, the Member States reached a number of key conclusions, agreeing, inter alia, that developing countries should significantly increase sustained investments in capacity building and the development and use of biotechnologies to maintain the natural resource base; that effective and enabling national biotechnology policies and science-based regulatory frameworks can facilitate the development and appropriate use of biotechnologies in developing countries; and that FAO and other relevant international organizations and donors should significantly increase their efforts to support the strengthening of national capacities in the development and appropriate use of pro-poor agricultural biotechnologies.

Transcriptome profiling of Pinus radiata juvenile wood with contrasting stiffness identifies putative candidate genes involved in microfibril orientation and cell wall mechanics

BACKGROUND

The mechanical properties of wood are largely determined by the orientation of cellulose microfibrils in secondary cell walls. Several genes and their allelic variants have previously been found to affect microfibril angle (MFA) and wood stiffness; however, the molecular mechanisms controlling microfibril orientation and mechanical strength are largely uncharacterised. In the present study, cDNA microarrays were used to compare gene expression in developing xylem with contrasting stiffness and MFA in juvenile Pinus radiata trees in order to gain further insights into the molecular mechanisms underlying microfibril orientation and cell wall mechanics.

RESULTS

Juvenile radiata pine trees with higher stiffness (HS) had lower MFA in the earlywood and latewood of each ring compared to low stiffness (LS) trees. Approximately 3.4 to 14.5% out of 3, 320 xylem unigenes on cDNA microarrays were differentially regulated in juvenile wood with contrasting stiffness and MFA. Greater variation in MFA and stiffness was observed in earlywood compared to latewood, suggesting earlywood contributes most to differences in stiffness; however, 3-4 times more genes were differentially regulated in latewood than in earlywood. A total of 108 xylem unigenes were differentially regulated in juvenile wood with HS and LS in at least two seasons, including 43 unigenes with unknown functions. Many genes involved in cytoskeleton development and secondary wall formation (cellulose and lignin biosynthesis) were preferentially transcribed in wood with HS and low MFA. In contrast, several genes involved in cell division and primary wall synthesis were more abundantly transcribed in LS wood with high MFA.

CONCLUSIONS

Microarray expression profiles in Pinus radiata juvenile wood with contrasting stiffness has shed more light on the transcriptional control of microfibril orientation and the mechanical properties of wood. The identified candidate genes provide an invaluable resource for further gene function and association genetics studies aimed at deepening our understanding of cell wall biomechanics with a view to improving the mechanical properties of wood.

Forest tree genomics: growing resources and applications

Over the past two decades, research in forest tree genomics has lagged behind that of model and agricultural systems. However, genomic research in forest trees is poised to enter into an important and productive phase owing to the advent of next-generation sequencing technologies, the enormous genetic diversity in forest trees and the need to mitigate the effects of climate change. Research on long-lived woody perennials is extending our molecular knowledge of complex life histories and adaptations to the environment - enriching a field that has traditionally drawn biological inference from a few short-lived herbaceous species.

Allelic variation in cell wall candidate genes affecting solid wood properties in natural populations and land races of Pinus radiata

Forest trees are ideally suited to association mapping due to their high levels of diversity and low genomic linkage disequilibrium. Using an association mapping approach, single-nucleotide polymorphism (SNP) markers influencing quantitative variation in wood quality were identified in a natural population of Pinus radiata. Of 149 sites examined, 10 demonstrated significant associations (P < 0.05, q < 0.1) with one or more traits after accounting for population structure and experimentwise error. Without accounting for marker interactions, phenotypic variation attributed to individual SNPs ranged from 2 to 6.5%. Undesirable negative correlations between wood quality and growth were not observed, indicating potential to break negative correlations by selecting for individual SNPs in breeding programs. Markers that yielded significant associations were reexamined in an Australian land race. SNPs from three genes (PAL1, PCBER, and SUSY) yielded significant associations. Importantly, associations with two of these genes validated associations with density previously observed in the discovery population. In both cases, decreased wood density was associated with the minor allele, suggesting that these SNPs may be under weak negative purifying selection for density in the natural populations. These results demonstrate the utility of LD mapping to detect associations, even when the power to detect SNPs with small effect is anticipated to be low.

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.

Forest-tree population genomics and adaptive evolution

Forest trees have gained much attention in recent years as nonclassical model eukaryotes for population, evolutionary and ecological genomic studies. Because of low domestication, large open-pollinated native populations, and high levels of both genetic and phenotypic variation, they are ideal organisms to unveil the molecular basis of population adaptive divergence in nature. Population genomics, in its broad-sense definition, is an emerging discipline that combines genome-wide sampling with traditional population genetic approaches to understanding evolution. Here we briefly review traditional methods of studying adaptive genetic variation in forest trees, and describe a new, integrated population genomics approach. First, alleles (haplotypes) at candidate genes for adaptive traits and their effects on phenotypes need to be characterized via sequencing and association mapping. At this stage, functional genomics can assist in understanding gene action and regulation by providing detailed transcriptional profiles. Second, frequencies of alleles in native populations for causative single-nucleotide polymorphisms are estimated to identify patterns of adaptive variation across heterogeneous environments. Population genomics, through deciphering allelic effects on phenotypes and identifying patterns of adaptive variation at the landscape level, will in the future constitute a useful tool, if cost-effective, to design conservation strategies for forest trees.

Evolutionary relationships among Pinus (Pinaceae) subsections inferred from multiple low-copy nuclear loci

Sequence data from nrITS and cpDNA have failed to fully resolve phylogenetic relationships among Pinus species. Four low-copy nuclear genes, developed from the screening of 73 mapped conifer anchor loci, were sequenced from 12 species representing all subsections. Individual loci do not uniformly support either the nrITS or cpDNA hypotheses and in some cases produce unique topologies. Combined analysis of low-copy nuclear loci produces a well-supported subsectional topology of two subgenera, each divided into two sections, congruent with prior hypotheses of deep divergence in Pinus. The placements of P. nelsonii, P. krempfii, and P. contorta have been of continued systematic interest. Results strongly support the placement of P. nelsonii as sister to the remaining members of sect. Parrya, suggest a moderately well-supported and consistent position of P. krempfii as sister to the remaining members of sect. Quinquefoliae, and are ambiguous about the placement of P. contorta. A successful phylogenetic strategy in Pinus will require many low-copy nuclear loci that include a high proportion of silent sites and derive from independent linkage groups. The locus screening and evaluation strategy presented here can be broadly applied to facilitate the development of phylogenetic markers from the increasing number of available genomic resources.

A composite linkage map from two crosses for the species complex Picea mariana x Picea rubens and analysis of synteny with other Pinaceae

Four individual linkage maps were constructed from two crosses for the species complex Picea mariana (Mill.) B.S.P. x Picea rubens Sarg in order to integrate their information into a composite map and to compare with other Pinaceae. For all individual linkage maps, 12 major linkage groups were recovered with 306 markers per map on average. Before building the composite linkage map, the common male parent between the two crosses made it possible to construct a reference linkage map to validate the relative position of homologous markers. The final composite map had a length of 2,319 cM (Haldane) and contained a total of 1,124 positioned markers, including 1,014 AFLPs, 3 RAPDs, 53 SSRs, and 54 ESTPs, assembled into 12 major linkage groups. Marker density of the composite map was statistically homogenous and was much higher (one marker every 2.1 cM) than that of the individual linkage maps (one marker every 5.7 to 7.1 cM). Synteny was well conserved between individual, reference, and composite linkage maps and 94% of homologous markers were colinear between the reference and composite maps. The combined information from the two crosses increased by about 24% the number of anchor markers compared to the information from any single cross. With a total number of 107 anchor markers (SSRs and ESTPs), the composite linkage map is a useful starting point for large-scale genome comparisons at the intergeneric level in the Pinaceae. Comparisons of this map with those in Pinus and Pseudotsuga allowed the identification of one breakdown in synteny where one linkage group homologous to both Picea and Pinus corresponded to two linkage groups in Pseudotsuga. Implications for the evolution of the Pinaceae genome are discussed.

Identification of quantitative trait loci for wood and fibre properties in two full-sib properties of Eucalyptus globulus

Regions of the genome influencing wood and fibre traits in Eucalyptus globulus Labill. have been identified in two full-sib pedigrees that share a common male parent. The first pedigree, cross A, contains 148 progeny, and the second pedigree, cross B, contains 135 progeny. Subsets of progeny of these two controlled crosses were planted at seven sites throughout Australia in 1990. Wood cores were taken at 0.9 m above ground in 1997, and wood and fibre traits were analysed for each individual. Three quantitative trait loci (QTL) affecting wood density, one QTL affecting pulp yield and one QTL affecting microfibril angle have been located in both pedigrees, using single-factor analysis of variance. Other QTLs affecting these traits, as well as fibre length and cellulose content, were located in cross A only.

Detection and verification of quantitative trait loci for resistance to Dothistroma needle blight in Pinus radiata

Six related radiata pine ( Pinus radiata) full-sib families were used to detect and independently verify quantitative trait loci (QTLs) for resistance to Dothistroma needle blight, caused by Dothistroma septospora. The detection families had from 26 to 30 individuals each, and had either a common maternal (31053) or paternal (31032) parent; one family (cross 4) consisted of progeny from both parents, 31053 x 31032. Approximately 200 additional progeny from cross 4 were clonally replicated and planted at two sites, with at least five to seven ramets of each individual per site. Marker segregation data were collected from a total of 250 RFLP and microsatellite markers, and single factor ANOVAs were conducted separately for each family and marker. A number of putative associations were observed, some across more than one family. Permutation tests were used to confirm expected probabilities of multiple associations based on chance alone. Seven markers representing at least four QTLs for resistance to Dothistroma were identified as being significant in more than one family; one of these was significant at P<0.05 in three families and highly significant at P<0.01 in a fourth. Further confirmation was obtained by testing those markers that were significant in more than one of the detection families (or highly significant in cross 4) in the clonally replicated progeny from cross 4. Four QTL positions were verified in the clonal populations, with a total percent variation accounted for of 12.5.