PhySortR: a fast, flexible tool for sorting phylogenetic trees in R

Genomes of nitrogen-fixing eukaryotes reveal a non-canonical model of organellogenesis

Endosymbiont gene transfer and import of host-encoded proteins are considered hallmarks of organelles necessary for stable integration of two cells. However, newer endosymbiotic models have challenged the origin and timing of such genetic integration during organellogenesis. Epithemia diatoms contain diazoplasts, closely related to recently-described nitrogen-fixing organelles, that are also stably integrated and co-speciating with their host algae. We report genomic analyses of two species, freshwater E.clementina and marine E.pelagica, which are highly divergent but share a common endosymbiotic origin. We found minimal evidence of genetic integration: nonfunctional diazoplast-to-nuclear DNA transfers in the E.clementina genome and 6 host-encoded proteins of unknown function in the E.clementina diazoplast proteome, far fewer than in other recently-acquired organelles. Epithemia diazoplasts are a valuable counterpoint to existing organellogenesis models, demonstrating that endosymbionts can be stably integrated and inherited absent significant genetic integration. The minimal genetic integration makes diazoplasts valuable blueprints for bioengineering endosymbiotic compartments de novo.

A k-mer-Based Approach for Phylogenetic Classification of Taxa in Environmental Genomic Data

In the age of genome sequencing, whole-genome data is readily and frequently generated, leading to a wealth of new information that can be used to advance various fields of research. New approaches, such as alignment-free phylogenetic methods that utilize k-mer-based distance scoring, are becoming increasingly popular given their ability to rapidly generate phylogenetic information from whole-genome data. However, these methods have not yet been tested using environmental data, which often tends to be highly fragmented and incomplete. Here, we compare the results of one alignment-free approach (which utilizes the D2 statistic) to traditional multi-gene maximum likelihood trees in 3 algal groups that have high-quality genome data available. In addition, we simulate lower-quality, fragmented genome data using these algae to test method robustness to genome quality and completeness. Finally, we apply the alignment-free approach to environmental metagenome assembled genome data of unclassified Saccharibacteria and Trebouxiophyte algae, and single-cell amplified data from uncultured marine stramenopiles to demonstrate its utility with real datasets. We find that in all instances, the alignment-free method produces phylogenies that are comparable, and often more informative, than those created using the traditional multi-gene approach. The k-mer-based method performs well even when there are significant missing data that include marker genes traditionally used for tree reconstruction. Our results demonstrate the value of alignment-free approaches for classifying novel, often cryptic or rare, species, that may not be culturable or are difficult to access using single-cell methods, but fill important gaps in the tree of life.

Multiple waves of viral invasions in Symbiodiniaceae algal genomes

Dinoflagellates from the family Symbiodiniaceae are phototrophic marine protists that engage in symbiosis with diverse hosts. Their large and distinct genomes are characterized by pervasive gene duplication and large-scale retroposition events. However, little is known about the role and scale of horizontal gene transfer (HGT) in the evolution of this algal family. In other dinoflagellates, high levels of HGTs have been observed, linked to major genomic transitions, such as the appearance of a viral-acquired nucleoprotein that originated via HGT from a large DNA algal virus. Previous work showed that Symbiodiniaceae from different hosts are actively infected by viral groups, such as giant DNA viruses and ssRNA viruses, that may play an important role in coral health. Latent viral infections may also occur, whereby viruses could persist in the cytoplasm or integrate into the host genome as a provirus. This hypothesis received experimental support; however, the cellular localization of putative latent viruses and their taxonomic affiliation are still unknown. In addition, despite the finding of viral sequences in some genomes of Symbiodiniaceae, viral origin, taxonomic breadth, and metabolic potential have not been explored. To address these questions, we searched for putative viral-derived proteins in thirteen Symbiodiniaceae genomes. We found fifty-nine candidate viral-derived HGTs that gave rise to twelve phylogenies across ten genomes. We also describe the taxonomic affiliation of these virus-related sequences, their structure, and their genomic context. These results lead us to propose a model to explain the origin and fate of Symbiodiniaceae viral acquisitions.

Improved Cladocopium goreaui Genome Assembly Reveals Features of a Facultative Coral Symbiont and the Complex Evolutionary History of Dinoflagellate Genes

Dinoflagellates of the family Symbiodiniaceae are crucial photosymbionts in corals and other marine organisms. Of these, Cladocopium goreaui is one of the most dominant symbiont species in the Indo-Pacific. Here, we present an improved genome assembly of C. goreaui combining new long-read sequence data with previously generated short-read data. Incorporating new full-length transcripts to guide gene prediction, the C. goreaui genome (1.2 Gb) exhibits a high extent of completeness (82.4% based on BUSCO protein recovery) and better resolution of repetitive sequence regions; 45,322 gene models were predicted, and 327 putative, topologically associated domains of the chromosomes were identified. Comparison with other Symbiodiniaceae genomes revealed a prevalence of repeats and duplicated genes in C. goreaui, and lineage-specific genes indicating functional innovation. Incorporating 2,841,408 protein sequences from 96 taxonomically diverse eukaryotes and representative prokaryotes in a phylogenomic approach, we assessed the evolutionary history of C. goreaui genes. Of the 5246 phylogenetic trees inferred from homologous protein sets containing two or more phyla, 35–36% have putatively originated via horizontal gene transfer (HGT), predominantly (19–23%) via an ancestral Archaeplastida lineage implicated in the endosymbiotic origin of plastids: 10–11% are of green algal origin, including genes encoding photosynthetic functions. Our results demonstrate the utility of long-read sequence data in resolving structural features of a dinoflagellate genome, and highlight how genetic transfer has shaped genome evolution of a facultative symbiont, and more broadly of dinoflagellates.

Insights into agar and secondary metabolite pathways from the genome of the red alga Gracilaria domingensis (Rhodophyta, Gracilariales)

Gracilariales is a clade of florideophycean red macroalgae known for being the main source of agar. We present a de novo genome assembly and annotation of Gracilaria domingensis, an agarophyte alga with flattened thallus widely distributed along Central and South American Atlantic intertidal zones. In addition to structural analysis, an organizational comparison was done with other Rhodophyta genomes. The nuclear genome has 78 Mbp, with 11,437 predicted coding genes, 4,075 of which did not have hits in sequence databases. We also predicted 1,567 noncoding RNAs, distributed in 14 classes. The plastid and mitochondrion genome structures were also obtained. Genes related to agar synthesis were identified. Genes for type II galactose sulfurylases could not be found. Genes related to ascorbate synthesis were found. These results suggest an intricate connection of cell wall polysaccharide synthesis and the redox systems through the use of L-galactose in Rhodophyta. The genome of G. domingensis should be valuable to phycological and aquacultural research, as it is the first tropical and Western Atlantic red macroalgal genome to be sequenced.

The evolution of huge Y chromosomes in Coccinia grandis and its sister, Coccinia schimperi

Microscopically dimorphic sex chromosomes in plants are rare, reducing our ability to study them. One difficulty has been the paucity of cultivatable species pairs for cytogenetic, genomic and experimental work. Here, we study the newly recognized sisters Coccinia grandis and Coccinia schimperi, both with large Y chromosomes as we here show for Co. schimperi. We built genetic maps for male and female Co. grandis using a full-sibling family, inferred gene sex-linkage, and, with Co. schimperi transcriptome data, tested whether X- and Y-alleles group by species or by sex. Most sex-linked genes for which we could include outgroups grouped the X- and Y-alleles by species, but some 10% instead grouped the two species' X-alleles. There was no relationship between XY synonymous-site divergences in these genes and gene position on the non-recombining part of the X, suggesting recombination arrest shortly before or after species divergence, here dated to about 3.6 Ma. Coccinia grandis and Co. schimperi are the species pair with the most heteromorphic sex chromosomes in vascular plants (the condition in their sister remains unknown), and future work could use them to study mechanisms of Y chromosome enlargement and parallel degeneration, or to test Haldane's rule about lower hybrid fitness in the heterogametic sex. This article is part of the theme issue 'Sex determination and sex chromosome evolution in land plants'.

Colocality to Cofunctionality: Eukaryotic Gene Neighborhoods as a Resource for Function Discovery

Diverging from the classic paradigm of random gene order in eukaryotes, gene proximity can be leveraged to systematically identify functionally related gene neighborhoods in eukaryotes, utilizing techniques pioneered in bacteria. Current methods of identifying gene neighborhoods typically rely on sequence similarity to characterized gene products. However, this approach is not robust for nonmodel organisms like algae, which are evolutionarily distant from well-characterized model organisms. Here, we utilize a comparative genomic approach to identify evolutionarily conserved proximal orthologous gene pairs conserved across at least two taxonomic classes of green algae. A total of 317 gene neighborhoods were identified. In some cases, gene proximity appears to have been conserved since before the streptophyte–chlorophyte split, 1,000 Ma. Using functional inferences derived from reconstructed evolutionary relationships, we identified several novel functional clusters. A putative mycosporine-like amino acid, “sunscreen,” neighborhood contains genes similar to either vertebrate or cyanobacterial pathways, suggesting a novel mosaic biosynthetic pathway in green algae. One of two putative arsenic-detoxification neighborhoods includes an organoarsenical transporter (ArsJ), a glyceraldehyde 3-phosphate dehydrogenase-like gene, homologs of which are involved in arsenic detoxification in bacteria, and a novel algal-specific phosphoglycerate kinase-like gene. Mutants of the ArsJ-like transporter and phosphoglycerate kinase-like genes in Chlamydomonas reinhardtii were found to be sensitive to arsenate, providing experimental support for the role of these identified neighbors in resistance to arsenate. Potential evolutionary origins of neighborhoods are discussed, and updated annotations for formerly poorly annotated genes are presented, highlighting the potential of this strategy for functional annotation.

Amoeba Genome Reveals Dominant Host Contribution to Plastid Endosymbiosis

Eukaryotic photosynthetic organelles, plastids, are the powerhouses of many aquatic and terrestrial ecosystems. The canonical plastid in algae and plants originated >1 Ga and therefore offers limited insights into the initial stages of organelle evolution. To address this issue, we focus here on the photosynthetic amoeba Paulinella micropora strain KR01 (hereafter, KR01) that underwent a more recent (∼124 Ma) primary endosymbiosis, resulting in a photosynthetic organelle termed the chromatophore. Analysis of genomic and transcriptomic data resulted in a high-quality draft assembly of size 707 Mb and 32,361 predicted gene models. A total of 291 chromatophore-targeted proteins were predicted in silico, 208 of which comprise the ancestral organelle proteome in photosynthetic Paulinella species with functions, among others, in nucleotide metabolism and oxidative stress response. Gene coexpression analysis identified networks containing known high light stress response genes as well as a variety of genes of unknown function (“dark” genes). We characterized diurnally rhythmic genes in this species and found that over 49% are dark. It was recently hypothesized that large double-stranded DNA viruses may have driven gene transfer to the nucleus in Paulinella and facilitated endosymbiosis. Our analyses do not support this idea, but rather suggest that these viruses in the KR01 and closely related P. micropora MYN1 genomes resulted from a more recent invasion.

Inferring Phylogenomic Relationship of Microbes Using Scalable Alignment-Free Methods

Inferring phylogenetic relationships among hundreds or thousands of microbial genomes is an increasingly common task. The conventional phylogenetic approach adopts multiple sequence alignment to compare gene-by-gene, concatenated multigene or whole-genome sequences, from which a phylogenetic tree would be inferred. These alignments follow the implicit assumption of full-length contiguity among homologous sequences. However, common events in microbial genome evolution (e.g., structural rearrangements and genetic recombination) violate this assumption. Moreover, aligning hundreds or thousands of sequences is computationally intensive and not scalable to the rate at which genome data are generated. Therefore, alignment-free methods present an attractive alternative strategy. Here we describe a scalable alignment-free strategy to infer phylogenetic relationships using complete genome sequences of bacteria and archaea, based on short, subsequences of length k (k-mers). We describe how this strategy can be extended to infer evolutionary relationships beyond a tree-like structure, to better capture both vertical and lateral signals of microbial evolution.

Genome Analyses of the Microalga Picochlorum Provide Insights into the Evolution of Thermotolerance in the Green Lineage

While the molecular events involved in cell responses to heat stress have been extensively studied, our understanding of the genetic basis of basal thermotolerance, and particularly its evolution within the green lineage, remains limited. Here, we present the 13.3-Mb haploid genome and transcriptomes of a halotolerant and thermotolerant unicellular green alga, Picochlorum costavermella (Trebouxiophyceae) to investigate the evolution of the genomic basis of thermotolerance. Differential gene expression at high and standard temperatures revealed that more of the gene families containing up-regulated genes at high temperature were recently evolved, and less originated at the ancestor of green plants. Inversely, there was an excess of ancient gene families containing transcriptionally repressed genes. Interestingly, there is a striking overlap between the thermotolerance and halotolerance transcriptional rewiring, as more than one-third of the gene families up-regulated at 35 °C were also up-regulated under variable salt concentrations in Picochlorum SE3. Moreover, phylogenetic analysis of the 9,304 protein coding genes revealed 26 genes of horizontally transferred origin in P. costavermella, of which five were differentially expressed at higher temperature. Altogether, these results provide new insights about how the genomic basis of adaptation to halo- and thermotolerance evolved in the green lineage.

Plastid and mitochondrial genomes of Coccophora langsdorfii (Fucales, Phaeophyceae) and the utility of molecular markers

Coccophora langsdorfii (Turner) Greville (Fucales) is an intertidal brown alga that is endemic to Northeast Asia and increasingly endangered by habitat loss and climate change. We sequenced the complete circular plastid and mitochondrial genomes of C. langsdorfii. The circular plastid genome is 124,450 bp and contains 139 protein-coding, 28 tRNA and 6 rRNA genes. The circular mitochondrial genome is 35,660 bp and contains 38 protein-coding, 25 tRNA and 3 rRNA genes. The structure and gene content of the C. langsdorfii plastid genome is similar to those of other species in the Fucales. The plastid genomes of brown algae in other orders share similar gene content but exhibit large structural recombination. The large in-frame insert in the cox2 gene in the mitochondrial genome of C. langsdorfii is typical of other brown algae. We explored the effect of this insertion on the structure and function of the cox2 protein. We estimated the usefulness of 135 plastid genes and 35 mitochondrial genes for developing molecular markers. This study shows that 29 organellar genes will prove efficient for resolving brown algal phylogeny. In addition, we propose a new molecular marker suitable for the study of intraspecific genetic diversity that should be tested in a large survey of populations of C. langsdorfii.

Evolution of novel mimicry rings facilitated by adaptive introgression in tropical butterflies

Understanding the genetic basis of phenotypic variation and the mechanisms involved in the evolution of adaptive novelty, especially in adaptive radiations, is a major goal in evolutionary biology. Here, we used whole-genome sequence data to investigate the origin of the yellow hindwing bar in the Heliconius cydno radiation. We found modular variation associated with hindwing phenotype in two narrow noncoding regions upstream and downstream of the cortex gene, which was recently identified as a pigmentation pattern controller in multiple species of Heliconius. Genetic variation at each of these modules suggests an independent control of the dorsal and ventral hindwing patterning, with the upstream module associated with the ventral phenotype and the downstream module with the dorsal one. Furthermore, we detected introgression between H. cydno and its closely related species Heliconius melpomene in these modules, likely allowing both species to participate in novel mimicry rings. In sum, our findings support the role of regulatory modularity coupled with adaptive introgression as an elegant mechanism by which novel phenotypic combinations can evolve and fuel an adaptive radiation.

Discovery and genomic analyses of hybridization between divergent lineages of Trypanosoma congolense, causative agent of Animal African Trypanosomiasis

Hybrid populations and introgressive hybridization remain poorly documented in pathogenic micro-organisms, as such that genetic exchange has been argued to play a minor role in their evolution. Recent work demonstrated the existence of hybrid microsatellite profiles in Trypanosoma congolense, a parasitic protozoan with detrimental effects on livestock productivity in sub-Saharan Africa. Here, we present the first population genomic study of T. congolense, revealing a remarkable number of single nucleotide polymorphisms (SNPs), small insertions/deletions (indels) and gene deletions among 56 parasite genomes from ten African countries. One group of parasites from Zambia was particularly diverse, displaying a substantial number of heterozygous SNP and indel sites compared to T. congolense parasites from the nine other sub-Saharan countries. Genomewide 5-kb phylogenetic analyses based on phased SNP data revealed that these parasites were the product of hybridization between phylogenetically distinct T. congolense lineages. Other parasites within the same region in Zambia presented a mosaic of haplotypic ancestry and genetic variability, indicating that hybrid parasites persisted and recombined beyond the initial hybridization event. Our observations challenge traditional views of trypanosome population biology and encourage future research on the role of hybridization in spreading genes for drug resistance, pathogenicity and virulence.

Aspergillus is monophyletic: Evidence from multiple gene phylogenies and extrolites profiles

Aspergillus is one of the economically most important fungal genera. Recently, the ICN adopted the single name nomenclature which has forced mycologists to choose one name for fungi (e.g. Aspergillus, Fusarium, Penicillium, etc.). Previously two proposals for the single name nomenclature in Aspergillus were presented: one attributes the name "Aspergillus" to clades comprising seven different teleomorphic names, by supporting the monophyly of this genus; the other proposes that Aspergillus is a non-monophyletic genus, by preserving the Aspergillus name only to species belonging to subgenus Circumdati and maintaining the sexual names in the other clades. The aim of our study was to test the monophyly of Aspergilli by two independent phylogenetic analyses using a multilocus phylogenetic approach. One test was run on the publicly available coding regions of six genes (RPB1, RPB2, Tsr1, Cct8, BenA, CaM), using 96 species of Penicillium, Aspergillus and related taxa. Bayesian (MrBayes) and Ultrafast Maximum Likelihood (IQ-Tree) and Rapid Maximum Likelihood (RaxML) analyses gave the same conclusion highly supporting the monophyly of Aspergillus. The other analyses were also performed by using publicly available data of the coding sequences of nine loci (18S rRNA, 5,8S rRNA, 28S rRNA (D1-D2), RPB1, RPB2, CaM, BenA, Tsr1, Cct8) of 204 different species. Both Bayesian (MrBayes) and Maximum Likelihood (RAxML) trees obtained by this second round of independent analyses strongly supported the monophyly of the genus Aspergillus. The stability test also confirmed the robustness of the results obtained. In conclusion, statistical analyses have rejected the hypothesis that the Aspergilli are non-monophyletic, and provided robust arguments that the genus is monophyletic and clearly separated from the monophyletic genus Penicillium. There is no phylogenetic evidence to split Aspergillus into several genera and the name Aspergillus can be used for all the species belonging to Aspergillus i.e. the clade comprising the subgenera Aspergillus, Circumdati, Fumigati, Nidulantes, section Cremei and certain species which were formerly part of the genera Phialosimplex and Polypaecilum. Section Cremei and the clade containing Polypaecilum and Phialosimplex are proposed as new subgenera of Aspergillus. The phylogenetic analysis also clearly shows that Aspergillus clavatoflavus and A. zonatus do not belong to the genus Aspergillus. Aspergillus clavatoflavus is therefore transferred to a new genus Aspergillago as Aspergillago clavatoflavus and A. zonatus was transferred to Penicilliopsis as P. zonata. The subgenera of Aspergillus share similar extrolite profiles indicating that the genus is one large genus from a chemotaxonomical point of view. Morphological and ecophysiological characteristics of the species also strongly indicate that Aspergillus is a polythetic class in phenotypic characters.

SICLE: a high-throughput tool for extracting evolutionary relationships from phylogenetic trees

We present the phylogeny analysis software SICLE (Sister Clade Extractor), an easy-to-use, high-throughput tool to describe the nearest neighbors to a node of interest in a phylogenetic tree as well as the support value for the relationship. The application is a command line utility that can be embedded into a phylogenetic analysis pipeline or can be used as a subroutine within another C++ program. As a test case, we applied this new tool to the published phylome of Salinibacter ruber, a species of halophilic Bacteriodetes, identifying 13 unique sister relationships to S. ruber across the 4,589 gene phylogenies. S. ruber grouped with bacteria, most often other Bacteriodetes, in the majority of phylogenies, but 91 phylogenies showed a branch-supported sister association between S. ruber and Archaea, an evolutionarily intriguing relationship indicative of horizontal gene transfer. This test case demonstrates how SICLE makes it possible to summarize the phylogenetic information produced by automated phylogenetic pipelines to rapidly identify and quantify the possible evolutionary relationships that merit further investigation. SICLE is available for free for noncommercial use at http://eebweb.arizona.edu/sicle/.