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.

Identification of gene isoforms and their switching events between male and female embryos of the parthenogenetic crustacean Daphnia magna.

The cladoceran crustacean Daphnia exhibits phenotypic plasticity, a phenomenon that leads to diverse phenotypes from one genome. Alternative usage of gene isoforms has been considered a key gene regulation mechanism for controlling different phenotypes. However, to understand the phenotypic plasticity of Daphnia, gene isoforms have not been comprehensively analyzed. Here we identified 25,654 transcripts derived from the 9710 genes expressed during environmental sex determination of Daphnia magna using the long-read RNA-Seq with PacBio Iso-Seq. We found that 14,924 transcripts were previously unidentified and 5713 genes produced two or more isoforms. By a combination of Illumina short-read RNA-Seq, we detected 824 genes that implemented switching of the highest expressed isoform between females and males. Among the 824 genes, we found isoform switching of an ortholog of CREB-regulated transcription coactivator, a major regulator of carbohydrate metabolism in animals, and a correlation of this switching event with the sexually dimorphic expression of carbohydrate metabolic genes. These results suggest that a comprehensive catalog of isoforms may lead to understanding the molecular basis for environmental sex determination of Daphnia. We also infer the applicability of the full-length isoform analyses to the elucidation of phenotypic plasticity in Daphnia.

Temporal and fertilizer-dependent dynamics of soil bacterial communities in buckwheat fields under long-term management.

This study integrated bacterial community and soil chemicals to characterize the soil ecosystem in an open upland field managed by six controlled fertilizer programs using the minimum amount of pesticides. Amplicon sequencing the 16S rRNA gene revealed that inorganic nitrogen fertilizer and compost altered the diversity and structure of the soil bacterial community throughout buckwheat (Fagopyrum esculentum Moench 'Hitachiakisoba') cultivation. The bacterial community comprised three clusters that contained bacteria that are prevalent in soils fertilized with nitrogen (cluster 1, 340 taxa), without nitrogen and compost (cluster 2, 234 taxa), and with compost-fertilized (cluster 3, 296 taxa). Cluster 2 contained more taxa in Actinobacteriota and less in Acidobacteriota, and cluster 3 contained more taxa in Gemmatimonadota compared with the other clusters. The most frequent taxa in cluster 1 were within the Chloroflexi phylum. The bacterial community structure correlated with soil chemical properties including pH, total organic carbon, SO42-, soluble Ca2+. A co-occurrence network of bacterial taxa and chemicals identified key bacterial groups comprising the center of a community network that determined topology and dynamics of the network. Temporal dynamics of the bacterial community structure indicated that Burkholderiales were associated with buckwheat ripening, indicating plant-bacteria interaction in the ecosystem.

Uncovering microbiomes of the rice phyllosphere using long-read metagenomic sequencing.

The plant microbiome is crucial for plant growth, yet many important questions remain, such as the identification of specific bacterial species in plants, their genetic content, and location of these genes on chromosomes or plasmids. To gain insights into the genetic makeup of the rice-phyllosphere, we perform a metagenomic analysis using long-read sequences. Here, 1.8 Gb reads are assembled into 26,067 contigs including 142 circular sequences. Within these contigs, 669 complete 16S rRNA genes are clustered into 166 bacterial species, 121 of which show low identity (<97%) to defined sequences, suggesting novel species. The circular contigs contain novel chromosomes and a megaplasmid, and most of the smaller circular contigs are defined as novel plasmids or bacteriophages. One circular contig represents the complete chromosome of a difficult-to-culture bacterium Candidatus Saccharibacteria. Our findings demonstrate the efficacy of long-read-based metagenomics for profiling microbial communities and discovering novel sequences in plant-microbiome studies.

Quantification of caffeine in coffee cans using electrochemical measurements, machine learning, and boron-doped diamond electrodes.

Electrochemical measurements, which exhibit high accuracy and sensitivity under low contamination, controlled electrolyte concentration, and pH conditions, have been used in determining various compounds. The electrochemical quantification capability decreases with an increase in the complexity of the measurement object. Therefore, solvent pretreatment and electrolyte addition are crucial in performing electrochemical measurements of specific compounds directly from beverages owing to the poor measurement quality caused by unspecified noise signals from foreign substances and unstable electrolyte concentrations. To prevent such signal disturbances from affecting quantitative analysis, spectral data of voltage-current values from electrochemical measurements must be used for principal component analysis (PCA). Moreover, this method enables highly accurate quantification even though numerical data alone are challenging to analyze. This study utilized boron-doped diamond (BDD) single-chip electrochemical detection to quantify caffeine content in commercial beverages without dilution. By applying PCA, we integrated electrochemical signals with known caffeine contents and subsequently utilized principal component regression to predict the caffeine content in unknown beverages. Consequently, we addressed existing research problems, such as the high quantification cost and the long measurement time required to obtain results after quantification. The average prediction accuracy was 93.8% compared to the actual content values. Electrochemical measurements are helpful in medical care and indirectly support our lives.

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.

Seq2Phase: language model-based accurate prediction of client proteins in liquid-liquid phase separation.

Motivation: Liquid-liquid phase separation (LLPS) enables compartmentalization in cells without biological membranes. LLPS plays essential roles in membraneless organelles such as nucleoli and p-bodies, helps regulate cellular physiology, and is linked to amyloid formation. Two types of proteins, scaffolds and clients, are involved in LLPS. However, computational methods for predicting LLPS client proteins from amino-acid sequences remain underdeveloped. Results: Here, we present Seq2Phase, an accurate predictor of LLPS client proteins. Information-rich features are extracted from amino-acid sequences by a deep-learning technique, Transformer, and fed into supervised machine learning. Predicted client proteins contained known LLPS regulators and showed localization enrichment into membraneless organelles, confirming the validity of the prediction. Feature analysis revealed that scaffolds and clients have different sequence properties and that textbook knowledge of LLPS-related proteins is biased and incomplete. Seq2Phase achieved high accuracies across human, mouse, yeast, and plant, showing that the method is not overfitted to specific species and has broad applicability. We predict that more than hundreds or thousands of LLPS client proteins remain undiscovered in each species and that Seq2Phase will advance our understanding of still enigmatic molecular and physiological bases of LLPS as well as its roles in disease. Availability and implementation: The software codes in Python underlying this article are available at https://github.com/IwasakiLab/Seq2Phase.

Bacteria can maintain rRNA operons solely on plasmids for hundreds of millions of years.

It is generally assumed that all bacteria must have at least one rRNA operon (rrn operon) on the chromosome, but some strains of the genera Aureimonas and Oecophyllibacter carry their sole rrn operon on a plasmid. However, other related strains and species have chromosomal rrn loci, suggesting that the exclusive presence of rrn operons on a plasmid is rare and unlikely to be stably maintained over long evolutionary periods. Here, we report the results of a systematic search for additional bacteria without chromosomal rrn operons. We find that at least four bacterial clades in the phyla Bacteroidota, Spirochaetota, and Pseudomonadota (Proteobacteria) lost chromosomal rrn operons independently. Remarkably, Persicobacteraceae have apparently maintained this peculiar genome organization for hundreds of millions of years. In our study, all the rrn-carrying plasmids in bacteria lacking chromosomal rrn loci possess replication initiator genes of the Rep_3 family. Furthermore, the lack of chromosomal rrn operons is associated with differences in copy numbers of rrn operons, plasmids, and chromosomal tRNA genes. Thus, our findings indicate that the absence of rrn loci in bacterial chromosomes can be stably maintained over long evolutionary periods.

Arbuscular mycorrhizal fungi heterokaryons have two nuclear populations with distinct roles in host-plant interactions.

Arbuscular mycorrhizal fungi (AMF) are prominent root symbionts that can carry thousands of nuclei deriving from two parental strains in a large syncytium. These co-existing genomes can also vary in abundance with changing environmental conditions. Here we assemble the nuclear genomes of all four publicly available AMF heterokaryons using PacBio high-fidelity and Hi-C sequencing. We find that the two co-existing genomes of these strains are phylogenetically related but differ in structure, content and epigenetics. We confirm that AMF heterokaryon genomes vary in relative abundance across conditions and show this can lead to nucleus-specific differences in expression during interactions with plants. Population analyses also reveal signatures of genetic exchange indicative of past events of sexual reproduction in these strains. This work uncovers the origin and contribution of two nuclear genomes in AMF heterokaryons and opens avenues for the improvement and environmental application of these strains.

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.

A fungal sesquiterpene biosynthesis gene cluster critical for mutualist-pathogen transition in Colletotrichum tofieldiae.

Plant-associated fungi show diverse lifestyles from pathogenic to mutualistic to the host; however, the principles and mechanisms through which they shift the lifestyles require elucidation. The root fungus Colletotrichum tofieldiae (Ct) promotes Arabidopsis thaliana growth under phosphate limiting conditions. Here we describe a Ct strain, designated Ct3, that severely inhibits plant growth. Ct3 pathogenesis occurs through activation of host abscisic acid pathways via a fungal secondary metabolism gene cluster related to the biosynthesis of sesquiterpene metabolites, including botrydial. Cluster activation during root infection suppresses host nutrient uptake-related genes and changes mineral contents, suggesting a role in manipulating host nutrition state. Conversely, disruption or environmental suppression of the cluster renders Ct3 beneficial for plant growth, in a manner dependent on host phosphate starvation response regulators. Our findings indicate that a fungal metabolism cluster provides a means by which infectious fungi modulate lifestyles along the parasitic-mutualistic continuum in fluctuating environments.

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.

Chromosome genome assembly and annotation of the spiny red gurnard (Chelidonichthys spinosus).

Chelidonichthys spinosus, a secondary economic fish, is increasingly being exploited and valued in China. However, overfishing has led to it being recognized as one of the most depleted marine species in China. In this study, we generated a chromosome-level genome of C. spinosus using PacBio, Illumina, and Hi-C sequencing data. Ultimately, we assembled a 624.7 Mb genome of C. spinosus, with a contig N50 of 13.77 Mb and scaffold N50 of 28.11 Mb. We further anchored and oriented the assembled sequences onto 24 pseudo-chromosomes using Hi-C techniques. In total, 25,358 protein-coding genes were predicted, of which 24,072 (94.93%) genes were functionally annotated. The dot plot reveals a prominent co-linearity between C. spinosus and Cyclopterus lumpus, indicating a remarkably close phylogenetic relationship between these two species. The assembled genome sequences provide valuable information for elucidating the genetic adaptation and potential molecular basis of C. spinosus. They also have the potential to provide insight into the evolutionary investigation of teleost fish and vertebrates.

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.

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/.

Machine learning enables prediction of metabolic system evolution in bacteria.

Evolution prediction is a long-standing goal in evolutionary biology, with potential impacts on strategic pathogen control, genome engineering, and synthetic biology. While laboratory evolution studies have shown the predictability of short-term and sequence-level evolution, that of long-term and system-level evolution has not been systematically examined. Here, we show that the gene content evolution of metabolic systems is generally predictable by applying ancestral gene content reconstruction and machine learning techniques to ~3000 bacterial genomes. Our framework, Evodictor, successfully predicted gene gain and loss evolution at the branches of the reference phylogenetic tree, suggesting that evolutionary pressures and constraints on metabolic systems are universally shared. Investigation of pathway architectures and meta-analysis of metagenomic datasets confirmed that these evolutionary patterns have physiological and ecological bases as functional dependencies among metabolic reactions and bacterial habitat changes. Last, pan-genomic analysis of intraspecies gene content variations proved that even "ongoing" evolution in extant bacterial species is predictable in our framework.

Nondestructive microbial discrimination using single-cell Raman spectra and random forest machine learning algorithm.

Raman microspectroscopy is a powerful tool for obtaining biomolecular information from single microbial cells in a nondestructive manner. Here, we detail steps to discriminate prokaryotic species using single-cell Raman spectra acquisitions followed by data preprocessing and random forest model tuning. In addition, we describe the steps required to evaluate the model. This protocol requires minimal preprocessing of Raman spectral data, making it accessible to non-spectroscopists, yet allows intuitive visualization of feature importance. For complete details on the use and execution of this protocol, please refer to Kanno et al. (2021).