Phylogenomics Supports the Monophyly of Aphelids and Fungi and Identifies New Molecular Synapomorphies.

The supergroup Holomycota, composed of Fungi and several related lineages of unicellular organisms (Nucleariida, Rozellida, Microsporidia, and Aphelida), represents one of the major branches in the phylogeny of eukaryotes. Nevertheless, except for the well-established position of Nucleariida as the first holomycotan branch to diverge, the relationships among the other lineages have so far remained unresolved largely owing to the lack of molecular data for some groups. This was notably the case aphelids, a poorly known group of endobiotic phagotrophic protists that feed on algae with cellulose walls. The first molecular phylogenies including aphelids supported their sister relationship with Rozellida and Microsporidia which, collectively, formed a new group called Opisthosporidia (the "Opisthosporidia hypothesis"). However, recent phylogenomic analyses including massive sequence data from two aphelid genera, Paraphelidium and Amoeboaphelidium, suggested that the aphelids are sister to fungi (the "Aphelida $+$ Fungi hypothesis"). Should this position be confirmed, aphelids would be key to understanding the early evolution of Holomycota and the origin of Fungi. Here, we carry out phylogenomic analyses with an expanded taxonomic sampling for aphelids after sequencing the transcriptomes of two species of the genus Aphelidium (Aphelidium insulamus and Aphelidium tribonematis) in order to test these competing hypotheses. Our new phylogenomic analyses including species from the three known aphelid genera strongly rejected the Opisthosporidia hypothesis. Furthermore, comparative genomic analyses further supported the Aphelida $+$ Fungi hypothesis via the identification of 19 orthologous genes exclusively shared by these two lineages. Seven of them originated from ancient horizontal gene transfer events predating the aphelid-fungal split and the remaining 12 likely evolved de novo, constituting additional molecular synapomorphies for this clade. Ancestral trait reconstruction based on our well-resolved phylogeny of Holomycota suggests that the progenitor of both fungi and rozellids, was aphelid-like, having an amoeboflagellate state and likely preying endobiotically on cellulose-containing, cell-walled organisms. Two lineages, which we propose to call Phytophagea and Opisthophagea, evolved from this ancestor. Phytophagea, grouping aphelids and classical fungi, mainly specialized in endobiotic predation of algal cells. Fungi emerged from this lineage after losing phagotrophy in favor of osmotrophy. Opisthophagea, grouping rozellids and Microsporidia, became parasites, mostly of chitin-containing hosts. This lineage entered a progressive reductive process that resulted in a unique lifestyle, especially in the highly derived Microsporidia. [Aphelida, fungi, Holomycota, horizontal gene transfer, phylogenomics, synapomorphy.].

Unveiling microbial interactions in stratified mat communities from a warm saline shallow pond.

Modern phototrophic microbial mats are complex communities often used as analogs of major Precambrian ecosystems. Characterizing biotic, notably metabolic, interactions among different microbial mat members is essential to gain insights into the ecology and biogeochemistry of these systems. We applied 16S/18S rRNA metabarcoding approaches to characterize the structure of archaea, bacteria and protist communities from microbial mats collected along strong physicochemical (oxygen, salinity, temperature, depth) gradients in a shallow pond at the salar de Llamara (Chile). All mats were highly diverse, including members of virtually all known high-rank eukaryotic and prokaryotic taxa but also many novel lineages. Bacterial candidate divisions accounted for almost 50% of sequences in deeper mats, while Archaea represented up to 40% of sequences in some mat layers. Molecular phylogenetic analyses revealed six novel deeply divergent archaeal groups, along abundant and diverse Pacearchaeota and Woesearchaeota. Multivariate statistical analyses showed that local environmental conditions strongly influenced community composition. Co-occurrence network structure was markedly different between surface mats located in the oxygenated zone and mats located in transition and anoxic water layers. We identified potential biotic interactions between various high- and low-rank taxa. Notably, a strong positive correlation was observed between Lokiarchaeota and the poorly known candidate bacterial division TA06.

Extremely acidic proteomes and metabolic flexibility in bacteria and highly diversified archaea thriving in geothermal chaotropic brines.

Few described archaeal, and fewer bacterial, lineages thrive under salt-saturating conditions, such as solar saltern crystallizers (salinity above 30% w/v). They accumulate molar K+ cytoplasmic concentrations to maintain osmotic balance ('salt-in' strategy) and have proteins adaptively enriched in negatively charged acidic amino acids. Here we analysed metagenomes and metagenome-assembled genomes from geothermally influenced hypersaline ecosystems with increasing chaotropicity in the Danakil Depression. Normalized abundances of universal single-copy genes confirmed that haloarchaea and Nanohaloarchaeota encompass 99% of microbial communities in the near-life-limiting conditions of the Western-Canyon Lakes. Danakil metagenome- and metagenome-assembled-genome-inferred proteomes, compared with those of freshwater, seawater and solar saltern ponds up to saturation (6-14-32% salinity), showed that Western-Canyon Lake archaea encode the most acidic proteomes ever observed (median protein isoelectric points ≤4.4). We identified previously undescribed haloarchaeal families as well as an Aenigmatarchaeota family and a bacterial phylum independently adapted to extreme halophily. Despite phylum-level diversity decreasing with increasing salinity-chaotropicity, and unlike in solar salterns, adapted archaea exceedingly diversified in Danakil ecosystems, challenging the notion of decreasing diversity under extreme conditions. Metabolic flexibility to utilize multiple energy and carbon resources generated by local hydrothermalism along feast-and-famine strategies seemingly shapes microbial diversity in these ecosystems near life limits.

Expanded phylogeny of extremely halophilic archaea shows multiple independent adaptations to hypersaline environments.

Extremely halophilic archaea (Haloarchaea, Nanohaloarchaeota, Methanonatronarchaeia and Halarchaeoplasmatales) thrive in saturating salt concentrations where they must maintain osmotic equilibrium with their environment. The evolutionary history of adaptations enabling salt tolerance remains poorly understood, in particular because the phylogeny of several lineages is conflicting. Here we present a resolved phylogeny of extremely halophilic archaea obtained using improved taxon sampling and state-of-the-art phylogenetic approaches designed to cope with the strong compositional biases of their proteomes. We describe two uncultured lineages, Afararchaeaceae and Asbonarchaeaceae, which break the long branches at the base of Haloarchaea and Nanohaloarchaeota, respectively. We obtained 13 metagenome-assembled genomes (MAGs) of these archaea from metagenomes of hypersaline aquatic systems of the Danakil Depression (Ethiopia). Our phylogenomic analyses including these taxa show that at least four independent adaptations to extreme halophily occurred during archaeal evolution. Gene-tree/species-tree reconciliation suggests that gene duplication and horizontal gene transfer played an important role in this process, for example, by spreading key genes (such as those encoding potassium transporters) across extremely halophilic lineages.

Metagenome-derived virus-microbe ratios across ecosystems.

It is generally assumed that viruses outnumber cells on Earth by at least tenfold. Virus-to-microbe ratios (VMR) are largely based on counts of fluorescently labelled virus-like particles. However, these exclude intracellular viruses and potentially include false positives (DNA-containing vesicles, gene-transfer agents, unspecifically stained inert particles). Here, we develop a metagenome-based VMR estimate (mVRM) that accounts for DNA viruses across all stages of their replication cycles (virion, intracellular lytic and lysogenic) by using normalised RPKM (reads per kilobase of gene sequence per million of mapped metagenome reads) counts of the major capsid protein (MCP) genes and cellular universal single-copy genes (USCGs) as proxies for virus and cell counts, respectively. After benchmarking this strategy using mock metagenomes with increasing VMR, we inferred mVMR across different biomes. To properly estimate mVMR in aquatic ecosystems, we generated metagenomes from co-occurring cellular and viral fractions (>50 kDa-200 µm size-range) in freshwater, seawater and solar saltern ponds (10 metagenomes, 2 control metaviromes). Viruses outnumbered cells in freshwater by ~13 fold and in plankton from marine and saline waters by ~2-4 fold. However, across an additional set of 121 diverse non-aquatic metagenomes including microbial mats, microbialites, soils, freshwater and marine sediments and metazoan-associated microbiomes, viruses, on average, outnumbered cells by barely two-fold. Although viruses likely are the most diverse biological entities on Earth, their global numbers might be closer to those of cells than previously estimated.

Comparison of tryptophan biosynthetic operon regulation in different Gram-positive bacterial species.

The tryptophan biosynthetic operon has been widely used as a model system for studying transcription regulation. In Bacillus subtilis, the trp operon is primarily regulated by a tryptophan-activated RNA-binding protein, TRAP. Here we show that in many other Gram-positive species the trp operon is regulated differently, by tRNA(Trp) sensing by the RNA-based T-box mechanism, with T-boxes arranged in tandem. Our analyses reveal an apparent relationship between trp operon organization and the specific regulatory mechanism(s) used.

Functional shifts in microbial mats recapitulate early Earth metabolic transitions.

Phototrophic microbial mats dominated terrestrial ecosystems for billions of years, largely causing, through cyanobacterial oxygenic photosynthesis, but also undergoing, the Great Oxidation Event approximately 2.5 billion years ago. Taking a space-for-time approach based on the universality of core metabolic pathways expressed at ecosystem level, we studied gene content and co-occurrence networks in high-diversity metagenomes from spatially close microbial mats along a steep redox gradient. The observed functional shifts suggest that anoxygenic photosynthesis was present but not predominant under early Precambrian conditions, being accompanied by other autotrophic processes. Our data also suggest that, in contrast to general assumptions, anoxygenic photosynthesis largely expanded in parallel with the subsequent evolution of oxygenic photosynthesis and aerobic respiration. Finally, our observations might represent space-for-time evidence that the Wood-Ljungdahl carbon fixation pathway dominated phototrophic mats in early ecosystems, whereas the Calvin cycle probably evolved from pre-existing variants before becoming the dominant contemporary form of carbon fixation.

The Skin Microbiome of the Neotropical Frog Craugastor fitzingeri: Inferring Potential Bacterial-Host-Pathogen Interactions From Metagenomic Data.

Skin symbiotic bacteria on amphibians can play a role in protecting their host against pathogens. Chytridiomycosis, the disease caused by Batrachochytrium dendrobatidis, Bd, has caused dramatic population declines and extinctions of amphibians worldwide. Anti-Bd bacteria from amphibian skin have been cultured, and skin bacterial communities have been described through 16S rRNA gene amplicon sequencing. Here, we present a shotgun metagenomic analysis of skin bacterial communities from a Neotropical frog, Craugastor fitzingeri. We sequenced the metagenome of six frogs from two different sites in Panamá: three frogs from Soberanía (Sob), a Bd-endemic site, and three frogs from Serranía del Sapo (Sapo), a Bd-naïve site. We described the taxonomic composition of skin microbiomes and found that Pseudomonas was a major component of these communities. We also identified that Sob communities were enriched in Actinobacteria while Sapo communities were enriched in Gammaproteobacteria. We described gene abundances within the main functional classes and found genes enriched either in Sapo or Sob. We then focused our study on five functional classes of genes: biosynthesis of secondary metabolites, metabolism of terpenoids and polyketides, membrane transport, cellular communication and antimicrobial drug resistance. These gene classes are potentially involved in bacterial communication, bacterial-host and bacterial-pathogen interactions among other functions. We found that C. fitzingeri metagenomes have a wide array of genes that code for secondary metabolites, including antibiotics and bacterial toxins, which may be involved in bacterial communication, but could also have a defensive role against pathogens. Several genes involved in bacterial communication and bacterial-host interactions, such as biofilm formation and bacterial secretion systems were found. We identified specific genes and pathways enriched at the different sites and determined that gene co-occurrence networks differed between sites. Our results suggest that skin microbiomes are composed of distinct bacterial taxa with a wide range of metabolic capabilities involved in bacterial defense and communication. Differences in taxonomic composition and pathway enrichments suggest that skin microbiomes from different sites have unique functional properties. This study strongly supports the need for shotgun metagenomic analyses to describe the functional capacities of skin microbiomes and to tease apart their role in host defense against pathogens.

The genomic sequence of Exiguobacterium chiriqhucha str. N139 reveals a species that thrives in cold waters and extreme environmental conditions.

We report the genome sequence of Exiguobacterium chiriqhucha str. N139, isolated from a high-altitude Andean lake. Comparative genomic analyses of the Exiguobacterium genomes available suggest that our strain belongs to the same species as the previously reported E. pavilionensis str. RW-2 and Exiguobacterium str. GIC 31. We describe this species and propose the chiriqhucha name to group them. 'Chiri qhucha' in Quechua means 'cold lake', which is a common origin of these three cosmopolitan Exiguobacteria. The 2,952,588-bp E. chiriqhucha str. N139 genome contains one chromosome and three megaplasmids. The genome analysis of the Andean strain suggests the presence of enzymes that confer E. chiriqhucha str. N139 the ability to grow under multiple environmental extreme conditions, including high concentrations of different metals, high ultraviolet B radiation, scavenging for phosphorous and coping with high salinity. Moreover, the regulation of its tryptophan biosynthesis suggests that novel pathways remain to be discovered, and that these pathways might be fundamental in the amino acid metabolism of the microbial community from Laguna Negra, Argentina.

Extensive Identification of Bacterial Riboflavin Transporters and Their Distribution across Bacterial Species.

Riboflavin, the precursor for the cofactors flavin mononucleotide (FMN) and flavin adenine dinucleotide, is an essential metabolite in all organisms. While the functions for de novo riboflavin biosynthesis and riboflavin import may coexist in bacteria, the extent of this co-occurrence is undetermined. The RibM, RibN, RfuABCD and the energy-coupling factor-RibU bacterial riboflavin transporters have been experimentally characterized. In addition, ImpX, RfnT and RibXY are proposed as riboflavin transporters based on positional clustering with riboflavin biosynthetic pathway (RBP) genes or conservation of the FMN riboswitch regulatory element. Here, we searched for the FMN riboswitch in bacterial genomes to identify genes encoding riboflavin transporters and assessed their distribution among bacteria. Two new putative riboflavin transporters were identified: RibZ in Clostridium and RibV in Mesoplasma florum. Trans-complementation of an Escherichia coli riboflavin auxotroph strain confirmed the riboflavin transport activity of RibZ from Clostridium difficile, RibXY from Chloroflexus aurantiacus, ImpX from Fusobacterium nucleatum and RfnT from Ochrobactrum anthropi. The analysis of the genomic distribution of all known bacterial riboflavin transporters revealed that most occur in species possessing the RBP and that some bacteria may even encode functional riboflavin transporters from two different families. Our results indicate that some species possess ancestral riboflavin transporters, while others possess transporters that appear to have evolved recently. Moreover, our data suggest that unidentified riboflavin transporters also exist. The present study doubles the number of experimentally characterized riboflavin transporters and suggests a specific, non-accessory role for these proteins in riboflavin-prototrophic bacteria.

Biochemical features and functional implications of the RNA-based T-box regulatory mechanism.

The T-box mechanism is a common regulatory strategy used for modulating the expression of genes of amino acid metabolism-related operons in gram-positive bacteria, especially members of the Firmicutes. T-box regulation is usually based on a transcription attenuation mechanism in which an interaction between a specific uncharged tRNA and the 5' region of the transcript stabilizes an antiterminator structure in preference to a terminator structure, thereby preventing transcription termination. Although single T-box regulatory elements are common, double or triple T-box arrangements are also observed, expanding the regulatory range of these elements. In the present study, we predict the functional implications of T-box regulation in genes encoding aminoacyl-tRNA synthetases, proteins of amino acid biosynthetic pathways, transporters, and regulatory proteins. We also consider the global impact of the use of this regulatory mechanism on cell physiology. Novel biochemical relationships between regulated genes and their corresponding metabolic pathways were revealed. Some of the genes identified, such as the quorum-sensing gene luxS, in members of the Lactobacillaceae were not previously predicted to be regulated by the T-box mechanism. Our analyses also predict an imbalance in tRNA sensing during the regulation of operons containing multiple aminoacyl-tRNA synthetase genes or biosynthetic genes involved in pathways common to more than one amino acid. Based on the distribution of T-box regulatory elements, we propose that this regulatory mechanism originated in a common ancestor of members of the Firmicutes, Chloroflexi, Deinococcus-Thermus group, and Actinobacteria and was transferred into the Deltaproteobacteria by horizontal gene transfer.