Araucaria angustifolia has high economic, social, and ecological importance in Brazil, although it is critically threatened with extinction. To understand araucaria’s adaptation, we used a genome-wide association studies (GWAS) to identify markers with signatures of selection associating genomic variation to phenotypic and climatic variables. We also used landscape genomics to identify geographic regions at the highest risk of extinction for the species due to climate change. We used phenotypic and genotypic data of 859 adult trees from a provenance-progeny trial (15 populations), 1,304 SNPs, climatic variables, and growth traits. The GWAS analyses were performed using a general linear model, the Wald test, and a Bayesian method based on population divergence. BLAST techniques were used to gather information about the selected markers. We estimated the proportion of variance explained by regression of genomic data against phenotypic and climatic variables. To estimate vulnerability to climate change, we used the gradient forests. We identified outlier SNPs associated with the climatic and phenotypic traits. Considering the climatic features as drivers of araucaria adaptation, we see that precipitation in the dry season is the leading and most predictable adaptation trait for araucaria. Genomic offset (Goff) for the most optimistic scenario shows that the main critical area is the transition between the tropical and temperate climates in Brazil. In Goff’s most pessimistic scenario, the entire temperate region presents a change in allele turnover. In this context, we propose strategies like assisted migration and targeted reforestation management to accelerate the adaptation of araucaria to the predicted scenarios.

Chloroplast DNA (cpDNA) is a valuable tool for studying plant population genetics and gene flow by pollen in conifers, particularly through the application of chloroplast DNA microsatellite markers (cpSSRs). This study focuses on Chilgoza pine (Pinus gerardiana L.), an economically and ecologically significant tree species in Afghanistan. Despite its importance, comprehensive genetic research on Chilgoza pine has been limited. To address this gap, we developed novel cpSSR markers based on the Chilgoza pine’s chloroplast genome to assess genetic diversity, population structure, and pollen dispersal in a population from Gardiz, Afghanistan. Needle samples from 199 trees across four subpopulations and two age cohorts (young and old) were collected and analyzed. Our findings revealed 27 chloroplast microsatellite markers, of which six exhibited polymorphism. Haplotype analysis identified 32 unique haplotypes, with one most prevalent haplotype. Genetic diversity analysis showed comparatively high levels of diversity, with no genetic differentiation between young and old tree cohorts. Fine-scale spatial genetic structure (FSGS) analysis revealed significant but weak family structure and relatedness in young cohorts, suggesting distance-dependent gene flow. Our study underscores the utility of cpSSRs in characterizing genetic diversity and structure, which is crucial for the conservation and sustainable management of Chilgoza pine forests. These findings provide insights for developing conservation strategies and highlight the importance of genetic marker studies to aid the preservation of biodiversity that supports local livelihoods.

Hymenoscyphus fraxineus, the causal agent of Ash Dieback (ADB), has been introduced to eastern Europe in the 1990s from where it spread causing decline in European ash populations. However, the genetic basis of the molecular response in tolerant and susceptible ash trees to this disease is still largely unknown. We performed RNA-sequencing to study the transcriptomic response to the disease in four ash genotypes (ADB-tolerant FAR3 and FS36, and ADB-susceptible UW1 and UW2), during a time-course of 7, 14, 21, and 28 days post-inoculation, including mock-infected trees as control samples for each sampling time point. The analysis yielded 395 and 500 Differentially Expressed Genes (DEGs) along the response for ADB-tolerant FAR3 and FS36, respectively, while ADB-susceptible UW1 and UW2 revealed 194 and 571 DEGs, respectively, with most DEGs found exclusively in just one of the genotypes. DEGs shared between tolerant genotypes FAR3 and FS36, included genes involved in the production of phytoalexins and other secondary metabolites with roles in plant defense. Moreover, we identified an earlier expression of genes involved in both pattern- and effector-triggered immunity (PTI and ETI) in ADB-tolerant genotypes, while in ADB-susceptible genotypes both responses were delayed (late response). Overall, these results revealed different transcriptomic expression patterns not only between ADB-tolerant and ADB-susceptible genotypes, but also within these two groups. This hints to individual response in the natural tolerance to ADB, possibly revealing diversified strategies across ash genotypes.

High-throughput sequencing (HTS) provides an efficient and cost-effective way to generate large amounts of sequence data. However, marker-based methods and the resulting datasets come with a range of challenges and disputes, including incomplete reference databases, controversial sequence similarity thresholds for delineating taxa, and downstream compositional data analysis. Here, we use HTS data from a soil nematode biodiversity experiment to address the following questions: (1) how the choice of reference database affects HTS data analysis, (2) whether the same ecological patterns are detected with ASV (100% similarity) versus classical OTU (97% similarity), and (3) how different data normalization methods affect the recovery of beta diversity patterns and identification of differentially abundant taxa. At this time, the SILVA database performed better than PR2, assigning more reads to family level and providing higher phylogenetic resolution. ASV- and OTU-based alpha and beta diversity of nematodes correlated closely, indicating that OTU-based studies represent useful reference points. For downstream data analyses, our results indicate that rarefaction-based methods are more vulnerable to missed findings, while clr-transformation based methods may overestimate tested effects. ANCOM-BC retains all data and accounts for uneven sampling fractions for each sample, suggesting that this is currently the optimal method to analyze compositional data. Overall, our study highlights the importance of comparing and selecting taxonomic reference databases before data analyses, and provides solid evidence for the similarity and comparability between OTU- and ASV-based nematode studies. Further, the results highlight the potential weakness of rarefaction-based and clr-transformation based methods. We recommend future studies use ASV and that both the taxonomic reference databases and normalization strategies are carefully tested and selected before analyzing the data.

Northern red oak (Quercus rubra L.) is an ecologically and economically important forest tree native to the northeastern United States. We present a chromosome-scale, haplotype-resolved genome of Q. rubra, a representative red oak species, 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 Megabases (Mb) with 95.27% of the genome sequences scaffolded into 12 chromosomes and 33,333 protein-coding genes. Comparisons to the genomes of Q. lobata and Q. mongolica reveal high collinearity, with intrachromosomal structural variants present. Orthologous gene family analysis with other oak and rosid 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 nine species analyzed. Terpene synthase gene family comparisons further reveal tandem gene duplications in TPS-b subfamily, similar to Q. robur. 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 a family of eight 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 genomics community which will further supplement the knowledge of Quercus genomics.