Search

Accurate inference of orthologous genes constitutes a prerequisite for comparative and evolutionary genomics. SonicParanoid is one of the fastest tools for orthology inference; however, its scalability and accuracy have been hampered by time-consuming all-versus-all alignments and the existence of proteins with complex domain architectures. Here, we present a substantial update of SonicParanoid, where a gradient boosting predictor halves the execution time and a language model doubles the recall. Application to empirical large-scale and standardized benchmark datasets shows that SonicParanoid2 is much faster than comparable methods and also the most accurate. SonicParanoid2 is available at https://gitlab.com/salvo981/sonicparanoid2 and https://zenodo.org/doi/10.5281/zenodo.11371108 ., (© 2024. The Author(s).)

Convergence provides clues to unveil the non-random nature of evolution. Intermediate paths toward convergence inform us of the stochasticity and the constraint of evolutionary processes. Although previous studies have suggested that substantial constraints exist in microevolutionary paths, it remains unclear whether macroevolutionary convergence follows stochastic or constrained paths. Here, we performed comparative genomics for hundreds of lactic acid bacteria (LAB) species, including clades showing a convergent gene repertoire and sharing fructose-rich habitats. By adopting phylogenetic comparative methods we showed that the genomic convergence of distinct fructophilic LAB (FLAB) lineages was caused by parallel losses of more than a hundred orthologs and the gene losses followed significantly similar orders. Our results further suggested that the loss of adhE, a key gene for phenotypic convergence to FLAB, follows a specific evolutionary path of domain architecture decay and amino acid substitutions in multiple LAB lineages sharing fructose-rich habitats. These findings unveiled the constrained evolutionary paths toward the convergence of free-living bacterial clades at the genomic and molecular levels., (© 2024. The Author(s).)

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., (© The Author(s) 2024. Published by Oxford University Press on behalf of Nucleic Acids Research.)

Soil plays a crucial role in human health through its impact on food and habitation. However, it often contains toxic heavy metals, with mercury being particularly hazardous when methylated. Currently, high-sensitivity, rapid detection of mercury is achievable only through electrochemical measurements. These measurements require pretreatment of the soil sample and the preparation of a calibration curve tailored to the sample's condition. In this study, we developed a method to determine the environmental standard value of mercury content in soil by significantly reducing the pretreatment process. Our approach involves analyzing current peaks from electrodeposition times using specific electrodes and solvent settings. This method demonstrates low error rates under low concentration conditions and can detect mercury levels as low as 0.5 ppb in soil leachate and reagent dilution series. This research facilitates the determination of low mercury concentrations in solutions containing various soil micro-compounds without the need for calibration curves.

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, SO 4 2- , soluble Ca 2+ . 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., (© 2024. The Author(s).)

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., (© 2024. The Author(s).)

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., (© 2024. The Author(s).)

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., Competing Interests: The authors have declared that no competing interests exist., (Copyright: © 2024 Honda et al. This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.)

RecA protein and RecA/Rad51 orthologues are required for homologous recombination and DNA repair in all living creatures. RecA/Rad51 catalyzes formation of the D-loop, an obligatory recombination intermediate, through an ATP-dependent reaction consisting of two phases: homology recognition between double-stranded (ds)DNA and single-stranded (ss)DNA to form a hybrid-duplex core of 6-8 base pairs and subsequent hybrid-duplex/D-loop processing. How dsDNA recognizes homologous ssDNA is controversial. The aromatic residue at the tip of the β-hairpin loop (L2) was shown to stabilize dsDNA-strand separation. We tested a model in which dsDNA strands were separated by the aromatic residue before homology recognition and found that the aromatic residue was not essential to homology recognition, but was required for D-loop processing. Contrary to the model, we found that the double helix was not unwound even a single turn during search for sequence homology, but rather was unwound only after the homologous sequence was recognized. These results suggest that dsDNA recognizes its homologous ssDNA before strand separation. The search for homologous sequence with homologous ssDNA without dsDNA-strand separation does not generate stress within the dsDNA; this would be an advantage for dsDNA to express homology-dependent functions in vivo and also in vitro., (© The Author(s) 2024. Published by Oxford University Press on behalf of Nucleic Acids Research.)

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 MJ1 T 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 MJ1 T were different from those of the host archaeon strain MJ1HA. Genomic and transcriptomic analyses highlighted the most highly expressed MJ1 T 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 MJ1 T and MJ1HA were motile, but sole MJ1 T 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., Competing Interests: The authors declare no conflict of interest.

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., Competing Interests: W.I. is an Editorial Board member of Bioinformatics Advances, but was not involved in the editorial process of this manuscript., (© The Author(s) 2023. Published by Oxford University Press.)

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&#95;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., (© 2023. The Author(s).)

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., (© 2023. BioMed Central Ltd., part of Springer Nature.)

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., (© 2023. Springer Nature Limited.)

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 mm 3 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.

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., (© 2023. The Author(s).)

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 Author(s) 2023. Published by Oxford University Press on behalf of Society for Molecular Biology and Evolution.)

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.

Metagenome-assembled genomes (MAGs) have revealed the hidden diversity and functions of uncultivated microbes, but their reconstruction from metagenomes remains a computationally difficult task. Repetitive or exogenous sequences, such as ribosomal RNA and horizontally transferred genes, are frequently absent from MAGs because of misassembly and binning errors. Here, we report that ribosomal protein genes are also often absent from MAGs, although they are neither repetitive nor exogenous. Comprehensive analyses of more than 190,000 MAGs revealed that these genes could be missing in more than 20-40% of near-complete (i.e., with completeness of 90% or higher) MAGs. While some uncultivated environmental microbes intrinsically lack some ribosomal protein genes, we found that this unexpected absence is largely due to special evolutionary patterns of codon usage bias in ribosomal protein genes and algorithmic characteristics of metagenomic binning, which is dependent on tetranucleotide frequencies of contigs. This problem reflects the microbial life-history strategy. Fast-growing microbes tend to have this difficulty, likely because of strong evolutionary pressures on ribosomal protein genes toward the efficient assembly of ribosomes. Our observations caution those who study genomics and phylogeny of uncultivated microbes, the diversity and evolution of microbial genes in the central dogma, and bioinformatics in metagenomics., (© 2022. The Author(s).)

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)., Competing Interests: The authors declare no competing interests., (© 2022 The Author(s).)

In our previous study, enterohepatic 12α-hydroxylated (12α) bile acid (BA) levels were found to be correlated with hepatic triacylglycerol concentration in rats fed high-fat (HF) diet. Since BA composition is diverse depending on animal species, we evaluated whether such a relationship is applicable in mice in response to an HF diet. C57BL/6JJmsSLC (B6) male mice were fed HF diet for 13 weeks and analyzed for triacylglycerol, cholesterol, oxysterols, and other metabolites in the liver. The BA composition was determined in the liver, small intestinal contents, portal plasma, aortic plasma, and feces. Neutral sterols were also measured in the feces. The ratio of 12α BA/non-12 BA increased in the liver, portal plasma, small intestinal contents, and feces of HF-fed B6 mice. Moreover, a positive correlation was observed between the ratio of fecal 12α BAs/non-12 BAs and hepatic triacylglycerol concentration. The concentration of 7α-hydroxycholesterol was increased in the liver of HF-fed B6 mice, whereas no increase was observed in the hepatic expression of cytochrome P450 family 7 subfamily A member 1. The present study showed that the ratio of 12α BA/non-12 BA in feces is closely associated with hepatic triacylglycerol accumulation in B6 mice fed HF diet., (© 2022. The Author(s).)

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., (© The Author(s) 2022. Published by Oxford University Press on behalf of Society for Molecular Biology and Evolution.)

Ferns, with about 12,000 species, are the second most diverse lineage of vascular plants after angiosperms. They have been the subject of numerous molecular phylogenetic studies, resulting in the publication of trees for every major clade and DNA sequences from nearly half of all species. Global fern phylogenies have been published periodically, but as molecular systematics research continues at a rapid pace, these become quickly outdated. Here, we develop a mostly automated, reproducible, open pipeline to generate a continuously updated fern tree of life (FTOL) from DNA sequence data available in GenBank. Our tailored sampling strategy combines whole plastomes (few taxa, many loci) with commonly sequenced plastid regions (many taxa, few loci) to obtain a global, species-level fern phylogeny with high resolution along the backbone and maximal sampling across the tips. We use a curated reference taxonomy to resolve synonyms in general compliance with the community-driven Pteridophyte Phylogeny Group I classification. The current FTOL includes 5,582 species, an increase of ca. 40% relative to the most recently published global fern phylogeny. Using an updated and expanded list of 51 fern fossil constraints, we find estimated ages for most families and deeper clades to be considerably older than earlier studies. FTOL and its accompanying datasets, including the fossil list and taxonomic database, will be updated on a regular basis and are available via a web portal (https://fernphy.github.io) and R packages, enabling immediate access to the most up-to-date, comprehensively sampled fern phylogeny. FTOL will be useful for anyone studying this important group of plants over a wide range of taxonomic scales, from smaller clades to the entire tree. We anticipate FTOL will be particularly relevant for macroecological studies at regional to global scales and will inform future taxonomic systems with the most recent hypothesis of fern phylogeny., Competing Interests: The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest., (Copyright © 2022 Nitta, Schuettpelz, Ramírez-Barahona and Iwasaki.)

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., (© The Author(s) 2022. Published by Oxford University Press.)

The use of metal deposition has been limited to a limited number of applicable samples due to the increased temperature caused by accelerated electron impact on the substrate surface. The surfaces of various biological samples have a nanoscale structure with specific properties, which have been simulated in numerous studies. However, no examples of nano/microscale reproductions of biological surface features have used moulds. In this study, a mould that imitates the surface shape of a cellular-level biological material was fabricated, for the first time, and the shape was successfully reproduced using the mould. Al thin films were deposited on bovine sperm using magnetron sputtering without thermal denaturation with a cathode operating at a biological temperature. It is difficult to deposit films used as metal coatings on pre-treated biological materials at temperatures below 40 °C during evaporation. The Al thin film was peeled off and used as a mould to reproduce the shape of the sperm with high accuracy using a polymer. The results of this study represent a major innovation in reproducible biomimetic moulding technology, demonstrating biological temperature sputtering. We expect our non-destructive metal deposition and metal nano-moulding methods for biological samples to be the basis for the effective utilization of various biological structures., (© 2022. The Author(s).)

Copy number variation (CNV) can cause phenotypic changes. However, in contrast to amino acid substitutions and cis -regulatory changes, little is known about the functional categories of genes in which CNV is important for adaptation to novel environments. It is also unclear whether the same genes repeatedly change the copy numbers for adapting to similar environments. Here, we investigate CNV associated with freshwater colonization in fishes, which was observed multiple times across different lineages. Using 48 ray-finned fishes across diverse orders, we identified 23 genes whose copy number increases were associated with freshwater colonization. These genes showed enrichment for peptide receptor activity, hexosyltransferase activity and unsaturated fatty acid metabolism. We further revealed that three of the genes showed copy number increases in freshwater populations compared to marine ancestral populations of the stickleback genus Gasterosteus . These results indicate that copy number increases of genes involved in fatty acid metabolism ( FADS2 ), immune function ( PSMB8a ) and thyroid hormone metabolism ( UGT2 ) may be important for freshwater colonization at both the inter-order macroevolutionary scale and at the intra-genus microevolutionary scale. Further analysis across diverse taxa will help to understand the role of CNV in the adaptation to novel environments. This article is part of the theme issue 'Genetic basis of adaptation and speciation: from loci to causative mutations'.

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., (© The Author(s) 2022. Published by Oxford University Press on behalf of Nucleic Acids Research.)

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.

Supporting cells of oocytes, i.e., cumulus cells, control oocyte quality, which determines fertilization success. Therefore, the transformation of mature and immature cumulus cells (MCCs and ICCs, respectively) into dysmature cumulus cells (DCCs) with dead characteristics deteriorates oocyte quality. However, the molecular basis for this transformation remains unclear. Here, we explored the link between autophagic decline and cumulus transformation using cumulus cells from patients with infertility, female mice, and human granulosa cell-derived KGN cell lines. When human cumulus cells were labeled with LysoTracker probes, fluorescence corresponding to lysosomes was enhanced in DCCs compared to that in MCCs and ICCs. Similarly, treatment with the autophagy inhibitor chloroquine elevated LysoTracker fluorescence in both mouse cumulus cells and KGN cells, subsequently suppressing ovulation in female mice. Electron microscopy analysis revealed the proliferation of abnormal lysosomes in chloroquine-treated KGN cells. Conversely, the addition of an autophagy inducer, trehalose, suppressed chloroquine-driven problematic lysosomal anomalies and ameliorated ovulation problems. Our results suggest that autophagy maintains the healthy state of the supporting cells of human oocytes by suppressing the formation of lysosomes. Thus, our results provide insights into the therapeutic effects of trehalose on female fertility.

Owing to its simplicity and sensitivity, electrochemical analysis is of high significance in the detection of pollutants and highly toxic substances in the environment. In electrochemical analysis, the sensitivity of the sensor and reliability of the obtained signal are especially dependent on the electrode characteristics. Electrodes with a high density of nanomaterials, which exhibit excellent activity, are useful as sensor substrates for pollutant detection. However, the effective placement of high-density nanomaterials requires a high degree of control over the particle size, particle shape, and distance between the particles on the substrate. In this study, we exploited the properties of boron-doped diamond (BDD) electrodes, which have a wide potential window, and succeeded in coating a highly dense layer of gold nanoparticles (AuNPs) at high potential. The AuNP-modified BDD (AuNP-BDD) electrodes comprising less than 100 nm AuNPs at a density of 125 particles/µm were electrochemically synthesized over a short period of 30-60 s. The AuNP-BDD electrodes were applied for detecting arsenic, which is one of the most abundant elements, and exhibited a limit of detection of 0.473 ppb in solution.

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., (© 2022 The Authors. American Journal of Botany published by Wiley Periodicals LLC on behalf of Botanical Society of America.)

Fluid control on a paper channel is necessary for analysis with multiple reagents, such as enzyme-linked immunosorbent assay (ELISA) in microfluidic paper-based analytical devices (µPADs). In this study, a thermo-responsive valve was fabricated by polymerizing N-isopropylacrylamide on a PVDF porous membrane by plasma-induced graft polymerization. The polymerized membrane was observed by scanning electron microscopy (SEM), and it was confirmed that more pores were closed at temperatures below 32 °C and more pores were opened at temperatures above 32 °C. Valve permeability tests confirmed that the proposed polymerized membrane was impermeable to water and proteins at temperatures below 32 °C and permeable to water at temperatures above 32 °C. The valve could also be reversibly and repeatedly opened and closed by changing the temperature near 32 °C. These results suggest that plasma-induced graft polymerization may be used to produce thermo-responsive valves that can be opened and closed without subsequent loss of performance. These results indicate that the thermo-responsive valve fabricated by plasma-induced graft polymerization could potentially be applied to ELISA with µPADs.

Phylogeny estimation (the reconstruction of evolutionary trees) has recently been applied to CRISPR-based cell lineage tracing, allowing the developmental history of an individual tissue or organism to be inferred from a large number of mutated sequences in somatic cells. However, current computational methods are not able to construct phylogenetic trees from extremely large numbers of input sequences. Here, we present a deep distributed computing framework to comprehensively trace accurate large lineages (FRACTAL) that substantially enhances the scalability of current lineage estimation software tools. FRACTAL first reconstructs only an upstream lineage of the input sequences and recursively iterates the same produce for its downstream lineages using independent computing nodes. We demonstrate the utility of FRACTAL by reconstructing lineages from >235 million simulated sequences and from >16 million cells from a simulated experiment with a CRISPR system that accumulates mutations during cell proliferation. We also successfully applied FRACTAL to evolutionary tree reconstructions and to an experiment using error-prone PCR (EP-PCR) for large-scale sequence diversification., (© 2022. The Author(s), under exclusive licence to Springer Nature America, Inc.)

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., (© The Author(s) 2022. Published by Oxford University Press.)

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., (© 2021 The Authors. New Phytologist © 2021 New Phytologist Foundation.)

A characteristic clinical feature of COVID-19 is the frequent incidence of microvascular thrombosis. In fact, COVID-19 autopsy reports have shown widespread thrombotic microangiopathy characterized by extensive diffuse microthrombi within peripheral capillaries and arterioles in lungs, hearts, and other organs, resulting in multiorgan failure. However, the underlying process of COVID-19-associated microvascular thrombosis remains elusive due to the lack of tools to statistically examine platelet aggregation (i.e., the initiation of microthrombus formation) in detail. Here we report the landscape of circulating platelet aggregates in COVID-19 obtained by massive single-cell image-based profiling and temporal monitoring of the blood of COVID-19 patients (n = 110). Surprisingly, our analysis of the big image data shows the anomalous presence of excessive platelet aggregates in nearly 90% of all COVID-19 patients. Furthermore, results indicate strong links between the concentration of platelet aggregates and the severity, mortality, respiratory condition, and vascular endothelial dysfunction level of COVID-19 patients., (© 2021. The Author(s).)

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., (© 2021 The Authors. Published under the terms of the CC BY 4.0 license.)

Accessing enormous uncultivated microorganisms (microbial dark matter) in various Earth environments requires accurate, nondestructive classification, and molecular understanding of the microorganisms in in situ and at the single-cell level. Here we demonstrate a combined approach of random forest (RF) machine learning and single-cell Raman microspectroscopy for accurate classification of phylogenetically diverse prokaryotes (three bacterial and three archaeal species from different phyla). Our RF classifier achieved a 98.8 ± 1.9% classification accuracy among the six species in pure populations and 98.4% for three species in an artificially mixed population. Feature importance scores against each wavenumber reveal that the presence of carotenoids and structure of membrane lipids play key roles in distinguishing the prokaryotic species. We also find unique Raman markers for an ammonia-oxidizing archaeon. Our approach with moderate data pretreatment and intuitive visualization of feature importance is easy to use for non-spectroscopists, and thus offers microbiologists a new single-cell tool for shedding light on microbial dark matter., Competing Interests: The authors declare no competing interests., (© 2021 The Author(s).)

Motivation: Reconstruction of gene copy number evolution is an essential approach for understanding how complex biological systems have been organized. Although various models have been proposed for gene copy number evolution, existing evolutionary models have not appropriately addressed the fact that different gene families can have very different gene gain/loss rates., Results: In this study, we developed Mirage (MIxtuRe model for Ancestral Genome Estimation), which allows different gene families to have flexible gene gain/loss rates. Mirage can use three models for formulating heterogeneous evolution among gene families: the discretized Γ model, probability distribution-free model and pattern mixture (PM) model. Simulation analysis showed that Mirage can accurately estimate heterogeneous gene gain/loss rates and reconstruct gene-content evolutionary history. Application to empirical datasets demonstrated that the PM model fits genome data from various taxonomic groups better than the other heterogeneous models. Using Mirage, we revealed that metabolic function-related gene families displayed frequent gene gains and losses in all taxa investigated., Availability and Implementation: The source code of Mirage is freely available at https://github.com/fukunagatsu/Mirage., Supplementary Information: Supplementary data are available at Bioinformatics Advances online., (© The Author(s) 2021. Published by Oxford University Press.)