Genome sequencing and population genomic analyses provide insights into the adaptive landscape of silver birch

Chromosome-level genome assemblies reveal genome evolution of an invasive plant Phragmites australis

Biological invasions pose a significant threat to ecosystems, disrupting local biodiversity and ecosystem functions. The genomic underpinnings of invasiveness, however, are still largely unknown, making it difficult to predict and manage invasive species effectively. The common reed (Phragmites australis) is a dominant grass species in wetland ecosystems and has become particularly invasive when transferred from Europe to North America. Here, we present a high-quality gap-free, telomere-to-telomere genome assembly of Phragmites australis consisting of 24 pseudochromosomes and a B chromosome. Fully phased subgenomes demonstrated considerable subgenome dominance and revealed the divergence of diploid progenitors approximately 30.9 million years ago. Comparative genomics using chromosome-level scaffolds for three other lineages and a previously published draft genome assembly of an invasive lineage revealed that gene family expansions in the form of tandem duplications may have contributed to the invasiveness of the lineage. This study sheds light on the genome evolution of Arundinoideae grasses and suggests that genetic drivers, such as gene family expansions and tandem duplications, may underly the processes of biological invasion in plants. These findings provide a crucial step toward understanding and managing the genetic basis of invasiveness in plant species.

Complex Polyploids: Origins, Genomic Composition, and Role of Introgressed Alleles

Introgression allows polyploid species to acquire new genomic content from diploid progenitors or from other unrelated diploid or polyploid lineages, contributing to genetic diversity and facilitating adaptive allele discovery. In some cases, high levels of introgression elicit the replacement of large numbers of alleles inherited from the polyploid's ancestral species, profoundly reshaping the polyploid's genomic composition. In such complex polyploids, it is often difficult to determine which taxa were the progenitor species and which taxa provided additional introgressive blocks through subsequent hybridization. Here, we use population-level genomic data to reconstruct the phylogenetic history of Betula pubescens (downy birch), a tetraploid species often assumed to be of allopolyploid origin and which is known to hybridize with at least four other birch species. This was achieved by modeling polyploidization and introgression events under the multispecies coalescent and then using an approximate Bayesian computation rejection algorithm to evaluate and compare competing polyploidization models. We provide evidence that B. pubescens is the outcome of an autoploid genome doubling event in the common ancestor of B. pendula and its extant sister species, B. platyphylla, that took place approximately 178,000-188,000 generations ago. Extensive hybridization with B. pendula, B. nana, and B. humilis followed in the aftermath of autopolyploidization, with the relative contribution of each of these species to the B. pubescens genome varying markedly across the species' range. Functional analysis of B. pubescens loci containing alleles introgressed from B. nana identified multiple genes involved in climate adaptation, while loci containing alleles derived from B. humilis revealed several genes involved in the regulation of meiotic stability and pollen viability in plant species.

Cryptic diploid lineage of Betula ermanii at its southern boundary populations in Japan

Polyploidy is thought to enable species diversification and adaptation to extreme environments. Resolving the ecological differences between a taxon's ploidy levels would therefore provide important insights into local adaptation and speciation. The genus Betula includes many polyploids, but estimates of their phylogenetic relationships and evolutionary history are uncertain because of cryptic lineages and species. As one of the southern boundary populations of Betula ermanii in Japan has been shown to have distinctive genetic characteristics and traits, the differences in ploidy levels between three southern boundary and various other Japanese B. ermanii populations were investigated using flow cytometry. Leaf and seed morphologies were also compared. Apart from individuals in southern boundary populations, all those sampled were tetraploid. Individuals from the southern boundary populations were mostly diploid, apart from a few from lower altitude Shikoku populations, which were tetraploid. Leaf and seed morphologies differed between tetraploids and diploids. Diploid individuals were characterized by leaves with a heart-shaped base and many leaf teeth, and seeds with relatively longer wings. The diploid populations could be considered a cryptic relict lineage of B. ermanii, and there is a possibility that this lineage is a diploid ancestor of B. ermanii and a relict population of the Sohayaki element. Further investigation of the Japanese Betula phylogenetic relationships would enable an informed discussion of taxonomic revisions.

The Whole Mitochondrial Genome Sequence of Dendrobium loddigesii Rolfe, an Endangered Orchid Species in China, Reveals a Complex Multi-Chromosome Structure

Dendrobium loddigesii is a precious traditional Chinese medicine with high medicinal and ornamental value. However, the characterization of its mitochondrial genome is still pending. Here, we assembled the complete mitochondrial genome of D. loddigesii and discovered that its genome possessed a complex multi-chromosome structure. The mitogenome of D. loddigesii consisted of 17 circular subgenomes, ranging in size from 16,323 bp to 56,781 bp. The total length of the mitogenome was 513,356 bp, with a GC content of 43.41%. The mitogenome contained 70 genes, comprising 36 protein-coding genes (PCGs), 31 tRNA genes, and 3 rRNA genes. Furthermore, we detected 403 repeat sequences as well as identified 482 RNA-editing sites and 8154 codons across all PCGs. Following the sequence similarity analysis, 27 fragments exhibiting homology to both the mitogenome and chloroplast genome were discovered, accounting for 9.86% mitogenome of D. loddigesii. Synteny analysis revealed numerous sequence rearrangements in D. loddigesii and the mitogenomes of related species. Phylogenetic analysis strongly supported that D. loddigesii and D. Amplum formed a single clade with 100% bootstrap support. The outcomes will significantly augment the orchid mitochondrial genome database, offering profound insights into Dendrobium's intricate mitochondrial genome architecture.

A haplotype-resolved chromosome-level assembly and annotation of European hazelnut (C. avellana cv. Jefferson) provides insight into mechanisms of eastern filbert blight resistance

European hazelnut (Corylus avellana L.) is an important tree nut crop. Hazelnut production in North America is currently limited in scalability due to Anisogramma anomala, a fungal pathogen that causes Eastern Filbert Blight (EFB) disease in hazelnut. Successful deployment of EFB resistant cultivars has been limited to the state of Oregon, where the breeding program at Oregon State University (OSU) has released cultivars with a dominant allele at a single resistance locus identified by classical breeding, linkage mapping, and molecular markers. C. avellana cultivar "Jefferson" is resistant to the predominant EFB biotype in Oregon and has been selected by the OSU breeding program as a model for hazelnut genetic and genomic research. Here, we present a near complete, haplotype-resolved chromosome-level hazelnut genome assembly for "Jefferson". This new assembly is a significant improvement over a previously published genome draft. Analysis of genomic regions linked to EFB resistance and self-incompatibility confirmed haplotype splitting and identified new gene candidates that are essential for downstream molecular marker development, thereby facilitating breeding efforts.

Chromosome-scale genomes of Quercus sichourensis and Quercus rex provide insights into the evolution and adaptation of Fagaceae

The Fagaceae, a plant family with a wide distribution and diverse adaptability, has garnered significant interest as a subject of study in plant speciation and adaptation. Meanwhile, certain Fagaceae species are regarded as highly valuable wood resources due to the exceptional quality of their wood. In this study, we present two high-quality, chromosome-scale genome sequences for Quercus sichourensis (848.75 Mb) and Quercus rex (883.46 Mb). Comparative genomics analysis reveals that the difference in the number of plant disease resistance genes and the nonsynonymous and synonymous substitution ratio (Ka/Ks) of protein-coding genes among Fagaceae species are related to different environmental adaptations. Interestingly, most genes related to starch synthesis in the investigated Quercoideae species are located on a single chromosome, as compared to the outgroup species, Fagus sylvatica. Furthermore, resequencing and population analysis of Q. sichourensis and Q. rex reveal that Q. sichourensis has lower genetic diversity and higher deleterious mutations compared to Q. rex. The high-quality, chromosome-level genomes and the population genomic analysis of the critically endangered Q. sichourensis and Q. rex will provide an invaluable resource as well as insights for future study in these two species, even the genus Quercus, to facilitate their conservation.

Conservation genomics provides insights into genetic resilience and adaptation of the endangered Chinese hazelnut, Corylus chinensis

Global climate change has increased concerns regarding biodiversity loss. However, many key conservation issues still required further research, including demographic history, deleterious mutation load, adaptive evolution, and putative introgression. Here we generated the first chromosome-level genome of the endangered Chinese hazelnut, Corylus chinensis, and compared the genomic signatures with its sympatric widespread C. kwechowensis-C. yunnanensis complex. We found large genome rearrangements across all Corylus species and identified species-specific expanded gene families that may be involved in adaptation. Population genomics revealed that both C. chinensis and the C. kwechowensis-C. yunnanensis complex had diverged into two genetic lineages, forming a consistent pattern of southwestern-northern differentiation. Population size of the narrow southwestern lineages of both species have decreased continuously since the late Miocene, whereas the widespread northern lineages have remained stable (C. chinensis) or have even recovered from population bottlenecks (C. kwechowensis-C. yunnanensis complex) during the Quaternary. Compared with C. kwechowensis-C. yunnanensis complex, C. chinensis showed significantly lower genomic diversity and higher inbreeding level. However, C. chinensis carried significantly fewer deleterious mutations than C. kwechowensis-C. yunnanensis complex, as more effective purging selection reduced the accumulation of homozygous variants. We also detected signals of positive selection and adaptive introgression in different lineages, which facilitated the accumulation of favorable variants and formation of local adaptation. Hence, both types of selection and exogenous introgression could have mitigated inbreeding and facilitated survival and persistence of C. chinensis. Overall, our study provides critical insights into lineage differentiation, local adaptation, and the potential for future recovery of endangered trees.

Allopolyploid origin and diversification of the Hawaiian endemic mints

Island systems provide important contexts for studying processes underlying lineage migration, species diversification, and organismal extinction. The Hawaiian endemic mints (Lamiaceae family) are the second largest plant radiation on the isolated Hawaiian Islands. We generated a chromosome-scale reference genome for one Hawaiian species, Stenogyne calaminthoides, and resequenced 45 relatives, representing 34 species, to uncover the continental origins of this group and their subsequent diversification. We further resequenced 109 individuals of two Stenogyne species, and their purported hybrids, found high on the Mauna Kea volcano on the island of Hawai'i. The three distinct Hawaiian genera, Haplostachys, Phyllostegia, and Stenogyne, are nested inside a fourth genus, Stachys. We uncovered four independent polyploidy events within Stachys, including one allopolyploidy event underlying the Hawaiian mints and their direct western North American ancestors. While the Hawaiian taxa may have principally diversified by parapatry and drift in small and fragmented populations, localized admixture may have played an important role early in lineage diversification. Our genomic analyses provide a view into how organisms may have radiated on isolated island chains, settings that provided one of the principal natural laboratories for Darwin's thinking about the evolutionary process.

Ancient Duplication and Lineage-Specific Transposition Determine Evolutionary Trajectory of ERF Subfamily across Angiosperms

AP2/ERF transcription factor family plays an important role in plant development and stress responses. Previous studies have shed light on the evolutionary trajectory of the AP2 and DREB subfamilies. However, knowledge about the evolutionary history of the ERF subfamily in angiosperms still remains limited. In this study, we performed a comprehensive analysis of the ERF subfamily from 107 representative angiosperm species by combining phylogenomic and synteny network approaches. We observed that the expansion of the ERF subfamily was driven not only by whole-genome duplication (WGD) but also by tandem duplication (TD) and transposition duplication events. We also found multiple transposition events in Poaceae, Brassicaceae, Poales, Brassicales, and Commelinids. These events may have had notable impacts on copy number variation and subsequent functional divergence of the ERF subfamily. Moreover, we observed a number of ancient tandem duplications occurred in the ERF subfamily across angiosperms, e.g., in Subgroup IX, IXb originated from ancient tandem duplication events within IXa. These findings together provide novel insights into the evolution of this important transcription factor family.

The genome and population genomics of allopolyploid Coffea arabica reveal the diversification history of modern coffee cultivars

Coffea arabica, an allotetraploid hybrid of Coffea eugenioides and Coffea canephora, is the source of approximately 60% of coffee products worldwide, and its cultivated accessions have undergone several population bottlenecks. We present chromosome-level assemblies of a di-haploid C. arabica accession and modern representatives of its diploid progenitors, C. eugenioides and C. canephora. The three species exhibit largely conserved genome structures between diploid parents and descendant subgenomes, with no obvious global subgenome dominance. We find evidence for a founding polyploidy event 350,000-610,000 years ago, followed by several pre-domestication bottlenecks, resulting in narrow genetic variation. A split between wild accessions and cultivar progenitors occurred ~30.5 thousand years ago, followed by a period of migration between the two populations. Analysis of modern varieties, including lines historically introgressed with C. canephora, highlights their breeding histories and loci that may contribute to pathogen resistance, laying the groundwork for future genomics-based breeding of C. arabica.

A cryptic syngameon within Betula shrubs revealed: Implications for conservation in changing subarctic environments

Arctic and subarctic ecosystems are rapidly transforming due to global warming, emphasizing the need to understand the genetic diversity and adaptive strategies of northern plant species for effective conservation. This study focuses on Betula glandulosa, a native North American tundra shrub known as dwarf birch, which demonstrates an apparent capacity to adapt to changing climate conditions. To address the taxonomic challenges associated with shrub birches and logistical complexities of sampling in the northernmost areas where species' ranges overlap, we adopted a multicriteria approach. Incorporating molecular data, ploidy level assessment and leaf morphology, we aimed to distinguish B. glandulosa individuals from other shrub birch species sampled. Our results revealed three distinct species and their hybrids within the 537 collected samples, suggesting the existence of a shrub birch syngameon, a reproductive network of interconnected species. Additionally, we identified two discrete genetic clusters within the core species, B. glandulosa, that likely correspond to two different glacial lineages. A comparison between the nuclear and chloroplast SNP data emphasizes a long history of gene exchange between different birch species and genetic clusters. Furthermore, our results highlight the significance of incorporating interfertile congeneric species in conservation strategies and underscores the need for a holistic approach to conservation in the context of climate change, considering the complex dynamics of species interactions. While further research will be needed to describe this shrub birches syngameon and its constituents, this study is a first step in recognizing its existence and disseminating awareness among ecologists and conservation practitioners. This biological phenomenon, which offers evolutionary flexibility and resilience beyond what its constituent species can achieve individually, may have significant ecological implications.

Dissection of figured wood trait in curly birch (Betula pendula Roth var. carelica (Mercklin) Hämet-Ahti) using high-throughput genotyping

Curly (Karelian) birch is a special variety of Betula pendula Roth distributed in the northwestern part of Europe. Karelian birch is well-known for its valuable figured curly wood also known as "wooden marble". The genetic basis underlying curly wood formation has been debated since last century, however, there was no data about loci responsible for the curly wood trait. In the present study, we analyzed two full-sibs populations derived from experimental crosses of curly birches and segregating for the trait. RADseq genotyping was applied to reveal how many loci are involved in 'curliness' formation and to search for genetic variants associated with this trait. One single interval on chromosome 10 was detected containing possible candidate genes. InDel marker BpCW1 was suggested for the first time for marker-assisted selection of trees with curly wood at their earliest stages of development.

Anthropogenic disturbance driving population decline of a dominant tree in East Asia evergreen broadleaved forests over the last 11,000 years

Current biodiversity loss is generally considered to have been caused by anthropogenic disturbance, but it is unclear when anthropogenic activities began to affect biodiversity loss. One hypothesis suggests it began with the Industrial Revolution, whereas others propose that anthropogenic disturbance has been associated with biodiversity decline since the early Holocene. To test these hypotheses, we examined the unique vegetation of evergreen broadleaved forests (EBLFs) in East Asia, where humans have affected landscapes since the early Holocene. We adopted a genomic approach to infer the demographic history of a dominant tree (Litsea elongata) of EBLFs. We used Holocene temperature and anthropogenic disturbance factors to calculate the correlation between these variables and the historical effective population size of L. elongata with Spearman statistics and integrated the maximum-entropy niche model to determine the impact of climate change and anthropogenic disturbance on fluctuation in its effective population size. We identified 9 well-defined geographic clades for the populations of L. elongata. Based on the estimated historical population sizes of these clades, all the populations contracted, indicating persistent population decline over the last 11,000 years. Demographic history of L. elongata and human population change, change in cropland use, and change in irrigated rice area were significantly negatively correlated, whereas climate change in the Holocene was not correlated with demographic history. Our results support the early human impact hypothesis and provide comprehensive evidence that early anthropogenic disturbance may contribute to the current biodiversity crisis in East Asia.

Genome-wide analysis of FRF gene family and functional identification of HvFRF9 under drought stress in barley

FHY3 and its homologous protein FAR1 are the founding members of FRS family. They exhibited diverse and powerful physiological functions during evolution, and participated in the response to multiple abiotic stresses. FRF genes are considered to be truncated FRS family proteins. They competed with FRS for DNA binding sites to regulate gene expression. However, only few studies are available on FRF genes in plants participating in the regulation of abiotic stress. With wide adaptability and high stress-resistance, barley is an excellent candidate for the identification of stress-resistance-related genes. In this study, 22 HvFRFs were detected in barley using bioinformatic analysis from whole genome. According to evolution and conserved motif analysis, the 22 HvFRFs could be divided into subfamilies I and II. Most promoters of subfamily I members contained abscisic acid and methyl jasmonate response elements; however, a large number promoters of subfamily II contained gibberellin and salicylic acid response elements. HvFRF9, one of the members of subfamily II, exhibited a expression advantage in different tissues, and it was most significantly upregulated under drought stress. In-situ PCR revealed that HvFRF9 is mainly expressed in the root epidermal cells, as well as xylem and phloem of roots and leaves, indicating that HvFRF9 may be related to absorption and transportation of water and nutrients. The results of subcellular localization indicated that HvFRF9 was mainly expressed in the nuclei of tobacco epidermal cells and protoplast of arabidopsis. Further, transgenic arabidopsis plants with HvFRF9 overexpression were generated to verify the role of HvFRF9 in drought resistance. Under drought stress, leaf chlorosis and wilting, MDA and O2 - contents were significantly lower, meanwhile, fresh weight, root length, PRO content, and SOD, CAT and POD activities were significantly higher in HvFRF9-overexpressing arabidopsis plants than in wild-type plants. Therefore, overexpression of HvFRF9 could significantly enhance the drought resistance in arabidopsis. These results suggested that HvFRF9 may play a key role in drought resistance in barley by increasing the absorption and transportation of water and the activity of antioxidant enzymes. This study provided a theoretical basis for drought resistance in barley and provided new genes for drought resistance breeding.

Advancing tree genomics to future proof next generation orchard production

The challenges facing tree orchard production in the coming years will be largely driven by changes in the climate affecting the sustainability of farming practices in specific geographical regions. Identifying key traits that enable tree crops to modify their growth to varying environmental conditions and taking advantage of new crop improvement opportunities and technologies will ensure the tree crop industry remains viable and profitable into the future. In this review article we 1) outline climate and sustainability challenges relevant to horticultural tree crop industries, 2) describe key tree crop traits targeted for improvement in agroecosystem productivity and resilience to environmental change, and 3) discuss existing and emerging genomic technologies that provide opportunities for industries to future proof the next generation of orchards.

Peach LAZY1 and DRO1 protein-protein interactions and co-expression with PRAF/RLD family support conserved gravity-related protein interactions across plants

IGT/LAZY proteins play a central role in determining gravitropic set point angle and orientation of lateral organs across plant species. Recent work in model systems has demonstrated that interactions between IGT/LAZY proteins and BREVIS RADIX (BRX)-domain containing proteins, such as PH, RCC1, AND FYVE/RCC1-LIKE DOMAIN (PRAF/RLD), and BREVIS RADIX LIKE (BRXL) family members, are mechanistically important for setting gravitropic set point angle. Here, we identified peach PRAF/RLD proteins as interactors of the peach IGT/LAZY proteins PpeLAZY1 and DEEPER ROOTING 1 (PpeDRO1) from a yeast-two-hybrid screen. We also show that the BRX domains of these interacting proteins have high sequence similarity with PRAF/RLD and BRX family proteins from rice and Arabidopsis. Further, PpeLAZY1 and the peach PRAF/RLD interactors are all expressed at relatively high levels in leaf, meristem, and shoot tip tissues. Together, this evidence supports the importance and conservation of IGT/LAZY-BRX-domain interactions, which underlie setting gravitropic set point angle across angiosperms.

The chromosome-scale genome of Magnolia sinica (Magnoliaceae) provides insights into the conservation of plant species with extremely small populations (PSESP)

Magnolia sinica (Magnoliaceae) is a highly threatened tree endemic to southeast Yunnan, China. In this study, we generated for the first time a high-quality chromosome-scale genome sequence from M. sinica, by combining Illumina and ONT data with Hi-C scaffolding methods. The final assembled genome size of M. sinica was 1.84 Gb, with a contig N50 of ca. 45 Mb and scaffold N50 of 92 Mb. Identified repeats constituted approximately 57% of the genome, and 43,473 protein-coding genes were predicted. Phylogenetic analysis shows that the magnolias form a sister clade with the eudicots and the order Ceratophyllales, while the monocots are sister to the other core angiosperms. In our study, a total of 21 individuals from the 5 remnant populations of M. sinica, as well as 22 specimens belonging to 8 related Magnoliaceae species, were resequenced. The results showed that M. sinica had higher genetic diversity (θw = 0.01126 and θπ = 0.01158) than other related species in the Magnoliaceae. However, population structure analysis suggested that the genetic differentiation among the 5 M. sinica populations was very low. Analyses of the demographic history of the species using different models consistently revealed that 2 bottleneck events occurred. The contemporary effective population size of M. sinica was estimated to be 10.9. The different patterns of genetic loads (inbreeding and numbers of deleterious mutations) suggested constructive strategies for the conservation of these 5 different populations of M. sinica. Overall, this high-quality genome will be a valuable genomic resource for conservation of M. sinica.

UPBEAT1-ROS-POD-PAL System under Different Xylogenesis Scenarios in Karelian Birch (Betula pendula Roth var. carelica (Mercl.) Hämet-Ahti)

BACKGROUND

We studied UPBEAT1 (UPB1) which regulated superoxide radical / hydrogen peroxide ratio together with peroxidase (POD) activity and PAL genes expression under different ways of apical meristem development during the xylem structural elements' formation in unique woody plants B. pendula var. pendula with straight-grained wood and B. pendula var. carelica with figured wood. The differentiation process predominanced in straight-grained wood (B. pendula var. pendula) or proliferation - in the figured wood. The investigation was conducted in the radial row (cambial zone - differentiating xylem - mature xylem) during the active cambial growth period.

OBJECTIVE

The study aimed to study the xylogenesis processes occurring in the 16-year-old straight-grained silver birch (Betula pendula Roth) and Karelian birch (Betula pendula Roth var. carelica (Mercl.) Hämet-Ahti) with figured wood.

METHODS

Hydrogen peroxide and superoxide radical contents and peroxidase activity were determined spectrophotometrically. Gene expression for PAL family genes and the UPBEAT1 gene was assessed using qRT-PCR.

RESULTS

Principal component analysis has confirmed trees with straight-grained and figured wood to be different according to UPBEAT1-ROS-POD-PAL system functioning.

CONCLUSION

The higher superoxide radical/hydrogen peroxide ratio in figured Karelian birch, along with UPBEAT1 transcription factor and PAL genes upregulation, distinguished it from straight-grained silver birch. This metabolic picture confirmed the shift of Karelian birch xylogenesis towards proliferation processes, accompanied by ROS and phenolic compounds' flow and POD activity.

Divergent mechanisms of reduced growth performance in Betula ermanii saplings from high-altitude and low-latitude range edges

The reduced growth performance of individuals from range edges is a common phenomenon in various taxa, and considered to be an evolutionary factor that limits the species' range. However, most studies did not distinguish between two mechanisms that can lead to this reduction: genetic load and adaptive selection to harsh conditions. To address this lack of understanding, we investigated the climatic and genetic factors underlying the growth performance of Betula ermanii saplings transplanted from 11 populations including high-altitude edge and low-latitude edge population. We estimated the climatic position of the populations within the overall B. ermanii's distribution, and the genetic composition and diversity using restriction-site associated DNA sequencing, and measured survival, growth rates and individual size of the saplings. The high-altitude edge population (APW) was located below the 95% significance interval for the mean annual temperature range, but did not show any distinctive genetic characteristics. In contrast, the low-latitude edge population (SHK) exhibited a high level of linkage disequilibrium, low genetic diversity, a distinct genetic composition from the other populations, and a high relatedness coefficient. Both APW and SHK saplings displayed lower survival rates, heights and diameters, while SHK saplings also exhibited lower growth rates than the other populations' saplings. The low heights and diameters of APW saplings was likely the result of adaptive selection to harsh conditions, while the low survival and growth rates of SHK saplings was likely the result of genetic load. Our findings shed light on the mechanisms underlying the reduced growth performance of range-edge populations.

Tree architecture: A strigolactone-deficient mutant reveals a connection between branching order and auxin gradient along the tree stem

Due to their long lifespan, trees and bushes develop higher order of branches in a perennial manner. In contrast to a tall tree, with a clearly defined main stem and branching order, a bush is shorter and has a less apparent main stem and branching pattern. To address the developmental basis of these two forms, we studied several naturally occurring architectural variants in silver birch (Betula pendula). Using a candidate gene approach, we identified a bushy kanttarelli variant with a loss-of-function mutation in the BpMAX1 gene required for strigolactone (SL) biosynthesis. While kanttarelli is shorter than the wild type (WT), it has the same number of primary branches, whereas the number of secondary branches is increased, contributing to its bush-like phenotype. To confirm that the identified mutation was responsible for the phenotype, we phenocopied kanttarelli in transgenic BpMAX1::RNAi birch lines. SL profiling confirmed that both kanttarelli and the transgenic lines produced very limited amounts of SL. Interestingly, the auxin (IAA) distribution along the main stem differed between WT and BpMAX1::RNAi. In the WT, the auxin concentration formed a gradient, being higher in the uppermost internodes and decreasing toward the basal part of the stem, whereas in the transgenic line, this gradient was not observed. Through modeling, we showed that the different IAA distribution patterns may result from the difference in the number of higher-order branches and plant height. Future studies will determine whether the IAA gradient itself regulates aspects of plant architecture.

Telomere-to-telomere and haplotype-resolved genome assembly of the Chinese cork oak (Quercus variabilis)

The Quercus variabilis, a deciduous broadleaved tree species, holds significant ecological and economical value. While a chromosome-level genome for this species has been made available, it remains riddled with unanchored sequences and gaps. In this study, we present a nearly complete comprehensive telomere-to-telomere (T2T) and haplotype-resolved reference genome for Q. variabilis. This was achieved through the integration of ONT ultra-long reads, PacBio HiFi long reads, and Hi-C data. The resultant two haplotype genomes measure 789 Mb and 768 Mb in length, with a contig N50 of 65 Mb and 56 Mb, and were anchored to 12 allelic chromosomes. Within this T2T haplotype-resolved assembly, we predicted 36,830 and 36,370 protein-coding genes, with 95.9% and 96.0% functional annotation for each haplotype genome. The availability of the T2T and haplotype-resolved reference genome lays a solid foundation, not only for illustrating genome structure and functional genomics studies but also to inform and facilitate genetic breeding and improvement of cultivated Quercus species.

Heritable epigenetic modification of BpPIN1 is associated with leaf shapes in Betula pendula

The new variety Betula pendula 'Dalecarlica', selected from Betula pendula, shows high ornamental value owing to its lobed leaf shape. In this study, to identify the genetic components of leaf shape formation, we performed bulked segregant analysis and molecular marker-based fine mapping to identify the causal gene responsible for lobed leaves in B. pendula 'Dalecarlica'. The most significant variations associated with leaf shape were identified within the gene BpPIN1 encoding a member of the PIN-FORMED family, responsible for the auxin efflux carrier. We further confirmed the hypomethylation at the promoter region promoting the expression level of BpPIN1, which causes stronger and longer veins and lobed leaf shape in B. pendula 'Dalecarlica'. These results indicated that DNA methylation at the BpPIN1 promoter region is associated with leaf shapes in B. pendula. Our findings revealed an epigenetic mechanism of BpPIN1 in the regulation of leaf shape in Betula  Linn. (birch), which could help in the molecular breeding of ornamental traits.

Genomic insights into selection for heterozygous alleles and woody traits in Populus tomentosa

Heterozygous alleles are widespread in outcrossing and clonally propagated woody plants. The variation in heterozygosity that underlies population adaptive evolution and phenotypic variation, however, remains largely unknown. Here, we describe a de novo chromosome-level genome assembly of Populus tomentosa, an economic and ecologically important native tree in northern China. By resequencing 302 natural accessions, we determined that the South subpopulation (Pop_S) encompasses the ancestral strains of P. tomentosa, while the Northwest subpopulation (Pop_NW) and Northeast subpopulation (Pop_NE) experienced different selection pressures during population evolution, resulting in significant population differentiation and a decrease in the extent of heterozygosity. Analysis of heterozygous selective sweep regions (HSSR) suggested that selection for lower heterozygosity contributed to the local adaptation of P. tomentosa by dwindling gene expression and genetic load in the Pop_NW and Pop_NE subpopulations. Genome-wide association studies (GWAS) revealed that 88 single nucleotide polymorphisms (SNPs) within 63 genes are associated with nine wood composition traits. Among them, the selection for the homozygous AA allele in PtoARF8 is associated with reductions in cellulose and hemicellulose contents by attenuating PtoARF8 expression, and the increase in lignin content is attributable to the selection for decreases in exon heterozygosity in PtoLOX3 during adaptive evolution of natural populations. This study provides novel insights into allelic variations in heterozygosity associated with adaptive evolution of P. tomentosa in response to the local environment and identifies a series of key genes for wood component traits, thereby facilitating genomic-based breeding of important traits in perennial woody plants.

A single amino acid substitution in MdLAZY1A dominantly impairs shoot gravitropism in Malus

Plant architecture is 1 of the most important factors that determines crop yield potential and productivity. In apple (Malus domestica), genetic improvement of tree architecture has been challenging due to a long juvenile phase and growth as complex trees composed of a distinct scion and a rootstock. To better understand the genetic control of apple tree architecture, the dominant weeping growth phenotype was investigated. We report the identification of MdLAZY1A (MD13G1122400) as the genetic determinant underpinning the Weeping (W) locus that largely controls weeping growth in Malus. MdLAZY1A is 1 of the 4 paralogs in apple that are most closely related to AtLAZY1 involved in gravitropism in Arabidopsis (Arabidopsis thaliana). The weeping allele (MdLAZY1A-W) contains a single nucleotide mutation c.584T>C that leads to a leucine to proline (L195P) substitution within a predicted transmembrane domain that colocalizes with Region III, 1 of the 5 conserved regions in LAZY1-like proteins. Subcellular localization revealed that MdLAZY1A localizes to the plasma membrane and nucleus in plant cells. Overexpressing the weeping allele in apple cultivar Royal Gala (RG) with standard growth habit impaired its gravitropic response and altered the growth to weeping-like. Suppressing the standard allele (MdLAZY1A-S) by RNA interference (RNAi) in RG similarly changed the branch growth direction to downward. Overall, the L195P mutation in MdLAZY1A is genetically causal for weeping growth, underscoring not only the crucial roles of residue L195 and Region III in MdLAZY1A-mediated gravitropic response but also a potential DNA base editing target for tree architecture improvement in Malus and other crops.

Cell polarity linked to gravity sensing is generated by LZY translocation from statoliths to the plasma membrane

Organisms have evolved under gravitational force, and many sense the direction of gravity by means of statoliths in specialized cells. In flowering plants, starch-accumulating plastids, known as amyloplasts, act as statoliths to facilitate downstream gravitropism. The gravity-sensing mechanism has long been considered a mechanosensing process by which amyloplasts transmit forces to intracellular structures, but the molecular mechanism underlying this has not been elucidated. We show here that LAZY1-LIKE (LZY) family proteins involved in statocyte gravity signaling associate with amyloplasts and the proximal plasma membrane. This results in polar localization according to the direction of gravity. We propose a gravity-sensing mechanism by which LZY translocation to the plasma membrane signals the direction of gravity by transmitting information on the position of amyloplasts.

Genomic Analysis of Plastid-Nuclear Interactions and Differential Evolution Rates in Coevolved Genes across Juglandaceae Species

The interaction between the nuclear and chloroplast genomes in plants is crucial for preserving essential cellular functions in the face of varying rates of mutation, levels of selection, and modes of transmission. Despite this, identifying nuclear genes that coevolve with chloroplast genomes at a genome-wide level has remained a challenge. In this study, we conducted an evolutionary rate covariation analysis to identify candidate nuclear genes coevolving with chloroplast genomes in Juglandaceae. Our analysis was based on 4,894 orthologous nuclear genes and 76 genes across seven chloroplast partitions in nine Juglandaceae species. Our results indicated that 1,369 (27.97%) of the nuclear genes demonstrated signatures of coevolution, with the Ycf1/2 partition yielding the largest number of hits (765) and the ClpP1 partition yielding the fewest (13). These hits were found to be significantly enriched in biological processes related to leaf development, photoperiodism, and response to abiotic stress. Among the seven partitions, AccD, ClpP1, MatK, and RNA polymerase partitions and their respective hits exhibited a narrow range, characterized by dN/dS values below 1. In contrast, the Ribosomal, Photosynthesis, Ycf1/2 partitions and their corresponding hits, displayed a broader range of dN/dS values, with certain values exceeding 1. Our findings highlight the differences in the number of candidate nuclear genes coevolving with the seven chloroplast partitions in Juglandaceae species and the correlation between the evolution rates of these genes and their corresponding chloroplast partitions.

Betula mcallisteri sp. nov. (sect. Acuminatae, Betulaceae), a new diploid species overlooked in the wild and in cultivation, and its relation to the widespread B. luminifera

Taxa are traditionally identified using morphological proxies for groups of evolutionarily isolated populations. These proxies are common characters deemed by taxonomists as significant. However, there is no general rule on which character or sets of characters are appropriate to circumscribe taxa, leading to discussions and uncertainty. Birch species are notoriously hard to identify due to strong morphological variability and factors such as hybridization and the existence of several ploidy levels. Here, we present evidence for an evolutionarily isolated line of birches from China that are not distinguishable by traditionally assumed taxon recognition proxies, such as fruit or leaf characters. We have discovered that some wild material in China and some cultivated in the Royal Botanic Gardens Edinburgh, formerly recognized as Betula luminifera, differ from other individuals by having a peeling bark and a lack of cambial fragrance. We use restriction site-associated DNA sequencing and flow cytometry to study the evolutionary status of the unidentified Betula samples to assess the extent of hybridization between the unidentified Betula samples and typical B. luminifera in natural populations. Molecular analyses show the unidentified Betula samples as a distinct lineage and reveal very little genetic admixture between the unidentified samples and B. luminifera. This may also be facilitated by the finding that B. luminifera is tetraploid, while the unidentified samples turned out to be diploid. We therefore conclude that the samples represent a yet unrecognized species, which is here described as Betula mcallisteri.

Differentially Expressed Conserved Plant Vegetative Phase-Change-Related microRNAs in Mature and Rejuvenated Silver Birch In Vitro Propagated Tissues

In plants, phase change from the juvenile stage to maturity involves physiological and anatomical changes, which are initiated and controlled by evolutionary highly conserved microRNAs. This process is of particular significance for the in vitro propagation of woody plant species, as individuals or tissues that have undergone the transition to vegetative maturity are recalcitrant to propagation. Conserved miRNAs differentially expressed between juvenile (including rejuvenated) and mature silver birch tissues were identified using high-throughput sequencing of small RNA libraries. Expression of some miR156 isoforms was high in juvenile tissues and has been previously reported to regulate phase transitions in a range of species. Additional miRNAs, such as miR394 and miR396, that were previously reported to be highly expressed in juvenile woody plant tissues were also differentially expressed in this study. However, expression of miR172, previously reported to be highly expressed in mature tissues, was low in all sample types in this study. The obtained results will provide insight for further investigation of the molecular mechanisms regulating vegetative phase change in silver birch and other perennial woody plant species, by analysing a wider range of genotypes, tissue types and maturation stages. This knowledge can potentially assist in identification of rejuvenated material at an earlier stage than currently possible, increasing the efficiency of silver birch in vitro propagation.

Integrative multi-omics analysis of three early diverged rosid species reveals an ancient hierarchical cold-responsive regulatory network

Elucidating regulators, including transcription factors (TFs) and RNA-binding proteins (RBPs), underlying gene transcriptional and post-transcriptional co-regulatory network is key to understand plant cold responses. Previous studies were mainly conducted on single species, and whether the regulators are conserved across different species remains elusive. Here, we selected three species that diverged at the early evolution of rosids (~99-113 million years ago), performed cold-responsive phylotranscriptome experiments, and integrated chromatin immunoprecipitation- and DNA affinity purification-sequencing (ChIP/DAP-seq) analysis to explore cold-responsive regulators and their regulatory networks. First, we detected over 10,000 cold-induced differentially expressed genes (DEGs) and alternative splicing genes (DASGs) in each species. Among the DEGs, a set of TFs and RBPs were conserved in rosid cold response. Compared to TFs, RBPs displayed a delayed cold-responsive pattern, implying a hierarchical regulation of DEGs and DASGs. By integrating DEGs and DASGs, we identified 259 overlapping DE-DASG orthogroups (closely-related homologs) that were cold-regulated at both transcriptional and post-transcriptional levels in all three studied species. Notably, pathway analysis on each of the DEGs, DASGs, and DE-DASGs in the three species showed a common enrichment connected to the circadian rhythm. Evidently, 226 cold-responsive genes were directly targeted by at least two circadian rhythm components (CCA1, LHY, RV4, RVE7, and RVE8). Finally, we revealed an ancient hierarchy of cold-responsive regulatory networks at transcriptional and post-transcriptional levels launched by circadian components in rosids. Altogether, this study sheds light on conserved regulators underlying cold-responsive regulatory networks across rosid species, despite a long evolutionary history after their divergence.

Fragmented habitat compensates for the adverse effects of genetic bottleneck

In the face of the human-caused biodiversity crisis, understanding the theoretical basis of conservation efforts of endangered species and populations has become increasingly important. According to population genetics theory, population subdivision helps organisms retain genetic diversity, crucial for adaptation in a changing environment. Habitat topography is thought to be important for generating and maintaining population subdivision, but empirical cases are needed to test this assumption. We studied Saimaa ringed seals, landlocked in a labyrinthine lake and recovering from a drastic bottleneck, with additional samples from three other ringed seal subspecies. Using whole-genome sequences of 145 seals, we analyzed the distribution of variation and genetic relatedness among the individuals in relation to the habitat shape. Despite a severe history of genetic bottlenecks with prevalent homozygosity in Saimaa ringed seals, we found evidence for the population structure mirroring the subregions of the lake. Our genome-wide analyses showed that the subpopulations had retained unique variation and largely complementary patterns of homozygosity, highlighting the significance of habitat connectivity in conservation biology and the power of genomic tools in understanding its impact. The central role of the population substructure in preserving genetic diversity at the metapopulation level was confirmed by simulations. Integration of genetic analyses in conservation decisions gives hope to Saimaa ringed seals and other endangered species in fragmented habitats.

Genome structure-based Juglandaceae phylogenies contradict alignment-based phylogenies and substitution rates vary with DNA repair genes

In lineages of allopolyploid origin, sets of homoeologous chromosomes may coexist that differ in gene content and syntenic structure. Presence or absence of genes and microsynteny along chromosomal blocks can serve to differentiate subgenomes and to infer phylogenies. We here apply genome-structural data to infer relationships in an ancient allopolyploid lineage, the walnut family (Juglandaceae), by using seven chromosome-level genomes, two of them newly assembled. Microsynteny and gene-content analyses yield identical topologies that place Platycarya with Engelhardia as did a 1980s morphological-cladistic study. DNA-alignment-based topologies here and in numerous earlier studies instead group Platycarya with Carya and Juglans, perhaps misled by past hybridization. All available data support a hybrid origin of Juglandaceae from extinct or unsampled progenitors nested within, or sister to, Myricaceae. Rhoiptelea chiliantha, sister to all other Juglandaceae, contains proportionally more DNA repair genes and appears to evolve at a rate 2.6- to 3.5-times slower than the remaining species.

Genomic signals of local adaptation and hybridization in Asian white birch

Disentangling the numerous processes that affect patterns of genome-wide diversity in widespread tree species has important implications for taxonomy, conservation, and forestry. Here, we investigate the population genomic structure of Asian white birch (Betula platyphylla) in China and seek to explain it in terms of hybridization, demography and adaptation. We generate whole genome sequence data from 83 individuals across the species range in China. Combining this with an existing data set for 79 European and Russian white birches, we show a clear distinction between B. pendula and B. platyphylla, which have sometimes been lumped taxonomically. Genomic diversity of B. platyphylla in north-western China and Central Russia is affected greatly by hybridization with B. pendula. Excluding these hybridized populations, B. platyphylla in China has a linear distribution from north-eastern to south-western China, along the edge of the inland mountainous region. Within this distribution, three genetic clusters are found, which we model as long diverged with subsequent episodes of gene flow. Patterns of covariation between allele frequencies and environmental variables in B. platyphylla suggest the role of natural selection in the distribution of diversity at 7609 SNPs of which 3767 were significantly differentiated among the genetic clusters. The putative adaptive SNPs are distributed throughout the genome and span 1633 genic regions. Of these genic regions, 87 were previously identified as candidates for selective sweeps in Eurasian B. pendula. We use the 7609 environmentally associated SNPs to estimate the risk of nonadaptedness for each sequenced B. platyphylla individual under a scenario of future climate change, highlighting areas where populations may be under future threat from rising temperatures.

Leaf physiology variations are modulated by natural variations that underlie stomatal morphology in Populus

Stomata are essential for photosynthesis and abiotic stress tolerance. Here, we used multiomics approaches to dissect the genetic architecture and adaptive mechanisms that underlie stomatal morphology in Populus tomentosa juvenile natural population (303 accessions). We detected 46 candidate genes and 15 epistatic gene-pairs, associated with 5 stomatal morphologies and 18 leaf development and photosynthesis traits, through genome-wide association studies. Expression quantitative trait locus mapping revealed that stomata-associated gene loci were significantly associated with the expression of leaf-related genes; selective sweep analysis uncovered significant differentiation in the allele frequencies of genes that underlie stomatal variations. An allelic regulatory network operating under drought stress and adequate precipitation conditions, with three key regulators (DUF538, TRA2 and AbFH2) and eight interacting genes, was identified that might regulate leaf physiology via modulation of stomatal shape and density. Validation of candidate gene variations in drought-tolerant and F₁ hybrid populations of P. tomentosa showed that the DUF538, TRA2 and AbFH2 loci cause functional stabilisation of spatiotemporal regulatory, whose favourable alleles can be faithfully transmitted to offspring. This study provides insights concerning leaf physiology and stress tolerance via the regulation of stomatal determination in perennial plants.

Climate-responsive DNA methylation is involved in the biosynthesis of lignin in birch

Lignin is one of the most important secondary metabolites and essential to the formation of cell walls. Changes in lignin biosynthesis have been reported to be associated with environmental variations and can influence plant fitness and their adaptation to abiotic stresses. However, the molecular mechanisms underlying this association remain unclear. In this study, we evaluated the relations between the lignin biosynthesis and environmental factors and explored the role of epigenetic modification (DNA methylation) in contributing to these relations if any in natural birch. Significantly negative correlations were observed between the lignin content and temperature ranges. Analyzing the transcriptomes of birches in two habitats with different temperature ranges showed that the expressions of genes and transcription factors (TFs) involving lignin biosynthesis were significantly reduced at higher temperature ranges. Whole-genome bisulfite sequencing revealed that promoter DNA methylation of two NAC-domain TFs, BpNST1/2 and BpSND1, may be involved in the inhibition of these gene expressions, and thereby reduced the content of lignin. Based on these results we proposed a DNA methylation-mediated lignin biosynthesis model which responds to environmental factors. Overall, this study suggests the possibility of environmental signals to induce epigenetic variations that result in changes in lignin content, which can aid to develop resilient plants to combat ongoing climate changes or to manipulate secondary metabolite biosynthesis for agricultural, medicinal, or industrial values.

Gravity sensing and responses in the coordination of the shoot gravitropic setpoint angle

Gravity is one of the fundamental environmental cues that affect plant development. Indeed, the plant architecture in the shoots and roots is modulated by gravity. Stems grow vertically upward, whereas lateral organs, such as the lateral branches in shoots, tend to grow at a specific angle according to a gravity vector known as the gravitropic setpoint angle (GSA). During this process, gravity is sensed in specialised gravity-sensing cells named statocytes, which convert gravity information into biochemical signals, leading to asymmetric auxin distribution and driving asymmetric cell division/expansion in the organs to achieve gravitropism. As a hypothetical offset mechanism against gravitropism to determine the GSA, the anti-gravitropic offset (AGO) has been proposed. According to this concept, the GSA is a balance of two antagonistic growth components, that is gravitropism and the AGO. Although the nature of the AGO has not been clarified, studies have suggested that gravitropism and the AGO share a common gravity-sensing mechanism in statocytes. This review discusses the molecular mechanisms underlying gravitropism as well as the hypothetical AGO in the control of the GSA.

Convergent evolution of AP2/ERF III and IX subfamilies through recurrent polyploidization and tandem duplication during eudicot adaptation to paleoenvironmental changes

Whole-genome duplication (WGD or polyploidization) has been suggested as a genetic contributor to angiosperm adaptation to environmental changes. However, many eudicot lineages did not undergo recent WGD (R-WGD) around and/or after the Cretaceous-Paleogene (K-Pg) boundary, times of severe environmental changes; how those plants survived has been largely ignored. Here, we collected 22 plants from major branches of the eudicot phylogeny and classified them into two groups according to the occurrence or absence of R-WGD: 12 R-WGD-containing plants (R-WGD-Y) and 10 R-WGD-lacking plants (R-WGD-N). Subsequently, we identified 496 gene-rich families in R-WGD-Y and revealed that members of the AP2/ERF transcription factor family were convergently over-retained after R-WGDs and showed exceptional cold stimulation. The evolutionary trajectories of the AP2/ERF family were then compared between R-WGD-Y and R-WGD-N to reveal convergent expansions of the AP2/ERF III and IX subfamilies through recurrent independent WGDs and tandem duplications (TDs) after the radiation of the plants. The expansions showed coincident enrichments in- times around and/or after the K-Pg boundary, when global cooling was a major environmental stressor. Consequently, convergent expansions and co-retentions of AP2/ERF III C-repeat binding factor (CBF) duplicates and their regulons in different eudicot lineages contributed to the rewiring of cold-specific regulatory networks. Moreover, promoter analysis of cold-responsive AP2/ERF genes revealed an underlying cis-regulatory code (G-box: CACGTG). We propose a seesaw model of WGDs and TDs in the convergent expansion of AP2/ERF III and IX genes that has contributed to eudicot adaptation during paleoenvironmental changes, and we suggest that TD may be a reciprocal/alternative mechanism for genetic innovation in plants that lack WGD.

Population structure and genetic variation of fragmented mountain birch forests in Iceland

Betula pubescens Ehrh. (mountain birch) is the only forest-forming tree in Iceland. Since human settlement (874 AD), the continuous 25,000-30,000 km 2 forest has shrunk to 1.200 km 2 of fragmented patches, making it a good object to study population genetic consequences of habitat fragmentation and disturbance. Further, genetic studies have also shown that hybridization between the tetraploid (2n=56) B. pubescens and the diploid (2n=28) Betula nana L. (dwarf birch) occurs among Iceland's natural populations. This study assessed the genetic variation within and among eleven birch forests remaining across Iceland. Genotype-by-sequencing methodology (GBS) provided a total of 24,585 SNPs, with a minor allele frequency > 5% for genetic analyses. The analysis showed similar diversity within forests, suggesting that fragmentation and hybridization have had a limited effect on the genetic variation within sites. A clear genetic divergence is found among forests from the different regions of Iceland that may reflect historical isolation; the differentiation between forests increased with geographic distances reflecting isolation by distance. Information on the distribution of genetic variation of birch in Iceland is essential for its conservation and to establish genotype-phenotype associations to predict responses to new environmental conditions imposed by climate change and novel biotic/abiotic stressors.

Genome-wide identification and characterization of OVATE family proteins in Betula luminifera reveals involvement of BlOFP3 and BlOFP5 genes in leaf development

Ovate family proteins (OFP) are plant-specific transcription factors involved in regulating morphologies of the lateral organs, plant growth and development. However, the functional roles of OFP genes in Betula luminifera, an important timber tree species, are not well studied. In this study, we identified 20 BlOFP genes and analyzed their phylogenetic relationship, gene structure, conserved motifs, and cis-elements. Further, expression analysis indicates that BlOFP genes were up-regulated in leaves on the one-year-old branch compared to leaves on the current-year branch and bract, except BlOFP7, BlOFP11, BlOFP14 and BlOFP12. The overexpression of BlOFP3 and BlOFP5 in Arabidopsis thaliana not only resulted in a slower growth rate but also produced sawtooth shape, flatter and darker green rosette leaves. Further investigation showed that the leaf thickness of the transgenic plants was more than double that of the wild type, which was caused by the increasement in the number and size of palisade tissue cells. Furthermore, the expression analysis also indicated that the expressions of several genes related to leaf development were significantly changed in the transgene plants. These results suggested the significant roles of BlOFP3 and BlOFP5 in leaf development. Moreover, protein-protein interaction studies showed that BlOFP3 interacts with BlKNAT5, and BlOFP5 interacts with BlKNAT5, BlBLH6 and BlBLH7. In conclusion, our study demonstrates that BlOFP3 and BlOFP5 were involved in leaf shape and thickness regulation by forming a complex with BlKNAT5, BlBLH6 and BlBLH7. In addition, our study serves as a guide for future functional genomic studies of OFP genes of the B. luminifera.

Genomic insights into rapid speciation within the world's largest tree genus Syzygium

Species radiations, despite immense phenotypic variation, can be difficult to resolve phylogenetically when genetic change poorly matches the rapidity of diversification. Genomic potential furnished by palaeopolyploidy, and relative roles for adaptation, random drift and hybridisation in the apportionment of genetic variation, remain poorly understood factors. Here, we study these aspects in a model radiation, Syzygium, the most species-rich tree genus worldwide. Genomes of 182 distinct species and 58 unidentified taxa are compared against a chromosome-level reference genome of the sea apple, Syzygium grande. We show that while Syzygium shares an ancient genome doubling event with other Myrtales, little evidence exists for recent polyploidy events. Phylogenomics confirms that Syzygium originated in Australia-New Guinea and diversified in multiple migrations, eastward to the Pacific and westward to India and Africa, in bursts of speciation visible as poorly resolved branches on phylogenies. Furthermore, some sublineages demonstrate genomic clines that recapitulate cladogenetic events, suggesting that stepwise geographic speciation, a neutral process, has been important in Syzygium diversification.

Strategic roadmap to assess forest vulnerability under air pollution and climate change

Although it is an integral part of global change, most of the research addressing the effects of climate change on forests have overlooked the role of environmental pollution. Similarly, most studies investigating the effects of air pollutants on forests have generally neglected the impacts of climate change. We review the current knowledge on combined air pollution and climate change effects on global forest ecosystems and identify several key research priorities as a roadmap for the future. Specifically, we recommend (1) the establishment of much denser array of monitoring sites, particularly in the South Hemisphere; (2) further integration of ground and satellite monitoring; (3) generation of flux-based standards and critical levels taking into account the sensitivity of dominant forest tree species; (4) long-term monitoring of N, S, P cycles and base cations deposition together at global scale; (5) intensification of experimental studies, addressing the combined effects of different abiotic factors on forests by assuring a better representation of taxonomic and functional diversity across the ~73,000 tree species on Earth; (6) more experimental focus on phenomics and genomics; (7) improved knowledge on key processes regulating the dynamics of radionuclides in forest systems; and (8) development of models integrating air pollution and climate change data from long-term monitoring programs.

Distinct Taphrina strains from the phyllosphere of birch exhibiting a range of witches' broom disease symptoms

The phyllosphere is an important microbial habitat and reservoir of organisms that modify plant health. Taphrina betulina is the causal agent of birch witches' broom disease. Taphrina species are dimorphic, infecting hosts in the filamentous form and residing in the host phyllosphere as non-infectious yeast. As such, they are expected to be found as resident yeasts on their hosts, even on healthy tissues; however, there is little experimental data supporting this supposition. With the aim of exploring the local infection ecology of T. betulina, we isolated yeasts from the phyllosphere of birch leaves, using three sample classes; infected leaves inside symptom-bearing branches, healthy leaves from symptom-free branches on symptom-bearing trees and leaves from symptom-free branches on symptom-free trees. Isolations yielded 224 yeast strains, representing 11 taxa, including T. betulina, which was the most common isolate and was found in all sample classes, including symptom-free samples. Genotyping revealed genetic diversity among these T. betulina isolates, with seven distinct genotypes differentiated by the markers used. Twenty-two representative T. betulina strains were selected for further study, revealing further phenotypic differences. These findings support that T. betulina is ubiquitous on birch and that individual trees host a diversity of T. betulina strains.

Frozen in time: rangewide genomic diversity, structure, and demographic history of relict American chestnut populations

American chestnut (Castanea dentata) was once the most economically and ecologically important hardwood species in the eastern United States. In the first half of the 20th century, an exotic fungal pathogen - Cryphonectria parasitica - decimated the species, killing billions of chestnut trees. Two approaches to developing blight resistant American chestnut populations show promise, but both will require introduction of adaptive genomic diversity from wild germplasm to produce diverse, locally adapted restoration populations. Here we characterize population structure, demographic history, and genomic diversity in a range-wide sample of 384 wild American chestnuts to inform conservation and breeding with blight resistant varieties. Population structure analyses suggest that the chestnut range can be roughly divided into northeast, central, and southwest populations. Within-population genomic diversity estimates revealed a clinal pattern with the highest diversity in the southwest, which likely reflects bottleneck events associated with Quaternary glaciation. Finally, we identified genomic regions under positive selection within each population, which suggests that defense against fungal pathogens is a common target of selection across all populations. Taken together, these results show that American chestnut underwent a postglacial expansion from the southern portion of its range leading to three extant genetic populations. These populations will serve as management units for breeding adaptive genetic variation into the blight-resistant tree populations for targeted reintroduction efforts.

Complex Physical Structure of Complete Mitochondrial Genome of Quercus acutissima (Fagaceae): A Significant Energy Plant

Quercus acutissima Carruth. is a Chinese important energy plant with high ecological and economic values. While the species chloroplast genome has been reported, its mitochondrial genome (mitogenome) is still unexplored. Here, we assembled and annotated the Q. acutissima mitogenome, and we compared its characteristic differences with several closely related species. The Q. acutissima mitogenome’s main structure is branched with three distinguished contigs (linear molecule 1, circular molecule 2, and circular molecule 3) with 448,982 bp total length and 45.72% GC content. The mitogenome contained 51 genes, including 32 protein-coding, 16 tRNA and 3 rRNA genes. We examined codon usage, repeated sequences, genome recombination, chloroplast to mitochondrion DNA transformation, RNA editing, and synteny in the Q. acutissima mitogenome. Phylogenetic trees based on 29 species mitogenomes clarified the species classification. Our results provided comprehensive information of Q. acutissima mitogenome, and they are expected to provide valuable information for Fagaceae evolutionary biology and to promote the species germplasm utilization.

Changes in the Differentiation Program of Birch Cambial Derivatives following Trunk Girdling

The mechanisms regulating the tree trunk radial growth can be studied in original experiments. One technique for studying cambium activity (the meristem involved in radial growth) under conditions of an increased photoassimilate level is trunk girdling. We girdled the trunks of 17- to 22-year-old silver birch plants (Betula pendula Roth var. pendula) during the active growth period and collected xylem and phloem samples at two height levels (1 cm and 35 cm) above girdle, 10, 20, and 30 days after girdling. We investigated the changes that occurred at the anatomical level, as well as the activities of sucrose-metabolizing enzymes and antioxidant-system enzymes and the expression of genes that encode proteins involved in sucrose and auxin transport and metabolism. A moderate increase in photoassimilates (35 cm above the girdle) resulted in a change in the ratio of phloem to xylem increments and an increase in the proportion of parenchyma in the conducting tissues. The increase of photoassimilates above the level at which they can be used in the processes of normal tissue growth and development (1 cm above the girdle) led to xylogenesis suppression and the stimulation of phloem formation, a significant increase in the parenchyma proportion in the conducting tissues, and formation of large sclereid complexes. The differentiation of parenchyma and sclereid cells coincided with biochemical and molecular markers of abnormal conducting tissue formation in Karelian birch, which are also characterized by high proportions of parenchyma and sclereid near the cambium. The results obtained are important in understanding the cambium responses to the photoassimilate distribution changes and estimating tree productivity and survival under changing environmental conditions.

Chromosome-level genome assembly of the diploid blueberry Vaccinium darrowii provides insights into its subtropical adaptation and cuticle synthesis

Vaccinium darrowii is a subtropical wild blueberry species that has been used to breed economically important southern highbush cultivars. The adaptive traits of V. darrowii to subtropical climates can provide valuable information for breeding blueberry and perhaps other plants, especially against the background of global warming. Here, we assembled the V. darrowii genome into 12 pseudochromosomes using Oxford Nanopore long reads complemented with Hi-C scaffolding technologies, and we predicted 41 815 genes using RNA-sequencing evidence. Syntenic analysis across three Vaccinium species revealed a highly conserved genome structure, with the highest collinearity between V. darrowii and Vaccinium corymbosum. This conserved genome structure may explain the high fertility observed during crossbreeding of V. darrowii with other blueberry cultivars. Analysis of gene expansion and tandem duplication indicated possible roles for defense- and flowering-associated genes in the adaptation of V. darrowii to the subtropics. Putative SOC1 genes in V. darrowii were identified based on phylogeny and expression analysis. Blueberries are covered in a thick cuticle layer and contain anthocyanins, which confer their powdery blue color. Using RNA sequencing, we delineated the cuticle biosynthesis pathways of Vaccinium species in V. darrowii. This result can serve as a reference for breeding berries whose colors are appealing to customers. The V. darrowii reference genome, together with the unique traits of this species, including its diploid genome, short vegetative phase, and high compatibility in hybridization with other blueberries, make V. darrowii a potential research model for blueberry species.

Changes in the Activity of the CLE41/PXY/WOX Signaling Pathway in the Birch Cambial Zone under Different Xylogenesis Patterns

The balance between cell proliferation and differentiation into other cell types is crucial for meristem indeterminacy, and both growth aspects are under genetic control. The peptide-receptor signaling module regulates the activity of the cambial stem cells and the differentiation of their derivatives, along with cytokinins and auxin. We identified the genes encoding the signaling module CLE41-PXY and the regulator of vascular cambium division WOX4 and studied their expression during the period of cambial growth in the radial row: the conducting phloem/cambial zone and the differentiating xylem in two forms of Betula pendula, silver birch and Karelian birch. We have shown that the expression maximum of the BpCLE41/44a gene precedes the expression maximum of the BpPXY gene. Non-figured Karelian birch plants with straight-grained wood are characterized by a more intensive growth and the high expression of CLE41/44-PXY-WOX4. Figured Karelian birch plants, where the disturbed ratio and spatial orientation of structural elements characterizes the wood, have high levels of BpWOX4 expression and a decrease in xylem growth as well as the formation of xylem with a lower vessel density. The mutual influences of CLE41-PXY signaling and auxin signaling on WOX4 gene activity and the proliferation of cambium stem cells are discussed.

BpEIL1 negatively regulates resistance to Rhizoctonia solani and Alternaria alternata in birch

Pathogen attacks affect tree health, causing considerable economic losses as well as serious damage to the surrounding environment. Understanding the disease resistance mechanisms of trees is important for tree breeding. In previous studies on birch (Betula platyphylla × B. pendula), we identified a lesion mimic mutant called lmd. We found that reduced expression of BpEIL1 was responsible for the phenotype in lmd. Following cloning, we acquired several BpEIL1 overexpression and suppression lines in birch. In this study, we cloned the BpEIL1 promoter and found that BpEIL1 was primarily expressed in leaves, particularly in veins. We further studied the traits of transgenic lines and the function of BpEIL1 in disease resistance in birch using the BpEIL1 overexpression line OE9, the suppression line SE13 and the non-transgenic line NT. We found that hydrogen peroxide accumulated in SE13 leaves. Ascorbate peroxidase and catalase activity significantly increased in SE13. SE13 was more resistant to the fungal pathogens Alternaria alternata and Rhizoctonia solani than were the OE9 and NT lines. RNA-seq indicated that pathways related to signal transduction, disease resistance and plant immunity were enriched in SE13. BpEIL1 is thus a negative regulatory transcription factor for disease resistance in birch. This study provides a reference for disease resistance of birch and other trees.

Overexpression of BpERF1.1 in Betula Platyphylla enhanced tolerance to multiple abiotic stresses

Cold, salt and drought are the main environmental stresses affecting plant growth and crop yield. ERF1, a member of the well-studied ERF gene family, has been proven to play an important role in abiotic stress in plants. In this study, an ERF1, BpERF1.1, was overexpressed in birch (Betula Platyphylla Suk.) and 11 transgenic lines were obtained. Effects of cold, salt and drought treatments were  checked  on the BpERF1.1 overexpressing lines and the corresponding WT line. The transgenic lines showed improved tolerance against these abiotic stresses than the WT line. Further RNA-seq analysis identified 689 differentially expressed genes (DEGs) in the transgenic birch compared with WT, including 228 up-regulated genes and 461 down-regulated genes. The result of gene ontology enrichment analysis showed that among these DEGs, 273 genes were involved in various plant biological process, and 83% of them were involved in cellular process, metabolic process, biological regulation and response to stimulus (11%). All these results demonstrate that BpERF1.1 gene improves the tolerance and resistance of birch against cold, salt and drought stress, probably by interconnecting with other genes involved in plant response to abiotic stresses.

Supplementary Information

The online version contains supplementary material available at 10.1007/s12298-022-01206-3.

High-quality genome and methylomes illustrate features underlying evolutionary success of oaks

The genus Quercus, which emerged ∼55 million years ago during globally warm temperatures, diversified into ∼450 extant species. We present a high-quality de novo genome assembly of a California endemic oak, Quercus lobata, revealing features consistent with oak evolutionary success. Effective population size remained large throughout history despite declining since early Miocene. Analysis of 39,373 mapped protein-coding genes outlined copious duplications consistent with genetic and phenotypic diversity, both by retention of genes created during the ancient γ whole genome hexaploid duplication event and by tandem duplication within families, including numerous resistance genes and a very large block of duplicated DUF247 genes, which have been found to be associated with self-incompatibility in grasses. An additional surprising finding is that subcontext-specific patterns of DNA methylation associated with transposable elements reveal broadly-distributed heterochromatin in intergenic regions, similar to grasses. Collectively, these features promote genetic and phenotypic variation that would facilitate adaptability to changing environments.

Genomic evidence for homoploid hybrid speciation between ancestors of two different genera

Homoploid hybrid speciation (HHS) has been increasingly recognized as occurring widely during species diversification of both plants and animals. However, previous studies on HHS have mostly focused on closely-related species while it has been rarely reported or tested between ancestors of different genera. Here, we explore the likely HHS origin of Carpinus sect. Distegocarpus between sect. Carpinus and Ostrya in the family Betulaceae. We generate a chromosome-level reference genome for C. viminea of sect. Carpinus and re-sequence genomes of 44 individuals from the genera Carpinus and Ostrya. Our integrated analyses of all genomic data suggest that sect. Distegocarpus, which has three species, likely originates through HHS during the early divergence between Carpinus and Ostrya. Our study highlights the likelihood of an HHS event between ancestors of the extant genera during their initial divergences, which may have led to reticulate phylogenies at higher taxonomic levels.