Inferring duplications, losses, transfers and incomplete lineage sorting with nonbinary species trees

ERCnet: Phylogenomic Prediction of Interaction Networks in the Presence of Gene Duplication

Assigning gene function from genome sequences is a rate-limiting step in molecular biology research. A protein's position within an interaction network can potentially provide insights into its molecular mechanisms. Phylogenetic analysis of evolutionary rate covariation (ERC) in protein sequence has been shown to be effective for large-scale prediction of functional relationships and interactions. However, gene duplication, gene loss, and other sources of phylogenetic incongruence are barriers for analyzing ERC on a genome-wide basis. Here, we developed ERCnet, a bioinformatic program designed to overcome these challenges, facilitating efficient all-versus-all ERC analyses for large protein sequence datasets. We simulated proteome datasets and found that ERCnet achieves combined false positive and negative error rates well below 10% and that our novel "branch-by-branch" length measurements outperforms "root-to-tip" approaches in most cases, offering a valuable new strategy for performing ERC. We also compiled a sample set of 35 angiosperm genomes to test the performance of ERCnet on empirical data, including its sensitivity to user-defined analysis parameters such as input dataset size and branch-length measurement strategy. We investigated the overlap between ERCnet runs with different species samples to understand how species number and composition affect predicted interactions and to identify the protein sets that consistently exhibit ERC across angiosperms. Our systematic exploration of the performance of ERCnet provides a roadmap for design of future ERC analyses to predict functional interactions in a wide array of genomic datasets. ERCnet code is freely available at https://github.com/EvanForsythe/ERCnet.

Adaptive evolution of sesquiterpene deoxyphomenone in mycoparasitism by Hansfordia pulvinata associated with horizontal gene transfer from Aspergillus species

Leaf mold caused by the ascomycete fungus Cladosporium fulvum is a devastating disease of tomato plants. The mycoparasitic fungus Hansfordia pulvinata is an effective biocontrol agent that parasitizes C. fulvum hyphae on leaves and secretes 13-deoxyphomenone, an eremophilane-type sesquiterpene, which was also identified as a sporulation-inducing factor in Aspergillus oryzae. Here, we identified deoxyphomenone biosynthesis (DPH) gene clusters conserved in both H. pulvinata and Aspergillus section Flavi, including A. oryzae and A. flavus. Functional disruption of DPH1 orthologous genes encoding sesquiterpene cyclase in H. pulvinata, A. oryzae, and its close relative A. flavus revealed that deoxyphomenone in H. pulvinata had exogenic antifungal activity against C. fulvum and controlled endogenic sporulation in Aspergillus species. Complete DPH clusters, highly similar to those in H. pulvinata, were exclusive to Aspergillus section Flavi, while species in other Aspergillus sections contained fragmented DPH clusters. A comparative genomics analysis revealed that these DPH gene clusters share a common origin and are horizontally transferred from an ancestor of Aspergillus to H. pulvinata. Our results suggest that after horizontal transfer, H. pulvinata maintained the DPH cluster as the inhibitory effect of deoxyphomenone on spore germination and mycelial growth contributed to its mycoparasitism on the host fungus C. fulvum.

IMPORTANCE

Tomato leaf mold disease caused by C. fulvum poses a significant economic threat to tomato production globally. Breeders have developed tomato cultivars with Cf resistance genes. C. fulvum frequently evolves new races that overcome these genetic defenses, complicating control efforts. Additionally, the pathogen has developed resistance to chemical fungicides, prompting the need for sustainable alternatives like biocontrol agents. The mycoparasitic fungus H. pulvinata is crucial as an effective agent against C. fulvum. Clarifying the mechanism of mycoparasitism is significant, as it enhances its application as a biocontrol agent against plant pathogens. This study revealed how H. pulvinata produces deoxyphomenone, an antifungal compound, through horizontal gene transfer from Aspergillus species. It is hypothesized that mycoparasitism could be one of the mechanisms that facilitated horizontal gene transfer between fungi. These insights facilitate the development of eco-friendly, sustainable agricultural practices by reducing dependence on chemical fungicides and promoting natural pathogen control methods.

Identification and expression responses of TCP gene family in Opisthopappus taihangensis under abiotic stress

The TCP gene family plays pivotal roles in the development and abiotic stress responses of plants; however, no data has been provided for this gene family in Opisthopappus taihangensis. Based on O. taihangensis genome, 14 TCP genes were identified and divided into two classes (I and II). After tandem and segmental duplication/whole-genome duplication (WGD), more loss and less gain events of OtTCPs occurred, which might be related with the underwent purifying selection during the evolution. The conserved motifs and structures of OtTCP genes contained light response, growth and development, hormone response, and stress-related cis-acting elements. Different OtTCP genes, even duplicated gene pairs, could be expressed in different tissues, which implied that OtTCP genes had diverse function. Among OtTCPs, OtTCP4, 9 and 11 of CYC clade (Class II) presented a relative wide expression pattern with no or one intron. The three TCP genes could be regarded as important candidate factors for O. taihangensis in growth, development and stress response. These results provided some clues and references for the further in-depth exploration of O. taihangensis resistance mechanisms, as well as those of other unique eco-environment plants.

Systematic characterization of the bZIP gene family in Colletotrichum siamense and functional analysis of three family members

The basic leucine zipper (bZIP) transcription factors (TFs) play important roles in many physiological processes of plant-pathogenic fungi, especially concerning fungal development, fungicide resistance, and pathogenicity. Colletotrichum siamense is the predominant species causing Colletotrichum leaf disease (CLD) in rubber trees. However, little is known about the bZIP genes in C. siamense. In this study, 25 bZIP genes were systematically identified in the genome of C. siamense, and molecular features were characterized. Evolutionarily, the CsbZIP genes were divided into 11 groups, with the members in the same group sharing similar gene structures and conserved protein motif organizations. Furthermore, protein-protein interaction (PPI) analysis revealed that 15 bZIP proteins had functional partners in common or interacted with other CsbZIP proteins. Additionally, the expression of 23 CsbZIP genes changed in response to the antifungal chemicals melatonin, prochloraz, and thymol, and the genes could be divided into three clusters based on their expression patterns. Finally, gene deletion mutants of CsbZIP01/09/17 were constructed and functional analysis indicated that these genes operated as important regulators of mycelial growth, fungicide resistance, ergosterol biosynthesis, and virulence in C. siamense. This study provided the foundations crucial for further investigation of the functions of CsbZIP TFs in fungicide resistance and virulence.

Deciphering the landscape and evolutionary trajectory of NLR immune receptors in Dioscorea alata

Dioscorea alata, a key tuber crop for global food security, is threatened by anthracnose disease caused by Colletotrichum gloeosporioides. However, identification of functional resistance genes against C. gloeosporioides in D. alata is challenging due to low flowering and hybridization efficiency of this plant. Nucleotide-binding leucine-rich repeat (NLR) genes constitute the largest group of plant disease resistance genes, from which functional genes against diverse pathogens across various crops have been cloned. In this study, a comprehensive genome-wide analysis identified 346 NLR genes from D. alata, including one RNL and 345 CNLs. These NLRs were unequally distributed on 20 chromosomes, with chromosome 3 harboring the highest number (78 NLR genes). The majority of NLR genes (91%) were located in multigene clusters, implying that tandem or proximal duplication was the primary driving force for NLR gene expansion in D. alata. Comparative analysis of Dioscoreaceae species revealed high variability and differential expansion patterns of NLR genes. In addition, transcriptome profiling of D. alata post-infection with C. gloeosporioides identified 12 differentially expressed NLR genes. In summary, this study sheds new light on the genetic architecture and evolutionary dynamics of D. alata NLR genes, offering valuable insights for cloning functional genes against C. gloeosporioides.

Evolutionary analysis of TIR- and non-TIR-NBS-LRR disease resistance genes in wild strawberries

Introduction

NBS-LRR genes (NLRs) are the most extensive category of plant resistance genes (R genes) and play a crucial role in pathogen defense. Understanding the diversity and evolutionary dynamics of NLRs in different plant species is essential for improving disease resistance. This study investigates the NLR gene family in eight diploid wild strawberry species to explore their structural characteristics, evolutionary relationships, and potential for enhancing disease resistance.

Methods

We conducted a comprehensive genome-wide identification and structural analysis of NLRs across eight diploid wild strawberry species. Phylogenetic analysis was performed to examine the relationships between TIR-NLRs (TNLs), Non-TIR-NLRs (non-TNLs), CC-NLRs (CNLs), and RPW8-NLRs (RNLs). Gene structures were compared, and gene expression was profiled across different NLR subfamilies. Additionally, in vitro leaf inoculation assays with Botrytis cinerea were performed to assess the resistance of various strawberry species.

Results

Our analysis revealed that non-TNLs constitute over 50% of the NLR gene family in all eight strawberry species, surpassing the proportion of TNLs. Phylogenetic analysis showed that TNLs diverged into two subclades: one grouping with CNLs and the other closely related to RNLs. A significantly higher number of non-TNLs were under positive selection compared to TNLs, indicating their rapid diversification. Gene structure analysis demonstrated that non-TNLs have shorter gene structures than TNLs and exhibit higher expression levels, particularly RNLs. Notably, non-TNLs showed dominant expression under both normal and infected conditions. In vitro leaf inoculation assays revealed that Fragaria pentaphylla and Fragaria nilgerrensis, which have the highest proportion of non-TNLs, exhibited significantly greater resistance to Botrytis cinerea compared to Fragaria vesca, which has the lowest proportion of non-TNLs.

Discussion

The findings of this study provide important insights into the evolutionary dynamics of NLRs in strawberries, particularly the significant role of non-TNLs in pathogen defense. The rapid diversification and higher expression levels of non-TNLs suggest their potential contribution to enhanced disease resistance. This research highlights the value of non-TNLs in strawberry breeding programs aimed at improving resistance to pathogens such as Botrytis cinerea.

R2R3-MYB Gene Family in Coptis teeta Wall.: Genome-Wide Identification, Phylogeny, Evolutionary Expansion, and Expression Analyses during Floral Development

The R2R3-MYB gene family represents a widely distributed class of plant transcription factors. This gene family plays an important role in many aspects of plant growth and development. However, the characterization of R2R3-MYB genes present in the genome of Coptis teeta has not been reported. Here, we describe the bioinformatic identification and characterization of 88 R2R3-MYB genes in this species, and the identification of members of the R2R3-MYB gene family in species within the order Ranales most closely related to Coptis teeta. The CteR2R3-MYB genes were shown to exhibit a higher degree of conservation compared to those of A. thaliana, as evidenced by phylogeny, conserved motifs, gene structure, and replication event analyses. Cis-acting element analysis confirmed the involvement of CteR2R3-MYB genes in a variety of developmental processes, including growth, cell differentiation, and reproduction mediated by hormone synthesis. In addition, through homology comparisons with the equivalent gene family in A. thaliana, protein regulatory network prediction and transcriptome data analysis of floral organs across three time periods of flower development, 17 candidate genes were shown to exhibit biased expression in two floral phenotypes of C. teeta. This suggests their potential involvement in floral development (anther development) in this species.

Reclassification of Botryococcus braunii chemical races into separate species based on a comparative genomics analysis

The colonial green microalga Botryococcus braunii is well known for producing liquid hydrocarbons that can be utilized as biofuel feedstocks. B. braunii is taxonomically classified as a single species made up of three chemical races, A, B, and L, that are mainly distinguished by the hydrocarbons produced. We previously reported a B race draft nuclear genome, and here we report the draft nuclear genomes for the A and L races. A comparative genomic study of the three B. braunii races and 14 other algal species within Chlorophyta revealed significant differences in the genomes of each race of B. braunii. Phylogenomically, there was a clear divergence of the three races with the A race diverging earlier than both the B and L races, and the B and L races diverging from a later common ancestor not shared by the A race. DNA repeat content analysis suggested the B race had more repeat content than the A or L races. Orthogroup analysis revealed the B. braunii races displayed more gene orthogroup diversity than three closely related Chlamydomonas species, with nearly 24-36% of all genes in each B. braunii race being specific to each race. This analysis suggests the three races are distinct species based on sufficient differences in their respective genomes. We propose reclassification of the three chemical races to the following species names: Botryococcus alkenealis (A race), Botryococcus braunii (B race), and Botryococcus lycopadienor (L race).

The Effect of Copy Number Hemiplasy on Gene Family Evolution

The evolution of gene families is complex, involving gene-level evolutionary events such as gene duplication, horizontal gene transfer, and gene loss, and other processes such as incomplete lineage sorting (ILS). Because of this, topological differences often exist between gene trees and species trees. A number of models have been recently developed to explain these discrepancies, the most realistic of which attempts to consider both gene-level events and ILS. When unified in a single model, the interaction between ILS and gene-level events can cause polymorphism in gene copy number, which we refer to as copy number hemiplasy (CNH). In this paper, we extend the Wright-Fisher process to include duplications and losses over several species, and show that the probability of CNH for this process can be significant. We study how well two unified models-multilocus multispecies coalescent (MLMSC), which models CNH, and duplication, loss, and coalescence (DLCoal), which does not-approximate the Wright-Fisher process with duplication and loss. We then study the effect of CNH on gene family evolution by comparing MLMSC and DLCoal. We generate comparable gene trees under both models, showing significant differences in various summary statistics; most importantly, CNH reduces the number of gene copies greatly. If this is not taken into account, the traditional method of estimating duplication rates (by counting the number of gene copies) becomes inaccurate. The simulated gene trees are also used for species tree inference with the summary methods ASTRAL and ASTRAL-Pro, demonstrating that their accuracy, based on CNH-unaware simulations calibrated on real data, may have been overestimated.

Inferences on the evolution of the ascorbic acid synthesis pathway in insects using Phylogenetic Tree Collapser (PTC), a tool for the automated collapsing of phylogenetic trees using taxonomic information

When inferring the evolution of a gene/gene family, it is advisable to use all available coding sequences (CDS) from as many species genomes as possible in order to infer and date all gene duplications and losses. Nowadays, this means using hundreds or even thousands of CDSs, which makes the inferred phylogenetic trees difficult to visualize and interpret. Therefore, it is useful to have an automated way of collapsing large phylogenetic trees according to a taxonomic term decided by the user (family, class, or order, for instance), in order to highlight the minimal set of sequences that should be used to recapitulate the full history of the gene/gene family being studied at that taxonomic level, that can be refined using additional software. Here we present the Phylogenetic Tree Collapser (PTC) program (https://github.com/pegi3s/phylogenetic-tree-collapser), a flexible tool for automated tree collapsing using taxonomic information, that can be easily used by researchers without a background in informatics, since it only requires the installation of Docker, Podman or Singularity. The utility of PTC is demonstrated by addressing the evolution of the ascorbic acid synthesis pathway in insects. A Docker image is available at Docker Hub (https://hub.docker.com/r/pegi3s/phylogenetic-tree-collapser) with PTC installed and ready-to-run.

The promise and pitfalls of synteny in phylogenomics

Reconstructing the tree of life remains a central goal in biology. Early methods, which relied on small numbers of morphological or genetic characters, often yielded conflicting evolutionary histories, undermining confidence in the results. Investigations based on phylogenomics, which use hundreds to thousands of loci for phylogenetic inquiry, have provided a clearer picture of life's history, but certain branches remain problematic. To resolve difficult nodes on the tree of life, 2 recent studies tested the utility of synteny, the conserved collinearity of orthologous genetic loci in 2 or more organisms, for phylogenetics. Synteny exhibits compelling phylogenomic potential while also raising new challenges. This Essay identifies and discusses specific opportunities and challenges that bear on the value of synteny data and other rare genomic changes for phylogenomic studies. Synteny-based analyses of highly contiguous genome assemblies mark a new chapter in the phylogenomic era and the quest to reconstruct the tree of life.

Adaptation in Unstable Environments and Global Gene Losses: Small but Stable Gene Networks by the May-Wigner Theory

Although gene loss is common in evolution, it remains unclear whether it is an adaptive process. In a survey of seven major mangrove clades that are woody plants in the intertidal zones of daily environmental perturbations, we noticed that they generally evolved reduced gene numbers. We then focused on the largest clade of Rhizophoreae and observed the continual gene set reduction in each of the eight species. A great majority of gene losses are concentrated on environmental interaction processes, presumably to cope with the constant fluctuations in the tidal environments. Genes of the general processes for woody plants are largely retained. In particular, fewer gene losses are found in physiological traits such as viviparous seeds, high salinity, and high tannin content. Given the broad and continual genome reductions, we propose the May-Wigner theory (MWT) of system stability as a possible mechanism. In MWT, the most effective solution for buffering continual perturbations is to reduce the size of the system (or to weaken the total genic interactions). Mangroves are unique as immovable inhabitants of the compound environments in the land-sea interface, where environmental gradients (such as salinity) fluctuate constantly, often drastically. Extending MWT to gene regulatory network (GRN), computer simulations and transcriptome analyses support the stabilizing effects of smaller gene sets in mangroves vis-à-vis inland plants. In summary, we show the adaptive significance of gene losses in mangrove plants, including the specific role of promoting phenotype innovation and a general role in stabilizing GRN in unstable environments as predicted by MWT.

Combining extracellular matrix proteome and phosphoproteome of chickpea and meta-analysis reveal novel proteoforms and evolutionary significance of clade-specific wall-associated events in plant

Extracellular matrix (ECM) plays central roles in cell architecture, innate defense and cell wall integrity (CWI) signaling. During transition to multicellularity, modular domain structures of ECM proteins and proteoforms have evolved due to continuous adaptation across taxonomic clades under different ecological niche. Although this incredible diversity has to some extent been investigated at protein level, extracellular phosphorylation events and molecular evolution of ECM proteoform families remains unexplored. We developed matrisome proteoform atlas in a grain legume, chickpea and performed meta-analyses of 74 plant matrisomes. MS/MS analysis identified 1,424 proteins and 315 phosphoproteins involved in diverse functions. Cross-species ECM protein network identified proteoforms associated with CWI maintenance system. Phylogenetic characterization of eighteen matrix protein families highlighted the role of taxon-specific paralogs and orthologs. Novel information was acquired on gene expansion and loss, co-divergence, sub functionalization and neofunctionalization during evolution. Modular networks of matrix protein families and hub proteins showed higher diversity across taxonomic clades than among organs. Furthermore, protein families differ in nonsynonymous to synonymous substitution rates. Our study pointed towards the matrix proteoform functionality, sequence divergence variation, interactions between wall remodelers and molecular evolution using a phylogenetic framework. This is the first report on comprehensive matrisome proteoform network illustrating presence of CWI signaling proteins in land plants.

Genomic evidence for rediploidization and adaptive evolution following the whole-genome triplication

Whole-genome duplication (WGD), or polyploidy, events are widespread and significant in the evolutionary history of angiosperms. However, empirical evidence for rediploidization, the major process where polyploids give rise to diploid descendants, is still lacking at the genomic level. Here we present chromosome-scale genomes of the mangrove tree Sonneratia alba and the related inland plant Lagerstroemia speciosa. Their common ancestor has experienced a whole-genome triplication (WGT) approximately 64 million years ago coinciding with a period of dramatic global climate change. Sonneratia, adapting mangrove habitats, experienced extensive chromosome rearrangements post-WGT. We observe the WGT retentions display sequence and expression divergence, suggesting potential neo- and sub-functionalization. Strong selection acting on three-copy retentions indicates adaptive value in response to new environments. To elucidate the role of ploidy changes in genome evolution, we improve a model of the polyploidization-rediploidization process based on genomic evidence, contributing to the understanding of adaptive evolution during climate change.

A chromosome-level genome assembly of Agave hybrid NO.11648 provides insights into the CAM photosynthesis

The subfamily Agavoideae comprises crassulacean acid metabolism (CAM), C3, and C4 plants with a young age of speciation and slower mutation accumulation, making it a model crop for studying CAM evolution. However, the genetic mechanism underlying CAM evolution remains unclear because of lacking genomic information. This study assembled the genome of Agave hybrid NO.11648, a constitutive CAM plant belonging to subfamily Agavoideae, at the chromosome level using data generated from high-throughput chromosome conformation capture, Nanopore, and Illumina techniques, resulting in 30 pseudo-chromosomes with a size of 4.87 Gb and scaffold N50 of 186.42 Mb. The genome annotation revealed 58 841 protein-coding genes and 76.91% repetitive sequences, with the dominant repetitive sequences being the I-type repeats (Copia and Gypsy accounting for 18.34% and 13.5% of the genome, respectively). Our findings also provide support for a whole genome duplication event in the lineage leading to A. hybrid, which occurred after its divergence from subfamily Asparagoideae. Moreover, we identified a gene duplication event in the phosphoenolpyruvate carboxylase kinase (PEPCK) gene family and revealed that three PEPCK genes (PEPCK3, PEPCK5, and PEPCK12) were involved in the CAM pathway. More importantly, we identified transcription factors enriched in the circadian rhythm, MAPK signaling, and plant hormone signal pathway that regulate the PEPCK3 expression by analysing the transcriptome and using yeast one-hybrid assays. Our results shed light on CAM evolution and offer an essential resource for the molecular breeding program of Agave spp.

Genome-wide identification and characterization of NBLRR genes in finger millet (Eleusine coracana L.) and their expression in response to Magnaporthe grisea infection

BACKGROUND

The nucleotide binding site leucine rich repeat (NBLRR) genes significantly regulate defences against phytopathogens in plants. The genome-wide identification and analysis of NBLRR genes have been performed in several species. However, the detailed evolution, structure, expression of NBLRRs and functional response to Magnaporthe grisea are unknown in finger millet (Eleusine coracana (L.) Gaertn.).

RESULTS

The genome-wide scanning of the finger millet genome resulted in 116 NBLRR (EcNBLRRs1-116) encompassing 64 CC-NB-LRR, 47 NB-LRR and 5 CCR-NB-LRR types. The evolutionary studies among the NBLRRs of five Gramineae species, viz., purple false brome (Brachypodium distachyon (L.) P.Beauv.), finger millet (E. coracana), rice (Oryza sativa L.), sorghum (Sorghum bicolor L. (Moench)) and foxtail millet (Setaria italica (L.) P.Beauv.) showed the evolution of NBLRRs in the ancestral lineage of the target species and subsequent divergence through gene-loss events. The purifying selection (Ka/Ks < 1) shaped the expansions of NBLRRs paralogs in finger millet and orthologs among the target Gramineae species. The promoter sequence analysis showed various stress- and phytohormone-responsive cis-acting elements besides growth and development, indicating their potential role in disease defence and regulatory mechanisms. The expression analysis of 22 EcNBLRRs in the genotypes showing contrasting responses to Magnaporthe grisea infection revealed four and five EcNBLRRs in early and late infection stages, respectively. The six of these nine candidate EcNBLRRs proteins, viz., EcNBLRR21, EcNBLRR26, EcNBLRR30, EcNBLRR45, EcNBLRR55 and EcNBLRR76 showed CC, NB and LRR domains, whereas the EcNBLRR23, EcNBLRR32 and EcNBLRR83 showed NB and LRR somains.

CONCLUSION

The identification and expression analysis of EcNBLRRs showed the role of EcNBLRR genes in assigning blast resistance in finger millet. These results pave the foundation for in-depth and targeted functional analysis of EcNBLRRs through genome editing and transgenic approaches.

Phylogenetic reconciliation: making the most of genomes to understand microbial ecology and evolution

In recent years, phylogenetic reconciliation has emerged as a promising approach for studying microbial ecology and evolution. The core idea is to model how gene trees evolve along a species tree and to explain differences between them via evolutionary events including gene duplications, transfers, and losses. Here, we describe how phylogenetic reconciliation provides a natural framework for studying genome evolution and highlight recent applications including ancestral gene content inference, the rooting of species trees, and the insights into metabolic evolution and ecological transitions they yield. Reconciliation analyses have elucidated the evolution of diverse microbial lineages, from Chlamydiae to Asgard archaea, shedding light on ecological adaptation, host-microbe interactions, and symbiotic relationships. However, there are many opportunities for broader application of the approach in microbiology. Continuing improvements to make reconciliation models more realistic and scalable, and integration of ecological metadata such as habitat, pH, temperature, and oxygen use offer enormous potential for understanding the rich tapestry of microbial life.

The Theory of Gene Family Histories

Most genes are part of larger families of evolutionary-related genes. The history of gene families typically involves duplications and losses of genes as well as horizontal transfers into other organisms. The reconstruction of detailed gene family histories, i.e., the precise dating of evolutionary events relative to phylogenetic tree of the underlying species has remained a challenging topic despite their importance as a basis for detailed investigations into adaptation and functional evolution of individual members of the gene family. The identification of orthologs, moreover, is a particularly important subproblem of the more general setting considered here. In the last few years, an extensive body of mathematical results has appeared that tightly links orthology, a formal notion of best matches among genes, and horizontal gene transfer. The purpose of this chapter is to broadly outline some of the key mathematical insights and to discuss their implication for practical applications. In particular, we focus on tree-free methods, i.e., methods to infer orthology or horizontal gene transfer as well as gene trees, species trees, and reconciliations between them without using a priori knowledge of the underlying trees or statistical models for the inference of phylogenetic trees. Instead, the initial step aims to extract binary relations among genes.

Cophylogeny Reconstruction Allowing for Multiple Associations Through Approximate Bayesian Computation

Phylogenetic tree reconciliation is extensively employed for the examination of coevolution between host and symbiont species. An important concern is the requirement for dependable cost values when selecting event-based parsimonious reconciliation. Although certain approaches deduce event probabilities unique to each pair of host and symbiont trees, which can subsequently be converted into cost values, a significant limitation lies in their inability to model the invasion of diverse host species by the same symbiont species (termed as a spread event), which is believed to occur in symbiotic relationships. Invasions lead to the observation of multiple associations between symbionts and their hosts (indicating that a symbiont is no longer exclusive to a single host), which are incompatible with the existing methods of coevolution. Here, we present a method called AmoCoala (an enhanced version of the tool Coala) that provides a more realistic estimation of cophylogeny event probabilities for a given pair of host and symbiont trees, even in the presence of spread events. We expand the classical 4-event coevolutionary model to include 2 additional outcomes, vertical and horizontal spreads, that lead to multiple associations. In the initial step, we estimate the probabilities of spread events using heuristic frequencies. Subsequently, in the second step, we employ an approximate Bayesian computation approach to infer the probabilities of the remaining 4 classical events (cospeciation, duplication, host switch, and loss) based on these values. By incorporating spread events, our reconciliation model enables a more accurate consideration of multiple associations. This improvement enhances the precision of estimated cost sets, paving the way to a more reliable reconciliation of host and symbiont trees. To validate our method, we conducted experiments on synthetic datasets and demonstrated its efficacy using real-world examples. Our results showcase that AmoCoala produces biologically plausible reconciliation scenarios, further emphasizing its effectiveness.

Gene Transfer-Based Phylogenetics: Analytical Expressions and Additivity via Birth-Death Theory

The genomic era has opened up vast opportunities in molecular systematics, one of which is deciphering the evolutionary history in fine detail. Under this mass of data, analyzing the point mutations of standard markers is often too crude and slow for fine-scale phylogenetics. Nevertheless, genome dynamics (GD) events provide alternative, often richer information. The synteny index (SI) between a pair of genomes combines gene order and gene content information, allowing the comparison of genomes of unequal gene content, together with order considerations of their common genes. Recently, genome dynamics has been modeled as a continuous-time Markov process, and gene distance in the genome as a birth-death-immigration process. Nevertheless, due to complexities arising in this setting, no precise and provably consistent estimators could be derived, resulting in heuristic solutions. Here, we extend this modeling approach by using techniques from birth-death theory to derive explicit expressions of the system's probabilistic dynamics in the form of rational functions of the model parameters. This, in turn, allows us to infer analytically accurate distances between organisms based on their SI. Subsequently, we establish additivity of this estimated evolutionary distance (a desirable property yielding phylogenetic consistency). Applying the new measure in simulation studies shows that it provides accurate results in realistic settings and even under model extensions such as gene gain/loss or over a tree structure. In the real-data realm, we applied the new formulation to unique data structure that we constructed-the ordered orthology DB-based on a new version of the EggNOG database, to construct a tree with more than 4.5K taxa. To the best of our knowledge, this is the largest gene-order-based tree constructed and it overcomes shortcomings found in previous approaches. Constructing a GD-based tree allows to confirm and contrast findings based on other phylogenetic approaches, as we show.

Origin, evolution, and diversification of the wall-associated kinase gene family in plants

KEY MESSAGE

The study of the origin, evolution, and diversification of the wall-associated kinase gene family in plants facilitates their functional investigations in the future. Wall-associated kinases (WAKs) make up one subfamily of receptor-like kinases (RLKs), and function directly in plant cell elongation and responses to biotic and abiotic stresses. The biological functions of WAKs have been extensively characterized in angiosperms; however, the origin and evolutionary history of the WAK family in green plants remain unclear. Here, we performed a comprehensive analysis of the WAK family to reveal its origin, evolution, and diversification in green plants. In total, 1061 WAK genes were identified in 37 species from unicellular algae to multicellular plants, and the results showed that WAK genes probably originated before bryophyte differentiation and were widely distributed in land plants, especially angiosperms. The phylogeny indicated that the land plant WAKs gave rise to five clades and underwent lineage-specific expansion after species differentiation. Cis-acting elements and expression patterns analyses of WAK genes in Arabidopsis and rice demonstrated the functional diversity of WAK genes in these two species. Many gene gains and losses have occurred in angiosperms, leading to an increase in the number of gene copies. The evolutionary trajectory of the WAK family during polyploidization was uncovered using Gossypium species. Our results provide insights into the evolution of WAK genes in green plants, facilitating their functional investigations in the future.

Incongruence in the phylogenomics era

Genome-scale data and the development of novel statistical phylogenetic approaches have greatly aided the reconstruction of a broad sketch of the tree of life and resolved many of its branches. However, incongruence - the inference of conflicting evolutionary histories - remains pervasive in phylogenomic data, hampering our ability to reconstruct and interpret the tree of life. Biological factors, such as incomplete lineage sorting, horizontal gene transfer, hybridization, introgression, recombination and convergent molecular evolution, can lead to gene phylogenies that differ from the species tree. In addition, analytical factors, including stochastic, systematic and treatment errors, can drive incongruence. Here, we review these factors, discuss methodological advances to identify and handle incongruence, and highlight avenues for future research.

Relative timing information and orthology in evolutionary scenarios

BACKGROUND

Evolutionary scenarios describing the evolution of a family of genes within a collection of species comprise the mapping of the vertices of a gene tree T to vertices and edges of a species tree S. The relative timing of the last common ancestors of two extant genes (leaves of T) and the last common ancestors of the two species (leaves of S) in which they reside is indicative of horizontal gene transfers (HGT) and ancient duplications. Orthologous gene pairs, on the other hand, require that their last common ancestors coincides with a corresponding speciation event. The relative timing information of gene and species divergences is captured by three colored graphs that have the extant genes as vertices and the species in which the genes are found as vertex colors: the equal-divergence-time (EDT) graph, the later-divergence-time (LDT) graph and the prior-divergence-time (PDT) graph, which together form an edge partition of the complete graph.

RESULTS

Here we give a complete characterization in terms of informative and forbidden triples that can be read off the three graphs and provide a polynomial time algorithm for constructing an evolutionary scenario that explains the graphs, provided such a scenario exists. While both LDT and PDT graphs are cographs, this is not true for the EDT graph in general. We show that every EDT graph is perfect. While the information about LDT and PDT graphs is necessary to recognize EDT graphs in polynomial-time for general scenarios, this extra information can be dropped in the HGT-free case. However, recognition of EDT graphs without knowledge of putative LDT and PDT graphs is NP-complete for general scenarios. In contrast, PDT graphs can be recognized in polynomial-time. We finally connect the EDT graph to the alternative definitions of orthology that have been proposed for scenarios with horizontal gene transfer. With one exception, the corresponding graphs are shown to be colored cographs.

What can be lost? Genomic perspective on the lipid metabolism of Mucoromycota

Mucoromycota is a phylum of early diverging fungal (EDF) lineages, of mostly plant-associated terrestrial fungi. Some strains have been selected as promising biotechnological organisms due to their ability to produce polyunsaturated fatty acids and efficient conversion of nutrients into lipids. Others get their lipids from the host plant and are unable to produce even the essential ones on their own. Following the advancement in EDF genome sequencing, we carried out a systematic survey of lipid metabolism protein families across different EDF lineages. This enabled us to explore the genomic basis of the previously documented ability to produce several types of lipids within the fungal tree of life. The core lipid metabolism genes showed no significant diversity in distribution, however specialized lipid metabolic pathways differed in this regard among different fungal lineages. In total 165 out of 202 genes involved in lipid metabolism were present in all tested fungal lineages, while remaining 37 genes were found to be absent in some of fungal lineages. Duplications were observed for 69 genes. For the first time we demonstrate that ergosterol is not being produced by several independent groups of plant-associated fungi due to the losses of different ERG genes. Instead, they possess an ancestral pathway leading to the synthesis of cholesterol, which is absent in other fungal lineages. The lack of diacylglycerol kinase in both Mortierellomycotina and Blastocladiomycota opens the question on sterol equilibrium regulation in these organisms. Early diverging fungi retained most of beta oxidation components common with animals including Nudt7, Nudt12 and Nudt19 pointing at peroxisome divergence in Dikarya. Finally, Glomeromycotina and Mortierellomycotina representatives have a similar set of desaturases and elongases related to the synthesis of complex, polyunsaturated fatty acids pointing at an ancient expansion of fatty acid metabolism currently being explored by biotechnological studies.

Molecular Evolution of Aryl Hydrocarbon Receptor Signaling Pathway Genes

The Aryl hydrocarbon receptor is an ancient transcriptional factor originally discovered as a sensor of dioxin. In addition to its function as a receptor of environmental toxicants, it plays an important role in development. Although a significant amount of research has been carried out to understand the AHR signal transduction pathway and its involvement in species' susceptibility to environmental toxicants, none of them to date has comprehensively studied its evolutionary origins. Studying the evolutionary origins of molecules can inform ancestral relationships of genes. The vertebrate genome has been shaped by two rounds of whole-genome duplications (WGD) at the base of vertebrate evolution approximately 600 million years ago, followed by lineage-specific gene losses, which often complicate the assignment of orthology. It is crucial to understand the evolutionary origins of this transcription factor and its partners, to distinguish orthologs from ancient non-orthologous homologs. In this study, we have investigated the evolutionary origins of proteins involved in the AHR pathway. Our results provide evidence of gene loss and duplications, crucial for understanding the functional connectivity of humans and model species. Multiple studies have shown that 2R-ohnologs (genes and proteins that have survived from the 2R-WGD) are enriched in signaling components relevant to developmental disorders and cancer. Our findings provide a link between the AHR pathway's evolutionary trajectory and its potential mechanistic involvement in pathogenesis.

Genome-Wide Comparison and Functional Characterization of HMGR Gene Family Associated with Shikonin Biosynthesis in Lithospermum erythrorhizon

3-hydroxy-3-methylglutaryl-CoA reductase (HMGR), as the rate-limiting enzyme in the mevalonate pathway, is essential for the biosynthesis of shikonin in Lithospermum erythrorhizon. However, in the absence of sufficient data, the principles of a genome-wide in-depth evolutionary exploration of HMGR family members in plants, as well as key members related to shikonin biosynthesis, remain unidentified. In this study, 124 HMGRs were identified and characterized from 36 representative plants, including L. erythrorhizon. Vascular plants were found to have more HMGR family genes than nonvascular plants. The phylogenetic tree revealed that during lineage and species diversification, the HMGRs evolved independently and intronless LerHMGRs emerged from multi-intron HMGR in land plants. Among them, Pinus tabuliformis and L. erythrorhizon had the most HMGR gene duplications, with 11 LerHMGRs most likely expanded through WGD/segmental and tandem duplications. In seedling roots and M9 cultured cells/hairy roots, where shikonin biosynthesis occurs, LerHMGR1 and LerHMGR2 were expressed significantly more than other genes. The enzymatic activities of LerHMGR1 and LerHMGR2 further supported their roles in catalyzing the conversion of HMG-CoA to mevalonate. Our findings provide insight into the molecular evolutionary properties and function of the HMGR family in plants and a basis for the genetic improvement of efficiently produced secondary metabolites in L. erythrorhizon.

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

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

Extreme genome diversity and cryptic speciation in a harmful algal-bloom-forming eukaryote

Harmful algal blooms of the toxic haptophyte Prymnesium parvum are a recurrent problem in many inland and estuarine waters around the world. Strains of P. parvum vary in the toxins they produce and in other physiological traits associated with harmful algal blooms, but the genetic basis for this variation is unknown. To investigate genome diversity in this morphospecies, we generated genome assemblies for 15 phylogenetically and geographically diverse strains of P. parvum, including Hi-C guided, near-chromosome-level assemblies for two strains. Comparative analysis revealed considerable DNA content variation between strains, ranging from 115 to 845 Mbp. Strains included haploids, diploids, and polyploids, but not all differences in DNA content were due to variation in genome copy number. Haploid genome size between strains of different chemotypes differed by as much as 243 Mbp. Syntenic and phylogenetic analyses indicate that UTEX 2797, a common laboratory strain from Texas, is a hybrid that retains two phylogenetically distinct haplotypes. Investigation of gene families variably present across the strains identified several functional categories associated with metabolic and genome size variation in P. parvum, including genes for the biosynthesis of toxic metabolites and proliferation of transposable elements. Together, our results indicate that P. parvum comprises multiple cryptic species. These genomes provide a robust phylogenetic and genomic framework for investigations into the eco-physiological consequences of the intra- and inter-specific genetic variation present in P. parvum and demonstrate the need for similar resources for other harmful algal-bloom-forming morphospecies.

The Gynandropsis gynandra genome provides insights into whole-genome duplications and the evolution of C4 photosynthesis in Cleomaceae

Gynandropsis gynandra (Cleomaceae) is a cosmopolitan leafy vegetable and medicinal plant, which has also been used as a model to study C4 photosynthesis due to its evolutionary proximity to C3 Arabidopsis (Arabidopsis thaliana). Here, we present the genome sequence of G. gynandra, anchored onto 17 main pseudomolecules with a total length of 740 Mb, an N50 of 42 Mb and 30,933 well-supported gene models. The G. gynandra genome and previously released genomes of C3 relatives in the Cleomaceae and Brassicaceae make an excellent model for studying the role of genome evolution in the transition from C3 to C4 photosynthesis. Our analyses revealed that G. gynandra and its C3 relative Tarenaya hassleriana shared a whole-genome duplication event (Gg-α), then an addition of a third genome (Th-α, +1×) took place in T. hassleriana but not in G. gynandra. Analysis of syntenic copy number of C4 photosynthesis-related gene families indicates that G. gynandra generally retained more duplicated copies of these genes than C3T. hassleriana, and also that the G. gynandra C4 genes might have been under positive selection pressure. Both whole-genome and single-gene duplication were found to contribute to the expansion of the aforementioned gene families in G. gynandra. Collectively, this study enhances our understanding of the polyploidy history, gene duplication and retention, as well as their impact on the evolution of C4 photosynthesis in Cleomaceae.

Exon shuffling potentiates a diverse repertoire of brown algal NB-ARC-TPR candidate immune receptor proteins via alternative splicing

Like other organisms, brown algae are subject to diseases caused by bacteria, fungi, and viruses. Brown algal immunity mechanisms are not well characterized; however, there is evidence suggesting that pathogen receptors exist in brown algae. One key protein family likely associated with brown algal innate immunity possesses an NB-ARC domain analogous to innate immune proteins in plants and animals. In this study, we conducted an extensive survey of NB-ARC genes in brown algae and obtained insights into the domain organization and evolutionary history of the encoded proteins. Our data show that brown algae possess an ancient NB-ARC-tetratricopeptide repeat (NB-TPR) domain architecture. We identified an N-terminal effector domain, the four-helix bundle, which was not previously found associated with NB-ARC domains. The phylogenetic tree including NB-ARC domains from all kingdoms of life suggests the three clades of brown algal NB-TPRs are likely monophyletic, whereas their TPRs seem to have distinct origins. One group of TPRs exhibit intense exon shuffling, with various alternative splicing and diversifying selection acting on them, suggesting exon shuffling is an important mechanism for evolving ligand-binding specificities. The reconciliation of gene duplication and loss events of the NB-ARC genes reveals that more independent gene gains than losses have occurred during brown algal evolution, and that tandem duplication has played a major role in the expansion of NB-ARC genes. Our results substantially enhance our understanding of the evolutionary history and exon shuffling mechanisms of the candidate innate immune repertoire of brown algae.

Evolutionary history of two evergreen Rhododendron species as revealed by chromosome-level genome assembly

Background

The genus Rhododendron (Ericaceae), a species-rich and widely distributed genus of woody plants, is distinguished for the beautiful and diverse flowers. Rhododendron delavayi Franch. and Rhododendron irroratum Franch., are highly attractive species widely distributed in south-west China and abundant new varieties have been selected from their genetic resources.

Methods

We constructed chromosome-scale genome assemblies for Rhododendron delavayi and Rhododendron irroratum. Phylogenetic and whole-genome duplication analyses were performed to elucidate the evolutionary history of Rhododendron. Further, different types of gene duplications were identified and their contributions to gene family expansion were investigated. Finally, comprehensive characterization and evolutionary analysis of R2R3-MYB and NBS-encoding genes were conducted to explore their evolutionary patterns.

Results

The phylogenetic analysis classified Rhododendron species into two sister clades, 'rhododendrons' and 'azaleas'. Whole-genome duplication (WGD) analysis unveiled only one WGD event that occurred in Rhododendron after the ancestral γ triplication. Gene duplication and gene family expansion analyses suggested that the younger tandem and proximal duplications contributed greatly to the expansion of gene families involved in secondary metabolite biosynthesis and stress response. The candidate R2R3-MYB genes likely regulating anthocyanin biosynthesis and stress tolerance in Rhododendron will facilitate the breeding for ornamental use. NBS-encoding genes had undergone significant expansion and experienced species-specific gain and loss events in Rhododendron plants.

Conclusions

The reference genomes presented here will provide important genetic resources for molecular breeding and genetic improvement of plants in this economically important Rhododendron genus.

virDTL: Viral Recombination Analysis Through Phylogenetic Reconciliation and Its Application to Sarbecoviruses and SARS-CoV-2

An accurate understanding of the evolutionary history of rapidly-evolving viruses like SARS-CoV-2, responsible for the COVID-19 pandemic, is crucial to tracking and preventing the spread of emerging pathogens. However, viruses undergo frequent recombination, which makes it difficult to trace their evolutionary history using traditional phylogenetic methods. In this study, we present a phylogenetic workflow, virDTL, for analyzing viral evolution in the presence of recombination. Our approach leverages reconciliation methods developed for inferring horizontal gene transfer in prokaryotes and, compared to existing tools, is uniquely able to identify ancestral recombinations while accounting for several sources of inference uncertainty, including in the construction of a strain tree, estimation and rooting of gene family trees, and reconciliation itself. We apply this workflow to the Sarbecovirus subgenus and demonstrate how a principled analysis of predicted recombination gives insight into the evolution of SARS-CoV-2. In addition to providing confirming evidence for the horseshoe bat as its zoonotic origin, we identify several ancestral recombination events that merit further study.

Comparative-genomic analysis reveals dynamic NLR gene loss and gain across Apiaceae species

Introduction: Nucleotide-binding leucine-rich repeat (NLR) genes play a crucial role in green plants' responding to various pathogens. Genome-scale evolutionary studies of NLR genes are important for discovering and applying functional NLR genes. However, little is known about the evolution of NLR genes in the Apiaceae family including agricultural and medical plants. Methods: In this study, comparative genomic analysis was performed in four Apiaceae species to trace the dynamic evolutionary patterns of NLR genes during speciation in this family. Results: The results revealed different number of NLR genes in these four Apiaceae species, namely, Angelica sinensis (95), Coriandrum sativum (183), Apium graveolens (153) and Daucus carota (149). Phylogenetic analysis demonstrated that NLR genes in these four species were derived from 183 ancestral NLR lineages and experienced different levels of gene-loss and gain events. The contraction pattern of the ancestral NLR lineages was discovered during the evolution of D. carota, whereas a different pattern of contraction after first expansion of NLR genes was observed for A. sinensis, C. sativum and A. graveolens. Discussion: Taken together, rapid and dynamic gene content variation has shaped evolutionary history of NLR genes in Apiaceae species.

A Novel Trichothecene Toxin Phenotype Associated with Horizontal Gene Transfer and a Change in Gene Function in Fusarium

Fusarium trichothecenes are among the mycotoxins of most concern to food and feed safety. Production of these mycotoxins and presence of the trichothecene biosynthetic gene (TRI) cluster have been confirmed in only two multispecies lineages of Fusarium: the Fusarium incarnatum-equiseti (Incarnatum) and F. sambucinum (Sambucinum) species complexes. Here, we identified and characterized a TRI cluster in a species that has not been formally described and is represented by Fusarium sp. NRRL 66739. This fungus is reported to be a member of a third Fusarium lineage: the F. buharicum species complex. Cultures of NRRL 66739 accumulated only two trichothecenes, 7-hydroxyisotrichodermin and 7-hydroxyisotrichodermol. Although these are not novel trichothecenes, the production profile of NRRL 66739 is novel, because in previous reports 7-hydroxyisotrichodermin and 7-hydroxyisotrichodermol were components of mixtures of 6-8 trichothecenes produced by several Fusarium species in Sambucinum. Heterologous expression analysis indicated that the TRI13 gene in NRRL 66739 confers trichothecene 7-hydroxylation. This contrasts the trichothecene 4-hydroxylation function of TRI13 in other Fusarium species. Phylogenetic analyses suggest that NRRL 66739 acquired the TRI cluster via horizontal gene transfer from a close relative of Incarnatum and Sambucinum. These findings provide insights into evolutionary processes that have shaped the distribution of trichothecene production among Fusarium species and the structural diversity of the toxins.

Karyotype evolution of the Asterids insights from the first genome sequences of the family Cornaceae

Cornaceae is a core representative family in Cornales, the earliest branching lineage in the Asterids on the life tree of angiosperms. This family includes the only genus Cornus, a group of ~55 species. These species occur widely in Northern Hemisphere and have been used as resources for horticultural ornaments, medicinal and industrial manufacturing. However, no any genome sequences are available for this family. Here, we reported a chromosome-level genome for Cornus controversa. This was generated using high-fidelity plus Hi-C sequencing, and totally ~771.80 Mb assembled sequences and 39,886 protein-coding genes were obtained. We provided evidence for a whole-genome duplication event (WGD) unique to C. controversa. The evolutionary features of this genome indicated that the expanded and unique genes might have contributed to response to stress, stimulus and defense. By using chromosome-level syntenic blocks shared between eight living genomes, we found high degrees of genomic diversification from the ancestral core-eudicot genome to the present-day genomes, suggesting an important role of WGD in genomic plasticity that leads to speciation and diversification. These results provide foundational insights on the evolutionary history of Cornaceae, as well as on the Asterids diversification.

Tree2GD: a phylogenomic method to detect large-scale gene duplication events

MOTIVATION

Whole-genome duplication events have long been discovered throughout the evolution of eukaryotes, contributing to genome complexity and biodiversity and leaving traces in the descending organisms. Therefore, an accurate and rapid phylogenomic method is needed to identify the retained duplicated genes on various lineages across the target taxonomy.

RESULTS

Here, we present Tree2GD, an integrated method to identify large-scale gene duplication events by automatically perform multiple procedures, including sequence alignment, recognition of homolog, gene tree/species tree reconciliation, Ks distribution of gene duplicates and synteny analyses. Application of Tree2GD on 2 datasets, 12 metazoan genomes and 68 angiosperms, successfully identifies all reported whole-genome duplication events exhibited by these species, showing effectiveness and efficiency of Tree2GD on phylogenomic analyses of large-scale gene duplications.

AVAILABILITY AND IMPLEMENTATION

Tree2GD is written in Python and C++ and is available at https://github.com/Dee-chen/Tree2gd.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Inferred ancestry of scytonemin biosynthesis proteins in cyanobacteria indicates a response to Paleoproterozoic oxygenation

Protection from radiation damage is an important adaptation for phototrophic microbes. Living in surface, shallow water, and peritidal environments, cyanobacteria are especially exposed to long-wavelength ultraviolet (UVA) radiation. Several groups of cyanobacteria within these environments are protected from UVA damage by the production of the pigment scytonemin. Paleontological evidence of cyanobacteria in UVA-exposed environments from the Proterozoic, and possibly as early as the Archaean, suggests a long evolutionary history of radiation protection within this group. We show that phylogenetic analyses of enzymes in the scytonemin biosynthesis pathway support this hypothesis and reveal a deep history of vertical inheritance of this pathway within extant cyanobacterial diversity. Referencing this phylogeny to cyanobacterial molecular clocks suggests that scytonemin production likely appeared during the early Proterozoic, soon after the Great Oxygenation Event. This timing is consistent with an adaptive scenario for the evolution of scytonemin production, wherein the threat of UVA-generated reactive oxygen species becomes significantly greater once molecular oxygen is more pervasive across photosynthetic environments.

Comprehensive Identification and Analyses of the GRF Gene Family in the Whole-Genome of Four Juglandaceae Species

The GRF gene family plays an important role in plant growth and development as regulators involved in plant hormone signaling and metabolism. However, the Juglandaceae GRF gene family remains to be studied. Here, we identified 15, 15, 19, and 20 GRF genes in J. regia, C. illinoinensis, J. sigillata, and J. mandshurica, respectively. The phylogeny shows that the Juglandaceae family GRF is divided into two subfamilies, the ε-group and the non-ε-group, and that selection pressure analysis did not detect amino acid loci subject to positive selection pressure. In addition, we found that the duplications of the Juglandaceae family GRF genes were all segmental duplication events, and a total of 79 orthologous gene pairs and one paralogous homologous gene pair were identified in four Juglandaceae families. The Ka/KS ratios between these homologous gene pairs were further analyzed, and the Ka/KS values were all less than 1, indicating that purifying selection plays an important role in the evolution of the Juglandaceae family GRF genes. The codon bias of genes in the GRF family of Juglandaceae species is weak, and is affected by both natural selection pressure and base mutation, and translation selection plays a dominant role in the mutation pressure in codon usage. Finally, expression analysis showed that GRF genes play important roles in pecan embryo development and walnut male and female flower bud development, but with different expression patterns. In conclusion, this study will serve as a rich genetic resource for exploring the molecular mechanisms of flower bud differentiation and embryo development in Juglandaceae. In addition, this is the first study to report the GRF gene family in the Juglandaceae family; therefore, our study will provide guidance for future comparative and functional genomic studies of the GRF gene family in the Juglandaceae specie.

Microbial communities of stratified aquatic ecosystems of Kandalaksha Bay (White Sea) shed light on the evolutionary history of green and brown morphotypes of Chlorobiota

Anoxygenic photoautotrophic metabolism of green sulfur bacteria of the family Chlorobiaceae played a significant role in establishing the Earth's biosphere. Two known major ecological forms of these phototrophs differ in their pigment composition and, therefore, in color: the green and brown forms. The latter form often occurs in low-light environments and is specialized to harvest blue light, which can penetrate to the greatest depth in the water column. In the present work, metagenomic sequencing was used to investigate the natural population of brown Chl. phaeovibrioides ZM in a marine stratified Zeleny Mys lagoon in the Kandalaksha Bay (the White Sea) to supplement the previously obtained genomes of brown Chlorobiaceae. The genomes of brown and green Chlorobiaceae were investigated using comparative genome analysis and phylogenetic and reconciliation analysis to reconstruct the evolution of these ecological forms. Our results support the suggestion that the last common ancestor of Chlorobiaceae belonged to the brown form, i.e. it was adapted to the conditions of low illumination. However, despite the vertical inheritance of these characteristics, among modern Chlorobiaceae populations, the genes responsible for synthesizing the pigments of the brown form are subject to active horizontal transfer.

A Polynomial-Time Algorithm for Minimizing the Deep Coalescence Cost for Level-1 Species Networks

Phylogenetic analyses commonly assume that the species history can be represented as a tree. However, in the presence of hybridization, the species history is more accurately captured as a network. Despite several advances in modeling phylogenetic networks, there is no known polynomial-time algorithm for parsimoniously reconciling gene trees with species networks while accounting for incomplete lineage sorting. To address this issue, we present a polynomial-time algorithm for the case of level-1 networks, in which no hybrid species is the direct ancestor of another hybrid species. This work enables more efficient reconciliation of gene trees with species networks, which in turn, enables more efficient reconstruction of species networks.

Dynamic genome evolution in a model fern

The large size and complexity of most fern genomes have hampered efforts to elucidate fundamental aspects of fern biology and land plant evolution through genome-enabled research. Here we present a chromosomal genome assembly and associated methylome, transcriptome and metabolome analyses for the model fern species Ceratopteris richardii. The assembly reveals a history of remarkably dynamic genome evolution including rapid changes in genome content and structure following the most recent whole-genome duplication approximately 60 million years ago. These changes include massive gene loss, rampant tandem duplications and multiple horizontal gene transfers from bacteria, contributing to the diversification of defence-related gene families. The insertion of transposable elements into introns has led to the large size of the Ceratopteris genome and to exceptionally long genes relative to other plants. Gene family analyses indicate that genes directing seed development were co-opted from those controlling the development of fern sporangia, providing insights into seed plant evolution. Our findings and annotated genome assembly extend the utility of Ceratopteris as a model for investigating and teaching plant biology.

Evolution of zygomycete secretomes and the origins of terrestrial fungal ecologies

Fungi survive in diverse ecological niches by secreting proteins and other molecules into the environment to acquire food and interact with various biotic and abiotic stressors. Fungal secretome content is, therefore, believed to be tightly linked to fungal ecologies. We sampled 132 genomes from the early-diverging terrestrial fungal lineage zygomycetes (Mucoromycota and Zoopagomycota) and characterized their secretome composition. Our analyses revealed that phylogeny played an important role in shaping the secretome composition of zygomycete fungi with trophic mode contributing a smaller amount. Reconstruction of the evolution of secreted digestive enzymes revealed lineage-specific expansions, indicating that Mucoromycota and Zoopagomycota followed different trajectories early in their evolutionary history. We identified the presence of multiple pathogenicity-related proteins in the lineages known as saprotrophs, suggesting that either the ecologies of these fungi are incompletely known, and/or that these pathogenicity-related proteins have important functions associated with saprotrophic ecologies, both of which invite further investigation.

Detection of interphylum transfers of the magnetosome gene cluster in magnetotactic bacteria

Magnetosome synthesis in magnetotactic bacteria (MTB) is regarded as a very ancient evolutionary process that dates back to deep-branching phyla. Magnetotactic bacteria belonging to one of such phyla, Nitrospirota, contain the classical genes for the magnetosome synthesis (e.g., mam, mms) and man genes, which were considered to be specific for this group. However, the recent discovery of man genes in MTB from the Thermodesulfobacteriota phylum has raised several questions about the inheritance of these genes in MTB. In this work, three new man genes containing MTB genomes affiliated with Nitrospirota and Thermodesulfobacteriota, were obtained. By applying reconciliation with these and the previously published MTB genomes, we demonstrate that the last common ancestor of all Nitrospirota was most likely not magnetotactic as assumed previously. Instead, our findings suggest that the genes for magnetosome synthesis were transmitted to the phylum Nitrospirota by horizontal gene transfer (HGT), which is the first case of the interphylum transfer of magnetosome genes detected to date. Furthermore, we provide evidence for the HGT of magnetosome genes from the Magnetobacteriaceae to the Dissulfurispiraceae family within Nitrospirota. Thus, our results imply a more significant role of HGT in the MTB evolution than deemed before and challenge the hypothesis of the ancient origin of magnetosome synthesis.

Spodoptera littoralis genome mining brings insights on the dynamic of expansion of gustatory receptors in polyphagous noctuidae

The bitter taste, triggered via gustatory receptors, serves as an important natural defense against the ingestion of poisonous foods in animals, and the increased host breadth is usually linked to an increase in the number of gustatory receptor genes. This has been especially observed in polyphagous insect species, such as noctuid species from the Spodoptera genus. However, the dynamic and physical mechanisms leading to these gene expansions and the evolutionary pressures behind them remain elusive. Among major drivers of genome dynamics are the transposable elements but, surprisingly, their potential role in insect gustatory receptor expansion has not been considered yet. In this work, we hypothesized that transposable elements and possibly positive selection would be involved in the highly dynamic evolution of gustatory receptor in Spodoptera spp. We first sequenced de novo the full 465 Mb genome of S. littoralis, and manually annotated the main chemosensory genes, including a large repertoire of 373 gustatory receptor genes (including 19 pseudogenes). We also improved the completeness of S. frugiperda and S. litura gustatory receptor gene repertoires. Then, we annotated transposable elements and revealed that a particular category of class I retrotransposons, the SINE transposons, was significantly enriched in the vicinity of gustatory receptor gene clusters, suggesting a transposon-mediated mechanism for the formation of these clusters. Selection pressure analyses indicated that positive selection within the gustatory receptor gene family is cryptic, only 7 receptors being identified as positively selected. Altogether, our data provide a new good quality Spodoptera genome, pinpoint interesting gustatory receptor candidates for further functional studies and bring valuable genomic information on the mechanisms of gustatory receptor expansions in polyphagous insect species.

The chromosome-level genome of a free-floating aquatic weed Pistia stratiotes provides insights into its rapid invasion

Pistia stratiotes (Araceae), commonly referred to as water lettuce, is one of the most notorious weeds that cause severe damage to the economy and natural ecosystems of infested areas. In order to explore the mechanism of its rapid invasion, here, we assembled a high-quality chromosome-level genome for P. stratiotes based on the Illumina sequencing, PacBio sequencing, and Hi-C scaffolding technology. The assembled genome is 311.87 Mb in size with a contig N50 of 1.08 Mb. The contigs were further anchored on 14 pseudochromosomes with a scaffold N50 of 21.21 Mb. A total of 20,356 protein-coding genes were predicted, of which 79.35% were functionally annotated here. Evolutionary analysis showed that P. stratiotes and Colocasia esculenta were clustered together as sister lineages that diverged approximately 61 Mya. The synteny analyses indicated that two whole-genome duplication (WGD) events occurred within a short period in P. stratiotes. Moreover, comparative genome analysis indicated that the expansion of gene families corresponding to disease resistance might contribute to rapid invasion in P. stratiotes. Also, we analyzed the disease-resistance gene family (NBS-LRR) involved in plant defense. A genome-wide search in P. stratiotes genome identified 85 NBS-LRR genes in this study. In conclusion, our present study provides some new insights into the evolution of the invasive aquatic plant P. stratiotes. Our reference genome will also provide valuable resources for future invasion genomics research programs.

The chromosome-scale genomes of Dipterocarpus turbinatus and Hopea hainanensis (Dipterocarpaceae) provide insights into fragrant oleoresin biosynthesis and hardwood formation

Dipterocarpaceae are typical tropical plants (dipterocarp forests) that are famous for their high economic value because of their production of fragrant oleoresins, top-quality timber and usage in traditional Chinese medicine. Currently, the lack of Dipterocarpaceae genomes has been a limiting factor to decipher the fragrant oleoresin biosynthesis and gain evolutionary insights into high-quality wood formation in Dipterocarpaceae. We generated chromosome-level genome assemblies for two representative Dipterocarpaceae species viz. Dipterocarpus turbinatus Gaertn. f. and Hopea hainanensis Merr. et Chun. Our whole-genome duplication (WGD) analysis revealed that Dipterocarpaceae underwent a shared WGD event, which showed significant impacts on increased copy numbers of genes related to the biosynthesis of terpene, BAHD acyltransferases, fatty acid and benzenoid/phenylpropanoid, which probably confer to the formation of their characteristic fragrant oleoresin. Additionally, compared with common soft wood plants, the expansion of gene families was also found to be associated with wood formation, such as in CESA (cellulose synthase), CSLE (cellulose synthase-like protein E), laccase and peroxidase in Dipterocarpaceae genomes, which might also contribute to the formation of harder, stronger and high-density timbers. Finally, an integrative analysis on a combination of genomic, transcriptomic and metabolic data from different tissues provided further insights into the molecular basis of fragrant oleoresins biosynthesis and high-quality wood formation of Dipterocarpaceae. Our study contributes the first two representative genomes for Dipterocarpaceae, which are valuable genetic resources for further researches on the fragrant oleoresins and superior-quality timber, genome-assisted breeding and improvement, and conservation biology of this family.

Comparative genomics reveal the convergent evolution of CYP82D and CYP706X members related to flavone biosynthesis in Lamiaceae and Asteraceae

Distant species producing the same secondary metabolites is an interesting and common phenomenon in nature. A classic example of this is scutellarein whose derivatives have been used clinically for more than 30 years. Scutellarein occurs in significant amounts in species of two different orders, Scutellaria baicalensis and Erigeron breviscapus, which diverged more than 100 million years ago. Here, according to the genome-wide selection and functional identification of 39 CYP450 genes from various angiosperms, we confirmed that only seven Scutellaria-specific CYP82D genes and one Erigeron CYP706X gene could perform the catalytic activity of flavone 6-hydroxylase (F6H), suggesting that the convergent evolution of scutellarein production in these two distant species was caused by two independently evolved CYP450 families. We also identified seven Scutellaria-specific CYP82D genes encoding flavone 8-hydroxylase (F8H). The evolutionary patterns of CYP82 and CYP706 families via kingdom-wide comparative genomics highlighted the evolutionary diversity of CYP82D and the specificity of CYP706X in angiosperms. Multi-collinearity and phylogenetic analysis of CYP82D in Scutellaria confirmed that the function of F6H evolved from F8H. Furthermore, the SbaiCYP82D1^A319D , EbreCYP706X^R130A , EbreCYP706X^F312D and EbreCYP706X^A318D mutants can significantly decrease the catalytic activity of F6H, revealing the contribution of crucial F6H amino acids to the scutellarein biosynthesis of distant species. This study provides important insights into the multi-origin evolution of the same secondary metabolite biosynthesis in the plant kingdom.

Large-scale analyses of heat shock transcription factors and database construction based on whole-genome genes in horticultural and representative plants

Heat shock transcription factor (Hsf) plays a critical role in regulating heat resistance. Here, 2950 Hsf family genes were identified from 111 horticultural and representative plants. More Hsf genes were detected in higher plants than lower plants. Based on all Hsf genes, we constructed a phylogenetic tree, which indicated that Hsf genes of each branch evolved independently after species differentiation. Furthermore, we uncovered the evolutionary trajectories of Hsf genes by motif analysis. There were only 6 motifs (M1 to M6) in lower plants, and then 4 novel motifs (M7-M10) appeared in higher plants. However, the motifs of some Hsf genes were lost in higher plant, indicating that Hsf genes have undergone sequence variation during the evolution. The number of Hsf gene loss was more than duplication after whole-genome duplication in higher plants. The heat response network was constructed using 24 Hsf genes, 2421 downstream, and 222 upstream genes of Arabidopsis. Further enrichment analysis revealed that Hsf genes and other transcription factors interacted with each other to response heat resistance. The global expression maps were illustrated for Hsf genes under various abiotic, biotic stresses, and several developmental stages in Arabidopsis. The syntenic and phylogenetic analyses were conducted using Hsf genes of Arabidopsis and Pan-genome of 18 Brassica rapa accessions. We also performed the expression pattern analysis of Hsf and six Hsp family genes using expression values from different tissues and heat treatments in B. rapa. The interaction network between Hsf and Hsp gene families was constructed in B. rapa, and several core genes were detected in the network. Finally, we constructed a Hsf database (http://hsfdb.bio2db.com) for researchers to retrieve Hsf gene family information. Therefore, our study will provide rich resources for the evolution and functional study of Hsf genes.

Efficiently sparse listing of classes of optimal cophylogeny reconciliations

Background

Cophylogeny reconciliation is a powerful method for analyzing host-parasite (or host-symbiont) co-evolution. It models co-evolution as an optimization problem where the set of all optimal solutions may represent different biological scenarios which thus need to be analyzed separately. Despite the significant research done in the area, few approaches have addressed the problem of helping the biologist deal with the often huge space of optimal solutions.

Results

In this paper, we propose a new approach to tackle this problem. We introduce three different criteria under which two solutions may be considered biologically equivalent, and then we propose polynomial-delay algorithms that enumerate only one representative per equivalence class (without listing all the solutions).

Conclusions

Our results are of both theoretical and practical importance. Indeed, as shown by the experiments, we are able to significantly reduce the space of optimal solutions while still maintaining important biological information about the whole space.