Iroki: automatic customization and visualization of phylogenetic trees

Genetic Diversity within a Collection of Italian Maize Inbred Lines: A Resource for Maize Genomics and Breeding

Genetic diversity is fundamental for studying the complex architecture of the traits of agronomic importance, controlled by major and minor loci. Moreover, well-characterized germplasm collections are essential tools for dissecting and analyzing genetic and phenotypic diversity in crops. A panel of 360 entries, a subset of a larger collection maintained within the GenBank at CREA Bergamo, which includes the inbreds derived from traditional Italian maize open-pollinated (OP) varieties and advanced breeding ones (Elite Inbreds), was analyzed to identify SNP markers using the tGBS® genotyping-by-sequencing technology. A total of 797,368 SNPs were found during the initial analysis. Imputation and filtering processes were carried out based on the percentage of missing data, redundant markers, and rarest allele frequencies, resulting in a final dataset of 15,872 SNP markers for which a physical map position was identified. Using this dataset, the inbred panel was characterized for linkage disequilibrium (LD), genetic diversity, population structure, and genetic relationships. LD decay at a genome-wide level indicates that the collection is a suitable resource for association mapping. Population structure analyses, which were carried out with different clustering methods, showed stable grouping statistics for four groups, broadly corresponding to 'Insubria', 'Microsperma', and 'Scagliolino' genotypes, with a fourth group composed prevalently of elite accessions derived from Italian and US breeding programs. Based on these results, the CREA Italian maize collection, genetically characterized in this study, can be considered an important tool for the mapping and characterization of useful traits and associated loci/alleles, to be used in maize breeding programs.

On the Origin and Evolution of the Mosquito Male-determining Factor Nix

The mosquito family Culicidae is divided into 2 subfamilies named the Culicinae and Anophelinae. Nix, the dominant male-determining factor, has only been found in the culicines Aedes aegypti and Aedes albopictus, 2 important arboviral vectors that belong to the subgenus Stegomyia. Here we performed sex-specific whole-genome sequencing and RNAseq of divergent mosquito species and explored additional male-inclusive datasets to investigate the distribution of Nix. Except for the Culex genus, Nix homologs were found in all species surveyed from the Culicinae subfamily, including 12 additional species from 3 highly divergent tribes comprising 4 genera, suggesting Nix originated at least 133 to 165 million years ago (MYA). Heterologous expression of 1 of 3 divergent Nix open reading frames (ORFs) in Ae. aegypti resulted in partial masculinization of genetic females as evidenced by morphology and doublesex splicing. Phylogenetic analysis suggests Nix is related to femaleless (fle), a recently described intermediate sex-determining factor found exclusively in anopheline mosquitoes. Nix from all species has a conserved structure, including 3 RNA-recognition motifs (RRMs), as does fle. However, Nix has evolved at a much faster rate than fle. The RRM3 of both Nix and fle are distantly related to the single RRM of a widely distributed and conserved splicing factor transformer-2 (tra2). The RRM3-based phylogenetic analysis suggests this domain in Nix and fle may have evolved from tra2 or a tra2-related gene in a common ancestor of mosquitoes. Our results provide insights into the evolution of sex determination in mosquitoes and will inform broad applications of mosquito-control strategies based on manipulating sex ratios toward nonbiting males.

Gallionellaceae pangenomic analysis reveals insight into phylogeny, metabolic flexibility, and iron oxidation mechanisms

IMPORTANCE

Neutrophilic iron-oxidizing bacteria (FeOB) produce copious iron (oxyhydr)oxides that can profoundly influence biogeochemical cycles, notably the fate of carbon and many metals. To fully understand environmental microbial iron oxidation, we need a thorough accounting of iron oxidation mechanisms. In this study, we show the Gallionellaceae FeOB genomes encode both characterized iron oxidases as well as uncharacterized multiheme cytochromes (MHCs). MHCs are predicted to transfer electrons from extracellular substrates and likely confer metabolic capabilities that help Gallionellaceae occupy a range of different iron- and mineral-rich niches. Gallionellaceae appear to specialize in iron oxidation, so it would be advantageous for them to have multiple mechanisms to oxidize various forms of iron, given the many iron minerals on Earth, as well as the physiological and kinetic challenges faced by FeOB. The multiple iron/mineral oxidation mechanisms may help drive the widespread ecological success of Gallionellaceae.

Gallionellaceae in rice root plaque: metabolic roles in iron oxidation, nutrient cycling, and plant interactions

IMPORTANCE

In waterlogged soils, iron plaque forms a reactive barrier between the root and soil, collecting phosphate and metals such as arsenic and cadmium. It is well established that iron-reducing bacteria solubilize iron, releasing these associated elements. In contrast, microbial roles in plaque formation have not been clear. Here, we show that there is a substantial population of iron oxidizers in plaque, and furthermore, that these organisms (Sideroxydans and Gallionella) are distinguished by genes for plant colonization and nutrient fixation. Our results suggest that iron-oxidizing and iron-reducing bacteria form and remodel iron plaque, making it a dynamic system that represents both a temporary sink for elements (P, As, Cd, C, etc.) as well as a source. In contrast to abiotic iron oxidation, microbial iron oxidation results in coupled Fe-C-N cycling, as well as microbe-microbe and microbe-plant ecological interactions that need to be considered in soil biogeochemistry, ecosystem dynamics, and crop management.

Accurate prediction of HLA class II antigen presentation across all loci using tailored data acquisition and refined machine learning

Accurate prediction of antigen presentation by human leukocyte antigen (HLA) class II molecules is crucial for rational development of immunotherapies and vaccines targeting CD4+ T cell activation. So far, most prediction methods for HLA class II antigen presentation have focused on HLA-DR because of limited availability of immunopeptidomics data for HLA-DQ and HLA-DP while not taking into account alternative peptide binding modes. We present an update to the NetMHCIIpan prediction method, which closes the performance gap between all three HLA class II loci. We accomplish this by first integrating large immunopeptidomics datasets describing the HLA class II specificity space across all loci using a refined machine learning framework that accommodates inverted peptide binders. Next, we apply targeted immunopeptidomics assays to generate data that covers additional HLA-DP specificities. The final method, NetMHCIIpan-4.3, achieves high accuracy and molecular coverage across all HLA class II allotypes.

Gut microbiomes of cycad-feeding insects tolerant to β-methylamino-L-alanine (BMAA) are rich in siderophore biosynthesis

Ingestion of the cycad toxins β-methylamino-L-alanine (BMAA) and azoxyglycosides is harmful to diverse organisms. However, some insects are specialized to feed on toxin-rich cycads with apparent immunity. Some cycad-feeding insects possess a common set of gut bacteria, which might play a role in detoxifying cycad toxins. Here, we investigated the composition of gut microbiota from a worldwide sample of cycadivorous insects and characterized the biosynthetic potential of selected bacteria. Cycadivorous insects shared a core gut microbiome consisting of six bacterial taxa, mainly belonging to the Proteobacteria, which we were able to isolate. To further investigate selected taxa from diverging lineages, we performed shotgun metagenomic sequencing of co-cultured bacterial sub-communities. We characterized the biosynthetic potential of four bacteria from Serratia, Pantoea, and two different Stenotrophomonas lineages, and discovered a suite of biosynthetic gene clusters notably rich in siderophores. Siderophore semi-untargeted metabolomics revealed a broad range of chemically related yet diverse iron-chelating metabolites, including desferrioxamine B, suggesting the occurrence of an unprecedented desferrioxamine-like biosynthetic pathway that remains to be identified. These results provide a foundation for future investigations into how cycadivorous insects tolerate diets rich in azoxyglycosides, BMAA, and other cycad toxins, including a possible role for bacterial siderophores.

The complete mitochondrial genome of Cyphocaris challengeri (Amphipoda: Cyphocarididae)

The amphipod Cyphocaris challengeri is a globally distributed, highly abundant species of zooplankton. Here, we report the complete mitochondrial genome of C. challengeri obtained using the Illumina sequencing platform from a specimen collected from Puget Sound, Washington. The mitogenome is a circular DNA molecule with a size of 14,338 bp and 26.7% GC content, with 13 protein-encoding genes, 2 rRNAs, and 22 tRNAs annotated. A maximum likelihood phylogenetic analysis including C. challengeri and all other available mitogenomes from Amphipoda places our mitogenome firmly within the Lysianassoidea superfamily, as expected. The newly described mitochondrial genome of C. challengeri fills a gap in valuable reference data for detecting this organism using molecular methods such as environmental DNA.

Ubiquitous, B12-dependent virioplankton utilizing ribonucleotide-triphosphate reductase demonstrate interseasonal dynamics and associate with a diverse range of bacterial hosts in the pelagic ocean

Through infection and lysis of their coexisting bacterial hosts, viruses impact the biogeochemical cycles sustaining globally significant pelagic oceanic ecosystems. Currently, little is known of the ecological interactions between lytic viruses and their bacterial hosts underlying these biogeochemical impacts at ecosystem scales. This study focused on populations of lytic viruses carrying the B12-dependent Class II monomeric ribonucleotide reductase (RNR) gene, ribonucleotide-triphosphate reductase (Class II RTPR), documenting seasonal changes in pelagic virioplankton and bacterioplankton using amplicon sequences of Class II RTPR and the 16S rRNA gene, respectively. Amplicon sequence libraries were analyzed using compositional data analysis tools that account for the compositional nature of these data. Both virio- and bacterioplankton communities responded to environmental changes typically seen across seasonal cycles as well as shorter term upwelling-downwelling events. Defining Class II RTPR-carrying viral populations according to major phylogenetic clades proved a more robust means of exploring virioplankton ecology than operational taxonomic units defined by percent sequence homology. Virioplankton Class II RTPR populations showed positive associations with a broad phylogenetic diversity of bacterioplankton including dominant taxa within pelagic oceanic ecosystems such as Prochlorococcus and SAR11. Temporal changes in Class II RTPR virioplankton, occurring as both free viruses and within infected cells, indicated possible viral-host pairs undergoing sustained infection and lysis cycles throughout the seasonal study. Phylogenetic relationships inferred from Class II RTPR sequences mirrored ecological patterns in virio- and bacterioplankton populations demonstrating possible genome to phenome associations for an essential viral replication gene.

Genome-wide analysis of horizontal transfer in non-model wild species from a natural ecosystem reveals new insights into genetic exchange in plants

Horizontal transfer (HT) refers to the exchange of genetic material between divergent species by mechanisms other than reproduction. In recent years, several studies have demonstrated HTs in eukaryotes, particularly in the context of parasitic relationships and in model species. However, very little is known about HT in natural ecosystems, especially those involving non-parasitic wild species, and the nature of the ecological relationships that promote these HTs. In this work, we conducted a pilot study investigating HTs by sequencing the genomes of 17 wild non-model species from a natural ecosystem, the Massane forest, located in southern France. To this end, we developed a new computational pipeline called INTERCHANGE that is able to characterize HTs at the whole genome level without prior annotation and directly in the raw sequencing reads. Using this pipeline, we identified 12 HT events, half of which occurred between lianas and trees. We found that mainly low copy number LTR-retrotransposons from the Copia superfamily were transferred between these wild plant species, especially those of the Ivana and Ale lineages. This study revealed a possible new route for HTs between non-parasitic plants and provides new insights into the genomic characteristics of horizontally transferred DNA in plant genomes.

Genetic Basis of Acinetobacter sp. K1 Adaptation Mechanisms to Extreme Environmental Conditions

Anthropogenic pollution often leads to the generation of technosols, technogenic soils with inhospitable conditions for all living organisms including microbiota. Aluminum production near Ziar nad Hronom (Slovakia) resulted in the creation of a highly alkaline and heavy-metal-rich brown mud landfill, from which a bacterial strain of a likely new species of the genus Acinetobacter, Acinetobacter sp. K1, was isolated. The whole-genome sequence analysis of this strain confirmed the presence of operon units enabling tolerance to the heavy metals copper, zinc, cobalt, cadmium, chromium, and metalloid arsenic, which are functionally active. Despite the predominance of plasmid-related sequences in the K1 genome, the results indicate that most of the resistance genes are chromosomally encoded. No significant alkali tolerance of Acinetobacter sp. K1 was observed in vitro, suggesting that community level mechanisms are responsible for the survival of this strain in the highly alkaline, brown mud bacterial community.

Machine learning reveals limited contribution of trans-only encoded variants to the HLA-DQ immunopeptidome

Human leukocyte antigen (HLA) class II antigen presentation is key for controlling and triggering T cell immune responses. HLA-DQ molecules, which are believed to play a major role in autoimmune diseases, are heterodimers that can be formed as both cis and trans variants depending on whether the α- and β-chains are encoded on the same (cis) or opposite (trans) chromosomes. So far, limited progress has been made for predicting HLA-DQ antigen presentation. In addition, the contribution of trans-only variants (i.e. variants not observed in the population as cis) in shaping the HLA-DQ immunopeptidome remains largely unresolved. Here, we seek to address these issues by integrating state-of-the-art immunoinformatics data mining models with large volumes of high-quality HLA-DQ specific mass spectrometry immunopeptidomics data. The analysis demonstrates highly improved predictive power and molecular coverage for models trained including these novel HLA-DQ data. More importantly, investigating the role of trans-only HLA-DQ variants reveals a limited to no contribution to the overall HLA-DQ immunopeptidome. In conclusion, this study furthers our understanding of HLA-DQ specificities and casts light on the relative role of cis versus trans-only HLA-DQ variants in the HLA class II antigen presentation space. The developed method, NetMHCIIpan-4.2, is available at https://services.healthtech.dtu.dk/services/NetMHCIIpan-4.2 .

A zero-agnostic model for copy number evolution in cancer

Motivation

New low-coverage single-cell DNA sequencing technologies enable the measurement of copy number profiles from thousands of individual cells within tumors. From this data, one can infer the evolutionary history of the tumor by modeling transformations of the genome via copy number aberrations. A widely used model to infer such copy number phylogenies is the copy number transformation (CNT) model in which a genome is represented by an integer vector and a copy number aberration is an event that either increases or decreases the number of copies of a contiguous segment of the genome. The CNT distance between a pair of copy number profiles is the minimum number of events required to transform one profile to another. While this distance can be computed efficiently, no efficient algorithm has been developed to find the most parsimonious phylogeny under the CNT model.

Results

We introduce the zero-agnostic copy number transformation (ZCNT) model, a simplification of the CNT model that allows the amplification or deletion of regions with zero copies. We derive a closed form expression for the ZCNT distance between two copy number profiles and show that, unlike the CNT distance, the ZCNT distance forms a metric. We leverage the closed-form expression for the ZCNT distance and an alternative characterization of copy number profiles to derive polynomial time algorithms for two natural relaxations of the small parsimony problem on copy number profiles. While the alteration of zero copy number regions allowed under the ZCNT model is not biologically realistic, we show on both simulated and real datasets that the ZCNT distance is a close approximation to the CNT distance. Extending our polynomial time algorithm for the ZCNT small parsimony problem, we develop an algorithm, Lazac, for solving the large parsimony problem on copy number profiles. We demonstrate that Lazac outperforms existing methods for inferring copy number phylogenies on both simulated and real data.

Yersiniomics, a Multi-Omics Interactive Database for Yersinia Species

The genus Yersinia includes a large variety of nonpathogenic and life-threatening pathogenic bacteria, which cause a broad spectrum of diseases in humans and animals, such as plague, enteritis, Far East scarlet-like fever (FESLF), and enteric redmouth disease. Like most clinically relevant microorganisms, Yersinia spp. are currently subjected to intense multi-omics investigations whose numbers have increased extensively in recent years, generating massive amounts of data useful for diagnostic and therapeutic developments. The lack of a simple and centralized way to exploit these data led us to design Yersiniomics, a web-based platform allowing straightforward analysis of Yersinia omics data. Yersiniomics contains a curated multi-omics database at its core, gathering 200 genomic, 317 transcriptomic, and 62 proteomic data sets for Yersinia species. It integrates genomic, transcriptomic, and proteomic browsers, a genome viewer, and a heatmap viewer to navigate within genomes and experimental conditions. For streamlined access to structural and functional properties, it directly links each gene to GenBank, the Kyoto Encyclopedia of Genes and Genomes (KEGG), UniProt, InterPro, IntAct, and the Search Tool for the Retrieval of Interacting Genes/Proteins (STRING) and each experiment to Gene Expression Omnibus (GEO), the European Nucleotide Archive (ENA), or the Proteomics Identifications Database (PRIDE). Yersiniomics provides a powerful tool for microbiologists to assist with investigations ranging from specific gene studies to systems biology studies. IMPORTANCE The expanding genus Yersinia is composed of multiple nonpathogenic species and a few pathogenic species, including the deadly etiologic agent of plague, Yersinia pestis. In 2 decades, the number of genomic, transcriptomic, and proteomic studies on Yersinia grew massively, delivering a wealth of data. We developed Yersiniomics, an interactive web-based platform, to centralize and analyze omics data sets on Yersinia species. The platform allows user-friendly navigation between genomic data, expression data, and experimental conditions. Yersiniomics will be a valuable tool to microbiologists.

Gallionellaceae pangenomic analysis reveals insight into phylogeny, metabolic flexibility, and iron oxidation mechanisms

The iron-oxidizing Gallionellaceae drive a wide variety of biogeochemical cycles through their metabolisms and biominerals. To better understand the environmental impacts of Gallionellaceae, we need to improve our knowledge of their diversity and metabolisms, especially any novel iron oxidation mechanisms. Here, we used a pangenomic analysis of 103 genomes to resolve Gallionellaceae phylogeny and explore the range of genomic potential. Using a concatenated ribosomal protein tree and key gene patterns, we determined Gallionellaceae has four genera, divided into two groups-iron-oxidizing bacteria (FeOB) Gallionella, Sideroxydans, and Ferriphaselus with known iron oxidases (Cyc2, MtoA) and nitrite-oxidizing bacteria (NOB) Candidatus Nitrotoga with nitrite oxidase (Nxr). The FeOB and NOB have similar electron transport chains, including genes for reverse electron transport and carbon fixation. Auxiliary energy metabolisms including S oxidation, denitrification, and organotrophy were scattered throughout the Gallionellaceae FeOB. Within FeOB, we found genes that may represent adaptations for iron oxidation, including a variety of extracellular electron uptake (EEU) mechanisms. FeOB genomes encoded more predicted c-type cytochromes overall, notably more multiheme c-type cytochromes (MHCs) with >10 CXXCH motifs. These include homologs of several predicted outer membrane porin-MHC complexes, including MtoAB and Uet. MHCs are known to efficiently conduct electrons across longer distances and function across a wide range of redox potentials that overlap with mineral redox potentials, which can help expand the range of usable iron substrates. Overall, the results of pangenome analyses suggest that the Gallionellaceae genera Gallionella, Sideroxydans, and Ferriphaselus are primarily iron oxidizers, capable of oxidizing dissolved Fe2+ as well as a range of solid iron or other mineral substrates.

Tomato brown rugose fruit virus Mo gene is a novel microbial source tracking marker

Microbial source tracking (MST) identifies sources of fecal contamination in the environment using fecal host-associated markers. While there are numerous bacterial MST markers, there are few viral markers. Here we design and test novel viral MST markers based on tomato brown rugose fruit virus (ToBRFV) genomes. We assembled eight nearly complete genomes of ToBRFV from wastewater and stool samples from the San Francisco Bay Area in the United States of America. Next, we developed two novel probe-based RT-PCR assays based on conserved regions of the ToBRFV genome, and tested the markers’ sensitivities and specificities using human and non-human animal stool as well as wastewater. TheToBRFV markers are sensitive and specific; in human stool and wastewater, they are more prevalent and abundant than a currently used marker, the pepper mild mottle virus (PMMoV) coat protein (CP) gene. We applied the assays to detect fecal contamination in urban stormwater samples and found that the ToBRFV markers matched cross-assembly phage (crAssphage), an established viral MST marker, in prevalence across samples. Taken together, ToBRFV is a promising viral human-associated MST marker.

Importance

Human exposure to fecal contamination in the environment can cause transmission of infectious diseases. Microbial source tracking (MST) can identify sources of fecal contamination so that contamination can be remediated and human exposures can be reduced. MST requires the use of fecal host-associated MST markers. Here we design and test novel MST markers from genomes of tomato brown rugose fruit virus (ToBRFV). The markers are sensitive and specific to human stool, and highly abundant in human stool and wastewater samples.

Orchid mycorrhizal fungi and ascomycetous fungi in epiphytic Vanda falcata roots occupy different niches during growth and development

Epiphytic orchids are commonly found in exposed environments, which plausibly lead to different root fungal community structures from terrestrial orchids. Until recently, few studies have been conducted to show the fungal community structure during the growth of a photosynthetic and epiphytic orchid in its natural growing site. In this study, the Vanda falcata (commonly known as Neofinetia falcata), one of Japan's ornamental orchids, was used to characterize the fungal community structure at different developmental stages. Amplicon sequencing analysis showed that all development stages contain a similar fungal community: Ascomycota dominate half of the community while one-third of the community belongs to Basidiomycota. Rhizoctonia-like fungi, a polyphyletic basidiomycetous fungal group forming mycorrhizas in many orchids, exist even in a smaller portion (around one-quarter) compared to other Basidiomycota members. While ascomycetous fungi exhibit pathogenicity, two Ceratobasidium strains isolated from young and adult plants could initiate seed germination in vitro. It was also found that the colonization of mycorrhizal fungi was concentrated in a part of the root where it directly attaches to the phorophyte bark, while ascomycetous fungi were distributed in the velamen but never colonized cortical cells. Additionally, the root parts attached to the bark have denser exodermal passage cells, and these cells were only colonized by mycorrhizal fungi that further penetrated into the cortical area. Therefore, we confirmed a process that physical regulation of fungal entry to partition the ascomycetes and mycorrhizal fungi results in the balanced mycorrhizal symbiosis in this orchid.

Conservation of Energetic Pathways for Electroautotrophy in the Uncultivated Candidate Order Tenderiales

Electromicrobiology can be used to understand extracellular electron uptake in previously undescribed chemolithotrophs. Enrichment and characterization of the uncultivated electroautotroph "Candidatus Tenderia electrophaga" using electromicrobiology led to the designation of the order Tenderiales. Representative Tenderiales metagenome-assembled genomes (MAGs) have been identified in a number of environmental surveys, yet a comprehensive characterization of conserved genes for extracellular electron uptake has thus far not been conducted. Using comparative genomics, we identified conserved orthologous genes within the Tenderiales and nearest-neighbor orders important for extracellular electron uptake based on a previously proposed pathway from "Ca. Tenderia electrophaga." The Tenderiales contained a conserved cluster we designated uetABCDEFGHIJ, which encodes proteins containing features that would enable transport of extracellular electrons to cytoplasmic membrane-bound energy-transducing complexes such as two conserved cytochrome cbb₃ oxidases. For example, UetJ is predicted to be an extracellular undecaheme c-type cytochrome that forms a heme wire. We also identified clusters of genes predicted to facilitate assembly and maturation of electron transport proteins, as well as cellular attachment to surfaces. Autotrophy among the Tenderiales is supported by the presence of carbon fixation and stress response pathways that could allow cellular growth by extracellular electron uptake. Key differences between the Tenderiales and other known neutrophilic iron oxidizers were revealed, including very few Cyc2 genes in the Tenderiales. Our results reveal a possible conserved pathway for extracellular electron uptake and suggest that the Tenderiales have an ecological role in coupling metal or mineral redox chemistry and the carbon cycle in marine and brackish sediments. IMPORTANCE Chemolithotrophic bacteria capable of extracellular electron uptake to drive energy metabolism and CO₂ fixation are known as electroautotrophs. The recently described order Tenderiales contains the uncultivated electroautotroph "Ca. Tenderia electrophaga." The "Ca. Tenderia electrophaga" genome contains genes proposed to make up a previously undescribed extracellular electron uptake pathway. Here, we use comparative genomics to show that this pathway is well conserved among Tenderiales spp. recovered by metagenome-assembled genomes. This conservation extends to near neighbors of the Tenderiales but not to other well-studied chemolithotrophs, including iron and sulfur oxidizers, indicating that these genes may be useful markers of growth using insoluble extracellular electron donors. Our findings suggest that extracellular electron uptake and electroautotrophy may be pervasive among the Tenderiales, and the geographic locations from which metagenome-assembled genomes were recovered offer clues to their natural ecological niche.

Bioremediation of copper in sediments from a constructed wetland ex situ with the novel bacterium Cupriavidus basilensis SRS

The H-02 constructed wetland was designed to remove metals (primarily copper and zinc) to treat building process water and storm water runoff from multiple sources associated with the Tritium Facility at the DOE-Savannah River Site, Aiken, SC. The concentration of Cu and Zn in the sediments has increased over the lifetime of the wetland and is a concern. A bioremediation option was investigated at the laboratory scale utilizing a newly isolated bacterium of the copper metabolizing genus Cupriavidus isolated from Tim's Branch Creek, a second-order stream that eventually serves as a tributary to the Savannah River, contaminated with uranium and other metals including copper, nickel, and mercury. Cupriavidus basilensis SRS is a rod-shaped, gram-negative bacterium which has been shown to have predatory tendencies. The isolate displayed resistance to the antibiotics ofloxacin, tetracycline, ciprofloxacin, select fungi, as well as Cu²⁺ and Zn²⁺. Subsequent ribosomal sequencing demonstrated a 100% confidence for placement in the genus Cupriavidus and a 99.014% match to the C. basilensis type strain. When H-02 wetland samples were inoculated with Cupriavidus basilensis SRS samples showed significant (p < 0.05) decrease in Cu²⁺ concentrations and variability in Zn²⁺ concentrations. Over the 72-h incubation there were no significant changes in the inoculate densities (10⁶-10⁸ cells/ML) indicating Cupriavidus basilensis SRS resiliency in this environment. This research expands our understanding of the Cupriavidus genus and demonstrates the potential for Cupriavidus basilensis SRS to bioremediate sites impacted with heavy metals, most notably copper.

The rate and role of pseudogenes of the Mycobacterium tuberculosis complex

Whole-genome sequence analyses have significantly contributed to the understanding of virulence and evolution of the Mycobacterium tuberculosis complex (MTBC), the causative pathogens of tuberculosis. Most MTBC evolutionary studies are focused on single nucleotide polymorphisms and deletions, but rare studies have evaluated gene content, whereas none has comprehensively evaluated pseudogenes. Accordingly, we describe an extensive study focused on quantifying and predicting possible functions of MTBC and Mycobacterium canettii pseudogenes. Using NCBI's PGAP-detected pseudogenes, we analysed 25 837 pseudogenes from 158 MTBC and M. canetii strains and combined transcriptomics and proteomics of M. tuberculosis H37Rv to gain insights about pseudogenes' expression. Our results indicate significant variability concerning rate and conservancy of in silico predicted pseudogenes among different ecotypes and lineages of tuberculous mycobacteria and pseudogenization of important virulence factors and genes of the metabolism and antimicrobial resistance/tolerance. We show that in silico predicted pseudogenes contribute considerably to MTBC genetic diversity at the population level. Moreover, the transcription machinery of M. tuberculosis can fully transcribe most pseudogenes, indicating intact promoters and recent pseudogene evolutionary emergence. Proteomics of M. tuberculosis and close evaluation of mutational lesions driving pseudogenization suggest that few in silico predicted pseudogenes are likely capable of neofunctionalization, nonsense mutation reversal, or phase variation, contradicting the classical definition of pseudogenes. Such findings indicate that genome annotation should be accompanied by proteomics and protein function assays to improve its accuracy. While indels and insertion sequences are the main drivers of the observed mutational lesions in these species, population bottlenecks and genetic drift are likely the evolutionary processes acting on pseudogenes' emergence over time. Our findings unveil a new perspective on MTBC's evolution and genetic diversity.

Rice apoplastic CBM1-interacting protein counters blast pathogen invasion by binding conserved carbohydrate binding module 1 motif of fungal proteins

When infecting plants, fungal pathogens secrete cell wall-degrading enzymes (CWDEs) that break down cellulose and hemicellulose, the primary components of plant cell walls. Some fungal CWDEs contain a unique domain, named the carbohydrate binding module (CBM), that facilitates their access to polysaccharides. However, little is known about how plants counteract pathogen degradation of their cell walls. Here, we show that the rice cysteine-rich repeat secretion protein OsRMC binds to and inhibits xylanase MoCel10A of the blast fungus pathogen Magnaporthe oryzae, interfering with its access to the rice cell wall and degradation of rice xylan. We found binding of OsRMC to various CBM1-containing enzymes, suggesting that it has a general role in inhibiting the action of CBM1. OsRMC is localized to the apoplast, and its expression is strongly induced in leaves infected with M. oryzae. Remarkably, knockdown and overexpression of OsRMC reduced and enhanced rice defense against M. oryzae, respectively, demonstrating that inhibition of CBM1-containing fungal enzymes by OsRMC is crucial for rice defense. We also identified additional CBM-interacting proteins (CBMIPs) from Arabidopsis thaliana and Setaria italica, indicating that a wide range of plants counteract pathogens through this mechanism.

Plants have evolved various activity-inhibiting proteins as a defense against fungal cell wall-degrading enzymes (CWDEs), but how plants counteract the function of fungal enzymes containing carbohydrate binding modules (CBMs) remains unknown. Here, we demonstrate that OsRMC, a member of the cysteine-rich repeat secretion protein family, interacts with fungal CBM1. OsRMC binding to CBM1 of a blast fungal xylanase blocks access to cellulose, resulting in the inhibition of xylanase enzymatic activity. Our study provides significant insights into plant countermeasures against CWDEs in the apoplastic space during plant-fungal pathogen interactions. It also reveals a molecular function of the DUF26 domain widely distributed in plant proteins.

Broadly neutralizing antibodies to SARS-related viruses can be readily induced in rhesus macaques

To prepare for future coronavirus (CoV) pandemics, it is desirable to generate vaccines capable of eliciting broadly neutralizing antibody responses to CoVs. Here, we show that immunization of macaques with SARS-CoV-2 spike (S) protein with a two-shot protocol generated potent serum receptor binding domain cross-neutralizing antibody responses to both SARS-CoV-2 and SARS-CoV-1. Furthermore, responses were equally effective against most SARS-CoV-2 variants of concern (VOCs) and some were highly effective against Omicron. This result contrasts with human infection or many two-shot vaccination protocols where responses were typically more SARS-CoV-2 specific and where VOCs were less well neutralized. Structural studies showed that cloned macaque neutralizing antibodies, particularly using a given heavy chain germline gene, recognized a relatively conserved region proximal to the angiotensin converting enzyme 2 receptor binding site (RBS), whereas many frequently elicited human neutralizing antibodies targeted more variable epitopes overlapping the RBS. B cell repertoire differences between humans and macaques appeared to influence the vaccine response. The macaque neutralizing antibodies identified a pan-SARS-related virus epitope region less well targeted by human antibodies that could be exploited in rational vaccine design.

Advancing Precision Vaccinology by Molecular and Genomic Surveillance of Severe Acute Respiratory Syndrome Coronavirus 2 in Germany, 2021

BACKGROUND

Comprehensive pathogen genomic surveillance represents a powerful tool to complement and advance precision vaccinology. The emergence of the Alpha variant in December 2020 and the resulting efforts to track the spread of this and other severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) variants of concern led to an expansion of genomic sequencing activities in Germany.

METHODS

At Robert Koch Institute (RKI), the German National Institute of Public Health, we established the Integrated Molecular Surveillance for SARS-CoV-2 (IMS-SC2) network to perform SARS-CoV-2 genomic surveillance at the national scale, SARS-CoV-2-positive samples from laboratories distributed across Germany regularly undergo whole-genome sequencing at RKI.

RESULTS

We report analyses of 3623 SARS-CoV-2 genomes collected between December 2020 and December 2021, of which 3282 were randomly sampled. All variants of concern were identified in the sequenced sample set, at ratios equivalent to those in the 100-fold larger German GISAID sequence dataset from the same time period. Phylogenetic analysis confirmed variant assignments. Multiple mutations of concern emerged during the observation period. To model vaccine effectiveness in vitro, we employed authentic-virus neutralization assays, confirming that both the Beta and Zeta variants are capable of immune evasion. The IMS-SC2 sequence dataset facilitated an estimate of the SARS-CoV-2 incidence based on genetic evolution rates. Together with modeled vaccine efficacies, Delta-specific incidence estimation indicated that the German vaccination campaign contributed substantially to a deceleration of the nascent German Delta wave.

CONCLUSIONS

SARS-CoV-2 molecular and genomic surveillance may inform public health policies including vaccination strategies and enable a proactive approach to controlling coronavirus disease 2019 spread as the virus evolves.

Tree House Explorer: A Novel Genome Browser for Phylogenomics

Tree House Explorer (THEx) is a genome browser that integrates phylogenomic data and genomic annotations into a single interactive platform for combined analysis. THEx allows users to visualize genome-wide variation in evolutionary histories and genetic divergence on a chromosome-by-chromosome basis, with continuous sliding window comparisons to gene annotations (GFF/GTF), recombination rates, and other user-specified, highly customizable feature annotations. THEx provides a new platform for interactive phylogenomic data visualization to analyze and interpret the diverse evolutionary histories woven throughout genomes. Hosted on Conda, THEx integrates seamlessly into new or pre-existing workflows.

Intradiol ring cleavage dioxygenases from herbivorous spider mites as a new detoxification enzyme family in animals

Background

Generalist herbivores such as the two-spotted spider mite Tetranychus urticae thrive on a wide variety of plants and can rapidly adapt to novel hosts. What traits enable polyphagous herbivores to cope with the diversity of secondary metabolites in their variable plant diet is unclear. Genome sequencing of T. urticae revealed the presence of 17 genes that code for secreted proteins with strong homology to “intradiol ring cleavage dioxygenases (DOGs)” from bacteria and fungi, and phylogenetic analyses show that they have been acquired by horizontal gene transfer from fungi. In bacteria and fungi, DOGs have been well characterized and cleave aromatic rings in catecholic compounds between adjacent hydroxyl groups. Such compounds are found in high amounts in solanaceous plants like tomato, where they protect against herbivory. To better understand the role of this gene family in spider mites, we used a multi-disciplinary approach to functionally characterize the various T. urticae DOG genes.

Results

We confirmed that DOG genes were present in the T. urticae genome and performed a phylogenetic reconstruction using transcriptomic and genomic data to advance our understanding of the evolutionary history of spider mite DOG genes. We found that DOG expression differed between mites from different plant hosts and was induced in response to jasmonic acid defense signaling. In consonance with a presumed role in detoxification, expression was localized in the mite’s gut region. Silencing selected DOGs expression by dsRNA injection reduced the mites’ survival rate on tomato, further supporting a role in mitigating the plant defense response. Recombinant purified DOGs displayed a broad substrate promiscuity, cleaving a surprisingly wide array of aromatic plant metabolites, greatly exceeding the metabolic capacity of previously characterized microbial DOGs.

Conclusion

Our findings suggest that the laterally acquired spider mite DOGs function as detoxification enzymes in the gut, disarming plant metabolites before they reach toxic levels. We provide experimental evidence to support the hypothesis that this proliferated gene family in T. urticae is causally linked to its ability to feed on an extremely wide range of host plants.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12915-022-01323-1.

Targeted isolation of diverse human protective broadly neutralizing antibodies against SARS-like viruses

The emergence of current severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) variants of concern (VOCs) and potential future spillovers of SARS-like coronaviruses into humans pose a major threat to human health and the global economy. Development of broadly effective coronavirus vaccines that can mitigate these threats is needed. Here, we utilized a targeted donor selection strategy to isolate a large panel of human broadly neutralizing antibodies (bnAbs) to sarbecoviruses. Many of these bnAbs are remarkably effective in neutralizing a diversity of sarbecoviruses and against most SARS-CoV-2 VOCs, including the Omicron variant. Neutralization breadth is achieved by bnAb binding to epitopes on a relatively conserved face of the receptor-binding domain (RBD). Consistent with targeting of conserved sites, select RBD bnAbs exhibited protective efficacy against diverse SARS-like coronaviruses in a prophylaxis challenge model in vivo. These bnAbs provide new opportunities and choices for next-generation antibody prophylactic and therapeutic applications and provide a molecular basis for effective design of pan-sarbecovirus vaccines.

Gene prediction in the immunoglobulin loci

The V(D)J recombination process rearranges the variable (V), diversity (D), and joining (J) genes in the immunoglobulin loci to generate antibody repertoires. Annotation of these loci across various species and predicting the V, D, and J genes (IG genes) is critical for studies of the adaptive immune system. However, since the standard gene finding algorithms are not suitable for predicting IG genes, they have been semi-manually annotated in very few species. We developed the IGDetective algorithm for predicting IG genes and applied it to species with the assembled IG loci. IGDetective generated the first large collection of IG genes across many species and enabled their evolutionary analysis, including the analysis of the "bat IG diversity" hypothesis. This analysis revealed extremely conserved V genes in evolutionary distant species indicating that these genes may be subjected to the same selective pressure, e.g., pressure driven by common pathogens. IGDetective also revealed extremely diverged V genes and a new family of evolutionary conserved V genes in bats with unusual noncanonical cysteines. Moreover, in difference from all other previously reported antibodies, these cysteines are located within complementarity-determining regions. Since cysteines form disulfide bonds, we hypothesize that these cysteine-rich V genes might generate antibodies with noncanonical conformations and could potentially form a unique part of the immune repertoire in bats. We also analyzed the diversity landscape of the recombination signal sequences and revealed their features that trigger the high/low usage of the IG genes.

Novel Viral DNA Polymerases From Metagenomes Suggest Genomic Sources of Strand-Displacing Biochemical Phenotypes

Viruses are the most abundant and diverse biological entities on the planet and constitute a significant proportion of Earth's genetic diversity. Most of this diversity is not represented by isolated viral-host systems and has only been observed through sequencing of viral metagenomes (viromes) from environmental samples. Viromes provide snapshots of viral genetic potential, and a wealth of information on viral community ecology. These data also provide opportunities for exploring the biochemistry of novel viral enzymes. The in vitro biochemical characteristics of novel viral DNA polymerases were explored, testing hypothesized differences in polymerase biochemistry according to protein sequence phylogeny. Forty-eight viral DNA Polymerase I (PolA) proteins from estuarine viromes, hot spring metagenomes, and reference viruses, encompassing a broad representation of currently known diversity, were synthesized, expressed, and purified. Novel functionality was shown in multiple PolAs. Intriguingly, some of the estuarine viral polymerases demonstrated moderate to strong innate DNA strand displacement activity at high enzyme concentration. Strand-displacing polymerases have important technological applications where isothermal reactions are desirable. Bioinformatic investigation of genes neighboring these strand displacing polymerases found associations with SNF2 helicase-associated proteins. The specific function of SNF2 family enzymes is unknown for prokaryotes and viruses. In eukaryotes, SNF2 enzymes have chromatin remodeling functions but do not separate nucleic acid strands. This suggests the strand separation function may be fulfilled by the DNA polymerase for viruses carrying SNF2 helicase-associated proteins. Biochemical data elucidated from this study expands understanding of the biology and ecological behavior of unknown viruses. Moreover, given the numerous biotechnological applications of viral DNA polymerases, novel viral polymerases discovered within viromes may be a rich source of biological material for further in vitro DNA amplification advancements.

Broadly neutralizing anti-S2 antibodies protect against all three human betacoronaviruses that cause severe disease

Pan-betacoronavirus neutralizing antibodies may hold the key to developing broadly protective vaccines against coronaviruses that cause severe disease, for anticipating novel pandemic-causing viruses, and to respond more effectively to SARS-CoV-2 variants. The emergence of the Omicron variant of SARS-CoV-2 has illustrated the limitations of solely targeting the receptor binding domain (RBD) of the envelope Spike (S)-protein. Here, we isolated a large panel of broadly neutralizing antibodies (bnAbs) from SARS-CoV-2 recovered-vaccinated donors that target a conserved S2 region in the fusion machinery on betacoronavirus spikes. Select bnAbs show broad in vivo protection against all three pathogenic betacoronaviruses, SARS-CoV-1, SARS-CoV-2 and MERS-CoV, that have spilled over into humans in the past 20 years to cause severe disease. The bnAbs provide new opportunities for antibody-based interventions and key insights for developing pan-betacoronavirus vaccines.

Targeted isolation of panels of diverse human protective broadly neutralizing antibodies against SARS-like viruses

The emergence of current SARS-CoV-2 variants of concern (VOCs) and potential future spillovers of SARS-like coronaviruses into humans pose a major threat to human health and the global economy 1-7 . Development of broadly effective coronavirus vaccines that can mitigate these threats is needed 8, 9 . Notably, several recent studies have revealed that vaccination of recovered COVID-19 donors results in enhanced nAb responses compared to SARS-CoV-2 infection or vaccination alone 10-13 . Here, we utilized a targeted donor selection strategy to isolate a large panel of broadly neutralizing antibodies (bnAbs) to sarbecoviruses from two such donors. Many of the bnAbs are remarkably effective in neutralization against sarbecoviruses that use ACE2 for viral entry and a substantial fraction also show notable binding to non-ACE2-using sarbecoviruses. The bnAbs are equally effective against most SARS-CoV-2 VOCs and many neutralize the Omicron variant. Neutralization breadth is achieved by bnAb binding to epitopes on a relatively conserved face of the receptor binding domain (RBD) as opposed to strain-specific nAbs to the receptor binding site that are commonly elicited in SARS-CoV-2 infection and vaccination 14-18 . Consistent with targeting of conserved sites, select RBD bnAbs exhibited in vivo protective efficacy against diverse SARS-like coronaviruses in a prophylaxis challenge model. The generation of a large panel of potent bnAbs provides new opportunities and choices for next-generation antibody prophylactic and therapeutic applications and, importantly, provides a molecular basis for effective design of pan-sarbecovirus vaccines.

Patient health records and whole viral genomes from an early SARS-CoV-2 outbreak in a Quebec hospital reveal features associated with favorable outcomes

The first confirmed case of COVID-19 in Quebec, Canada, occurred at Verdun Hospital on February 25, 2020. A month later, a localized outbreak was observed at this hospital. We performed tiled amplicon whole genome nanopore sequencing on nasopharyngeal swabs from all SARS-CoV-2 positive samples from 31 March to 17 April 2020 in 2 local hospitals to assess viral diversity (unknown at the time in Quebec) and potential associations with clinical outcomes. We report 264 viral genomes from 242 individuals-both staff and patients-with associated clinical features and outcomes, as well as longitudinal samples and technical replicates. Viral lineage assessment identified multiple subclades in both hospitals, with a predominant subclade in the Verdun outbreak, indicative of hospital-acquired transmission. Dimensionality reduction identified two subclades with mutations of clinical interest, namely in the Spike protein, that evaded supervised lineage assignment methods-including Pangolin and NextClade supervised lineage assignment tools. We also report that certain symptoms (headache, myalgia and sore throat) are significantly associated with favorable patient outcomes. Our findings demonstrate the strength of unsupervised, data-driven analyses whilst suggesting that caution should be used when employing supervised genomic workflows, particularly during the early stages of a pandemic.

A comprehensive in silico investigation into the nsSNPs of Drd2 gene predicts significant functional consequences in dopamine signaling and pharmacotherapy

DRD2 is a neuronal cell surface protein involved in brain development and function. Variations in the Drd2 gene have clinical significance since DRD2 is a pharmacotherapeutic target for treating psychiatric disorders like ADHD and schizophrenia. Despite numerous studies on the disease association of single nucleotide polymorphisms (SNPs) in the intronic regions, investigation into the coding regions is surprisingly limited. In this study, we aimed at identifying potential functionally and pharmaco-therapeutically deleterious non-synonymous SNPs of Drd2. A wide array of bioinformatics tools was used to evaluate the impact of nsSNPs on protein structure and functionality. Out of 260 nsSNPs retrieved from the dbSNP database, initially 9 were predicted as deleterious by 15 tools. Upon further assessment of their domain association, conservation profile, homology models and inter-atomic interaction, the mutant F389V was considered as the most impactful. In-depth analysis of F389V through Molecular Docking and Dynamics Simulation revealed a decline in affinity for its native agonist dopamine and an increase in affinity for the antipsychotic drug risperidone. Remarkable alterations in binding interactions and stability of the protein-ligand complex in simulated physiological conditions were also noted. These findings will improve our understanding of the consequence of nsSNPs in disease-susceptibility and therapeutic efficacy.

Micronutrient supplements can promote disruptive protozoan and fungal communities in the developing infant gut

Supplementation with micronutrients, including vitamins, iron and zinc, is a key strategy to alleviate child malnutrition. However, association of gastrointestinal disorders with iron has led to ongoing debate over their administration. To better understand their impact on gut microbiota, we analyse the bacterial, protozoal, fungal and helminth communities of stool samples collected from a subset of 80 children at 12 and 24 months of age, previously enrolled into a large cluster randomized controlled trial of micronutrient supplementation in Pakistan (ClinicalTrials.gov identifier NCT00705445). We show that while bacterial diversity is reduced in supplemented children, vitamins and iron (as well as residence in a rural setting) may promote colonization with distinct protozoa and mucormycetes, whereas the addition of zinc appears to ameliorate this effect. We suggest that the risks and benefits of micronutrient interventions may depend on eukaryotic communities, potentially exacerbated by exposure to a rural setting. Larger studies are needed to evaluate the clinical significance of these findings and their impact on health outcomes.

NgAgo possesses guided DNA nicking activity

Prokaryotic Argonautes (pAgos) have been proposed as more flexible tools for gene-editing as they do not require sequence motifs adjacent to their targets for function, unlike popular CRISPR/Cas systems. One promising pAgo candidate, from the halophilic archaeon Natronobacterium gregoryi (NgAgo), has been the subject of debate regarding its potential in eukaryotic systems. Here, we revisit this enzyme and characterize its function in prokaryotes. NgAgo expresses poorly in non-halophilic hosts with most of the protein being insoluble and inactive even after refolding. However, we report that the soluble fraction does indeed act as a nicking DNA endonuclease. NgAgo shares canonical domains with other catalytically active pAgos but also contains a previously unrecognized single-stranded DNA binding domain (repA). Both repA and the canonical PIWI domains participate in DNA cleavage activities of NgAgo. NgAgo can be programmed with guides to nick targeted DNA in Escherichia coli and in vitro 1 nt outside the 3' end of the guide sequence. We also found that these endonuclease activities are essential for enhanced NgAgo-guided homologous recombination, or gene-editing, in E. coli. Collectively, our results demonstrate the potential of NgAgo for gene-editing and provide new insight into seemingly contradictory reports.

Aspergillus fumigatus In-Host HOG Pathway Mutation for Cystic Fibrosis Lung Microenvironment Persistence

The prevalence of Aspergillus fumigatus colonization in individuals with cystic fibrosis (CF) and subsequent fungal persistence in the lung is increasingly recognized. However, there is no consensus for clinical management of A. fumigatus in CF individuals, due largely to uncertainty surrounding A. fumigatus CF pathogenesis and virulence mechanisms. To address this gap in knowledge, a longitudinal series of A. fumigatus isolates from an individual with CF were collected over 4.5 years. Isolate genotypes were defined with whole-genome sequencing that revealed both transitory and persistent A. fumigatus in the lung. Persistent lineage isolates grew most readily in a low-oxygen culture environment, and conidia were more sensitive to oxidative stress-inducing conditions than those from nonpersistent isolates. Closely related persistent isolates harbored a unique allele of the high-osmolarity glycerol (HOG) pathway mitogen-activated protein kinase kinase, Pbs2 (pbs2C2). Data suggest this novel pbs2C2 allele arose in vivo and is necessary for the fungal response to osmotic stress in a low-oxygen environment through hyperactivation of the HOG (SakA) signaling pathway. Hyperactivation of the HOG pathway through pbs2C2 comes at the cost of decreased conidial stress resistance in the presence of atmospheric oxygen levels. These novel findings shed light on pathoadaptive mechanisms of A. fumigatus in CF, lay the foundation for identifying persistent A. fumigatus isolates that may require antifungal therapy, and highlight considerations for successful culture of persistent Aspergillus CF isolates. IMPORTANCE Aspergillus fumigatus infection causes a spectrum of clinical manifestations. For individuals with cystic fibrosis (CF), allergic bronchopulmonary aspergillosis (ABPA) is an established complication, but there is a growing appreciation for A. fumigatus airway persistence in CF disease progression. There currently is little consensus for clinical management of A. fumigatus long-term culture positivity in CF. A better understanding of A. fumigatus pathogenesis mechanisms in CF is expected to yield insights into when antifungal therapies are warranted. Here, a 4.5-year longitudinal collection of A. fumigatus isolates from a patient with CF identified a persistent lineage that harbors a unique allele of the Pbs2 mitogen-activated protein kinase kinase (MAPKK) necessary for unique CF-relevant stress phenotypes. Importantly for A. fumigatus CF patient diagnostics, this allele provides increased fitness under CF lung-like conditions at a cost of reduced in vitro growth under standard laboratory conditions. These data illustrate a molecular mechanism for A. fumigatus CF lung persistence with implications for diagnostics and antifungal therapy.

Phylogeny and abiotic conditions shape the diel floral emission patterns of desert Brassicaceae species

A key facet of floral scent is diel fluctuations in emission, often studied in the context of plant-pollinator interactions, while contributions of environment and phylogeny remain overlooked. Here, we ask if these factors are involved in shaping temporal variations in scent emission. To that end, we coupled light/dark floral emission measurements of 17 desert Brassicaceae species with environmental and phylogenetic data to explore the individual/combined impacts of these predictors on diel emission patterns. We further investigated these patterns by conducting high-resolution emission measurements in a subset of genetically distant species with contrasting temporal dynamics. While diel shifts in magnitude and richness of emission were strongly affected by genetic relatedness, they also reflect the environmental conditions under which the species grow. Specifically, light/dark emission ratios were negatively affected by an increase in winter temperatures, known to impact both plant physiology and insect locomotion, and sandy soil fractions, previously shown to exert stress that tempers with diel metabolic rhythms. Additionally, the biosynthetic origins of the compounds were associated with their corresponding production patterns, possibly to maximize emission efficacy. Using a multidisciplinary chemical/ecological approach, we uncover and differentiate the main factors shaping floral scent diel fluctuations, highlighting their consequences under changing global climate.

Iron Oxidation by a Fused Cytochrome-Porin Common to Diverse Iron-Oxidizing Bacteria

Iron (Fe) oxidation is one of Earth’s major biogeochemical processes, key to weathering, soil formation, water quality, and corrosion. However, our understanding of microbial contribution is limited by incomplete knowledge of microbial iron oxidation mechanisms, particularly in neutrophilic iron oxidizers. The genomes of many diverse iron oxidizers encode a homolog to an outer membrane cytochrome (Cyc2) shown to oxidize iron in two acidophiles. Phylogenetic analyses show Cyc2 sequences from neutrophiles cluster together, suggesting a common function, though this function has not been verified in these organisms. Therefore, we investigated the iron oxidase function of heterologously expressed Cyc2 from a neutrophilic iron oxidizer Mariprofundus ferrooxydans PV-1. Cyc2_PV-1 is capable of oxidizing iron, and its redox potential is 208 ± 20 mV, consistent with the ability to accept electrons from Fe²⁺ at neutral pH. These results support the hypothesis that Cyc2 functions as an iron oxidase in neutrophilic iron-oxidizing organisms. The results of sequence analysis and modeling reveal that the entire Cyc2 family shares a unique fused cytochrome-porin structure, with a defining consensus motif in the cytochrome region. On the basis of results from structural analyses, we predict that the monoheme cytochrome Cyc2 specifically oxidizes dissolved Fe²⁺, in contrast to multiheme iron oxidases, which may oxidize solid Fe(II). With our results, there is now functional validation for diverse representatives of Cyc2 sequences. We present a comprehensive Cyc2 phylogenetic tree and offer a roadmap for identifying cyc2/Cyc2 homologs and interpreting their function. The occurrence of cyc2 in many genomes beyond known iron oxidizers presents the possibility that microbial iron oxidation may be a widespread metabolism.

Electrophysiological and Transcriptomic Features Reveal a Circular Taxonomy of Cortical Neurons

The complete understanding of the mammalian brain requires exact knowledge of the function of each neuron subpopulation composing its parts. To achieve this goal, an exhaustive, precise, reproducible, and robust neuronal taxonomy should be defined. In this paper, a new circular taxonomy based on transcriptomic features and novel electrophysiological features is proposed. The approach is validated by analysing more than 1850 electrophysiological signals of different mouse visual cortex neurons proceeding from the Allen Cell Types database. The study is conducted on two different levels: neurons and their cell-type aggregation into Cre lines. At the neuronal level, electrophysiological features have been extracted with a promising model that has already proved its worth in neuronal dynamics. At the Cre line level, electrophysiological and transcriptomic features are joined on cell types with available genetic information. A taxonomy with a circular order is revealed by a simple transformation of the first two principal components that allow the characterization of the different Cre lines. Moreover, the proposed methodology locates other Cre lines in the taxonomy that do not have transcriptomic features available. Finally, the taxonomy is validated by Machine Learning methods which are able to discriminate the different neuron types with the proposed electrophysiological features.

Chasing breeding footprints through structural variations in Cucumis melo and wild relatives

Cucumis melo (melon or muskmelon) is an important crop in the family of the Cucurbitaceae. Melon is cross pollinated and domesticated at several locations throughout the breeding history, resulting in highly diverse genetic structure in the germplasm. Yet, the relations among the groups and cultivars are still incomplete. We shed light on the melonbreeding history, analyzing structural variations ranging from 50 bp up to 100 kb, identified from whole genome sequences of 100 selected melon accessions and wild relatives. Phylogenetic trees based on SV types completely resolve cultivars and wild accessions into two monophyletic groups and clustering of cultivars largely correlates with their geographic origin. Taking into account morphology, we found six mis-categorized cultivars. Unique inversions are more often shared between cultivars, carrying advantageous genes and do not directly originate from wild species. Approximately 60% of the inversion breaks carry a long poly A/T motif, and following observations in other plant species, suggest that inversions in melon likely resulted from meiotic recombination events. We show that resistance genes in the linkage V region are expanded in the cultivar genomes compared to wild relatives. Furthermore, particular agronomic traits such as fruit ripening, fragrance, and stress response are specifically selected for in the melon subspecies. These results represent distinctive footprints of selective breeding that shaped today's melon. The sequences and genomic relations between land races, wild relatives, and cultivars will serve the community to identify genetic diversity, optimize experimental designs, and enhance crop development.

Genomic epidemiology and associated clinical outcomes of a SARS-CoV-2 outbreak in a general adult hospital in Quebec

The first confirmed case of COVID-19 in Quebec, Canada, occurred at Verdun Hospital on February 25, 2020. A month later, a localized outbreak was observed at this hospital. We performed tiled amplicon whole genome nanopore sequencing on nasopharyngeal swabs from all SARS-CoV-2 positive samples from 31 March to 17 April 2020 in 2 local hospitals to assess the viral diversity of the outbreak. We report 264 viral genomes from 242 individuals (both staff and patients) with associated clinical features and outcomes, as well as longitudinal samples, technical replicates and the first publicly disseminated SARS-CoV-2 genomes in Quebec. Viral lineage assessment identified multiple subclades in both hospitals, with a predominant subclade in the Verdun outbreak, indicative of hospital-acquired transmission. Dimensionality reduction identified two subclades that evaded supervised lineage assignment methods, including Pangolin, and identified certain symptoms (headache, myalgia and sore throat) that are significantly associated with favorable patient outcomes. We also address certain limitations of standard SARS-CoV-2 bioinformatics procedures, notably when presented with multiple viral haplotypes.

Functional Diversity of the Litter-Associated Fungi from an Oxalate-Carbonate Pathway Ecosystem in Madagascar

The oxalate-carbonate pathway (OCP) is a biogeochemical process linking oxalate oxidation and carbonate precipitation. Currently, this pathway is described as a tripartite association involving oxalogenic plants, oxalogenic fungi, and oxalotrophic bacteria. While the OCP has recently received increasing interest given its potential for capturing carbon in soils, there are still many unknowns, especially regarding the taxonomic and functional diversity of the fungi involved in this pathway. To fill this gap, we described an active OCP site in Madagascar, under the influence of the oxalogenic tree Tamarindus indica, and isolated, identified, and characterized 50 fungal strains from the leaf litter. The fungal diversity encompassed three phyla, namely Mucoromycota, Ascomycota, and Basidiomycota, and 23 genera. Using various media, we further investigated their functional potential. Most of the fungal strains produced siderophores and presented proteolytic activities. The majority were also able to decompose cellulose and xylan, but only a few were able to solubilize inorganic phosphate. Regarding oxalate metabolism, several strains were able to produce calcium oxalate crystals while others decomposed calcium oxalate. These results challenge the current view of the OCP by indicating that fungi are both oxalate producers and degraders. Moreover, they strengthen the importance of the role of fungi in C, N, Ca, and Fe cycles.

Genetic diversity of porcine circovirus 3 strains and the first detection of two different PCV3 strains coinfecting the same host in Minas Gerais, Brazil

Porcine circovirus 3 (PCV3) is a recently emerged circovirus discovered in 2016 that has drawn the attention of the swine industry worldwide. In this study, we evaluated the genetic diversity of PCV3 strains on pig farms. A total of 261 samples from sows, weaning pigs, growing pigs, and stillborn/mummified fetuses were analyzed by quantitative real-time PCR. The results revealed that at least two main lineages of PCV3 are circulating in Brazil. For the first time, it was possible to detect the presence of two different PCV3 strains in the same host.

Acetan and Acetan-Like Polysaccharides: Genetics, Biosynthesis, Structure, and Viscoelasticity

Bacteria produce a variety of multifunctional polysaccharides, including structural, intracellular, and extracellular polysaccharides. They are attractive for the industrial sector due to their natural origin, sustainability, biodegradability, low toxicity, stability, unique viscoelastic properties, stable cost, and supply. When incorporated into different matrices, they may control emulsification, stabilization, crystallization, water release, and encapsulation. Acetan is an important extracellular water-soluble polysaccharide produced mainly by bacterial species of the genera Komagataeibacter and Acetobacter. Since its original description in Komagataeibacter xylinus, acetan-like polysaccharides have also been described in other species of acetic acid bacteria. Our knowledge on chemical composition of different acetan-like polysaccharides, their viscoelasticity, and the genetic basis for their production has expanded during the last years. Here, we review data on acetan biosynthesis, its molecular structure, genetic organization, and mechanical properties. In addition, we have performed an extended bioinformatic analysis on acetan-like polysaccharide genetic clusters in the genomes of Komagataeibacter and Acetobacter species. The analysis revealed for the first time a second acetan-like polysaccharide genetic cluster, that is widespread in both genera. All species of the Komagataeibacter possess at least one acetan genetic cluster, while it is present in only one third of the Acetobacter species surveyed.

Genome size versus geographic range size in birds

Why do some species occur in small, restricted areas, while others are distributed globally? Environmental heterogeneity increases with area and so does the number of species. Hence, diverse biotic and abiotic conditions across large ranges may lead to specific adaptations that are often linked to a species' genome size and chromosome number. Therefore, a positive association between genome size and geographic range is anticipated. Moreover, high cognitive ability in organisms would be favored by natural selection to cope with the dynamic conditions within large geographic ranges. Here, we tested these hypotheses in birds-the most mobile terrestrial vertebrates-and accounted for the effects of various confounding variables, such as body mass, relative brain mass, and geographic latitude. Using phylogenetic generalized least squares and phylogenetic confirmatory path analysis, we demonstrated that range size is positively associated with bird genome size but probably not with chromosome number. Moreover, relative brain mass had no effect on range size, whereas body mass had a possible weak and negative effect, and range size was larger at higher geographic latitudes. However, our models did not fully explain the overall variation in range size. Hence, natural selection may impose larger genomes in birds with larger geographic ranges, although there may be additional explanations for this phenomenon.

Active-site loop variations adjust activity and selectivity of the cumene dioxygenase

Active-site loops play essential roles in various catalytically important enzyme properties like activity, selectivity, and substrate scope. However, their high flexibility and diversity makes them challenging to incorporate into rational enzyme engineering strategies. Here, we report the engineering of hot-spots in loops of the cumene dioxygenase from Pseudomonas fluorescens IP01 with high impact on activity, regio- and enantioselectivity. Libraries based on alanine scan, sequence alignments, and deletions along with a novel insertion approach result in up to 16-fold increases in activity and the formation of novel products and enantiomers. CAVER analysis suggests possible increases in the active pocket volume and formation of new active-site tunnels, suggesting additional degrees of freedom of the substrate in the pocket. The combination of identified hot-spots with the Linker In Loop Insertion approach proves to be a valuable addition to future loop engineering approaches for enhanced biocatalysts.

Soil Texture, Sampling Depth and Root Hairs Shape the Structure of ACC Deaminase Bacterial Community Composition in Maize Rhizosphere

Preservation of the phytostimulatory functions of plant growth-promoting bacteria relies on the adaptation of their community to the rhizosphere environment. Here, an amplicon sequencing approach was implemented to specifically target microorganisms with 1-aminocyclopropane-1-carboxylate deaminase activity, carrying the acdS gene. We stated the hypothesis that the relative phylogenetic distribution of acdS carrying microorganisms is affected by the presence or absence of root hairs, soil type, and depth. To this end, a standardized soil column experiment was conducted with maize wild type and root hair defective rth3 mutant in the substrates loam and sand, and harvest was implemented from three depths. Most acdS sequences (99%) were affiliated to Actinobacteria and Proteobacteria, and the strongest influence on the relative abundances of sequences were exerted by the substrate. Variovorax, Acidovorax, and Ralstonia sequences dominated in loam, whereas Streptomyces and Agromyces were more abundant in sand. Soil depth caused strong variations in acdS sequence distribution, with differential levels in the relative abundances of acdS sequences affiliated to Tetrasphaera, Amycolatopsis, and Streptomyces in loam, but Burkholderia, Paraburkholderia, and Variovorax in sand. Maize genotype influenced the distribution of acdS sequences mainly in loam and only in the uppermost depth. Variovorax acdS sequences were more abundant in WT, but Streptomyces, Microbacterium, and Modestobacter in rth3 rhizosphere. Substrate and soil depth were strong and plant genotype a further significant single and interacting drivers of acdS carrying microbial community composition in the rhizosphere of maize. This suggests that maize rhizosphere acdS carrying bacterial community establishes according to the environmental constraints, and that root hairs possess a minor but significant impact on acdS carrying bacterial populations.

Non-oral Prevotella stepping into the spotlight

Species now affiliated to genus Prevotella have been known for decades as an integral part of human oral cavity microbiota. They were frequently isolated from patients with periodontitis or from dental root canals but also from healthy subjects. With the exception of Prevotella intermedia, they were considered opportunistic pathogens, as they were isolated also from various bacterial abscesses from the head, neck, breast, skin and various other body sites. Consequently, Prevotella were not in the focus of research activities. On the other hand, the four species found in the rumen never caused any disease and seemed early on to be numerous and important part of the rumen ecosystem indicating this genus harbored bacteria with enormously diverse habitats and lifestyles. The purpose of this review is to illustrate the main research themes performed in Prevotella on a path from less noted oral bacteria and from hard to cultivate and study rumen organisms to important mutualistic bacteria in guts of various mammals warranting major research efforts.

Aerobic and anaerobic iron oxidizers together drive denitrification and carbon cycling at marine iron-rich hydrothermal vents

In principle, iron oxidation can fuel significant primary productivity and nutrient cycling in dark environments such as the deep sea. However, we have an extremely limited understanding of the ecology of iron-based ecosystems, and thus the linkages between iron oxidation, carbon cycling, and nitrate reduction. Here we investigate iron microbial mats from hydrothermal vents at Lōʻihi Seamount, Hawaiʻi, using genome-resolved metagenomics and metatranscriptomics to reconstruct potential microbial roles and interactions. Our results show that the aerobic iron-oxidizing Zetaproteobacteria are the primary producers, concentrated at the oxic mat surface. Their fixed carbon supports heterotrophs deeper in the mat, notably the second most abundant organism, Candidatus Ferristratum sp. (uncultivated gen. nov.) from the uncharacterized DTB120 phylum. Candidatus Ferristratum sp., described using nine high-quality metagenome-assembled genomes with similar distributions of genes, expressed nitrate reduction genes narGH and the iron oxidation gene cyc2 in situ and in response to Fe(II) in a shipboard incubation, suggesting it is an anaerobic nitrate-reducing iron oxidizer. Candidatus Ferristratum sp. lacks a full denitrification pathway, relying on Zetaproteobacteria to remove intermediates like nitrite. Thus, at Lōʻihi, anaerobic iron oxidizers coexist with and are dependent on aerobic iron oxidizers. In total, our work shows how key community members work together to connect iron oxidation with carbon and nitrogen cycling, thus driving the biogeochemistry of exported fluids.