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

Phylogenomic analyses highlight innovation and introgression in the continental radiations of Fagaceae across the Northern Hemisphere

Northern Hemisphere forests changed drastically in the early Eocene with the diversification of the oak family (Fagaceae). Cooling climates over the next 20 million years fostered the spread of temperate biomes that became increasingly dominated by oaks and their chestnut relatives. Here we use phylogenomic analyses of nuclear and plastid genomes to investigate the timing and pattern of major macroevolutionary events and ancient genome-wide signatures of hybridization across Fagaceae. Innovation related to seed dispersal is implicated in triggering waves of continental radiations beginning with the rapid diversification of major lineages and resulting in unparalleled transformation of forest dynamics within 15 million years following the K-Pg extinction. We detect introgression at multiple time scales, including ancient events predating the origination of genus-level diversity. As oak lineages moved into newly available temperate habitats in the early Miocene, secondary contact between previously isolated species occurred. This resulted in adaptive introgression, which may have further amplified the diversification of white oaks across Eurasia.

Fagaceae are diverse family including trees of ecological and economic importance. This phylogenomic analysis of nuclear and plastid genomes reconstructs evolutionary history and finds evidence of multiple adaptive introgression events in this important plant family.

Genome-Wide Identification and Analysis of Cell Cycle Genes in Birch

Research Highlights: This study identified the cell cycle genes in birch that likely play important roles during the plant’s growth and development. This analysis provides a basis for understanding the regulatory mechanism of various cell cycles in Betula pendula Roth. Background and Objectives: The cell cycle factors not only influence cell cycles progression together, but also regulate accretion, division, and differentiation of cells, and then regulate growth and development of the plant. In this study, we identified the putative cell cycle genes in the B. pendula genome, based on the annotated cell cycle genes in Arabidopsis thaliana (L.) Heynh. It can be used as a basis for further functional research. Materials and Methods: RNA-seq technology was used to determine the transcription abundance of all cell cycle genes in xylem, roots, leaves, and floral tissues. Results: We identified 59 cell cycle gene models in the genome of B. pendula, with 17 highly expression genes among them. These genes were BpCDKA.1, BpCDKB1.1, BpCDKB2.1, BpCKS1.2, BpCYCB1.1, BpCYCB1.2, BpCYCB2.1, BpCYCD3.1, BpCYCD3.5, BpDEL1, BpDpa2, BpE2Fa, BpE2Fb, BpKRP1, BpKRP2, BpRb1, and BpWEE1. Conclusions: By combining phylogenetic analysis and tissue-specific expression data, we identified 17 core cell cycle genes in the Betulapendula genome.

Improved chromosome-level genome assembly of the Glanville fritillary butterfly (Melitaea cinxia) integrating Pacific Biosciences long reads and a high-density linkage map

BACKGROUND

The Glanville fritillary (Melitaea cinxia) butterfly is a model system for metapopulation dynamics research in fragmented landscapes. Here, we provide a chromosome-level assembly of the butterfly's genome produced from Pacific Biosciences sequencing of a pool of males, combined with a linkage map from population crosses.

RESULTS

The final assembly size of 484 Mb is an increase of 94 Mb on the previously published genome. Estimation of the completeness of the genome with BUSCO indicates that the genome contains 92-94% of the BUSCO genes in complete and single copies. We predicted 14,810 genes using the MAKER pipeline and manually curated 1,232 of these gene models.

CONCLUSIONS

The genome and its annotated gene models are a valuable resource for future comparative genomics, molecular biology, transcriptome, and genetics studies on this species.

Two divergent haplotypes from a highly heterozygous lychee genome suggest independent domestication events for early and late-maturing cultivars

Lychee is an exotic tropical fruit with a distinct flavor. The genome of cultivar ‘Feizixiao’ was assembled into 15 pseudochromosomes, totaling ~470 Mb. High heterozygosity (2.27%) resulted in two complete haplotypic assemblies. A total of 13,517 allelic genes (42.4%) were differentially expressed in diverse tissues. Analyses of 72 resequenced lychee accessions revealed two independent domestication events. The extremely early maturing cultivars preferentially aligned to one haplotype were domesticated from a wild population in Yunnan, whereas the late-maturing cultivars that mapped mostly to the second haplotype were domesticated independently from a wild population in Hainan. Early maturing cultivars were probably developed in Guangdong via hybridization between extremely early maturing cultivar and late-maturing cultivar individuals. Variable deletions of a 3.7 kb region encompassed by a pair of CONSTANS-like genes probably regulate fruit maturation differences among lychee cultivars. These genomic resources provide insights into the natural history of lychee domestication and will accelerate the improvement of lychee and related crops.

Two divergent haplotypes from a highly heterozygous lychee genome of the cultivar ‘Feizixiao’ and resequencing of 72 lychee accessions provide insights into the genome evolution and domestication history of lychee.

Can Forest Trees Cope with Climate Change?-Effects of DNA Methylation on Gene Expression and Adaptation to Environmental Change

Epigenetic modifications, including chromatin modifications and DNA methylation, play key roles in regulating gene expression in both plants and animals. Transmission of epigenetic markers is important for some genes to maintain specific expression patterns and preserve the status quo of the cell. This article provides a review of existing research and the current state of knowledge about DNA methylation in trees in the context of global climate change, along with references to the potential of epigenome editing tools and the possibility of their use for forest tree research. Epigenetic modifications, including DNA methylation, are involved in evolutionary processes, developmental processes, and environmental interactions. Thus, the implications of epigenetics are important for adaptation and phenotypic plasticity because they provide the potential for tree conservation in forest ecosystems exposed to adverse conditions resulting from global warming and regional climate fluctuations.

Silver birch pollen-derived microRNAs promote NF-κB-mediated inflammation in human lung cells

The pollen of Betula pendula Roth (silver birch) is considered to be the main cause of allergy-related rhinitis in Europe and its protein-based allergens such as Bet v 1 are well characterized. However, little is known about non-protein components of birch pollen, e.g., small RNAs and their proinflammatory activity. In the present study, next-generation sequencing (NGS) and bioinformatic approaches were used for silver birch pollen (SBP)-derived microRNA profiling and evaluation of microRNA target genes and pathways in human. Human lung cells, namely WI-38 fibroblasts and A549 alveolar epithelial cells were then stimulated with SBP microRNA in vitro and imaging cytometry-based analysis of the levels of proinflammatory cytokines, autophagy parameters and small RNA processing regulators was conducted. Bioinformatic analysis revealed that SBP microRNA may interfere with autophagy, inflammation and allergy pathways in human. SBP and SBP-derived microRNA induced NF-κB-mediated proinflammatory response in human lung cells as judged by increased levels of NF-κB p65, IL-8 and TNFα. NSUN2 and NSUN5 were involved in pollen-derived microRNA processing. Pollen-derived microRNA also modulated autophagic pathway by changes in the pools of LC3B and p62 that may affect autophagy-based adaptive responses during allergic lung inflammation. We postulate that SBP-derived microRNAs can be considered as novel proinflammatory environmental agents.

Genetic architecture of the metabolic pathway of salicylic acid biosynthesis in Populus

Salicylic acid (SA) is a vital hormone for adaptive responses to biotic and abiotic stresses, which facilitates growth-immunity trade-offs in plants. However, the genetic regulatory networks underlying the metabolic pathway of SA biosynthesis in perennial species remain unclear. Here, we integrated genome-wide association study (GWAS) with metabolite and expression profiling methodologies to dissect the genetic architecture of SA biosynthesis in Populus. First, we quantified nine intermediate metabolites of SA biosynthesis in 300 unrelated Populus tomentosa Carr. individuals. Then, we used a systematic genetic strategy to identify candidate genes for constructing the genetic regulatory network of SA biosynthesis. We focused on WRKY70, an efficient transcription factor, as the key causal gene in the regulatory network, and combined the novel genes coordinating the accumulation of SA. Finally, we identified eight GWAS signals and eight expression quantitative trait loci situated in a selective sweep, and showed the presence of large allele frequency differences among the three geographic populations, revealing that candidate genes subject to selection were involved in SA biosynthesis. This study provides an integrated strategy for dissecting the genetic architecture of the metabolic pathway of SA biosynthesis in Populus, thereby enhancing our understanding of genetic regulation of SA biosynthesis in trees, and accelerating marker-assisted breeding efforts toward high-resistance elite varieties of Populus.

LAZY1-LIKE-mediated gravity signaling pathway in root gravitropic set-point angle control

Gravity signaling components contribute to the control of root gravitropic set-point angle through protein polarization relay within columella.

Transcriptome analysis of North American sweet birch (Betula lenta) revealed a higher expression of genes involved in the biosynthesis of secondary metabolites than European silver birch (B. pendula)

The North American Betula lenta L. (sweet birch) has been used for medicinal reasons for centuries by native Americans. Although sophisticated technologies have rapidly been developed, a large information gap has been observed regarding genetic regulators of medicinally important compounds in sweet birch. Very little is known on the different genes involved in secondary metabolic biosynthesis in sweet birch. To gain a deeper insight into genetic factors, we performed a transcriptome analysis of each three biological samples from different independent trees of sweet and European silver birch (B. pendula Roth). This allowed us to precisely quantify the transcripts of about 24,000 expressed genes including 29 prominent candidate genes putatively involved in the biosynthesis of secondary metabolites like terpenoids, and aromatic benzoic acids. A total number of 597 genes were differentially expressed between B. lenta and B. pendula, while 264 and 210 genes showed upregulation in the bark and leaf of B. lenta, respectively. Moreover, we identified 39 transcriptional regulatory elements, involved in secondary metabolite biosynthesis, upregulated in B. lenta. Our study demonstrated the potential of RNA sequencing to identify candidate genes interacting in secondary metabolite biosynthesis in sweet birch. The candidate genes identified in this study could be subjected to genetic engineering to functionally characterize them in sweet birch. This knowledge can be beneficial to the increase of therapeutically important compounds.

Molecular Genetic Characteristics of Different Scenarios of Xylogenesis on the Example of Two Forms of Silver Birch Differing in the Ratio of Structural Elements in the Xylem

Silver birch (Betula pendula Roth) is an economically important species in Northern Europe. The current research focused on the molecular background of different xylogenesis scenarios in the birch trunks. The study objects were two forms of silver birch, silver birch trees, and Karelian birch trees; the latter form is characterized by the formation of two types of wood, non-figured (straight-grained) and figured, respectively, while it is currently not clear which factors cause this difference. We identified VND/NST/SND genes that regulate secondary cell wall biosynthesis in the birch genome and revealed differences in their expression in association with the formation of xylem with different ratios of structural elements. High expression levels of BpVND7 accompanied differentiation of the type of xylem which is characteristic of the species. At the same time, the appearance of figured wood was accompanied by the low expression levels of the VND genes and increased levels of expression of NST and SND genes. We identified BpARF5 as a crucial regulator of auxin-dependent vascular patterning and its direct target—BpHB8. A decrease in the BpARF5 level expression in differentiating xylem was a specific characteristic of both Karelian birch with figured and non-figured wood. Decreased BpARF5 level expression in non-figured trees accompanied by decreased BpHB8 and VND/NST/SND expression levels compared to figured Karelian birch trees. According to the results obtained, we suggested silver birch forms differing in wood anatomy as valuable objects in studying the regulation of xylogenesis.

Genotypic traits and tradeoffs of fast growth in silver birch, a pioneer tree

Fast-growing and slow-growing plant species are suggested to show integrated economics spectrums and the tradeoffs of fast growth are predicted to emerge as susceptibility to herbivory and resource competition. We tested if these predictions also hold for fast-growing and slow-growing genotypes within a silver birch, Betula pendula population. We exposed cloned saplings of 17 genotypes with slow, medium or fast height growth to reduced insect herbivory, using an insecticide, and to increasing resource competition, using naturally varying field plot grass cover. We measured shoot and root growth, ectomycorrhizal (EM) fungal production using ergosterol analysis and soil N transfer to leaves using ¹⁵N-labelled pulse of NH₄⁺. We found that fast-growing genotypes grew on average 78% faster, produced 56% and 16% more leaf mass and ergosterol, and showed 78% higher leaf N uptake than slow-growing genotypes. The insecticide decreased leaf damage by 83% and increased shoot growth, leaf growth and leaf N uptake by 38%, 52% and 76%, without differences between the responses of fast-growing and slow-growing genotypes, whereas root mass decreased with increasing grass cover. Shoot and leaf growth of fast-growing genotypes decreased and EM fungal production of slow-growing genotypes increased with increasing grass cover. Our results suggest that fast growth is genotypically associated with higher allocation to EM fungi, better soil N capture and greater leaf production, and that the tradeoff of fast growth is sensitivity to competition, but not to insect herbivory. EM fungi may have a dual role: to support growth of fast-growing genotypes under low grass competition and to maintain growth of slow-growing genotypes under intensifying competition.

Genome-Wide Identification, Structure Characterization, and Expression Pattern Profiling of the Aquaporin Gene Family in Betula pendula

Aquaporin water channels (AQPs) constitute a large family of transmembrane proteins present throughout all kingdoms of life. They play key roles in the flux of water and many solutes across the membranes. The AQP diversity, protein features, and biological functions of silver birch are still unknown. A genome analysis of Betula pendula identified 33 putative genes encoding full-length AQP sequences (BpeAQPs). They are grouped into five subfamilies, representing ten plasma membrane intrinsic proteins (PIPs), eight tonoplast intrinsic proteins (TIPs), eight NOD26-like intrinsic proteins (NIPs), four X intrinsic proteins (XIPs), and three small basic intrinsic proteins (SIPs). The BpeAQP gene structure is conserved within each subfamily, with exon numbers ranging from one to five. The predictions of the aromatic/arginine selectivity filter (ar/R), Froger's positions, specificity-determining positions, and 2D and 3D biochemical properties indicate noticeable transport specificities to various non-aqueous substrates between members and/or subfamilies. Nevertheless, overall, the BpePIPs display mostly hydrophilic ar/R selective filter and lining-pore residues, whereas the BpeTIP, BpeNIP, BpeSIP, and BpeXIP subfamilies mostly contain hydrophobic permeation signatures. Transcriptional expression analyses indicate that 23 BpeAQP genes are transcribed, including five organ-related expressions. Surprisingly, no significant transcriptional expression is monitored in leaves in response to cold stress (6 °C), although interesting trends can be distinguished and will be discussed, notably in relation to the plasticity of this pioneer species, B. pendula. The current study presents the first detailed genome-wide analysis of the AQP gene family in a Betulaceae species, and our results lay a foundation for a better understanding of the specific functions of the BpeAQP genes in the responses of the silver birch trees to cold stress.

Strategy by latitude? Higher photosynthetic capacity and root mass fraction in northern than southern silver birch (Betula pendula Roth) in uniform growing conditions

Growth of northern trees is limited by short growing seasons. In multi-year trials, northern trees usually grow less than southern ones but can have higher gas exchange, whereas differences in biomass allocation and its relation to photosynthesis are less known. We characterized silver birch (Betula pendula Roth) provenances from southern (latitude 61°) and northern (latitude 67°) Finland in uniform chamber conditions. In a time-series experiment, we measured traits related to growth, biomass allocation and photosynthesis, and determined gas exchange responses to temperature and light. We found provenance differences in photosynthetic capacity and growth. The northern provenance allocated relatively more to roots, having a higher root mass fraction and lower shoot:root ratio than the southern provenance. On the other hand, the northern provenance had fewer leaves and lower total leaf dry weight (DW) than the southern provenance. The northern provenance attained higher rates of net photosynthesis (Anet) and higher stomatal conductance (gs) in all measured temperatures and higher photosynthesis at the optimum temperature (Aopt) than the southern provenance, but there was no difference in the optimum temperature of photosynthesis (Topt, 18.3 °C for the southern provenance vs 18.9 °C for the northern one). Photosynthetic light response curves showed no between-provenance differences. In a time-series, the northern provenance had higher Anet than the southern provenance, but gs was similar. The northern provenance had higher maximum quantum yield of photosystem II photochemistry (Fv/Fm) than the southern provenance. There were no differences between provenances in height, total plant DW, shoot DW, root DW or shoot mass fraction. Our results suggest that the provenances occupy a common thermal niche, or can at least relatively quickly acclimate to a common growth temperature. Thus, carbon assimilation of these northern trees may not be significantly affected by rising temperatures alone. In an equal photoperiod and optimal conditions, we found different one-season biomass accumulation strategies: southern trees grow with more leaves, while northern trees reach similar total assimilation (total DW, height) with more efficient photosynthetic capacity per leaf area (higher gas exchange, higher Fv/Fm) and relatively more investment in the below-ground fraction of the plant.

Genome-wide identification and analysis of class III peroxidases in Betula pendula

BACKGROUND

Class III peroxidases (POD) proteins are widely present in the plant kingdom that are involved in a broad range of physiological processes including stress responses and lignin polymerization throughout the plant life cycle. At present, POD genes have been studied in Arabidopsis, rice, poplar, maize and Chinese pear, but there are no reports on the identification and function of POD gene family in Betula pendula.

RESULTS

We identified 90 nonredundant POD genes in Betula pendula. (designated BpPODs). According to phylogenetic relationships, these POD genes were classified into 12 groups. The BpPODs are distributed in different numbers on the 14 chromosomes, and some BpPODs were located sequentially in tandem on chromosomes. In addition, we analyzed the conserved domains of BpPOD proteins and found that they contain highly conserved motifs. We also investigated their expression patterns in different tissues, the results showed that some BpPODs might play an important role in xylem, leaf, root and flower. Furthermore, under low temperature conditions, some BpPODs showed different expression patterns at different times.

CONCLUSIONS

The research on the structure and function of the POD genes in Betula pendula plays a very important role in understanding the growth and development process and the molecular mechanism of stress resistance. These results lay the theoretical foundation for the genetic improvement of Betula pendula.

ITS secondary structure reconstruction to resolve taxonomy and phylogeny of the Betula L. genus

The taxonomy and phylogeny of the Betula L. genus remain unresolved and are very difficult to assess due to several factors, especially because of frequent hybridization among different species. In the current study, we used nucleotide sequences of two internal transcribed spacer regions (ITS1 and ITS2), which are commonly used as phylogenetic markers. In addition to their nucleotide variation we reconstructed their secondary structure and used it to resolve phylogenetic relationships of some birch species. We explored whether consideration of secondary structure in phylogenetic analyses based on neighbor-joining, maximum parsimony, maximum likelihood, and Bayesian inference methods would help us obtain more solid support of the reconstructed phylogenetic trees. The results were not unambiguous. There were only a few clades with higher support when secondary structure was included into analysis. The phylogenetic trees generated using different methods were mostly in agreement with each other. However, the resolving power of these markers is still insufficient to reliably discriminate some closely related species. To achieve this aim more reliably there is a need for application of modern genomic approaches in combination with traditional ones.

A chromosome-scale genome assembly of European hazel (Corylus avellana L.) reveals targets for crop improvement

The European hazelnut (Corylus avellana L.) is a tree crop of economic importance worldwide, but especially for northern Turkey, where the majority of production takes place. Hazelnut production is currently challenged by environmental stresses, such as a recent outbreak of severe powdery mildew disease; furthermore, allergy to hazelnuts is an increasing health concern in some regions. In order to provide a foundation for using the available hazelnut genetic resources for crop improvement, we produced a fully assembled genome sequence and annotation for a hazelnut species, from C. avellana cv. 'Tombul', one of the most important Turkish varieties. A hybrid sequencing strategy, combining short reads, long reads and proximity ligation methods, enabled us to resolve heterozygous regions and produce a high-quality 370-Mb assembly that agrees closely with cytogenetic studies and genetic maps of the 11 C. avellana chromosomes, and covers 97.8% of the estimated genome size. The genome includes 27 270 high-confidence protein-coding genes, over 20 000 of which were functionally annotated based on homology with known plant proteins. We focused particularly on gene families encoding hazelnut allergens, and the Mildew resistance Locus O (MLO) proteins that are an important susceptibility factor for powdery mildew. The complete assembly enabled us to differentiate between members of these families and to identify homologues that may be important in mildew disease and hazelnut allergy. These findings provide examples of how the genome can be used to guide research and to develop effective strategies for crop improvement in C. avellana.

Holocene chloroplast genetic variation of shrubs (Alnus alnobetula, Betula nana, Salix sp.) at the siberian tundra-taiga ecotone inferred from modern chloroplast genome assembly and sedimentary ancient DNA analyses

Climate warming alters plant composition and population dynamics of arctic ecosystems. In particular, an increase in relative abundance and cover of deciduous shrub species (shrubification) has been recorded. We inferred genetic variation of common shrub species (Alnus alnobetula, Betula nana, Salix sp.) through time. Chloroplast genomes were assembled from modern plants (n = 15) from the Siberian forest-tundra ecotone. Sedimentary ancient DNA (sedaDNA; n = 4) was retrieved from a lake on the southern Taymyr Peninsula and analyzed by metagenomics shotgun sequencing and a hybridization capture approach. For A. alnobetula, analyses of modern DNA showed low intraspecies genetic variability and a clear geographical structure in haplotype distribution. In contrast, B. nana showed high intraspecies genetic diversity and weak geographical structure. Analyses of sedaDNA revealed a decreasing relative abundance of Alnus since 5,400 cal yr BP, whereas Betula and Salix increased. A comparison between genetic variations identified in modern DNA and sedaDNA showed that Alnus variants were maintained over the last 6,700 years in the Taymyr region. In accordance with modern individuals, the variants retrieved from Betula and Salix sedaDNA showed higher genetic diversity. The success of the hybridization capture in retrieving diverged sequences demonstrates the high potential for future studies of plant biodiversity as well as specific genetic variation on ancient DNA from lake sediments. Overall, our results suggest that shrubification has species-specific trajectories. The low genetic diversity in A. alnobetula suggests a local population recruitment and growth response of the already present communities, whereas the higher genetic variability and lack of geographical structure in B. nana may indicate a recruitment from different populations due to more efficient seed dispersal, increasing the genetic connectivity over long distances.

The Corylus mandshurica genome provides insights into the evolution of Betulaceae genomes and hazelnut breeding

Hazelnut is popular for its flavor, and it has also been suggested that hazelnut is beneficial to cardiovascular health because it is rich in oleic acid. Here, we report the first high-quality chromosome-scale genome for the hazelnut species Corylus mandshurica (2n = 22), which has a high concentration of oleic acid in its nuts. The assembled genome is 367.67 Mb in length, and the contig N50 is 14.85 Mb. All contigs were assembled into 11 chromosomes, and 28,409 protein-coding genes were annotated. We reconstructed the evolutionary trajectories of the genomes of Betulaceae species and revealed that the 11 chromosomes of the hazelnut genus were derived from the most ancestral karyotype in Betula pendula, which has 14 protochromosomes, by inferring homology among five Betulaceae genomes. We identified 96 candidate genes involved in oleic acid biosynthesis, and 10 showed rapid evolution or positive selection. These findings will help us to understand the mechanisms of lipid synthesis and storage in hazelnuts. Several gene families related to salicylic acid metabolism and stress responses experienced rapid expansion in this hazelnut species, which may have increased its stress tolerance. The reference genome presented here constitutes a valuable resource for molecular breeding and genetic improvement of the important agronomic properties of hazelnut.

Resolving phylogeny and polyploid parentage using genus-wide genome-wide sequence data from birch trees

Numerous plant genera have a history including frequent hybridisation and polyploidisation (allopolyploidisation), which means that their phylogeny is a network of reticulate evolution that cannot be accurately depicted as a bifurcating tree with a single tip per species. The genus Betula, which contains many ecologically important tree species, is a case in point. We generated genome-wide sequence reads for 27 diploid and 36 polyploid Betula species or subspecies using restriction site associated DNA (RAD) sequences. These reads were assembled into contigs with a mean length of 675 bp. We reconstructed the evolutionary relationships among diploid Betula species using both supermatrix (concatenation) and species tree methods. We identified the closest diploid relatives of the polyploids according to the relative rates at which reads from polyploids mapped to contigs from different diploid species within a concatenated reference sequence. By mapping reads from allopolyploids to their different putative diploid relatives we assembled contigs from the putative sub-genomes of allopolyploid taxa. We used these to build new phylogenies that included allopolyploid sub-genomes as separate tips. This approach yielded a highly evidenced phylogenetic hypothesis for the genus Betula, including the complex reticulate origins of the majority of its polyploid taxa. Our phylogeny divides the genus into two well supported clades, which, interestingly, differ in their seed-wing morphology. We therefore propose to split Betula into two subgenera.

Genome sequence and evolution of Betula platyphylla

Betula L. (birch) is a pioneer hardwood tree species with ecological, economic, and evolutionary importance in the Northern Hemisphere. We sequenced the Betula platyphylla genome and assembled the sequences into 14 chromosomes. The Betula genome lacks evidence of recent whole-genome duplication and has the same paleoploidy level as Vitis vinifera and Prunus mume. Phylogenetic analysis of lignin pathway genes coupled with tissue-specific expression patterns provided clues for understanding the formation of higher ratios of syringyl to guaiacyl lignin observed in Betula species. Our transcriptome analysis of leaf tissues under a time-series cold stress experiment revealed the presence of the MEKK1-MKK2-MPK4 cascade and six additional mitogen-activated protein kinases that can be linked to a gene regulatory network involving many transcription factors and cold tolerance genes. Our genomic and transcriptome analyses provide insight into the structures, features, and evolution of the B. platyphylla genome. The chromosome-level genome and gene resources of B. platyphylla obtained in this study will facilitate the identification of important and essential genes governing important traits of trees and genetic improvement of B. platyphylla.

Genetic modification of the flavonoid pathway alters growth and reveals flexible responses to enhanced UVB - Role of foliar condensed tannins

Accumulation of certain phenolics is a well-known response of plants to enhanced UVB radiation (280-315 nm), but few experiments have compared the relative importance of different phenolic groups for UVB resilience. To study how an altered phenolic profile affects the responses and resilience of silver birch (Betula pendula) to enhanced UVB, we used RNA interference (RNAi) targeting dihydroflavonol reductase (DFR), anthocyanidin synthase (ANS), or anthocyanidin reductase (ANR) to change the accumulation of phenolics. The unmodified control line and RNAi-modified plants were grown for 51 days under ambient or +32% enhanced UVB dose in a greenhouse. RNAi greatly affected phenolic profile and plant growth. There were no interactive effects of RNAi and UVB on growth or photosynthesis, which indicates that the RNAi and unmodified control plants were equally resilient. UVB enhancement led to an accumulation of foliar flavonoids and condensed tannins, and an increase in the density of stem glands and glandular trichomes on upper leaf surfaces in both the control and RNAi-modified plants. Our results do not indicate a photoprotective role for condensed tannins. However, decreased growth of high-flavonoid low-tannin DFRi and ANRi plants implies that the balance of flavonoids and condensed tannins might be important for normal plant growth.

The roles of the IGT gene family in plant architecture: past, present, and future

Genetic improvement of architectural traits offers tremendous opportunities to dramatically improve crop densities, productivity, and ultimately sustainability. Among these, the orientation, or gravitropic set point angle (GSA), of plant organs is critical to optimize crop profiles, light capture, and nutrient acquisition. Mutant GSA phenotypes have been studied in plants since the 1930's but only recently have the underlying genes been identified. Many of these genes have turned out to fall within the IGT (LAZY1/DRO1/TAC1) family, which initially was not previously recognized due to the lack of sequence conservation of homologous genes across species. Here we discuss recent progress on IGT family genes in various plant species over the past century, review possible functional mechanisms, and provide further analysis of their evolution in land plants and their past and future roles in crop domestication.

Hybrid speciation via inheritance of alternate alleles of parental isolating genes

It is increasingly realized that homoploid hybrid speciation (HHS), which involves no change in chromosome number, is an important mechanism of speciation. HHS will likely increase in frequency as ecological and geographical barriers between species are continuing to be disrupted by human activities. HHS requires the establishment of reproductive isolation between a hybrid and its parents, but the underlying genes and genetic mechanisms remain largely unknown. In this study, we reveal by integrated approaches that reproductive isolation originates in one homoploid hybrid plant species through the inheritance of alternate alleles at genes that determine parental premating isolation. The parent species of this hybrid species are reproductively isolated by differences in flowering time and survivorship on soils containing high concentrations of iron. We found that the hybrid species inherits alleles of parental isolating major genes related to flowering time from one parent and alleles of major genes related to iron tolerance from the other parent. In this way, it became reproductively isolated from one parent by the difference in flowering time and from the other by habitat adaptation (iron tolerance). These findings and further modeling results suggest that HHS may occur relatively easily via the inheritance of alternate parental premating isolating genes and barriers.

Comparative gene expression analysis reveals that multiple mechanisms regulate the weeping trait in Prunus mume

Prunus mume (also known as Mei) is an important ornamental plant that is popular with Asians. The weeping trait in P. mume has attracted the attention of researchers for its high ornamental value. However, the formation of the weeping trait of woody plants is a complex process and the molecular basis of weeping stem development is unclear. Here, the morphological and histochemical characteristics and transcriptome profiles of upright and weeping stems from P. mume were studied. Significant alterations in the histochemical characteristics of upright and weeping stems were observed, and the absence of phloem fibres and less xylem in weeping stems might be responsible for their inability to resist gravity and to grow downward. Transcriptome analysis showed that differentially expressed genes (DEGs) were enriched in phenylpropanoid biosynthesis and phytohormone signal transduction pathways. To investigate the differential responses to hormones, upright and weeping stems were treated with IAA (auxin) and GA₃ (gibberellin A3), respectively, and the results revealed that weeping stems had a weaker IAA response ability and reduced upward bending angles than upright stems. On the contrary, weeping stems had increased upward bending angles than upright stems with GA₃ treatment. Compared to upright stems, interestingly, DEGs associated with diterpenoid biosynthesis and phenylpropanoid biosynthesis were significantly enriched after being treated with IAA, and expression levels of genes associated with phenylpropanoid biosynthesis, ABC transporters, glycosylphosphatidylinositol (GPI)—anchor biosynthesis were altered after being treated with GA₃ in weeping stems. Those results reveal that multiple molecular mechanisms regulate the formation of weeping trait in P. mume, which lays a theoretical foundation for the cultivation of new varieties.

Mitochondrial Phylogenomics of Fagales Provides Insights Into Plant Mitogenome Mosaic Evolution

Fagales are an order of woody plants and comprise more than 1,100 species, most of which produce economically important timbers, nuts, and fruits. Their nuclear and plastid genomes are well-sequenced and provided valuable resources to study their phylogeny, breeding, resistance, etc. However, little is known about the mitochondrial genomes (mitogenomes), which hinder a full understanding of their genome evolution. In this study, we assembled complete mitogenomes of 23 species, covering five of the seven families of Fagales. These mitogenomes had similar gene sets but varied 2.4 times in size. The mitochondrial genes were highly conserved, and their capacity in phylogeny was challenging. The mitogenomic structure was extremely dynamic, and synteny among species was poor. Further analyses of the Fagales mitogenomes revealed extremely mosaic characteristics, with horizontal transfer (HGT)-like sequences from almost all seed plant taxa and even mitoviruses. The largest mitogenome, Carpinus cordata, did not have large amounts of specific sequences but instead contained a high proportion of sequences homologous to other Fagales. Independent and unequal transfers of third-party DNA, including nuclear genome and other resources, may partially account for the HGT-like fragments and unbalanced size expansions observed in Fagales mitogenomes. Supporting this, a mitochondrial plasmid-like of nuclear origin was found in Carpinus. Overall, we deciphered the last genetic materials of Fagales, and our large-scale analyses provide new insights into plant mitogenome evolution and size variation.

Genomic Selection for Forest Tree Improvement: Methods, Achievements and Perspectives

The breeding of forest trees is only a few decades old, and is a much more complicated, longer, and expensive endeavor than the breeding of agricultural crops. One breeding cycle for forest trees can take 20–30 years. Recent advances in genomics and molecular biology have revolutionized traditional plant breeding based on visual phenotype assessment: the development of different types of molecular markers has made genotype selection possible. Marker-assisted breeding can significantly accelerate the breeding process, but this method has not been shown to be effective for selection of complex traits on forest trees. This new method of genomic selection is based on the analysis of all effects of quantitative trait loci (QTLs) using a large number of molecular markers distributed throughout the genome, which makes it possible to assess the genomic estimated breeding value (GEBV) of an individual. This approach is expected to be much more efficient for forest tree improvement than traditional breeding. Here, we review the current state of the art in the application of genomic selection in forest tree breeding and discuss different methods of genotyping and phenotyping. We also compare the accuracies of genomic prediction models and highlight the importance of a prior cost-benefit analysis before implementing genomic selection. Perspectives for the further development of this approach in forest breeding are also discussed: expanding the range of species and the list of valuable traits, the application of high-throughput phenotyping methods, and the possibility of using epigenetic variance to improve of forest trees.

Betula platyphylla BpHOX2 transcription factor binds to different cis-acting elements and confers osmotic tolerance

The homeodomain-leucine zipper (HD-Zip) proteins play crucial roles in plant developmental and environmental responses. However, how they mediate gene expression to facilitate abiotic stress tolerance remains unknown. In the present study, we characterized BpHOX2 (encoding a HD-Zip I family protein) from birch (Betula platyphylla). BpHOX2 is predominately expressed in mature stems and leaves, but expressed at a low level in apical buds and roots, suggesting that it has tissue-specific characteristics. BpHOX2 expression was highly induced by osmotic and salt, but only slightly induced by abscisic acid. Overexpression of BpHOX2 markedly improved osmotic tolerance, while knockdown of BpHOX2 increased sensitivity to osmotic stress. BpHOX2 could induce the expression of pyrroline-5-carboxylate synthase, peroxidase, and superoxide dismutase genes to improve proline levels and the reactive oxygen species scavenging capability. Chromatin immunoprecipitation sequencing combined with RNA sequencing showed that BpHOX2 could bind to at least four cis-acting elements, including dehydration-responsive element "RCCGAC", Myb-p binding box "CCWACC," and two novel cis-acting elements with the sequences of "AAGAAG" and "TACGTG" (termed HBS1 and HBS2, respectively) to regulate gene expression. Our results suggested that BpHOX2 is a transcription factor that binds to different cis-acting elements to regulate gene expression, ultimately improving osmotic tolerance in birch.

Expression Analysis of Key Auxin Biosynthesis, Transport, and Metabolism Genes of Betula pendula with Special Emphasis on Figured Wood Formation in Karelian Birch

Auxin status in woody plants is believed to be a critical factor for the quantity and quality of the wood formed. It has been previously demonstrated that figured wood formation in Karelian birch (Betula pendula Roth var. carelica (Merckl.) Hämet-Ahti) is associated with a reduced auxin level and elevated sugar content in the differentiating xylem, but the molecular mechanisms of the abnormal xylogenesis remained largely unclear. We have identified genes involved in auxin biosynthesis (Yucca), polar auxin transport (PIN) and the conjugation of auxin with amino acids (GH3) and UDP-glucose (UGT84B1) in the B. pendula genome, and analysed their expression in trunk tissues of trees differing in wood structure. Almost all the investigated genes were overexpressed in Karelian birch trunks. Although Yucca genes were overexpressed, trunk tissues in areas developing figured grain had traits of an auxin-deficient phenotype. Overexpression of GH3s and UGT84B1 appears to have a greater effect on figured wood formation. Analysis of promoters of the differentially expressed genes revealed a large number of binding sites with various transcription factors associated with auxin and sugar signalling. These data agree with the hypothesis that anomalous figured wood formation in Karelian birch may be associated with the sugar induction of auxin conjugation.

Identification and analysis of key genes involved in methyl salicylate biosynthesis in different birch species

Species of the perennial woody plant genus Betula dominate subalpine forests and play a significant role in preserving biological diversity. In addition to their conventional benefits, birches synthesize a wide range of secondary metabolites having pharmacological significance. Methyl salicylate (MeSA) is one of these naturally occurring compounds constitutively produced by different birch species. MeSA is therapeutically important in human medicine for muscle injuries and joint pain. However, MeSA is now mainly produced synthetically due to a lack of information relating to MeSA biosynthesis and regulation. In this study, we performed a comprehensive bioinformatics analysis of two candidate genes mediating MeSA biosynthesis, SALICYLIC ACID METHYLTRANSFERASE (SAMT) and SALICYLIC ACID-BINDING PROTEIN 2 (SABP2), of high (B. lenta, B. alleghaniensis, B. medwediewii, and B. grossa) and low (B. pendula, B. utilis, B. alnoides, and B. nana) MeSA-producing birch species. Phylogenetic analyses of SAMT and SABP2 genes and homologous genes from other plant species confirmed their evolutionary relationships. Multiple sequence alignments of the amino acid revealed the occurrence of important residues for substrate specificity in SAMT and SABP2. The analysis of cis elements in different birches indicated a functional multiplicity of SAMT and SABP2 and provided insights into the regulation of both genes. We successfully developed six prominent single nucleotide substitution markers that were validated with 38 additional birch individuals to differentiate high and low MeSA-producing birch species. Relative tissue-specific expression analysis of SAMT in leaf and bark tissue of two high and two low MeSA-synthesizing birches revealed a high expression in the bark of both high MeSA-synthesizing birches. In contrast, SABP2 expression in tissues revealed indifferent levels of expression between species belonging to the two groups. The comparative expression and bioinformatics analyses provided vital information that could be used to apply plant genetic engineering technology in the mass production of organic MeSA.

The genome sequence of Mesua ferrea and comparative demographic histories of forest trees

Mesua ferrea (Family: Calophyllaceae) is a tropical forest plant used for timber, biofuel, and traditional medicine. Colloquially, it is known as Nagkesar (Cobra saffron) and is the state flower of Tripura (India). In this study, we perform the whole-genome assembly of Mesua ferrea using ~180X coverage paired-end Illumina data. Our de novo assembly is 614 Mega-base pair (Mbp), has an N50 of 392 Kilo-base pairs (Kbp), and an assembly quality comparable to other published Malpighiales genomes. Further, we collate the genomic datasets of 14 additional forest tree species to compare the temporal dynamics of Effective Population Size (N_e) and find evidence of a substantial bottleneck in all tropical forest plants during Mid-Pleistocene glaciations. The availability of this high-quality draft genome assembly will prove to be a useful resource for functional and comparative genomic studies.

Genomic analysis reveals the genetic diversity, population structure, evolutionary history and relationships of Chinese pepper

Chinese pepper, mainly including Zanthoxylum bungeanum and Zanthoxylum armatum, is an economically important crop popular in Asian countries due to its unique taste characteristics and potential medical uses. Numerous cultivars of Chinese pepper have been developed in China through long-term domestication. To better understand the population structure, demographic history, and speciation of Chinese pepper, we performed a comprehensive analysis at a genome-wide level by analyzing 38,395 genomic SNPs that were identified in 112 cultivated and wild accessions using a high-throughput genome-wide genotyping-by-sequencing (GBS) approach. Our analysis provides genetic evidence of multiple splitting events occurring between and within species, resulting in at least four clades in Z. bungeanum and two clades in Z. armatum. Despite no evidence of recent admixture between species, we detected substantial gene flow within species. Estimates of demographic dynamics and species distribution modeling suggest that climatic oscillations during the Pleistocene (including the Penultimate Glaciation and the Last Glacial Maximum) and recent domestication events together shaped the demography and evolution of Chinese pepper. Our analyses also suggest that southeastern Gansu province is the most likely origin of Z. bungeanum in China. These findings provide comprehensive insights into genetic diversity, population structure, demography, and adaptation in Zanthoxylum.

Mitochondrial Genome of Fagus sylvatica L. as a Source for Taxonomic Marker Development in the Fagales

European beech, Fagus sylvatica L., is one of the most important and widespread deciduous tree species in Central Europe and is widely managed for its hard wood. The complete DNA sequence of the mitochondrial genome of Fagus sylvatica L. was assembled and annotated based on Illumina MiSeq reads and validated using long reads from nanopore MinION sequencing. The genome assembled into a single DNA sequence of 504,715 bp in length containing 58 genes with predicted function, including 35 protein-coding, 20 tRNA and three rRNA genes. Additionally, 23 putative protein-coding genes were predicted supported by RNA-Seq data. Aiming at the development of taxon-specific mitochondrial genetic markers, the tool SNPtax was developed and applied to select genic SNPs potentially specific for different taxa within the Fagales. Further validation of a small SNP set resulted in the development of four CAPS markers specific for Fagus, Fagaceae, or Fagales, respectively, when considering over 100 individuals from a total of 69 species of deciduous trees and conifers from up to 15 families included in the marker validation. The CAPS marker set is suitable to identify the genus Fagus in DNA samples from tree tissues or wood products, including wood composite products.

Population sequencing enhances understanding of tea plant evolution

Tea is an economically important plant characterized by a large genome, high heterozygosity, and high species diversity. In this study, we assemble a 3.26-Gb high-quality chromosome-scale genome for the 'Longjing 43' cultivar of Camellia sinensis var. sinensis. Genomic resequencing of 139 tea accessions from around the world is used to investigate the evolution and phylogenetic relationships of tea accessions. We find that hybridization has increased the heterozygosity and wide-ranging gene flow among tea populations with the spread of tea cultivation. Population genetic and transcriptomic analyses reveal that during domestication, selection for disease resistance and flavor in C. sinensis var. sinensis populations has been stronger than that in C. sinensis var. assamica populations. This study provides resources for marker-assisted breeding of tea and sets the foundation for further research on tea genetics and evolution.

Transcriptome Sequencing Analysis of Birch (Betula platyphylla Sukaczev) under Low-Temperature Stress

Low temperature is one of the common abiotic stresses that adversely affect the growth and development of plants. In this study, we used RNA-Seq to identify low-temperature-responsive genes in birch and further analyzed the underlying molecular mechanism. Birch seedlings were treated by the low temperature (6 °C) for 0, 1, 1.5, 2, 2.5, and 3 h, respectively. A total of 3491 genes were differentially expressed after low-temperature stress. Gene Ontology (GO) and Kyoto Encyclopedia of Gene and Genomes (KEGG) functional enrichment analysis were performed for the differentially expressed genes (DEGs). GO analysis indicated that 3491 DEGs were distributed into 1002 categories, and these DEGs were enriched in “cell process”, “metabolic process”, and “stimulus response”, under the “biological process” category; in “organelles” and “cell components”, under the “cell component” category; and in “catalytic activity” and “adhesion”, under the “molecular function” category. The KEGG enrichment indicated that 119 DEGs were involved in Ca2+ and plant hormone signal transduction; 205 DEGs were involved in secondary metabolic processes, such as lipid metabolism and phenylpropanoid biosynthesis pathway; and 20 DEGs were involved in photosynthesis. In addition, a total of 362 transcription factors (TFs) were differentially expressed under low-temperature stress, including AP2/ERF, C2H2, MYB-HB-like, WRKY, bHLH, WD40-like, and GRAS families. Gene Bpev01.c0480.g0081 (calmodulin-like CML38), Bpev01.c1074.g0005 (calmodulin-like CML25), Bpev01.c1074.g0001 (Calcium-binding EF-hand family protein), Bpev01.c2029.g0005 (calmodulin-like protein), Bpev01.c0154.g0008 (POD), Bpev01.c0015.g0143 (N-acetyl-1-glutamate synthase), and Bpev01.c0148.g0010 (branched chain amino acid transferase) were up-regulated at a high level, under low-temperature stress.

Genomic insights into historical population dynamics, local adaptation, and climate change vulnerability of the East Asian Tertiary relict Euptelea (Eupteleaceae)

The warm-temperate and subtropical climate zones of East Asia are a hotspot of plant species richness and endemism, including a noticeable number of species-poor Tertiary relict tree genera. However, little is understood about when East Asian Tertiary relict plants diversified, how they responded demographically to past environmental change, and to what extent their current genomic composition (and adaptive capacity) might mitigate the effects of global warming. Here, we obtained genomic (RAD-SNP) data for 171 samples from two extant species of Euptelea in China (24 E. pleiosperma populations) and Japan (11 E. polyandra populations) to elucidate their divergence and demographic histories, genome-wide associations with current environmental variables, and genomic vulnerability to future climate change. Our results indicate that Late Miocene changes in climate and/or sea level promoted species divergence, whereas Late Pliocene uplifting in southwest China likely fostered lineage divergence within E. pleiosperma. Its subsequent range expansion into central/east (CE) China bears genomic signatures of climate-driven selection, yet extant CE populations are predicted to be most vulnerable to future climate change. For E. polyandra, geography was the only significant predictor of genomic variation. Our findings indicate a profound impact of Late Neogene geological and climate change on the evolutionary history of Euptelea, with much stronger signals of local adaptation left in China than in Japan. This study deepens our understanding of the complex evolutionary forces that influence the distribution of genetic variation of Tertiary relict trees, and provides insights into their susceptibility to global change and potential for adaptive responses. Our results lay the groundwork for future conservation and restoration programs for Euptelea.

Tandem gene duplications drive divergent evolution of caffeine and crocin biosynthetic pathways in plants

Background

Plants have evolved a panoply of specialized metabolites that increase their environmental fitness. Two examples are caffeine, a purine psychotropic alkaloid, and crocins, a group of glycosylated apocarotenoid pigments. Both classes of compounds are found in a handful of distantly related plant genera (Coffea, Camellia, Paullinia, and Ilex for caffeine; Crocus, Buddleja, and Gardenia for crocins) wherein they presumably evolved through convergent evolution. The closely related Coffea and Gardenia genera belong to the Rubiaceae family and synthesize, respectively, caffeine and crocins in their fruits.

Results

Here, we report a chromosomal-level genome assembly of Gardenia jasminoides, a crocin-producing species, obtained using Oxford Nanopore sequencing and Hi-C technology. Through genomic and functional assays, we completely deciphered for the first time in any plant the dedicated pathway of crocin biosynthesis. Through comparative analyses with Coffea canephora and other eudicot genomes, we show that Coffea caffeine synthases and the first dedicated gene in the Gardenia crocin pathway, GjCCD4a, evolved through recent tandem gene duplications in the two different genera, respectively. In contrast, genes encoding later steps of the Gardenia crocin pathway, ALDH and UGT, evolved through more ancient gene duplications and were presumably recruited into the crocin biosynthetic pathway only after the evolution of the GjCCD4a gene.

Conclusions

This study shows duplication-based divergent evolution within the coffee family (Rubiaceae) of two characteristic secondary metabolic pathways, caffeine and crocin biosynthesis, from a common ancestor that possessed neither complete pathway. These findings provide significant insights on the role of tandem duplications in the evolution of plant specialized metabolism.

Genetic Analysis by nuSSR Markers of Silver Birch (Betula pendula Roth) Populations in Their Southern European Distribution Range

In the main distribution area the genetic pattern of silver birch is dominated by two haplotypes: haplotype A located in the western and north-western Europe, and haplotype C in eastern and southeastern Europe, characterized by high levels of neutral genetic variability within populations, and low differentiation among populations. Information about the amount and structure of genetic variation in the southern marginal areas, representing rear populations left during the expansion of this species from southern glacial refugia, are lacking. The general aim of the study was to investigate the existence of the climatic characteristics typical of the environmental niche of the species, jointly to genetic organization, variation and gene flow, in marginal populations on the Italian Apennines and Greek Southern Rhodope and compare them with populations of the southern part of the main distribution range on the Alps and Balkans. Genetic analysis was performed using nuclear microsatellites loci on 311 trees sampled from 14 populations. Environmental analysis was performed on the multivariate analysis of derived climatic variables. The allelic pattern was analyzed to assess genetic diversity, population diversity and differentiation, population structure and gene flow. The geographic and environmental peripherality did not always match, with some Apennine sites at higher elevation enveloped in the environmental niche. In the peripheral populations on the Apennines, we observed a lower genetic diversity and higher differentiation, with evident genetic barriers detected around these sites. These characteristics were not shown in the marginal Greek populations. Unexpectedly, the southern Italian marginal populations showed genetic links with the Greek and central area of the distribution range. The Greek populations also showed evident gene flow with the Alpine and Balkan areas. The disparity of results in these two marginal areas show that it is not the geographic peripherality or even the ecological marginality that may shape the genetic diversity and structure of marginal populations, but primarily their position as part of the continuous range or as disjunct populations. This outcome suggests different considerations on how to manage their gene pools and the role that these rear populations can play in maintaining the biodiversity of this species.

Overexpression of an AP2/ERF family gene, BpERF13, in birch enhances cold tolerance through upregulating CBF genes and mitigating reactive oxygen species

The AP2/ERF (APETALA2/ethylene-responsive factor) family of transcription factors (TF) is involved in regulating biotic and abiotic stress responses in plants. To explore the role of AP2/ERFs in cold tolerance in woody plants, BpERF13 was cloned and characterized in Betula platyphylla (white birch), a species primarily found in Asia in temperate and boreal climates. Based on phylogenetic analysis, BpERF13 is a member of the IXb subfamily of ERFs. Using qRT-PCR, we found that BpERF13 was differentially expressed in different tissues, and its expression could be induced by cold treatment (4 °C). BpERF13 protein, fused with GFP, was exclusively localized to nuclei. To further assess the role of BpERF13 in cold tolerance, BpERF13 overexpression (OE) transgenic lines were generated in B. platyphylla and used for cold stress treatment and biochemical/physiological studies. BpERF13 overexpression lines had significantly increased tolerance to subfreezing treatment and reduced reactive oxygen species. Using a TF-centered yeast one-hybrid (Y1H) experimental system, we showed that BpERF13 could bind to LTRECOREATCOR15 and MYBCORE cis-elements to activate a reporter gene. ChIP-seq and ChIP-PCR experiments further demonstrated that BpERF13 bound to these cis-elements when present in the 5' proximal regions of superoxide dismutase (SOD), peroxidase (POD), and C-repeat-binding factor (CBF) genes. qRT-PCR was employed to examine the expression levels of these genes in response to cold stress; SOD, POD, and CBF genes were significantly upregulated in BpERF13 transgenic lines compared to wild-type plants in response to cold stress. These results indicate that the transcription factor BpERF13 regulates physiological processes underlying cold tolerance in woody plants.

Differences in growth and gas exchange between southern and northern provenances of silver birch (Betula pendula Roth) in northern Europe

Due to its ubiquity across northern latitudes, silver birch (Betula pendula Roth) is an attractive model species for studying geographical trait variation and acclimation capacity. Six birch provenances from 60 to 67°N across Finland were grown in a common garden and studied for provenance and genotype variation. We looked for differences in height growth, photosynthetic gas exchange and chlorophyll content index (CCI) and compared the gas exchange of early and late leaves on short and long shoots, respectively. The provenances stratified into southern and northern groups. Northern provenances attained less height growth increment and had higher stomatal conductance (gs) and lower intrinsic water-use efficiency (WUE, Anet/gs) than southern provenances, whereas net photosynthesis (Anet) or CCI did not show clear grouping. Short shoot leaves had lower gs and higher WUE than long shoot leaves in all provenances, but there was no difference in Anet between shoot types. The separation of the provenances into two groups according to their physiological responses might reflect the evolutionary history of B. pendula. Latitudinal differences in gas exchange and water use traits can have plausible consequences for global carbon and water fluxes in a warming climate.

Chromosomal-level assembly of Juglans sigillata genome using Nanopore, BioNano, and Hi-C analysis

BACKGROUND

Juglans sigillata, or iron walnut, belonging to the order Juglandales, is an economically important tree species in Asia, especially in the Yunnan province of China. However, little research has been conducted on J. sigillata at the molecular level, which hinders understanding of its evolution, speciation, and synthesis of secondary metabolites, as well as its wide adaptability to its plateau environment. To address these issues, a high-quality reference genome of J. sigillata would be useful.

FINDINGS

To construct a high-quality reference genome for J. sigillata, we first generated 38.0 Gb short reads and 66.31 Gb long reads using Illumina and Nanopore sequencing platforms, respectively. The sequencing data were assembled into a 536.50-Mb genome assembly with a contig N50 length of 4.31 Mb. Additionally, we applied BioNano technology to identify contacts among contigs, which were then used to assemble contigs into scaffolds, resulting in a genome assembly with scaffold N50 length of 16.43 Mb and contig N50 length of 4.34 Mb. To obtain a chromosome-level genome assembly, we constructed 1 Hi-C library and sequenced 79.97 Gb raw reads using the Illumina HiSeq platform. We anchored ∼93% of the scaffold sequences into 16 chromosomes and evaluated the quality of our assembly using the high contact frequency heat map. Repetitive elements account for 50.06% of the genome, and 30,387 protein-coding genes were predicted from the genome, of which 99.8% have been functionally annotated. The genome-wide phylogenetic tree indicated an estimated divergence time between J. sigillata and Juglans regia of 49 million years ago on the basis of single-copy orthologous genes.

CONCLUSIONS

We provide the first chromosome-level genome for J. sigillata. It will lay a valuable foundation for future research on the genetic improvement of J. sigillata.

ELIMÄKI Locus Is Required for Vertical Proprioceptive Response in Birch Trees

Tree architecture has evolved to support a top-heavy above-ground biomass, but this integral feature poses a weight-induced challenge to trunk stability. Maintaining an upright stem is expected to require vertical proprioception through feedback between sensing stem weight and responding with radial growth. Despite its apparent importance, the principle by which plant stems respond to vertical loading forces remains largely unknown. Here, by manipulating the stem weight of downy birch (Betula pubescens) trees, we show that cambial development is modulated systemically along the stem. We carried out a genetic study on the underlying regulation by combining an accelerated birch flowering program with a recessive mutation at the ELIMÄKI locus (EKI), which causes a mechanically defective response to weight stimulus resulting in stem collapse after just 3 months. We observed delayed wood morphogenesis in eki compared with WT, along with a more mechanically elastic cambial zone and radial compression of xylem cell size, indicating that rapid tissue differentiation is critical for cambial growth under mechanical stress. Furthermore, the touch-induced mechanosensory pathway was transcriptionally misregulated in eki, indicating that the ELIMÄKI locus is required to integrate the weight-growth feedback regulation. By studying this birch mutant, we were able to dissect vertical proprioception from the gravitropic response associated with reaction wood formation. Our study provides evidence for both local and systemic responses to mechanical stimuli during secondary plant development.

Phylogeny and highland adaptation of Chinese species in Allium section Daghestanica (Amaryllidaceae) revealed by transcriptome sequencing

Allium L. is one of the largest monocotyledonous genera with extensive distribution in the Northern Hemisphere. The fundamental phylogenies of Allium have been investigated using many morphological and molecular characters. However, the morphological characters may not agree with the molecular results in some Allium groups or sections (such as the Chinese Allium section Daghestanica), which may result in ambiguous species relationships and hinder further evolutionary and adaptive researches. Here, transcriptome sequences of the six Chinese endemics from Allium section Daghestanica were collected, with their single-copy genes (SCGs) were extracted. The interspecies relationships were analyzed using concatenation and coalescent methods. The branch-site model (BSM) was conducted to detect the positively selected genes (PSGs) in five highland species of this section. Based on 1644, 1281 and 1580 SCGs in flowers, leaves, and flowers-leaves combination respectively, a robust consistent and well-resolved phylogeny was generated from the concatenation method. Strong conflicts among individual gene trees were detected in the coalescent method, and morphological characters were incongruent with molecular relationships to some degree. Many PSGs were involved in responses of various stresses and stimuli (e.g. hypoxia, low temperature, aridity), DNA repair, metabolism, nutrient or energy intake, photosynthesis, and signal transduction. Our study revealed a clear interspecies relationship of Chinese endemics in Allium section Daghestanica and suggested that the discordance between morphological characters and molecular relationships might result from that the former are more susceptible to convergence compared with the latter. PSGs detected in our study may provide some insights into highland adaptation in Allium species.

A chromosome-level reference genome of the hornbeam, Carpinus fangiana

Betulaceae, the birch family, comprises six living genera and over 160 species, many of which are economically valuable. To deepen our knowledge of Betulaceae species, we have sequenced the genome of a hornbeam, Carpinus fangiana, which belongs to the most species-rich genus of the Betulaceae subfamily Coryloideae. Based on over 75 Gb (~200x) of high-quality next-generation sequencing data, we assembled a 386.19 Mb C. fangiana genome with contig N50 and scaffold N50 sizes of 35.32 kb and 1.91 Mb, respectively. Furthermore, 357.84 Mb of the genome was anchored to eight chromosomes using over 50 Gb (~130x) Hi-C sequencing data. Transcriptomes representing six tissues were sequenced to facilitate gene annotation, and over 5.50 Gb high-quality data were generated for each tissue. The structural annotation identified a total of 27,381 protein-coding genes in the assembled genome, of which 94.36% were functionally annotated. Additionally, 4,440 non-coding genes were predicted.

Measurement(s)	DNA • RNA • sequence_assembly • sequence feature annotation
Technology Type(s)	DNA sequencing • RNA sequencing • genome assembly • sequence annotation
Sample Characteristic - Organism	Carpinus fangiana

Machine-accessible metadata file describing the reported data: 10.6084/m9.figshare.11382390

The Terrestrial Carnivorous Plant Utricularia reniformis Sheds Light on Environmental and Life-Form Genome Plasticity

Utricularia belongs to Lentibulariaceae, a widespread family of carnivorous plants that possess ultra-small and highly dynamic nuclear genomes. It has been shown that the Lentibulariaceae genomes have been shaped by transposable elements expansion and loss, and multiple rounds of whole-genome duplications (WGD), making the family a platform for evolutionary and comparative genomics studies. To explore the evolution of Utricularia, we estimated the chromosome number and genome size, as well as sequenced the terrestrial bladderwort Utricularia reniformis (2n = 40, 1C = 317.1-Mpb). Here, we report a high quality 304 Mb draft genome, with a scaffold NG50 of 466-Kb, a BUSCO completeness of 87.8%, and 42,582 predicted genes. Compared to the smaller and aquatic U. gibba genome (101 Mb) that has a 32% repetitive sequence, the U. reniformis genome is highly repetitive (56%). The structural differences between the two genomes are the result of distinct fractionation and rearrangements after WGD, and massive proliferation of LTR-retrotransposons. Moreover, GO enrichment analyses suggest an ongoing gene birth-death-innovation process occurring among the tandem duplicated genes, shaping the evolution of carnivory-associated functions. We also identified unique patterns of developmentally related genes that support the terrestrial life-form and body plan of U. reniformis. Collectively, our results provided additional insights into the evolution of the plastic and specialized Lentibulariaceae genomes.

Negative feedback loop between BpAP1 and BpPI/BpDEF heterodimer in Betula platyphylla × B. pendula

MADS-box genes encode transcription factors involved in the control of many important developmental processes, especially the flower development of angiosperms. Analysis on gene regulatory relationship between MADS-box genes is useful for understanding the molecular mechanism of flower development. In this study, we focused on the regulatory relationship between MADS-box transcription factors APETALA1 (AP1) and PISTILLATA(PI)/DEFICIENS (DEF) in birch. We found that BpPI was an authentic target gene of BpAP1, and BpAP1 activated the expression of BpPI via directly binding to the CArG box motif. Functional analysis of BpPI showed that overexpression of BpPI may delay flowering via restricting flowering activators, in which BpAP1 was significantly down-regulated. We further investigated the regulatory of BpAP1 by BpPI, and found that BpPI could directly bind to the promoter of BpAP1 to restrict BpAP1 expression. In addition, we also found that BpPI could interact with its hypothetical partner BpDEF to co-regulate BpAP1 in birch. Our results suggested that overexpression of BpPI may delay flowering via restricting flowering activators, and there is a negative feedback loop between BpAP1 and BpPI/BpDEF heterodimer in birch. Our results will bring new evidences for further analysis of the molecular mechanism of flower formation in plants that produced unisexual flowers.

Characterization and T-DNA insertion sites identification of a multiple-branches mutant br in Betula platyphylla × Betula pendula

BACKGROUND

Plant architecture, which is mostly determined by shoot branching, plays an important role in plant growth and development. Thus, it is essential to explore the regulatory molecular mechanism of branching patterns based on the economic and ecological importance. In our previous work, a multiple-branches birch mutant br was identified from 19 CINNAMOYL-COENZYME A REDUCTASE 1 (CCR1)-overexpressed transgenic lines, and the expression patterns of differentially expressed genes in br were analyzed. In this study, we further explored some other characteristics of br, including plant architecture, wood properties, photosynthetic characteristics, and IAA and Zeatin contents. Meanwhile, the T-DNA insertion sites caused by the insertion of exogenous BpCCR1 in br were identified to explain the causes of the mutation phenotypes.

RESULTS

The mutant br exhibited slower growth, more abundant and weaker branches, and lower wood basic density and lignin content than BpCCR1 transgenic line (OE2) and wild type (WT). Compared to WT and OE2, br had high stomatal conductance (Gs), transpiration rate (Tr), but a low non-photochemical quenching coefficient (NPQ) and chlorophyll content. In addition, br displayed an equal IAA and Zeatin content ratio of main branches' apical buds to lateral branches' apical buds and high ratio of Zeatin to IAA content. Two T-DNA insertion sites caused by the insertion of exogenous BpCCR1 in br genome were found. On one site, chromosome 2 (Chr2), no known gene was detected on the flanking sequence. The other site was on Chr5, with an insertion of 388 bp T-DNA sequence, resulting in deletion of 107 bp 5' untranslated region (UTR) and 264 bp coding sequence (CDS) on CORONATINE INSENSITIVE 1 (BpCOII). In comparison with OE2 and WT, BpCOI1 was down-regulated in br, and the sensitivity of br to Methyl Jasmonate (MeJA) was abnormal.

CONCLUSIONS

Plant architecture, wood properties, photosynthetic characteristics, and IAA and Zeatin contents in main and lateral branches' apical buds changed in br over the study's time period. One T-DNA insertion was identified on the first exon of BpCOI1, which resulted in the reduction of BpCOI1 expression and abnormal perception to MeJA in br. These mutation phenotypes might be associated with a partial loss of BpCOI1 in birch.

Genome-Wide Analysis of NAC Gene Family in Betula pendula

NACs (NAM, ATAF1/2, and CUC2) are plant-specific transcription factors that play diverse roles in various plant developmental processes. In this study, we identified the NAC gene family in birch (Betula pendula) and further analyzed the function of BpNACs. Phylogenetic analysis reveals that the 114 BpNACs can be divided into seven subfamilies. We investigated the expression levels of these BpNACs in different tissues of birch including roots, xylem, leaves, and flowers, and the results showed that the BpNACs seem to be expressed higher in xylem and roots than leaves and flowers. In addition to tissue-specific expression analysis, we investigated the expression of BpNACs under low-temperature stress. A total of 21 BpNACs were differentially expressed under low-temperature stress, of which 17 were up-regulated, and four were down-regulated. Using the gene expression data, we reconstructed the gene co-expression network for the 21 low-temperature-responsive BpNACs. In conclusion, our results provide insight into the evolution of NAC genes in the B. pendula genome, and provide a basis for understanding the molecular mechanism for BpNAC-mediated cold responses in birch.

The avocado genome informs deep angiosperm phylogeny, highlights introgressive hybridization, and reveals pathogen-influenced gene space adaptation

The avocado is a nutritious, economically important fruit species that occupies an unresolved position near the earliest evolutionary branchings of flowering plants. Our nuclear genome sequences of Mexican and Hass variety avocados inform ancient evolutionary relationships and genome doublings and the admixed nature of Hass and provide a look at how pathogen interactions have shaped the avocado’s more recent genomic evolutionary history.

The avocado, Persea americana, is a fruit crop of immense importance to Mexican agriculture with an increasing demand worldwide. Avocado lies in the anciently diverged magnoliid clade of angiosperms, which has a controversial phylogenetic position relative to eudicots and monocots. We sequenced the nuclear genomes of the Mexican avocado race, P. americana var. drymifolia, and the most commercially popular hybrid cultivar, Hass, and anchored the latter to chromosomes using a genetic map. Resequencing of Guatemalan and West Indian varieties revealed that ∼39% of the Hass genome represents Guatemalan source regions introgressed into a Mexican race background. Some introgressed blocks are extremely large, consistent with the recent origin of the cultivar. The avocado lineage experienced 2 lineage-specific polyploidy events during its evolutionary history. Although gene-tree/species-tree phylogenomic results are inconclusive, syntenic ortholog distances to other species place avocado as sister to the enormous monocot and eudicot lineages combined. Duplicate genes descending from polyploidy augmented the transcription factor diversity of avocado, while tandem duplicates enhanced the secondary metabolism of the species. Phenylpropanoid biosynthesis, known to be elicited by Colletotrichum (anthracnose) pathogen infection in avocado, is one enriched function among tandems. Furthermore, transcriptome data show that tandem duplicates are significantly up- and down-regulated in response to anthracnose infection, whereas polyploid duplicates are not, supporting the general view that collections of tandem duplicates contribute evolutionarily recent “tuning knobs” in the genome adaptive landscapes of given species.