Proteinortho: detection of (co-)orthologs in large-scale analysis

Evolution of diapause in the African turquoise killifish by remodeling the ancient gene regulatory landscape

Suspended animation states allow organisms to survive extreme environments. The African turquoise killifish has evolved diapause as a form of suspended development to survive a complete drought. However, the mechanisms underlying the evolution of extreme survival states are unknown. To understand diapause evolution, we performed integrative multi-omics (gene expression, chromatin accessibility, and lipidomics) in the embryos of multiple killifish species. We find that diapause evolved by a recent remodeling of regulatory elements at very ancient gene duplicates (paralogs) present in all vertebrates. CRISPR-Cas9-based perturbations identify the transcription factors REST/NRSF and FOXOs as critical for the diapause gene expression program, including genes involved in lipid metabolism. Indeed, diapause shows a distinct lipid profile, with an increase in triglycerides with very-long-chain fatty acids. Our work suggests a mechanism for the evolution of complex adaptations and offers strategies to promote long-term survival by activating suspended animation programs in other species.

The genomic and phylogenetic analysis of Marseillevirus cajuinensis raises questions about the evolution of Marseilleviridae lineages and their taxonomical organization

Marseilleviruses (MsV) are a group of viruses that compose the Marseilleviridae family within the Nucleocytoviricota phylum. They have been found in different samples, mainly in freshwater. MsV are classically organized into five phylogenetic lineages (A/B/C/D/E), but the current taxonomy does not fully represent all the diversity of the MsV lineages. Here, we describe a novel strain isolated from a Brazilian saltwater sample named Marseillevirus cajuinensis. Based on genomics and phylogenetic analyses, M. cajuinensis exhibits a 380,653-bp genome that encodes 515 open reading frames. Additionally, M. cajuinensis encodes a transfer RNA, a feature that is rarely described for Marseilleviridae. Phylogeny suggests that M. cajuinensis forms a divergent branch within the MsV lineage A. Furthermore, our analysis suggests that the common ancestor for the five classical lineages of MsV diversified into three major groups. The organization of MsV into three main groups is reinforced by a comprehensive analysis of clusters of orthologous groups, sequence identities, and evolutionary distances considering several MsV isolates. Taken together, our results highlight the importance of discovering new viruses to expand the knowledge about known viruses that belong to the same lineages or families. This work proposes a new perspective on the Marseilleviridae lineages organization that could be helpful to a future update in the taxonomy of the Marseilleviridae family.

IMPORTANCE

Marseilleviridae is a family of viruses whose members were mostly isolated from freshwater samples. In this work, we describe the first Marseillevirus isolated from saltwater samples, which we called Marseillevirus cajuinensis. Most of M. cajuinensis genomic features are comparable to other Marseilleviridae members, such as its high number of unknown proteins. On the other hand, M. cajuinensis encodes a transfer RNA, which is a gene category involved in protein translation that is rarely described in this viral family. Additionally, our phylogenetic analyses suggested the existence of, at least, three major Marseilleviridae groups. These observations provide a new perspective on Marseilleviridae lineages organization, which will be valuable in future updates to the taxonomy of the family since the current official classification does not capture all the Marseilleviridae known diversity.

On the evolution of chromosomal regions with high gene strand bias in bacteria

On circular bacterial chromosomes, the majority of genes are coded on the leading strand. This gene strand bias (GSB) ranges from up to 85% in some Bacillota to a little more than 50% in other phyla. The factors determining the extent of the strand bias remain to be found. Here, we report that species in the phylum Gemmatimonadota share a unique chromosome architecture, distinct from neighboring phyla: in a conserved 600-kb region around the terminus of replication, almost all genes were located on the leading strands, while on the remaining part of the chromosome, the strand preference was more balanced. The high strand bias (HSB) region harbors the rRNA clusters, core, and highly expressed genes. Selective pressure for reduction of collisions with DNA replication to minimize detrimental mutations can explain the conservation of essential genes in this region. Repetitive and mobile elements are underrepresented, suggesting reduced recombination frequency by structural isolation from other parts of the chromosome. We propose that the HSB region forms a distinct chromosomal domain. Gemmatimonadota chromosomes evolved mainly by expansion through horizontal gene transfer and duplications outside of the ancient high strand bias region. In support of our hypothesis, we could further identify two Spiroplasma strains on a similar evolutionary path.IMPORTANCEOn bacterial chromosomes, a preferred location of genes on the leading strand has evolved to reduce conflicts between replication and transcription. Despite a vast body of research, the question why bacteria show large differences in their gene strand bias is still not solved. The discovery of "hybrid" chromosomes in different phyla, including Gemmatimonadota, in which a conserved high strand bias is found exclusively in a region at ter, points toward a role of nucleoid structure, additional to replication, in the evolution of strand preferences. A fine-grained structural analysis of the ever-increasing number of available bacterial genomes could help to better understand the forces that shape the sequential and spatial organization of the cell's information content.

Widespread occurrence and diverse origins of polintoviruses influence lineage-specific genome dynamics in stony corals

Stony corals (Order: Scleractinia) are central to vital marine habitats known as coral reefs. Numerous stressors in the Anthropocene are contributing to the ongoing decline in coral reef health and coverage. While viruses are established modulators of marine microbial dynamics, their interactions within the coral holobiont and impact on coral health and physiology remain unclear. To address this key knowledge gap, we investigated diverse stony coral genomes for 'endogenous' viruses. Our study uncovered a remarkable number of integrated viral elements recognized as 'Polintoviruses' (Class Polintoviricetes) in thirty Scleractinia genomes; with several species harboring hundreds to thousands of polintoviruses. We reveal massive paralogous expansion of polintoviruses in stony coral genomes, alongside the presence of integrated elements closely related to Polinton-like viruses (PLVs), a group of viruses that exist as free virions. These results suggest multiple integrations of polintoviruses and PLV-relatives, along with paralogous expansions, shaped stony coral genomes. Re-analysis of existing gene expression data reveals all polintovirus structural and non-structural hallmark genes are expressed, providing support for free virion production from polintoviruses. Our results, revealing a significant diversity of polintovirus across the Scleractinia order, open a new research avenue into polintovirus and their possible roles in disease, genomic plasticity, and environmental adaptation in this key group of organisms.

Jeotgalibacillus haloalkalitolerans sp. nov., a novel alkalitolerant and halotolerant bacterium, isolated from the confluence of the Fenhe River and the Yellow River

A Gram-stain positive, aerobic, alkalitolerant and halotolerant bacterium, designated HH7-29 T, was isolated from the confluence of the Fenhe River and the Yellow River in Shanxi Province, PR China. Growth occurred at pH 6.0-12.0 (optimum, pH 8.0-8.5) and 15-40℃ (optimum, 32℃) with 0.5-24% NaCl (optimum, 2-9%). The predominant fatty acids (> 10.0%) were iso-C15:0 and anteiso-C15:0. The major menaquinones were MK-7 and MK-8. The polar lipids were phosphatidylglycerol, diphosphatidylglycerol and two unidentified phospholipids. Phylogenetic analyses based on the 16S rRNA gene sequence revealed that strain HH7-29 T was a member of the genus Jeotgalibacillus, exhibiting high sequence similarity to the 16S rRNA gene sequences of Jeotgalibacillus alkaliphilus JC303T (98.4%), Jeotgalibacillus salarius ASL-1 T (98.1%) and Jeotgalibacillus alimentarius YKJ-13 T (98.1%). The genomic DNA G + C content was 43.0%. Gene annotation showed that strain HH7-29 T had lower protein isoelectric points (pIs) and possessed genes related to ion transport and organic osmoprotectant uptake, implying its potential tolerance to salt and alkali. The average nucleotide identity, digital DNA-DNA hybridization values, amino acid identity values, and percentage of conserved proteins values between strain HH7-29 T and its related species were 71.1-83.8%, 19.5-27.4%, 66.5-88.4% and 59.8-76.6%, respectively. Based on the analyses of phenotypic, chemotaxonomic, phylogenetic and genomic features, strain HH7-29 T represents a novel species of the genus Jeotgalibacillus, for which the name Jeotgalibacillus haloalkalitolerans sp. nov. is proposed. The type strain is HH7-29 T (= KCTC 43417 T = MCCC 1K07541T).

Genome-based taxonomic classification of the genus Sulfitobacter along with the proposal of a new genus Parasulfitobacter gen. nov. and exploring the gene clusters associated with sulfur oxidation

BACKGROUND

The genus Sulfitobacter, a member of the family Roseobacteraceae, is widely distributed in the ocean and is believed to play crucial roles in the global sulfur cycle. However, gene clusters associated with sulfur oxidation in genomes of the type strains of this genus have been poorly studied. Furthermore, taxonomic errors have been identified in this genus, potentially leading to significant confusion in ecological and evolutionary interpretations in subsequent studies of the genus Sulfitobacter. This study aims to investigate the taxonomic status of this genus and explore the metabolism associated with sulfur oxidation.

RESULTS

This study suggests that Sulfitobacter algicola does not belong to Sulfitobacter and should be reclassified into a novel genus, for which we propose the name Parasulfitobacter gen. nov., with Parasulfitobacter algicola comb. nov. as the type species. Additionally, enzymes involved in the sulfur oxidation process, such as the sulfur oxidization (Sox) system, the disulfide reductase protein family, and the sulfite dehydrogenase (SoeABC), were identified in almost all Sulfitobacter species. This finding implies that the majority of Sulfitobacter species can oxidize reduced sulfur compounds. Differences in the modular organization of sox gene clusters among Sulfitobacter species were identified, along with the presence of five genes with unknown function located in some of the sox gene clusters. Lastly, this study revealed the presence of the demethylation pathway and the cleavage pathway used by many Sulfitobacter species to degrade dimethylsulfoniopropionate (DMSP). These pathways enable these bacteria to utilize DMSP as important source of sulfur and carbon or as a defence strategy.

CONCLUSIONS

Our findings contribute to interpreting the mechanism by which Sulfitobacter species participate in the global sulfur cycle. The taxonomic rearrangement of S. algicola into the novel genus Parasulfitobacter will prevent confusion in ecological and evolutionary interpretations in future studies of the genus Sulfitobacter.

Elucidating human gut microbiota interactions that robustly inhibit diverse Clostridioides difficile strains across different nutrient landscapes

The human gut pathogen Clostridioides difficile displays extreme genetic variability and confronts a changeable nutrient landscape in the gut. We mapped gut microbiota inter-species interactions impacting the growth and toxin production of diverse C. difficile strains in different nutrient environments. Although negative interactions impacting C. difficile are prevalent in environments promoting resource competition, they are sparse in an environment containing C. difficile-preferred carbohydrates. C. difficile strains display differences in interactions with Clostridium scindens and the ability to compete for proline. C. difficile toxin production displays substantial community-context dependent variation and does not trend with growth-mediated inter-species interactions. C. difficile shows substantial differences in transcriptional profiles in the presence of the closely related species C. hiranonis or C. scindens. In co-culture with C. hiranonis, C. difficile exhibits massive alterations in metabolism and other cellular processes, consistent with their high metabolic overlap. Further, Clostridium hiranonis inhibits the growth and toxin production of diverse C. difficile strains across different nutrient environments and ameliorates the disease severity of a C. difficile challenge in a murine model. In sum, strain-level variability and nutrient environments are major variables shaping gut microbiota interactions with C. difficile.

Unveiling Gene Expression Dynamics during Early Embryogenesis in Cynoglossus semilaevis: A Transcriptomic Perspective

Commencing with sperm-egg fusion, the early stages of metazoan development include the cleavage and formation of blastula and gastrula. These early embryonic events play a crucial role in ontogeny and are accompanied by a dramatic remodeling of the gene network, particularly encompassing the maternal-to-zygotic transition. Nonetheless, the gene expression dynamics governing early embryogenesis remain unclear in most metazoan lineages. We conducted transcriptomic profiling on two types of gametes (oocytes and sperms) and early embryos (ranging from the four-cell to the gastrula stage) of an economically valuable flatfish-the Chinese tongue sole Cynoglossus semilaevis (Pleuronectiformes: Cynoglossidae). Comparative transcriptome analysis revealed that large-scale zygotic genome activation (ZGA) occurs in the blastula stage, aligning with previous findings in zebrafish. Through the comparison of the most abundant transcripts identified in each sample and the functional analysis of co-expression modules, we unveiled distinct functional enrichments across different gametes/developmental stages: actin- and immune-related functions in sperms; mitosis, transcription inhibition, and mitochondrial function in oocytes and in pre-ZGA embryos (four- to 1000-cell stage); and organ development in post-ZGA embryos (blastula and gastrula). These results provide insights into the intricate transcriptional regulation of early embryonic development in Cynoglossidae fish and expand our knowledge of developmental constraints in vertebrates.

Identification of Parthenogenesis-Inducing Effector Proteins in Wolbachia

Bacteria in the genus Wolbachia have evolved numerous strategies to manipulate arthropod sex, including the conversion of would-be male offspring to asexually reproducing females. This so-called "parthenogenesis induction" phenotype can be found in a number of Wolbachia strains that infect arthropods with haplodiploid sex determination systems, including parasitoid wasps. Despite the discovery of microbe-mediated parthenogenesis more than 30 yr ago, the underlying genetic mechanisms have remained elusive. We used a suite of genomic, computational, and molecular tools to identify and characterize two proteins that are uniquely found in parthenogenesis-inducing Wolbachia and have strong signatures of host-associated bacterial effector proteins. These putative parthenogenesis-inducing proteins have structural homology to eukaryotic protein domains including nucleoporins, the key insect sex determining factor Transformer, and a eukaryotic-like serine-threonine kinase with leucine-rich repeats. Furthermore, these proteins significantly impact eukaryotic cell biology in the model Saccharomyces cerevisiae. We suggest that these proteins are parthenogenesis-inducing factors and our results indicate that this would be made possible by a novel mechanism of bacterial-host interaction.

Marinobacter qingdaonensis sp. nov., a moderately halotolerant bacterium isolated from intertidal sediment

A Gram-stain-negative, aerobic, rod-shaped and halotolerant bacterium, designated as strain ASW11-75T, was isolated from intertidal sediments in Qingdao, PR China, and identified using a polyphasic taxonomic approach. Growth of strain ASW11-75T occurred at 10-45 °C (optimum, 37 °C), pH 6.5-9.0 (optimum, pH 8.0) and 0.5-18.0 % NaCl concentrations (optimum, 2.5 %). Phylogenetic analyses based on 16S rRNA gene sequences and 1179 single-copy orthologous clusters indicated that strain ASW11-75T is affiliated with the genus Marinobacter. Strain ASW11-75T showed highest 16S rRNA gene sequence similarity to 'Marinobacter arenosus' CAU 1620T (98.5 %). The digital DNA-DNA hybridization and average nucleotide identity values between strain ASW11-75T and its closely related strains (Marinobacter salarius R9SW1T, Marinobacter similis A3d10T, 'Marinobacter arenosus' CAU 1620T, Marinobacter sediminum R65T, Marinobacter salinus Hb8T, Marinobacter alexandrii LZ-8T and Marinobacter nauticus ATCC 49840T) were 19.8-24.5 % and 76.6-80.7 %, respectively. The predominant cellular fatty acids were C16 : 0, C18 : 1  ω9c and C16 : 0 N alcohol. The polar lipids were phosphatidylethanolamine, phosphatidylglycerol, diphosphatidylglycerol, one unidentified aminophospholipid and two unidentified lipids. The major isoprenoid quinone was ubiquinone-9. The genomic DNA G+C content was 62.2 mol%. Based on genomic and gene function analysis, strain ASW11-75T had lower protein isoelectric points with higher ratios of acidic residues to basic residues and possessed genes related to ion transport and organic osmoprotectant uptake, implying its potential tolerance to salt. The results of polyphasic characterization indicated strain ASW11-75T represents a novel Marinobacter species, for which the name Marinobacter qingdaonensis sp. nov. with the type strain ASW11-75T is proposed. The type strain is ASW11-75T (=KCTC 82497T=MCCC 1K05587T).

Patterns recovered in phylogenomic analysis of Candida auris and close relatives implicate broad environmental flexibility in Candida/Clavispora clade yeasts

Fungal pathogens commonly originate from benign or non-pathogenic strains living in the natural environment. The recently emerged human pathogen, Candida auris, is one example of a fungus believed to have originated in the environment and recently transitioned into a clinical setting. To date, however, there is limited evidence about the origins of this species in the natural environment and when it began associating with humans. One approach to overcome this gap is to reconstruct phylogenetic relationships between (1) strains isolated from clinical and non-clinical environments and (2) between species known to cause disease in humans and benign environmental saprobes. C. auris belongs to the Candida/Clavispora clade, a diverse group of 45 yeast species including human pathogens and environmental saprobes. We present a phylogenomic analysis of the Candida/Clavispora clade aimed at understanding the ecological breadth and evolutionary relationships between an expanded sample of environmentally and clinically isolated yeasts. To build a robust framework for investigating these relationships, we developed a whole-genome sequence dataset of 108 isolates representing 18 species, including four newly sequenced species and 18 environmentally isolated strains. Our phylogeny, based on 619 orthologous genes, shows environmentally isolated species and strains interspersed with clinically isolated counterparts, suggesting that there have been many transitions between humans and the natural environment in this clade. Our findings highlight the breadth of environments these yeasts inhabit and imply that many clinically isolated yeasts in this clade could just as easily live outside the human body in diverse natural environments and vice versa.

Cryptic diversity of cellulose-degrading gut bacteria in industrialized humans

Humans, like all mammals, depend on the gut microbiome for digestion of cellulose, the main component of plant fiber. However, evidence for cellulose fermentation in the human gut is scarce. We have identified ruminococcal species in the gut microbiota of human populations that assemble functional multienzymatic cellulosome structures capable of degrading plant cell wall polysaccharides. One of these species, which is strongly associated with humans, likely originated in the ruminant gut and was subsequently transferred to the human gut, potentially during domestication where it underwent diversification and diet-related adaptation through the acquisition of genes from other gut microbes. Collectively, these species are abundant and widespread among ancient humans, hunter-gatherers, and rural populations but are rare in populations from industrialized societies thus indicating potential disappearance in response to the westernized lifestyle.

Multi-Omics integration can be used to rescue metabolic information for some of the dark region of the Pseudomonas putida proteome

In every omics experiment, genes or their products are identified for which even state of the art tools are unable to assign a function. In the biotechnology chassis organism Pseudomonas putida, these proteins of unknown function make up 14% of the proteome. This missing information can bias analyses since these proteins can carry out functions which impact the engineering of organisms. As a consequence of predicting protein function across all organisms, function prediction tools generally fail to use all of the types of data available for any specific organism, including protein and transcript expression information. Additionally, the release of Alphafold predictions for all Uniprot proteins provides a novel opportunity for leveraging structural information. We constructed a bespoke machine learning model to predict the function of recalcitrant proteins of unknown function in Pseudomonas putida based on these sources of data, which annotated 1079 terms to 213 proteins. Among the predicted functions supplied by the model, we found evidence for a significant overrepresentation of nitrogen metabolism and macromolecule processing proteins. These findings were corroborated by manual analyses of selected proteins which identified, among others, a functionally unannotated operon that likely encodes a branch of the shikimate pathway.

Redefining the treponemal history through pre-Columbian genomes from Brazil

The origins of treponemal diseases have long remained unknown, especially considering the sudden onset of the first syphilis epidemic in the late 15th century in Europe and its hypothesized arrival from the Americas with Columbus' expeditions1,2. Recently, ancient DNA evidence has revealed various treponemal infections circulating in early modern Europe and colonial-era Mexico3-6. However, there has been to our knowledge no genomic evidence of treponematosis recovered from either the Americas or the Old World that can be reliably dated to the time before the first trans-Atlantic contacts. Here, we present treponemal genomes from nearly 2,000-year-old human remains from Brazil. We reconstruct four ancient genomes of a prehistoric treponemal pathogen, most closely related to the bejel-causing agent Treponema pallidum endemicum. Contradicting the modern day geographical niche of bejel in the arid regions of the world, the results call into question the previous palaeopathological characterization of treponeme subspecies and showcase their adaptive potential. A high-coverage genome is used to improve molecular clock date estimations, placing the divergence of modern T. pallidum subspecies firmly in pre-Columbian times. Overall, our study demonstrates the opportunities within archaeogenetics to uncover key events in pathogen evolution and emergence, paving the way to new hypotheses on the origin and spread of treponematoses.

Genome biology and evolution of mating-type loci in four cereal rust fungi

Permanent heterozygous loci, such as sex- or mating-compatibility regions, often display suppression of recombination and signals of genomic degeneration. In Basidiomycota, two distinct loci confer mating compatibility. These loci encode homeodomain (HD) transcription factors and pheromone receptor (Pra)-ligand allele pairs. To date, an analysis of genome level mating-type (MAT) loci is lacking for obligate biotrophic basidiomycetes in the Pucciniales, an order containing serious agricultural plant pathogens. Here, we focus on four species of Puccinia that infect oat and wheat, including P. coronata f. sp. avenae, P. graminis f. sp. tritici, P. triticina and P. striiformis f. sp. tritici. MAT loci are located on two separate chromosomes supporting previous hypotheses of a tetrapolar mating compatibility system in the Pucciniales. The HD genes are multiallelic in all four species while the PR locus appears biallelic, except for P. graminis f. sp. tritici, which potentially has multiple alleles. HD loci are largely conserved in their macrosynteny, both within and between species, without strong signals of recombination suppression. Regions proximal to the PR locus, however, displayed signs of recombination suppression and genomic degeneration in the three species with a biallelic PR locus. Our observations support a link between recombination suppression, genomic degeneration, and allele diversity of MAT loci that is consistent with recent mathematical modelling and simulations. Finally, we confirm that MAT genes are expressed during the asexual infection cycle, and we propose that this may support regulating nuclear maintenance and pairing during infection and spore formation. Our study provides insights into the evolution of MAT loci of key pathogenic Puccinia species. Understanding mating compatibility can help predict possible combinations of nuclear pairs, generated by sexual reproduction or somatic recombination, and the potential evolution of new virulent isolates of these important plant pathogens.

Informatic challenges and advances in illuminating the druggable proteome

The understudied members of the druggable proteomes offer promising prospects for drug discovery efforts. While large-scale initiatives have generated valuable functional information on understudied members of the druggable gene families, translating this information into actionable knowledge for drug discovery requires specialized informatics tools and resources. Here, we review the unique informatics challenges and advances in annotating understudied members of the druggable proteome. We demonstrate the application of statistical evolutionary inference tools, knowledge graph mining approaches, and protein language models in illuminating understudied protein kinases, pseudokinases, and ion channels.

Draft genome sequence of three hydrocarbon-degrading Pseudomonadota strains isolated from an abandoned century-old oil exploration well

We present genome sequences of three Pseudomonadota strains isolated from an abandoned century-old oil exploration well. A Pseudomonas sp. genome showed a size of 5,378,420 bp, while Acinetobacter genomes sized 3,522,593 and 3,864,311 bp. Genomes included catabolic genes for benzoate, 4-hydroxybenzoate, salicylate, vanillate, indoleacetate, anthranilate, n-alkanes, 4-hydroxyphenylacetate, phenylacetate, among others.

Temporal profiling of physiological, histological, and transcriptomic dissection during auxin-induced adventitious root formation in tetraploid Robinia pseudoacacia micro-cuttings

MAIN CONCLUSION

Optimal levels of indole-3-butyric acid (IBA) applied at the stem base promote adventitious root (AR) initiation and primordia formation, thus promoting the rooting of leafy micro-cuttings of tetraploid Robinia pseudoacacia. Tetraploid Robinia pseudoacacia L. is a widely cultivated tree in most regions of China that has a hard-rooting capability, propagated by stem cuttings. This study utilizes histological, physiological, and transcriptomic approaches to explore how root primordia are induced after indole butyric acid (IBA) treatment of micro-cuttings. IBA application promoted cell divisions in some cells within the vasculature, showing subcellular features associated with adventitious root (AR) founder cells. The anatomical structure explicitly showed that AR initiated from the cambium layer and instigate the inducible development of AR primordia. Meanwhile, the hormone data showed that similar to that of indole-3-acetic acid, the contents of trans-zeatin and abscisic acid peaked at early stages of AR formation and increased gradually in primordia formation across the subsequent stages, suggesting their indispensable roles in AR induction. On the contrary, 24-epibrassinolide roughly maintained at extremely high levels during primordium initiation thoroughly, indicating its presence was involved in cell-specific reorganization during AR development. Furthermore, antioxidant activities transiently increased in the basal region of micro-cuttings and may serve as biochemical indicators for distinct rooting phases, potentially aiding in AR formation. Transcriptomic analysis during the early stages of root formation shows significant downregulation of the abscisic acid and jasmonate signaling pathways, while ethylene and cytokinin signaling seems upregulated. Network analysis of genes involved in carbon metabolism and photosynthesis indicates that the basal region of the micro-cuttings undergoes rapid reprogramming, which results in the breakdown of sugars into pyruvate. This pyruvate is then utilized to fuel the tricarboxylic acid cycle, thereby sustaining growth through aerobic respiration. Collectively, our findings provide a time-course morphophysiological dissection and also suggest the regulatory role of a conserved auxin module in AR development in these species.

The genome of Salmacisia buchloëana, the parasitic puppet master pulling strings of sexual phenotypic monstrosities in buffalograss

To complete its parasitic lifecycle, Salmacisia buchloëana, a biotrophic fungus, manipulates reproductive organ development, meristem determinacy, and resource allocation in its dioecious plant host, buffalograss (Bouteloua dactyloides; Poaceae). To gain insight into S. buchloëana's ability to manipulate its host, we sequenced and assembled the 20.1 Mb genome of S. buchloëana into 22 chromosome-level pseudomolecules. Phylogenetic analysis suggests that S. buchloëana is nested within the genus Tilletia and diverged from Tilletia caries and Tilletia walkeri ∼40 MYA. We find that S. buchloëana contains a novel chromosome arm with no syntenic relationship to other publicly available Tilletia genomes, and that genes on the novel arm are upregulated upon infection, suggesting that this unique chromosomal segment may have played a critical role in S. buchloëana's evolution and host specificity. Salmacisia buchloëana has one of the largest fractions of serine peptidases (1.53% of the proteome) and one of the highest GC contents (62.3%) in all classified fungi. Analysis of codon base composition indicated that GC content is controlled more by selective constraints than directional mutation, and that S. buchloëana has a unique bias for the serine codon UCG. Finally, we identify 3 inteins within the S. buchloëana genome, 2 of which are located in a gene often used in fungal taxonomy. The genomic and transcriptomic resources generated here will aid plant pathologists and breeders by providing insight into the extracellular components contributing to sex determination in dioecious grasses.

Alteromonas arenosi sp. nov., a novel bioflocculant-producing bacterium, isolated from intertidal sand

A novel Gram-stain-negative, strictly aerobic and bioflocculant-producing bacterium, designated as ASW11-36T, was isolated from an intertidal sand collected from coastal areas of Qingdao, PR China. Growth occurred at 15-40 °C (optimum, 30 °C), pH 7.0-9.0 (optimum, pH 7.5) and with 1.5-7.0% (w/v) NaCl (optimum, 2.5-3.0%). In the whole-cell fatty acid pattern prevailed C16:0 and summed feature 3 (C16:1 ω7c and/or C16:1 ω6c). The major isoprenoid quinone was determined to be Q-8 and the major polar lipids were phosphatidylethanolamine (PE) and phosphatidylglycerol (PG), one unidentified aminolipid (AL), one unidentified glycolipid (GL), and two lipids (L1, L2). Based on the phylogenetic analyses of 16S rRNA gene sequences and 618 single-copy orthologous clusters, strain ASW11-36T could represent a novel member of the genus Alteromonas and was closely related to Alteromonas flava P0211T (98.4%) and Alteromonas facilis P0213T (98.3%). The pairwise average nucleotide identity and digital DNA-DNA hybridization values of the ASW11-36T genome assembly against the closely related species genomes were 71.8% and 21.7%, respectively, that clearly lower than the proposed thresholds for species. Based on phenotypic, phylogenetic, and chemotaxonomic analyses, strain ASW11-36T is considered to represent a novel species of the genus Alteromonas, for which the name Alteromonas arenosi sp. nov. is proposed. The type strain is ASW11-36T (= KCTC 82496T = MCCC 1K05585T). In addition, the strain yielded 65% of flocculating efficiency in kaolin suspension with CaCl2 addition. The draft genome of ASW11-36T shared abundant putative CAZy family related genes, especially involved in the biosynthesis of exopolysaccharides, implying its potential environmental and biological applications.

Insights into the differences related to the resistance mechanisms to the highly toxic fruit Hippomane mancinella (Malpighiales: Euphorbiaceae) between the larvae of the sister species Anastrepha acris and Anastrepha ludens (Diptera: Tephritidae) through comparative transcriptomics

The Manchineel, Hippomane mancinella ("Death Apple Tree") is one of the most toxic fruits worldwide and nevertheless is the host plant of the monophagous fruit fly species Anastrepha acris (Diptera: Tephritidae). Here we aimed at elucidating the detoxification mechanisms in larvae of A. acris reared on a diet enriched with the toxic fruit (6% lyophilizate) through comparative transcriptomics. We compared the performance of A. acris larvae with that of the sister species A. ludens, a highly polyphagous pest species that is unable to infest H. mancinella in nature. The transcriptional alterations in A. ludens were significantly greater than in A. acris. We mainly found two resistance mechanisms in both species: structural, activating cuticle protein biosynthesis (chitin-binding proteins likely reducing permeability to toxic compounds in the intestine), and metabolic, triggering biosynthesis of serine proteases and xenobiotic metabolism activation by glutathione-S-transferases and cytochrome P450 oxidoreductase. Some cuticle proteins and serine proteases were not orthologous between both species, suggesting that in A. acris, a structural resistance mechanism has been selected allowing specialization to the highly toxic host plant. Our results represent a nice example of how two phylogenetically close species diverged over recent evolutionary time related to resistance mechanisms to plant secondary metabolites.

Refining and illuminating acetogenic Eubacterium strains for reclassification and metabolic engineering

BACKGROUND

The genus Eubacterium is quite diverse and includes several acetogenic strains capable of fermenting C1-substrates into valuable products. Especially, Eubacterium limosum and closely related strains attract attention not only for their capability to ferment C1 gases and liquids, but also due to their ability to produce butyrate. Apart from its well-elucidated metabolism, E. limosum is also genetically accessible, which makes it an interesting candidate to be an industrial biocatalyst.

RESULTS

In this study, we examined genomic, phylogenetic, and physiologic features of E. limosum and the closest related species E. callanderi as well as E. maltosivorans. We sequenced the genomes of the six Eubacterium strains 'FD' (DSM 3662T), 'Marburg' (DSM 3468), '2A' (DSM 2593), '11A' (DSM 2594), 'G14' (DSM 107592), and '32' (DSM 20517) and subsequently compared these with previously available genomes of the E. limosum type strain (DSM 20543T) as well as the strains 'B2', 'KIST612', 'YI' (DSM 105863T), and 'SA11'. This comparison revealed a close relationship between all eleven Eubacterium strains, forming three distinct clades: E. limosum, E. callanderi, and E. maltosivorans. Moreover, we identified the gene clusters responsible for methanol utilization as well as genes mediating chain elongation in all analyzed strains. Subsequent growth experiments revealed that strains of all three clades can convert methanol and produce acetate, butyrate, and hexanoate via reverse β-oxidation. Additionally, we used a harmonized electroporation protocol and successfully transformed eight of these Eubacterium strains to enable recombinant plasmid-based expression of the gene encoding the fluorescence-activating and absorption shifting tag (FAST). Engineered Eubacterium strains were verified regarding their FAST-mediated fluorescence at a single-cell level using a flow cytometry approach. Eventually, strains 'FD' (DSM 3662T), '2A' (DSM 2593), '11A' (DSM 2594), and '32' (DSM 20517) were genetically engineered for the first time.

CONCLUSION

Strains of E. limosum, E. callanderi, and E. maltosivorans are outstanding candidates as biocatalysts for anaerobic C1-substrate conversion into valuable biocommodities. A large variety of strains is genetically accessible using a harmonized electroporation protocol, and FAST can serve as a reliable fluorescent reporter protein to characterize genetically engineered cells. In total eleven strains have been assigned to distinct clades, providing a clear and updated classification. Thus, the description of respective Eubacterium species has been emended, improved, aligned, and is requested to be implemented in respective databases.

Diversity and conservation of the genome architecture of phages infecting the Alphaproteobacteria

IMPORTANCE

This study reports the results of the largest analysis of genome sequences from phages that infect the Alphaproteobacteria class of bacterial hosts. We analyzed over 100 whole genome sequences of phages to construct dotplots, categorize them into genetically distinct clusters, generate a bootstrapped phylogenetic tree, compute protein orthologs, and predict packaging strategies. We determined that the phage sequences primarily cluster by the bacterial host family, phage morphotype, and genome size. We expect that the findings reported in this seminal study will facilitate future analyses that will improve our knowledge of the phages that infect these hosts.

Comparative genomics sheds light on transcription factor-mediated regulation in the extreme acidophilic Acidithiobacillia representatives

Extreme acidophiles thrive in acidic environments, confront a multitude of challenges, and demonstrate remarkable adaptability in their metabolism to cope with the ever-changing environmental fluctuations, which encompass variations in temperature, pH levels, and the availability of electron acceptors and donors. The survival and proliferation of members within the Acidithiobacillia class rely on the deployment of transcriptional regulatory systems linked to essential physiological traits. The study of these transcriptional regulatory systems provides valuable insights into critical processes, such as energy metabolism and nutrient assimilation, and how they integrate into major genetic-metabolic circuits. In this study, we examined the transcriptional regulatory repertoires and potential interactions of forty-three Acidithiobacillia complete and draft genomes, encompassing nine species. To investigate the function and diversity of Transcription Factors (TFs) and their DNA Binding Sites (DBSs), we conducted a genome-wide comparative analysis, which allowed us to identify these regulatory elements in representatives of Acidithiobacillia. We classified TFs into gene families and compared their occurrence among all representatives, revealing conservation patterns across the class. The results identified conserved regulators for several pathways, including iron and sulfur oxidation, the main pathways for energy acquisition, providing new evidence for viable regulatory interactions and branch-specific conservation in Acidithiobacillia. The identification of TFs and DBSs not only corroborates existing experimental information for selected species, but also introduces novel candidates for experimental validation. Moreover, these promising candidates have the potential for further extension to new representatives within the class.

The Theory of Gene Family Histories

Most genes are part of larger families of evolutionary-related genes. The history of gene families typically involves duplications and losses of genes as well as horizontal transfers into other organisms. The reconstruction of detailed gene family histories, i.e., the precise dating of evolutionary events relative to phylogenetic tree of the underlying species has remained a challenging topic despite their importance as a basis for detailed investigations into adaptation and functional evolution of individual members of the gene family. The identification of orthologs, moreover, is a particularly important subproblem of the more general setting considered here. In the last few years, an extensive body of mathematical results has appeared that tightly links orthology, a formal notion of best matches among genes, and horizontal gene transfer. The purpose of this chapter is to broadly outline some of the key mathematical insights and to discuss their implication for practical applications. In particular, we focus on tree-free methods, i.e., methods to infer orthology or horizontal gene transfer as well as gene trees, species trees, and reconciliations between them without using a priori knowledge of the underlying trees or statistical models for the inference of phylogenetic trees. Instead, the initial step aims to extract binary relations among genes.

Family-Free Genome Comparison

The comparison of large-scale genome structures across distinct species offers valuable insights into the species' phylogeny, genome organization, and gene associations. In this chapter, we review the family-free genome comparison tool FFGC that, relying on built-in interfaces with a sequence comparison tool (either BLAST+ or DIAMOND) and with an ILP solver (either CPLEX or Gurobi), provides several methods for analyses that do not require prior classification of genes across the studied genomes. Taking annotated genome sequences as input, FFGC is a complete workflow for genome comparison allowing not only the computation of measures of similarity and dissimilarity but also the inference of gene families, simultaneously based on sequence similarities and large-scale genomic features.

Plasmid-encoded toxin defence mediates mutualistic microbial interactions

Gut environments harbour dense microbial ecosystems in which plasmids are widely distributed. Plasmids facilitate the exchange of genetic material among microorganisms while enabling the transfer of a diverse array of accessory functions. However, their precise impact on microbial community composition and function remains largely unexplored. Here we identify a prevalent bacterial toxin and a plasmid-encoded resistance mechanism that mediates the interaction between Lactobacilli and Enterococci. This plasmid is widespread across ecosystems, including the rumen and human gut microbiota. Biochemical characterization of the plasmid revealed a defence mechanism against reuterin, a toxin produced by various gut microbes, such as Limosilactobacillus reuteri. Using a targeted metabolomic approach, we find reuterin to be prevalent across rumen ecosystems with impacts on microbial community structure. Enterococcus strains carrying the protective plasmid were isolated and their interactions with L. reuteri, the toxin producer, were studied in vitro. Interestingly, we found that by conferring resistance against reuterin, the plasmid mediates metabolic exchange between the defending and the attacking microbial species, resulting in a beneficial relationship or mutualism. Hence, we reveal here an ecological role for a plasmid-coded defence system in mediating a beneficial interaction.

Listeria monocytogenes, a silent foodborne pathogen in Ecuador

Listeria monocytogenes is a foodborne pathogen that can produce serious, even fatal, infections. Among other foods, it can be found in unpasteurized dairy and ready-to-eat products. Surveillance of L. monocytogenes is of great interest since sources of infection are difficult to determine due to the long incubation period, and because the symptoms of listeriosis are similar to other diseases. We performed a genomic study of L. monocytogenes isolated from fresh cheeses and clinical samples from Ecuador. Sixty-five isolates were evaluated and sequenced, 14 isolates from cheese samples and 20 from clinical listeriosis cases from the National Institute of National Institute of Public Health Research, and 31 isolates from artisanal cheese samples from 8 provinces. All isolates exhibited heterogeneous patterns of the presence of pathogenicity islands. All isolates exhibited at least 4 genes from LIPI-1, but all references (26 L. monocytogenes closed genomes available in the NCBI database) showed the complete island, which encompasses 5 genes but is present in only two Ecuadorian isolates. Most isolates lacked gene actA. Genes from LIPI-2 were absent in all isolates. LIPI-3 and LIPI-4 were present in only a few references and isolates. With respect to the stress survival islets, our samples either presented SSI-1 or SSI-F2365, except for one isolate that presented SSI-F2365 and also one gene from SSI-1. None of the samples presented SSI-2. The predominant ST (sequence type) was ST2 (84.62% 55/65), and the only ST found in food (93.33% 42/45) and clinical samples (65% 13/20). Isolates were not grouped according to their sampling origin, date, or place in a phylogenetic tree obtained from the core alignment. The presence of ST2 in food and clinical samples, with high genomic similarity, suggests a foodborne infection risk linked to the consumption of fresh cheeses in Ecuador.

Proteinortho6: pseudo-reciprocal best alignment heuristic for graph-based detection of (co-)orthologs

Proteinortho is a widely used tool to predict (co)-orthologous groups of genes for any set of species. It finds application in comparative and functional genomics, phylogenomics, and evolutionary reconstructions. With a rapidly increasing number of available genomes, the demand for large-scale predictions is also growing. In this contribution, we evaluate and implement major algorithmic improvements that significantly enhance the speed of the analysis without reducing precision. Graph-based detection of (co-)orthologs is typically based on a reciprocal best alignment heuristic that requires an all vs. all comparison of proteins from all species under study. The initial identification of similar proteins is accelerated by introducing an alternative search tool along with a revised search strategy-the pseudo-reciprocal best alignment heuristic-that reduces the number of required sequence comparisons by one-half. The clustering algorithm was reworked to efficiently decompose very large clusters and accelerate processing. Proteinortho6 reduces the overall processing time by an order of magnitude compared to its predecessor while maintaining its small memory footprint and good predictive quality.

A timeline of bacterial and archaeal diversification in the ocean

Microbial plankton play a central role in marine biogeochemical cycles, but the timing in which abundant lineages diversified into ocean environments remains unclear. Here, we reconstructed the timeline in which major clades of bacteria and archaea colonized the ocean using a high-resolution benchmarked phylogenetic tree that allows for simultaneous and direct comparison of the ages of multiple divergent lineages. Our findings show that the diversification of the most prevalent marine clades spans throughout a period of 2.2 Ga, with most clades colonizing the ocean during the last 800 million years. The oldest clades - SAR202, SAR324, Ca. Marinimicrobia, and Marine Group II - diversified around the time of the Great Oxidation Event, during which oxygen concentration increased but remained at microaerophilic levels throughout the Mid-Proterozoic, consistent with the prevalence of some clades within these groups in oxygen minimum zones today. We found the diversification of the prevalent heterotrophic marine clades SAR11, SAR116, SAR92, SAR86, and Roseobacter as well as the Marine Group I to occur near to the Neoproterozoic Oxygenation Event (0.8-0.4 Ga). The diversification of these clades is concomitant with an overall increase of oxygen and nutrients in the ocean at this time, as well as the diversification of eukaryotic algae, consistent with the previous hypothesis that the diversification of heterotrophic bacteria is linked to the emergence of large eukaryotic phytoplankton. The youngest clades correspond to the widespread phototrophic clades Prochlorococcus, Synechococcus, and Crocosphaera, whose diversification happened after the Phanerozoic Oxidation Event (0.45-0.4 Ga), in which oxygen concentrations had already reached their modern levels in the atmosphere and the ocean. Our work clarifies the timing at which abundant lineages of bacteria and archaea colonized the ocean, thereby providing key insights into the evolutionary history of lineages that comprise the majority of prokaryotic biomass in the modern ocean.

Identification of parthenogenesis-inducing effector proteins in Wolbachia

Bacteria in the genus Wolbachia have evolved numerous strategies to manipulate arthropod sex, including the conversion of would-be male offspring to asexually reproducing females. This so-called "parthenogenesis-induction" phenotype can be found in a number of Wolbachia strains that infect arthropods with haplodiploid sex determination systems, including parasitoid wasps. Despite the discovery of microbe-mediated parthenogenesis more than 30 years ago, the underlying genetic mechanisms have remained elusive. We used a suite of genomic, computational, and molecular tools to identify and characterize two proteins that are uniquely found in parthenogenesis-inducing Wolbachia and have strong signatures of host-associated bacterial effector proteins. These putative parthenogenesis-inducing proteins have structural homology to eukaryotic protein domains including nucleoporins, the key insect sex-determining factor Transformer, and a eukaryotic-like serine-threonine kinase with leucine rich repeats. Furthermore, these proteins significantly impact eukaryotic cell biology in the model, Saccharomyces cerevisiae. We suggest these proteins are parthenogenesis-inducing factors and our results indicate this would be made possible by a novel mechanism of bacterial-host interaction.

Genome-scale approach to reconstructing the phylogenetic tree of psyllids (superfamily Psylloidea) with account of systematic bias

Psyllids (class Insecta: order Hemiptera: superfamily Psylloidea) are a taxonomically and phylogenetically challenging clade. Recent studies have largely advanced the phylogeny of this group, yet the family-level relationships among Aphalaridae, Carsidaridae, and others remain unresolved. Genome-scale phylogenetic analysis is known to provide a finer resolution for problems like that. However, such phylogenomics also introduces new problems: incorrect trees with high confidence yielded due to systematic error (bias). Here we addressed these issues using hundreds of single-copy orthologous (SCO) genes in psyllid transcriptomes and genomes. Our analyses revealed conflicts between the nucleotide-based and amino-acid-based phylogenetic trees. While the nucleotide-based phylogeny strongly supported the (Aphalaridae + Carsidaridae) + Others relationship, the amino-acid-based one recovered Aphalaridae + (Carsidaridae + Others) with 100% support. Further inspection revealed significant compositional heterogeneity in nucleotide sequences for 67% of SCO genes, but not in the corresponding translated amino acid sequences. We then used different strategies to combat this compositional bias, and found that using the RY-coding strategy (coding the standard nucleotides as purines and pyrimidines) the nucleotide-based phylogeny became consistent with the amino-acid-based one. We further applied RY-coding to a published concatenated nucleotide dataset and recovered the Aphalaridae monophyly (which is refuted by the original literature on non-recoded sequences) at the base of psyllid tree. Moreover, it was found that variations in evolutionary rate could lead to errors in nucleotide-based phylogeny. The fast-evolving Heteropsylla cubana (Psyllidae: Ciriacreminae) was incorrectly placed within the subfamily Psyllinae. This bias can be avoided by using data removal or RY-coding strategies. Together, our results strongly support the family relationship of Aphalaridae + (Carsidaridae + Others), and show that the amino-acid-based concatenation analysis is more robust than nucleotide-based one. Future phylogenomic analysis of psyllid nucleotide sequences should take into account methods such as the RY-coding scheme to address potential systematic biases arising from composition and rate heterogeneities.

Species- and strain-level diversity of Corynebacteria isolated from human facial skin

BACKGROUND

Sequencing of the human skin microbiome revealed that Corynebacterium is an ubiquitous and abundant bacterial genus on human skin. Shotgun sequencing further highlighted the microbial "dark matter" of the skin microbiome, consisting of microorganisms, including corynebacterial species that were not cultivated and genome-sequenced so far. In this pilot project, facial human skin swabs of 13 persons were cultivated to selectively obtain corynebacteria. 54 isolates were collected and 15 of these were genome-sequenced and the pan-genome was determined. The strains were biochemically characterized and antibiotic susceptibility testing (AST) was performed.

RESULTS

Among the 15 sequenced strains, nine different corynebacterial species were found, including two so far undescribed species, tentatively named "Corynebacterium vikingii" and "Corynebacterium borealis", for which closed genome sequences were obtained. Strain variability beyond the species level was determined in biochemical tests, such as the variable presence of urease activity and the capacity to ferment different sugars. The ability to grow under anaerobic conditions on solid agar was found to be species-specific. AST revealed resistances to clindamycin in seven strains. A Corynebacterium pseudokroppenstedtii strain showed additional resistance towards beta-lactam and fluoroquinolone antibiotics; a chromosomally located 17 kb gene cluster with five antibiotic resistance genes was found in the closed genome of this strain.

CONCLUSIONS

Taken together, this pilot study identified an astonishing diversity of cutaneous corynebacterial species in a relatively small cohort and determined species- and strain-specific individualities regarding biochemical and resistance profiles. This further emphasizes the need for cultivation-based studies to be able to study these microorganisms in more detail, in particular regarding their host-interacting and, potentially, -beneficial and/or -detrimental properties.

Host macrocyclic acylcarnitines mediate symbiotic interactions between frogs and their skin microbiome

The host-microbiome associations occurring on the skin of vertebrates significantly influence hosts' health. However, the factors mediating their interactions remain largely unknown. Herein, we used integrated technical and ecological frameworks to investigate the skin metabolites sustaining a beneficial symbiosis between tree frogs and bacteria. We characterize macrocyclic acylcarnitines as the major metabolites secreted by the frogs' skin and trace their origin to an enzymatic unbalance of carnitine palmitoyltransferases. We found that these compounds colocalize with bacteria on the skin surface and are mostly represented by members of the Pseudomonas community. We showed that Pseudomonas sp. MPFS isolated from frogs' skin can exploit acylcarnitines as its sole carbon and nitrogen source, and this metabolic capability is widespread in Pseudomonas. We summarize frogs' multiple mechanisms to filter environmental bacteria and highlight that acylcarnitines likely evolved for another function but were co-opted to provide nutritional benefits to the symbionts.

De novo genome assembly and annotation of the medicinal plant Tinospora cordifolia (Willd.) Miers ex Hook. f. & Thom's

Indian natural climbing shrub Tinospora cordifolia, often known as "Guduchi" and "Amrita," is a highly esteemed medicinal plant in the Indian system of medicine (ISM). It is a member of the Menispermaceae family which consists of a rich source of protein, micronutrients, and rich source of bioactive components which are used in treating various systemic diseases. The current study was designed to know the biological characterization of the plant genome and biosynthesis of plant metabolites essential for its medicinal applications. Tinospora cordifolia's complete genome was sequenced using Illumina HiSeq2500 sequencing technology. The draft genome was assembled through a de novo method. An integrative genome annotation approach was used to perform functional gene prediction. The pathway analysis was carried out using the KEGG database. The total genome size obtained after genome assembly was 894 Mb with an N50 of 9148 bp. The integrative annotation approach resulted in 35,111 protein-coding genes. In addition, genes responsible for the synthesis of syringin, a secondary metabolite found in plants, were identified. In comparison to the standard drug (dopamine, rasagiline, and selegiline), syringin's molecular docking exhibited a greater binding affinity from the range of - 4.3 to - 6.6 kcal/mol for all the targets of Parkinson's disease and for Alzheimer's targets; it has shown the maximum potency from the range of - 6.5 to - 7.4 kcal/mol with respect to the standard drug (donepezil, galantamine, and rivastigmine). This study provides the genomic information of Tinospora cordifolia which is helpful in understanding genomic insights and metabolic pathways connected to the corresponding plant genome and predicts the possible useful effect for the molecular characterization of therapeutic drugs.

Comparative Population Transcriptomics Provide New Insight into the Evolutionary History and Adaptive Potential of World Ocean Krill

Genetic variation is instrumental for adaptation to changing environments but it is unclear how it is structured and contributes to adaptation in pelagic species lacking clear barriers to gene flow. Here, we applied comparative genomics to extensive transcriptome datasets from 20 krill species collected across the Atlantic, Indian, Pacific, and Southern Oceans. We compared genetic variation both within and between species to elucidate their evolutionary history and genomic bases of adaptation. We resolved phylogenetic interrelationships and uncovered genomic evidence to elevate the cryptic Euphausia similis var. armata into species. Levels of genetic variation and rates of adaptive protein evolution vary widely. Species endemic to the cold Southern Ocean, such as the Antarctic krill Euphausia superba, showed less genetic variation and lower evolutionary rates than other species. This could suggest a low adaptive potential to rapid climate change. We uncovered hundreds of candidate genes with signatures of adaptive evolution among Antarctic Euphausia but did not observe strong evidence of adaptive convergence with the predominantly Arctic Thysanoessa. We instead identified candidates for cold-adaptation that have also been detected in Antarctic fish, including genes that govern thermal reception such as TrpA1. Our results suggest parallel genetic responses to similar selection pressures across Antarctic taxa and provide new insights into the adaptive potential of important zooplankton already affected by climate change.

Key processes required for the different stages of fungal carnivory by a nematode-trapping fungus

Nutritional deprivation triggers a switch from a saprotrophic to predatory lifestyle in soil-dwelling nematode-trapping fungi (NTF). In particular, the NTF Arthrobotrys oligospora secretes food and sex cues to lure nematodes to its mycelium and is triggered to develop specialized trapping devices. Captured nematodes are then invaded and digested by the fungus, thus serving as a food source. In this study, we examined the transcriptomic response of A. oligospora across the stages of sensing, trap development, and digestion upon exposure to the model nematode Caenorhabditis elegans. A. oligospora enacts a dynamic transcriptomic response, especially of protein secretion-related genes, in the presence of prey. Two-thirds of the predicted secretome of A. oligospora was up-regulated in the presence of C. elegans at all time points examined, and among these secreted proteins, 38.5% are predicted to be effector proteins. Furthermore, functional studies disrupting the t-SNARE protein Sso2 resulted in impaired ability to capture nematodes. Additionally, genes of the DUF3129 family, which are expanded in the genomes of several NTF, were highly up-regulated upon nematode exposure. We observed the accumulation of highly expressed DUF3129 proteins in trap cells, leading us to name members of this gene family as Trap Enriched Proteins (TEPs). Gene deletion of the most highly expressed TEP gene, TEP1, impairs the function of traps and prevents the fungus from capturing prey efficiently. In late stages of predation, we observed up-regulation of a variety of proteases, including metalloproteases. Following penetration of nematodes, these metalloproteases facilitate hyphal growth required for colonization of prey. These findings provide insights into the biology of the predatory lifestyle switch in a carnivorous fungus and provide frameworks for other fungal-nematode predator-prey systems.

Symmetric expression of ohnologs encoding conserved antiviral responses in tetraploid common carp suggest absence of subgenome dominance after whole genome duplication

Allopolyploids often experience subgenome dominance, with one subgenome showing higher levels of gene expression and greater gene retention. Here, we address the functionality of both subgenomes of allotetraploid common carp (Cyprinus carpio) by analysing a functional network of interferon-stimulated genes (ISGs) crucial in anti-viral immune defence. As an indicator of subgenome dominance we investigated retainment of a core set of ohnologous ISGs. To facilitate our functional genomic analysis a high quality genome was assembled (WagV4.0). Transcriptome data from an in vitro experiment mimicking a viral infection was used to infer ISG expression. Transcriptome analysis confirmed induction of 88 ISG ohnologs on both subgenomes. In both control and infected states, average expression of ISG ohnologs was comparable between the two subgenomes. Also, the highest expressing and most inducible gene copies of an ohnolog pair could be derived from either subgenome. We found no strong evidence of subgenome dominance for common carp.

RNA localization during early development of the axolotl

The asymmetric localization of biomolecules is critical for body plan development. One of the most popular model organisms for early embryogenesis studies is Xenopus laevis but there is a lack of information in other animal species. Here, we compared the early development of two amphibian species-the frog X. laevis and the axolotl Ambystoma mexicanum. This study aimed to identify asymmetrically localized RNAs along the animal-vegetal axis during the early development of A. mexicanum. For that purpose, we performed spatial transcriptome-wide analysis at low resolution, which revealed dynamic changes along the animal-vegetal axis classified into the following categories: profile alteration, de novo synthesis and degradation. Surprisingly, our results showed that many of the vegetally localized genes, which are important for germ cell development, are degraded during early development. Furthermore, we assessed the motif presence in UTRs of degraded mRNAs and revealed the enrichment of several motifs in RNAs of germ cell markers. Our results suggest novel reorganization of the transcriptome during embryogenesis of A. mexicanum to converge to the similar developmental pattern as the X. laevis.

Assessing sequence heterogeneity in Chlorellaceae DNA barcode markers for phylogenetic inference

Phylogenetic inference is an important approach that allows the recovery of the evolutionary history and the origin of the Chlorellaceae species. Despite the species' potential for biofuel feedstock production, their high phenotypic plasticity and similar morphological structures among the species have muddled the taxonomy and identification of the Chlorellaceae species. This study aimed to decipher Chlorellaceae DNA barcode marker heterogeneity by examining the sequence divergence and genomic properties of 18S rRNA, ITS (ITS1-5.8S rRNA-ITS2-28S rRNA), and rbcL from 655 orthologous sequences of 64 species across 31 genera in the Chlorellaceae family. The study assessed the distinct evolutionary properties of the DNA markers that may have caused the discordance between individual trees in the phylogenetic inference using the Robinson-Foulds distance and the Shimodaira-Hasegawa test. Our findings suggest that using the supermatrix approach improves the congruency between trees by reducing stochastic error and increasing the confidence of the inferred Chlorellaceae phylogenetic tree. This study also found that the phylogenies inferred through the supermatrix approach might not always be well supported by all markers. The study highlights that assessing sequence heterogeneity prior to the phylogenetic inference could allow the approach to accommodate sequence evolutionary properties and support species identification from the most congruent phylogeny, which can better represent the evolution of Chlorellaceae species.

Parerythrobacter aestuarii sp. nov., Isolated from Seawater in the Tidal Region of Taizhou

A yellow, Gram-stain-negative, aerobic, and rod-shaped strain, designated as C18T, was isolated from seawater in the tidal region of Taizhou. Growth of strain C18T occurs at 20-45 °C, at pH 5.5-8.0 and with 1.0-8.0% (w/v) NaCl. The 16S rRNA gene sequence analysis showed that strain C18T shared sequence identities with the genera Erythrobacter (< 98.4%), Qipengyuania (< 98.0%), Altererythrobacter (< 96.4%), Parerythrobacter (< 96.2%), Aurantiacibacter (< 96.2%), Tsuneonella (< 96.0%), Pelagerythrobacter (< 96.0%), Alteriqipengyuania (< 95.9%), and Parapontixanthobacter (< 95.7%) type strains. While the phylogenomic tree based on single-copy orthologous clusters revealed that strain C18T was stably clustered into the genus Parerythrobacter. Average nucleotide identity and digital DNA-DNA hybridization values of strain C18T and Parerythrobacter type strains were 73.5-75.2% and 18.5-19.4%, respectively, which were lower than the species delineation thresholds. The sole respiratory quinones were identified as ubiquinone-10. The major fatty acids (> 10%) were C17:1ω6c and summed feature 8 (C18:1ω7c and/or C18:1ω6c). Polar lipids included diphosphatidylglycerol, phosphatidylethanolamine, phosphatidylglycerol, sphingoglycolipid, an unidentified phospholipid, and an unidentified aminophospholipid. Based on the genetic, chemotaxonomic and phenotypic results, strain C18T is concluded to represent a novel species in the genus Parerythrobacter, for which the name Parerythrobacter aestuarii sp. nov. is proposed. The type strain is C18T (= KCTC 82594T = MCCC 1K05109T).

Genome-based metabolic and phylogenomic analysis of three Terrisporobacter species

Acetogenic bacteria are of high interest for biotechnological applications as industrial platform organisms, however, acetogenic strains from the genus Terrisporobacter have hitherto been neglected. To date, three published type strains of the genus Terrisporobacter are only covered by draft genome sequences, and the genes and pathway responsible for acetogenesis have not been analyzed. Here, we report complete genome sequences of the bacterial type strains Terrisporobacter petrolearius JCM 19845T, Terrisporobacter mayombei DSM 6539T and Terrisporobacter glycolicus DSM 1288T. Functional annotation, KEGG pathway module reconstructions and screening for virulence factors were performed. Various species-specific vitamin, cofactor and amino acid auxotrophies were identified and a model for acetogenesis of Terrisporobacter was constructed. The complete genomes harbored a gene cluster for the reductive proline-dependent branch of the Stickland reaction located on an approximately 21 kb plasmid, which is exclusively found in the Terrisporobacter genus. Phylogenomic analysis of available Terrisporobacter genomes suggested a reclassification of most isolates as T. glycolicus into T. petrolearius.

AlexandrusPS: A User-Friendly Pipeline for the Automated Detection of Orthologous Gene Clusters and Subsequent Positive Selection Analysis

The detection of adaptive selection in a system approach considering all protein-coding genes allows for the identification of mechanisms and pathways that enabled adaptation to different environments. Currently, available programs for the estimation of positive selection signals can be divided into two groups. They are either easy to apply but can analyze only one gene family at a time, restricting system analysis; or they can handle larger cohorts of gene families, but require considerable prerequisite data such as orthology associations, codon alignments, phylogenetic trees, and proper configuration files. All these steps require extensive computational expertise, restricting this endeavor to specialists. Here, we introduce AlexandrusPS, a high-throughput pipeline that overcomes technical challenges when conducting transcriptome-wide positive selection analyses on large sets of nucleotide and protein sequences. The pipeline streamlines 1) the execution of an accurate orthology prediction as a precondition for positive selection analysis, 2) preparing and organizing configuration files for CodeML, 3) performing positive selection analysis using CodeML, and 4) generating an output that is easy to interpret, including all maximum likelihood and log-likelihood test results. The only input needed from the user is the CDS and peptide FASTA files of proteins of interest. The pipeline is provided in a Docker image, requiring no program or module installation, enabling the application of the pipeline in any computing environment. AlexandrusPS and its documentation are available via GitHub (https://github.com/alejocn5/AlexandrusPS).

Genome sequence analysis of Malayan pangolin (Manis javanica) forensic samples reveals the presence of Paraburkholderia fungorum sequences

Background

The Malayan pangolin (Manis javanica) is a placental mammal and is listed as Critically Endangered on the IUCN Red List of Threatened Species. Most previous attempts to breed pangolins in captivity have met with little success because of dietary issues, infections, and other complications, although a previous study reported breeding pangolins in captivity to the third generation. In our previous pangolin genome sequencing data analysis, we obtained a considerable amount of bacterial DNA from a pregnant female Malayan pangolin (named "UM3"), which was likely infected by Paraburkholderia fungorum-an agent of biodegradation and bioremediation in agriculture.

Methodology

Here, we further confirmed and characterized this bacterial species using PCR, histological staining, whole-genome sequencing, and bioinformatics approaches. PCR assays with in-house designed primer sets and 16S universal primers showed clear positive bands in the cerebrum, cerebellum, lung, and blood of UM3 suggesting that UM3 might have developed septicaemia. Histological staining showed the presence of Gram-negative rod-shaped bacteria in the pangolin brain and lungs, indicating the colonization of the bacteria in these two organs. In addition, PCR screening of UM3's fetal tissues revealed the presence of P. fungorum in the gastrocnemius muscle, but not in other tissues that we examined. We also sequenced and reconstructed the genome of pangolin P. fungorum, which has a genome size of 7.7 Mbps.

Conclusion

Our study is the first to present detailed evidence of the presence of P. fungorum in a pangolin and her fetus (although preliminary results were presented in our previous article). Here, we raise the concern that P. fungorum may potentially infect humans, especially YOPI (young, old, pregnant, and immunocompromised) people. Therefore, caution should be exercised when using this bacterial species as biodegradation or bioremediation agents in agriculture.

Alteromonas Aquimaris sp. nov., Isolated from Surface Seawater

A Gram-negative, aerobic, motile by flagellum, and rod-shaped bacterium, designated ASW11-7T, was isolated from coastal surface seawater sample collected from the Yellow Sea, PR China. Strain ASW11-7T grew optimally at 37℃, 4.0% (w/v) NaCl and pH 7.0. Phylogenetic analysis based on the 16S rRNA gene sequences revealed that strain ASW11-7T belongs to the genus Alteromonas and most closely related to Alteromonas ponticola MYP5T (99.6% similarity), followed by Alteromonas confluentis DSSK2-12T (98.2%), Alteromonas lipolytica JW12T (98.2%), and Alteromonas hispanica F-32T (98.0%). The draft genome of strain ASW11-7T had a length of 3,530,922 bp with a G + C content of 44.9%, predicting 3108 coding sequences, 5 rRNA, 4 ncRNAs, 49 tRNAs genes, and 18 pseudogenes. The average nucleotide identity and digital DNA-DNA hybridization values between genomic sequences of strain ASW11-7T and closely related species of Alteromonas were in ranges of 66.9-77.8% and 18.3-27.5%, respectively. The major fatty acids of strain ASW11-7T were C16:0, summed feature 3 (C16:1ω7c/C16:1ω6c), and summed feature 8 (C18:1ω7c/C18:1ω6c). The predominant respiratory quinone was Q-8 and the major polar lipids were phosphatidylethanolamine and phosphatidylglycerol. Based on the phenotypic properties, genotypic distinctiveness, and chemotaxonomic features, strain ASW11-7T is considered to represent a novel Alteromonas species, for which the name Alteromonas aquimaris sp. nov. is proposed. The type strain is ASW11-7T (= KCTC 92853T = MCCC 1K07240T).

Simvastatin induces human gut bacterial cell surface genes

Drugs intended to target mammalian cells can have broad off-target effects on the human gut microbiota with potential downstream consequences for drug efficacy and side effect profiles. Yet, despite a rich literature on antibiotic resistance, we still know very little about the mechanisms through which commensal bacteria evade non-antibiotic drugs. Here, we focus on statins, one of the most prescribed drug types in the world and an essential tool in the prevention and treatment of high circulating cholesterol levels. Prior work in humans, mice, and cell culture support an off-target effect of statins on human gut bacteria; however, the genetic determinants of statin sensitivity remain unknown. We confirmed that simvastatin inhibits the growth of diverse human gut bacterial strains grown in communities and in pure cultures. Drug sensitivity varied between phyla and was dose-dependent. We selected two representative simvastatin-sensitive species for more in-depth analysis: Eggerthella lenta (phylum: Actinobacteriota) and Bacteroides thetaiotaomicron (phylum: Bacteroidota). Transcriptomics revealed that both bacterial species upregulate genes in response to simvastatin that alter the cell membrane, including fatty acid biogenesis (E. lenta) and drug efflux systems (B. thetaiotaomicron). Transposon mutagenesis identified a key efflux system in B. thetaiotaomicron that enables growth in the presence of statins. Taken together, these results emphasize the importance of the bacterial cell membrane in countering the off-target effects of host-targeted drugs. Continued mechanistic dissection of the various mechanisms through which the human gut microbiota evades drugs will be essential to understand and predict the effects of drug administration in human cohorts and the potential downstream consequences for health and disease.

Evidence that the domesticated fungus Leucoagaricus gongylophorus recycles its cytoplasmic contents as nutritional rewards to feed its leafcutter ant farmers

Leafcutter ants farm a fungal cultivar (Leucoagaricus gongylophorus) that converts inedible vegetation into food that sustains colonies with up to millions of workers. Analogous to edible fruits of crops domesticated by humans, L. gongylophorus has evolved specialized nutritional rewards-swollen hyphal cells called gongylidia that package metabolites and are consumed by ant farmers. Yet, little is known about how gongylidia form, and thus how fungal physiology and ant provisioning collectively govern farming performance. We explored the process of gongylidium formation using advanced microscopy to image the cultivar at scales of nanometers, and both in vitro experiments and in silico analyses to examine the mechanisms of gongylidia formation when isolated from ant farmers. We first used transmission electron, fluorescence, and confocal microscopy imaging to see inside hyphal cells. This imaging showed that the cultivar uses a process called autophagy to recycle its own cellular material (e.g. cytosol, mitochondria) and then shuttles the resulting metabolites into a vacuole whose continual expansion displaces other organelles and causes the gongylidium cell's bulging bulb-like appearance. We next used scanning electron microscopy and light microscopy to link this intracellular rearrangement to the external branching patterns of gongylidium cells as they clump together into edible bundles called staphyla. We next confirmed that autophagy plays a critical role in gongylidium formation both: (1) in vitro as gongylidium suppression occurred when isolated fungal cultures were grown on media with autophagy inhibitors, and (2) in silico as differential transcript expression (RNA-seq) analyses showed upregulation of multiple autophagy gene isoforms in gongylidia relative to undifferentiated hyphae. While autophagy is a ubiquitous and often highly derived process across the tree of life, our study reveals a new role for autophagy as a mechanism of functional integration between ant farmers and their fungal crop, and potentially as a signifier of higher-level homeostasis between uniquely life-time committed ectosymbionts.

The origin of eukaryotes and rise in complexity were synchronous with the rise in oxygen

The origin of eukaryotes was among the most important events in the history of life, spawning a new evolutionary lineage that led to all complex multicellular organisms. However, the timing of this event, crucial for understanding its environmental context, has been difficult to establish. The fossil and biomarker records are sparse and molecular clocks have thus far not reached a consensus, with dates spanning 2.1-0.91 billion years ago (Ga) for critical nodes. Notably, molecular time estimates for the last common ancestor of eukaryotes are typically hundreds of millions of years younger than the Great Oxidation Event (GOE, 2.43-2.22 Ga), leading researchers to question the presumptive link between eukaryotes and oxygen. We obtained a new time estimate for the origin of eukaryotes using genetic data of both archaeal and bacterial origin, the latter rarely used in past studies. We also avoided potential calibration biases that may have affected earlier studies. We obtained a conservative interval of 2.2-1.5 Ga, with an even narrower core interval of 2.0-1.8 Ga, for the origin of eukaryotes, a period closely aligned with the rise in oxygen. We further reconstructed the history of biological complexity across the tree of life using three universal measures: cell types, genes, and genome size. We found that the rise in complexity was temporally consistent with and followed a pattern similar to the rise in oxygen. This suggests a causal relationship stemming from the increased energy needs of complex life fulfilled by oxygen.

Multiomic interpretation of fungus-infected ant metabolomes during manipulated summit disease

Camponotus floridanus ants show altered behaviors followed by a fatal summiting phenotype when infected with manipulating Ophiocordyceps camponoti-floridani fungi. Host summiting as a strategy to increase transmission is also observed with parasite taxa beyond fungi, including aquatic and terrestrial helminths and baculoviruses. The drastic phenotypic changes can sometimes reflect significant molecular changes in gene expression and metabolite concentrations measured in manipulated hosts. Nevertheless, the underlying mechanisms still need to be fully characterized. To investigate the small molecules producing summiting behavior, we infected C. floridanus ants with O. camponoti-floridani and sampled their heads for LC-MS/MS when we observed the characteristic summiting phenotype. We link this metabolomic data with our previous genomic and transcriptomic data to propose mechanisms that underlie manipulated summiting behavior in "zombie ants." This "multiomic" evidence points toward the dysregulation of neurotransmitter levels and neuronal signaling. We propose that these processes are altered during infection and manipulation based on (1) differential expression of neurotransmitter synthesis and receptor genes, (2) altered abundance of metabolites and neurotransmitters (or their precursors) with known behavioral effects in ants and other insects, and (3) possible suppression of a connected immunity pathway. We additionally report signals for metabolic activity during manipulation related to primary metabolism, detoxification, and anti-stress protectants. Taken together, these findings suggest that host manipulation is likely a multi-faceted phenomenon, with key processes changing at multiple levels of molecular organization.

Genome assemblies of the simultaneously hermaphroditic flatworms Macrostomum cliftonense and Macrostomum hystrix

The free-living, simultaneously hermaphroditic flatworms of the genus Macrostomum are increasingly used as model systems in various contexts. In particular, Macrostomum lignano, the only species of this group with a published genome assembly, has emerged as a model for the study of regeneration, reproduction, and stem-cell function. However, challenges have emerged due to M. lignano being a hidden polyploid, having recently undergone whole-genome duplication and chromosome fusion events. This complex genome architecture presents a significant roadblock to the application of many modern genetic tools. Hence, additional genomic resources for this genus are needed. Here, we present such resources for Macrostomum cliftonense and Macrostomum hystrix, which represent the contrasting mating behaviors of reciprocal copulation and hypodermic insemination found in the genus. We use a combination of PacBio long-read sequencing and Illumina shot-gun sequencing, along with several RNA-Seq data sets, to assemble and annotate highly contiguous genomes for both species. The assemblies span ∼227 and ∼220 Mb and are represented by 399 and 42 contigs for M. cliftonense and M. hystrix, respectively. Furthermore, high BUSCO completeness (∼84-85%), low BUSCO duplication rates (8.3-6.2%), and low k-mer multiplicity indicate that these assemblies do not suffer from the same assembly ambiguities of the M. lignano genome assembly, which can be attributed to the complex karyology of this species. We also show that these resources, in combination with the prior resources from M. lignano, offer an excellent foundation for comparative genomic research in this group of organisms.