Ggtree: A serialized data object for visualization of a phylogenetic tree and annotation data

Association analysis of antibiotic and disinfectant resistome in human and foodborne E. coli in Beijing, China

The widespread use of chemical disinfectants and antibiotics poses a major threat to food safety and human health. However, the mechanisms of co-transmission of antimicrobial resistance genes (ARGs) and biocides and metal resistance genes (BMRGs) of foodborne pathogens in the food chain is still unclear. This study isolated 343 E. coli strains from animal-derived foods in Beijing and incorporated online data of human-derived E. coli strains from NCBI. Our results demonstrated a relatively uniform distribution of strains from various regions in Beijing, indicating a lack of region-specific clustering. Additionally, predominant sequence types varied between food- and human-derived strains, suggesting a preference for different hosts and environments. Phenotypic association analysis showed that the chlorine disinfectants peroxides had a significant positive correlation with tetracyclines. Many more ARGs and BMRGs were enriched in human-associated E. coli compared with those in chicken- and pork-origin. The quaternary ammonium compounds (QACs) resistance gene qacEΔ1 had a strong correlation with aminoglycoside resistance gene aadA5, folate pathway antagonist resistance gene dfrA17, sul1 and macrolide resistance gene mph(A). The correlation results indicated a significant association between the copper resistance gene cluster pco and the silver resistance gene cluster sil. Coexistence of many resistance genes was observed within the qacEΔ1 gene structure, with qacEΔ1-sul1 being the most common combination. Our findings demonstrated that the epidemiological spread of resistance is affected by a combination of heavy metals, disinfectants and antibiotic use, suggesting that the prevention and control strategies of antimicrobial resistance need to be multifaceted and comprehensive.

Histones and histone variant families in prokaryotes

Histones are important chromatin-organizing proteins in eukaryotes and archaea. They form superhelical structures around which DNA is wrapped. Recent studies have shown that some archaea and bacteria contain alternative histones that exhibit different DNA binding properties, in addition to highly divergent sequences. However, the vast majority of these histones are identified in metagenomes and thus are difficult to study in vivo. The recent revolutionary breakthroughs in computational protein structure prediction by AlphaFold2 and RoseTTAfold allow for unprecedented insights into the potential function and structure of previously uncharacterized proteins. Here, we categorize the prokaryotic histone space into 17 distinct groups based on AlphaFold2 predictions. We identify a superfamily of histones, termed α3 histones, which are common in archaea and present in several bacteria. Importantly, we establish the existence of a large family of histones throughout archaea and in some bacteriophages that, instead of wrapping DNA, bridge DNA, thereby diverging from conventional nucleosomal histones.

Digital Microbe: a genome-informed data integration framework for team science on emerging model organisms

The remarkable pace of genomic data generation is rapidly transforming our understanding of life at the micron scale. Yet this data stream also creates challenges for team science. A single microbe can have multiple versions of genome architecture, functional gene annotations, and gene identifiers; additionally, the lack of mechanisms for collating and preserving advances in this knowledge raises barriers to community coalescence around shared datasets. "Digital Microbes" are frameworks for interoperable and reproducible collaborative science through open source, community-curated data packages built on a (pan)genomic foundation. Housed within an integrative software environment, Digital Microbes ensure real-time alignment of research efforts for collaborative teams and facilitate novel scientific insights as new layers of data are added. Here we describe two Digital Microbes: 1) the heterotrophic marine bacterium Ruegeria pomeroyi DSS-3 with > 100 transcriptomic datasets from lab and field studies, and 2) the pangenome of the cosmopolitan marine heterotroph Alteromonas containing 339 genomes. Examples demonstrate how an integrated framework collating public (pan)genome-informed data can generate novel and reproducible findings.

PTMoreR-enabled cross-species PTM mapping and comparative phosphoproteomics across mammals

To support PTM proteomic analysis and annotation in different species, we developed PTMoreR, a user-friendly tool that considers the surrounding amino acid sequences of PTM sites during BLAST, enabling a motif-centric analysis across species. By controlling sequence window similarity, PTMoreR can map phosphoproteomic results between any two species, perform site-level functional enrichment analysis, and generate kinase-substrate networks. We demonstrate that the majority of real P-sites in mice can be inferred from experimentally derived human P-sites with PTMoreR mapping. Furthermore, the compositions of 129 mammalian phosphoproteomes can also be predicted using PTMoreR. The method also identifies cross-species phosphorylation events that occur on proteins with an increased tendency to respond to the environmental factors. Moreover, the classic kinase motifs can be extracted across mammalian species, offering an evolutionary angle for refining current motifs. PTMoreR supports PTM proteomics in non-human species and facilitates quantitative phosphoproteomic analysis.

Spatial patterns of influenza A virus spread across compartments in commercial swine farms in the United States

AbstractMultiple genetic variants of H1 and H3 influenza A viruses (IAVs) circulate concurrently in US swine farms. Understanding the spatial transmission patterns of IAVs among these farms is crucial for developing effective control strategies and mitigating the emergence of novel IAVs. In this study, we analyzed 1,909 IAV genomic sequences from 785 US swine farms, representing 33 farming systems across 12 states, primarily in the Midwest from 2004 to 2023. Bayesian phylogeographic analyses were performed to identify the dispersal patterns of both H1 and H3 virus genetic lineages and to elucidate their spatial migration patterns within and between different systems. Our results showed that both intra-system and inter-system migrations occurred between the swine farms, with intra-system migrations being more frequent. However, migration rates for H1 and H3 IAVs were similar between intra-system and inter-system migration events. Spatial migration patterns aligned with expected pig movement across different compartments of swine farming systems. Sow-Farms were identified as key sources of viruses, with bi-directional migration observed between these farms and other parts of the system, including Wean-to-Finish and Gilt-Development-Units. High intra-system migration was detected across farms in the same region, while spread to geographically distant intra- and inter- system farms was less frequently. These findings suggest that prioritizing resources towards systems frequently confronting influenza problems and targeting pivotal source farms, such as sow farms, could be an effective strategy for controlling influenza in US commercial swine operations.

Reconstructing an ancient fish: Three-dimensional skeletal restoration of the head of Mawsonia (Sarcopterygii, Actinistia) using CT scan, and an adjusted model for body size estimation in fossil coelacanths

Mawsonia constitutes one of the most conspicuous fossil coelacanth taxa, due to its unique anatomy and possible maximum body size. It typifies Mesozoic coelacanth morphology, before the putative disappearance of the group in the fossil record. In this work, the three-dimensional cranial anatomy and body size estimations of this genus are re-evaluated from a recently described specimen from Upper Jurassic deposits of Uruguay. The 3D restoration was performed directly on the material based on anatomical information provided by the living coelacanth Latimeria and previous two-dimensional restorations of the head of Mawsonia. The montage was then scanned with computed tomography and virtually adjusted to generate an interactive online resource for future anatomical, taxonomic and biomechanical research. In general terms, the model constitutes a tool to improve both the anatomical knowledge of this genus and its comparison with other coelacanths. It also facilitates the evaluation of possible evolutionary trends and the discussion of particular features with potential palaeobiological implications, such as the anterior position of the eye and the development of the pseudomaxillary fold. Regarding the body size, a previous model for body size estimation based on the gular plate was submitted to OLS, RMA, segmented linear and PGLS regressions (including the evaluation of regression statistics, variance analysis, t-tests and residual analysis). The results point to a power relationship between gular and total lengths showing a better support than a simple linear relationship. The new resulting equations were applied to the studied individual and are provided for future estimates. Although an isometric evolutionary growth cannot be rejected with the available evidence, additional models developed with other bones will be necessary to evaluate possible hidden evolutionary allometric trends in this group of fishes, thus avoiding overestimates.

Confirmed Case of Autochthonous Human Babesiosis, Hungary

We report a case of autochthonous human babesiosis in Hungary, confirmed by PCR and partial sequencing of the Babesia spp. 18S rRNA gene. Babesiosis should be considered during the differential diagnosis of febrile illnesses, and peripheral blood smears to detect Babesia spp. should be part of the routine clinical workup.

Contribution of homoeologous exchange to domestication of polyploid Brassica

BACKGROUND

Polyploidy is widely recognized as a significant evolutionary force in the plant kingdom, contributing to the diversification of plants. One of the notable features of allopolyploidy is the occurrence of homoeologous exchange (HE) events between the subgenomes, causing changes in genomic composition, gene expression, and phenotypic variations. However, the role of HE in plant adaptation and domestication remains unclear.

RESULTS

Here we analyze the whole-genome resequencing data from Brassica napus accessions representing the different morphotypes and ecotypes, to investigate the role of HE in domestication. Our findings demonstrate frequent occurrence of HEs in Brassica napus, with substantial HE patterns shared across populations, indicating their potential role in promoting crop domestication. HE events are asymmetric, with the A genome more frequently replacing C genome segments. These events show a preference for specific genomic regions and vary among populations. We also identify candidate genes in HE regions specific to certain populations, which likely contribute to flowering-time diversification across diverse morphotypes and ecotypes. In addition, we assemble a new genome of a swede accession, confirming the HE signals on the genome and their potential involvement in root tuber development. By analyzing HE in another allopolyploid species, Brassica juncea, we characterize a potential broader role of HE in allopolyploid crop domestication.

CONCLUSIONS

Our results provide novel insights into the domestication of polyploid Brassica species and highlight homoeologous exchange as a crucial mechanism for generating variations that are selected for crop improvement in polyploid species.

Population genomics uncovers global distribution, antimicrobial resistance, and virulence genes of the opportunistic pathogen Klebsiella aerogenes

Klebsiella aerogenes is an understudied and clinically important pathogen. We therefore investigate its population structure by genome analysis aligned with metadata. We sequence 130 non-duplicated K. aerogenes clinical isolates and identify two inter-patient transmission events. We then retrieve all publicly available K. aerogenes genomes (n = 1,026, accessed by January 1, 2023) and analyze them with our 130 genomes. We develop a core-genome multi-locus sequence-typing scheme. We find that K. aerogenes is a species complex comprising four phylogroups undergoing evolutionary divergence, likely forming three species. We delineate remarkable clonal diversity and identify three worldwide-distributed carbapenemase-encoding clonal clusters, representing high-risk lineages. We uncover that K. aerogenes has an open genome equipped by a large arsenal of antimicrobial resistance genes. We identify two genetic regions specific for K. aerogenes, encoding a type VI secretion system and flagella/chemotaxis for motility, respectively, both contributing to the virulence. These results provide much-needed insights into the population structure and pan-genomes of K. aerogenes.

Genomic Characterization of a Dog-Mediated Rabies Outbreak in El Pedregal, Arequipa, Peru

Background

Rabies, a re-emerging zoonosis with the highest known human case fatality rate, has been largely absent from Peru, except for endemic circulation in the Puno region on the Bolivian border and re-emergence in Arequipa City in 2015, where it has persisted. In 2021, an outbreak occurred in the rapidly expanding city of El Pedregal near Arequipa, followed by more cases in 2022 after nearly a year of epidemiological silence. While currently under control, questions persist regarding the origin of the El Pedregal outbreak and implications for maintaining rabies control in Peru.

Methods

We sequenced 25 dog rabies virus (RABV) genomes from the El Pedregal outbreak (n=11) and Arequipa City (n=14) from 2021-2023 using Nanopore sequencing in Peru. Historical genomes from Puno (n=4, 2010-2012) and Arequipa (n=5, 2015-2019), were sequenced using an Illumina approach in the UK. In total, 34 RABV genomes were analyzed, including archived and newly obtained samples. The genomes were analyzed phylogenetically to understand the outbreak's context and origins.

Results

Phylogenomic analysis identified two genetic clusters in El Pedregal: 2021 cases stemmed from a single introduction unrelated to Arequipa cases, while the 2022 sequence suggested a new introduction from Arequipa rather than persistence. In relation to canine RABV diversity in Latin America, all new sequences belonged to a new minor clade, Cosmopolitan Am5, sharing relatives from Bolivia, Argentina, and Brazil.

Conclusion

Genomic insights into the El Pedregal outbreak revealed multiple introductions over a 2-year window. Eco-epidemiological conditions, including migratory worker patterns, suggest human-mediated movement drove introductions. Despite outbreak containment, El Pedregal remains at risk of dog-mediated rabies due to ongoing circulation in Arequipa, Puno, and Bolivia. Human-mediated movement of dogs presents a major risk for rabies re-emergence in Peru, jeopardizing regional dog-mediated rabies control. Additional sequence data is needed for comprehensive phylogenetic analyses.

Reducing language barriers, promoting information absorption, and communication using fanyi

ABSTRACT

Interpreting genes of interest is essential for identifying molecular mechanisms, but acquiring such information typically involves tedious manual retrieval. To streamline this process, the fanyi package offers tools to retrieve gene information from sources like National Center for Biotechnology Information (NCBI), significantly enhancing accessibility. Additionally, understanding the latest research advancements and sharing achievements are crucial for junior researchers. However, language barriers often restrict knowledge absorption and career development. To address these challenges, we developed the fanyi package, which leverages artificial intelligence (AI)-driven online translation services to accurately translate among multiple languages. This dual functionality allows researchers to quickly capture and comprehend information, promotes a multilingual environment, and fosters innovation in academic community. Meanwhile, the translation functions are versatile and applicable beyond biomedicine research to other domains as well. The fanyi package is freely available at https://github.com/YuLab-SMU/fanyi .

Total substitution and partial modification of the set of non-ribosomal peptide synthetases clusters lead to pyoverdine diversity in the Pseudomonas fluorescens complex

Pyoverdines are high affinity siderophores produced by most Pseudomonas with a wide role in microbial interspecies interactions. They are primarily composed of a conserved chromophore moiety, an acyl side chain and a peptide backbone which may be highly variable among strains. Upon ferric iron sequestration, pyoverdines are internalized through specialized receptors. The peptide precursor of pyoverdine, termed ferribactin, is synthesized by a set of non-ribosomal peptide synthetase (NRPS) enzymes and further modified by tailoring enzymes. While PvdL, the NRPS responsible for the synthesis of the peptide moiety that derives into the chromophore is conserved, the NRPSs for the peptide backbone are different across fluorescent Pseudomonas. Although the variation of pyoverdine is a widely recognized characteristic within the genus, the evolutionary events associated with the diversity and distribution of this trait remain mostly unknown. This study analyzed the NRPSs clusters for the biosynthesis of the peptide backbone of ferribactin in the genomes of a representative subset of strains of the Pseudomonas fluorescens complex. Bioinformatic analysis of the specificity of adenylation domains of the NRPSs allowed the prediction of 30 different pyoverdine variants. Phylogenetic reconstruction and mapping of the NRPS clusters pinpointed two different general levels of modifications. In the first level, a complete replacement of the set of NRPRs by horizontal transfer occurs. In the second level, the original set of NRPSs is modified through different mechanisms, including partial substitution of the NRPS genes by horizontal transfer, adenylation domain specificity change or NRPS accessory domain gain/loss.

Structural and genetic diversity in the secreted mucins MUC5AC and MUC5B

The secreted mucins MUC5AC and MUC5B are large glycoproteins that play critical defensive roles in pathogen entrapment and mucociliary clearance. Their respective genes contain polymorphic and degenerate protein-coding variable number tandem repeats (VNTRs) that make the loci difficult to investigate with short reads. We characterize the structural diversity of MUC5AC and MUC5B by long-read sequencing and assembly of 206 human and 20 nonhuman primate (NHP) haplotypes. We find that human MUC5B is largely invariant (5,761-5,762 amino acids [aa]); however, seven haplotypes have expanded VNTRs (6,291-7,019 aa). In contrast, 30 allelic variants of MUC5AC encode 16 distinct proteins (5,249-6,325 aa) with cysteine-rich domain and VNTR copy-number variation. We group MUC5AC alleles into three phylogenetic clades: H1 (46%, ∼5,654 aa), H2 (33%, ∼5,742 aa), and H3 (7%, ∼6,325 aa). The two most common human MUC5AC variants are smaller than NHP gene models, suggesting a reduction in protein length during recent human evolution. Linkage disequilibrium and Tajima's D analyses reveal that East Asians carry exceptionally large blocks with an excess of rare variation (p < 0.05) at MUC5AC. To validate this result, we use Locityper for genotyping MUC5AC haplogroups in 2,600 unrelated samples from the 1000 Genomes Project. We observe a signature of positive selection in H1 among East Asians and a depletion of the likely ancestral haplogroup (H3). In Europeans, H3 alleles show an excess of common variation and deviate from Hardy-Weinberg equilibrium (p < 0.05), consistent with heterozygote advantage and balancing selection. This study provides a generalizable strategy to characterize complex protein-coding VNTRs for improved disease associations.

Sacbrood viruses and select Lake Sinai virus variants dominated Apis mellifera colonies symptomatic for European foulbrood

European foulbrood (EFB) is a prevalent disease in the European honey bee (Apis mellifera) in the United States, which can lead to colony decline and collapse. The bacterial components of EFB are well-studied, but the diversity of viral infections within infected colonies has not been explored. In this study, we use meta-transcriptomics sequencing of 12 honey bee hives, symptomatic (+, n = 6) and asymptomatic (-, n = 6) for EFB, to investigate viral infection associated with the disease. We assembled 41 viral genomes, belonging to three families (Iflaviridae, Dicistroviridae, and Sinhaliviridae), all previously reported in honey bees, including Lake Sinai virus, deformed wing virus, sacbrood virus, Black queen cell virus, and Israeli acute paralysis virus. In colonies with severe EFB, we observed a higher occurrence of viral genomes (34 genomes) in contrast to fewer recovered from healthy colonies (seven genomes) and a complete absence of Dicistroviridae genomes.We observed specific Lake Sinai virus clades associated exclusively with EFB + or EFB - colonies, in addition to EFB-afflicted colonies that exhibited an increase in relative abundance of sacbrood viruses. Multivariate analyses highlighted that a combination of site and EFB disease status influenced RNA virome composition, while EFB status alone did not significantly impact it, presenting a challenge for comparisons between colonies kept in different yards. These findings contribute to the understanding of viral dynamics in honey bee colonies compromised by EFB and underscore the need for future investigations to consider viral composition when investigating EFB.IMPORTANCEThis study on the viromes of honey bee colonies affected by European foulbrood (EFB) sheds light on the dynamics of viral populations in bee colonies in the context of a prevalent bacterial brood disease. The identification of distinct Lake Sinai virus and sacbrood virus clades associated with colonies affected by severe EFB suggests a potential connection between viral composition and disease status, emphasizing the need for further investigation into the role of viruses during EFB infection. The observed increase in sacbrood viruses during EFB infection suggests a potential viral dysbiosis, with potential implications for honey bee brood health. These findings contribute valuable insights related to beekeeping practices, offering a foundation for future research aimed at understanding and mitigating the impact of bacterial and viral infection in commercial honey bee operations and the management of EFB.

On the wing: Morphological variation in the osteology of Mediterranean, Near Eastern, and European Anatidae (excluding Anserinae)

Accurate identification of waterfowl bones in archaeological and fossil assemblages has potential to unlock new methods of past environmental reconstruction, as species have differing habitat preferences and migration patterns. Therefore, identifying the presence of avian species with different ecological niches is key to determining past environments and ultimately how prehistoric people responded to climatic and environmental realignments. However, the identification of osteological remains of waterbirds such as ducks to species level is notoriously challenging. We address this by presenting a new two-dimensional geometric morphometric protocol on wing elements from over 20 duck species and test the utility of these shape data for correct species identification. This is an ideal starting point to expand utilization of these types of approaches in avifaunal research and test applicability to an extremely difficult taxonomic group.

Long-term surveillance of SARS-CoV-2 in the school community from Campo Grande, Brazil

BACKGROUND

The COVID-19 pandemic has significantly impacted education systems worldwide, with Brazil being one of the countries with the longest school closures. Over a million children and teenagers have been affected, leading to increased hunger and nutritional deficiencies. This study aimed to implement long-term surveillance of SARS-CoV-2 infections in public and private schools in Campo Grande, Brazil, after returning to in-person classes.

METHODS

The study involved testing and genomic surveillance at 23 public and private schools in Campo Grande, Mato Grosso do Sul, Brazil, from October 18, 2021 to November 21, 2022. The participants eligible for enrollment were students aged 6-17 years and staff members from school institutions. At the time of collection, participants were asked if they had symptoms in the last two weeks. Whole-genome sequencing of SARS-CoV-2 was conducted to identify circulating variants and to compare them with those detected in the municipality. The demographic data and clinical history of the participants were described, and a logistic regression model was used to understand how the RT-qPCR results could be related to different characteristics.

RESULTS

The study included 999 participants, most of whom were women. A total of 85 tests were positive, with an overall positivity rate of 3.2%. The dynamics of case frequency were consistent with those observed in the municipality during the study period. The most common symptoms reported were cough, rhinorrhea, headache, and sore throat. Symptoms were significantly associated with SARS-CoV-2 infection. Eleven lineages were identified in school community samples, with a frequency of occurrence per period similar to that found in the sequences available for the municipality. The most prevalent lineages within the sampling period were BA.2 (59.3%) and BA.5 (29.6%).

CONCLUSIONS

Our findings demonstrate that schools can play a crucial role in epidemiological surveillance, helping trigger rapid responses to pathogens such as SARS-CoV-2. Long-term surveillance can be used to track outbreaks and assess the role of children and adults in transmission. It can also contribute to pandemic preparedness, enabling a rapid response to emergencies, such as COVID-19.

Metabolic reprogramming and heterogeneity during the decidualization process of endometrial stromal cells

The human endometrial decidualization is a transformative event in the pregnant uterus that involves the differentiation of stromal cells into decidual cells. While crucial to the establishment of a successful pregnancy, the metabolic characteristics of decidual cells in vivo remain largely unexplored. Here, we integrated the single-cell RNA sequencing (scRNA-seq) datasets on the endometrium of the menstrual cycle and the maternal-fetal interface in the first trimester to comprehensively decrypt the metabolic characteristics of stromal fibroblast cells. Our results revealed that the differentiation of stromal cells into decidual cells is accompanied by increased amino acid and sphingolipid metabolism. Furthermore, metabolic heterogeneity exists in decidual cells with differentiation maturity disparities. Decidual cells with high metabolism exhibit higher cellular activity and show a strong propensity for signaling. In addition, significant metabolic reprogramming in amino acids and lipids also occurs during the transition from non-pregnancy to pregnancy in the uteri of pigs, cattle, and mice. Our analysis provides comprehensive insights into the dynamic landscape of stromal fibroblast cell metabolism, contributing to our understanding of the metabolism at the molecular dynamics underlying the decidualization process in the human endometrium.

Arginine catabolism is essential to polymyxin dependence in Acinetobacter baumannii

Polymyxins are often the only effective antibiotics against the "Critical" pathogen Acinetobacter baumannii. Worryingly, highly polymyxin-resistant A. baumannii displaying dependence on polymyxins has emerged in the clinic, leading to diagnosis and treatment failures. Here, we report that arginine metabolism is essential for polymyxin-dependent A. baumannii. Specifically, the arginine degradation pathway was significantly altered in polymyxin-dependent strains compared to wild-type strains, with critical metabolites (e.g., L-arginine and L-glutamate) severely depleted and expression of the astABCDE operon significantly increased. Supplementation of arginine increased bacterial metabolic activity and suppressed polymyxin dependence. Deletion of astA, the first gene in the arginine degradation pathway, decreased phosphatidylglycerol and increased phosphatidylethanolamine levels in the outer membrane, thereby reducing the interaction with polymyxins. This study elucidates the molecular mechanism by which arginine metabolism impacts polymyxin dependence in A. baumannii, underscoring its critical role in improving diagnosis and treatment of life-threatening infections caused by "undetectable" polymyxin-dependent A. baumannii.

Maternal gestational diabetes mellitus associates with altered gut microbiome composition and head circumference abnormalities in male offspring

The impact of gestational diabetes mellitus (GDM) on maternal or infant microbiome trajectory remains poorly understood. Utilizing large-scale longitudinal fecal samples from 264 mother-baby dyads, we present the gut microbiome trajectory of the mothers throughout pregnancy and infants during the first year of life. GDM mothers had a distinct microbiome diversity and composition during the gestation period. GDM leaves fingerprints on the infant's gut microbiome, which are confounded by delivery mode. Further, Clostridium species positively correlate with a larger head circumference at month 12 in male offspring but not females. The gut microbiome of GDM mothers with male fetuses displays depleted gut-brain modules, including acetate synthesis I and degradation and glutamate synthesis II. The gut microbiome of female infants of GDM mothers has higher histamine degradation and dopamine degradation. Together, our integrative analysis indicates that GDM affects maternal and infant gut composition, which is associated with sexually dimorphic infant head growth.

MAGinator enables accurate profiling of de novo MAGs with strain-level phylogenies

Metagenomic sequencing has provided great advantages in the characterisation of microbiomes, but currently available analysis tools lack the ability to combine subspecies-level taxonomic resolution and accurate abundance estimation with functional profiling of assembled genomes. To define the microbiome and its associations with human health, improved tools are needed to enable comprehensive understanding of the microbial composition and elucidation of the phylogenetic and functional relationships between the microbes. Here, we present MAGinator, a freely available tool, tailored for profiling of shotgun metagenomics datasets. MAGinator provides de novo identification of subspecies-level microbes and accurate abundance estimates of metagenome-assembled genomes (MAGs). MAGinator utilises the information from both gene- and contig-based methods yielding insight into both taxonomic profiles and the origin of genes and genetic content, used for inference of functional content of each sample by host organism. Additionally, MAGinator facilitates the reconstruction of phylogenetic relationships between the MAGs, providing a framework to identify clade-level differences.

Rapid Adaptation and Interspecific Introgression in the North American Crop Pest Helicoverpa zea

Insect crop pests threaten global food security. This threat is amplified through the spread of nonnative species and through adaptation of native pests to control measures. Adaptations such as pesticide resistance can result from selection on variation within a population, or through gene flow from another population. We investigate these processes in an economically important noctuid crop pest, Helicoverpa zea, which has evolved resistance to a wide range of pesticides. Its sister species Helicoverpa armigera, first detected as an invasive species in Brazil in 2013, introduced the pyrethroid-resistance gene CYP337B3 to South American H. zea via adaptive introgression. To understand whether this could contribute to pesticide resistance in North America, we sequenced 237 H. zea genomes across 10 sample sites. We report H. armigera introgression into the North American H. zea population. Two individuals sampled in Texas in 2019 carry H. armigera haplotypes in a 4 Mbp region containing CYP337B3. Next, we identify signatures of selection in the panmictic population of nonadmixed H. zea, identifying a selective sweep at a second cytochrome P450 gene: CYP333B3. We estimate that its derived allele conferred a ∼5% fitness advantage and show that this estimate explains independently observed rare nonsynonymous CYP333B3 mutations approaching fixation over a ∼20-year period. We also detect putative signatures of selection at a kinesin gene associated with Bt resistance. Overall, we document two mechanisms of rapid adaptation: the introduction of fitness-enhancing alleles through interspecific introgression, and selection on intraspecific variation.

shinyTempSignal: an R shiny application for exploring temporal and other phylogenetic signals

The molecular clock model is fundamental for inferring species divergence times from molecular sequences. However, its direct application may introduce significant biases due to sequencing errors, recombination events, and inaccurately labeled sampling times. Improving accuracy necessitates rigorous quality control measures to identify and remove potentially erroneous sequences. Furthermore, while not all branches of a phylogenetic tree may exhibit a clear temporal signal, specific branches may still adhere to the assumptions, with varying evolutionary rates. Supporting a relaxed molecular clock model better aligns with the complexities of evolution. The root-to-tip regression method has been widely used to analyze the temporal signal in phylogenetic studies and can be generalized for detecting other phylogenetic signals. Despite its utility, there remains a lack of corresponding software implementations for broader applications. To address this gap, we present shinyTempSignal, an interactive web application implemented with the shiny framework, available as an R package and publicly accessible at https://github.com/YuLab-SMU/shinyTempSignal. This tool facilitates the analysis of temporal and other phylogenetic signals under both strict and relaxed models. By extending the root-to-tip regression method to diverse signals, shinyTempSignal helps in the detection of evolving features or traits, thereby laying the foundation for deeper insights and subsequent analyses.

Protocol for fast clonal family inference and analysis from large-scale B cell receptor repertoire sequencing data

The expeditious identification and comprehensive analysis of clonal families from extensive B cell receptor (BCR) repertoire sequencing data are imperative for elucidating the intricacies of B cell immune responses. Here, we introduce a computational pipeline designed to swiftly deduce clonal families from bulk BCR heavy-chain sequencing data, accompanied by a suite of functional modules tailored to streamline post-clustering analysis. The outlined methodology encompasses guidelines for software installation, meticulous data preparation, and the systematic inference and analysis of clonal families. For complete details on the use and execution of this protocol, please refer to Wang et al.1.

Genetic and Phenotypic Evaluation of European Maize Landraces as a Tool for Conservation and Valorization of Agrobiodiversity

The ECPGR European Evaluation Network (EVA) for Maize involves genebanks, research institutions, and private breeding companies from nine countries focusing on the valorization of maize genetic resources across Europe. This study describes a diverse collection of 626 local landraces and traditional varieties of maize (Zea mays L.) from nine European genebanks, including criteria for selection of the collection and its genetic and phenotypic diversity. High-throughput pool genotyping grouped the landraces into nine genetic groups with a threshold of 0.6 admixture, while 277 accessions were designated admixed and likely to have resulted from previous breeding activities. The grouping correlated well with the geographic origins of the collection, also reflecting the various pathways of introduction of maize to Europe. Phenotypic evaluations of 588 accessions for flowering time and plant architecture in multilocation trials over three years confirmed the great diversity within the collection, although phenotypic clusters only partially correlated with the genetic grouping. The EVA approach promotes conservation of genetic resources and opens an opportunity to increase genetic variability for developing improved varieties and populations for farmers, with better adaptation to specific environments and greater tolerance to various stresses. As such, the EVA maize collection provides valuable sources of diversity for facing climate change due to the varieties' local adaptation.

The haplotype-resolved genome of diploid Chrysanthemum indicum unveils new acacetin synthases genes and their evolutionary history

Acacetin, a flavonoid compound, possesses a wide range of pharmacological effects, including antimicrobial, immune regulation, and anticancer effects. Some key steps in its biosynthetic pathway were largely unknown in flowering plants. Here, we present the first haplotype-resolved genome of Chrysanthemum indicum, whose dried flowers contain abundant flavonoids and have been utilized as traditional Chinese medicine. Various phylogenetic analyses revealed almost equal proportion of three tree topologies among three Chrysanthemum species (C. indicum, C. nankingense, and C. lavandulifolium), indicating that frequent gene flow among Chrysanthemum species or incomplete lineage sorting due to rapid speciation might contribute to conflict topologies. The expanded gene families in C. indicum were associated with oxidative functions. Through comprehensive candidate gene screening, we identified five flavonoid O-methyltransferase (FOMT) candidates, which were highly expressed in flowers and whose expressional levels were significantly correlated with the content of acacetin. Further experiments validated two FOMTs (CI02A009970 and CI03A006662) were capable of catalyzing the conversion of apigenin into acacetin, and these two genes are possibly responsible acacetin accumulation in disc florets and young leaves, respectively. Furthermore, combined analyses of ancestral chromosome reconstruction and phylogenetic trees revealed the distinct evolutionary fates of the two validated FOMT genes. Our study provides new insights into the biosynthetic pathway of flavonoid compounds in the Asteraceae family and offers a model for tracing the origin and evolutionary routes of single genes. These findings will facilitate in vitro biosynthetic production of flavonoid compounds through cellular and metabolic engineering and expedite molecular breeding of C. indicum cultivars.

Draft genome sequences of 13 Tenacibaculum maritimum isolates from farmed Norwegian Cyclopterus lumpus (lumpfish) and Scophthalmus maximus (turbot)

Thirteen bacterial isolates of Tenacibaculum maritimum were sequenced and assembled. The strains were isolated from four disease outbreaks in farmed marine fish in Norway. Eight isolates were from Cyclopterus lumpus (lumpfish), and five were from Scophthalmus maximus (turbot). Overall, sequence similarity did not correlate with host species or geographic location.

Unique trajectory of gene family evolution from genomic analysis of nearly all known species in an ancient yeast lineage

Gene gains and losses are a major driver of genome evolution; their precise characterization can provide insights into the origin and diversification of major lineages. Here, we examined gene family evolution of 1,154 genomes from nearly all known species in the medically and technologically important yeast subphylum Saccharomycotina. We found that yeast gene family and genome evolution are distinct from plants, animals, and filamentous ascomycetes and are characterized by small genome sizes and smaller gene numbers but larger gene family sizes. Faster-evolving lineages (FELs) in yeasts experienced significantly higher rates of gene losses-commensurate with a narrowing of metabolic niche breadth-but higher speciation rates than their slower-evolving sister lineages (SELs). Gene families most often lost are those involved in mRNA splicing, carbohydrate metabolism, and cell division and are likely associated with intron loss, metabolic breadth, and non-canonical cell cycle processes. Our results highlight the significant role of gene family contractions in the evolution of yeast metabolism, genome function, and speciation, and suggest that gene family evolutionary trajectories have differed markedly across major eukaryotic lineages.

Regional Comparisons of Sensitivities of Phytophthora citrophthora and P. syringae Causing Citrus Brown Rot in California to Four New and Two Older Fungicides

Isolates of the citrus brown rot pathogens Phytophthora citrophthora and P. syringae from the Inland Empire (IE) and Ventura Co. (VE) regions of southern California were evaluated for their sensitivity to ethaboxam, fluopicolide, mandipropamid, and oxathiapiprolin, and the previously published baselines that were generated for Central Valley (CV) isolates of California were expanded. Fungicides were generally more toxic to CV isolates of both species for all four fungicides. Specific differences were found in the toxicity of ethaboxam to P. syringae where CV isolates on average were 6.8 or 8.2 times more sensitive than those from the VE or IE regions, respectively. Based on the grouping of isolates in an unweighted pair-group method with arithmetic mean (UPGMA) dendrogram, as well as fastStructure analyses and plotting of principal component analyses (PCAs), differences in ethaboxam sensitivity could be related to differences in genetic background of the isolates. Isolates of P. citrophthora from the IE and VE had slightly reduced (i.e., 1.5×) sensitivity to mandipropamid as compared with isolates from the CV and were found on distinct branches in the UPGMA dendrogram. Differences in genetic background of less sensitive isolates within each species indicate that these two phenotypes emerged multiple times independently. IE and VE isolates of both species were sensitive to mefenoxam. Moderate resistance to potassium phosphite (EC50 values of 25 to 75 μg/ml) was present in IE and VE isolates of P. syringae, whereas some IE isolates of P. citrophthora were considered resistant with EC50 values of up to 113.69 μg/ml. Resistance to potassium phosphite did not relate to distinct genotypes.

Mucor thermorhizoides-A New Species from Post-mining Site in Sudety Mountains (Poland)

Mucor representatives are mostly rapidly growing cosmopolitan soil saprotrophs of early diverged Mucoromycotina subphylum. Although this is the most speciose genus within the group, some lineages are still understudied. In this study, new species of Mucor was isolated from the post-mining area in southwestern Poland, where soil chemical composition analysis revealed high concentration of hydrocarbons and heavy metals. Phylogenetic analysis based on multigene phylogeny showed that the new isolate clusters distinctly from other Mucor species as a sister group to Mucor microsporus. New species Mucor thermorhizoides Abramczyk (Mucorales, Mucoromycota) is characterized by the extensive rhizoid production in elevated temperatures and formation of two layers of sporangiophores. It also significantly differs from M. microsporus in the shape of spores and the size of sporangia. M. thermorhizoides was shown to be able to grow in oligotrophic conditions at low temperatures. Together with M. microsporus they represent understudied and highly variable lineage of the Mucor genus.

Expansion of pneumococcal serotype 23F and 14 lineages with genotypic changes in capsule polysaccharide locus and virulence gene profiles post introduction of pneumococcal conjugate vaccine in Blantyre, Malawi

Since the introduction of the 13-valent pneumococcal conjugate vaccine (PCV13) in Malawi in 2011, there has been persistent carriage of vaccine serotype (VT) Streptococcus pneumoniae, despite high vaccine coverage. To determine if there has been a genetic change within the VT capsule polysaccharide (cps) loci since the vaccine's introduction, we compared 1022 whole-genome-sequenced VT isolates from 1998 to 2019. We identified the clonal expansion of a multidrug-resistant, penicillin non-susceptible serotype 23F GPSC14-ST2059 lineage, a serotype 14 GPSC9-ST782 lineage and a novel serotype 14 sequence type GPSC9-ST18728 lineage. Serotype 23F GPSC14-ST2059 had an I253T mutation within the capsule oligosaccharide repeat unit polymerase Wzy protein, which is predicted in silico to alter the protein pocket cavity. Moreover, serotype 23F GPSC14-ST2059 had SNPs in the DNA binding sites for the cps transcriptional repressors CspR and SpxR. Serotype 14 GPSC9-ST782 harbours a non-truncated version of the large repetitive protein (Lrp), containing a Cna protein B-type domain which is also present in proteins associated with infection and colonisation. These emergent lineages also harboured genes associated with antibiotic resistance, and the promotion of colonisation and infection which were absent in other lineages of the same serotype. Together these data suggest that in addition to serotype replacement, modifications of the capsule locus associated with changes in virulence factor expression and antibiotic resistance may promote vaccine escape. In summary, the study highlights that the persistence of vaccine serotype carriage despite high vaccine coverage in Malawi may be partly caused by expansion of VT lineages post-PCV13 rollout.

Endosymbionts of citrus leafminer Phyllocnistis citrella Stainton among different citrus orchards in China

Endosymbionts regulate the behavior of pest species, which could provide insights into their control. The citrus leafminer (Phyllocnistis citrella Stainton) is a widely distributed pest associated with diseases of citrus, especially of young trees. Here, we determined the endosymbiont composition of P. citrella in citrus orchards across China. The resulting dataset comprised average 50,430 high-quality reads for bacterial 16S rRNA V3-V4 regions of endosymbionts from 36 P. citrella larvae sampled from 12 citrus orchards across China. The sequencing depth and sampling size of this dataset were sufficient to reveal most of the endosymbionts of P. citrella. In total, 2,875 bacterial amplicon sequence variants were obtained; taxonomic analysis revealed a total of 372 bacterial genera, most of which were Proteobacteria phylum with Undibacterium being the most abundant genus. This dataset provides the first evidence of P. citrella endosymbionts that could support the development of pest management approaches in citrus orchards.

Exploring the Frozen Armory: Antiphage Defense Systems in Cold-Adapted Bacteria with a Focus on CRISPR-Cas Systems

Our understanding of the antiphage defense system arsenal in bacteria is rapidly expanding, but little is known about its occurrence in cold-adapted bacteria. In this study, we aim to shed light on the prevalence and distribution of antiphage defense systems in cold-adapted bacteria, with a focus on CRISPR-Cas systems. Using bioinformatics tools, Prokaryotic Antiviral Defense LOCator (PADLOC) and CRISPRCasTyper, we mapped the presence and diversity of antiphage defense systems in 938 available genomes of cold-adapted bacteria from diverse habitats. We confirmed that CRISPR-Cas systems are less frequent in cold-adapted bacteria, compared to mesophilic and thermophilic species. In contrast, several antiphage defense systems, such as dXTPases and DRTs, appear to be more frequently compared to temperate bacteria. Additionally, our study provides Cas endonuclease candidates with a potential for further development into cold-active CRISPR-Cas genome editing tools. These candidates could have broad applications in research on cold-adapted organisms. Our study provides a first-time map of antiphage defense systems in cold-adapted bacteria and a detailed overview of CRISPR-Cas diversity.

Common Methods for Phylogenetic Tree Construction and Their Implementation in R

A phylogenetic tree can reflect the evolutionary relationships between species or gene families, and they play a critical role in modern biological research. In this review, we summarize common methods for constructing phylogenetic trees, including distance methods, maximum parsimony, maximum likelihood, Bayesian inference, and tree-integration methods (supermatrix and supertree). Here we discuss the advantages, shortcomings, and applications of each method and offer relevant codes to construct phylogenetic trees from molecular data using packages and algorithms in R. This review aims to provide comprehensive guidance and reference for researchers seeking to construct phylogenetic trees while also promoting further development and innovation in this field. By offering a clear and concise overview of the different methods available, we hope to enable researchers to select the most appropriate approach for their specific research questions and datasets.

Systemic metabolic depletion of gut microbiome undermines responsiveness to melanoma immunotherapy

Immunotherapy has proven to be a boon for patients battling metastatic melanoma, significantly improving their clinical condition and overall quality of life. A compelling link between the composition of the gut microbiome and the efficacy of immunotherapy has been established in both animal models and human patients. However, the precise biological mechanisms by which gut microbes influence treatment outcomes remain poorly understood. Using a robust dataset of 680 fecal metagenomes from melanoma patients, a detailed catalog of metagenome-assembled genomes (MAGs) was constructed to explore the compositional and functional properties of the gut microbiome. Our study uncovered significant findings that deepen the understanding of the intricate relationship between gut microbes and the efficacy of melanoma immunotherapy. In particular, we discovered the specific metagenomic profile of patients with favorable treatment outcomes, characterized by a prevalence of MAGs with increased overall metabolic potential and proficiency in polysaccharide utilization, along with those responsible for cobalamin and amino acid production. Furthermore, our investigation of the biosynthetic pathways of short-chain fatty acids, known for their immunomodulatory role, revealed a differential abundance of these pathways among the specific MAGs. Among others, the cobalamin-dependent Wood-Ljungdahl pathway of acetate synthesis was directly associated with responsiveness to melanoma immunotherapy.

Balancing selection and the functional effects of shared polymorphism in cryptic Daphnia species

The patterns of genetic variation within and between related taxa represent the genetic history of a species. Shared polymorphisms, loci with identical alleles across species, are of unique interest as they may represent cases of ancient selection maintaining functional variation post-speciation. In this study, we investigate the abundance of shared polymorphism in the Daphnia pulex species complex. We test whether shared mutations are consistent with the action of balancing selection or alternative hypotheses such as hybridization, incomplete lineage sorting, or convergent evolution. We analyzed over 2,000 genomes from North American and European D. pulex and several outgroup species to examine the prevalence and distribution of shared alleles between the focal species pair, North American and European D. pulex. We show that while North American and European D. pulex diverged over ten million years ago, they retained tens of thousands of shared alleles. We found that the number of shared polymorphisms between North American and European D. pulex cannot be explained by hybridization or incomplete lineage sorting alone. Instead, we show that most shared polymorphisms could be the product of convergent evolution, that a limited number appear to be old trans-specific polymorphisms, and that balancing selection is affecting young and ancient mutations alike. Finally, we provide evidence that a blue wavelength opsin gene with trans-specific polymorphisms has functional effects on behavior and fitness in the wild. Ultimately, our findings provide insights into the genetic basis of adaptation and the maintenance of genetic diversity between species.

Potential routes of plastics biotransformation involving novel plastizymes revealed by global multi-omic analysis of plastic associated microbes

Despite increasing efforts across various disciplines, the fate, transport, and impact of synthetic plastics on the environment and public health remain poorly understood. To better elucidate the microbial ecology of plastic waste and its potential for biotransformation, we conducted a large-scale analysis of all publicly available meta-omic studies investigating plastics (n = 27) in the environment. Notably, we observed low prevalence of known plastic degraders throughout most environments, except for substantial enrichment in riverine systems. This indicates rivers may be a highly promising environment for discovery of novel plastic bioremediation products. Ocean samples associated with degrading plastics showed clear differentiation from non-degrading polymers, showing enrichment of novel putative biodegrading taxa in the degraded samples. Regarding plastisphere pathogenicity, we observed significant enrichment of antimicrobial resistance genes on plastics but not of virulence factors. Additionally, we report a co-occurrence network analysis of 10 + million proteins associated with the plastisphere. This analysis revealed a localized sub-region enriched with known and putative plastizymes-these may be useful for deeper investigation of nature's ability to biodegrade man-made plastics. Finally, the combined data from our meta-analysis was used to construct a publicly available database, the Plastics Meta-omic Database (PMDB)-accessible at plasticmdb.org. These data should aid in the integrated exploration of the microbial plastisphere and facilitate research efforts investigating the fate and bioremediation potential of environmental plastic waste.

Natural mating ability is associated with gut microbiota composition and function in captive male giant pandas

The issue of poor sexual performance of some male giant pandas seriously impairs the growth and the genetic diversity of the captive population, yet there is still no clear understanding of the cause of the loss of this ability and its underlying mechanism. In this study, we analyzed the gut microbiota and its function in 72 fecal samples obtained from 20 captive male giant pandas, with an equal allocation between individuals capable and incapable of natural mating. Additionally, we investigated fecal hormone levels and behavioral differences between the two groups. A correlation analysis was then conducted among these factors to explore the influencing factors of their natural mating ability. The results showed significant differences in the composition of gut microbiota between the two groups of male pandas. The capable group had significantly higher abundance of Clostridium sensu stricto 1 (p adjusted = .0021, GLMM), which was positively correlated with fatty acid degradation and two-component system functions (Spearman, p adjusted < .05). Additionally, the capable group showed higher gene abundance in gut microbiota function including purine and pyrimidine metabolism and galactose metabolism, as well as pathways related to biological processes such as ribosome and homologous recombination (DEseq2, p adjusted < .05). We found no significant differences in fecal cortisol and testosterone levels between the two groups, and no difference was found in their behavior either. Our study provides a theoretical and practical basis for further studying the behavioral degradation mechanisms of giant pandas and other endangered mammal species.

Despite structural identity, ace-1 heterogenous duplication resistance alleles are quite diverse in Anopheles mosquitoes

Anopheles gambiae s.l. has been the target of intense insecticide treatment since the mid-20th century to try and control malaria. A substitution in the ace-1 locus has been rapidly selected for, allowing resistance to organophosphate and carbamate insecticides. Since then, two types of duplication of the ace-1 locus have been found in An. gambiae s.l. populations: homogeneous duplications that are composed of several resistance copies, or heterogeneous duplications that contain both resistance and susceptible copies. The substitution induces a trade-off between resistance in the presence of insecticides and disadvantages in their absence: the heterogeneous duplications allow the fixation of the intermediate heterozygote phenotype. So far, a single heterogeneous duplication has been described in An. gambiae s.l. populations (in contrast with the multiple duplicated alleles found in Culex pipiens mosquitoes). We used a new approach, combining long and short-read sequencing with Sanger sequencing to precisely identify and describe at least nine different heterogeneous duplications, in two populations of An. gambiae s.l. We show that these alleles share the same structure as the previously identified heterogeneous and homogeneous duplications, namely 203-kb tandem amplifications with conserved breakpoints. Our study sheds new light on the origin and maintenance of these alleles in An. gambiae s.l. populations, and their role in mosquito adaptation.

Diversity and Clonality in Populations of Phytophthora citrophthora and P. syringae Causing Brown Rot of Citrus in California

Phytophthora citrophthora and P. syringae are currently the primary causal organisms of brown rot of citrus fruits in California. To possibly find an explanation for the prevalence of the previously minor species P. syringae, we determined the population structures of both pathogens in California using next-generation sequencing and population genomics analyses. Whole-genome sequencing and aligning with newly assembled reference genomes identified 972,266 variants in 132 isolates of P. citrophthora and 422,208 variants in 154 isolates (including 24 from noncitrus tree crops) of P. syringae originating from three major growing regions. The resulting data sets were visualized using principal component analysis, discriminant analysis of principal components, unweighted pair-group method with arithmetic mean dendrograms, fastStructure, and minimum spanning networks, and we obtained the index of association, diversity summary statistics, and genetic distance statistics values GST, G''ST, and Jost's D. Subpopulations of both species were mostly defined by their geographic origin indicating restricted dispersal of inoculum. Except for five isolates, the population structure of P. citrophthora (that is heterothallic and unlikely to reproduce sexually) was clonal to semi-clonal, with very little genetic diversity within and among subgroups. In contrast, the population structure of P. syringae was also clonal to semi-clonal, but isolates were placed into four main clusters of much higher diversity. Clonality in both species can be explained by a high level of asexual reproduction. The higher diversity in the homothallic P. syringae is likely due to commonly occurring sexual reproduction. One distinct cluster of P. syringae consisted solely of isolates from noncitrus hosts; therefore, the origin of P. syringae in citrus could not be resolved.

GC Content Across Insect Genomes: Phylogenetic Patterns, Causes and Consequences

The proportions of A:T and G:C nucleotide pairs are often unequal and can vary greatly between animal species and along chromosomes. The causes and consequences of this variation are incompletely understood. The recent release of high-quality genome sequences from the Darwin Tree of Life and other large-scale genome projects provides an opportunity for GC heterogeneity to be compared across a large number of insect species. Here we analyse GC content along chromosomes, and within protein-coding genes and codons, of 150 insect species from four holometabolous orders: Coleoptera, Diptera, Hymenoptera, and Lepidoptera. We find that protein-coding sequences have higher GC content than the genome average, and that Lepidoptera generally have higher GC content than the other three insect orders examined. GC content is higher in small chromosomes in most Lepidoptera species, but this pattern is less consistent in other orders. GC content also increases towards subtelomeric regions within protein-coding genes in Diptera, Coleoptera and Lepidoptera. Two species of Diptera, Bombylius major and B. discolor, have very atypical genomes with ubiquitous increase in AT content, especially at third codon positions. Despite dramatic AT-biased codon usage, we find no evidence that this has driven divergent protein evolution. We argue that the GC landscape of Lepidoptera, Diptera and Coleoptera genomes is influenced by GC-biased gene conversion, strongest in Lepidoptera, with some outlier taxa affected drastically by counteracting processes.

Structural and genetic diversity in the secreted mucins, MUC5AC and MUC5B

The secreted mucins MUC5AC and MUC5B play critical defensive roles in airway pathogen entrapment and mucociliary clearance by encoding large glycoproteins with variable number tandem repeats (VNTRs). These polymorphic and degenerate protein coding VNTRs make the loci difficult to investigate with short reads. We characterize the structural diversity of MUC5AC and MUC5B by long-read sequencing and assembly of 206 human and 20 nonhuman primate (NHP) haplotypes. We find that human MUC5B is largely invariant (5761-5762aa); however, seven haplotypes have expanded VNTRs (6291-7019aa). In contrast, 30 allelic variants of MUC5AC encode 16 distinct proteins (5249-6325aa) with cysteine-rich domain and VNTR copy number variation. We grouped MUC5AC alleles into three phylogenetic clades: H1 (46%, ~5654aa), H2 (33%, ~5742aa), and H3 (7%, ~6325aa). The two most common human MUC5AC variants are smaller than NHP gene models, suggesting a reduction in protein length during recent human evolution. Linkage disequilibrium (LD) and Tajima's D analyses reveal that East Asians carry exceptionally large MUC5AC LD blocks with an excess of rare variation (p<0.05). To validate this result, we used Locityper for genotyping MUC5AC haplogroups in 2,600 unrelated samples from the 1000 Genomes Project. We observed signatures of positive selection in H1 and H2 among East Asians and a depletion of the likely ancestral haplogroup (H3). In Africans and Europeans, H3 alleles show an excess of common variation and deviate from Hardy-Weinberg equilibrium, consistent with heterozygote advantage and balancing selection. This study provides a generalizable strategy to characterize complex protein coding VNTRs for improved disease associations.

Population-level transposable element expression dynamics influence trait evolution in a fungal crop pathogen

The rapid adaptive evolution of microbes is driven by strong selection pressure acting on genetic variation. How adaptive genetic variation is generated within species and how such variation influences phenotypic trait expression is often not well understood though. We focused on the recent activity of transposable elements (TEs) using deep population genomics and transcriptomics analyses of a fungal plant pathogen with a highly active content of TEs in the genome. Zymoseptoria tritici causes one of the most damaging diseases on wheat, with recent adaptation to the host and environment being facilitated by TE-associated mutations. We obtained genomic and RNA-sequencing data from 146 isolates collected from a single wheat field. We established a genome-wide map of TE insertion polymorphisms in the population by analyzing recent TE insertions among individuals. We quantified the locus-specific transcription of individual TE copies and found considerable population variation at individual TE loci in the population. About 20% of all TE copies show transcription in the genome suggesting that genomic defenses such as repressive epigenetic marks and repeat-induced polymorphisms are at least partially ineffective at preventing the proliferation of TEs in the genome. A quarter of recent TE insertions are associated with expression variation of neighboring genes providing broad potential to influence trait expression. We indeed found that TE insertions are likely responsible for variation in virulence on the host and potentially diverse components of secondary metabolite production. Our large-scale transcriptomics study emphasizes how TE-derived polymorphisms segregate even in individual microbial populations and can broadly underpin trait variation in pathogens.IMPORTANCEPathogens can rapidly adapt to new hosts, antimicrobials, or changes in the environment. Adaptation arises often from mutations in the genome; however, how such variation is generated remains poorly understood. We investigated the most dynamic regions of the genome of Zymoseptoria tritici, a major fungal pathogen of wheat. We focused on the transcription of transposable elements. A large proportion of the transposable elements not only show signatures of potential activity but are also variable within a single population of the pathogen. We find that this variation in activity is likely influencing many important traits of the pathogen. Hence, our work provides insights into how a microbial species can adapt over the shortest time periods based on the activity of transposable elements.

Transmission and dynamics of mother-infant gut viruses during pregnancy and early life

Early development of the gut ecosystem is crucial for lifelong health. While infant gut bacterial communities have been studied extensively, the infant gut virome remains under-explored. To study the development of the infant gut virome over time and the factors that shape it, we longitudinally assess the composition of gut viruses and their bacterial hosts in 30 women during and after pregnancy and in their 32 infants during their first year of life. Using shotgun metagenomic sequencing applied to dsDNA extracted from Virus-Like Particles (VLPs) and bacteria, we generate 205 VLP metaviromes and 322 total metagenomes. With this data, we show that while the maternal gut virome composition remains stable during late pregnancy and after birth, the infant gut virome is dynamic in the first year of life. Notably, infant gut viromes contain a higher abundance of active temperate phages compared to maternal gut viromes, which decreases over the first year of life. Moreover, we show that the feeding mode and place of delivery influence the gut virome composition of infants. Lastly, we provide evidence of co-transmission of viral and bacterial strains from mothers to infants, demonstrating that infants acquire some of their virome from their mother's gut.

Whole Genome Duplication Events Likely Contributed to the Aquatic Adaptive Evolution of Parkerioideae

As the only aquatic lineage of Pteridaceae, Parkerioideae is distinct from many xeric-adapted species of the family and consists of the freshwater Ceratopteris species and the only mangrove ferns from the genus Acrostichum. Previous studies have shown that whole genome duplication (WGD) has occurred in Parkerioideae at least once and may have played a role in their adaptive evolution; however, more in-depth research regarding this is still required. In this study, comparative and evolutionary transcriptomics analyses were carried out to identify WGDs and explore their roles in the environmental adaptation of Parkerioideae. Three putative WGD events were identified within Parkerioideae, two of which were specific to Ceratopteris and Acrostichum, respectively. The functional enrichment analysis indicated that the lineage-specific WGD events have played a role in the adaptation of Parkerioideae to the low oxygen concentrations of aquatic habitats, as well as different aquatic environments of Ceratopteris and Acrostichum, such as the adaptation of Ceratopteris to reduced light levels and the adaptation of Acrostichum to high salinity. Positive selection analysis further provided evidence that the putative WGD events may have facilitated the adaptation of Parkerioideae to changes in habitat. Moreover, the gene family analysis indicated that the plasma membrane H+-ATPase (AHA), vacuolar H+-ATPase (VHA), and suppressor of K+ transport growth defect 1 (SKD1) may have been involved in the high salinity adaptation of Acrostichum. Our study provides new insights into the evolution and adaptations of Parkerioideae in different aquatic environments.

A phylogenomic perspective on interspecific competition

Evolutionary processes may have substantial impacts on community assembly, but evidence for phylogenetic relatedness as a determinant of interspecific interaction strength remains mixed. In this perspective, we consider a possible role for discordance between gene trees and species trees in the interpretation of phylogenetic signal in studies of community ecology. Modern genomic data show that the evolutionary histories of many taxa are better described by a patchwork of histories that vary along the genome rather than a single species tree. If a subset of genomic loci harbour trait-related genetic variation, then the phylogeny at these loci may be more informative of interspecific trait differences than the genome background. We develop a simple method to detect loci harbouring phylogenetic signal and demonstrate its application through a proof-of-principle analysis of Penicillium genomes and pairwise interaction strength. Our results show that phylogenetic signal that may be masked genome-wide could be detectable using phylogenomic techniques and may provide a window into the genetic basis for interspecific interactions.

TreeViewer: Flexible, modular software to visualise and manipulate phylogenetic trees

Phylogenetic trees illustrate evolutionary relationships between taxa or genes. Tree figures are crucial when presenting results and data, and by creating clear and effective plots, researchers can describe many kinds of evolutionary patterns. However, producing tree plots can be a time-consuming task, especially as multiple different programs are often needed to adjust and illustrate all data associated with a tree. We present TreeViewer, a new software to draw phylogenetic trees. TreeViewer is flexible, modular, and user-friendly. Plots are produced as the result of a user-defined pipeline, which can be finely customised and easily applied to different trees. Every feature of the program is documented and easily accessible, either in the online manual or within the program's interface. We show how TreeViewer can be used to produce publication-ready figures, saving time by not requiring additional graphical post-processing tools. TreeViewer is freely available for Windows, macOS, and Linux operating systems and distributed under an AGPLv3 licence from https://treeviewer.org. It has a graphical user interface (GUI), as well as a command-line interface, which is useful to work with very large trees and for automated pipelines. A detailed user manual with examples and tutorials is also available. TreeViewer is mainly aimed at users wishing to produce highly customised, publication-quality tree figures using a single GUI software tool. Compared to other GUI tools, TreeViewer offers a richer feature set and a finer degree of customisation. Compared to command-line-based tools and software libraries, TreeViewer's graphical interface is more accessible. The flexibility of TreeViewer's approach to phylogenetic tree plotting enables the program to produce a wide variety of publication-ready figures. Users are encouraged to create their own custom modules to expand the functionalities of the program. This sets the scene for an ever-expanding and ever-adapting software framework that can easily adjust to respond to new challenges.

Longitudinal quantification of Bifidobacterium longum subsp. infantis reveals late colonization in the infant gut independent of maternal milk HMO composition

Breast milk contains human milk oligosaccharides (HMOs) that cannot be digested by infants, yet nourish their developing gut microbiome. While Bifidobacterium are the best-known utilizers of individual HMOs, a longitudinal study examining the evolving microbial community at high-resolution coupled with mothers' milk HMO composition is lacking. Here, we developed a high-throughput method to quantify Bifidobacterium longum subsp. infantis (BL. infantis), a proficient HMO-utilizer, and applied it to a longitudinal cohort consisting of 21 mother-infant dyads. We observed substantial changes in the infant gut microbiome over the course of several months, while the HMO composition in mothers' milk remained relatively stable. Although Bifidobacterium species significantly influenced sample variation, no specific HMOs correlated with Bifidobacterium species abundance. Surprisingly, we found that BL. infantis colonization began late in the breastfeeding period both in our cohort and in other geographic locations, highlighting the importance of focusing on BL. infantis dynamics in the infant gut.

Physiology of γ-aminobutyric acid production by Akkermansia muciniphila

Gut bacteria hold the potential to produce a broad range of metabolites that can modulate human functions, including molecules with neuroactive potential. One such molecule is γ-aminobutyric acid (GABA), the main inhibitory neurotransmitter of the central nervous system in animals. Metagenomic analyses suggest that the genomes of many gut bacteria encode glutamate decarboxylase (GAD), the enzyme that catalyzes GABA production. The genome of Akkermansia muciniphila, a mucin specialist and potential next-generation probiotic from the human gut, is predicted to encode GAD, suggesting a contributing role in GABA production in the human gut. In this study, A. muciniphila was grown in batch cultures with and without pH control. In both experiments, A. muciniphila was found to produce GABA as a response to acid (pH <5.5), although only when GABA precursors, either glutamate or glutamine, were present in the medium. Proteomic analysis comparing A. muciniphila grown with and without precursors at pH 4 did not show a difference in GAD expression, suggesting that it is expressed regardless of the presence of GABA precursors. To further investigate the function of A. muciniphila GAD, we heterologously expressed the gad gene (encoded by locus tag Amuc_0372) with a His tag in Escherichia coli and purified the GAD protein. Enzyme assays showed GAD activity in a pH range between 4 and 6, with the highest specific activity at pH 5 of 144 ± 16 µM GABA/min/mg. Overall, our results demonstrate the ability of A. muciniphila to produce GABA as an acid response and unravel the conditions under which GABA production in A. muciniphila occurs.IMPORTANCEAkkermansia muciniphila is considered to be a beneficial bacterium from the human gut, but the exact mechanisms by which A. muciniphila influences its host are not yet fully understood. To this end, it is important to identify which metabolites are produced and consumed by A. muciniphila that may contribute to a healthy gut. In the present study, we demonstrate the ability of A. muciniphila to produce γ-aminobutyric acid (GABA) when grown in an acidic environment, which often occurs in the gut. GABA is the major inhibitory neurotransmitter in the central nervous system and is present in the human gut. For this reason, it is considered an important bacterial metabolite. Our finding that A. muciniphila produces GABA in acidic environments adds to the growing body of understanding of its relationship with host health and provides an explanation on how it can survive acid stress in the human gut.

Thermodynamic restrictions determine ammonia tolerance of functional floras during anaerobic digestion

Ammonia inhibition is a major challenge in anaerobic digestion processes, affecting the activity and performance of functional floras, including syntrophic butyrate oxidation (FSBO), syntrophic propionate oxidation (FSPO), acetoclastic methanogenesis (FAMs), and hydrogenotrophic methanogenesis (FHMs). FHMs was more tolerant to ammonia, with a half maximal inhibitory concentration (IC50) of 18.80 g/L, followed by FSBO (IC50 = 14.26 g/L) and FSPO (IC50 = 10.47 g/L), and FAMs was the most sensitive to ammonia with the lowest IC50 (1.74 g/L). The order of ammonia tolerance (IC50) was found to be correlated with the Gibbs free energy of metabolic reactions in each functional flora. The results suggested the functional flora with higher energy availability (exergonic capacity) exhibited enhanced ammonia resistance ability. These findings provide insights into the thermodynamic restrictions and ammonia tolerance mechanisms of functional floras, which can guide the optimization and operation of anaerobic digestion systems for efficient methane production.

Host genetic regulation of human gut microbial structural variation

Although the impact of host genetics on gut microbial diversity and the abundance of specific taxa is well established1-6, little is known about how host genetics regulates the genetic diversity of gut microorganisms. Here we conducted a meta-analysis of associations between human genetic variation and gut microbial structural variation in 9,015 individuals from four Dutch cohorts. Strikingly, the presence rate of a structural variation segment in Faecalibacterium prausnitzii that harbours an N-acetylgalactosamine (GalNAc) utilization gene cluster is higher in individuals who secrete the type A oligosaccharide antigen terminating in GalNAc, a feature that is jointly determined by human ABO and FUT2 genotypes, and we could replicate this association in a Tanzanian cohort. In vitro experiments demonstrated that GalNAc can be used as the sole carbohydrate source for F. prausnitzii strains that carry the GalNAc-metabolizing pathway. Further in silico and in vitro studies demonstrated that other ABO-associated species can also utilize GalNAc, particularly Collinsella aerofaciens. The GalNAc utilization genes are also associated with the host's cardiometabolic health, particularly in individuals with mucosal A-antigen. Together, the findings of our study demonstrate that genetic associations across the human genome and bacterial metagenome can provide functional insights into the reciprocal host-microbiome relationship.