Integration of Mitoflash and Time-Series Transcriptomics Facilitates Energy Dynamics Tracking and Substrate Supply Analysis of Floral Thermogenesis in Lotus.

The high biosynthetic and energetic demands of floral thermogenesis render thermogenic plants the ideal systems to characterize energy metabolism in plants, but real-time tracking of energy metabolism in plant cells remains challenging. In this study, a new method was developed for tracking the mitochondrial energy metabolism at the single mitochondria level by real-time imaging of mitochondrial superoxide production (i.e., mitoflash). Using this method, we observed the increased mitoflash frequencies in the receptacles of Nelumbo nucifera Gaertn. at the thermogenic stages. This increase, combined with the higher expression of antioxidant response-related genes identified through time-series transcriptomics at the same stages, shows us a new regulatory mechanism for plant redox balance. Furthermore, we found that the upregulation of respiratory metabolism-related genes during the thermogenic stages not only correlates with changes in mitoflash frequency but also underscores the critical roles of these pathways in ensuring adequate substrate supply for thermogenesis. Metabolite analysis revealed that sugars are likely one of the substrates for thermogenesis and may be transported over long distances by sugar transporters. Taken together, our findings demonstrate that mitoflash is a reliable tool for tracking energy metabolism in thermogenic plants and contributes to our understanding of the regulatory mechanisms underlying floral thermogenesis.

Genetics of flight in spongy moths (Lymantria dispar ssp.): functionally integrated profiling of a complex invasive trait.

BACKGROUND: Flight can drastically enhance dispersal capacity and is a key trait defining the potential of exotic insect species to spread and invade new habitats. The phytophagous European spongy moths (ESM, Lymantria dispar dispar) and Asian spongy moths (ASM; a multi-species group represented here by L. d. asiatica and L. d. japonica), are globally invasive species that vary in adult female flight capability-female ASM are typically flight capable, whereas female ESM are typically flightless. Genetic markers of flight capability would supply a powerful tool for flight profiling of these species at any intercepted life stage. To assess the functional complexity of spongy moth flight and to identify potential markers of flight capability, we used multiple genetic approaches aimed at capturing complementary signals of putative flight-relevant genetic divergence between ESM and ASM: reduced representation genome-wide association studies, whole genome sequence comparisons, and developmental transcriptomics. We then judged the candidacy of flight-associated genes through functional analyses aimed at addressing the proximate demands of flight and salient features of the ecological context of spongy moth flight evolution. RESULTS: Candidate gene sets were typically non-overlapping across different genetic approaches, with only nine gene annotations shared between any pair of approaches. We detected an array of flight-relevant functional themes across gene sets that collectively suggest divergence in flight capability between European and Asian spongy moth lineages has coincided with evolutionary differentiation in multiple aspects of flight development, execution, and surrounding life history. Overall, our results indicate that spongy moth flight evolution has shaped or been influenced by a large and functionally broad network of traits. CONCLUSIONS: Our study identified a suite of flight-associated genes in spongy moths suited to exploration of the genetic architecture and evolution of flight, or validation for flight profiling purposes. This work illustrates how complementary genetic approaches combined with phenotypically targeted functional analyses can help to characterize genetically complex traits.

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.

Does urbanisation lead to parallel demographic shifts across the world in a cosmopolitan plant?

Urbanisation is occurring globally, leading to dramatic environmental changes that are altering the ecology and evolution of species. In particular, the expansion of human infrastructure and the loss and fragmentation of natural habitats in cities is predicted to increase genetic drift and reduce gene flow by reducing the size and connectivity of populations. Alternatively, the 'urban facilitation model' suggests that some species will have greater gene flow into and within cities leading to higher diversity and lower differentiation in urban populations. These alternative hypotheses have not been contrasted across multiple cities. Here, we used the genomic data from the GLobal Urban Evolution project (GLUE), to study the effects of urbanisation on non-adaptive evolutionary processes of white clover (Trifolium repens) at a global scale. We found that white clover populations presented high genetic diversity and no evidence of reduced Ne linked to urbanisation. On the contrary, we found that urban populations were less likely to experience a recent decrease in effective population size than rural ones. In addition, we found little genetic structure among populations both globally and between urban and rural populations, which showed extensive gene flow between habitats. Interestingly, white clover displayed overall higher gene flow within urban areas than within rural habitats. Our study provides the largest comprehensive test of the demographic effects of urbanisation. Our results contrast with the common perception that heavily altered and fragmented urban environments will reduce the effective population size and genetic diversity of populations and contribute to their isolation.

High-quality genome assembly enables prediction of allele-specific gene expression in hybrid poplar.

Poplar (Populus) is a well-established model system for tree genomics and molecular breeding, and hybrid poplar is widely used in forest plantations. However, distinguishing its diploid homologous chromosomes is difficult, complicating advanced functional studies on specific alleles. In this study, we applied a trio-binning design and PacBio high-fidelity long-read sequencing to obtain haplotype-phased telomere-to-telomere genome assemblies for the 2 parents of the well-studied F1 hybrid "84K" (Populus alba × Populus tremula var. glandulosa). Almost all chromosomes, including the telomeres and centromeres, were completely assembled for each haplotype subgenome apart from 2 small gaps on one chromosome. By incorporating information from these haplotype assemblies and extensive RNA-seq data, we analyzed gene expression patterns between the 2 subgenomes and alleles. Transcription bias at the subgenome level was not uncovered, but extensive-expression differences were detected between alleles. We developed machine-learning (ML) models to predict allele-specific expression (ASE) with high accuracy and identified underlying genome features most highly influencing ASE. One of our models with 15 predictor variables achieved 77% accuracy on the training set and 74% accuracy on the testing set. ML models identified gene body CHG methylation, sequence divergence, and transposon occupancy both upstream and downstream of alleles as important factors for ASE. Our haplotype-phased genome assemblies and ML strategy highlight an avenue for functional studies in Populus and provide additional tools for studying ASE and heterosis in hybrids.

Efficient genomics-based 'end-to-end' selective tree breeding framework.

Since their initiation in the 1950s, worldwide selective tree breeding programs followed the recurrent selection scheme of repeated cycles of selection, breeding (mating), and testing phases and essentially remained unchanged to accelerate this process or address environmental contingencies and concerns. Here, we introduce an "end-to-end" selective tree breeding framework that: (1) leverages strategically preselected GWAS-based sequence data capturing trait architecture information, (2) generates unprecedented resolution of genealogical relationships among tested individuals, and (3) leads to the elimination of the breeding phase through the utilization of readily available wind-pollinated (OP) families. Individuals' breeding values generated from multi-trait multi-site analysis were also used in an optimum contribution selection protocol to effectively manage genetic gain/co-ancestry trade-offs and traits' correlated response to selection. The proof-of-concept study involved a 40-year-old spruce OP testing population growing on three sites in British Columbia, Canada, clearly demonstrating our method's superiority in capturing most of the available genetic gains in a substantially reduced timeline relative to the traditional approach. The proposed framework is expected to increase the efficiency of existing selective breeding programs, accelerate the start of new programs for ecologically and environmentally important tree species, and address climate-change caused biotic and abiotic stress concerns more effectively.

Unraveling the transcriptional features and gene expression networks of pathogenic and saprotrophic Ophiostoma species during the infection of Ulmus americana.

American elm (Ulmus americana), highly prized for its ornamental value, has suffered two successive outbreaks of Dutch elm disease (DED) caused by ascomycete fungi belonging to the genus Ophiostoma. To identify the genes linked to the pathogenicity of different species and lineages of Ophiostoma, we inoculated 2-year-old U. americana saplings with six strains representing three species of DED fungi, and one strain of the saprotroph Ophiostoma quercus. Differential expression analyses were performed following RNA sequencing of fungal transcripts recovered at 3- and 10-days post-infection. Based on a total of 8,640 Ophiostoma genes, we observed a difference in fungal gene expression depending on the strain inoculated and the time of incubation in host tissue. Some genes overexpressed in the more virulent strains of Ophiostoma encode hydrolases that possibly act synergistically. A mutant of Ophiostoma novo-ulmi in which the gene encoding the ogf1 transcription factor had been deleted did not produce transcripts for the gene encoding the hydrophobin cerato-ulmin and was less virulent. Weighted gene correlation network analyses identified several candidate pathogenicity genes distributed among 13 modules of interconnected genes.IMPORTANCEOphiostoma is a genus of cosmopolitan fungi that belongs to the family Ophiostomataceae and includes the pathogens responsible for two devastating pandemics of Dutch elm disease (DED). As the mechanisms of action of DED agents remain unclear, we carried out the first comparative transcriptomic study including representative strains of the three Ophiostoma species causing DED, along with the phylogenetically close saprotrophic species Ophiostoma quercus. Statistical analyses of the fungal transcriptomes recovered at 3 and 10 days following infection of Ulmus americana saplings highlighted several candidate genes associated with virulence and host-pathogen interactions wherein each strain showed a distinct transcriptome. The results of this research underscore the importance of investigating the transcriptional behavior of different fungal taxa to understand their pathogenicity and virulence in relation to the timeline of infection.

Is there hybridization between 2 species of the same genus in sympatry?-The genetic relationships between Anoplophora glabripennis, Anoplophora chinensis, and putative hybrids.

Anoplophora glabripennis (Asian longhorn beetle, ALB) and Anoplophora chinensis (Citrus longhorn beetle, CLB) are native forest pests in China; they have become important international quarantine pests. They are found using the same Salix aureo-pendula host tree of Cixi, Zhejiang province, China. On this host tree, we collected additional beetles that appeared to be morphologically intermediate between ALB and CLB. By using a stereoscope, we observed that there were several bumps on the base of the elytra, which was inconsistent with ALB, which typically has a smooth elytral base, but was more like CLB, which has numerous short tubercles on the elytral base. Given their sympatry and intermediate morphology, we hypothesized that these may represent ALB × CLB hybrids. We studied the genomic profiles for 46 samples (ALB, CLB, and putative hybrids) using genotyping-by-sequencing (GBS) providing a reduced representation of the entire genome. Employing principal component analyses on the 163 GBS-derived single nucleotide polymorphism data, we found putative hybrids tightly clustered with ALB, but genetically distinct from the CLB individuals. Therefore, our initial hybrid hypothesis was not supported by genomic data. Further, while mating experiments between adult ALB and CLB were successful in 4 separate years (2017, 2018, 2020, and 2021), and oviposition behavior was observed, no progeny was produced. Having employed population genomic analysis and biological hybridization experiments, we conclude that the putative hybrids represent newly discovered morphological variants within ALB. Our approach further confirmed the advantage of genome-wide information for Anoplophora species assignment in certain ambiguous classification cases.

Unraveling the evolutionary dynamics of the TPS gene family in land plants.

Terpenes and terpenoids are key natural compounds for plant defense, development, and composition of plant oil. The synthesis and accumulation of a myriad of volatile terpenoid compounds in these plants may dramatically alter the quality and flavor of the oils, which provide great commercial utilization value for oil-producing plants. Terpene synthases (TPSs) are important enzymes responsible for terpenic diversity. Investigating the differentiation of the TPS gene family could provide valuable theoretical support for the genetic improvement of oil-producing plants. While the origin and function of TPS genes have been extensively studied, the exact origin of the initial gene fusion event - it occurred in plants or microbes - remains uncertain. Furthermore, a comprehensive exploration of the TPS gene differentiation is still pending. Here, phylogenetic analysis revealed that the fusion of the TPS gene likely occurred in the ancestor of land plants, following the acquisition of individual C- and N- terminal domains. Potential mutual transfer of TPS genes was observed among microbes and plants. Gene synteny analysis disclosed a differential divergence pattern between TPS-c and TPS-e/f subfamilies involved in primary metabolism and those (TPS-a/b/d/g/h subfamilies) crucial for secondary metabolites. Biosynthetic gene clusters (BGCs) analysis suggested a correlation between lineage divergence and potential natural selection in structuring terpene diversities. This study provides fresh perspectives on the origin and evolution of the TPS gene family.

Haplotype-resolved genome assembly of Coriaria nepalensis a non-legume nitrogen-fixing shrub.

Coriaria nepalensis Wall. (Coriariaceae) is a nitrogen-fixing shrub which forms root nodules with the actinomycete Frankia. Oils and extracts of C. nepalensis have been reported to be bacteriostatic and insecticidal, and C. nepalensis bark provides a valuable tannin resource. Here, by combining PacBio HiFi sequencing and Hi-C scaffolding techniques, we generated a haplotype-resolved chromosome-scale genome assembly for C. nepalensis. This genome assembly is approximately 620 Mb in size with a contig N50 of 11 Mb, with 99.9% of the total assembled sequences anchored to 40 pseudochromosomes. We predicted 60,862 protein-coding genes of which 99.5% were annotated from databases. We further identified 939 tRNAs, 7,297 rRNAs, and 982 ncRNAs. The chromosome-scale genome of C. nepalensis is expected to be a significant resource for understanding the genetic basis of root nodulation with Frankia, toxicity, and tannin biosynthesis.

Chromosome-scale genome assembly and insights into the metabolome and gene regulation of leaf color transition in an important oak species, Quercus dentata.

Quercus dentata Thunb., a dominant forest tree species in northern China, has significant ecological and ornamental value due to its adaptability and beautiful autumn coloration, with color changes from green to yellow into red resulting from the autumnal shifts in leaf pigmentation. However, the key genes and molecular regulatory mechanisms for leaf color transition remain to be investigated. First, we presented a high-quality chromosome-scale assembly for Q. dentata. This 893.54 Mb sized genome (contig N50 = 4.21 Mb, scaffold N50 = 75.55 Mb; 2n = 24) harbors 31 584 protein-coding genes. Second, our metabolome analyses uncovered pelargonidin-3-O-glucoside, cyanidin-3-O-arabinoside, and cyanidin-3-O-glucoside as the main pigments involved in leaf color transition. Third, gene co-expression further identified the MYB-bHLH-WD40 (MBW) transcription activation complex as central to anthocyanin biosynthesis regulation. Notably, transcription factor (TF) QdNAC (QD08G038820) was highly co-expressed with this MBW complex and may regulate anthocyanin accumulation and chlorophyll degradation during leaf senescence through direct interaction with another TF, QdMYB (QD01G020890), as revealed by our further protein-protein and DNA-protein interaction assays. Our high-quality genome assembly, metabolome, and transcriptome resources further enrich Quercus genomics and will facilitate upcoming exploration of ornamental values and environmental adaptability in this important genus.

Gapless genome assembly of azalea and multi-omics investigation into divergence between two species with distinct flower color.

The genus Rhododendron (Ericaceae), with more than 1000 species highly diverse in flower color, is providing distinct ornamental values and a model system for flower color studies. Here, we investigated the divergence between two parental species with different flower color widely used for azalea breeding. Gapless genome assembly was generated for the yellow-flowered azalea, Rhododendron molle. Comparative genomics found recent proliferation of long terminal repeat retrotransposons (LTR-RTs), especially Gypsy, has resulted in a 125 Mb (19%) genome size increase in species-specific regions, and a significant amount of dispersed gene duplicates (13 402) and pseudogenes (17 437). Metabolomic assessment revealed that yellow flower coloration is attributed to the dynamic changes of carotenoids/flavonols biosynthesis and chlorophyll degradation. Time-ordered gene co-expression networks (TO-GCNs) and the comparison confirmed the metabolome and uncovered the specific gene regulatory changes underpinning the distinct flower pigmentation. B3 and ERF TFs were found dominating the gene regulation of carotenoids/flavonols characterized pigmentation in R. molle, while WRKY, ERF, WD40, C2H2, and NAC TFs collectively regulated the anthocyanins characterized pigmentation in the red-flowered R simsii. This study employed a multi-omics strategy in disentangling the complex divergence between two important azaleas and provided references for further functional genetics and molecular breeding.

Building resiliency in conifer forests: Interior spruce crosses among weevil resistant and susceptible parents produce hybrids appropriate for multi-trait selection.

Tree planting programs now need to consider climate change increasingly, therefore, the resistance to pests plays an essential role in enabling tree adaptation to new ranges through tree population movement. The weevil Pissodes strobi (Peck) is a major pest of spruces and substantially reduces lumber quality. We revisited a large Interior spruce provenance/progeny trial (2,964 genotypes, 42 families) of varying susceptibility, established in British Columbia. We employed multivariate mixed linear models to estimate covariances between, and genetic control of, juvenile height growth and resistance traits. We performed linear regressions and ordinal logistic regressions to test for impact of parental origin on growth and susceptibility to the pest, respectively. A significant environmental component affected the correlations between resistance and height, with outcomes dependent on families. Parents sourced from above 950 m a.s.l. elevation negatively influenced host resistance to attacks, probably due to higher P. engelmannii proportion. For the genetic contribution of parents sourced from above 1,200 m a.s.l., however, we found less attack severity, probably due to a marked mismatch in phenologies. This clearly highlights that interspecific hybrid status might be a good predictor for weevil attacks and delineates the boundaries of successful spruce population movement. Families resulting from crossing susceptible parents generally showed fast-growing trees were the most affected by weevil attacks. Such results indicate that interspecific 'hybrids' with a higher P. glauca ancestry might be genetically better equipped with an optimized resource allocation between defence and growth and might provide the solution for concurrent improvement in resistance against weevil attacks, whilst maintaining tree productivity.

Potential allopolyploid origin of Ericales revealed with gene-tree reconciliation.

Few incidents of ancient allopolyploidization (polyploidization by hybridization or merging diverged genomes) were previously revealed, although there is significant evidence for the accumulation of whole genome duplications (WGD) in plants. Here, we focused on Ericales, one of the largest and most diverse angiosperm orders with significant ornamental and economic value. Through integrating 24 high-quality whole genome data selected from ~ 200 Superasterids genomes/species and an algorithm of topology-based gene-tree reconciliation, we explored the evolutionary history of in Ericales with ancient complex. We unraveled the allopolyploid origin of Ericales and detected extensive lineage-specific gene loss following the polyploidization. Our study provided a new hypothesis regarding the origin of Ericales and revealed an instructive perspective of gene loss as a pervasive source of genetic variation and adaptive phenotypic diversity in Ericales.

Genome-scale phylogeography resolves the native population structure of the Asian longhorned beetle, Anoplophora glabripennis (Motschulsky).

Human-assisted movement has allowed the Asian longhorned beetle (ALB, Anoplophora glabripennis (Motschulsky)) to spread beyond its native range and become a globally regulated invasive pest. Within its native range of China and the Korean peninsula, human-mediated dispersal has also caused cryptic translocation of insects, resulting in population structure complexity. Previous studies used genetic methods to detangle this complexity but were unable to clearly delimit native populations which is needed to develop downstream biosurveillance tools. We used genome-wide markers to define historical population structure in native ALB populations and contemporary movement between regions. We used genotyping-by-sequencing to generate 6102 single-nucleotide polymorphisms (SNPs) and amplicon sequencing to genotype 53 microsatellites. In total, we genotyped 712 individuals from ALB's native distribution. We observed six distinct population clusters among native ALB populations, with a clear delineation between northern and southern groups. Most of the individuals from South Korea were distinct from populations in China. Our results also indicate historical divergence among populations and suggest limited large-scale admixture, but we did identify a restricted number of cases of contemporary movement between regions. We identified SNPs under selection and describe a clinal allele frequency pattern in a missense variant associated with glycerol kinase, an important enzyme in the utilization of an insect cryoprotectant. We further demonstrate that small numbers of SNPs can assign individuals to geographic regions with high probability, paving the way for novel ALB biosurveillance tools.

Lilac (Syringa oblata) genome provides insights into its evolution and molecular mechanism of petal color change.

Color change during flower opening is common; however, little is understood on the biochemical and molecular basis related. Lilac (Syringa oblata), a well-known woody ornamental plant with obvious petal color changes, is an ideal model. Here, we presented chromosome-scale genome assembly for lilac, resolved the flavonoids metabolism, and identified key genes and potential regulatory networks related to petal color change. The genome assembly is 1.05 Gb anchored onto 23 chromosomes, with a BUSCO score of 96.6%. Whole-genome duplication (WGD) event shared within Oleaceae was revealed. Metabolome quantification identified delphinidin-3-O-rutinoside (Dp3Ru) and cyanidin-3-O-rutinoside (Cy3Ru) as the major pigments; gene co-expression networks indicated WRKY an essential regulation factor at the early flowering stage, ERF more important in the color transition period (from violet to light nearly white), while the MBW complex participated in the entire process. Our results provide a foundation for functional study and molecular breeding in lilac.

SubPhaser: a robust allopolyploid subgenome phasing method based on subgenome-specific k-mers.

With advanced sequencing technology, dozens of complex polyploid plant genomes have been characterized. However, for many polyploid species, their diploid ancestors are unknown or extinct, making it impossible to unravel the subgenomes and genome evolution directly. We developed a novel subgenome-phasing algorithm, SubPhaser, specifically designed for a neoallopolyploid or a homoploid hybrid. SubPhaser first searches for the subgenome-specific sequence (k-mer), then assigns homoeologous chromosomes into subgenomes, and further provides tools to annotate and investigate specific sequences. SubPhaser works well on neoallopolyploids and homoploid hybrids containing subgenome-specific sequences like wheat, but fails on autopolyploids lacking subgenome-specific sequences like alfalfa, indicating that SubPhaser can phase neoallopolyploid/homoploid hybrids with high accuracy, sensitivity and performance. This highly accurate, highly sensitive, ancestral data free chromosome phasing algorithm, SubPhaser, offers significant application value for subgenome phasing in neoallopolyploids and homoploid hybrids, and for the subsequent exploration of genome evolution and related genetic/epigenetic mechanisms.

UV-B and UV-C radiation trigger both common and distinctive signal perceptions and transmissions in Pinus tabuliformis Carr.

In plants, ultraviolet (UV)-light is an important driver for growth and natural distribution, and is also a valuable tool for manipulating productivity as well as biotic interactions. Understanding of plant responses to different UV radiation is sparse, especially from a systems biology perspective and particularly for conifers. Here, we evaluated the physiological and transcriptomic responses to the short-term application of high-irradiance UV-B and UV-C waves on Pinus tabuliformis Carr., a major conifer in Northern China. By undertaking time-ordered gene coexpression network analyses and network comparisons incorporating physiological traits and gene expression variation, we uncovered communalities but also differences in P. tabuliformis responses to UV-B and UV-C. Both types of spectral bands caused a significant inhibition of photosynthesis, and conversely, the improvement of antioxidant capacity, flavonoid production and signaling pathways related to stress resistance, indicating a clear switch from predominantly primary metabolism to enhanced defensive metabolism in pine. We isolated distinct subnetworks for photoreceptor-mediated signal transduction, maximum quantum efficiency of photosystem II (Fv/Fm) regulation and flavonoid biosynthesis in response to UV-B and UV-C radiation. From these subnetworks, we further identified phototropins as potentially important elements in both UV-B and UV-C signaling and, for the first time, suggesting peptide hormones to be involved in promoting flavonoid biosynthesis against UV-B, while these hormones seem not to be implicated in the defense against UV-C exposure. The present study employed an effective strategy for disentangling the complex physiological and genetic regulatory mechanisms in a nonmodel plant species, and thus, provides a suitable reference for future functional evaluations and artificial UV-light mediated growing strategies in plant production.

Independent Evolution Has Led to Distinct Genomic Signatures in Dutch Elm Disease-Causing Fungi and Other Vascular Wilts-Causing Fungal Pathogens.

Vascular wilts are important diseases caused by plant pathogenic fungi that result in the rapid death of their plant hosts. This is due to a systemic defense mechanism whereby the plant induces the compartmentalization of the infected vascular system in order to reduce the propagation of the fungus. The ascomycete class Sordariomycetes contains several species that cause vascular wilts in diverse plant hosts, and they can be classified into four taxonomic orders. The genetic mechanisms of pathogenesis have already been investigated in Fusarium and Verticillium species, but they have not yet been compared with other well-known wilt-causing species, especially fungi causing oak wilt or Dutch elm disease (DED). Here we analyzed 20 whole genome assemblies of wilt-causing fungi together with 56 other species using phylogenetic approaches to trace expansions and contractions of orthologous gene families and gene classes related to pathogenicity. We found that the wilt-causing pathogens evolved seven times, experiencing the largest fold changes in different classes of genes almost every time. However, some similarities exist across groups of wilt pathogens, particularly in Microascales and Ophiostomatales, and these include the common gains and losses of genes that make up secondary metabolite clusters (SMC). DED pathogens do not experience large-scale gene expansions, with most of the gene classes, except for some SMC families, reducing in number. We also found that gene family expansions in the most recent common ancestors of wilt pathogen groups are enriched for carbohydrate metabolic processes. Our study shows that wilt-causing species evolve primarily through distinct changes in their repertoires of pathogenicity-related genes and that there is the potential importance of carbohydrate metabolism genes for regulating osmosis in those pathogens that penetrate the plant vascular system.

Centromere-Specific Retrotransposons and Very-Long-Chain Fatty Acid Biosynthesis in the Genome of Yellowhorn (Xanthoceras sorbifolium, Sapindaceae), an Oil-Producing Tree With Significant Drought Resistance.

In-depth genome characterization is still lacking for most of biofuel crops, especially for centromeres, which play a fundamental role during nuclear division and in the maintenance of genome stability. This study applied long-read sequencing technologies to assemble a highly contiguous genome for yellowhorn (Xanthoceras sorbifolium), an oil-producing tree, and conducted extensive comparative analyses to understand centromere structure and evolution, and fatty acid biosynthesis. We produced a reference-level genome of yellowhorn, ∼470 Mb in length with ∼95% of contigs anchored onto 15 chromosomes. Genome annotation identified 22,049 protein-coding genes and 65.7% of the genome sequence as repetitive elements. Long terminal repeat retrotransposons (LTR-RTs) account for ∼30% of the yellowhorn genome, which is maintained by a moderate birth rate and a low removal rate. We identified the centromeric regions on each chromosome and found enrichment of centromere-specific retrotransposons of LINE1 and Gypsy in these regions, which have evolved recently (∼0.7 MYA). We compared the genomes of three cultivars and found frequent inversions. We analyzed the transcriptomes from different tissues and identified the candidate genes involved in very-long-chain fatty acid biosynthesis and their expression profiles. Collinear block analysis showed that yellowhorn shared the gamma (γ) hexaploidy event with Vitis vinifera but did not undergo any further whole-genome duplication. This study provides excellent genomic resources for understanding centromere structure and evolution and for functional studies in this important oil-producing plant.