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

Proteomics, phylogenetics, and co-expression analyses indicate novel interactions in the plastid CLP chaperone-protease system

The chloroplast chaperone CLPC1 unfolds and delivers substrates to the stromal CLPPRT protease complex for degradation. We previously used an in vivo trapping approach to identify interactors with CLPC1 in Arabidopsis thaliana by expressing a STREPII-tagged copy of CLPC1 mutated in its Walker B domains (CLPC1-TRAP) followed by affinity purification and mass spectrometry. To create a larger pool of candidate substrates, adaptors, or regulators, we carried out a far more sensitive and comprehensive in vivo protein trapping analysis. We identified 59 highly enriched CLPC1 protein interactors, in particular proteins belonging to families of unknown functions (DUF760, DUF179, DUF3143, UVR-DUF151, HugZ/DUF2470), as well as the UVR domain proteins EXE1 and EXE2 implicated in singlet oxygen damage and signaling. Phylogenetic and functional domain analyses identified other members of these families that appear to localize (nearly) exclusively to plastids. In addition, several of these DUF proteins are of very low abundance as determined through the Arabidopsis PeptideAtlas http://www.peptideatlas.org/builds/arabidopsis/ showing that enrichment in the CLPC1-TRAP was extremely selective. Evolutionary rate covariation indicated that the HugZ/DUF2470 family coevolved with the plastid CLP machinery suggesting functional and/or physical interactions. Finally, mRNA-based coexpression networks showed that all 12 CLP protease subunits tightly coexpressed as a single cluster with deep connections to DUF760-3. Coexpression modules for other trapped proteins suggested specific functions in biological processes, e.g., UVR2 and UVR3 were associated with extraplastidic degradation, whereas DUF760-6 is likely involved in senescence. This study provides a strong foundation for discovery of substrate selection by the chloroplast CLP protease system.

Deciphering Microbial Gene Family Evolution Using Duplication-Transfer-Loss Reconciliation and RANGER-DTL

Phylogenetic reconciliation has emerged as a principled, highly effective technique for investigating the origin, spread, and evolutionary history of microbial gene families. Proper application of phylogenetic reconciliation requires a clear understanding of potential pitfalls and sources of error, and knowledge of the most effective reconciliation-based tools and protocols to use to maximize accuracy. In this book chapter, we provide a brief overview of Duplication-Transfer-Loss (DTL) reconciliation, the standard reconciliation model used to study microbial gene families and provide a step-by-step computational protocol to maximize the accuracy of DTL reconciliation and minimize false-positive evolutionary inferences.

Genome-wide Identification and Evolutionary Analysis of NBS-LRR Genes From Secale cereale

Secale cereale is an important crop in the Triticeae tribe of the Poaceae family, and it has unique agronomic characteristics and genome properties. It possesses resistance to many diseases and serves as an important resource for the breeding of other Triticeae crops. We performed a genome-wide study on S. cereale to identify the largest group of plant disease resistance genes (R genes), the nucleotide-binding site-leucine-rich repeat receptor (NBS-LRR) genes. In its genome, 582 NBS-LRR genes were identified, including one from the RNL subclass and 581 from the CNL subclass. The NBS-LRR gene number in the S. cereale genome is greater than that in barley and the diploid wheat genomes. S. cereale chromosome 4 contains the largest number of NBS-LRR genes among the seven chromosomes, which is different from the pattern in barley and the genomes B and D of wheat but similar to that in the genome A of wheat. Further synteny analysis suggests that more NBS-LRR genes on chromosome 4 have been inherited from a common ancestor by S. cereale and the wheat genome A than the wheat genomes B and D. Phylogenetic analysis revealed that at least 740 NBS-LRR lineages are present in the common ancestor of S. cereale, Hordeum vulgare and Triticum urartu. However, most of them have only been inherited by one or two species, with only 65 of them preserved in all three species. The S. cereale genome inherited 382 of these ancestral NBS-LRR lineages, but 120 of them have been lost in both H. vulgare and T. urartu. This study provides the full NBS-LRR profile of the S. cereale genome, which is a resource for S. cereale breeding and indicates that S. cereale can be an important material for the molecular breeding of other Triticeae crops.

Frequent Gene Duplication/Loss Shapes Distinct Evolutionary Patterns of NLR Genes in Arecaceae Species

Nucleotide-binding leucine-rich repeat (NLR) genes play a key role in plant immune responses and have co-evolved with pathogens since the origin of green plants. Comparative genomic studies on the evolution of NLR genes have been carried out in several angiosperm lineages. However, most of these lineages come from the dicot clade. In this study, comparative analysis was performed on NLR genes from five Arecaceae species to trace the dynamic evolutionary pattern of the gene family during species speciation in this monocot lineage. The results showed that NLR genes from the genomes of Elaeis guineensis (262), Phoenix dactylifera (85), Daemonorops jenkinsiana (536), Cocos nucifera (135) and Calamus simplicifolius (399) are highly variable. Frequent domain loss and alien domain integration have occurred to shape the NLR protein structures. Phylogenetic analysis revealed that NLR genes from the five genomes were derived from dozens of ancestral genes. D. jenkinsiana and E. guineensis genomes have experienced “consistent expansion” of the ancestral NLR lineages, whereas a pattern of “first expansion and then contraction” of NLR genes was observed for P. dactylifera, C. nucifera and C. simplicifolius. The results suggest that rapid and dynamic gene content and structure variation have shaped the NLR profiles of Arecaceae species.

Multiple Optimal Reconciliations Under the Duplication-Loss-Coalescence Model

Gene trees can differ from species trees due to a variety of biological phenomena, the most prevalent being gene duplication, horizontal gene transfer, gene loss, and coalescence. To explain topological incongruence between the two trees, researchers apply reconciliation methods, often relying on a maximum parsimony framework. However, while several studies have investigated the space of maximum parsimony reconciliations (MPRs) under the duplication-loss and duplication-transfer-loss models, the space of MPRs under the duplication-loss-coalescence (DLC) model remains poorly understood. To address this problem, we present new algorithms for computing the size of MPR space under the DLC model and sampling from this space uniformly at random. Our algorithms are efficient in practice, with runtime polynomial in the size of the species and gene tree when the number of genes that map to any given species is fixed, thus proving that the MPR problem is fixed-parameter tractable. We have applied our methods to a biological data set of 16 fungal species to provide the first key insights in the space of MPRs under the DLC model. Our results show that a plurality reconciliation, and underlying events, are likely to be representative of MPR space.

Reconciliation Reconsidered: In Search of a Most Representative Reconciliation in the Duplication-Transfer-Loss Model

Maximum parsimony reconciliation is a fundamental technique for studying the evolutionary histories of pairs of entities such as genes and species, parasites and hosts, and species and their biogeographical habitats. In these contexts, reconciliation is generally performed using the duplication-transfer-loss (DTL) model in a maximum parsimony framework. While efficient maximum parsimony reconciliation algorithms are known for the DTL model, the number of such reconciliations can grow exponentially with the sizes of the two phylogenetic trees. Choosing a maximum parsimony reconciliation arbitrarily may lead to conclusions that are not supported, and may even be contradicted, by other equally optimal reconciliations. This paper addresses the fundamental problem of how well a single reconciliation can represent the entire space of optimal reconciliations.

Odorant-binding proteins in canine anal sac glands indicate an evolutionarily conserved role in mammalian chemical communication

Background

Chemical communication is an important aspect of the behavioural ecology of a wide range of mammals. In dogs and other carnivores, anal sac glands are thought to convey information to conspecifics by secreting a pallet of small volatile molecules produced by symbiotic bacteria. Because these glands are unique to carnivores, it is unclear how their secretions relate to those of other placental mammals that make use of different tissues and secretions for chemical communication. Here we analyse the anal sac glands of domestic dogs to verify the secretion of proteins and infer their evolutionary relationship to those involved in the chemical communication of non-carnivoran mammals.

Results

Proteomic analysis of anal sac gland secretions of 17 dogs revealed the consistently abundant presence of three related proteins. Homology searches against online databases indicate that these proteins are evolutionary related to ‘odorant binding proteins’ (OBPs) found in a wide range of mammalian secretions and known to contribute to chemical communication. Screening of the dog’s genome sequence show that the newly discovered OBPs are encoded by a single cluster of three genes in the pseudoautosomal region of the X-chromosome. Comparative genomic screening indicates that the same locus is shared by a wide range of placental mammals and that it originated at least before the radiation of extant placental orders. Phylogenetic analyses suggest a dynamic evolution of gene duplication and loss, resulting in large gene clusters in some placental taxa and recurrent loss of this locus in others. The homology of OBPs in canid anal sac glands and those found in other mammalian secretions implies that these proteins maintained a function in chemical communication throughout mammalian evolutionary history by multiple shifts in expression between secretory tissues involved in signal release and nasal mucosa involved in signal reception.

Conclusions

Our study elucidates a poorly understood part of the biology of a species that lives in close association with humans. In addition, it shows that the protein repertoire underlying chemical communication in mammals is more evolutionarily stable than the variation of involved glands and tissues would suggest.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12862-021-01910-w.

De Novo Gene Birth, Horizontal Gene Transfer, and Gene Duplication as Sources of New Gene Families Associated with the Origin of Symbiosis in Amanita

By introducing novel capacities and functions, new genes and gene families may play a crucial role in ecological transitions. Mechanisms generating new gene families include de novo gene birth, horizontal gene transfer, and neofunctionalization following a duplication event. The ectomycorrhizal (ECM) symbiosis is a ubiquitous mutualism and the association has evolved repeatedly and independently many times among the fungi, but the evolutionary dynamics enabling its emergence remain elusive. We developed a phylogenetic workflow to first understand if gene families unique to ECM Amanita fungi and absent from closely related asymbiotic species are functionally relevant to the symbiosis, and then to systematically infer their origins. We identified 109 gene families unique to ECM Amanita species. Genes belonging to unique gene families are under strong purifying selection and are upregulated during symbiosis, compared with genes of conserved or orphan gene families. The origins of seven of the unique gene families are strongly supported as either de novo gene birth (two gene families), horizontal gene transfer (four), or gene duplication (one). An additional 34 families appear new because of their selective retention within symbiotic species. Among the 109 unique gene families, the most upregulated gene in symbiotic cultures encodes a 1-aminocyclopropane-1-carboxylate deaminase, an enzyme capable of downregulating the synthesis of the plant hormone ethylene, a common negative regulator of plant-microbial mutualisms.

Improved Duplication-Transfer-Loss Reconciliation with Extinct and Unsampled Lineages

Duplication-Transfer-Loss (DTL) reconciliation is a widely used computational technique for understanding gene family evolution and inferring horizontal gene transfer (transfer for short) in microbes. However, most existing models and implementations of DTL reconciliation cannot account for the effect of unsampled or extinct species lineages on the evolution of gene families, likely affecting their accuracy. Accounting for the presence and possible impact of any unsampled species lineages, including those that are extinct, is especially important for inferring and studying horizontal transfer since many genes in the species lineages represented in the reconciliation analysis are likely to have been acquired through horizontal transfer from unsampled lineages. While models of DTL reconciliation that account for transfer from unsampled lineages have already been proposed, they use a relatively simple framework for transfer from unsampled lineages and cannot explicitly infer the location on the species tree of each unsampled or extinct lineage associated with an identified transfer event. Furthermore, there does not yet exist any systematic studies to assess the impact of accounting for unsampled lineages on the accuracy of DTL reconciliation. In this work, we address these deficiencies by (i) introducing an extended DTL reconciliation model, called the DTLx reconciliation model, that accounts for unsampled and extinct species lineages in a new, more functional manner compared to existing models, (ii) showing that optimal reconciliations under the new DTLx reconciliation model can be computed just as efficiently as under the fastest DTL reconciliation model, (iii) providing an efficient algorithm for sampling optimal DTLx reconciliations uniformly at random, (iv) performing the first systematic simulation study to assess the impact of accounting for unsampled lineages on the accuracy of DTL reconciliation, and (v) comparing the accuracies of inferring transfers from unsampled lineages under our new model and the only other previously proposed parsimony-based model for this problem.

Maximum parsimony reconciliation in the DTLOR model

BACKGROUND

Analyses of microbial evolution often use reconciliation methods. However, the standard duplication-transfer-loss (DTL) model does not account for the fact that species trees are often not fully sampled and thus, from the perspective of reconciliation, a gene family may enter the species tree from the outside. Moreover, within the genome, genes are often rearranged, causing them to move to new syntenic regions.

RESULTS

We extend the DTL model to account for two events that commonly arise in the evolution of microbes: origin of a gene from outside the sampled species tree and rearrangement of gene syntenic regions. We describe an efficient algorithm for maximum parsimony reconciliation in this new DTLOR model and then show how it can be extended to account for non-binary gene trees to handle uncertainty in gene tree topologies. Finally, we describe preliminary experimental results from the integration of our algorithm into the existing xenoGI tool for reconstructing the histories of genomic islands in closely related bacteria.

CONCLUSIONS

Reconciliation in the DTLOR model can offer new insights into the evolution of microbes that is not currently possible under the DTL model.

Gene tree and species tree reconciliation with endosymbiotic gene transfer

MOTIVATION

It is largely established that all extant mitochondria originated from a unique endosymbiotic event integrating an α-proteobacterial genome into an eukaryotic cell. Subsequently, eukaryote evolution has been marked by episodes of gene transfer, mainly from the mitochondria to the nucleus, resulting in a significant reduction of the mitochondrial genome, eventually completely disappearing in some lineages. However, in other lineages such as in land plants, a high variability in gene repertoire distribution, including genes encoded in both the nuclear and mitochondrial genome, is an indication of an ongoing process of Endosymbiotic Gene Transfer (EGT). Understanding how both nuclear and mitochondrial genomes have been shaped by gene loss, duplication and transfer is expected to shed light on a number of open questions regarding the evolution of eukaryotes, including rooting of the eukaryotic tree.

RESULTS

We address the problem of inferring the evolution of a gene family through duplication, loss and EGT events, the latter considered as a special case of horizontal gene transfer occurring between the mitochondrial and nuclear genomes of the same species (in one direction or the other). We consider both EGT events resulting in maintaining (EGTcopy) or removing (EGTcut) the gene copy in the source genome. We present a linear-time algorithm for computing the DLE (Duplication, Loss and EGT) distance, as well as an optimal reconciled tree, for the unitary cost, and a dynamic programming algorithm allowing to output all optimal reconciliations for an arbitrary cost of operations. We illustrate the application of our EndoRex software and analyze different costs settings parameters on a plant dataset and discuss the resulting reconciled trees.

AVAILABILITY AND IMPLEMENTATION

EndoRex implementation and supporting data are available on the GitHub repository via https://github.com/AEVO-lab/EndoRex.

Comprehensive analysis of SSRs and database construction using all complete gene-coding sequences in major horticultural and representative plants

Simple sequence repeats (SSRs) are one of the most important genetic markers and widely exist in most species. Here, we identified 249,822 SSRs from 3,951,919 genes in 112 plants. Then, we conducted a comprehensive analysis of these SSRs and constructed a plant SSR database (PSSRD). Interestingly, more SSRs were found in lower plants than in higher plants, showing that lower plants needed to adapt to early extreme environments. Four specific enriched functional terms in the lower plant Chlamydomonas reinhardtii were detected when it was compared with seven other higher plants. In addition, Guanylate_cyc existed in more genes of lower plants than of higher plants. In our PSSRD, we constructed an interactive plotting function in the chart interface, and users can easily view the detailed information of SSRs. All SSR information, including sequences, primers, and annotations, can be downloaded from our database. Moreover, we developed Web SSR Finder and Batch SSR Finder tools, which can be easily used for identifying SSRs. Our database was developed using PHP, HTML, JavaScript, and MySQL, which are freely available at http://www.pssrd.info/ . We conducted an analysis of the Myb gene families and flowering genes as two applications of the PSSRD. Further analysis indicated that whole-genome duplication and whole-genome triplication played a major role in the expansion of the Myb gene families. These SSR markers in our database will greatly facilitate comparative genomics and functional genomics studies in the future.

Genome-wide signatures of plastid-nuclear coevolution point to repeated perturbations of plastid proteostasis systems across angiosperms

Nuclear and plastid (chloroplast) genomes experience different mutation rates, levels of selection, and transmission modes, yet key cellular functions depend on their coordinated interactions. Functionally related proteins often show correlated changes in rates of sequence evolution across a phylogeny [evolutionary rate covariation (ERC)], offering a means to detect previously unidentified suites of coevolving and cofunctional genes. We performed phylogenomic analyses across angiosperm diversity, scanning the nuclear genome for genes that exhibit ERC with plastid genes. As expected, the strongest hits were highly enriched for genes encoding plastid-targeted proteins, providing evidence that cytonuclear interactions affect rates of molecular evolution at genome-wide scales. Many identified nuclear genes functioned in post-transcriptional regulation and the maintenance of protein homeostasis (proteostasis), including protein translation (in both the plastid and cytosol), import, quality control, and turnover. We also identified nuclear genes that exhibit strong signatures of coevolution with the plastid genome, but their encoded proteins lack organellar-targeting annotations, making them candidates for having previously undescribed roles in plastids. In sum, our genome-wide analyses reveal that plastid-nuclear coevolution extends beyond the intimate molecular interactions within chloroplast enzyme complexes and may be driven by frequent rewiring of the machinery responsible for maintenance of plastid proteostasis in angiosperms.

Comprehensive identification and characterization of simple sequence repeats based on the whole-genome sequences of 14 forest and fruit trees

Simple sequence repeats (SSRs) are popular and important molecular markers that exist widely in plants. Here, we conducted a comprehensive identification and comparative analysis of SSRs in 14 tree species. A total of 16, 298 SSRs were identified from 429, 449 genes, and primers were successfully designed for 99.44% of the identified SSRs. Our analysis indicated that tri-nucleotide SSRs were the most abundant, with an average of ~834 per species. Functional enrichment analysis by combining SSR-containing genes in all species, revealed 50 significantly enriched terms, with most belonging to transcription factor families associated with plant development and abiotic stresses such as Myeloblastosis_DNA-bind_4 (Myb_DNA-bind_4), APETALA2 (AP2), and Fantastic Four meristem regulator (FAF). Further functional enrichment analysis showed that 48 terms related to abiotic stress regulation and floral development were significantly enriched in ten species, whereas no significantly enriched terms were found in four species. Interestingly, the largest number of enriched terms was detected in Citrus sinensis (L.) Osbeck, accounting for 54.17% of all significantly enriched functional terms. Finally, we analyzed AP2 and trihelix gene families (Myb_DNA-bind_4) due to their significant enrichment in SSR-containing genes. The results indicated that whole-genome duplication (WGD) and whole genome triplication (WGT) might have played major roles in the expansion of the AP2 gene family but only slightly affected the expansion of the trihelix gene family during evolution. In conclusion, the identification and comprehensive characterization of SSR markers will greatly facilitate future comparative genomics and functional genomics studies.

Ten Years of Collaborative Progress in the Quest for Orthologs

Accurate determination of the evolutionary relationships between genes is a foundational challenge in biology. Homology-evolutionary relatedness-is in many cases readily determined based on sequence similarity analysis. By contrast, whether or not two genes directly descended from a common ancestor by a speciation event (orthologs) or duplication event (paralogs) is more challenging, yet provides critical information on the history of a gene. Since 2009, this task has been the focus of the Quest for Orthologs (QFO) Consortium. The sixth QFO meeting took place in Okazaki, Japan in conjunction with the 67th National Institute for Basic Biology conference. Here, we report recent advances, applications, and oncoming challenges that were discussed during the conference. Steady progress has been made toward standardization and scalability of new and existing tools. A feature of the conference was the presentation of a panel of accessible tools for phylogenetic profiling and several developments to bring orthology beyond the gene unit-from domains to networks. This meeting brought into light several challenges to come: leveraging orthology computations to get the most of the incoming avalanche of genomic data, integrating orthology from domain to biological network levels, building better gene models, and adapting orthology approaches to the broad evolutionary and genomic diversity recognized in different forms of life and viruses.

Genomic insights into the host specific adaptation of the Pneumocystis genus

Pneumocystis jirovecii, the fungal agent of human Pneumocystis pneumonia, is closely related to macaque Pneumocystis. Little is known about other Pneumocystis species in distantly related mammals, none of which are capable of establishing infection in humans. The molecular basis of host specificity in Pneumocystis remains unknown as experiments are limited due to an inability to culture any species in vitro. To explore Pneumocystis evolutionary adaptations, we have sequenced the genomes of species infecting macaques, rabbits, dogs and rats and compared them to available genomes of species infecting humans, mice and rats. Complete whole genome sequence data enables analysis and robust phylogeny, identification of important genetic features of the host adaptation, and estimation of speciation timing relative to the rise of their mammalian hosts. Our data reveals insights into the evolution of P. jirovecii, the sole member of the genus able to infect humans.

Cissé, Ma et al. utilize genomic data from Pneumocystis species infecting macaques, rabbit, dogs and rats to investigate the molecular basis of host specificity in Pneumocystis. Their analyses provide insight to the specific adaptations enabling the infection of humans by P. jirovecii.

Treerecs: an integrated phylogenetic tool, from sequences to reconciliations

MOTIVATION

Gene and species tree reconciliation methods are used to interpret gene trees, root them and correct uncertainties that are due to scarcity of signal in multiple sequence alignments. So far, reconciliation tools have not been integrated in standard phylogenetic software and they either lack performance on certain functions, or usability for biologists.

RESULTS

We present Treerecs, a phylogenetic software based on duplication-loss reconciliation. Treerecs is simple to install and to use. It is fast and versatile, has a graphic output, and can be used along with methods for phylogenetic inference on multiple alignments like PLL and Seaview.

AVAILABILITY AND IMPLEMENTATION

Treerecs is open-source. Its source code (C++, AGPLv3) and manuals are available from https://project.inria.fr/treerecs/.

Light environment drives evolution of color vision genes in butterflies and moths

Opsins, combined with a chromophore, are the primary light-sensing molecules in animals and are crucial for color vision. Throughout animal evolution, duplications and losses of opsin proteins are common, but it is unclear what is driving these gains and losses. Light availability is implicated, and dim environments are often associated with low opsin diversity and loss. Correlations between high opsin diversity and bright environments, however, are tenuous. To test if increased light availability is associated with opsin diversification, we examined diel niche and identified opsins using transcriptomes and genomes of 175 butterflies and moths (Lepidoptera). We found 14 independent opsin duplications associated with bright environments. Estimating their rates of evolution revealed that opsins from diurnal taxa evolve faster-at least 13 amino acids were identified with higher dN/dS rates, with a subset close enough to the chromophore to tune the opsin. These results demonstrate that high light availability increases opsin diversity and evolution rate in Lepidoptera.

Capybara: equivalence ClAss enumeration of coPhylogenY event-BAsed ReconciliAtions

MOTIVATION

Phylogenetic tree reconciliation is the method of choice in analyzing host-symbiont systems. Despite the many reconciliation tools that have been proposed in the literature, two main issues remain unresolved: (i) listing suboptimal solutions (i.e. whose score is 'close' to the optimal ones) and (ii) listing only solutions that are biologically different 'enough'. The first issue arises because the optimal solutions are not always the ones biologically most significant; providing many suboptimal solutions as alternatives for the optimal ones is thus very useful. The second one is related to the difficulty to analyze an often huge number of optimal solutions. In this article, we propose Capybara that addresses both of these problems in an efficient way. Furthermore, it includes a tool for visualizing the solutions that significantly helps the user in the process of analyzing the results.

AVAILABILITY AND IMPLEMENTATION

The source code, documentation and binaries for all platforms are freely available at https://capybara-doc.readthedocs.io/.

CONTACT

yishu.wang@univ-lyon1.fr or blerina.sinaimeri@inria.fr.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Chromosomal scale assembly of parasitic wasp genome reveals symbiotic virus colonization

Endogenous viruses form an important proportion of eukaryote genomes and a source of novel functions. How large DNA viruses integrated into a genome evolve when they confer a benefit to their host, however, remains unknown. Bracoviruses are essential for the parasitism success of parasitoid wasps, into whose genomes they integrated ~103 million years ago. Here we show, from the assembly of a parasitoid wasp genome at a chromosomal scale, that bracovirus genes colonized all ten chromosomes of Cotesia congregata. Most form clusters of genes involved in particle production or parasitism success. Genomic comparison with another wasp, Microplitis demolitor, revealed that these clusters were already established ~53 mya and thus belong to remarkably stable genomic structures, the architectures of which are evolutionary constrained. Transcriptomic analyses highlight temporal synchronization of viral gene expression without resulting in immune gene induction, suggesting that no conflicts remain between ancient symbiotic partners when benefits to them converge.

Genome-Wide Identification and Comparative Analysis of ARF Family Genes in Three Apiaceae Species

The family Apiaceae includes many important vegetables and medicinal plants. Auxin response factors (ARFs) play critical roles in regulating plant growth and development. Here, we performed a comprehensive analysis of the ARF gene family in three Apiaceae species, celery, coriander, and carrot, and compared the results with the ARF gene family of lettuce, Arabidopsis, and grape. We identified 156 ARF genes in all six species and 89 genes in the three Apiaceae species, including 28, 34, and 27 in celery, coriander, and carrot, respectively. The paralogous gene number in coriander was far greater than that in carrot and celery. Our analysis revealed that ARF genes of the three Apiaceae species in 34 branches of the phylogenetic tree underwent significant positive selection. Additionally, our findings indicated that whole-genome duplication played an important role in ARF gene family expansion. Coriander contained a greater number of ARF genes than celery and carrot because of more gene duplications and less gene losses. We also analyzed the expression of ARF genes in three tissues by RNA-seq and verified the results by quantitative real-time PCR. Furthermore, we found that several paralogous genes exhibited divergent expression patterns. Overall, this study provides a valuable resource for exploring how ARF family genes regulate plant growth and development in other plants. Since this is the first report of the ARF gene family in Apiaceae, our results will serve as a guide for comparative and functional analyses of ARF and other gene families in Apiaceae.

Diversification and selection pattern of CYP6B genes in Japanese Papilio butterflies and their association with host plant spectra

Herbivorous insects are thought to have evolved counteradaptations to conquer chemical defenses in their host plants in a stepwise co-evolutionary process. Papilio butterflies use CYP6B gene family members to metabolize furanocoumarins in their Rutaceae or Apiaceae host plants. CYP6Bs have functionally diverged among Papilio species to be able to metabolite diverse types of furanocoumarins in their host plants. In this study, we examined the diversification and selection patterns of CYP6B among nine Papilio species in Japan (eight Rutaceae specialists and one Apiaceae specialist) and their association with host plant spectra and furanocoumarin profiles. We compared host plant spectrum of eight Rutaceae feeding Papilio species and also performed a furanocoumarin profiling of their host plants. In addition, we reconstructed CYP6B gene phylogeny and performed selection analysis based on the transcriptome data of those nine Papilio species. Among Rutaceae-feeding Papilio species, host plant spectrum differences were correlated with their furanocoumarin profiles. However, all tested Papilio species had similar duplicated sets of CYP6B, with no apparent lineage-specific or host plant-specific pattern of CYP6B diversification. Selection analysis showed a signature of positive selection on a CYP6B branch. The positively selected sites located at predicted substrate recognition sites and we also found that these CYP6B genes were observed only in Rutaceae-feeding species. These findings indicate that most CYP6B diversification occurred in ancestral species of these Papilio species, possibly in association with specific host plant chemical defenses and subsequent gene loss due to host specialization. These processes would have shaped the complex diversification patterns of the CYP6B gene family in Papilio butterflies. Our results also show potentially important CYP6B clades among Papilio species which likely to have diverged functions and associated with host plant phytochemicals in ancestral Papilio species.

Evolution of NLR Resistance Genes in Magnoliids: Dramatic Expansions of CNLs and Multiple Losses of TNLs

Magnoliids are the third-largest group of angiosperms and occupy a critical position in angiosperm evolution. In the past years, due to the lack of sequenced genomes, the disease resistance gene (R gene) profile of magnoliids remains poorly understood. By the genome-wide identification of 1,832 NLR genes from seven magnoliid genomes, we built a framework for the evolution of magnoliid R genes. TNL genes were completely absent from five magnoliids, presumably due to immune pathway deficiencies. A total of 74 ancestral R genes (70 CNLs, 3 TNLs, and 1 RNL) were recovered in a common ancestor of magnoliids, from which all current NLR gene repertoires were derived. Tandem duplication served as the major drive for NLR genes expansion in seven magnoliid genomes, as most surveyed angiosperms. Due to recent rapid expansions, most magnoliids exhibited "a first expansion followed by a slight contraction and a further stronger expansion" evolutionary pattern, while both Litsea cubeba and Persea americana showed a two-times-repeated pattern of "expansion followed by contraction." The transcriptome analysis of seven different tissues of Saururus chinensis revealed a low expression of most NLR genes, with some R genes displaying a relatively higher expression in roots and fruits. Overall, our study sheds light on the evolution of NLR genes in magnoliids, compensates for insufficiency in major angiosperm lineages, and provides an important reference for a better understanding of angiosperm NLR genes.

Suppression of wheat blast resistance by an effector of Pyricularia oryzae is counteracted by a host specificity resistance gene in wheat

Wheat blast caused by the Triticum pathotype of Pyricularia oryzae poses a serious threat to wheat production in South America and Asia and is now becoming a pandemic disease. Here, we show that Rmg8, a promising wheat gene for resistance breeding, is suppressed by PWT4, an effector gene of P. oryzae, and in turn that the suppression is counteracted by Rwt4, a wheat gene recognizing PWT4. When PWT4 was introduced into a wheat blast isolate carrying AVR-Rmg8 (an avirulence gene corresponding to Rmg8), PWT4 suppressed wheat resistance conferred by Rmg8. PWT4 did not alter the expression of AVR-Rmg8, but higher expression of PWT4 led to more efficient suppression. This suppression was observed in rwt4 carriers, but not in Rwt4 carriers, indicating that it is counteracted by Rwt4. PWT4 was assumed to have been horizontally transferred from a weed-associated cryptic species, P. pennisetigena, to an Avena isolate of P. oryzae in Brazil. This implies a potential risk of the acquisition of PWT4 by the wheat blast fungus and the 'breakdown' of Rmg8. We suggest that Rmg8 should be introduced together with Rwt4 into a wheat cultivar when it is used for resistance breeding.

Evolution of MicroRNA Biogenesis Genes in the Sterlet (Acipenser ruthenus) and Other Polyploid Vertebrates

MicroRNAs play a crucial role in eukaryotic gene regulation. For a long time, only little was known about microRNA-based gene regulatory mechanisms in polyploid animal genomes due to difficulties of polyploid genome assembly. However, in recent years, several polyploid genomes of fish, amphibian, and even invertebrate species have been sequenced and assembled. Here we investigated several key microRNA-associated genes in the recently sequenced sterlet (Acipenser ruthenus) genome, whose lineage has undergone a whole genome duplication around 180 MYA. We show that two paralogs of drosha, dgcr8, xpo1, and xpo5 as well as most ago genes have been retained after the acipenserid-specific whole genome duplication, while ago1 and ago3 genes have lost one paralog. While most diploid vertebrates possess only a single copy of dicer1, we strikingly found four paralogs of this gene in the sterlet genome, derived from a tandem segmental duplication that occurred prior to the last whole genome duplication. ago1,3,4 and exportins1,5 look to be prone to additional segment duplications producing up to four-five paralog copies in ray-finned fishes. We demonstrate for the first time exon microsatellite amplification in the acipenserid drosha2 gene, resulting in a highly variable protein product, which may indicate sub- or neofunctionalization. Paralogous copies of most microRNA metabolism genes exhibit different expression profiles in various tissues and remain functional despite the rediploidization process. Subfunctionalization of microRNA processing gene paralogs may be beneficial for different pathways of microRNA metabolism. Genetic variability of microRNA processing genes may represent a substrate for natural selection, and, by increasing genetic plasticity, could facilitate adaptations to changing environments.

An investigation of irreproducibility in maximum likelihood phylogenetic inference

Phylogenetic trees are essential for studying biology, but their reproducibility under identical parameter settings remains unexplored. Here, we find that 3515 (18.11%) IQ-TREE-inferred and 1813 (9.34%) RAxML-NG-inferred maximum likelihood (ML) gene trees are topologically irreproducible when executing two replicates (Run1 and Run2) for each of 19,414 gene alignments in 15 animal, plant, and fungal phylogenomic datasets. Notably, coalescent-based ASTRAL species phylogenies inferred from Run1 and Run2 sets of individual gene trees are topologically irreproducible for 9/15 phylogenomic datasets, whereas concatenation-based phylogenies inferred twice from the same supermatrix are reproducible. Our simulations further show that irreproducible phylogenies are more likely to be incorrect than reproducible phylogenies. These results suggest that a considerable fraction of single-gene ML trees may be irreproducible. Increasing reproducibility in ML inference will benefit from providing analyses’ log files, which contain typically reported parameters (e.g., program, substitution model, number of tree searches) but also typically unreported ones (e.g., random starting seed number, number of threads, processor type).

Replicate runs of maximum likelihood phylogenetic analyses can generate different tree topologies due to differences in parameters, such as random seeds. Here, Shen et al. demonstrate that replicate runs can generate substantially different tree topologies even with identical data and parameters.

The Multilocus Multispecies Coalescent: A Flexible New Model of Gene Family Evolution

Incomplete lineage sorting (ILS), the interaction between coalescence and speciation, can generate incongruence between gene trees and species trees, as can gene duplication (D), transfer (T), and loss (L). These processes are usually modeled independently, but in reality, ILS can affect gene copy number polymorphism, that is, interfere with DTL. This has been previously recognized, but not treated in a satisfactory way, mainly because DTL events are naturally modeled forward-in-time, while ILS is naturally modeled backward-in-time with the coalescent. Here, we consider the joint action of ILS and DTL on the gene tree/species tree problem in all its complexity. In particular, we show that the interaction between ILS and duplications/transfers (without losses) can result in patterns usually interpreted as resulting from gene loss, and that the realized rate of D, T, and L becomes nonhomogeneous in time when ILS is taken into account. We introduce algorithmic solutions to these problems. Our new model, the multilocus multispecies coalescent, which also accounts for any level of linkage between loci, generalizes the multispecies coalescent (MSC) model and offers a versatile, powerful framework for proper simulation, and inference of gene family evolution. [Gene duplication; gene loss; horizontal gene transfer; incomplete lineage sorting; multispecies coalescent; hemiplasy; recombination.].

Comprehensive analysis of polygalacturonase genes offers new insights into their origin and functional evolution in land plants

Polygalacturonase (PG) is a hydrolase that participates in pectin degradation, pod shattering and fruit softening. Here, we identified 2786 PG genes across 54 plants, which could be divided into three groups. Evolutionary analysis suggested that PG family originated from the charophyte green algae, and Subgroups A2-A4 evolved from the Subgroup A1 after the tracheophyte-angiosperm split. Whole-genome duplication was the major force leading to PG gene expansion. Interestingly, the PG genes continuously expanded in eudicots, whereas it contracted in monocots after the eudicot-monocot split. PG genes in Group A are expressed at high levels in floral organs, whereas genes in Groups B and C are expressed at high levels in various tissues. Moreover, three BnaPG15 members were found for their potential possibility in pod shattering in Brassica napus. Our results provide new insight into the evolutionary history of PG family, and their potentially functional role in plants.

Phylogenomics of the Andean Tetraploid Clade of the American Amaryllidaceae (Subfamily Amaryllidoideae): Unlocking a Polyploid Generic Radiation Abetted by Continental Geodynamics

One of the two major clades of the endemic American Amaryllidaceae subfam. Amaryllidoideae constitutes the tetraploid-derived (n = 23) Andean-centered tribes, most of which have 46 chromosomes. Despite progress in resolving phylogenetic relationships of the group with plastid and nrDNA, certain subclades were poorly resolved or weakly supported in those previous studies. Sequence capture using anchored hybrid enrichment was employed across 95 species of the clade along with five outgroups and generated sequences of 524 nuclear genes and a partial plastome. Maximum likelihood phylogenetic analyses were conducted on concatenated supermatrices, and coalescent-based species tree analyses were run on the gene trees, followed by hybridization network, age diversification and biogeographic analyses. The four tribes Clinantheae, Eucharideae, Eustephieae, and Hymenocallideae (sister to Clinantheae) are resolved in all analyses with > 90 and mostly 100% support, as are almost all genera within them. Nuclear gene supermatrix and species tree results were largely in concordance; however, some instances of cytonuclear discordance were evident. Hybridization network analysis identified significant reticulation in Clinanthus, Hymenocallis, Stenomesson and the subclade of Eucharideae comprising Eucharis, Caliphruria, and Urceolina. Our data support a previous treatment of the latter as a single genus, Urceolina, with the addition of Eucrosia dodsonii. Biogeographic analysis and penalized likelihood age estimation suggests an origin in the Cauca, Desert and Puna Neotropical bioprovinces for the complex in the mid-Oligocene, with more dispersals than vicariances in its history, but no extinctions. Hymenocallis represents the only instance of long-distance vicariance from the tropical Andean origin of its tribe Hymenocallideae. The absence of extinctions correlates with the lack of diversification rate shifts within the clade. The Eucharideae experienced a sudden lineage radiation ca. 10 Mya. We tie much of the divergences in the Andean-centered lineages to the rise of the Andes, and suggest that the Amotape-Huancabamba Zone functioned as both a corridor (dispersal) and a barrier to migration (vicariance). Several taxonomic changes are made. This is the largest DNA sequence data set to be applied within Amaryllidaceae to date.

Diversity and evolutionary dynamics of spore-coat proteins in spore-forming species of Bacillales

Among members of the Bacillales order, there are several species capable of forming a structure called an endospore. Endospores enable bacteria to survive under unfavourable growth conditions and germinate when environmental conditions are favourable again. Spore-coat proteins are found in a multilayered proteinaceous structure encasing the spore core and the cortex. They are involved in coat assembly, cortex synthesis and germination. Here, we aimed to determine the diversity and evolutionary processes that have influenced spore-coat genes in various spore-forming species of Bacillales using an in silico approach. For this, we used sequence similarity searching algorithms to determine the diversity of coat genes across 161 genomes of Bacillales. The results suggest that among Bacillales, there is a well-conserved core genome, composed mainly by morphogenetic coat proteins and spore-coat proteins involved in germination. However, some spore-coat proteins are taxa-specific. The best-conserved genes among different species may promote adaptation to changeable environmental conditions. Because most of the Bacillus species harbour complete or almost complete sets of spore-coat genes, we focused on this genus in greater depth. Phylogenetic reconstruction revealed eight monophyletic groups in the Bacillus genus, of which three are newly discovered. We estimated the selection pressures acting over spore-coat genes in these monophyletic groups using classical and modern approaches and detected horizontal gene transfer (HGT) events, which have been further confirmed by scanning the genomes to find traces of insertion sequences. Although most of the genes are under purifying selection, there are several cases with individual sites evolving under positive selection. Finally, the HGT results confirm that sporulation is an ancestral feature in Bacillus.

Origin and Evolutionary Dynamics of the miR2119 and ADH1 Regulatory Module in Legumes

MicroRNAs are important regulators of gene expression in eukaryotes. Previously, we reported that in Phaseolus vulgaris, the precursor for miR2119 is located in the same gene as miR398a, conceiving a dicistronic MIR gene. Both miRNA precursors are transcribed and processed from a single transcript resulting in two mature microRNAs that regulate the mRNAs encoding ALCOHOL DEHYDROGENASE 1 (ADH1) and COPPER-ZINC SUPEROXIDE DISMUTASE 1 (CSD1). Genes for miR398 are distributed throughout the spermatophytes; however, miR2119 is only found in Leguminosae species, indicating its recent emergence. Here, we used public databases to explore the presence of the miR2119 sequence in several plant species. We found that miR2119 is present only in specific clades within the Papilionoideae subfamily, including important crops used for human consumption and forage. Within this subfamily, MIR2119 and MIR398a are found together as a single gene in the genomes of the Millettioids and Hologalegina. In contrast, in the Dalbergioids MIR2119 is located in a different locus from MIR398a, suggesting this as the ancestral genomic organization. To our knowledge, this is a unique example where two separate MIRNA genes have merged to generate a single polycistronic gene. Phylogenetic analysis of ADH1 gene sequences in the Papilionoideae subfamily revealed duplication events resulting in up to four ADH1 genes in certain species. Notably, the presence of MIR2119 correlates with the conservation of target sites in particular ADH1 genes in each clade. Our results suggest that post-transcriptional regulation of ADH1 genes by miR2119 has contributed to shaping the expansion and divergence of this gene family in the Papilionoideae. Future experimental work on ADH1 regulation by miR2119 in more legume species will help to further understand the evolutionary history of the ADH1 gene family and the relevance of miRNA regulation in this process.

Minimizing genomic duplication episodes

BACKGROUND

The genomic duplication study is fundamental to understand the process of evolution. In evolutionary molecular biology, many approaches focus on discovering the occurrences of gene duplications and multiple gene duplication episodes and their locations in the Tree of Life. To reconstruct such episodes, one can cluster single gene duplications inferred by reconciling a set of gene trees with a species tree.

RESULTS

We propose an efficient quadratic time algorithm to solve the problem of genomic duplication clustering, in which input gene trees are rooted, episode locations are restricted to preserve the minimal number of single gene duplications, clustering rules are described by minimum episodes method, and the goal is based on the recently introduced new approach to minimize the maximal number of duplication episodes on a single path, called here the MP score. Based on our theoretical results, we show new algorithmic relationships between the MP score and the minimum episodes (ME) score, defined as the minimal number of duplication episodes.

CONCLUSIONS

Our evaluation analysis on three empirical datasets demonstrates, that under the model in which the minimal number of duplications is preserved, the duplication clusterings with minimal MP score support the clusterings with the minimal total number of duplication episodes.

AVAILABILITY

The software is available at https://bitbucket.org/pgor17/rmp.

Gene Loss and Evolution of the Plastome

Chloroplasts are unique organelles within the plant cells and are responsible for sustaining life forms on the earth due to their ability to conduct photosynthesis. Multiple functional genes within the chloroplast are responsible for a variety of metabolic processes that occur in the chloroplast. Considering its fundamental role in sustaining life on the earth, it is important to identify the level of diversity present in the chloroplast genome, what genes and genomic content have been lost, what genes have been transferred to the nuclear genome, duplication events, and the overall origin and evolution of the chloroplast genome. Our analysis of 2511 chloroplast genomes indicated that the genome size and number of coding DNA sequences (CDS) in the chloroplasts genome of algae are higher relative to other lineages. Approximately 10.31% of the examined species have lost the inverted repeats (IR) in the chloroplast genome that span across all the lineages. Genome-wide analyses revealed the loss of the Rbcl gene in parasitic and heterotrophic plants occurred approximately 56 Ma ago. PsaM, Psb30, ChlB, ChlL, ChlN, and Rpl21 were found to be characteristic signature genes of the chloroplast genome of algae, bryophytes, pteridophytes, and gymnosperms; however, none of these genes were found in the angiosperm or magnoliid lineage which appeared to have lost them approximately 203-156 Ma ago. A variety of chloroplast-encoded genes were lost across different species lineages throughout the evolutionary process. The Rpl20 gene, however, was found to be the most stable and intact gene in the chloroplast genome and was not lost in any of the analyzed species, suggesting that it is a signature gene of the plastome. Our evolutionary analysis indicated that chloroplast genomes evolved from multiple common ancestors ~1293 Ma ago and have undergone vivid recombination events across different taxonomic lineages.

Biased Gene Retention in the Face of Introgression Obscures Species Relationships

Phylogenomic analyses are recovering previously hidden histories of hybridization, revealing the genomic consequences of these events on the architecture of extant genomes. We applied phylogenomic techniques and several complementary statistical tests to show that introgressive hybridization appears to have occurred between close relatives of Arabidopsis, resulting in cytonuclear discordance and impacting our understanding of species relationships in the group. The composition of introgressed and retained genes indicates that selection against incompatible cytonuclear and nuclear-nuclear interactions likely acted during introgression, whereas linkage also contributed to genome composition through the retention of ancient haplotype blocks. We also applied divergence-based tests to determine the species branching order and distinguish donor from recipient lineages. Surprisingly, these analyses suggest that cytonuclear discordance arose via extensive nuclear, rather than cytoplasmic, introgression. If true, this would mean that most of the nuclear genome was displaced during introgression whereas only a small proportion of native alleles were retained.

The Origin of the Legumes is a Complex Paleopolyploid Phylogenomic Tangle Closely Associated with the Cretaceous-Paleogene (K-Pg) Mass Extinction Event

The consequences of the Cretaceous–Paleogene (K–Pg) boundary (KPB) mass extinction for the evolution of plant diversity remain poorly understood, even though evolutionary turnover of plant lineages at the KPB is central to understanding assembly of the Cenozoic biota. The apparent concentration of whole genome duplication (WGD) events around the KPB may have played a role in survival and subsequent diversification of plant lineages. To gain new insights into the origins of Cenozoic biodiversity, we examine the origin and early evolution of the globally diverse legume family (Leguminosae or Fabaceae). Legumes are ecologically (co-)dominant across many vegetation types, and the fossil record suggests that they rose to such prominence after the KPB in parallel with several well-studied animal clades including Placentalia and Neoaves. Furthermore, multiple WGD events are hypothesized to have occurred early in legume evolution. Using a recently inferred phylogenomic framework, we investigate the placement of WGDs during early legume evolution using gene tree reconciliation methods, gene count data and phylogenetic supernetwork reconstruction. Using 20 fossil calibrations we estimate a revised timeline of legume evolution based on 36 nuclear genes selected as informative and evolving in an approximately clock-like fashion. To establish the timing of WGDs we also date duplication nodes in gene trees. Results suggest either a pan-legume WGD event on the stem lineage of the family, or an allopolyploid event involving (some of) the earliest lineages within the crown group, with additional nested WGDs subtending subfamilies Papilionoideae and Detarioideae. Gene tree reconciliation methods that do not account for allopolyploidy may be misleading in inferring an earlier WGD event at the time of divergence of the two parental lineages of the polyploid, suggesting that the allopolyploid scenario is more likely. We show that the crown age of the legumes dates to the Maastrichtian or early Paleocene and that, apart from the Detarioideae WGD, paleopolyploidy occurred close to the KPB. We conclude that the early evolution of the legumes followed a complex history, in which multiple auto- and/or allopolyploidy events coincided with rapid diversification and in association with the mass extinction event at the KPB, ultimately underpinning the evolutionary success of the Leguminosae in the Cenozoic. [Allopolyploidy; Cretaceous–Paleogene (K–Pg) boundary; Fabaceae, Leguminosae; paleopolyploidy; phylogenomics; whole genome duplication events]

Distinct Evolutionary Patterns of NBS-Encoding Genes in Three Soapberry Family (Sapindaceae) Species

Nucleotide-binding site (NBS)-type disease resistance genes (R genes) play key roles in plant immune responses and have co-evolved with pathogens over the course of plant lifecycles. Comparative genomic studies tracing the dynamic evolution of NBS-encoding genes have been conducted using many important plant lineages. However, studies on Sapindaceae species have not been performed. In this study, a discrepant number of NBS-encoding genes were identified in the genomes of Xanthoceras sorbifolium (180), Dinnocarpus longan (568), and Acer yangbiense (252). These genes were unevenly distributed and usually clustered as tandem arrays on chromosomes, with few existed as singletons. The phylogenetic analysis revealed that NBS-encoding genes formed three monophyletic clades, RPW8-NBS-LRR (RNL), TIR-NBS-LRR (TNL), and CC-NBS-LRR (CNL), which were distinguished by amino acid motifs. The NBS-encoding genes of the X. sorbifolium, D. longan, and A. yangbiense genomes were derived from 181 ancestral genes (three RNL, 23 TNL, and 155 CNL), which exhibited dynamic and distinct evolutionary patterns due to independent gene duplication/loss events. Specifically, X. sorbifolium exhibited a “first expansion and then contraction” evolutionary pattern, while A. yangbiense and D. longan exhibited a “first expansion followed by contraction and further expansion” evolutionary pattern. However, further expansion in D. longan was stronger than in A. yangbiense after divergence, suggesting that D. longan gained more genes in response to various pathogens. Additionally, the ancient and recent expansion of CNL genes generated the dominance of this subclass in terms of gene numbers, while the low copy number status of RNL genes was attributed to their conserved functions.

Identification and distribution of gene clusters required for synthesis of sphingolipid metabolism inhibitors in diverse species of the filamentous fungus Fusarium

BACKGROUND

Sphingolipids are structural components and signaling molecules in eukaryotic membranes, and many organisms produce compounds that inhibit sphingolipid metabolism. Some of the inhibitors are structurally similar to the sphingolipid biosynthetic intermediate sphinganine and are referred to as sphinganine-analog metabolites (SAMs). The mycotoxins fumonisins, which are frequent contaminants in maize, are one family of SAMs. Due to food and feed safety concerns, fumonisin biosynthesis has been investigated extensively, including characterization of the fumonisin biosynthetic gene cluster in the agriculturally important fungi Aspergillus and Fusarium. Production of several other SAMs has also been reported in fungi, but there is almost no information on their biosynthesis. There is also little information on how widely SAM production occurs in fungi or on the extent of structural variation of fungal SAMs.

RESULTS

Using fumonisin biosynthesis as a model, we predicted that SAM biosynthetic gene clusters in fungi should include a polyketide synthase (PKS), an aminotransferase and a dehydrogenase gene. Surveys of genome sequences identified five putative clusters with this three-gene combination in 92 of 186 Fusarium species examined. Collectively, the putative SAM clusters were distributed widely but discontinuously among the species. We propose that the SAM5 cluster confers production of a previously reported Fusarium SAM, 2-amino-14,16-dimethyloctadecan-3-ol (AOD), based on the occurrence of AOD production only in species with the cluster and on deletion analysis of the SAM5 cluster PKS gene. We also identified SAM clusters in 24 species of other fungal genera, and propose that one of the clusters confers production of sphingofungin, a previously reported Aspergillus SAM.

CONCLUSION

Our results provide a genomics approach to identify novel SAM biosynthetic gene clusters in fungi, which should in turn contribute to identification of novel SAMs with applications in medicine and other fields. Information about novel SAMs could also provide insights into the role of SAMs in the ecology of fungi. Such insights have potential to contribute to strategies to reduce fumonisin contamination in crops and to control crop diseases caused by SAM-producing fungi.

Repeated horizontal gene transfers triggered parallel evolution of magnetotaxis in two evolutionary divergent lineages of magnetotactic bacteria

Under the same selection pressures, two genetically divergent populations may evolve in parallel toward the same adaptive solutions. Here, we hypothesized that magnetotaxis (i.e., magnetically guided chemotaxis) represents a key adaptation to micro-oxic habitats in aquatic sediments and that its parallel evolution homogenized the phenotypes of two evolutionary divergent clusters of freshwater spirilla. All magnetotactic bacteria affiliated to the Magnetospirillum genus (Alphaproteobacteria class) biomineralize the same magnetic particle chains and share highly similar physiological and ultrastructural features. We looked for the processes that could have contributed at shaping such an evolutionary pattern by reconciling species and gene trees using newly sequenced genomes of Magnetospirillum related bacteria. We showed that repeated horizontal gene transfers and homologous recombination of entire operons contributed to the parallel evolution of magnetotaxis. We propose that such processes could represent a more parsimonious and rapid solution for adaptation compared with independent and repeated de novo mutations, especially in the case of traits as complex as magnetotaxis involving tens of interacting proteins. Besides strengthening the idea about the importance of such a function in micro-oxic habitats, these results reinforce previous observations in experimental evolution suggesting that gene flow could alleviate clonal interference and speed up adaptation under some circumstances.

Phylogenetics Identifies Two Eumetazoan TRPM Clades and an Eighth TRP Family, TRP Soromelastatin (TRPS)

Transient receptor potential melastatins (TRPMs) are most well known as cold and menthol sensors, but are in fact broadly critical for life, from ion homeostasis to reproduction. Yet, the evolutionary relationship between TRPM channels remains largely unresolved, particularly with respect to the placement of several highly divergent members. To characterize the evolution of TRPM and like channels, we performed a large-scale phylogenetic analysis of >1,300 TRPM-like sequences from 14 phyla (Annelida, Arthropoda, Brachiopoda, Chordata, Cnidaria, Echinodermata, Hemichordata, Mollusca, Nematoda, Nemertea, Phoronida, Priapulida, Tardigrada, and Xenacoelomorpha), including sequences from a variety of recently sequenced genomes that fill what would otherwise be substantial taxonomic gaps. These findings suggest: 1) the previously recognized TRPM family is in fact two distinct families, including canonical TRPM channels and an eighth major previously undescribed family of animal TRP channel, TRP soromelastatin; 2) two TRPM clades predate the last bilaterian-cnidarian ancestor; and 3) the vertebrate-centric trend of categorizing TRPM channels as 1-8 is inappropriate for most phyla, including other chordates.

Convergent Loss of an EDS1/PAD4 Signaling Pathway in Several Plant Lineages Reveals Coevolved Components of Plant Immunity and Drought Response

Plant innate immunity relies on nucleotide binding leucine-rich repeat receptors (NLRs) that recognize pathogen-derived molecules and activate downstream signaling pathways. We analyzed the variation in NLR gene copy number and identified plants with a low number of NLR genes relative to sister species. We specifically focused on four plants from two distinct lineages, one monocot lineage (Alismatales) and one eudicot lineage (Lentibulariaceae). In these lineages, the loss of NLR genes coincides with loss of the well-known downstream immune signaling complex ENHANCED DISEASE SUSCEPTIBILITY 1 (EDS1)/PHYTOALEXIN DEFICIENT 4 (PAD4). We expanded our analysis across whole proteomes and found that other characterized immune genes were absent only in Lentibulariaceae and Alismatales. Additionally, we identified genes of unknown function that were convergently lost together with EDS1/PAD4 in five plant species. Gene expression analyses in Arabidopsis (Arabidopsis thaliana) and Oryza sativa revealed that several homologs of the candidates are differentially expressed during pathogen infection, drought, and abscisic acid treatment. Our analysis provides evolutionary evidence for the rewiring of plant immunity in some plant lineages, as well as the coevolution of the EDS1/PAD4 pathway and drought responses.

Cryo-Electron Tomography Reveals the Complex Ultrastructural Organization of Multicellular Filamentous Chloroflexota (Chloroflexi) Bacteria

The cell biology of Chloroflexota is poorly studied. We applied cryo-focused ion beam milling and cryo-electron tomography to study the ultrastructural organization of thermophilic Roseiflexus castenholzii and Chloroflexus aggregans, and mesophilic “Ca. Viridilinea mediisalina.” These species represent the three main lineages within a group of multicellular filamentous anoxygenic phototrophic Chloroflexota bacteria belonging to the Chloroflexales order. We found surprising structural complexity in the Chloroflexales. As with filamentous cyanobacteria, cells of C. aggregans and “Ca. Viridilinea mediisalina” share the outer membrane-like layers of their intricate multilayer cell envelope. Additionally, cells of R. castenholzii and “Ca. Viridilinea mediisalina” are connected by septal channels that resemble cyanobacterial septal junctions. All three strains possess long pili anchored close to cell-to-cell junctions, a morphological feature comparable to that observed in cyanobacteria. The cytoplasm of the Chloroflexales bacteria is crowded with intracellular organelles such as different types of storage granules, membrane vesicles, chlorosomes, gas vesicles, chemoreceptor-like arrays, and cytoplasmic filaments. We observed a higher level of complexity in the mesophilic strain compared to the thermophilic strains with regards to the composition of intracellular bodies and the organization of the cell envelope. The ultrastructural details that we describe in these Chloroflexales bacteria will motivate further cell biological studies, given that the function and evolution of the many discovered morphological traits remain enigmatic in this diverse and widespread bacterial group.

Diversity and Evolutionary Dynamics of Antiphage Defense Systems in Ralstonia solanacearum Species Complex

Over the years, many researchers have reported a great diversity of bacteriophages infecting members of the Ralstonia solanacearum species complex (RSSC). This diversity has driven bacterial evolution by leading the emergence and maintenance of bacterial defense systems to combat phage infection. In this work, we present an in silico study of the arsenal of defense systems that RSSC harbors and their evolutionary history. For this purpose, we used a combination of genomic, phylogenetic and associative methods. We found that in addition to the CRISPR-Cas system already reported, there are eight other antiphage defense systems including the well-known Restriction-Modification and Toxin-Antitoxin systems. Furthermore, we found a tenth defense system, which is dedicated to reducing the incidence of plasmid transformation in bacteria. We undertook an analysis of the gene gain and loss patterns of the defense systems in 15 genomes of RSSC. Results indicate that the dynamics are inclined toward the gain of defense genes as opposed to the rest of the genes that were preferably lost throughout evolution. This was confirmed by evidence on independent gene acquisition that has occurred by profuse horizontal transfer. The mutation and recombination rates were calculated as a proxy of evolutionary rates. Again, genes encoding the defense systems follow different rates of evolution respect to the rest of the genes. These results lead us to conclude that the evolution of RSSC defense systems is highly dynamic and responds to a different evolutionary regime than the rest of the genes in the genomes of RSSC.

Ancestral Genomes: a resource for reconstructed ancestral genes and genomes across the tree of life

A growing number of whole genome sequencing projects, in combination with development of phylogenetic methods for reconstructing gene evolution, have provided us with a window into genomes that existed millions, and even billions, of years ago. Ancestral Genomes (http://ancestralgenomes.org) is a resource for comprehensive reconstructions of these ‘fossil genomes’. Comprehensive sets of protein-coding genes have been reconstructed for 78 genomes of now-extinct species that were the common ancestors of extant species from across the tree of life. The reconstructed genes are based on the extensive library of over 15 000 gene family trees from the PANTHER database, and are updated on a yearly basis. For each ancestral gene, we assign a stable identifier, and provide additional information designed to facilitate analysis: an inferred name, a reconstructed protein sequence, a set of inferred Gene Ontology (GO) annotations, and a ‘proxy gene’ for each ancestral gene, defined as the least-diverged descendant of the ancestral gene in a given extant genome. On the Ancestral Genomes website, users can browse the Ancestral Genomes by selecting nodes in a species tree, and can compare an extant genome with any of its reconstructed ancestors to understand how the genome evolved.

Horizontal Gene Transfer Phylogenetics: A Random Walk Approach

The dramatic decrease in time and cost for generating genetic sequence data has opened up vast opportunities in molecular systematics, one of which is the ability to decipher the evolutionary history of strains of a species. Under this fine systematic resolution, the standard markers are too crude to provide a phylogenetic signal. Nevertheless, among prokaryotes, genome dynamics in the form of horizontal gene transfer (HGT) between organisms and gene loss seem to provide far richer information by affecting both gene order and gene content. The "synteny index" (SI) between a pair of genomes combines these latter two factors, allowing comparison of genomes with unequal gene content, together with order considerations of their common genes. Although this approach is useful for classifying close relatives, no rigorous statistical modeling for it has been suggested. Such modeling is valuable, as it allows observed measures to be transformed into estimates of time periods during evolution, yielding the "additivity" of the measure. To the best of our knowledge, there is no other additivity proof for other gene order/content measures under HGT. Here, we provide a first statistical model and analysis for the SI measure. We model the "gene neighborhood" as a "birth-death-immigration" process affected by the HGT activity over the genome, and analytically relate the HGT rate and time to the expected SI. This model is asymptotic and thus provides accurate results, assuming infinite size genomes. Therefore, we also developed a heuristic model following an "exponential decay" function, accounting for biologically realistic values, which performed well in simulations. Applying this model to 1,133 prokaryotes partitioned to 39 clusters by the rank of genus yields that the average number of genome dynamics events per gene in the phylogenetic depth of genus is around half with significant variability between genera. This result extends and confirms similar results obtained for individual genera in different manners.

A New Paradigm for Identifying Reconciliation-Scenario Altering Mutations Conferring Environmental Adaptation

An important goal in microbial computational genomics is to identify crucial events in the evolution of a gene that severely alter the duplication, loss, and mobilization patterns of the gene within the genomes in which it disseminates. In this article, we formalize this microbiological goal as a new pattern-matching problem in the domain of gene tree and species tree reconciliation, denoted "Reconciliation-Scenario Altering Mutation (RSAM) Discovery." We propose an [Formula: see text] time algorithm to solve this new problem, where m and n are the number of vertices of the input gene tree and species tree, respectively, and k is a user-specified parameter that bounds from above the number of optimal solutions of interest. The algorithm first constructs a hypergraph representing the k highest scoring reconciliation scenarios between the given gene tree and species tree, and then interrogates this hypergraph for subtrees matching a prespecified RSAM pattern. Our algorithm is optimal in the sense that the number of hypernodes in the hypergraph can be lower bounded by [Formula: see text]. We implement the new algorithm as a tool, called RSAM-finder, and demonstrate its application to the identification of RSAMs in toxins and drug resistance elements across a data set spanning hundreds of species.

Two cases of type-a Haemophilus influenzae meningitis within the same week in the same hospital are phylogenetically unrelated but recently exchanged capsule genes

Haemophilus influenzae causes common and sometimes severe adult and pediatric disease including chronic obstructive respiratory disease, otitis media and infections of the central nervous system. Serotype b strains, with a b-type capsule, have been the historical cause of invasive disease, and the introduction of a serotype b-specific vaccine has led to their decline. However, unencapsulated or non-b-type H. influenzae infections are not prevented by the vaccine and appear to be increasing in frequency. Here we report two pediatric cases of severe central nervous system H. influenzae infection presenting to the same hospital in San Diego, California during the same week in January 2016. Due to good vaccine coverage in this part of the world, H. influenzae cases are normally rare and seeing two cases in the same week was unexpected. We thus suspected a recent transmission chain, and possible local outbreak. To test this hypothesis, we isolated and sequenced whole genomes from each patient and placed them in a phylogenetic tree spanning the known diversity of H. influenzae. Surprisingly, we found that the two isolates (SD2016_1 and SD2016_2) belonged to distantly related lineages, suggesting two independent transmission events and ruling out a local outbreak. Despite being distantly related, the two isolates belong to two different lineages that have exchanged capsule loci in the recent past. Therefore, as in other bacterial pathogens, capsule switching by horizontal gene transfer may be an important evolutionary mechanism of vaccine evasion in H. influenzae.

Counting and sampling gene family evolutionary histories in the duplication-loss and duplication-loss-transfer models

Given a set of species whose evolution is represented by a species tree, a gene family is a group of genes having evolved from a single ancestral gene. A gene family evolves along the branches of a species tree through various mechanisms, including-but not limited to-speciation ([Formula: see text]), gene duplication ([Formula: see text]), gene loss ([Formula: see text]), and horizontal gene transfer ([Formula: see text]). The reconstruction of a gene tree representing the evolution of a gene family constrained by a species tree is an important problem in phylogenomics. However, unlike in the multispecies coalescent evolutionary model that considers only speciation and incomplete lineage sorting events, very little is known about the search space for gene family histories accounting for gene duplication, gene loss and horizontal gene transfer (the [Formula: see text]-model). In this work, we introduce the notion of evolutionary histories defined as a binary ordered rooted tree describing the evolution of a gene family, constrained by a species tree in the [Formula: see text]-model. We provide formal grammars describing the set of all evolutionary histories that are compatible with a given species tree, whether it is ranked or unranked. These grammars allow us, using either analytic combinatorics or dynamic programming, to efficiently compute the number of histories of a given size, and also to generate random histories of a given size under the uniform distribution. We apply these tools to obtain exact asymptotics for the number of gene family histories for two species trees, the rooted caterpillar and complete binary tree, as well as estimates of the range of the exponential growth factor of the number of histories for random species trees of size up to 25. Our results show that including horizontal gene transfers induce a dramatic increase of the number of evolutionary histories. We also show that, within ranked species trees, the number of evolutionary histories in the [Formula: see text]-model is almost independent of the species tree topology. These results establish firm foundations for the development of ensemble methods for the prediction of reconciliations.

Evolutionary Histories of Type III Polyketide Synthases in Fungi

Type III polyketide synthases (PKSs) produce secondary metabolites with diverse biological activities, including antimicrobials. While they have been extensively studied in plants and bacteria, only a handful of type III PKSs from fungi has been characterized in the last 15 years. The exploitation of fungal type III PKSs to produce novel bioactive compounds requires understanding the diversity of these enzymes, as well as of their biosynthetic pathways. Here, phylogenetic and reconciliation analyses of 522 type III PKSs from 1,193 fungal genomes revealed complex evolutionary histories with massive gene duplications and losses, explaining their discontinuous distribution in the fungal tree of life. In addition, horizontal gene transfer events from bacteria to fungi and, to a lower extent, between fungi, could be inferred. Ancestral gene duplication events have resulted in the divergence of eight phylogenetic clades. Especially, two clades show ancestral linkage and functional co-evolution between a type III PKS and a reducing PKS genes. Investigation of the occurrence of protein domains in fungal type III PKS predicted gene clusters highlighted the diversity of biosynthetic pathways, likely reflecting a large chemical landscape. Type III PKS genes are most often located next to genes encoding cytochrome P450s, MFS transporters and transcription factors, defining ancestral core gene clusters. This analysis also allowed predicting gene clusters for the characterized fungal type III PKSs and provides working hypotheses for the elucidation of the full biosynthetic pathways. Altogether, our analyses provide the fundamental knowledge to motivate further characterization and exploitation of fungal type III PKS biosynthetic pathways.

Evolution of the extracytoplasmic function σ factor protein family

Understanding transcription has been a central goal of the scientific community for decades. However, much is still unknown, especially concerning how it is regulated. In bacteria, a single DNA-directed RNA-polymerase performs the whole of transcription. It contains multiple subunits, among which the σ factor that confers promoter specificity. Besides the housekeeping σ factor, bacteria encode several alternative σ factors. The most abundant and diverse family of alternative σ factors, the extracytoplasmic function (ECF) family, regulates transcription of genes associated with stressful scenarios, making them key elements of adaptation to specific environmental changes. Despite this, the evolutionary history of ECF σ factors has never been investigated. Here, we report on our analysis of thousands of members of this family. We show that single events are in the origin of alternative modes of regulation of ECF σ factor activity that require partner proteins, but that multiple events resulted in acquisition of regulatory extensions. Moreover, in Bacteroidetes there is a recent duplication of an ecologically relevant gene cluster that includes an ECF σ factor, whereas in Planctomycetes duplication generates distinct C-terminal extensions after fortuitous insertion of the duplicated σ factor. At last, we also demonstrate horizontal transfer of ECF σ factors between soil bacteria.