Using terrestrial laser scanning to evaluate non-destructive aboveground biomass allometries in diverse Northern California forests

A crucial part of carbon accounting is quantifying a tree’s aboveground biomass (AGB) using allometric equations, but species-specific equations are limited because data to inform these equations requires destructive harvesting of many trees which is difficult and time-consuming. Here, we used terrestrial laser scanning (TLS) to non-destructively estimate AGB for 282 trees from 5 species at 3 locations in Northern California using stem and branch volume estimates from quantitative structure models (QSMs) and wood density from the literature. We then compared TLS QSM estimates of AGB with published allometric equations and used TLS-based AGB, diameter at breast height (DBH), and height to derive new species-specific allometric AGB equations for our study species. To validate the use of TLS, we used traditional forestry approaches to collect DBH (n = 550) and height (n = 291) data on individual trees. TLS-based DBH and height were not significantly different from field inventory data (R2 = 0.98 for DBH, R2 = 0.95 for height). Across all species, AGB calculated from TLS QSM volumes were approximately 30% greater than AGB estimates using published Forest Service’s Forest Inventory and Analysis Program equations, and TLS QSM AGB estimates were 10% greater than AGB calculated with existing equations, although this variation was species-dependent. In particular, TLS AGB estimates for Quercus agrifolia and Sequoia sempervirens differed the most from AGB estimates calculated using published equations. New allometric equations created using TLS data with DBH and height performed better than equations that only included DBH and matched most closely with AGB estimates generated from QSMs. Our results support the use of TLS as a method to rapidly estimate height, DBH, and AGB of multiple trees at a plot-level when species are identified and wood density is known. In addition, the creation of new TLS-based non-destructive allometric equations for our 5 study species may have important applications and implications for carbon quantification over larger spatial scales, especially since our equations estimated greater AGB than previous approaches

Multi-trait ensemble genomic prediction and simulations of recurrent selection highlight importance of complex trait genetic architecture for long-term genetic gains in wheat

Cereal crop breeders have achieved considerable genetic gain in genetically complex traits, such as grain yield, while maintaining genetic diversity. However, focus on selection for yield has negatively impacted other important traits. To better understand multi-trait selection within a breeding context, and how it might be optimized, we analysed genotypic and phenotypic data from a genetically diverse, 16-founder wheat multi-parent advanced generation inter-cross population. Compared to single-trait models, multi-trait ensemble genomic prediction models increased prediction accuracy for almost 90 % of traits, improving grain yield prediction accuracy by 3–52 %. For complex traits, non-parametric models (Random Forest) also outperformed simplified, additive models (LASSO), increasing grain yield prediction accuracy by 10–36 %. Simulations of recurrent genomic selection then showed that sustained greater forward prediction accuracy optimized long-term genetic gains. Simulations of selection on grain yield found indirect responses in related traits, involving optimized antagonistic trait relationships. We found multi-trait selection indices could effectively optimize undesirable relationships, such as the trade-off between grain yield and protein content, or combine traits of interest, such as yield and weed competitive ability. Simulations of phenotypic selection found that including Random Forest rather than LASSO genetic models, and multi-trait rather than single-trait models as the true genetic model accelerated and extended long-term genetic gain whilst maintaining genetic diversity. These results (i) suggest important roles of pleiotropy and epistasis in the wider context of wheat breeding programmes, and (ii) provide insights into mechanisms for continued genetic gain in a limited genepool and optimization of multiple traits for crop improvement

Limited haplotype diversity underlies polygenic trait architecture across 70 years of wheat breeding

Background Selection has dramatically shaped genetic and phenotypic variation in bread wheat. We can assess the genomic basis of historical phenotypic changes, and the potential for future improvement, using experimental populations that attempt to undo selection through the randomizing effects of recombination. Results We bred the NIAB Diverse MAGIC multi-parent population comprising over 500 recombinant inbred lines, descended from sixteen historical UK bread wheat varieties released between 1935 and 2004. We sequence the founders’ genes and promoters by capture, and the MAGIC population by low-coverage whole-genome sequencing. We impute 1.1 M high-quality SNPs that are over 99% concordant with array genotypes. Imputation accuracy only marginally improves when including the founders’ genomes as a haplotype reference panel. Despite capturing 73% of global wheat genetic polymorphism, 83% of genes cluster into no more than three haplotypes. We phenotype 47 agronomic traits over 2 years and map 136 genome-wide significant associations, concentrated at 42 genetic loci with large and often pleiotropic effects. Around half of these overlap known quantitative trait loci. Most traits exhibit extensive polygenicity, as revealed by multi-locus shrinkage modelling. Conclusions Our results are consistent with a gene pool of low haplotypic diversity, containing few novel loci of large effect. Most past, and projected future, phenotypic changes arising from existing variation involve fine-scale shuffling of a few haplotypes to recombine dozens of polygenic alleles of small effect. Moreover, extensive pleiotropy means selection on one trait will have unintended consequences, exemplified by the negative trade-off between yield and protein content, unless selection and recombination can break unfavorable trait-trait associations

An eight-parent multiparent advanced generation inter-cross population for winter-sown wheat: creation, properties, and validation

MAGIC populations represent one of a new generation of crop genetic mapping resources combining high genetic recombination and diversity. We describe the creation and validation of an eight-parent MAGIC population consisting of 1091 F7 lines of winter-sown wheat (Triticum aestivum L.). Analyses based on genotypes from a 90,000-single nucleotide polymorphism (SNP) array find the population to be well-suited as a platform for fine-mapping quantitative trait loci (QTL) and gene isolation. Patterns of linkage disequilibrium (LD) show the population to be highly recombined; genetic marker diversity among the founders was 74% of that captured in a larger set of 64 wheat varieties, and 54% of SNPs segregating among the 64 lines also segregated among the eight founder lines. In contrast, a commonly used reference bi-parental population had only 54% of the diversity of the 64 varieties with 27% of SNPs segregating. We demonstrate the potential of this MAGIC resource by identifying a highly diagnostic marker for the morphological character "awn presence/absence" and independently validate it in an association-mapping panel. These analyses show this large, diverse, and highly recombined MAGIC population to be a powerful resource for the genetic dissection of target traits in wheat, and it is well-placed to efficiently exploit ongoing advances in phenomics and genomics. Genetic marker and trait data, together with instructions for access to seed, are available at http://www.niab.com/MAGIC/

Multi-trait ensemble genomic prediction and simulations of recurrent selection highlight importance of complex trait genetic architecture for long-term genetic gains in wheat

Cereal crop breeders have achieved considerable genetic gain in genetically complex traits, such as grain yield, while maintaining genetic diversity. However, focus on selection for yield has negatively impacted other important traits. To better understand multi-trait selection within a breeding context, and how it might be optimized, we analysed genotypic and phenotypic data from a genetically diverse, 16-founder wheat multi-parent advanced generation inter-cross population. Compared to single-trait models, multi-trait ensemble genomic prediction models increased prediction accuracy for almost 90 % of traits, improving grain yield prediction accuracy by 3–52 %. For complex traits, non-parametric models (Random Forest) also outperformed simplified, additive models (LASSO), increasing grain yield prediction accuracy by 10–36 %. Simulations of recurrent genomic selection then showed that sustained greater forward prediction accuracy optimized long-term genetic gains. Simulations of selection on grain yield found indirect responses in related traits, involving optimized antagonistic trait relationships. We found multi-trait selection indices could effectively optimize undesirable relationships, such as the trade-off between grain yield and protein content, or combine traits of interest, such as yield and weed competitive ability. Simulations of phenotypic selection found that including Random Forest rather than LASSO genetic models, and multi-trait rather than single-trait models as the true genetic model accelerated and extended long-term genetic gain whilst maintaining genetic diversity. These results (i) suggest important roles of pleiotropy and epistasis in the wider context of wheat breeding programmes, and (ii) provide insights into mechanisms for continued genetic gain in a limited genepool and optimization of multiple traits for crop improvement

TLS2trees: A scalable tree segmentation pipeline for TLS data

1. Above-ground biomass (AGB) is an important metric used to quantify the mass of carbon stored in terrestrial ecosystems. For forests, this is routinely estimated at the plot scale (typically 1 ha) using inventory measurements and allometry. In recent years, terrestrial laser scanning (TLS) has appeared as a disruptive technology that can generate a more accurate assessment of tree and plot scale AGB; however, operationalising TLS methods has had to overcome a number of challenges. One such challenge is the segmentation of individual trees from plot level point clouds that are required to estimate woody volume, this is often done manually (e.g. with interactive point cloud editing software) and can be very time consuming. / 2. Here we present TLS2trees, an automated processing pipeline and set of Python command line tools that aims to redress this processing bottleneck. TLS2trees consists of existing and new methods and is specifically designed to be horizontally scalable. The processing pipeline is demonstrated on 7.5 ha of TLS data captured across 10 plots of seven forest types; from open savanna to dense tropical rainforest. / 3. A total of 10,557 trees are segmented with TLS2trees: these are compared to 1281 manually segmented trees. Results indicate that TLS2trees performs well, particularly for larger trees (i.e. the cohort of largest trees that comprise 50% of total plot volume), where plot-wise tree volume bias is ±0.4 m3 and %RMSE is 60%. Segmentation performance decreases for smaller trees, for example where DBH ≤10 cm; a number of reasons are suggested including performance of semantic segmentation step. / 4. The volume and scale of TLS data captured in forest plots is increasing. It is suggested that to fully utilise this data for activities such as monitoring, reporting and verification or as reference data for satellite missions an automated processing pipeline, such as TLS2trees, is required. To facilitate improvements to TLS2trees, as well as modification for other laser scanning modes (e.g. mobile and UAV laser scanning), TLS2trees is a free and open-source software

SHRiMP: Accurate Mapping of Short Color-space Reads

The development of Next Generation Sequencing technologies, capable of sequencing hundreds of millions of short reads (25–70 bp each) in a single run, is opening the door to population genomic studies of non-model species. In this paper we present SHRiMP - the SHort Read Mapping Package: a set of algorithms and methods to map short reads to a genome, even in the presence of a large amount of polymorphism. Our method is based upon a fast read mapping technique, separate thorough alignment methods for regular letter-space as well as AB SOLiD (color-space) reads, and a statistical model for false positive hits. We use SHRiMP to map reads from a newly sequenced Ciona savignyi individual to the reference genome. We demonstrate that SHRiMP can accurately map reads to this highly polymorphic genome, while confirming high heterozygosity of C. savignyi in this second individual. SHRiMP is freely available at http://compbio.cs.toronto.edu/shrimp

Using next-generation sequencing for high resolution multiplex analysis of copy number variation from nanogram quantities of DNA from formalin-fixed paraffin-embedded specimens

The use of next-generation sequencing technologies to produce genomic copy number data has recently been described. Most approaches, however, reply on optimal starting DNA, and are therefore unsuitable for the analysis of formalin-fixed paraffin-embedded (FFPE) samples, which largely precludes the analysis of many tumour series. We have sought to challenge the limits of this technique with regards to quality and quantity of starting material and the depth of sequencing required. We confirm that the technique can be used to interrogate DNA from cell lines, fresh frozen material and FFPE samples to assess copy number variation. We show that as little as 5 ng of DNA is needed to generate a copy number karyogram, and follow this up with data from a series of FFPE biopsies and surgical samples. We have used various levels of sample multiplexing to demonstrate the adjustable resolution of the methodology, depending on the number of samples and available resources. We also demonstrate reproducibility by use of replicate samples and comparison with microarray-based comparative genomic hybridization (aCGH) and digital PCR. This technique can be valuable in both the analysis of routine diagnostic samples and in examining large repositories of fixed archival material