QTL analysis of leaf morphological characters in a Quercus robur full-sib family (Q. robur x Q. robur ssp. slavonica)

Genome Wide Association Study Identifies Candidate Genes Related to the Earlywood Tracheid Properties in Picea crassifolia Kom

Picea crassifolia Kom. is one of the timber and ecological conifers in China and its wood tracheid traits directly affect wood formation and adaptability under harsh environment. Molecular studies on P. crassifolia remain inadequate because relatively few genes have been associated with these traits. To identify markers and candidate genes that can potentially be used for genetic improvement of wood tracheid traits, we examined 106 clones of P. crassifolia, and investigated phenotypic data for 14 wood tracheid traits before specific-locus amplified fragment sequencing (SLAF-seq) was employed to perform a genome wide association study (GWAS). Subsequently, the results were used to screen single nucleotide polymorphism (SNP) loci and candidate genes that exhibited a significant correlation with the studied traits. We developed 4,058,883 SLAF-tags and 12,275,765 SNP loci, and our analyses identified a total of 96 SNP loci that showed significant correlations with three earlywood tracheid traits using a mixed linear model (MLM). Next, candidate genes were screened in the 100 kb zone (50 kb upstream, 50 kb downstream) of each of the SNP loci, whereby 67 candidate genes were obtained in earlywood tracheid traits, including 34 genes of known function and 33 genes of unknown function. We provide the most significant SNP for each trait-locus combination and candidate genes occurring within the GWAS hits. These resources provide a foundation for the development of markers that could be used in wood traits improvement and candidate genes for the development of earlywood tracheid in P. crassifolia.

Quercus species divergence is driven by natural selection on evolutionarily less integrated traits

Functional traits are organismal attributes that can respond to environmental cues, thereby providing important ecological functions. In addition, an organism's potential for adaptation is defined by the patterns of covariation among groups of functionally related traits. Whether an organism is evolutionarily constrained or has the potential for adaptation is based on the phenotypic integration or modularity of these traits. Here, we revisited leaf morphology in two European sympatric white oaks (Quercus petraea (Matt.) Liebl. and Quercus robur L.), sampling 2098 individuals, across much of their geographical distribution ranges. At the phenotypic level, leaf morphology traditionally encompasses discriminant attributes among different oak species. Here, we estimated in situ heritability, genetic correlation, and integration across such attributes. Also, we performed Selection Response Decomposition to test these traits for potential differences in oak species' evolutionary responses. Based on the uncovered functional units of traits (modules) in our study, the morphological module "leaf size gradient" was highlighted among functionally integrated traits. Equally, this module was defined in both oaks as being under "global regulation" in vegetative bud establishment and development. Lamina basal shape and intercalary veins' number were not, or, less integrated within the initially defined leaf functional unit, suggesting more than one module within the leaf traits' ensemble. Since these traits generally show the greatest species discriminatory power, they potentially underwent effective differential response to selection among oaks. Indeed, the selection of these traits could have driven the ecological preferences between the two sympatric oaks growing under different microclimates.

Terpene Synthase Genes in Quercus robur - Gene Characterization, Expression and Resulting Terpenes Due to Cockchafer Feeding

Root herbivory caused by larvae of the forest cockchafer (Melolontha hippocastani) enhances the impact of drought on trees, particularly in oak forest rejuvenations. In Germany, geographically distant oak stands show differences in infestation strength by the forest cockchafer. While in Southwestern Germany this insect causes severe damage, oak forests in northern Germany are rarely infested. It is known that root-released volatile organic compounds (VOCs) are perceived by soil herbivores, thus guiding the larvae toward the host roots. In this work, we exposed seedlings of two distant oak provenances to forest cockchafer larvae and studied their population genetic properties, their root-based VOC chemotypes, their attraction for larvae and terpene synthase gene expression. Based on nuclear and chloroplast marker analysis, we found both oak populations to be genetically highly variable while showing typical patterns of migration from different refugial regions. However, no clear association between genetic constitution of the different provenances and the abundance of cockchafer populations on site was observed. In contrast to observations in the field, bioassays revealed a preference of the larvae for the northeastern oak provenance. The behavior of larvae was most likely related to root-released volatile terpenes and benzenoids since their composition and quantity differed between oak populations. We assume repellent effects of these compounds because the populations attractive to insects showed low abundance of these compounds. Five different oak terpene synthase (TPS) genes were identified at the genomic level which can be responsible for biosynthesis of the released terpenes. TPS gene expression patterns in response to larval feeding revealed geographic variation rather than genotypic variation. Our results support the assumption that root-released VOC are influencing the perception of roots by herbivores.

Oak genome reveals facets of long lifespan

Oaks are an important part of our natural and cultural heritage. Not only are they ubiquitous in our most common landscapes¹ but they have also supplied human societies with invaluable services, including food and shelter, since prehistoric times². With 450 species spread throughout Asia, Europe and America³, oaks constitute a critical global renewable resource. The longevity of oaks (several hundred years) probably underlies their emblematic cultural and historical importance. Such long-lived sessile organisms must persist in the face of a wide range of abiotic and biotic threats over their lifespans. We investigated the genomic features associated with such a long lifespan by sequencing, assembling and annotating the oak genome. We then used the growing number of whole-genome sequences for plants (including tree and herbaceous species) to investigate the parallel evolution of genomic characteristics potentially underpinning tree longevity. A further consequence of the long lifespan of trees is their accumulation of somatic mutations during mitotic divisions of stem cells present in the shoot apical meristems. Empirical⁴ and modelling⁵ approaches have shown that intra-organismal genetic heterogeneity can be selected for⁶ and provides direct fitness benefits in the arms race with short-lived pests and pathogens through a patchwork of intra-organismal phenotypes⁷. However, there is no clear proof that large-statured trees consist of a genetic mosaic of clonally distinct cell lineages within and between branches. Through this case study of oak, we demonstrate the accumulation and transmission of somatic mutations and the expansion of disease-resistance gene families in trees.

Leaf morphological variability and intraspecific taxonomic units for pedunculate oak and grayish oak (genus Quercus L., series Pedunculatae Schwz.) in Southern Carpathian Region (Romania)

Even though pedunculate oak (Quercus robur L.) and grayish oak (Quercus pedunculiflora K. Koch) have different ecological requirements, they have been considered as having low differentiation at the level of morphological traits and genetic variation. The leaf morphology for 862 trees has been assessed in 16 natural populations, seven of Q. robur, eight of Q. pedunculiflora and a mixed forest were both taxa coexist. In total, fifteen descriptors have been analysed by using discriminant analysis, while it was found that with only four out of the fifteen leaf traits (abaxial pubescence, abaxial colour of the leaf, petiole length and basal shape of lamina) the two taxa could be clearly differentiated. A dendrogram has been constructed on the basis of these traits, where the populations of each taxon have been clustered together. PU and CL traits of Q. pedunculiflora were discussed for their adaptive value for drought resistance in the steppe habitats occupied by this taxon. Using the leaves' morphological descriptors and data from the literature, intra-taxonomic units (varieties, forms and sub-forms) have been identified in all analysed populations. Eight intraspecific units for Q. robur and six for Q. pedunculiflora have been identified in the investigated area. An analysis of spatial distribution of the two taxa and of their intraspecific units has been performed using maps of ecoregions for the study area.

Screening of Candidate Leaf Morphology Genes by Integration of QTL Mapping and RNA Sequencing Technologies in Oilseed Rape (Brassica napus L.)

Leaf size and shape play important roles in agronomic traits, such as yield, quality and stress responses. Wide variations in leaf morphological traits exist in cultivated varieties of many plant species. By now, the genetics of leaf shape and size have not been characterized in Brassica napus. In this study, a population of 172 recombinant inbred lines (RILs) was used for quantitative trait locus (QTL) analysis of leaf morphology traits. Furthermore, fresh young leaves of extreme lines with more leaf lobes (referred to as ‘A’) and extreme lines with fewer lobes (referred to as ‘B’) selected from the RIL population and leaves of dissected lines (referred to as ‘P’) were used for transcriptional analysis. A total of 31 QTLs for the leaf morphological traits tested in this study were identified on 12 chromosomes, explaining 5.32–39.34% of the phenotypic variation. There were 8, 6, 2, 5, 8, and 2 QTLs for PL (petiole length), PN (lobe number), LW (lamina width), LL (Lamina length), LL/LTL (the lamina size ratio) and LTL (leaf total length), respectively. In addition, 74, 1,166 and 1,272 differentially expressed genes (DEGs) were identified in ‘A vs B’, ‘A vs P’ and ‘B vs P’ comparisons, respectively. The Gene ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) databases were used to predict the functions of these DEGs. Gene regulators of leaf shape and size, such as ASYMMETRIC LEAVES 2, gibberellin 20-oxidase 3, genes encoding gibberellin-regulated family protein, genes encoding growth-regulating factor and KNOTTED1-like homeobox were also detected in DEGs. After integrating the QTL mapping and RNA sequencing data, 33 genes, including a gene encoding auxin-responsive GH3 family protein and a gene encoding sphere organelles protein-related gene, were selected as candidates that may control leaf shape. Our findings should be valuable for studies of the genetic control of leaf morphological trait regulation in B. napus.

Limited Pollen Dispersal Contributes to Population Genetic Structure but Not Local Adaptation in Quercus oleoides Forests of Costa Rica

Background

Quercus oleoides Cham. and Schlect., tropical live oak, is a species of conservation importance in its southern range limit of northwestern Costa Rica. It occurs in high-density stands across a fragmented landscape spanning a contrasting elevation and precipitation gradient. We examined genetic diversity and spatial genetic structure in this geographically isolated and genetically distinct population. We characterized population genetic diversity at 11 nuclear microsatellite loci in 260 individuals from 13 sites. We monitored flowering time at 10 sites, and characterized the local environment in order to compare observed spatial genetic structure to hypotheses of isolation-by-distance and isolation-by-environment. Finally, we quantified pollen dispersal distances and tested for local adaptation through a reciprocal transplant experiment in order to experimentally address these hypotheses.

Results

High genetic diversity is maintained in the population and the genetic variation is significantly structured among sampled sites. We identified 5 distinct genetic clusters and average pollen dispersal predominately occurred over short distances. Differences among sites in flowering phenology and environmental factors, however, were not strictly associated with genetic differentiation. Growth and survival of upland and lowland progeny in their native and foreign environments was expected to exhibit evidence of local adaptation due to the more extreme dry season in the lowlands. Seedlings planted in the lowland garden experienced much higher mortality than seedlings in the upland garden, but we did not identify evidence for local adaptation.

Conclusion

Overall, this study indicates that the Costa Rican Q. oleoides population has a rich population genetic history. Despite environmental heterogeneity and habitat fragmentation, isolation-by-distance and isolation-by-environment alone do not explain spatial genetic structure. These results add to studies of genetic structure by examining a common, tropical tree over multiple habitats and provide information for managers of a successional forest in a protected area.

Genetical genomics of Populus leaf shape variation

BACKGROUND

Leaf morphology varies extensively among plant species and is under strong genetic control. Mutagenic screens in model systems have identified genes and established molecular mechanisms regulating leaf initiation, development, and shape. However, it is not known whether this diversity across plant species is related to naturally occurring variation at these genes. Quantitative trait locus (QTL) analysis has revealed a polygenic control for leaf shape variation in different species suggesting that loci discovered by mutagenesis may only explain part of the naturally occurring variation in leaf shape. Here we undertook a genetical genomics study in a poplar intersectional pseudo-backcross pedigree to identify genetic factors controlling leaf shape. The approach combined QTL discovery in a genetic linkage map anchored to the Populus trichocarpa reference genome sequence and transcriptome analysis.

RESULTS

A major QTL for leaf lamina width and length:width ratio was identified in multiple experiments that confirmed its stability. A transcriptome analysis of expanding leaf tissue contrasted gene expression between individuals with alternative QTL alleles, and identified an ADP-ribosylation factor (ARF) GTPase (PtARF1) as a candidate gene for regulating leaf morphology in this pedigree. ARF GTPases are critical elements in the vesicular trafficking machinery. Disruption of the vesicular trafficking function of ARF by the pharmacological agent Brefeldin A (BFA) altered leaf lateral growth in the narrow-leaf P. trichocarpa suggesting a molecular mechanism of leaf shape determination. Inhibition of the vesicular trafficking processes by BFA interferes with cycling of PIN proteins and causes their accumulation in intercellular compartments abolishing polar localization and disrupting normal auxin flux with potential effects on leaf expansion.

CONCLUSIONS

In other model systems, ARF proteins have been shown to control the localization of auxin efflux carriers, which function to establish auxin gradients and apical-basal cell polarity in developing plant organs. Our results support a model where PtARF1 transcript abundance changes the dynamics of endocytosis-mediated PIN localization in leaf cells, thus affecting lateral auxin flux and subsequently lamina leaf expansion. This suggests that evolution of differential cellular polarity plays a significant role in leaf morphological variation observed in subgenera of genus Populus.

Genetic and genomic approaches to assess adaptive genetic variation in plants: forest trees as a model

With the increasing availability of sequence information at putatively important genes or regulatory regions, the characterization of adaptive genetic diversity and their association with phenotypic trait variation becomes feasible for many non-model organisms such as forest trees. Especially in predominantly outcrossing forest tree populations with large effective size, a high genetic variation in relevant genes is maintained, that is the raw material for the adaptation to changing and variable environments, and likewise for plant breeding. Oaks (Quercus spp.) are excellent model species to study the adaptation of forest trees to changing environments. They show a wide geographic distribution in Europe as dominant tree species in many forests and grow under a wide range of climatic and edaphic conditions. With the availability of a growing amount of functional and expressional candidate genes, we are now able to test the functional importance of single nucleotide polymorphisms (SNPs) by associating nucleotide variation in these genes with phenotypic variation in adaptive traits in segregating or natural populations. Here, we report on quantitative trait locus (QTL), candidate gene and association mapping approaches that are applicable to characterize gene markers and SNPs associated with variation in adaptive traits, such as bud burst, drought resistance and other traits showing selective responses to environmental change and stress. Because genome-wide association mapping studies are not feasible because of the enormous amount of SNP markers required in outcrossing trees with high recombination rates, the success of such an approach depends largely on the reasonable selection of candidate genes.