Search

The selection of drought-tolerant genotypes is globally recognized as an effective strategy to maintain the growth and survival of commercial tree species exposed to future drought periods. New genomic selection tools that reduce the time of progeny trials are required to substitute traditional tree breeding programs. We investigated the genetic variation of water stress tolerance in New Zealand-grown Pinus radiata D. Don using 622 commercially-used genotypes from 63 families. We used quantitative pedigree-based (Genomic Best Linear Unbiased Prediction or ABLUP) and genomic-based (Genomic Best Linear Unbiased Prediction or GBLUP) approaches to examine the heritability estimates associated with water stress tolerance in P. radiata. Tree seedling growth traits, foliar carbon isotope composition (δ13C), and dark-adapted chlorophyll fluorescence (Y) were monitored before, during and after 10 months of water stress. Height growth showed a constant and moderate heritability level, while the heritability estimate for diameter growth and δ13C decreased with water stress. In contrast, chlorophyll fluorescence exhibited low heritability after 5 and 10 months of water stress. The GBLUP approach provided less breeding value accuracy than ABLUP, however, the relative selection efficiency of GBLUP was greater compared with ABLUP selection techniques. Although there was no significant relationship directly between δ13C and Y, the genetic correlations were significant and stronger for GBLUP. The positive genetic correlations between δ13C and tree biomass traits under water stress indicated that intraspecific variation in δ13C was likely driven by differences in the genotype’s photosynthetic capacity. The results show that foliar δ13C can predict P. radiata genotype tolerance to water stress using ABLUP and GBLUP approaches and that such approaches can provide a faster screening and selection of drought-tolerant genotypes for forestry breeding programs.

The recovery of genealogy in both natural and captive populations is critical for any decision in the management of genetic resources. It allows for the estimation of genetic parameters such as heritability and genetic correlations, as well as defining an optimal mating design that maintains a large effective population size. We utilised “genotyping-by-sequencing” (GBS) in combination with bioinformatics tools developed specifically for GBS data to recover genetic relatedness, with a focus on parent-offspring relationships in a Eucalyptus nitens breeding population as well as recognition of individuals representing other Eucalyptus species and putative hybrids. We found a clear advantage on using tools specifically designed for data of highly variable sequencing quality when recovering genetic relatedness. The parent-offspring relatedness showed a significant response to data filtering from 0.05 to 0.3 when the standard approach (G1) was used, while it oscillated around 0.4 when the specifically designed method (G5) was implemented. Additionally, comparisons with commonly used tools demonstrated vulnerability of the relatedness estimates to incorrect imputation of missing data when shallow sequencing information and genetically distant individuals are present in the population. In turn, these biased imputed genotypes negatively affected the estimation of genetic relatedness between parents and offspring. Careful filtering for both genetic outliers and shallowly sequenced markers led to improvements in estimations of genetic relatedness. Alternatively, a method that avoided missing data imputation and took sequence depth into consideration improved the accuracy of parent-offspring relationship coefficients where sequencing data quality was highly variable.

BackgroundNon-key traits (NKTs) in radiata pine (Pinus radiataD. Don) refer to traits other than growth, wood density and stiffness, but still of interest to breeders. Branch-cluster frequency, stem straightness, external resin bleeding and internal checking are examples of such traits and are targeted for improvement in radiata pine research programmes. Genomic selection can be conducted before the performance of selection candidates is available so that generation intervals can be reduced. Radiata pine is a species with a long generation interval, which if reduced could significantly increase genetic gain per unit of time. The aim of this study was to evaluate the accuracy and predictive ability of genomic selection and its efficiency over traditional forward selection in radiata pine for the following NKTs: branch-cluster frequency, stem straightness, internal checking, and external resin bleeding.ResultsNine hundred and eighty-eight individuals were genotyped using exome capture genotyping by sequencing (GBS) and 67,168 single nucleotide polymorphisms (SNPs) used to develop genomic estimated breeding values (GEBVs) with genomic best linear unbiased prediction (GBLUP). The documented pedigree was corrected using a subset of 704 SNPs. The percentage of trio parentage confirmed was about 49% and about 50% of parents were re-assigned. The accuracy of GEBVs was 0.55–0.75 when using the documented pedigree and 0.61–0.80 when using the SNP-corrected pedigree. A higher percentage of additive genetic variance was explained and a higher predictive ability was observed when using the SNP-corrected pedigree than using the documented pedigree. With the documented pedigree, genomic selection was similar to traditional forward selection when assuming a generation interval of 17 years, but worse than traditional forward selection when assuming a generation interval of 14 years. After the pedigree was corrected, genomic selection led to 37–115% and 13–77% additional genetic gain over traditional forward selection when generation intervals of 17 years and 14 years were assumed, respectively.ConclusionIt was concluded that genomic selection with a pedigree corrected by SNP information was an efficient way of improving non-key traits in radiata pine breeding.

Background Effective matching of genotypes and environments is required for the species to reach optimal productivity and act effectively for carbon sequestration. A common garden experiment across five different environments was undertaken to assess genotype x environment interaction (GxE) of coast redwood in order to understand the performance of genotypes across environments. Results The quantitative genetic analysis discovered no GxE between investigated environments for diameter at breast height (DBH). However, no genetic component was detected at one environment possibly due to stressful conditions. The implementation of universal response function allowed for the identification of important environmental factors affecting species productivity. Additionally, this approach enabled us to predict the performance of species across the New Zealand environmental conditions. Conclusions In combination with quantitative genetic analysis which identified genetically superior material, the URF model can directly identify the optimal geographical regions to maximize productivity. However, the finding of ideally uncorrelated climatic variables for species with narrow ecological amplitude is rather challenging, which complicates construction of informative URF model. This, along with a small number of tested environments, tended to overfit a prediction model which resulted in extreme predictions in untested environments.

Genomic selection is expected to enhance the genetic improvement of forest tree species by providing more accurate estimates of breeding values through marker-based relationship matrices compared with pedigree-based methodologies. When adequately robust genomic prediction models are available, an additional increase in genetic gains can be made possible with the shortening of the breeding cycle through elimination of the progeny testing phase and early selection of parental candidates. The potential of genomic selection was investigated in an advanced Eucalyptus nitens breeding population focused on improvement for solid wood production. A high-density SNP chip (EUChip60K) was used to genotype 691 individuals in the breeding population, which represented two seed orchards with different selection histories. Phenotypic records for growth and form traits at age six, and for wood quality traits at age seven were available to build genomic prediction models using GBLUP, which were compared to the traditional pedigree-based alternative using BLUP. GBLUP demonstrated that breeding value accuracy would be improved and substantial increases in genetic gains towards solid wood production would be achieved. Cross-validation within and across two different seed orchards indicated that genomic predictions would likely benefit in terms of higher predictive accuracy from increasing the size of the training data sets through higher relatedness and better utilization of LD

Background Many conifer breeding programs are paying increasing attention to breeding for resistance to needle disease due to the increasing importance of climate change. Phenotyping of traits related to resistance has many biological and temporal constraints that can often confound the ability to achieve reliable phenotypes and consequently, reliable genetic progress. The development of next generation sequencing platforms has also enabled implementation of genomic approaches in species lacking robust reference genomes. Genomic selection is, therefore, a promising strategy to overcome the constraints of needle disease phenotyping. Results We found high accuracy in the prediction of genomic breeding values in the disease-related traits that were well characterized, reaching 0.975 for genotyped individuals and 0.587 for non-genotyped individuals. This compared well with pedigree-based accuracies of up to 0.746. Surprisingly, poorly phenotyped disease traits also showed very high accuracy in terms of correlation of predicted genomic breeding values with pedigree-based counterparts. However, this was likely caused by the fact that both were clustered around the population mean, while deviations from the population mean caused by genetic effects did not appear to be well described. Caution should therefore be taken with the interpretation of results in poorly phenotyped traits. Conclusions Implementation of genomic selection in this test population of Pinus radiata resulted in a relatively high prediction accuracy of needle loss due to Dothistroma septosporum compared with a pedigree-based approach. Using genomics to avoid biological/temporal constraints where phenotyping is reliable appears promising. Unsurprisingly, reliable phenotyping, resulting in good heritability estimates, is a fundamental requirement for the development of a reliable prediction model. Furthermore, our results are also specific to the single pathogen mating-type that is present in New Zealand, and may change with future incursion of other pathogen varieties. There is no doubt, however, that once a robust genomic prediction model is built, it will be invaluable to not only select for host tolerance, but for other economically important traits simultaneously. This tool will thus future-proof our forests by mitigating the risk of disease outbreaks induced by future changes in climate.

BackgroundForest trees can occupy extensive geography and environmentally highly variable areas which result in high genetic variability in the direction of pressure from natural selection. At the same time, the majority of conifer species are wind-pollinated from both short and long distances, resulting in wide-spread gene flow, which can lead to maladaptation to local conditions. Quantitative analyses of provenance/progeny tests correct for genetic differences between populations to ensure unbiased genetic parameters are obtained. Commonly, the provenance effect is fitted as a fixed term or can be implemented as a contemporary group in the pedigree.ResultsThe use of a provenance effect, either as a fixed term or as the same contemporary groups in both maternal and paternal sides of the pedigree, resulted in fairly similar precision of genetic parameters in our case. However, when we developed a phantom contemporary group for the paternal side of the pedigree that considered a different genetic quality of pollen compared with the maternal contribution from trees in the local environment, the model fit and accuracy of breeding values increased.ConclusionConsideration of the mating dynamics and the vector of gene flow are important factors in modelling contemporary genetic groups, particularly when implementing pedigrees within a mixed model framework to obtain unbiased estimates of genetic parameters. This approach is especially important in traits involved in local adaptation.

There is a wide diversity of bioinformatic tools available for the assembly of next generation sequence and subsequence variant calling to identify genetic markers at scale. Integration of genomics tools such as genomic selection, association studies, pedigree analysis and analysis of genetic diversity, into operational breeding is a goal for New Zealand's most widely planted exotic tree species, Pinus radiata. In the absence of full reference genomes for large megagenomes such as in conifers, RNA sequencing in a range of genotypes and tissue types, offers a rich source of genetic markers for downstream application. We compared nine different assembler and variant calling software combinations in a single transcriptomic library and found that Single Nucleotide Polymorphism (SNPs) discovery could vary by as much as an order of magnitude (8,061 SNPs up to 86,815 SNPs). The assembler with the best realignment of the packages trialled, Trinity, in combination with several variant callers was then applied to a much larger multi-genotype, multi-tissue transcriptome and identified 683,135 in silico SNPs across a predicted 449,951 exons when mapped to the Pinus taeda ver 1.01e reference.

Open-pollinated (OP) mating is frequently used in forest tree breeding due to the relative temporal and financial efficiency of the approach. The trade-off is the lower precision of the estimated genetic parameters. Pedigree/sib-ship reconstruction has been proven as a tool to correct and complete pedigree information and to improve the precision of genetic parameter estimates. Our study analyzed an advanced generation Eucalyptus population from an OP breeding program using single-step genetic evaluation. The relationship matrix inferred from sib-ship reconstruction was used to rescale the marker-based relationship matrix (G matrix). This was compared with a second scenario that used rescaling based on the documented pedigree. The proposed single-step model performed better with respect to both model fit and the theoretical accuracy of breeding values. We found that the prediction accuracy was superior when using the pedigree information only when compared with using a combination of the pedigree and genomic information. This pattern appeared to be mainly a result of accumulated unrecognized relatedness over several breeding cycles, resulting in breeding values being shrunk toward the population mean. Using biased, pedigree-based breeding values as the base with which to correlate predicted GEBVs, resulted in the underestimation of prediction accuracies. Using breeding values estimated on the basis of sib-ship reconstruction resulted in increased prediction accuracies of the genotyped individuals. Therefore, selection of the correct base for estimation of prediction accuracy is critical. The beneficial impact of sib-ship reconstruction using G matrix rescaling was profound, especially in traits with inbreeding depression, such as stem diameter.

Phenotyping is the accurate and precise physical description of organisms. Accurate and quantitative phenotyping underpins the delivery of benefits from genetic improvement programs in agriculture. In forest trees, phenotyping at an equivalent precision has been impossible because trees and forests are large, long-lived, and highly variable. These facts have restricted the delivery of genetic gains in forestry compared to other agricultural sectors. We describe a landscape-scale phenotyping platform that integrates remote sensing, spatial information systems, and genomics to facilitate the delivery of greater gains enabling forestry to catch up with other sectors. Combining remote sensing at a range of spatial and temporal scales with genomics will ultimately impact on tree breeding globally.

Accurate inference of relatedness between individuals in breeding population contributes to the precision of genetic parameter estimates, effectiveness of inbreeding management and the amount of genetic progress delivered from breeding programs. Pedigree reconstruction has been proven to be an efficient tool to correct pedigree errors and recover hidden relatedness in open pollinated progeny tests but the method can be limited by the lack of parental genotypes and the high proportion of alien pollen from outside the breeding population. Our study investigates the efficiency of sib-ship reconstruction in an advanced breeding population of Eucalyptus nitens with only partially tracked pedigree. The sib-ship reconstruction allowed the identification of selfs (4% of the sample) and the exploration of their potential effect on inbreeding depression in the traits studied. We detected signs of inbreeding depression in diameter at breast height and growth strain while no indications were observed in wood density, wood stiffness and tangential air-dry shrinkage. After the application of a corrected sib-ship relationship matrix, additive genetic variance and heritability were observed to increase where signs of inbreeding depression were initially detected. Conversely, the same genetic parameters for traits that appeared to be free of inbreeding depression decreased in size. It therefore appeared that greater genetic variance may be due, at least in part, to contributions from inbreeding in these studied populations rather than a removal of inbreeding as is traditionally thought.

Genomic selection is expected to enhance the genetic improvement of forest tree species by providing more accurate estimates of breeding values through marker-based relationship matrices compared with pedigree-based methodologies. When adequately robust genomic prediction models are available, an additional increase in genetic gains can be made possible with the shortening of the breeding cycle through elimination of the progeny testing phase and early selection of parental candidates. The potential of genomic selection was investigated in an advanced Eucalyptus nitens breeding population focused on improvement for solid wood production. A high-density SNP chip (EUChip60K) was used to genotype 691 individuals in the breeding population, which represented two seed orchards with different selection histories. Phenotypic records for growth and form traits at age six, and for wood quality traits at age seven were available to build genomic prediction models using GBLUP, which were compared to the traditional pedigree-based alternative using BLUP. GBLUP demonstrated that breeding value accuracy would be improved and substantial increases in genetic gains towards solid wood production would be achieved. Cross-validation within and across two different seed orchards indicated that genomic predictions would likely benefit in terms of higher predictive accuracy from increasing the size of the training data sets through higher relatedness and better utilization of LD.

Pedigree reconstruction using molecular markers enables efficient management of inbreeding in open-pollinated breeding strategies, replacing expensive and time-consuming controlled pollination. This is particularly useful in preferentially outcrossed, insect pollinated Eucalypts known to suffer considerable inbreeding depression from related matings. A single nucleotide polymorphism (SNP) marker panel consisting of 106 markers was selected for pedigree reconstruction from the recently developed high-density Eucalyptus Infinium SNP chip (EuCHIP60K). The performance of this SNP panel for pedigree reconstruction in open-pollinated progenies of two Eucalyptus nitens seed orchards was compared with that of two microsatellite panels with 13 and 16 markers respectively. The SNP marker panel out-performed one of the microsatellite panels in the resolution power to reconstruct pedigrees and out-performed both panels with respect to data quality. Parentage of all but one offspring in each clonal seed orchard was correctly matched to the expected seed parent using the SNP marker panel, whereas parentage assignment to less than a third of the expected seed parents were supported using the 13-microsatellite panel. The 16-microsatellite panel supported all but one of the recorded seed parents, one better than the SNP panel, although there was still a considerable level of missing and inconsistent data. SNP marker data was considerably superior to microsatellite data in accuracy, reproducibility and robustness. Although microsatellites and SNPs data provide equivalent resolution for pedigree reconstruction, microsatellite analysis requires more time and experience to deal with the uncertainties of allele calling and faces challenges for data transferability across labs and over time. While microsatellite analysis will continue to be useful for some breeding tasks due to the high information content, existing infrastructure and low operating costs, the multi-species SNP resource available with the EuCHIP60k, opens a whole new array of opportunities for high-throughput, genome-wide or targeted genotyping in species of Eucalyptus.

The indole-diterpene paxilline is a potent tremorgenic mammalian mycotoxin and a known inhibitor of maxi-K ion channels. The gene cluster responsible for paxilline biosynthesis in Penicillium paxilli was identified by mapping four large plasmid-induced chromosome deletions. The cluster is predicted to lie within a 50 kb region of chromosome Va and to contain 17 genes, including a geranylgeranyl pyrophosphate (GGPP) synthase (paxG), two FAD-dependent monooxygenases (paxM and N), two cytochrome P450 monooxygenases (paxP and Q), a dimethylallyltryptophan (DMAT) synthase (paxD) and two possible transcription factors (paxR and paxS), which contain a Zn(II)2Cys6 DNA-binding motif. Targeted replacement of paxG confirmed that it is essential for paxilline biosynthesis but dispensable for growth. The GGPP for primary metabolism is predicted to be provided by a second GGPP synthase (ggs1) that was cloned, sequenced and mapped to chromosome IV. Semi-quantitative reverse transcriptase-polymerase chain reaction analysis demonstrated that the expression of paxG, paxM and paxP in submerged liquid cultures of P. paxilli increased dramatically with the onset of paxilline biosynthesis. In contrast, the expression of beta-tubulin (tub2) and ggs1 was not induced. This is the first description of the molecular cloning and genetic analysis of an indole-diterpene gene cluster.

Associations of single nucleotide polymorphisms (SNP) with the form traits, branch cluster frequency and stem straightness, were studied in three radiata pine (Pinus radiata D. Don) breeding populations. For branch cluster frequency, high genotype by environment (G×E) interactions were found between two sites of the POP1 trial series (Tarawera and Glenledi) and between two sites of the POP2 trial series (Tarawera and Woodhill). For stem straightness, high G×E interactions were found between two sites (Tarawera and Woodhill) in both the POP2 and POP3 trial series. Thirty-two and 26 SNPs were associated with branch cluster frequency and stem straightness respectively, across a number of sites and trial series in a genetic model assuming heterogeneous additive genetic variances across sites. These SNPs had significant associations with branch cluster frequency or stem straightness either in one site only, across sites within a trial series or across different trial series. Seven SNPs had significant associations with branch cluster frequency and seven SNPs with stem straightness across two of three trial series. One SNP had significant association with stem straightness across all three trial series. Strong G×E interactions resulted in fewer significant SNPs in common across the sites among which there were large effect differences. The percentage of genetic variance explained by SNPs showing significant associations ranged from 0.23 to 8.76 % for branch cluster frequency and from 0.37 to 12.75 % for stem straightness. This study showed that SNP effects for a trait change across environments if G×E interactions exist for that trait.

Standard protocols for extracting genomic DNA from Pinus radiata D. Don needles, such as CTAB-based methods, can yield large quantities of DNA. However, final DNA purity can be an issue due to carry over of contaminants that can impede accurate high throughput genotyping. This study evaluated eight DNA extraction and purification protocols to determine which method provided the greatest improvement in call rates and accuracy when using the Sequenom iPLEX® Gold MassARRAY® genotyping technology. Of the methods tested, genomic DNA extracted using the Machery-Nagel NucleoSpin®-96 Plant II kit performed the best overall, and was more efficiently and accurately genotyped than genomic DNA extracted using the standard CTAB method. This study also demonstrated that the quality and assay performance of CTAB-extracted genomic DNA is greatly improved by further purification with the Qiagen® QIAquick 96 PCR Purification kit. Using these improvements, the Sequenom iPLEX® Gold MassARRAY® genotyping technology is now a viable option for genotyping plant genomes such as Pinus radiata.