High-quality genetic mapping with ddRADseq in the non-model tree Quercus rubra

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.

A haplotype-resolved chromosome-scale genome for Quercus rubra L. provides insights into the genetics of adaptive traits for red oak species

Northern red oak (Quercus rubra L.) is an ecologically and economically important forest tree native to North America. We present a chromosome-scale genome of Q. rubra generated by the combination of PacBio sequences and chromatin conformation capture (Hi-C) scaffolding. This is the first reference genome from the red oak clade (section Lobatae). The Q. rubra assembly spans 739 Mb with 95.27% of the genome in 12 chromosomes and 33,333 protein-coding genes. Comparisons to the genomes of Quercus lobata and Quercus mongolica revealed high collinearity, with intrachromosomal structural variants present. Orthologous gene family analysis with other tree species revealed that gene families associated with defense response were expanding and contracting simultaneously across the Q. rubra genome. Quercus rubra had the most CC-NBS-LRR and TIR-NBS-LRR resistance genes out of the 9 species analyzed. Terpene synthase gene family comparisons further reveal tandem gene duplications in TPS-b subfamily, similar to Quercus robur. Phylogenetic analysis also identified 4 subfamilies of the IGT/LAZY gene family in Q. rubra important for plant structure. Single major QTL regions were identified for vegetative bud break and marcescence, which contain candidate genes for further research, including a putative ortholog of the circadian clock constituent cryptochrome (CRY2) and 8 tandemly duplicated genes for serine protease inhibitors, respectively. Genome-environment associations across natural populations identified candidate abiotic stress tolerance genes and predicted performance in a common garden. This high-quality red oak genome represents an essential resource to the oak genomic community, which will expedite comparative genomics and biological studies in Quercus species.

Key triggers of adaptive genetic variability of sessile oak [Q. petraea (Matt.) Liebl.] from the Balkan refugia: outlier detection and association of SNP loci from ddRAD-seq data

Knowledge on the genetic composition of Quercus petraea in south-eastern Europe is limited despite the species' significant role in the re-colonisation of Europe during the Holocene, and the diverse climate and physical geography of the region. Therefore, it is imperative to conduct research on adaptation in sessile oak to better understand its ecological significance in the region. While large sets of SNPs have been developed for the species, there is a continued need for smaller sets of SNPs that are highly informative about the possible adaptation to this varied landscape. By using double digest restriction site associated DNA sequencing data from our previous study, we mapped RAD-seq loci to the Quercus robur reference genome and identified a set of SNPs putatively related to drought stress-response. A total of 179 individuals from eighteen natural populations at sites covering heterogeneous climatic conditions in the southeastern natural distribution range of Q. petraea were genotyped. The detected highly polymorphic variant sites revealed three genetic clusters with a generally low level of genetic differentiation and balanced diversity among them but showed a north-southeast gradient. Selection tests showed nine outlier SNPs positioned in different functional regions. Genotype-environment association analysis of these markers yielded a total of 53 significant associations, explaining 2.4-16.6% of the total genetic variation. Our work exemplifies that adaptation to drought may be under natural selection in the examined Q. petraea populations.

Applying molecular and genetic methods to trees and their fungal communities

Forests provide invaluable economic, ecological, and social services. At the same time, they are exposed to several threats, such as fragmentation, changing climatic conditions, or increasingly destructive pests and pathogens. Trees, the inherent species of forests, cannot be viewed as isolated organisms. Manifold (micro)organisms are associated with trees playing a pivotal role in forest ecosystems. Of these organisms, fungi may have the greatest impact on the life of trees. A multitude of molecular and genetic methods are now available to investigate tree species and their associated organisms. Due to their smaller genome sizes compared to tree species, whole genomes of different fungi are routinely compared. Such studies have only recently started in forest tree species. Here, we summarize the application of molecular and genetic methods in forest conservation genetics, tree breeding, and association genetics as well as for the investigation of fungal communities and their interrelated ecological functions. These techniques provide valuable insights into the molecular basis of adaptive traits, the impacts of forest management, and changing environmental conditions on tree species and fungal communities and can enhance tree-breeding cycles due to reduced time for field testing. It becomes clear that there are multifaceted interactions among microbial species as well as between these organisms and trees. We demonstrate the versatility of the different approaches based on case studies on trees and fungi. KEY POINTS: • Current knowledge of genetic methods applied to forest trees and associated fungi. • Genomic methods are essential in conservation, breeding, management, and research. • Important role of phytobiomes for trees and their ecosystems.

Rapid plant trait evolution can alter coastal wetland resilience to sea level rise

Rapid evolution remains a largely unrecognized factor in models that forecast the fate of ecosystems under scenarios of global change. In this work, we quantified the roles of heritable variation in plant traits and of trait evolution in explaining variability in forecasts of the state of coastal wetland ecosystems. A common garden study of genotypes of the dominant sedge Schoenoplectus americanus, "resurrected" from time-stratified seed banks, revealed that heritable variation and evolution explained key ecosystem attributes such as the allocation and distribution of belowground biomass. Incorporating heritable trait variation and evolution into an ecosystem model altered predictions of carbon accumulation and soil surface accretion (a determinant of marsh resilience to sea level rise), demonstrating the importance of accounting for evolutionary processes when forecasting ecosystem dynamics.

Compatibility between snails and schistosomes: insights from new genetic resources, comparative genomics, and genetic mapping

The freshwater snail Biomphalaria glabrata is an important intermediate host of the parasite Schistosoma mansoni that causes human intestinal schistosomiasis. To better understand vector snail biology and help advance innovative snail control strategies, we have developed a new snail model consisting of two homozygous B. glabrata lines (iM line and iBS90) with sharply contrasting schistosome-resistance phenotypes. We produced and compared high-quality genome sequences for iM line and iBS90 which were assembled from 255 (N50 = 22.7 Mb) and 346 (N50 = 19.4 Mb) scaffolds, respectively. Using F2 offspring bred from the two lines and the newly generated iM line genome, we constructed 18 linkage groups (representing the 18 haploid chromosomes) covering 96% of the genome and identified three new QTLs (quantitative trait loci), two involved in snail resistance/susceptibility and one relating to body pigmentation. This study provides excellent genomic resources for unveiling complex vector snail biology, reveals genomic difference between resistant and susceptible lines, and offers novel insights into genetic mechanism of the compatibility between snail and schistosome.

Ectomycorrhizal fungal decay traits along a soil nitrogen gradient

The extent to which ectomycorrhizal (ECM) fungi decay soil organic matter (SOM) has implications for accurately predicting forest ecosystem response to climate change. Investigating the distribution of gene traits associated with SOM decay among ectomycorrhizal fungal communities could improve understanding of SOM dynamics and plant nutrition. We hypothesized that soil inorganic nitrogen (N) availability structures the distribution of ECM fungal genes associated with SOM decay and, specifically, that ECM fungal communities occurring in inorganic N-poor soils have greater SOM decay potential. To test this hypothesis, we paired amplicon and shotgun metagenomic sequencing of 60 ECM fungal communities associating with Quercus rubra along a natural soil inorganic N gradient. Ectomycorrhizal fungal communities occurring in low inorganic N soils were enriched in gene families involved in the decay of lignin, cellulose, and chitin. Ectomycorrhizal fungal community composition was the strongest driver of shifts in metagenomic estimates of fungal decay potential. Our study simultaneously illuminates the identity of key ECM fungal taxa and gene families potentially involved in the decay of SOM, and we link rhizomorphic and medium-distance hyphal morphologies with enhanced SOM decay potential. Coupled shifts in ECM fungal community composition and community-level decay gene frequencies are consistent with outcomes of trait-mediated community assembly processes.

Ectomycorrhizal access to organic nitrogen mediates CO2 fertilization response in a dominant temperate tree

Plant–mycorrhizal interactions mediate plant nitrogen (N) limitation and can inform model projections of the duration and strength of the effect of increasing CO₂ on plant growth. We present dendrochronological evidence of a positive, but context-dependent fertilization response of Quercus rubra L. to increasing ambient CO₂ (iCO₂) along a natural soil nutrient gradient in a mature temperate forest. We investigated this heterogeneous response by linking metagenomic measurements of ectomycorrhizal (ECM) fungal N-foraging traits and dendrochronological models of plant uptake of inorganic N and N bound in soil organic matter (N-SOM). N-SOM putatively enhanced tree growth under conditions of low inorganic N availability, soil conditions where ECM fungal communities possessed greater genomic potential to decay SOM and obtain N-SOM. These trees were fertilized by 38 years of iCO₂. In contrast, trees occupying inorganic N rich soils hosted ECM fungal communities with reduced SOM decay capacity and exhibited neutral growth responses to iCO₂. This study elucidates how the distribution of N-foraging traits among ECM fungal communities govern tree access to N-SOM and subsequent growth responses to iCO₂.

Root-mycorrhizal interactions could help explain the heterogeneity of plant responses to CO₂ fertilisation and nutrient availability. Here the authors combine tree-ring and metagenomic data to reveal that tree growth responses to increasing CO₂ along a soil nutrient gradient depend on the nitrogen foraging traits of ectomycorrhizal fungi.

A High-Density SNP Genetic Map Construction Using ddRAD-Seq and Mapping of Capsule Shattering Trait in Sesame

The seed-bearing capsule of sesame shatters at harvest. This wildish trait makes the crop unsuitable for mechanized harvesting and also restricts its commercial potential by limiting the cultivation for countries that have no access to low-cost labor. Therefore, the underlying genetic basis of the capsule shattering trait is highly important in order to develop mechanization-ready varieties for sustainable sesame farming. In the present study, we generated a sesame F₂ population derived from a cross between a capsule shattering cultivar (Muganli-57) and a non-shattering mutant (PI 599446), which was used to construct a genetic map based on double-digest restriction-site-associated DNA sequencing. The resulting high-density genetic map contained 782 single-nucleotide polymorphisms (SNPs) and spanned a length of 697.3 cM, with an average marker interval of 0.89 cM. Based on the reference genome, the capsule shattering trait was mapped onto SNP marker S8_5062843 (78.9 cM) near the distal end of LG8 (chromosome 8). In order to reveal genes potentially controlling the shattering trait, the marker region (S8_5062843) was examined, and a candidate gene including six CDSs was identified. Annotation showed that the gene encodes a protein with 440 amino acids, sharing ∼99% homology with transcription repressor KAN1. Compared with the capsule shattering allele, the SNP change and altered splicing in the flanking region of S8_5062843 caused a frameshift mutation in the mRNA, resulting in the loss of function of this gene in the mutant parent and thus in non-shattering capsules and leaf curling. With the use of genomic data, InDel and CAPS markers were developed to differentiate shattering and non-shattering capsule genotypes in marker-assisted selection studies. The obtained results in the study can be beneficial in breeding programs to improve the shattering trait and enhance sesame productivity.

Transcriptome profiles of Quercus rubra responding to increased O3 stress

BACKGROUND

Climate plays an essential role in forest health, and climate change may increase forest productivity losses due to abiotic and biotic stress. Increased temperature leads to the increased formation of ozone (O3). Ozone is formed by the interaction of sunlight, molecular oxygen and by the reactions of chemicals commonly found in industrial and automobile emissions such as nitrogen oxides and volatile organic compounds. Although it is well known that productivity of Northern red oak (Quercus rubra) (NRO), an ecologically and economically important species in the forests of eastern North America, is reduced by exposure to O3, limited information is available on its responses to exogenous stimuli at the level of gene expression.

RESULTS

RNA sequencing yielded more than 323 million high-quality raw sequence reads. De novo assembly generated 52,662 unigenes, of which more than 42,000 sequences could be annotated through homology-based searches. A total of 4140 differential expressed genes (DEGs) were detected in response to O3 stress, as compared to their respective controls. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analyses of the O3-response DEGs revealed perturbation of several biological pathways including energy, lipid, amino acid, carbohydrate and terpenoid metabolism as well as plant-pathogen interaction.

CONCLUSION

This study provides the first reference transcriptome for NRO and initial insights into the genomic responses of NRO to O3. Gene expression profiling reveals altered primary and secondary metabolism of NRO seedlings, including known defense responses such as terpenoid biosynthesis.

ddRAD Sequencing-Based Identification of Genomic Boundaries and Permeability in Quercus ilex and Q. suber Hybrids

Hybridization and its relevance is a hot topic in ecology and evolutionary biology. Interspecific gene flow may play a key role in species adaptation to environmental change, as well as in the survival of endangered populations. Despite the fact that hybridization is quite common in plants, many hybridizing species, such as Quercus spp., maintain their integrity, while precise determination of genomic boundaries between species remains elusive. Novel high throughput sequencing techniques have opened up new perspectives in the comparative analysis of genomes and in the study of historical and current interspecific gene flow. In this work, we applied ddRADseq technique and developed an ad hoc bioinformatics pipeline for the study of ongoing hybridization between two relevant Mediterranean oaks, Q. ilex and Q. suber. We adopted a local scale approach, analyzing adult hybrids (sensu lato) identified in a mixed stand and their open-pollinated progenies. We have identified up to 9,435 markers across the genome and have estimated individual introgression levels in adults and seedlings. Estimated contribution of Q. suber to the genome is higher, on average, in hybrid progenies than in hybrid adults, suggesting preferential backcrossing with this parental species, maybe followed by selection during juvenile stages against individuals with higher Q. suber genomic contribution. Most discriminating markers seem to be scattered throughout the genome, suggesting that a large number of small genomic regions underlie boundaries between these species. A noticeable proportion of the markers (26%) showed allelic frequencies in adult hybrids very similar to one of the parental species, and very different from the other; a finding that seems relevant for understanding the hybridization process and the occurrence of adaptive introgression. Candidate marker databases developed in this study constitute a valuable resource to design large scale re-sequencing experiments in Mediterranean sclerophyllous oak species and could provide insight in species boundaries and on adaptive introgression between Q. suber and Q. ilex.

Genome-Wide Variant Identification and High-Density Genetic Map Construction Using RADseq for Platycladus orientalis (Cupressaceae)

Platycladus orientalis is an ecologically important native conifer in Northern China and exotic species in many parts of the world; however, knowledge about the species’ genetics and genome are very limited. The availability of well-developed battery of genetic markers, with large genome coverage, is a prerequisite for the species genetic dissection of adaptive attributes and efficient selective breeding. Here, we present a genome-wide genotyping method with double-digestion restriction site associated DNA sequencing (ddRAD-seq) that is effective in generating large number of Mendelian markers for genome mapping and other genetic applications. Using 139 megagametophytes collected from a single mother tree, we assembled 397,226 loci, of which 108,683 (27.4%) were polymorphic. After stringent filtering for 1:1 segregation ratio and missing rate of <20%, the remaining 23,926 loci (22% of the polymorphic loci) were ordered into 11 linkage groups (LGs) and distributed across 7,559 unique positions, with a total map length of 1,443 cM and an average spacing of 0.2 cM between adjacent unique positions. The 11 LGs correspond to the species’ 11 haploid genome chromosome number. This genetic map is among few high-density maps available for conifers to date, and represents the first genetic map for P. orientalis. The information generated serves as a solid foundation not only for marker-assisted breeding efforts, but also for comparative conifer genomic studies.

Dissecting Resistance to Phytophthora cinnamomi in Interspecific Hybrid Chestnut Crosses Using Sequence-Based Genotyping and QTL Mapping

The soilborne oomycete Phytophthora cinnamomi-which causes root rot, trunk cankers, and stem lesions on an estimated 5,000 plant species worldwide-is a lethal pathogen of American chestnut (Castanea dentata) as well as many other woody plant species. P. cinnamomi is particularly damaging to chestnut and chinquapin trees (Castanea spp.) in the southern portion of its native range in the United States due to relatively mild climatic conditions that are conductive to disease development. Introduction of resistant genotypes is the most practical solution for disease management in forests because treatment with fungicides and eradication of the pathogen are neither practical nor economically feasible in natural ecosystems. Using backcross families derived from crosses of American chestnuts with two resistant Chinese chestnut cultivars Mahogany and Nanking, we constructed linkage maps and identified quantitative trait loci (QTLs) for resistance to P. cinnamomi that had been introgressed from these Chinese chestnut cultivars. In total, 957 plants representing five cohorts of three hybrid crosses were genotyped by sequencing and phenotyped by standardized inoculation and visual examination over a 6-year period from 2011 to 2016. Eight parental linkage maps comprising 7,715 markers were constructed, and 17 QTLs were identified on four linkage groups (LGs): LG_A, LG_C, LG_E, and LG_K. The most consistent QTLs were detected on LG_E in seedlings from crosses with both 'Mahogany' and 'Nanking' and LG_K in seedlings from 'Mahogany' crosses. Two consistent large and medium effect QTLs located ∼10 cM apart were present in the middle and at the lower end of LG_E; other QTLs were considered to have small effects. These results imply that the genetic architecture of resistance to P. cinnamomi in Chinese chestnut × American chestnut hybrid progeny may resemble the P. sojae-soybean pathosystem, with a few dominant QTLs along with quantitatively inherited partial resistance conferred by multiple small-effect QTLs.

Back to America: tracking the origin of European introduced populations of Quercus rubra L

Quercus rubra has been introduced in Europe since the end of the 17th century. It is widely distributed today across this continent and considered invasive in some countries. Here, we investigated the distribution of genetic diversity of both native and introduced populations with the aim of tracing the origin of introduced populations. A large sampling of 883 individuals from 73 native and 38 European locations were genotyped at 69 SNPs. In the natural range, we found a continuous geographic gradient of variation with a predominant latitudinal component. We explored the existence of ancestral populations by performing Bayesian clustering analysis and found support for two or three ancestral genetic clusters. Approximate Bayesian Computations analyses based on these two or three clusters support recent extensive secondary contacts between them, suggesting that present-day continuous genetic variation resulted from recent admixture. In the introduced range, one main genetic cluster was not recovered in Europe, suggesting that source populations were preferentially located in the northern part of the natural distribution. However, our results cannot refute the introduction of populations from the southern states that did not survive in Europe.