DendroBLAST: approximate phylogenetic trees in the absence of multiple sequence alignments

Characteristics and phylogenetic distribution of megaplasmids and prediction of a putative chromid in Pseudomonas aeruginosa

Research on megaplasmids that contribute to the spread of antimicrobial resistance (AMR) in Pseudomonas aeruginosa strains has grown in recent years due to the now widely used technologies allowing long-read sequencing. Here, we systematically analyzed distinct and consistent genetic characteristics of megaplasmids found in P. aeruginosa. Our data provide information on their phylogenetic distribution and hypotheses tracing the potential evolutionary paths of megaplasmids. Most of the megaplasmids we found belong to the IncP-2-type, with conserved and syntenic genetic backbones carrying modules of genes associated with chemotaxis apparatus, tellurite resistance and plasmid replication, segregation, and transmission. Extensively variable regions harbor abundant AMR genes, especially those encoding β-lactamases such as VIM-2, IMP-45, and KPC variants, which are high-risk elements in nosocomial infection. IncP-2 megaplasmids act as effective vehicles transmitting AMR genes to diverse regions. One evolutionary model of the origin of megaplasmids claims that chromids can develop from megaplasmids. These chromids have been characterized as an intermediate between a megaplasmid and a chromosome, also containing core genes that can be found on the chromosome but not on the megaplasmid. Using in silico prediction, we identified the "PABCH45 unnamed replicon" as a putative chromid in P. aeruginosa, which shows a much higher similarity and closer phylogenetic relationship to chromosomes than to megaplasmids while also encoding plasmid-like partition genes. We propose that such a chromid could facilitate genome expansion, allowing for more rapid adaptations to novel ecological niches or selective conditions, in comparison to megaplasmids.

Comparative analysis, diversification, and functional validation of plant nucleotide-binding site domain genes

Nucleotide-binding site (NBS) domain genes are one of the superfamily of resistance genes involved in plant responses to pathogens. The current study identified 12,820 NBS-domain-containing genes across 34 species covering from mosses to monocots and dicots. These identified genes are classified into 168 classes with several novel domain architecture patterns encompassing significant diversity among plant species. Several classical (NBS, NBS-LRR, TIR-NBS, TIR-NBS-LRR, etc.) and species-specific structural patterns (TIR-NBS-TIR-Cupin_1-Cupin_1, TIR-NBS-Prenyltransf, Sugar_tr-NBS etc.) were discovered. We observed 603 orthogroups (OGs) with some core (most common orthogroups; OG0, OG1, OG2, etc.) and unique (highly specific to species; OG80, OG82, etc.) OGs with tandem duplications. The expression profiling presented the putative upregulation of OG2, OG6, and OG15 in different tissues under various biotic and abiotic stresses in susceptible and tolerant plants to cotton leaf curl disease (CLCuD). The genetic variation between susceptible (Coker 312) and tolerant (Mac7) Gossypium hirsutum accessions identified several unique variants in NBS genes of Mac7 (6583 variants) and Coker312 (5173 variants). The protein-ligand and proteins-protein interaction showed a strong interaction of some putative NBS proteins with ADP/ATP and different core proteins of the cotton leaf curl disease virus. The silencing of GaNBS (OG2) in resistant cotton through virus-induced gene silencing (VIGS) demonstrated its putative role in virus tittering. The presented study will be further helpful in understanding the plant adaptation mechanism.

Genetic and Functional Diversity Help Explain Pathogenic, Weakly Pathogenic, and Commensal Lifestyles in the Genus Xanthomonas

The genus Xanthomonas has been primarily studied for pathogenic interactions with plants. However, besides host and tissue-specific pathogenic strains, this genus also comprises nonpathogenic strains isolated from a broad range of hosts, sometimes in association with pathogenic strains, and other environments, including rainwater. Based on their incapacity or limited capacity to cause symptoms on the host of isolation, nonpathogenic xanthomonads can be further characterized as commensal and weakly pathogenic. This study aimed to understand the diversity and evolution of nonpathogenic xanthomonads compared to their pathogenic counterparts based on their cooccurrence and phylogenetic relationship and to identify genomic traits that form the basis of a life history framework that groups xanthomonads by ecological strategies. We sequenced genomes of 83 strains spanning the genus phylogeny and identified eight novel species, indicating unexplored diversity. While some nonpathogenic species have experienced a recent loss of a type III secretion system, specifically the hrp2 cluster, we observed an apparent lack of association of the hrp2 cluster with lifestyles of diverse species. We performed association analysis on a large data set of 337 Xanthomonas strains to explain how xanthomonads may have established association with the plants across the continuum of lifestyles from commensals to weak pathogens to pathogens. Presence of distinct transcriptional regulators, distinct nutrient utilization and assimilation genes, transcriptional regulators, and chemotaxis genes may explain lifestyle-specific adaptations of xanthomonads.

Genomic and Transcriptomic Insights into the Genetic Basis of Foam Secretion in Rice Spittlebug Callitettix versicolor

Many animal species produce protective foams, the majority of which exhibit evolutionary adaptability. Although the function and composition of foams have been widely studied, the genetic basis of foam secretion remains unknown. Unlike most species that produce foam under specific situations, spittlebugs continuously secrete foams throughout all nymphal stages. Here, we capitalize on the rice spittlebug (Callitettix versicolor) to explore the genetic basis of foam secretion through genomic and transcriptomic approaches. Our comparative genomic analysis for C. versicolor and eight other insect species reveals 606 species-specific gene families and 66 expanded gene families, associated with carbohydrate and lipid metabolism. These functions are in accordance with the composition of foams secreted by spittlebugs. Transcriptomic analyses of malpighian tubules across developmental stages detected 3192 differentially expressed genes. Enrichment analysis of these genes highlights functions also revealed by our comparative genomic analysis and aligns with previous histochemical and morphological observations of foam secretion. This consistency suggests the important roles of these candidate genes in foam production. Our study not only provides novel insights into the genetic basis of foam secretion in rice spittlebugs but also contributes valuable knowledge for future evolutionary studies of spittlebugs and the development of pest control strategies for C. versicolor.

Strategies for describing myxozoan pathogens, dreadful fish diseases in aquaculture

Myxozoans are obligate endoparasites, cosmopolitan in distribution with both vertebrate and invertebrate hosts. Their myxospores consist of shell valves, polar capsules with coiled polar tubules that are extrudible, and infective amoeboid germs. Myxozoan parasites are most abundant, and due to their increasing number in recent years, they can pose an emerging threat to the fish industry worldwide. Hence, the immediate need is to devise a strategy to understand and detect parasites and parasitism. They may proliferate to different organs with the advancement of infection. This all warrants the development/devising of strategies and results of integrative studies in order to identify these dreadful parasites and resolve taxonomic issues. Different methods whether classical methods including gross morphology or advanced methods such as electron microscopy (SEM, TEM, STEM), Confocal laser scanning microscopy (CLSM), histopathological studies, site preference, host and tissue specificity, a molecular approach using new markers can be clubbed for identification because these parasites are hidden and are difficult to recognize. This group was earlier classified only on the basis of myxospores morphology, but due to the high structural variability of this group advanced methods and approaches have to be implied which can minimize the problems in assigning new species.

Identification of Incomplete Annotations of Biosynthesis Pathways in Rhodophytes Using a Multi-Omics Approach

Rhodophytes (red algae) are an important source of natural products and are, therefore, a current research focus in terms of metabolite production. The recent increase in publicly available Rhodophyte whole genome and transcriptome assemblies provides the resources needed for in silico metabolic pathway analysis. Thus, this study aimed to create a Rhodophyte multi-omics resource, utilising both genomes and transcriptome assemblies with functional annotations to explore Rhodophyte metabolism. The genomes and transcriptomes of 72 Rhodophytes were functionally annotated and integrated with metabolic reconstruction and phylogenetic inference, orthology prediction, and gene duplication analysis to analyse their metabolic pathways. This resource was utilised via two main investigations: the identification of bioactive sterol biosynthesis pathways and the evolutionary analysis of gene duplications for known enzymes. We report that sterol pathways, including campesterol, β-sitosterol, ergocalciferol and cholesterol biosynthesis pathways, all showed incomplete annotated pathways across all Rhodophytes despite prior in vivo studies showing otherwise. Gene duplication analysis revealed high rates of duplication of halide-associated haem peroxidases in Florideophyte algae, which are involved in the biosynthesis of drug-related halogenated secondary metabolites. In summary, this research revealed trends in Rhodophyte metabolic pathways that have been under-researched and require further functional analysis. Furthermore, the high duplication of haem peroxidases and other peroxidase enzymes offers insight into the potential drug development of Rhodophyte halogenated secondary metabolites.

Whole genome assembly and annotation of the endangered Caribbean coral Acropora cervicornis

Coral species in the genus Acropora are key ecological components of coral reefs worldwide and represent the most diverse genus of scleractinian corals. While key species of Indo-Pacific Acropora have annotated genomes, no annotated genome has been published for either of the two species of Caribbean Acropora. Here we present the first fully annotated genome of the endangered Caribbean staghorn coral, Acropora cervicornis. We assembled and annotated this genome using high-fidelity nanopore long-read sequencing with gene annotations validated with mRNA sequencing. The assembled genome size is 318 Mb, with 28,059 validated genes. Comparative genomic analyses with other Acropora revealed unique features in A. cervicornis, including contractions in immune pathways and expansions in signaling pathways. Phylogenetic analysis confirms previous findings showing that A. cervicornis diverged from Indo-Pacific relatives around 41 million years ago, with the closure of the western Tethys Sea, prior to the primary radiation of Indo-Pacific Acropora. This new A. cervicornis genome enriches our understanding of the speciose Acropora and addresses evolutionary inquiries concerning speciation and hybridization in this diverse clade.

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.

Inter-species association mapping links splice site evolution to METTL16 and SNRNP27K

Eukaryotic genes are interrupted by introns that are removed from transcribed RNAs by splicing. Patterns of splicing complexity differ between species, but it is unclear how these differences arise. We used inter-species association mapping with Saccharomycotina species to correlate splicing signal phenotypes with the presence or absence of splicing factors. Here, we show that variation in 5' splice site sequence preferences correlate with the presence of the U6 snRNA N6-methyladenosine methyltransferase METTL16 and the splicing factor SNRNP27K. The greatest variation in 5' splice site sequence occurred at the +4 position and involved a preference switch between adenosine and uridine. Loss of METTL16 and SNRNP27K orthologs, or a single SNRNP27K methionine residue, was associated with a preference for +4 U. These findings are consistent with splicing analyses of mutants defective in either METTL16 or SNRNP27K orthologs and models derived from spliceosome structures, demonstrating that inter-species association mapping is a powerful orthogonal approach to molecular studies. We identified variation between species in the occurrence of two major classes of 5' splice sites, defined by distinct interaction potentials with U5 and U6 snRNAs, that correlates with intron number. We conclude that variation in concerted processes of 5' splice site selection by U6 snRNA is associated with evolutionary changes in splicing signal phenotypes.

Whole genome assemblies of Zophobas morio and Tenebrio molitor

Zophobas morio (=Zophobas atratus) and Tenebrio molitor are darkling beetles with industrial importance due to their use as feeder insects and their apparent ability to biodegrade plastics. High quality genome assemblies were recently reported for both species. Here, we report additional independent Z. morio and T. molitor genome assemblies generated from Nanopore and Illumina data. Following scaffolding against the published genomes, haploid assemblies of 462 Mb (scaffold N90 of 16.8 Mb) and 258 Mb (scaffold N90 of 5.9 Mb) were produced for Z. morio and T. molitor, respectively. Gene prediction led to the prediction of 28,544 and 19,830 genes for Z. morio and T. molitor, respectively. Benchmarking Universal Single Copy Orthologs (BUSCO) analyses suggested that both assemblies have a high level of completeness; 91.5 and 89.0% of the BUSCO endopterygota marker genes were complete in the Z. morio assembly and proteome, respectively, while 99.1 and 92.8% were complete in the T. molitor assembly and proteome, respectively. Phylogenomic analyses of four genera from the family Tenebrionidae yielded phylogenies consistent with those previously constructed based on mitochondrial genomes. Synteny analyses revealed large stretches of macrosynteny across the family Tenebrionidae, as well as numerous within-chromosome rearrangements. Finally, orthogroup analysis identified ∼28,000 gene families across the family Tenebrionidae, of which 8,185 were identified in all five of the analyzed species, and 10,837 were conserved between Z. morio and T. molitor. We expect that the availability of multiple whole genome sequences for Z. morio and T. molitor will facilitate population genetics studies to identify genetic variation associated with industrially relevant phenotypes.

The Sugar Metabolic Model of Aspergillus niger Can Only Be Reliably Transferred to Fungi of Its Phylum

Fungi play a critical role in the global carbon cycle by degrading plant polysaccharides to small sugars and metabolizing them as carbon and energy sources. We mapped the well-established sugar metabolic network of Aspergillus niger to five taxonomically distant species (Aspergillus nidulans, Penicillium subrubescens, Trichoderma reesei, Phanerochaete chrysosporium and Dichomitus squalens) using an orthology-based approach. The diversity of sugar metabolism correlates well with the taxonomic distance of the fungi. The pathways are highly conserved between the three studied Eurotiomycetes (A. niger, A. nidulans, P. subrubescens). A higher level of diversity was observed between the T. reesei and A. niger, and even more so for the two Basidiomycetes. These results were confirmed by integrative analysis of transcriptome, proteome and metabolome, as well as growth profiles of the fungi growing on the corresponding sugars. In conclusion, the establishment of sugar pathway models in different fungi revealed the diversity of fungal sugar conversion and provided a valuable resource for the community, which would facilitate rational metabolic engineering of these fungi as microbial cell factories.

Multispecies transcriptomes reveal core fruit development genes

During angiosperm evolution there have been repeated transitions from an ancestral dry fruit to a derived fleshy fruit, often with dramatic ecological and economic consequences. Following the transition to fleshy fruits, domestication may also dramatically alter the fruit phenotype via artificial selection. Although the morphologies of these fruits are well documented, relatively less is known about the molecular basis of these developmental and evolutionary shifts. We generated RNA-seq libraries from pericarp tissue of desert tobacco and both cultivated and wild tomato species at common developmental time points and combined this with corresponding, publicly available data from Arabidopsis and melon. With this broadly sampled dataset consisting of dry/fleshy fruits and wild/domesticated species, we applied novel bioinformatic methods to investigate conserved and divergent patterns of gene expression during fruit development and evolution. A small set of 121 orthologous "core" fruit development genes show a common pattern of expression across all five species. These include key players in developmental patterning such as orthologs of KNOLLE, PERIANTHIA, and ARGONAUTE7. GO term enrichment suggests that these genes function in basic cell division processes, cell wall biosynthesis, and developmental patterning. We furthermore uncovered a number of "accessory" genes with conserved expression patterns within but not among fruit types, and whose functional enrichment highlights the conspicuous differences between these phenotypic classes. We observe striking conservation of gene expression patterns despite large evolutionary distances, and dramatic phenotypic shifts, suggesting a conserved function for a small subset of core fruit development genes.

Genome-guided prediction of acid resistance mechanisms in acidophilic methanotrophs of phylogenetically deep-rooted Verrucomicrobia isolated from geothermal environments

Verrucomicrobia are a group of microorganisms that have been proposed to be deeply rooted in the Tree of Life. Some are methanotrophs that oxidize the potent greenhouse gas methane and are thus important in decreasing atmospheric concentrations of the gas, potentially ameliorating climate change. They are widespread in various environments including soil and fresh or marine waters. Recently, a clade of extremely acidophilic Verrucomicrobia, flourishing at pH < 3, were described from high-temperature geothermal ecosystems. This novel group could be of interest for studies about the emergence of life on Earth and to astrobiologists as homologs for possible extraterrestrial life. In this paper, we describe predicted mechanisms for survival of this clade at low pH and suggest its possible evolutionary trajectory from an inferred neutrophilic ancestor. Extreme acidophiles are defined as organisms that thrive in extremely low pH environments (≤ pH 3). Many are polyextremophiles facing high temperatures and high salt as well as low pH. They are important to study for both providing fundamental insights into biological mechanisms of survival and evolution in such extreme environments and for understanding their roles in biotechnological applications such as industrial mineral recovery (bioleaching) and mitigation of acid mine drainage. They are also, potentially, a rich source of novel genes and pathways for the genetic engineering of microbial strains. Acidophiles of the Verrucomicrobia phylum are unique as they are the only known aerobic methanotrophs that can grow optimally under acidic (pH 2–3) and moderately thermophilic conditions (50–60°C). Three moderately thermophilic genera, namely Methylacidiphilum, Methylacidimicrobium, and Ca. Methylacidithermus, have been described in geothermal environments. Most of the investigations of these organisms have focused on their methane oxidizing capabilities (methanotrophy) and use of lanthanides as a protein cofactor, with no extensive study that sheds light on the mechanisms that they use to flourish at extremely low pH. In this paper, we extend the phylogenetic description of this group of acidophiles using whole genome information and we identify several mechanisms, potentially involved in acid resistance, including “first line of defense” mechanisms that impede the entry of protons into the cell. These include the presence of membrane-associated hopanoids, multiple copies of the outer membrane protein (Slp), and inner membrane potassium channels (kup, kdp) that generate a reversed membrane potential repelling the intrusion of protons. Acidophilic Verrucomicrobia also display a wide array of proteins potentially involved in the “second line of defense” where protons that evaded the first line of defense and entered the cell are expelled or neutralized, such as the glutamate decarboxylation (gadAB) and phosphate-uptake systems. An exclusive N-type ATPase F₀-F₁ was identified only in acidophiles of Verrucomicrobia and is predicted to be a specific adaptation in these organisms. Phylogenetic analyses suggest that many predicted mechanisms are evolutionarily conserved and most likely entered the acidophilic lineage of Verrucomicrobia by vertical descent from a common ancestor. However, it is likely that some defense mechanisms such as gadA and kup entered the acidophilic Verrucomicrobia lineage by horizontal gene transfer.

A genome-wide comparative evolutionary analysis of zinc finger-BED transcription factor genes in land plants

Zinc finger (Zf)-BED proteins are a novel superfamily of transcription factors that controls numerous activities in plants including growth, development, and cellular responses to biotic and abiotic stresses. Despite their important roles in gene regulation, little is known about the specific functions of Zf-BEDs in land plants. The current study identified a total of 750 Zf-BED-encoding genes in 35 land plant species including mosses, bryophytes, lycophytes, gymnosperms, and angiosperms. The gene family size was somewhat proportional to genome size. All identified genes were categorized into 22 classes based on their specific domain architectures. Of these, class I (Zf-BED_DUF-domain_Dimer_Tnp_hAT) was the most common in the majority of the land plants. However, some classes were family-specific, while the others were species-specific, demonstrating diversity at different classification levels. In addition, several novel functional domains were also predicated including WRKY and nucleotide-binding site (NBS). Comparative genomics, transcriptomics, and proteomics provided insights into the evolutionary history, duplication, divergence, gene gain and loss, species relationship, expression profiling, and structural diversity of Zf-BEDs in land plants. The comprehensive study of Zf-BEDs in Gossypium sp., (cotton) also demonstrated a clear footprint of polyploidization. Overall, this comprehensive evolutionary study of Zf-BEDs in land plants highlighted significant diversity among plant species.

Transcription factor lineages in plant-pathogenic fungi, connecting diversity with fungal virulence

Plant-pathogenic fungi span diverse taxonomic lineages. Their host-infection strategies are often specialised and require the coordinated regulation of molecular virulence factors. Transcription factors (TFs) are fundamental regulators of gene expression, yet relatively few virulence-specific regulators are characterised in detail and their evolutionary trajectories are not well understood. Hence, this study compared the full range of TFs across taxonomically-diverse fungal proteomes and classified their lineages through an orthology analysis. The primary aims were to characterise differences in the range and profile of TF lineages broadly linked to plant-host association or pathogenic lifestyles, and to better characterise the evolutionary origin and trajectory of experimentally-validated virulence regulators. We observed significantly fewer TFs among obligate, host-associated pathogens, largely attributed to contractions in several Zn2Cys6 TF-orthogroup lineages. We also present novel insight into the key virulence-regulating TFs Ste12, Pf2 and EBR1, providing evidence for their ancestral origins, expansion and/or loss. Ultimately, the analysis presented here provides both primary evidence for TF evolution in fungal phytopathogenicity, as well as a practical phylogenetic resource to guide further detailed investigation on the regulation of virulence within key pathogen lineages.

SHOOT: phylogenetic gene search and ortholog inference

Determining the evolutionary relationships between genes is fundamental to comparative biological research. Here, we present SHOOT. SHOOT searches a user query sequence against a database of phylogenetic trees and returns a tree with the query sequence correctly placed within it. We show that SHOOT performs this analysis with comparable speed to a BLAST search. We demonstrate that SHOOT phylogenetic placements are as accurate as conventional tree inference, and it can identify orthologs with high accuracy. In summary, SHOOT is a fast and accurate tool for phylogenetic analyses of novel query sequences. It is available online at www.shoot.bio.

Genome wide study of cysteine rich receptor like proteins in Gossypium sp

Cysteine-rich receptor-like-kinases (CRKs), a transmembrane subfamily of receptor-like kinase, play crucial roles in plant adaptation. As such cotton is the major source of fiber for the textile industry, but environmental stresses are limiting its growth and production. Here, we have performed a deep computational analysis of CRKs in five Gossypium species, including G. arboreum (60 genes), G. raimondii (74 genes), G. herbaceum (65 genes), G. hirsutum (118 genes), and G. barbadense (120 genes). All identified CRKs were classified into 11 major classes and 43 subclasses with the finding of several novel CRK-associated domains including ALMT, FUSC_2, Cript, FYVE, and Pkinase. Of these, DUF26_DUF26_Pkinase_Tyr was common and had elevated expression under different biotic and abiotic stresses. Moreover, the 35 land plants comparison identified several new CRKs domain-architectures. Likewise, several SNPs and InDels were observed in CLCuD resistant G. hirsutum. The miRNA target side prediction and their expression profiling in different tissues predicted miR172 as a major CRK regulating miR. The expression profiling of CRKs identified multiple clusters with co-expression under certain stress conditions. The expression analysis under CLCuD highlighted the role of GhCRK057, GhCRK059, GhCRK058, and GhCRK081 in resistant accession. Overall, these results provided primary data for future potential functional analysis as well as a reference study for other agronomically important crops.

The genome sequence of the avian vampire fly (Philornis downsi), an invasive nest parasite of Darwin’s finches in Galápagos

The invasive avian vampire fly (Philornis downsi, Diptera: Muscidae) is considered one of the greatest threats to the endemic avifauna of the Galápagos Islands. The fly larvae parasitize nearly every passerine species, including Darwin’s finches. Most P. downsi research to date has focused on the effects of the fly on avian host fitness and mitigation methods. A lag in research related to the genetics of this invasion demonstrates, in part, the need to develop full-scale genomic resources with which to address further questions within this system. In this study, an adult female P. downsi was sequenced to generate a high-quality genome assembly. We examined various features of the genome (e.g., coding regions and noncoding transposable elements) and carried out comparative genomics analysis against other dipteran genomes. We identified lists of gene families that are significantly expanding or contracting in P. downsi that are related to insecticide resistance, detoxification, and counter defense against host immune responses. The P. downsi genome assembly provides an important resource for studying the molecular basis of successful invasion in the Galápagos and the dynamics of its population across multiple islands. The findings of significantly changing gene families associated with insecticide resistance and immune responses highlight the need for further investigations into the role of different gene families in aiding the fly’s successful invasion. Furthermore, this genomic resource provides a necessary tool to better inform future research studies and mitigation strategies aimed at minimizing the fly’s impact on Galápagos birds.

A Large-Scale Multiple Genome Comparison of Acidophilic Archaea (pH ≤ 5.0) Extends Our Understanding of Oxidative Stress Responses in Polyextreme Environments

Acidophilic archaea thrive in anaerobic and aerobic low pH environments (pH < 5) rich in dissolved heavy metals that exacerbate stress caused by the production of reactive oxygen species (ROS) such as hydrogen peroxide (H2O2), hydroxyl radical (OH) and superoxide (O2−). ROS react with lipids, proteins and nucleic acids causing oxidative stress and damage that can lead to cell death. Herein, genes and mechanisms potentially involved in ROS mitigation are predicted in over 200 genomes of acidophilic archaea with sequenced genomes. These organisms are often be subjected to simultaneous multiple stresses such as high temperature, high salinity, low pH and high heavy metal loads. Some of the topics addressed include: (1) the phylogenomic distribution of these genes and what this can tell us about the evolution of these mechanisms in acidophilic archaea; (2) key differences in genes and mechanisms used by acidophilic versus non-acidophilic archaea and between acidophilic archaea and acidophilic bacteria and (3) how comparative genomic analysis predicts novel genes or pathways involved in oxidative stress responses in archaea and likely horizontal gene transfer (HGT) events.

Lanza V, Rodríguez-Beltrán J, Galán JC, San Millán A, Cantón R, Coque TM. Evolutionary Pathways and Trajectories in Antibiotic Resistance

Evolution is the hallmark of life. Descriptions of the evolution of microorganisms have provided a wealth of information, but knowledge regarding "what happened" has precluded a deeper understanding of "how" evolution has proceeded, as in the case of antimicrobial resistance. The difficulty in answering the "how" question lies in the multihierarchical dimensions of evolutionary processes, nested in complex networks, encompassing all units of selection, from genes to communities and ecosystems. At the simplest ontological level (as resistance genes), evolution proceeds by random (mutation and drift) and directional (natural selection) processes; however, sequential pathways of adaptive variation can occasionally be observed, and under fixed circumstances (particular fitness landscapes), evolution is predictable. At the highest level (such as that of plasmids, clones, species, microbiotas), the systems' degrees of freedom increase dramatically, related to the variable dispersal, fragmentation, relatedness, or coalescence of bacterial populations, depending on heterogeneous and changing niches and selective gradients in complex environments. Evolutionary trajectories of antibiotic resistance find their way in these changing landscapes subjected to random variations, becoming highly entropic and therefore unpredictable. However, experimental, phylogenetic, and ecogenetic analyses reveal preferential frequented paths (highways) where antibiotic resistance flows and propagates, allowing some understanding of evolutionary dynamics, modeling and designing interventions. Studies on antibiotic resistance have an applied aspect in improving individual health, One Health, and Global Health, as well as an academic value for understanding evolution. Most importantly, they have a heuristic significance as a model to reduce the negative influence of anthropogenic effects on the environment.

Phenotypic Characterization and Comparative Genomics of the Melanin-Producing Yeast Exophiala lecanii-corni Reveals a Distinct Stress Tolerance Profile and Reduced Ribosomal Genetic Content

The black yeast Exophiala lecanii-corni of the order Chaetothyriales is notable for its ability to produce abundant quantities of DHN-melanin. While many other Exophiala species are frequent causal agents of human infection, E. lecanii-corni CBS 102400 lacks the thermotolerance requirements that enable pathogenicity, making it appealing for use in targeted functional studies and biotechnological applications. Here, we report the stress tolerance characteristics of E. lecanii-corni, with an emphasis on the influence of melanin on its resistance to various forms of stress. We find that E. lecanii-corni has a distinct stress tolerance profile that includes variation in resistance to temperature, osmotic, and oxidative stress relative to the extremophilic and pathogenic black yeast Exophiala dermatitidis. Notably, the presence of melanin substantially impacts stress resistance in E. lecanii-corni, while this was not found to be the case in E. dermatitidis. The cellular context, therefore, influences the role of melanin in stress protection. In addition, we present a detailed analysis of the E. lecanii-corni genome, revealing key differences in functional genetic content relative to other ascomycetous species, including a significant decrease in abundance of genes encoding ribosomal proteins. In all, this study provides insight into how genetics and physiology may underlie stress tolerance and enhances understanding of the genetic diversity of black yeasts.

Genome sequence, transcriptome, and annotation of rodent malaria parasite Plasmodium yoelii nigeriensis N67

Background

Rodent malaria parasites are important models for studying host-malaria parasite interactions such as host immune response, mechanisms of parasite evasion of host killing, and vaccine development. One of the rodent malaria parasites is Plasmodium yoelii, and multiple P. yoelii strains or subspecies that cause different disease phenotypes have been widely employed in various studies. The genomes and transcriptomes of several P. yoelii strains have been analyzed and annotated, including the lethal strains of P. y. yoelii YM (or 17XL) and non-lethal strains of P. y. yoelii 17XNL/17X. Genomic DNA sequences and cDNA reads from another subspecies P. y. nigeriensis N67 have been reported for studies of genetic polymorphisms and parasite response to drugs, but its genome has not been assembled and annotated.

Results

We performed genome sequencing of the N67 parasite using the PacBio long-read sequencing technology, de novo assembled its genome and transcriptome, and predicted 5383 genes with high overall annotation quality. Comparison of the annotated genome of the N67 parasite with those of YM and 17X parasites revealed a set of genes with N67-specific orthology, expansion of gene families, particularly the homologs of the Plasmodium chabaudi erythrocyte membrane antigen, large numbers of SNPs and indels, and proteins predicted to interact with host immune responses based on their functional domains.

Conclusions

The genomes of N67 and 17X parasites are highly diverse, having approximately one polymorphic site per 50 base pairs of DNA. The annotated N67 genome and transcriptome provide searchable databases for fast retrieval of genes and proteins, which will greatly facilitate our efforts in studying the parasite biology and gene function and in developing effective control measures against malaria.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-07555-9.

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.

Evolutionary analyses of the major variant surface antigen-encoding genes reveal population structure of Plasmodium falciparum within and between continents

Malaria remains a major public health problem in many countries. Unlike influenza and HIV, where diversity in immunodominant surface antigens is understood geographically to inform disease surveillance, relatively little is known about the global population structure of PfEMP1, the major variant surface antigen of the malaria parasite Plasmodium falciparum. The complexity of the var multigene family that encodes PfEMP1 and that diversifies by recombination, has so far precluded its use in malaria surveillance. Recent studies have demonstrated that cost-effective deep sequencing of the region of var genes encoding the PfEMP1 DBLα domain and subsequent classification of within host sequences at 96% identity to define unique DBLα types, can reveal structure and strain dynamics within countries. However, to date there has not been a comprehensive comparison of these DBLα types between countries. By leveraging a bioinformatic approach (jumping hidden Markov model) designed specifically for the analysis of recombination within var genes and applying it to a dataset of DBLα types from 10 countries, we are able to describe population structure of DBLα types at the global scale. The sensitivity of the approach allows for the comparison of the global dataset to ape samples of Plasmodium Laverania species. Our analyses show that the evolution of the parasite population emerging out of Africa underlies current patterns of DBLα type diversity. Most importantly, we can distinguish geographic population structure within Africa between Gabon and Ghana in West Africa and Uganda in East Africa. Our evolutionary findings have translational implications in the context of globalization. Firstly, DBLα type diversity can provide a simple diagnostic framework for geographic surveillance of the rapidly evolving transmission dynamics of P. falciparum. It can also inform efforts to understand the presence or absence of global, regional and local population immunity to major surface antigen variants. Additionally, we identify a number of highly conserved DBLα types that are present globally that may be of biological significance and warrant further characterization.

De novo genome assembly of the tobacco hornworm moth (Manduca sexta)

The tobacco hornworm, Manduca sexta, is a lepidopteran insect that is used extensively as a model system for studying insect biology, development, neuroscience, and immunity. However, current studies rely on the highly fragmented reference genome Msex_1.0, which was created using now-outdated technologies and is hindered by a variety of deficiencies and inaccuracies. We present a new reference genome for M. sexta, JHU_Msex_v1.0, applying a combination of modern technologies in a de novo assembly to increase continuity, accuracy, and completeness. The assembly is 470 Mb and is ∼20× more continuous than the original assembly, with scaffold N50 > 14 Mb. We annotated the assembly by lifting over existing annotations and supplementing with additional supporting RNA-based data for a total of 25,256 genes. The new reference assembly is accessible in annotated form for public use. We demonstrate that improved continuity of the M. sexta genome improves resequencing studies and benefits future research on M. sexta as a model organism.

Genetic Structure, Function, and Evolution of Capsule Biosynthesis Loci in Vibrio parahaemolyticus

Capsule-forming extracellular polysaccharides are crucial for bacterial host colonization, invasion, immune evasion, and ultimately pathogenicity. Due to warming ocean waters and human encroachment of coastal ecosystems, Vibrio parahaemolyticus has emerged as a globally important foodborne enteropathogen implicated in acute gastroenteritis, wound infections, and septic shock. Conventionally, the antigenic properties of lipopolysaccharide (LPS, O antigen) and capsular polysaccharide (CPS, K antigen) have provided a basis for serotyping V. parahaemolyticus, whereas disclosure of genetic elements encoding 13 O-serogroups have allowed molecular serotyping methods to be developed. However, the genetic structure of CPS loci for 71 K-serogroups has remained unidentified, limiting progress in understanding its roles in V. parahaemolyticus pathophysiology. In this study, we identified and characterized the genetic structure and their evolutionary relationship of CPS loci of 40 K-serogroups through whole genome sequencing of 443 V. parahaemolyticus strains. We found a distinct pattern of CPS gene cluster across different K-serogroups and expanded its new 3'-border by identifying glpX as a key gene conserved across all K-serogroups. A total of 217 genes involved in CPS biosynthesis were annotated. Functional contents and genetic structure of the 40 K-serogroups were analyzed. Based on inferences from species trees and gene trees, we proposed an evolution model of the CPS gene clusters of 40 K-serogroups. Horizontal gene transfer by recombination from other Vibrio species, gene duplication is likely to play instrumental roles in the evolution of CPS in V. parahaemolyticus. This is the first time, to the best of our knowledge, that a large scale of CPS gene clusters of different K-serogroups in V. parahaemolyticus have been identified and characterized in evolutionary contexts. This work should help advance understanding on the variation of CPS in V. parahaemolyticus and provide a framework for developing diagnostically relevant serotyping methods.

Genomic profiling of bacterial and fungal communities and their predictive functionality during pulque fermentation by whole-genome shotgun sequencing

Pulque is a culturally important 4,000-year-old traditional Mexican fermented drink. Pulque is produced by adding fresh aguamiel (agave sap) to mature pulque, resulting in a mixture of microbial communities and chemical compositions. We performed shotgun metagenomic sequencing of five stages of pulque fermentation to characterize organismal and functional diversity. We identified 6 genera (Acinetobacter, Lactobacillus, Lactococcus, Leuconostoc, Saccharomyces and Zymomonas) and 10 species (Acinetobacter boissieri, Acinetobacter nectaris, Lactobacillus sanfranciscensis, Lactococcus lactis, Lactococcus piscium, Lactococcus plantarum, Leuconostoc citreum, Leuconostoc gelidum, Zymomonas mobilis and Saccharomyces cerevisiae) that were present ≥ 1% in at least one stage of pulque fermentation. The abundance of genera and species changed during fermentation and was associated with a decrease in sucrose and increases in ethanol and lactic acid, suggesting that resource competition shapes organismal diversity. We also predicted functional profiles, based on organismal gene content, for each fermentation stage and identified an abundance of genes associated with the biosynthesis of folate, an essential B-vitamin. Additionally, we investigated the evolutionary relationships of S. cerevisiae and Z. mobilis, two of the major microbial species found in pulque. For S. cerevisiae, we used a metagenomics assembly approach to identify S. cerevisiae scaffolds from pulque, and performed phylogenetic analysis of these sequences along with a collection of 158 S. cerevisiae strains. This analysis suggests that S. cerevisiae from pulque is most closely related to Asian strains isolated from sake and bioethanol. Lastly, we isolated and sequenced the whole-genomes of three strains of Z. mobilis from pulque and compared their relationship to seven previously sequenced isolates. Our results suggest pulque strains may represent a distinct lineage of Z. mobilis.

Gene gain and loss across the metazoan tree of life

Although recent research has revealed high genomic complexity in the earliest-splitting animals and their ancestors, the macroevolutionary trends orchestrating gene repertoire evolution throughout the animal phyla remain poorly understood. We used a phylogenomic approach to interrogate genome evolution across all animal phyla. Our analysis uncovered a bimodal distribution of recruitment of orthologous genes, with most genes gained very 'early' (that is, at deep nodes) or very 'late', representing lineage-specific acquisitions. The emergence of animals was characterized by high values of gene birth and duplications. Deuterostomes, ecdysozoans and Xenacoelomorpha were characterized by no gene gain but rampant differential gene loss. Genes considered as animal hallmarks, such as Notch/Delta, were convergently duplicated in all phyla and at different evolutionary depths. Genes duplicated in all nodes from Metazoa to phylum-specific levels were enriched in functions related to the neural system, suggesting that this system has been continuously and independently reshaped throughout evolution across animals. Our results indicate that animal genomes evolved by unparalleled gene duplication followed by differential gene loss, and provide an atlas of gene repertoire evolution throughout the animal tree of life to navigate how, when and how often each gene in each genome was gained, duplicated or lost.

Genome Insights of the Plant-Growth Promoting Bacterium Cronobacter muytjensii JZ38 With Volatile-Mediated Antagonistic Activity Against Phytophthora infestans

Salinity stress is a major challenge to agricultural productivity and global food security in light of a dramatic increase of human population and climate change. Plant growth promoting bacteria can be used as an additional solution to traditional crop breeding and genetic engineering. In the present work, the induction of plant salt tolerance by the desert plant endophyte Cronobacter sp. JZ38 was examined on the model plant Arabidopsis thaliana using different inoculation methods. JZ38 promoted plant growth under salinity stress via contact and emission of volatile compounds. Based on the 16S rRNA and whole genome phylogenetic analysis, fatty acid analysis and phenotypic identification, JZ38 was identified as Cronobacter muytjensii and clearly separated and differentiated from the pathogenic C. sakazakii. Full genome sequencing showed that JZ38 is composed of one chromosome and two plasmids. Bioinformatic analysis and bioassays revealed that JZ38 can grow under a range of abiotic stresses. JZ38 interaction with plants is correlated with an extensive set of genes involved in chemotaxis and motility. The presence of genes for plant nutrient acquisition and phytohormone production could explain the ability of JZ38 to colonize plants and sustain plant growth under stress conditions. Gas chromatography-mass spectrometry analysis of volatiles produced by JZ38 revealed the emission of indole and different sulfur volatile compounds that may play a role in contactless plant growth promotion and antagonistic activity against pathogenic microbes. Indeed, JZ38 was able to inhibit the growth of two strains of the phytopathogenic oomycete Phytophthora infestans via volatile emission. Genetic, transcriptomic and metabolomics analyses, combined with more in vitro assays will provide a better understanding the highlighted genes' involvement in JZ38's functional potential and its interaction with plants. Nevertheless, these results provide insight into the bioactivity of C. muytjensii JZ38 as a multi-stress tolerance promoting bacterium with a potential use in agriculture.

OrthoFinder: phylogenetic orthology inference for comparative genomics

Here, we present a major advance of the OrthoFinder method. This extends OrthoFinder's high accuracy orthogroup inference to provide phylogenetic inference of orthologs, rooted gene trees, gene duplication events, the rooted species tree, and comparative genomics statistics. Each output is benchmarked on appropriate real or simulated datasets, and where comparable methods exist, OrthoFinder is equivalent to or outperforms these methods. Furthermore, OrthoFinder is the most accurate ortholog inference method on the Quest for Orthologs benchmark test. Finally, OrthoFinder's comprehensive phylogenetic analysis is achieved with equivalent speed and scalability to the fastest, score-based heuristic methods. OrthoFinder is available at https://github.com/davidemms/OrthoFinder.

Primer Design for an Accurate View of Picocyanobacterial Community Structure by Using High-Throughput Sequencing

High-throughput sequencing (HTS) of the 16S rRNA gene has been used successfully to describe the structure and dynamics of microbial communities. Picocyanobacteria are important members of bacterioplankton communities, and, so far, they have predominantly been targeted using universal bacterial primers, providing a limited resolution of the picocyanobacterial community structure and dynamics. To increase such resolution, the study of a particular target group is best approached with the use of specific primers. Here, we aimed to design and evaluate specific primers for aquatic picocyanobacterial genera to be used with high-throughput sequencing. Since the various regions of the 16S rRNA gene have different degrees of conservation in different bacterial groups, we therefore first determined which hypervariable region of the 16S rRNA gene provides the highest taxonomic and phylogenetic resolution for the genera Synechococcus, Prochlorococcus, and Cyanobium An in silico analysis showed that the V5, V6, and V7 hypervariable regions appear to be the most informative for this group. We then designed primers flanking these hypervariable regions and tested them in natural marine and freshwater communities. We successfully detected that most (97%) of the obtained reads could be assigned to picocyanobacterial genera. We defined operational taxonomic units as exact sequence variants (zero-radius operational taxonomic units [zOTUs]), which allowed us to detect higher genetic diversity and infer ecologically relevant information about picocyanobacterial community composition and dynamics in different aquatic systems. Our results open the door to future studies investigating picocyanobacterial diversity in aquatic systems.IMPORTANCE The molecular diversity of the aquatic picocyanobacterial community cannot be accurately described using only the available universal 16S rRNA gene primers that target the whole bacterial and archaeal community. We show that the hypervariable regions V5, V6, and V7 of the 16S rRNA gene are better suited to study the diversity, community structure, and dynamics of picocyanobacterial communities at a fine scale using Illumina MiSeq sequencing. Due to its variability, it allows reconstructing phylogenies featuring topologies comparable to those generated when using the complete 16S rRNA gene sequence. Further, we successfully designed a new set of primers flanking the V5 to V7 region whose specificity for picocyanobacterial genera was tested in silico and validated in several freshwater and marine aquatic communities. This work represents a step forward for understanding the diversity and ecology of aquatic picocyanobacteria and sets the path for future studies on picocyanobacterial diversity.

A Malvaceae mystery: A mallow maelstrom of genome multiplications and maybe misleading methods?

Previous research suggests that Gossypium has undergone a 5- to 6-fold multiplication following its divergence from Theobroma. However, the number of events, or where they occurred in the Malvaceae phylogeny remains unknown. We analyzed transcriptomic and genomic data from representatives of eight of the nine Malvaceae subfamilies. Phylogenetic analysis of nuclear data placed Dombeya (Dombeyoideae) as sister to the rest of Malvadendrina clade, but the plastid DNA tree strongly supported Durio (Helicteroideae) in this position. Intraspecific Ks plots indicated that all sampled taxa, except Theobroma (Byttnerioideae), Corchorus (Grewioideae), and Dombeya (Dombeyoideae), have experienced whole genome multiplications (WGMs). Quartet analysis suggested WGMs were shared by Malvoideae-Bombacoideae and Sterculioideae-Tilioideae, but did not resolve whether these are shared with each other or Helicteroideae (Durio). Gene tree reconciliation and Bayesian concordance analysis suggested a complex history. Alternative hypotheses are suggested, each involving two independent autotetraploid and one allopolyploid event. They differ in that one entails an allopolyploid origin for the Durio lineage, whereas the other invokes an allopolyploid origin for Malvoideae-Bombacoideae. We highlight the need for more genomic information in the Malvaceae and improved methods to resolve complex evolutionary histories that may include allopolyploidy, incomplete lineage sorting, and variable rates of gene and genome evolution.

In silico exploration of Red Sea Bacillus genomes for natural product biosynthetic gene clusters

BACKGROUND

The increasing spectrum of multidrug-resistant bacteria is a major global public health concern, necessitating discovery of novel antimicrobial agents. Here, members of the genus Bacillus are investigated as a potentially attractive source of novel antibiotics due to their broad spectrum of antimicrobial activities. We specifically focus on a computational analysis of the distinctive biosynthetic potential of Bacillus paralicheniformis strains isolated from the Red Sea, an ecosystem exposed to adverse, highly saline and hot conditions.

RESULTS

We report the complete circular and annotated genomes of two Red Sea strains, B. paralicheniformis Bac48 isolated from mangrove mud and B. paralicheniformis Bac84 isolated from microbial mat collected from Rabigh Harbor Lagoon in Saudi Arabia. Comparing the genomes of B. paralicheniformis Bac48 and B. paralicheniformis Bac84 with nine publicly available complete genomes of B. licheniformis and three genomes of B. paralicheniformis, revealed that all of the B. paralicheniformis strains in this study are more enriched in nonribosomal peptides (NRPs). We further report the first computationally identified trans-acyltransferase (trans-AT) nonribosomal peptide synthetase/polyketide synthase (PKS/ NRPS) cluster in strains of this species.

CONCLUSIONS

B. paralicheniformis species have more genes associated with biosynthesis of antimicrobial bioactive compounds than other previously characterized species of B. licheniformis, which suggests that these species are better potential sources for novel antibiotics. Moreover, the genome of the Red Sea strain B. paralicheniformis Bac48 is more enriched in modular PKS genes compared to B. licheniformis strains and other B. paralicheniformis strains. This may be linked to adaptations that strains surviving in the Red Sea underwent to survive in the relatively hot and saline ecosystems.

Functional Annotations of Paralogs: A Blessing and a Curse

Gene duplication followed by mutation is a classic mechanism of neofunctionalization, producing gene families with functional diversity. In some cases, a single point mutation is sufficient to change the substrate specificity and/or the chemistry performed by an enzyme, making it difficult to accurately separate enzymes with identical functions from homologs with different functions. Because sequence similarity is often used as a basis for assigning functional annotations to genes, non-isofunctional gene families pose a great challenge for genome annotation pipelines. Here we describe how integrating evolutionary and functional information such as genome context, phylogeny, metabolic reconstruction and signature motifs may be required to correctly annotate multifunctional families. These integrative analyses can also lead to the discovery of novel gene functions, as hints from specific subgroups can guide the functional characterization of other members of the family. We demonstrate how careful manual curation processes using comparative genomics can disambiguate subgroups within large multifunctional families and discover their functions. We present the COG0720 protein family as a case study. We also discuss strategies to automate this process to improve the accuracy of genome functional annotation pipelines.

An Ancestral Role for CONSTITUTIVE TRIPLE RESPONSE1 Proteins in Both Ethylene and Abscisic Acid Signaling

Land plants have evolved adaptive regulatory mechanisms enabling the survival of environmental stresses associated with terrestrial life. Here, we focus on the evolution of the regulatory CONSTITUTIVE TRIPLE RESPONSE1 (CTR1) component of the ethylene signaling pathway that modulates stress-related changes in plant growth and development. First, we compare CTR1-like proteins from a bryophyte, Physcomitrella patens (representative of early divergent land plants), with those of more recently diverged lycophyte and angiosperm species (including Arabidopsis [Arabidopsis thaliana]) and identify a monophyletic CTR1 family. The fully sequenced P. patens genome encodes only a single member of this family (PpCTR1L). Next, we compare the functions of PpCTR1L with that of related angiosperm proteins. We show that, like angiosperm CTR1 proteins (e.g. AtCTR1 of Arabidopsis), PpCTR1L modulates downstream ethylene signaling via direct interaction with ethylene receptors. These functions, therefore, likely predate the divergence of the bryophytes from the land-plant lineage. However, we also show that PpCTR1L unexpectedly has dual functions and additionally modulates abscisic acid (ABA) signaling. In contrast, while AtCTR1 lacks detectable ABA signaling functions, Arabidopsis has during evolution acquired another homolog that is functionally distinct from AtCTR1. In conclusion, the roles of CTR1-related proteins appear to have functionally diversified during land-plant evolution, and angiosperm CTR1-related proteins appear to have lost an ancestral ABA signaling function. Our study provides new insights into how molecular events such as gene duplication and functional differentiation may have contributed to the adaptive evolution of regulatory mechanisms in plants.

OrthoFinder: solving fundamental biases in whole genome comparisons dramatically improves orthogroup inference accuracy

Identifying homology relationships between sequences is fundamental to biological research. Here we provide a novel orthogroup inference algorithm called OrthoFinder that solves a previously undetected gene length bias in orthogroup inference, resulting in significant improvements in accuracy. Using real benchmark datasets we demonstrate that OrthoFinder is more accurate than other orthogroup inference methods by between 8 % and 33 %. Furthermore, we demonstrate the utility of OrthoFinder by providing a complete classification of transcription factor gene families in plants revealing 6.9 million previously unobserved relationships.

OrthoFinder: solving fundamental biases in whole genome comparisons dramatically improves orthogroup inference accuracy

Identifying homology relationships between sequences is fundamental to biological research. Here we provide a novel orthogroup inference algorithm called OrthoFinder that solves a previously undetected gene length bias in orthogroup inference, resulting in significant improvements in accuracy. Using real benchmark datasets we demonstrate that OrthoFinder is more accurate than other orthogroup inference methods by between 8 % and 33 %. Furthermore, we demonstrate the utility of OrthoFinder by providing a complete classification of transcription factor gene families in plants revealing 6.9 million previously unobserved relationships.