Frequent nonhomologous replacement of replicative helicase loaders by viruses in Vibrionaceae.

Several microbial genomes lack textbook-defined essential genes. If an essential gene is absent from a genome, then an evolutionarily independent gene of unknown function complements its function. Here, we identified frequent nonhomologous replacement of an essential component of DNA replication initiation, a replicative helicase loader gene, in Vibrionaceae. Our analysis of Vibrionaceae genomes revealed two genes with unknown function, named vdhL1 and vdhL2, that were substantially enriched in genomes without the known helicase-loader genes. These genes showed no sequence similarities to genes with known function but encoded proteins structurally similar with a viral helicase loader. Analyses of genomic syntenies and coevolution with helicase genes suggested that vdhL1/2 encodes a helicase loader. The in vitro assay showed that Vibrio harveyi VdhL1 and Vibrio ezurae VdhL2 promote the helicase activity of DnaB. Furthermore, molecular phylogenetics suggested that vdhL1/2 were derived from phages and replaced an intrinsic helicase loader gene of Vibrionaceae over 20 times. This high replacement frequency implies the host's advantage in acquiring a viral helicase loader gene.

The ABCF proteins in Escherichia coli individually cope with 'hard-to-translate' nascent peptide sequences.

Organisms possess a wide variety of proteins with diverse amino acid sequences, and their synthesis relies on the ribosome. Empirical observations have led to the misconception that ribosomes are robust protein factories, but in reality, they have several weaknesses. For instance, ribosomes stall during the translation of the proline-rich sequences, but the elongation factor EF-P assists in synthesizing proteins containing the poly-proline sequences. Thus, living organisms have evolved to expand the translation capability of ribosomes through the acquisition of translation elongation factors. In this study, we have revealed that Escherichia coli ATP-Binding Cassette family-F (ABCF) proteins, YheS, YbiT, EttA and Uup, individually cope with various problematic nascent peptide sequences within the exit tunnel. The correspondence between noncanonical translations and ABCFs was YheS for the translational arrest by nascent SecM, YbiT for poly-basic sequence-dependent stalling and poly-acidic sequence-dependent intrinsic ribosome destabilization (IRD), EttA for IRD at the early stage of elongation, and Uup for poly-proline-dependent stalling. Our results suggest that ATP hydrolysis-coupled structural rearrangement and the interdomain linker sequence are pivotal for handling 'hard-to-translate' nascent peptides. Our study highlights a new aspect of ABCF proteins to reduce the potential risks that are encoded within the nascent peptide sequences.

Piwi-piRNA complexes induce stepwise changes in nuclear architecture at target loci.

PIWI-interacting RNAs (piRNAs) are germline-specific small RNAs that form effector complexes with PIWI proteins (Piwi-piRNA complexes) and play critical roles for preserving genomic integrity by repressing transposable elements (TEs). Drosophila Piwi transcriptionally silences specific targets through heterochromatin formation and increases histone H3K9 methylation (H3K9me3) and histone H1 deposition at these loci, with nuclear RNA export factor variant Nxf2 serving as a co-factor. Using ChEP and DamID-seq, we now uncover a Piwi/Nxf2-dependent target association with nuclear lamins. Hi-C analysis of Piwi or Nxf2-depleted cells reveals decreased intra-TAD and increased inter-TAD interactions in regions harboring Piwi-piRNA target TEs. Using a forced tethering system, we analyze the functional effects of Piwi-piRNA/Nxf2-mediated recruitment of piRNA target regions to the nuclear periphery. Removal of active histone marks is followed by transcriptional silencing, chromatin conformational changes, and H3K9me3 and H1 association. Our data show that the Piwi-piRNA pathway can induce stepwise changes in nuclear architecture and chromatin state at target loci for transcriptional silencing.

Whiskers as hydrodynamic prey sensors in foraging seals.

The darkness of the deep ocean limits the vision of diving predators, except when prey emit bioluminescence. It is hypothesized that deep-diving seals rely on highly developed whiskers to locate their prey. However, if and how seals use their whiskers while foraging in natural conditions remains unknown. We used animal-borne tags to show that free-ranging elephant seals use their whiskers for hydrodynamic prey sensing. Small, cheek-mounted video loggers documented seals actively protracting their whiskers in front of their mouths with rhythmic whisker movement, like terrestrial mammals exploring their environment. Seals focused their sensing effort at deep foraging depths, performing prolonged whisker protraction to detect, pursue, and capture prey. Feeding-event recorders with light sensors demonstrated that bioluminescence contributed to only about 20% of overall foraging success, confirming that whiskers play the primary role in sensing prey. Accordingly, visual prey detection complemented and enhanced prey capture. The whiskers' role highlights an evolutionary alternative to echolocation for adapting to the extreme dark of the deep ocean environment, revealing how sensory abilities shape foraging niche segregation in deep-diving mammals. Mammals typically have mobile facial whiskers, and our study reveals the significant function of whiskers in the natural foraging behavior of a marine predator. We demonstrate the importance of field-based sensory studies incorporating multimodality to better understand how multiple sensory systems are complementary in shaping the foraging success of predators.

Cell surface architecture of the cultivated DPANN archaeon Nanobdella aerobiophila.

The DPANN archaeal clade includes obligately ectosymbiotic species. Their cell surfaces potentially play an important role in the symbiotic interaction between the ectosymbionts and their hosts. However, little is known about the mechanism of ectosymbiosis. Here, we show cell surface structures of the cultivated DPANN archaeon Nanobdella aerobiophila strain MJ1T and its host Metallosphaera sedula strain MJ1HA, using a variety of electron microscopy techniques, i.e., negative-staining transmission electron microscopy, quick-freeze deep-etch TEM, and 3D electron tomography. The thickness, unit size, and lattice symmetry of the S-layer of strain MJ1T were different from those of the host archaeon strain MJ1HA. Genomic and transcriptomic analyses highlighted the most highly expressed MJ1T gene for a putative S-layer protein with multiple glycosylation sites and immunoglobulin-like folds, which has no sequence homology to known S-layer proteins. In addition, genes for putative pectin lyase- or lectin-like extracellular proteins, which are potentially involved in symbiotic interaction, were found in the MJ1T genome based on in silico 3D protein structure prediction. Live cell imaging at the optimum growth temperature of 65°C indicated that cell complexes of strains MJ1T and MJ1HA were motile, but sole MJ1T cells were not. Taken together, we propose a model of the symbiotic interaction and cell cycle of Nanobdella aerobiophila.IMPORTANCEDPANN archaea are widely distributed in a variety of natural and artificial environments and may play a considerable role in the microbial ecosystem. All of the cultivated DPANN archaea so far need host organisms for their growth, i.e., obligately ectosymbiotic. However, the mechanism of the ectosymbiosis by DPANN archaea is largely unknown. To this end, we performed a comprehensive analysis of the cultivated DPANN archaeon, Nanobdella aerobiophila, using electron microscopy, live cell imaging, transcriptomics, and genomics, including 3D protein structure prediction. Based on the results, we propose a reasonable model of the symbiotic interaction and cell cycle of Nanobdella aerobiophila, which will enhance our understanding of the enigmatic physiology and ecological significance of DPANN archaea.

MitoFish, MitoAnnotator, and MiFish Pipeline: Updates in 10 Years.

MitoFish, MitoAnnotator, and MiFish Pipeline are comprehensive databases of fish mitochondrial genomes (mitogenomes), accurate annotation software of fish mitogenomes, and a web platform for metabarcoding analysis of fish mitochondrial environmental DNA (eDNA), respectively. The MitoFish Suite currently receives over 48,000 visits worldwide every year; however, the performance and usefulness of the online platforms can still be improved. Here, we present essential updates on these platforms, including an enrichment of the reference data sets, an enhanced searching function, substantially faster genome annotation and eDNA analysis with the denoising of sequencing errors, and a multisample comparative analysis function. These updates have made our platform more intuitive, effective, and reliable. These updated platforms are freely available at http://mitofish.aori.u-tokyo.ac.jp/.

The Quest for Orthologs orthology benchmark service in 2022.

The Orthology Benchmark Service (https://orthology.benchmarkservice.org) is the gold standard for orthology inference evaluation, supported and maintained by the Quest for Orthologs consortium. It is an essential resource to compare existing and new methods of orthology inference (the bedrock for many comparative genomics and phylogenetic analysis) over a standard dataset and through common procedures. The Quest for Orthologs Consortium is dedicated to maintaining the resource up to date, through regular updates of the Reference Proteomes and increasingly accessible data through the OpenEBench platform. For this update, we have added a new benchmark based on curated orthology assertion from the Vertebrate Gene Nomenclature Committee, and provided an example meta-analysis of the public predictions present on the platform.

AlphaCutter: Efficient removal of non-globular regions from predicted protein structures.

A huge number of high-quality predicted protein structures are now publicly available. However, many of these structures contain non-globular regions, which diminish the performance of downstream structural bioinformatic applications. In this study, we develop AlphaCutter for the removal of non-globular regions from predicted protein structures. A large-scale cleaning of 542,380 predicted SwissProt structures highlights that AlphaCutter is able to (1) remove non-globular regions that are undetectable using pLDDT scores and (2) preserve high integrity of the cleaned domain regions. As useful applications, AlphaCutter improved the folding energy scores and sequence recovery rates in the re-design of domain regions. On average, AlphaCutter takes less than 3 s to clean a protein structure, enabling efficient cleaning of the exploding number of predicted protein structures. AlphaCutter is available at https://github.com/johnnytam100/AlphaCutter. AlphaCutter-cleaned SwissProt structures are available for download at https://doi.org/10.5281/zenodo.7944483.

Chromosome-level genome assemblies of Cutaneotrichosporon spp. (Trichosporonales, Basidiomycota) reveal imbalanced evolution between nucleotide sequences and chromosome synteny.

BACKGROUND: Since DNA information was first used in taxonomy, barcode sequences such as the internal transcribed spacer (ITS) region have greatly aided fungal identification; however, a barcode sequence alone is often insufficient. Thus, multi-gene- or whole-genome-based methods were developed. We previously isolated Basidiomycota yeasts classified in the Trichosporonales. Some strains were described as Cutaneotrichosporon cavernicola and C. spelunceum, whereas strain HIS471 remained unidentified. We analysed the genomes of these strains to elucidate their taxonomic relationship and genetic diversity. RESULTS: The long-read-based assembly resulted in chromosome-level draft genomes consisting of seven chromosomes and one mitochondrial genome. The genome of strain HIS471 has more than ten chromosome inversions or translocations compared to the type strain of C. cavernicola despite sharing identical ITS barcode sequences and displaying an average nucleotide identity (ANI) above 93%. Also, the chromosome synteny between C. cavernicola and the related species, C. spelunceum, showed significant rearrangements, whereas the ITS sequence identity exceeds 98.6% and the ANI is approximately 82%. Our results indicate that the relative evolutionary rates of barcode sequences, whole-genome nucleotide sequences, and chromosome synteny in Cutaneotrichosporon significantly differ from those in the model yeast Saccharomyces. CONCLUSIONS: Our results revealed that the relative evolutionary rates of nucleotide sequences and chromosome synteny are different among fungal clades, likely because different clades have diverse mutation/repair rates and distinct selection pressures on their genomic sequences and syntenic structures. Because diverse syntenic structures can be a barrier to meiotic recombination and may lead to speciation, the non-linear relationships between nucleotide and synteny diversification indicate that sequence-level distances at the barcode or whole-genome level are not sufficient for delineating species boundaries.

Origination of LTR Retroelement-Derived NYNRIN Coincides with Therian Placental Emergence.

The emergence of the placenta is a revolutionary event in the evolution of therian mammals, to which some LTR retroelement-derived genes, such as PEG10, RTL1, and syncytin, are known to contribute. However, therian genomes contain many more LTR retroelement-derived genes that may also have contributed to placental evolution. We conducted large-scale evolutionary genomic and transcriptomic analyses to comprehensively search for LTR retroelement-derived genes whose origination coincided with therian placental emergence and that became consistently expressed in therian placentae. We identified NYNRIN as another Ty3/Gypsy LTR retroelement-derived gene likely to contribute to placental emergence in the therian stem lineage. NYNRIN knockdown inhibited the invasion of HTR8/SVneo invasive-type trophoblasts, whereas the knockdown of its nonretroelement-derived homolog KHNYN did not. Functional enrichment analyses suggested that NYNRIN modulates trophoblast invasion by regulating epithelial-mesenchymal transition and extracellular matrix remodeling and that the ubiquitin-proteasome system is responsible for the functional differences between NYNRIN and KHNYN. These findings extend our knowledge of the roles of LTR retroelement-derived genes in the evolution of therian mammals.

Mirage 2.0: fast and memory-efficient reconstruction of gene-content evolution considering heterogeneous evolutionary patterns among gene families.

SUMMARY: We present Mirage 2.0, which accurately estimates gene-content evolutionary history by considering heterogeneous evolutionary patterns among gene families. Notably, we introduce a deterministic pattern mixture model, which makes Mirage substantially faster and more memory-efficient to be applicable to large datasets with thousands of genomes. AVAILABILITY AND IMPLEMENTATION: The source code is freely available at https://github.com/fukunagatsu/Mirage. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Inverse Potts model improves accuracy of phylogenetic profiling.

MOTIVATION: Phylogenetic profiling is a powerful computational method for revealing the functions of function-unknown genes. Although conventional similarity metrics in phylogenetic profiling achieved high prediction accuracy, they have two estimation biases: an evolutionary bias and a spurious correlation bias. While previous studies reduced the evolutionary bias by considering a phylogenetic tree, few studies have analyzed the spurious correlation bias. RESULTS: To reduce the spurious correlation bias, we developed metrics based on the inverse Potts model (IPM) for phylogenetic profiling. We also developed a metric based on both the IPM and a phylogenetic tree. In an empirical dataset analysis, we demonstrated that these IPM-based metrics improved the prediction performance of phylogenetic profiling. In addition, we found that the integration of several metrics, including the IPM-based metrics, had superior performance to a single metric. AVAILABILITY AND IMPLEMENTATION: The source code is freely available at https://github.com/fukunagatsu/Ipm. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Design and Fabrication of a Thin and Micro-Optical Sensor for Rapid Prototyping.

Optical sensing offers several advantages owing to its non-invasiveness and high sensitivity. The miniaturization of optical sensors will mitigate spatial and weight constraints, expanding their applications and extending the principal advantages of optical sensing to different fields, such as healthcare, Internet of Things, artificial intelligence, and other aspects of society. In this study, we present the development of a miniature optical sensor for monitoring thrombi in extracorporeal membrane oxygenation (ECMO). The sensor, based on a complementary metal-oxide semiconductor integrated circuit (CMOS-IC), also serves as a photodiode, amplifier, and light-emitting diode (LED)-mounting substrate. It is sized 3.8 × 4.8 × 0.75 mm3 and provides reflectance spectroscopy at three wavelengths. Based on semiconductor and microelectromechanical system (MEMS) processes, the design of the sensor achieves ultra-compact millimeter size, customizability, prototyping, and scalability for mass production, facilitating the development of miniature optical sensors for a variety of applications.

Long reads and Hi-C sequencing illuminate the two-compartment genome of the model arbuscular mycorrhizal symbiont Rhizophagus irregularis.

Chromosome folding links genome structure with gene function by generating distinct nuclear compartments and topologically associating domains. In mammals, these undergo preferential interactions and regulate gene expression. However, their role in fungal genome biology is unclear. Here, we combine Nanopore (ONT) sequencing with chromatin conformation capture sequencing (Hi-C) to reveal chromosome and epigenetic diversity in a group of obligate plant symbionts: the arbuscular mycorrhizal fungi (AMF). We find that five phylogenetically distinct strains of the model AMF Rhizophagus irregularis carry 33 chromosomes with substantial within-species variability in size, as well as in gene and repeat content. Strain-specific Hi-C contact maps reveal a 'checkerboard' pattern that underline two dominant euchromatin (A) and heterochromatin (B) compartments. Each compartment differs in the level of gene transcription, regulation of candidate effectors and methylation frequencies. The A-compartment is more gene-dense and contains most core genes, while the B-compartment is more repeat-rich and has higher rates of chromosomal rearrangement. While the B-compartment is transcriptionally repressed, it has significantly more secreted proteins and in planta upregulated candidate effectors, suggesting a possible host-induced change in chromosome conformation. Overall, this study provides a fine-scale view into the genome biology and evolution of model plant symbionts, and opens avenues to study the epigenetic mechanisms that modify chromosome folding during host-microbe interactions.

Spatial phylogenetics of Japanese ferns: Patterns, processes, and implications for conservation.

PREMISE: Biodiversity is often only measured with species richness; however, this metric ignores evolutionary history and is not sufficient for making conservation decisions. Here, we characterize multiple facets and drivers of biodiversity to understand how these relate to bioregions and conservation status in the ferns of Japan. METHODS: We compiled a community data set of 1239 grid cells (20 × 20 km each) including 672 taxa based on >300,000 specimen records. We combined the community data with a phylogeny and functional traits to analyze taxonomic, phylogenetic, and functional diversity and modeled biodiversity metrics in response to environmental factors and reproductive mode. Hierarchical clustering was used to delimit bioregions. Conservation status and threats were assessed by comparing the overlap of significantly diverse grid cells with conservation zones and range maps of native Japanese deer. RESULTS: Taxonomic richness was highest at mid-latitudes. Phylogenetic and functional diversity and phylogenetic endemism were highest in small southern islands. Relative phylogenetic and functional diversity were high at high and low latitudes, and low at mid-latitudes. Grid cells were grouped into three (phylogenetic) or four (taxonomic) major bioregions. Temperature and apomixis were identified as drivers of biodiversity patterns. Conservation status was generally high for grid cells with significantly high biodiversity, but the threat due to herbivory by deer was greater for taxonomic richness than other metrics. CONCLUSIONS: Our integrative approach reveals previously undetected patterns and drivers of biodiversity in the ferns of Japan. Future conservation efforts should recognize that threats can vary by biodiversity metric and consider multiple metrics when establishing conservation priorities.

The Quest for Orthologs benchmark service and consensus calls in 2020.

The identification of orthologs-genes in different species which descended from the same gene in their last common ancestor-is a prerequisite for many analyses in comparative genomics and molecular evolution. Numerous algorithms and resources have been conceived to address this problem, but benchmarking and interpreting them is fraught with difficulties (need to compare them on a common input dataset, absence of ground truth, computational cost of calling orthologs). To address this, the Quest for Orthologs consortium maintains a reference set of proteomes and provides a web server for continuous orthology benchmarking (http://orthology.benchmarkservice.org). Furthermore, consensus ortholog calls derived from public benchmark submissions are provided on the Alliance of Genome Resources website, the joint portal of NIH-funded model organism databases.

A distinct lineage of giant viruses brings a rhodopsin photosystem to unicellular marine predators.

Giant viruses are remarkable for their large genomes, often rivaling those of small bacteria, and for having genes thought exclusive to cellular life. Most isolated to date infect nonmarine protists, leaving their strategies and prevalence in marine environments largely unknown. Using eukaryotic single-cell metagenomics in the Pacific, we discovered a Mimiviridae lineage of giant viruses, which infects choanoflagellates, widespread protistan predators related to metazoans. The ChoanoVirus genomes are the largest yet from pelagic ecosystems, with 442 of 862 predicted proteins lacking known homologs. They are enriched in enzymes for modifying organic compounds, including degradation of chitin, an abundant polysaccharide in oceans, and they encode 3 divergent type-1 rhodopsins (VirR) with distinct evolutionary histories from those that capture sunlight in cellular organisms. One (VirRDTS) is similar to the only other putative rhodopsin from a virus (PgV) with a known host (a marine alga). Unlike the algal virus, ChoanoViruses encode the entire pigment biosynthesis pathway and cleavage enzyme for producing the required chromophore, retinal. We demonstrate that the rhodopsin shared by ChoanoViruses and PgV binds retinal and pumps protons. Moreover, our 1.65-Å resolved VirRDTS crystal structure and mutational analyses exposed differences from previously characterized type-1 rhodopsins, all of which come from cellular organisms. Multiple VirR types are present in metagenomes from across surface oceans, where they are correlated with and nearly as abundant as a canonical marker gene from Mimiviridae Our findings indicate that light-dependent energy transfer systems are likely common components of giant viruses of photosynthetic and phagotrophic unicellular marine eukaryotes.

Graph Splitting: A Graph-Based Approach for Superfamily-Scale Phylogenetic Tree Reconstruction.

A protein superfamily contains distantly related proteins that have acquired diverse biological functions through a long evolutionary history. Phylogenetic analysis of the early evolution of protein superfamilies is a key challenge because existing phylogenetic methods show poor performance when protein sequences are too diverged to construct an informative multiple sequence alignment (MSA). Here, we propose the Graph Splitting (GS) method, which rapidly reconstructs a protein superfamily-scale phylogenetic tree using a graph-based approach. Evolutionary simulation showed that the GS method can accurately reconstruct phylogenetic trees and be robust to major problems in phylogenetic estimation, such as biased taxon sampling, heterogeneous evolutionary rates, and long-branch attraction when sequences are substantially diverge. Its application to an empirical data set of the triosephosphate isomerase (TIM)-barrel superfamily suggests rapid evolution of protein-mediated pyrimidine biosynthesis, likely taking place after the RNA world. Furthermore, the GS method can also substantially improve performance of widely used MSA methods by providing accurate guide trees.

A Fast Likelihood Method to Reconstruct and Visualize Ancestral Scenarios.

The reconstruction of ancestral scenarios is widely used to study the evolution of characters along phylogenetic trees. One commonly uses the marginal posterior probabilities of the character states, or the joint reconstruction of the most likely scenario. However, marginal reconstructions provide users with state probabilities, which are difficult to interpret and visualize, whereas joint reconstructions select a unique state for every tree node and thus do not reflect the uncertainty of inferences. We propose a simple and fast approach, which is in between these two extremes. We use decision-theory concepts (namely, the Brier score) to associate each node in the tree to a set of likely states. A unique state is predicted in tree regions with low uncertainty, whereas several states are predicted in uncertain regions, typically around the tree root. To visualize the results, we cluster the neighboring nodes associated with the same states and use graph visualization tools. The method is implemented in the PastML program and web server. The results on simulated data demonstrate the accuracy and robustness of the approach. PastML was applied to the phylogeography of Dengue serotype 2 (DENV2), and the evolution of drug resistances in a large HIV data set. These analyses took a few minutes and provided convincing results. PastML retrieved the main transmission routes of human DENV2 and showed the uncertainty of the human-sylvatic DENV2 geographic origin. With HIV, the results show that resistance mutations mostly emerge independently under treatment pressure, but resistance clusters are found, corresponding to transmissions among untreated patients.

SonicParanoid: fast, accurate and easy orthology inference.

Motivation: Orthology inference constitutes a common base of many genome-based studies, as a pre-requisite for annotating new genomes, finding target genes for biotechnological applications and revealing the evolutionary history of life. Although its importance keeps rising with the ever-growing number of sequenced genomes, existing tools are computationally demanding and difficult to employ. Results: Here, we present SonicParanoid, which is faster than, but comparably accurate to, the well-established tools with a balanced precision-recall trade-off. Furthermore, SonicParanoid substantially relieves the difficulties of orthology inference for those who need to construct and maintain their own genomic datasets. Availability and implementation: SonicParanoid is available with a GNU GPLv3 license on the Python Package Index and BitBucket. Documentation is available at http://iwasakilab.bs.s.u-tokyo.ac.jp/sonicparanoid. Supplementary information: Supplementary data are available at Bioinformatics online.