Identification of a Cis-acting regulatory polymorphism in a Eucalypt COBRA-like gene affecting cellulose content

Uncovering epigenetic and transcriptional regulation of growth in Douglas-fir: identification of differential methylation regions in mega-sized introns

Tree growth performance can be partly explained by genetics, while a large proportion of growth variation is thought to be controlled by environmental factors. However, to what extent DNA methylation, a stable epigenetic modification, contributes to phenotypic plasticity in the growth performance of long-lived trees remains unclear. In this study, a comparative analysis of targeted DNA genotyping, DNA methylation and mRNAseq profiling for needles of 44-year-old Douglas-fir trees (Pseudotsuga menziesii (Mirb.) Franco) having contrasting growth characteristics was performed. In total, we identified 195 differentially expressed genes (DEGs) and 115 differentially methylated loci (DML) that are associated with genes involved in fitness-related processes such as growth, stress management, plant development and energy resources. Interestingly, all four intronic DML were identified in mega-sized (between 100 and 180 kbp in length) and highly expressed genes, suggesting specialized regulation mechanisms of these long intron genes in gymnosperms. DNA repetitive sequences mainly comprising long-terminal repeats of retroelements are involved in growth-associated DNA methylation regulation (both hyper- and hypomethylation) of 99 DML (86.1% of total DML). Furthermore, nearly 14% of the DML was not tagged by single nucleotide polymorphisms, suggesting a unique contribution of the epigenetic variation in tree growth.

Genome-wide association studies: an intuitive solution for SNP identification and gene mapping in trees

Analysis of natural diversity in wild/cultivated plants can be used to understand the genetic basis for plant breeding programs. Recent advancements in DNA sequencing have expanded the possibilities for genetically altering essential features. There have been several recently disclosed statistical genetic methods for discovering the genes impacting target qualities. One of these useful methods is the genome-wide association study (GWAS), which effectively identifies candidate genes for a variety of plant properties by examining the relationship between a molecular marker (such as SNP) and a target trait. Conventional QTL mapping with highly structured populations has major limitations. The limited number of recombination events results in poor resolution for quantitative traits. Only two alleles at any given locus can be studied simultaneously. Conventional mapping approach fails to work in perennial plants and vegetatively propagated crops. These limitations are sidestepped by association mapping or GWAS. The flexibility of GWAS comes from the fact that the individuals being examined need not be linked to one another, allowing for the use of all meiotic and recombination events to increase resolution. Phenotyping, genotyping, population structure analysis, kinship analysis, and marker-trait association analysis are the fundamental phases of GWAS. With the rapid development of sequencing technologies and computational methods, GWAS is becoming a potent tool for identifying the natural variations that underlie complex characteristics in crops. The use of high-throughput sequencing technologies along with genotyping approaches like genotyping-by-sequencing (GBS) and restriction site associated DNA (RAD) sequencing may be highly useful in fast-forward mapping approach like GWAS. Breeders may use GWAS to quickly unravel the genomes through QTL and association mapping by taking advantage of natural variances. The drawbacks of conventional linkage mapping can be successfully overcome with the use of high-resolution mapping and the inclusion of multiple alleles in GWAS.

Engineering traits through CRISPR/cas genome editing in woody species to improve forest diversity and yield

Dangers confronting forest ecosystems are many and the strength of these biological systems is deteriorating, thus substantially affecting tree physiology, phenology, and growth. The establishment of genetically engineered trees into degraded woodlands, which would be adaptive to changing climate, could help in subsiding ecological threats and bring new prospects. This should not be resisted due to the apprehension of transgene dispersal in forests. Consequently, it is important to have a deep insight into the genetic structure and phenotypic limits of the reproductive capability of tree stands/population(s) to endure tolerance and survival. Importantly, for a better understanding of genes and their functional mechanisms, gene editing (GeEd) technology is an excellent molecular tool to unravel adaptation progressions. Therefore, GeEd could be harnessed for resolving the allelic interactions for the creation of gene diversity, and transgene dispersal may be alleviated among the population or species in different bioclimatic zones around the globe. This review highlights the potential of the CRISPR/Cas tools in genomic, transcriptomic, and epigenomic-based assorted and programmable alterations of genes in trees that might be able to fix the trait-specific gene function. Also, we have discussed the application of diverse forms of GeEd to genetically improve several traits, such as wood density, phytochemical constituents, biotic and abiotic stress tolerance, and photosynthetic efficiency in trees. We believe that the technology encourages fundamental research in the forestry sector besides addressing key aspects, which might fasten tree breeding and germplasm improvement programs worldwide.

CAGEs are Golgi-localized GT31 enzymes involved in cellulose biosynthesis in Arabidopsis

Cellulose is the main structural component in the plant cell walls. We show that two glycosyltransferase family 31 (GT31) enzymes of Arabidopsis thaliana, here named cellulose synthesis associated glycosyltransferases 1 and 2 (CAGE1 and 2), influence both primary and secondary cell wall cellulose biosynthesis. cage1cage2 mutants show primary cell wall defects manifesting as impaired growth and cell expansion in seedlings and etiolated hypocotyls, along with secondary cell wall defects, apparent as collapsed xylem vessels and reduced xylem wall thickness in the inflorescence stem. Single and double cage mutants also show increased sensitivity to the cellulose biosynthesis inhibitor isoxaben. The cage1cage2 phenotypes were associated with an approximately 30% reduction in cellulose content, an approximately 50% reduction in secondary cell wall CELLULOSE SYNTHASE (CESA) protein levels in stems and reduced cellulose biosynthesis rate in seedlings. CESA transcript levels were not significantly altered in cage1cage2 mutants, suggesting that the reduction in CESA levels was caused by a post-transcriptional mechanism. Both CAGE1 and 2 localize to the Golgi apparatus and are predicted to synthesize β-1,3-galactans on arabinogalactan proteins. In line with this, the cage1cage2 mutants exhibit reduced levels of β-Yariv binding to arabinogalactan protein linked β-1,3-galactan. This leads us to hypothesize that defects in arabinogalactan biosynthesis underlie the cellulose deficiency of the mutants.

Targeted re-sequencing and genome-wide association analysis for wood property traits in breeding population of Eucalyptus tereticornis × E. grandis

Globally, Eucalyptus plantations occupy 22 million ha area and is one of the preferred hardwood species due to their short rotation, rapid growth, adaptability and wood properties. In this study, we present results of GWAS in parents and 100 hybrids of Eucalyptus tereticornis × E. grandis using 762 genes presumably involved in wood formation. Comparative analysis between parents predicted 32,202 polymorphic SNPs with high average read depth of 269-562× per individual per nucleotide. Seventeen wood related traits were phenotyped across three diverse environments and GWAS was conducted using 13,610 SNPs. A total of 45 SNP-trait associations were predicted across two locations. Seven large effect markers were identified which explained more than 80% of phenotypic variation for fibre area. This study has provided an array of candidate genes which may govern fibre morphology in this genus and has predicted potential SNPs which can guide future breeding programs in tropical Eucalyptus.

Identification of Candidate Gene-Based Markers for Girth Growth in Rubber Trees

Girth growth is an important factor in both latex and timber production of the rubber tree. In this study, we performed candidate gene association mapping for girth growth in rubber trees using intron length polymorphism markers (ILP) in identifying the candidate genes responsible for girth growth. The COBL064_1 marker developed from the candidate gene (COBL4) regulating cellulose deposition and oriented cell expansion in the plant cell wall showed the strongest association with girth growth across two seasons in the Amazonian population and was validated in the breeding lines. We then applied single molecule real-time (SMRT) circular consensus sequencing (CCS) to analyze a wider gene region of the COBL4 to pinpoint the single nucleotide polymorphism (SNP) that best explains the association with the traits. A SNP in the 3’ UTR showing linkage disequilibrium with the COBL064_1 most associated with girth growth. This study showed that the cost-effective method of ILP gene-based markers can assist in identification of SNPs in the candidate gene associated with girth growth. The SNP markers identified in this study added useful markers for the improvement of girth growth in rubber tree breeding programs.

First record of off-season flowering in Populus deltoides from India: paradigm of climate change indicator

Populus deltoides is a fast-growing woody species possessing plethora of industrial applications. This species evolutionarily developed unisexual male and female catkin inflorescence on separate trees. Flowering usually occurs during early spring before the development of foliage, where buds appear near axils or at the extending shoots. In 2019, surveys were undertaken to study the flowering pattern of P. deltoides in the states of Punjab, Haryana, Uttar Pradesh and Uttarakhand in northern India. Interestingly, an anomalous flowering behaviour (appearance of off-season male catkins during autumn, i.e. October) was observed in a plantation trial at Kapurthala, Punjab. The male catkins were 2.7-3.1 ± 0.07 cm long and 0.3-0.5 ± 0.03 cm wide, which is significant for flowering and liberation of pollen grains. Preliminary results suggested that climatic factors, such as episodes of high or low temperature and the precipitation variation forcing the tree species to behave differently. Unearthing the climate-driven off-season flowering in other tree species alluded the stimulation of phytohormones, such as gibberellic and salicylic acid concentrations influencing the flowering time, therefore, needs further investigation in case of P. deltoides. Overall, this work provides early clues of changing climatic scenario altering the flowering pattern of a tropical forestry tree species.

Non-Destructive Evaluation Techniques and What They Tell Us about Wood Property Variation

To maximize utilization of our forest resources, detailed knowledge of wood property variation and the impacts this has on end-product performance is required at multiple scales (within and among trees, regionally). As many wood properties are difficult and time-consuming to measure our knowledge regarding their variation is often inadequate as is our understanding of their responses to genetic and silvicultural manipulation. The emergence of many non-destructive evaluation (NDE) methodologies offers the potential to greatly enhance our understanding of the forest resource; however, it is critical to recognize that any technique has its limitations and it is important to select the appropriate technique for a given application. In this review, we will discuss the following technologies for assessing wood properties both in the field: acoustics, Pilodyn, Resistograph and Rigidimeter and the lab: computer tomography (CT) scanning, DiscBot, near infrared (NIR) spectroscopy, radial sample acoustics and SilviScan. We will discuss these techniques, explore their utilization, and list applications that best suit each methodology. As an end goal, NDE technologies will help researchers worldwide characterize wood properties, develop accurate models for prediction, and utilize field equipment that can validate the predictions. The continued advancement of NDE technologies will also allow researchers to better understand the impact on wood properties on product performance.

A Role for Epigenetic Regulation in the Adaptation and Stress Responses of Non-model Plants

In recent years enormous progress has been made in understanding the role of epigenetic regulation response to environmental stimuli, especially in response to stresses. Molecular mechanisms involved in chromatin dynamics and silencing have been explained, leading to an appreciation of how new phenotypes can be generated quickly in response to environmental modifications. In some cases, it has also been shown that epigenetic modifications can be stably transmitted to the next generations. Despite this, the vast majority of studies have been carried out with model plants, particularly with Arabidopsis, and very little is known on how native plants in their natural habitat react to changes in their environment. Climate change has been affecting, sometimes drastically, the conditions of numerous ecosystems around the world, forcing populations of native species to adapt quickly. Although part of the adaptation can be explained by the preexisting genetic variation in the populations, recent studies have shown that new stable phenotypes can be generated through epigenetic modifications in few generations, contributing to the stability and survival of the plants in their natural habitat. Here, we review the recent data that suggest that epigenetic variation can help natural populations to cope to with change in their environments.

Genome-Wide Association Studies to Improve Wood Properties: Challenges and Prospects

Wood formation is an excellent model system for quantitative trait analysis due to the strong associations between the transcriptional and metabolic traits that contribute to this complex process. Investigating the genetic architecture and regulatory mechanisms underlying wood formation will enhance our understanding of the quantitative genetics and genomics of complex phenotypic variation. Genome-wide association studies (GWASs) represent an ideal statistical strategy for dissecting the genetic basis of complex quantitative traits. However, elucidating the molecular mechanisms underlying many favorable loci that contribute to wood formation and optimizing GWAS design remain challenging in this omics era. In this review, we summarize the recent progress in GWAS-based functional genomics of wood property traits in major timber species such as Eucalyptus, Populus, and various coniferous species. We discuss several appropriate experimental designs for extensive GWAS in a given undomesticated tree population, such as omics-wide association studies and high-throughput phenotyping technologies. We also explain why more attention should be paid to rare allelic and major structural variation. Finally, we explore the potential use of GWAS for the molecular breeding of trees. Such studies will help provide an integrated understanding of complex quantitative traits and should enable the molecular design of new cultivars.

Thirteen years under arid conditions: exploring marker-trait associations in Eucalyptus cladocalyx for complex traits related to flowering, stem form and growth

We present an association analysis for seven key traits related to flowering, stem form and growth in Eucalyptus cladocalyx, a tree species suitable for low rainfall sites, using a long-term progeny trial with 49 open-pollinated maternal families in the southern Atacama Desert, Chile. The progeny trial was carried out in an arid environment with a mean annual rainfall of 152 mm. Simple sequence repeats (SSR) from a full consensus map of Eucalyptus were used for genotyping 245 individual trees. Twenty-three significant marker-trait associations were identified, explaining between 5.9 and 23.7% of the phenotypic variance. The marker EMBRA101 located on LG10 at 56.5 cM was concomitantly associated with diameter at breast height and tree height. Nine SSR were significantly associated with stem forking and stem straightness, explaining between 5.9 and 14.8% of the phenotypic variation. To our knowledge, this is the first study reporting a SSR-based association mapping analysis for stem form traits in Eucalyptus. These results provide novel and valuable information for understanding the genetic base of key traits in E. cladocalyx for breeding purposes under arid conditions.

Increased fire frequency promotes stronger spatial genetic structure and natural selection at regional and local scales in Pinus halepensis Mill

BACKGROUND AND AIMS

The recurrence of wildfires is predicted to increase due to global climate change, resulting in severe impacts on biodiversity and ecosystem functioning. Recurrent fires can drive plant adaptation and reduce genetic diversity; however, the underlying population genetic processes have not been studied in detail. In this study, the neutral and adaptive evolutionary effects of contrasting fire regimes were examined in the keystone tree species Pinus halepensis Mill. (Aleppo pine), a fire-adapted conifer. The genetic diversity, demographic history and spatial genetic structure were assessed at local (within-population) and regional scales for populations exposed to different crown fire frequencies.

METHODS

Eight natural P. halepensis stands were sampled in the east of the Iberian Peninsula, five of them in a region exposed to frequent crown fires (HiFi) and three of them in an adjacent region with a low frequency of crown fires (LoFi). Samples were genotyped at nine neutral simple sequence repeats (SSRs) and at 251 single nucleotide polymorphisms (SNPs) from coding regions, some of them potentially important for fire adaptation.

KEY RESULTS

Fire regime had no effects on genetic diversity or demographic history. Three high-differentiation outlier SNPs were identified between HiFi and LoFi stands, suggesting fire-related selection at the regional scale. At the local scale, fine-scale spatial genetic structure (SGS) was overall weak as expected for a wind-pollinated and wind-dispersed tree species. HiFi stands displayed a stronger SGS than LoFi stands at SNPs, which probably reflected the simultaneous post-fire recruitment of co-dispersed related seeds. SNPs with exceptionally strong SGS, a proxy for microenvironmental selection, were only reliably identified under the HiFi regime.

CONCLUSIONS

An increasing fire frequency as predicted due to global change can promote increased SGS with stronger family structures and alter natural selection in P. halepensis and in plants with similar life history traits.

Uncovering SNP and indel variations of tetraploid cottons by SLAF-seq

Background

Cotton (Gossypium spp.), as the world’s most utilized textile fibre source, is an important, economically valuable crop worldwide. Understanding the genomic variation of tetraploid cotton species is important for exploitation of the excellent characteristics of wild cotton and for improving the diversity of cotton in breeding. However, the discovery of DNA polymorphisms in tetraploid cotton genomes has lagged behind other important crops.

Results

A total of 111,795,823 reads, 467,735 specific length amplified fragment (SLAF) tags and 139,176 high-quality DNA polymorphisms were identified using specific length amplified fragment sequencing (SLAF-seq), including 132,880 SNPs and 6,296 InDels between the reference genome (TM-1) and the five tetraploid cotton species. Intriguingly, gene ontology (GO) enrichment analysis revealed that a number of significant terms were related to reproduction in G. barbadense acc. 3–79. Based on the new data sets, we reconstructed phylogenetic trees that showed a high concordance to the phylogeny of diploid and polyploid cottons. A large amount of interspecific genetic variations were identified, and some of them were validated by the single-strand conformation polymorphism (SSCP) method, which will be applied in introgression genetics and breeding with G. hirsutum cv. Emian22 as the receptor and the other species as donors.

Conclusions

Using SLAF-seq, a large number of DNA polymorphisms were identified. The comprehensive analysis of DNA polymorphisms provided invaluable insights into the different tetraploid cotton species. More importantly, the identification of numerous interspecific genetic variations provides the basis and is very practical for future introgression breeding. The results presented herein provide a valuable genomic resource for new insights into the genetics and breeding of cotton.

Electronic supplementary material

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

MLPA-Based Analysis of Copy Number Variation in Plant Populations

Copy number variants (CNVs) are intraspecies duplications/deletions of large DNA segments (>1 kb). A growing number of reports highlight the functional and evolutionary impact of CNV in plants, increasing the need for appropriate tools that enable locus-specific CNV genotyping on a population scale. Multiplex ligation-dependent probe amplification (MLPA) is considered a gold standard in genotyping CNV in humans. Consequently, numerous commercial MLPA assays for CNV-related human diseases have been created. We routinely genotype complex multiallelic CNVs in human and plant genomes using the modified MLPA procedure based on fully synthesized oligonucleotide probes (90-200 nt), which greatly simplifies the design process and allows for the development of custom assays. Here, we present a step-by-step protocol for gene-specific MLPA probe design, multiplexed assay setup and data analysis in a copy number genotyping experiment in plants. As a case study, we present the results of a custom assay designed to genotype the copy number status of 12 protein coding genes in a population of 80 Arabidopsis accessions. The genes were pre-selected based on whole genome sequencing data and are localized in the genomic regions that display different levels of population-scale variation (non-variable, biallelic, or multiallelic, as well as CNVs overlapping whole genes or their fragments). The presented approach is suitable for population-scale validation of the CNV regions inferred from whole genome sequencing data analysis and for focused analysis of selected genes of interest. It can also be very easily adopted for any plant species, following optimization of the template amount and design of the appropriate control probes, according to the general guidelines presented in this paper.

Regional heritability mapping and genome-wide association identify loci for complex growth, wood and disease resistance traits in Eucalyptus

Although genome-wide association studies (GWAS) have provided valuable insights into the decoding of the relationships between sequence variation and complex phenotypes, they have explained little heritability. Regional heritability mapping (RHM) provides heritability estimates for genomic segments containing both common and rare allelic effects that individually contribute too little variance to be detected by GWAS. We carried out GWAS and RHM for seven growth, wood and disease resistance traits in a breeding population of 768 Eucalyptus hybrid trees using EuCHIP60K. Total genomic heritabilities accounted for large proportions (64-89%) of pedigree-based trait heritabilities, providing additional evidence that complex traits in eucalypts are controlled by many sequence variants across the frequency spectrum, each with small contributions to the phenotypic variance. RHM detected 26 quantitative trait loci (QTLs) encompassing 2191 single nucleotide polymorphisms (SNPs), whereas GWAS detected 13 single SNP-trait associations. RHM and GWAS QTLs individually explained 5-15% and 4-6% of the genomic heritability, respectively. RHM was superior to GWAS in capturing larger proportions of genomic heritability. Equated to previously mapped QTLs, our results highlighted genomic regions for further examination towards gene discovery. RHM-QTLs bearing a combination of common and rare variants could be useful enhancements to incorporate prior knowledge of the underlying genetic architecture in genomic prediction models.

Dissection of Insertion-Deletion Variants within Differentially Expressed Genes Involved in Wood Formation in Populus

Short insertions and deletions (InDels) are one of the major genetic variants and are distributed widely across the genome; however, few investigations of InDels have been conducted in long-lived perennial plants. Here, we employed a combination of RNA-seq and population resequencing to identify InDels within differentially expressed (DE) genes underlying wood formation in a natural population of Populus tomentosa (435 individuals) and utilized InDel-based association mapping to detect the causal variants under additive, dominance, and epistasis underlying growth and wood properties. In the present paper, 5,482 InDels detected from 629 DE genes showed uneven distributions throughout all 19 chromosomes, and 95.9% of these loci were diallelic InDels. Seventy-four InDels (positive false discovery rate q ≤ 0.10) from 68 genes exhibited significant additive/dominant effects on 10 growth and wood-properties, with an average of 14.7% phenotypic variance explained. Potential pleiotropy was observed in one-third of the InDels (representing 24 genes). Seven genes exhibited significantly differential expression among the genotypic classes of associated InDels, indicating possible important roles for these InDels. Epistasis analysis showed that overlapping interacting genes formed unique interconnected networks for each trait, supporting the putative biochemical links that control quantitative traits. Therefore, the identification and utilization of InDels in trees will be recognized as an effective marker system for molecular marker-assisted breeding applications, and further facilitate our understanding of quantitative genomics.

Allelic variation in PtoPsbW associated with photosynthesis, growth, and wood properties in Populus tomentosa

Photosynthesis is one of the most important reactions on earth. PsbW, a nuclear-encoded subunit of photosystem II (PSII), stabilizes PSII structure and plays an important role in photosynthesis. Here, we used candidate gene-based linkage disequilibrium (LD) mapping to detect significant associations between allelic variations of PtoPsbW and traits related to photosynthesis, growth, and wood properties in Populus tomentosa. PtoPsbW showed the highest expression in leaves and it increased during the development of these leaves, suggesting that PtoPsbW may play an important role in plant growth and development. Analysis of nucleotide diversity and LD revealed that PtoPsbW has low single-nucleotide polymorphism (SNP) diversity (π _tot = 0.0048 and θ _w = 0.0050) and relatively low average value of LD (0.1500), indicating that PtoPsbW is conserved due to its indispensable function. Using single-SNP associations in an association population of 435 individuals, we identified five significant associations at the threshold of P ≤ 0.05, explaining 3.28-15.98 % of the phenotypic variation. Haplotype-based association analyses indicated that 13 haplotypes (P ≤ 0.05) from six blocks were associated with photosynthesis, growth, and wood properties. Our work shows that identifying allelic variation and LD can help to decipher the genetic basis of photosynthesis and could potentially be applied for molecular marker-assisted selection in Populus.

Genetic architecture of wood properties based on association analysis and co-expression networks in white spruce

Association studies are widely utilized to analyze complex traits but their ability to disclose genetic architectures is often limited by statistical constraints, and functional insights are usually minimal in nonmodel organisms like forest trees. We developed an approach to integrate association mapping results with co-expression networks. We tested single nucleotide polymorphisms (SNPs) in 2652 candidate genes for statistical associations with wood density, stiffness, microfibril angle and ring width in a population of 1694 white spruce trees (Picea glauca). Associations mapping identified 229-292 genes per wood trait using a statistical significance level of P < 0.05 to maximize discovery. Over-representation of genes associated for nearly all traits was found in a xylem preferential co-expression group developed in independent experiments. A xylem co-expression network was reconstructed with 180 wood associated genes and several known MYB and NAC regulators were identified as network hubs. The network revealed a link between the gene PgNAC8, wood stiffness and microfibril angle, as well as considerable within-season variation for both genetic control of wood traits and gene expression. Trait associations were distributed throughout the network suggesting complex interactions and pleiotropic effects. Our findings indicate that integration of association mapping and co-expression networks enhances our understanding of complex wood traits.

A Genome-Wide Association Study for Culm Cellulose Content in Barley Reveals Candidate Genes Co-Expressed with Members of the CELLULOSE SYNTHASE A Gene Family

Cellulose is a fundamentally important component of cell walls of higher plants. It provides a scaffold that allows the development and growth of the plant to occur in an ordered fashion. Cellulose also provides mechanical strength, which is crucial for both normal development and to enable the plant to withstand both abiotic and biotic stresses. We quantified the cellulose concentration in the culm of 288 two – rowed and 288 six – rowed spring type barley accessions that were part of the USDA funded barley Coordinated Agricultural Project (CAP) program in the USA. When the population structure of these accessions was analysed we identified six distinct populations, four of which we considered to be comprised of a sufficient number of accessions to be suitable for genome-wide association studies (GWAS). These lines had been genotyped with 3072 SNPs so we combined the trait and genetic data to carry out GWAS. The analysis allowed us to identify regions of the genome containing significant associations between molecular markers and cellulose concentration data, including one region cross-validated in multiple populations. To identify candidate genes we assembled the gene content of these regions and used these to query a comprehensive RNA-seq based gene expression atlas. This provided us with gene annotations and associated expression data across multiple tissues, which allowed us to formulate a supported list of candidate genes that regulate cellulose biosynthesis. Several regions identified by our analysis contain genes that are co-expressed with CELLULOSE SYNTHASE A (HvCesA) across a range of tissues and developmental stages. These genes are involved in both primary and secondary cell wall development. In addition, genes that have been previously linked with cellulose synthesis by biochemical methods, such as HvCOBRA, a gene of unknown function, were also associated with cellulose levels in the association panel. Our analyses provide new insights into the genes that contribute to cellulose content in cereal culms and to a greater understanding of the interactions between them.

Single nucleotide polymorphisms in a cellulose synthase gene (PtoCesA3) are associated with growth and wood properties in Populus tomentosa

In plants, the composition and organization of the cell wall determine cell shape, enable cell expansion, and affect the properties of woody tissues. Cellulose synthase (CesA) genes encode the enzymes involved in the synthesis of cellulose which is the major component of plant primary and secondary cell walls. Here, we isolated a full-length PtoCesA3 cDNA from the stem cambium tissue of Populus tomentosa. Tissue-specific expression profiling showed that PtoCesA3 is highly expressed during primary cell wall formation. Estimation of single nucleotide polymorphism (SNP) diversity and linkage disequilibrium (LD) revealed that PtoCesA3 harbors high SNP diversity (π(T) = 0.00995 and θ(w) = 0.0102) and low LD (r(2) ≥ 0.1, within 1,280 bp). Association analysis in a P. tomentosa association population (460 individuals) showed that seven SNPs (false discovery rate Q < 0.10) and five haplotypes (Q < 0.10) were significantly associated with growth and wood properties, explaining 4.09-7.02% of the phenotypic variance. All significant marker-trait associations were validated in at least one of the three smaller subsets (climatic regions) while five associations were repeated in the linkage population. Variation in RNA transcript abundance among genotypic classes of significant loci was also confirmed in the association or linkage populations. Identification of PtoCesA3 and examining its allelic polymorphisms using association studies open an avenue to understand the mechanism of cellulose synthesis in the primary cell wall and its effects on the properties of woody tissues.

Characterisation of adaptive genetic diversity in environmentally contrasted populations of Eucalyptus camaldulensis Dehnh. (river red gum)

As an increasing number of ecosystems face departures from long standing environmental conditions under climate change, our understanding of the capacity of species to adapt will become important for directing conservation and management of biodiversity. Insights into the potential for genetic adaptation might be gained by assessing genomic signatures of adaptation to historic or prevailing environmental conditions. The river red gum (Eucalyptus camaldulensis Dehnh.) is a widespread Australian eucalypt inhabiting riverine and floodplain habitats which spans strong environmental gradients. We investigated the effects of adaptation to environment on population level genetic diversity of E. camaldulensis, examining SNP variation in candidate gene loci sampled across 20 climatically diverse populations approximating the species natural distribution. Genetic differentiation among populations was high (F(ST) = 17%), exceeding previous estimates based on neutral markers. Complementary statistical approaches identified 6 SNP loci in four genes (COMT, Dehydrin, ERECTA and PIP2) which, after accounting for demographic effects, exhibited higher than expected levels of genetic differentiation among populations and whose allelic variation was associated with local environment. While this study employs but a small proportion of available diversity in the eucalyptus genome, it draws our attention to the potential for application of wide spread eucalypt species to test adaptive hypotheses.

RNA-Seq using two populations reveals genes and alleles controlling wood traits and growth in Eucalyptus nitens

Eucalyptus nitens is a perennial forest tree species grown mainly for kraft pulp production in many parts of the world. Kraft pulp yield (KPY) is a key determinant of plantation profitability and increasing the KPY of trees grown in plantations is a major breeding objective. To speed up the breeding process, molecular markers that can predict KPY are desirable. To achieve this goal, we carried out RNA-Seq studies on trees at extremes of KPY in two different trials to identify genes and alleles whose expression correlated with KPY. KPY is positively correlated with growth measured as diameter at breast height (DBH) in both trials. In total, six RNA bulks from two treatments were sequenced on an Illumina HiSeq platform. At 5% false discovery rate level, 3953 transcripts showed differential expression in the same direction in both trials; 2551 (65%) were down-regulated and 1402 (35%) were up-regulated in low KPY samples. The genes up-regulated in low KPY trees were largely involved in biotic and abiotic stress response reflecting the low growth among low KPY trees. Genes down-regulated in low KPY trees mainly belonged to gene categories involved in wood formation and growth. Differential allelic expression was observed in 2103 SNPs (in 1068 genes) and of these 640 SNPs (30%) occurred in 313 unique genes that were also differentially expressed. These SNPs may represent the cis-acting regulatory variants that influence total gene expression. In addition we also identified 196 genes which had Ka/Ks ratios greater than 1.5, suggesting that these genes are under positive selection. Candidate genes and alleles identified in this study will provide a valuable resource for future association studies aimed at identifying molecular markers for KPY and growth.

The involvement of ADAMTS-5 genetic polymorphisms in predisposition and diffusion tensor imaging alterations of lumbar disc degeneration

Low back pain is a global health problem in which more than 40% is caused by lumbar intervertebral disc degeneration (LDD). ADAMTS-5 (A disintegrin and metalloproteinase with thrombospondin motifs-5) was shown to be involved in LDD by functional analyses. To identify whether there is an association between ADAMTS-5 and LDD, and what is the contribution of ADAMTS-5 genetic polymorphisms to MD (Mean diffusivity) changes in lumbar IVD (Intervertebral disc). We firstly genotyped selected ADAMTS-5 SNPs (Single nucleotide polymorphisms) in a Chinese Han population. After the primary analyses of allelic, genotypic, and haplotypic association, we performed SNP-SNP interaction analysis. We subsequently genotyped another 50 participants and acquired the corresponding MD values from individual lumbar IVDs. The association analysis between the genotypic groups divided by the above positive SNPs and the corresponding MD values were also performed. Significant associations were identified in rs151058, rs229052, and rs162502. None of the 2-SNP haplotypic analysis survived the 10,000 permutation test. The following interaction analysis demonstrated that rs151058 was strong associated with LDD when conditioning on rs162502. Significant difference of MD values between AA and G+ carriers was identified in rs162502. This is the first study indicating that the SNPs of ADAMTS-5 may contribute to predisposition of LDD. An interaction between rs151058 and rs229052 may exist in ADAMTS-5 with LDD. The rs162502 might be associated with altered MD values.

Single-nucleotide polymorphisms in PtoCesA7 and their association with growth and wood properties in Populus tomentosa

Cellulose synthase (CesA) genes encode the enzymes that synthesize cellulose; therefore, CesAs play central roles in plant development and affect the yield and quality of wood, essential properties for industrial applications of plant biomass. To effectively manipulate wood biosynthesis in trees and improve wood quality, we thus require a better understanding of the natural variation in CesAs. Association studies have emerged as a powerful tool for identification of variation associated with quantitative traits. Here, we used a candidate gene-based association mapping approach to identify PtoCesA7 allelic variants that associate with growth and wood quality traits in Populus tomentosa. We isolated a full-length PtoCesA7 cDNA and observed high PtoCesA7 expression in xylem, consistent with the xylem-specific expression of CesA7. Nucleotide diversity and linkage disequilibrium (LD) in PtoCesA7, sampled from the P. tomentosa natural distribution, revealed that PtoCesA7 harbors high nucleotide diversity (π(T) = 0.0091) and low LD (r(2) ≥ 0.1, within 800 bp). By association analysis, we identified seven single-nucleotide polymorphisms (SNPs) (false discovery rate Q < 0.10) and 12 haplotypes (Q < 0.10) that associated with growth and wood properties, explaining 3.62-10.59 % of the phenotypic variance. We also validated 9 of the 10 significant marker-trait associations in at least one of three smaller subsets (climatic regions) or in a linkage-mapping population. Thus, our study identified functional PtoCesA7 allelic variants associated with growth and wood quality traits, giving new insights into genes affecting wood quality and quantity. From an applied perspective, the SNPs revealed in this study have potential applications in marker-assisted breeding.

Genomewide discovery of DNA polymorphisms in rice cultivars with contrasting drought and salinity stress response and their functional relevance

Next-generation sequencing technologies provide opportunities to understand the genetic basis of phenotypic differences, such as abiotic stress response, even in the closely related cultivars via identification of large number of DNA polymorphisms. We performed whole-genome resequencing of three rice cultivars with contrasting responses to drought and salinity stress (sensitive IR64, drought-tolerant Nagina 22 and salinity-tolerant Pokkali). More than 356 million 90-bp paired-end reads were generated, which provided about 85% coverage of the rice genome. Applying stringent parameters, we identified a total of 1 784 583 nonredundant single-nucleotide polymorphisms (SNPs) and 154 275 InDels between reference (Nipponbare) and the three resequenced cultivars. We detected 401 683 and 662 509 SNPs between IR64 and Pokkali, and IR64 and N22 cultivars, respectively. The distribution of DNA polymorphisms was found to be uneven across and within the rice chromosomes. One-fourth of the SNPs and InDels were detected in genic regions, and about 3.5% of the total SNPs resulted in nonsynonymous changes. Large-effect SNPs and InDels, which affect the integrity of the encoded protein, were also identified. Further, we identified DNA polymorphisms present in the differentially expressed genes within the known quantitative trait loci. Among these, a total of 548 SNPs in 232 genes, located in the conserved functional domains, were identified. The data presented in this study provide functional markers and promising target genes for salinity and drought tolerance and present a valuable resource for high-throughput genotyping and molecular breeding for abiotic stress traits in rice.

Genomics and the Evolution of Phenotypic Traits

Evolutionary genetics has entered an unprecedented era of discovery, catalyzed in large part by the development of technologies that provide information about genome sequence and function. An important benefit is the ability to move beyond a handful of model organisms in lab settings to identify the genetic basis for evolutionarily interesting traits in many organisms in natural settings. Other benefits are the abilities to identify causal mutations and validate their phenotypic consequences more readily and in many more species. Genomic technologies have reinvigorated interest in some of the most fundamental and persistent questions in evolutionary genetics, revealed previously unsuspected evolutionary phenomena, and opened the door to a wide range of new questions.

Demonstration of genome-wide association studies for identifying markers for wood property and male strobili traits in Cryptomeria japonica

Genome-wide association studies (GWAS) are an alternative to bi-parental QTL mapping in long-lived perennials. In the present study, we examined the potential of GWAS in conifers using 367 unrelated plus trees of Cryptomeria japonica D. Don, which is the most widely planted and commercially important tree species in Japan, and tried to detect significant associations between wood property traits and quantity of male strobili on the one hand, and 1,032 single nucleotide polymorphisms (SNPs) assigned to 1,032 genes on the other. Association analysis was performed with the mixed linear model taking into account kinship relationships and subpopulation structure. In total, 6 SNPs were found to have significant associations with the variations in phenotype. These SNPs were not associated with the positions of known genes and QTLs that have been reported to date, thus they may identify novel QTLs. These 6 SNPs were all found in sequences showing similarities with known genes, although further analysis is required to dissect the ways in which they affect wood property traits and abundance of male strobili. These presumptive QTL loci provide opportunities for improvement of C. japonica, based on a marker approach. The results suggest that GWAS has potential for use in future breeding programs in C. japonica.

Allelic variation in a cellulose synthase gene (PtoCesA4) associated with growth and wood properties in Populus tomentosa

Lignocellulosic biomass from trees provides a renewable feedstock for biofuels, lumber, pulp, paper, and other uses. Dissecting the mechanism underlying natural variation of the complex traits controlling growth and lignocellulose biosynthesis in trees can enable marker-assisted breeding to improve wood quality and yield. Here, we combined linkage disequilibrium (LD)-based association analysis with traditional linkage analysis to detect the genetic effect of a Populus tomentosa cellulose synthase gene, PtoCesA4. PtoCesA4 is strongly expressed in developing xylem and leaves. Nucleotide diversity and LD in PtoCesA4, sampled from the P. tomentosa natural distribution, revealed that PtoCesA4 harbors high single nucleotide polymorphism (SNP) diversity (πT = 0.0080 and θw = 0.0098) and low LD (r(2) ≥ 0.1, within 1400 bp), demonstrating that the potential of a candidate-gene-based LD approach in understanding the molecular basis underlying quantitative variation in this species. By combining single SNP, multi-SNP, and haplotype-based associations in an association population of 460 individuals with single SNP linkage analysis in a family-based linkage populations (1200 individuals), we identified three strong associations (false discovery rate Q < 0.05) in both populations. These include two nonsynonymous markers (SNP49 associated with α-cellulose content and SNP59 associated with fiber width) and a noncoding marker (SNP18 associated with α-cellulose content). Variation in RNA transcript abundance among genotypic classes of SNP49 was confirmed in these two populations. Therefore, combining different methods allowed us to examine functional PtoCesA4 allelic variation underlying natural variation in complex quantitative traits related to growth and lignocellulosic biosynthesis.

Network analysis reveals the relationship among wood properties, gene expression levels and genotypes of natural Populus trichocarpa accessions

High-throughput approaches have been widely applied to elucidate the genetic underpinnings of industrially important wood properties. Wood traits are polygenic in nature, but gene hierarchies can be assessed to identify the most important gene variants controlling specific traits within complex networks defining the overall wood phenotype. We tested a large set of genetic, genomic, and phenotypic information in an integrative approach to predict wood properties in Populus trichocarpa. Nine-yr-old natural P. trichocarpa trees including accessions with high contrasts in six traits related to wood chemistry and ultrastructure were profiled for gene expression on 49k Nimblegen (Roche NimbleGen Inc., Madison, WI, USA) array elements and for 28,831 polymorphic single nucleotide polymorphisms (SNPs). Pre-selected transcripts and SNPs with high statistical dependence on phenotypic traits were used in Bayesian network learning procedures with a stepwise K2 algorithm to infer phenotype-centric networks. Transcripts were pre-selected at a much lower logarithm of Bayes factor (logBF) threshold than SNPs and were not accommodated in the networks. Using persistent variables, we constructed cross-validated networks for variability in wood attributes, which contained four to six variables with 94-100% predictive accuracy. Accommodated gene variants revealed the hierarchy in the genetic architecture that underpins substantial phenotypic variability, and represent new tools to support the maximization of response to selection.

Genome-wide association mapping for wood characteristics in Populus identifies an array of candidate single nucleotide polymorphisms

Establishing links between phenotypes and molecular variants is of central importance to accelerate genetic improvement of economically important plant species. Our work represents the first genome-wide association study to the inherently complex and currently poorly understood genetic architecture of industrially relevant wood traits. Here, we employed an Illumina Infinium 34K single nucleotide polymorphism (SNP) genotyping array that generated 29,233 high-quality SNPs in c. 3500 broad-based candidate genes within a population of 334 unrelated Populus trichocarpa individuals to establish genome-wide associations. The analysis revealed 141 significant SNPs (α ≤ 0.05) associated with 16 wood chemistry/ultrastructure traits, individually explaining 3-7% of the phenotypic variance. A large set of associations (41% of all hits) occurred in candidate genes preselected for their suggested a priori involvement with secondary growth. For example, an allelic variant in the FRA8 ortholog explained 21% of the total genetic variance in fiber length, when the trait's heritability estimate was considered. The remaining associations identified SNPs in genes not previously implicated in wood or secondary wall formation. Our findings provide unique insights into wood trait architecture and support efforts for population improvement based on desirable allelic variants.

Expression of a conifer COBRA-like gene ClCOBL1 from Chinese fir (Cunninghamia lanceolata) alters the leaf architecture in tobacco

The cell wall plays crucial roles in establishing the morphology of the plant cell, defence response to biotic and abiotic stresses, and mechanical properties of organs. The COBRA gene encodes a putative glycosylphosphatidylinositol (GPI)-anchored protein that possesses the ability to modulate cellulose deposition and orient cell expansion in plant cell. We reported here the functional characterization of ClCOBL1, a conifer COBRA-like gene from the differentiating xylem of Chinese fir (Cunninghamia lanceolata (Lamb.) Hook). ClCOBL1 belonged to a woody plant-specific clade of the COBRA protein family with several conserved motifs. Expression pattern demonstrated that ClCOBL1 was constitutively expressed but with high level in cambium region. ClCOBL1 protein was mainly located in the cell wall and plasma membrane. Overexpression of ClCOBL1 in tobacco plants yielded altered leaf adaxial-abaxial patterning and short, swollen corolla tubes. The changed leaf architecture in the ClCOBL1 overexpressors was associated with the differential expression of leaf adaxial-abaxial identity genes. Our results indicated that ClCOBL1 was involved in the determination of leaf dorsoventrality and anisotropic expansion possibly by affecting the expression of adaxial and abaxial identity genes.

The UDP-glucuronate decarboxylase gene family in Populus: structure, expression, and association genetics

In woody crop plants, the oligosaccharide components of the cell wall are essential for important traits such as bioenergy content, growth, and structural wood properties. UDP-glucuronate decarboxylase (UXS) is a key enzyme in the synthesis of UDP-xylose for the formation of xylans during cell wall biosynthesis. Here, we isolated a multigene family of seven members (PtUXS1-7) encoding UXS from Populus tomentosa, the first investigation of UXSs in a tree species. Analysis of gene structure and phylogeny showed that the PtUXS family could be divided into three groups (PtUXS1/4, PtUXS2/5, and PtUXS3/6/7), consistent with the tissue-specific expression patterns of each PtUXS. We further evaluated the functional consequences of nucleotide polymorphisms in PtUXS1. In total, 243 single-nucleotide polymorphisms (SNPs) were identified, with a high frequency of SNPs (1/18 bp) and nucleotide diversity (π_T = 0.01033, θ_w =0.01280). Linkage disequilibrium (LD) analysis showed that LD did not extend over the entire gene (r²<0.1, P<0.001, within 700 bp). SNP- and haplotype-based association analysis showed that nine SNPs (Q <0.10) and 12 haplotypes (P<0.05) were significantly associated with growth and wood property traits in the association population (426 individuals), with 2.70% to 12.37% of the phenotypic variation explained. Four significant single-marker associations (Q <0.10) were validated in a linkage mapping population of 1200 individuals. Also, RNA transcript accumulation varies among genotypic classes of SNP10 was further confirmed in the association population. This is the first comprehensive study of the UXS gene family in woody plants, and lays the foundation for genetic improvements of wood properties and growth in trees using genetic engineering or marker-assisted breeding.

Epigenetic regulation of adaptive responses of forest tree species to the environment

Epigenetic variation is likely to contribute to the phenotypic plasticity and adaptative capacity of plant species, and may be especially important for long-lived organisms with complex life cycles, including forest trees. Diverse environmental stresses and hybridization/polyploidization events can create reversible heritable epigenetic marks that can be transmitted to subsequent generations as a form of molecular "memory". Epigenetic changes might also contribute to the ability of plants to colonize or persist in variable environments. In this review, we provide an overview of recent data on epigenetic mechanisms involved in developmental processes and responses to environmental cues in plant, with a focus on forest tree species. We consider the possible role of forest tree epigenetics as a new source of adaptive traits in plant breeding, biotechnology, and ecosystem conservation under rapid climate change.

Polymorphic simple sequence repeat (SSR) loci within cellulose synthase (PtoCesA) genes are associated with growth and wood properties in Populus tomentosa

Chinese white poplar (Populus tomentosa), an important commercial tree species for timber and pulp production in northern China, has been used to examine the individual genes and allelic diversity responsible for complex traits controlling growth and lignocellulosic biosynthesis. Taking advantage of the low degree of linkage disequilibrium (LD) within P. tomentosa association populations, we examined associations between 15 cellulose synthase (PtoCesA) genes and traits including growth and wood properties. Thirty-six novel simple sequence repeat (SSR) markers within PtoCesA genes were detected by re-sequencing and genotyped in an association population (460 individuals). Single-marker and haplotype-based LD approaches were used to identify significant marker-trait associations. Family-based linkage studies and real-time PCR testing were conducted to validate the functional significance of SSR variation. Fifteen single-marker associations from seven PtoCesA genes and nine haplotype-based associations within six genes were identified in the association population (false discovery rate Q < 0.05). Next, five SSR marker-trait associations (Q < 0.05) from four PtoCesA genes were successfully validated in a linkage mapping population (1200 individuals). The results imply a functional role for these genes in mediating wood properties, demonstrating the potential of combining single-marker and haplotype-based LD approaches to detect functional allelic variation underlying quantitative traits in a low-LD population.

Allelic variation in PtGA20Ox associates with growth and wood properties in Populus spp

Populus tomentosa is an economically important tree crop that produces wood for lumber, pulp, paper, and biofuels. Wood quality traits are likely to be strongly affected by the plant hormone gibberellic acid (GA), which regulates growth. GA20Ox encodes one of the major regulatory enzymes of GA biosynthesis and may therefore play a large role in growth and wood quality. Here, linkage disequilibrium (LD) studies were used to identify significant associations between single nucleotide polymorphisms (SNPs) within PtGA20Ox and growth and wood-quality traits of P. tomentosa. We isolated a full-length GA20Ox cDNA from Populus tomentosa by reverse transcription (RT)-PCR; this 1401 bp cDNA clone had an open reading frame of 1158 bp and encoded a protein of 385 amino acids. PtGA20Ox transcripts were maximally expressed in the mature xylem of vascular tissues, suggesting that PtGA20Ox is highly expressed and specifically associated with secondary xylem formation. Resequencing the PtGA20Ox locus of 36 individuals identified 55 SNPs, and the frequency of SNPs was 1/31 bp. The 29 most common SNPs (frequency>0.1) were genotyped in an association population (426 individuals) that was also phenotyped for key growth and wood quality traits. LD did not extend over the entire gene (r(2)<0.1, within 500 bp), demonstrating that a candidate-gene-based LD approach may the best way to understand the molecular basis underlying quantitative variation in this species. SNP- and haplotype-based association analyses indicated that four SNPs (false discovery rate Q<0.05) and 14 haplotypes (P<0.05) were significantly associated with growth and wood properties. The phenotypic variance explained by each SNP ranged from 3.44% to 14.47%. The SNP markers identified in this study can be applied to breeding programs for the improvement of growth and wood-property traits by marker-assisted selection.

High-throughput and parallel SNP discovery in selected candidate genes in Eucalyptus camaldulensis using Illumina NGS platform

Next generation sequencing (NGS) technologies have revolutionized the pace and scale of genomics- and transcriptomics-based SNP discovery across different plant and animal species. Herein, 72-base paired-end Illumina sequencing was employed for high-throughput, parallel and large-scale SNP discovery in 41 growth-related candidate genes in Eucalyptus camaldulensis. Approximately 100 kb of genome from 96 individuals was amplified and sequenced using a hierarchical DNA/PCR pooling strategy and assembled over corresponding E. grandis reference. A total of 1191 SNPs (minimum 5% other allele frequency) were identified with an average frequency of 1 SNP/83.9 bp, whereas in exons and introns, it was 1 SNP/108.4 bp and 1 SNP/65.6 bp, respectively. A total of 75 insertions and 89 deletions were detected of which approximately 15% were exonic. Transitions (Tr) were in excess than transversions (Tv) (Tr/Tv: 1.89), but exceeded in exons (Tr/Tv: 2.73). In exons, synonymous SNPs (Ka) prevailed over the non-synonymous SNPs (Ks; average Ka/Ks ratio: 0.72, range: 0-3.00 across genes). Many of the exonic SNPs/indels had potential to change amino acid sequence of respective genes. Transcription factors appeared more conserved, whereas enzyme coding genes appeared under relaxed control. Further, 541 SNPs were classified into 196 'equal frequency' (EF) blocks with almost similar minor allele frequencies to facilitate selection of one tag-SNP/EF-block. There were 241 (approximately 20%) 'zero-SNP' blocks with absence of SNPs in surrounding ±60 bp windows. The data thus indicated enormous extant and unexplored diversity in E. camaldulensis in the studied genes with potential applications for marker-trait associations.

Transcriptome sequencing of Eucalyptus camaldulensis seedlings subjected to water stress reveals functional single nucleotide polymorphisms and genes under selection

Background

Water stress limits plant survival and production in many parts of the world. Identification of genes and alleles responding to water stress conditions is important in breeding plants better adapted to drought. Currently there are no studies examining the transcriptome wide gene and allelic expression patterns under water stress conditions. We used RNA sequencing (RNA-seq) to identify the candidate genes and alleles and to explore the evolutionary signatures of selection.

Results

We studied the effect of water stress on gene expression in Eucalyptus camaldulensis seedlings derived from three natural populations. We used reference-guided transcriptome mapping to study gene expression. Several genes showed differential expression between control and stress conditions. Gene ontology (GO) enrichment tests revealed up-regulation of 140 stress-related gene categories and down-regulation of 35 metabolic and cell wall organisation gene categories. More than 190,000 single nucleotide polymorphisms (SNPs) were detected and 2737 of these showed differential allelic expression. Allelic expression of 52% of these variants was correlated with differential gene expression. Signatures of selection patterns were studied by estimating the proportion of nonsynonymous to synonymous substitution rates (Ka/Ks). The average Ka/Ks ratio among the 13,719 genes was 0.39 indicating that most of the genes are under purifying selection. Among the positively selected genes (Ka/Ks > 1.5) apoptosis and cell death categories were enriched. Of the 287 positively selected genes, ninety genes showed differential expression and 27 SNPs from 17 positively selected genes showed differential allelic expression between treatments.

Conclusions

Correlation of allelic expression of several SNPs with total gene expression indicates that these variants may be the cis-acting variants or in linkage disequilibrium with such variants. Enrichment of apoptosis and cell death gene categories among the positively selected genes reveals the past selection pressures experienced by the populations used in this study.

Association genetics in Corymbia citriodora subsp. variegata identifies single nucleotide polymorphisms affecting wood growth and cellulosic pulp yield

Wood is an important biological resource which contributes to nutrient and hydrology cycles through ecosystems, and provides structural support at the plant level. Thousands of genes are involved in wood development, yet their effects on phenotype are not well understood. We have exploited the low genomic linkage disequilibrium (LD) and abundant phenotypic variation of forest trees to explore allelic diversity underlying wood traits in an association study. Candidate gene allelic diversity was modelled against quantitative variation to identify SNPs influencing wood properties, growth and disease resistance across three populations of Corymbia citriodora subsp. variegata, a forest tree of eastern Australia. Nine single nucleotide polymorphism (SNP) associations from six genes were identified in a discovery population (833 individuals). Associations were subsequently tested in two smaller populations (130-160 individuals), 'validating' our findings in three cases for actin 7 (ACT7) and COP1 interacting protein 7 (CIP7). The results imply a functional role for these genes in mediating wood chemical composition and growth, respectively. A flip in the effect of ACT7 on pulp yield between populations suggests gene by environment interactions are at play. Existing evidence of gene function lends strength to the observed associations, and in the case of CIP7 supports a role in cortical photosynthesis.

Association mapping in forest trees and fruit crops

Association mapping (AM), also known as linkage disequilibrium (LD) mapping, is a viable approach to overcome limitations of pedigree-based quantitative trait loci (QTL) mapping. In AM, genotypic and phenotypic correlations are investigated in unrelated individuals. Unlike QTL mapping, AM takes advantage of both LD and historical recombination present within the gene pool of an organism, thus utilizing a broader reference population. In plants, AM has been used in model species with available genomic resources. Pursuing AM in tree species requires both genotyping and phenotyping of large populations with unique architectures. Recently, genome sequences and genomic resources for forest and fruit crops have become available. Due to abundance of single nucleotide polymorphisms (SNPs) within a genome, along with availability of high-throughput resequencing methods, SNPs can be effectively used for genotyping trees. In addition to DNA polymorphisms, copy number variations (CNVs) in the form of deletions, duplications, and insertions also play major roles in control of expression of phenotypic traits. Thus, CNVs could provide yet another valuable resource, beyond those of microsatellite and SNP variations, for pursuing genomic studies. As genome-wide SNP data are generated from high-throughput sequencing efforts, these could be readily reanalysed to identify CNVs, and subsequently used for AM studies. However, forest and fruit crops possess unique architectural and biological features that ought to be taken into consideration when collecting genotyping and phenotyping data, as these will also dictate which AM strategies should be pursued. These unique features as well as their impact on undertaking AM studies are outlined and discussed.

Genomic selection for growth and wood quality in Eucalyptus: capturing the missing heritability and accelerating breeding for complex traits in forest trees

• Genomic selection (GS) is expected to cause a paradigm shift in tree breeding by improving its speed and efficiency. By fitting all the genome-wide markers concurrently, GS can capture most of the 'missing heritability' of complex traits that quantitative trait locus (QTL) and association mapping classically fail to explain. Experimental support of GS is now required. • The effectiveness of GS was assessed in two unrelated Eucalyptus breeding populations with contrasting effective population sizes (N(e) = 11 and 51) genotyped with > 3000 DArT markers. Prediction models were developed for tree circumference and height growth, wood specific gravity and pulp yield using random regression best linear unbiased predictor (BLUP). • Accuracies of GS varied between 0.55 and 0.88, matching the accuracies achieved by conventional phenotypic selection. Substantial proportions (74-97%) of trait heritability were captured by fitting all genome-wide markers simultaneously. Genomic regions explaining trait variation largely coincided between populations, although GS models predicted poorly across populations, likely as a result of variable patterns of linkage disequilibrium, inconsistent allelic effects and genotype × environment interaction. • GS brings a new perspective to the understanding of quantitative trait variation in forest trees and provides a revolutionary tool for applied tree improvement. Nevertheless population-specific predictive models will likely drive the initial applications of GS in forest tree breeding.

Cellulose factories: advancing bioenergy production from forest trees

Fast-growing, short-rotation forest trees, such as Populus and Eucalyptus, produce large amounts of cellulose-rich biomass that could be utilized for bioenergy and biopolymer production. Major obstacles need to be overcome before the deployment of these genera as energy crops, including the effective removal of lignin and the subsequent liberation of carbohydrate constituents from wood cell walls. However, significant opportunities exist to both select for and engineer the structure and interaction of cell wall biopolymers, which could afford a means to improve processing and product development. The molecular underpinnings and regulation of cell wall carbohydrate biosynthesis are rapidly being elucidated, and are providing tools to strategically develop and guide the targeted modification required to adapt forest trees for the emerging bioeconomy. Much insight has already been gained from the perturbation of individual genes and pathways, but it is not known to what extent the natural variation in the sequence and expression of these same genes underlies the inherent variation in wood properties of field-grown trees. The integration of data from next-generation genomic technologies applied in natural and experimental populations will enable a systems genetics approach to study cell wall carbohydrate production in trees, and should advance the development of future woody bioenergy and biopolymer crops.

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.

Pectin Methylesterase genes influence solid wood properties of Eucalyptus pilularis

This association study of Eucalyptus pilularis populations provides empirical evidence for the role of Pectin Methylesterase (PME) in influencing solid wood characteristics of Eucalyptus. PME6 was primarily associated with the shrinkage and collapse of drying timber, which are phenotypic traits consistent with the role of pectin as a hydrophilic polysaccharide. PME7 was primarily associated with cellulose and pulp yield traits and had an inverse correlation with lignin content. Selection of specific alleles in these genes may be important for improving trees as sources of high-quality wood products. A heterozygote advantage was postulated for the PME7 loci and, in combination with haplotype blocks, may explain the absence of a homozygous class at all single-nucleotide polymorphisms investigated in this gene.

Association genetics of wood physical traits in the conifer white spruce and relationships with gene expression

Marker-assisted selection holds promise for highly influencing tree breeding, especially for wood traits, by considerably reducing breeding cycles and increasing selection accuracy. In this study, we used a candidate gene approach to test for associations between 944 single-nucleotide polymorphism markers from 549 candidate genes and 25 wood quality traits in white spruce. A mixed-linear model approach, including a weak but nonsignificant population structure, was implemented for each marker-trait combination. Relatedness among individuals was controlled using a kinship matrix estimated either from the known half-sib structure or from the markers. Both additive and dominance effect models were tested. Between 8 and 21 single-nucleotide polymorphisms (SNPs) were found to be significantly associated (P ≤ 0.01) with each of earlywood, latewood, or total wood traits. After controlling for multiple testing (Q ≤ 0.10), 13 SNPs were still significant across as many genes belonging to different families, each accounting for between 3 and 5% of the phenotypic variance in 10 wood characters. Transcript accumulation was determined for genes containing SNPs associated with these traits. Significantly different transcript levels (P ≤ 0.05) were found among the SNP genotypes of a 1-aminocyclopropane-1-carboxylate oxidase, a β-tonoplast intrinsic protein, and a long-chain acyl-CoA synthetase 9. These results should contribute toward the development of efficient marker-assisted selection in an economically important tree 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.