Adapterama I: universal stubs and primers for 384 unique dual-indexed or 147,456 combinatorially-indexed Illumina libraries (iTru & iNext)

Gut Microbiome Compositional and Functional Features Associate with Alzheimer's Disease Pathology

BACKGROUND

The gut microbiome is a potentially modifiable factor in Alzheimer's disease (AD); however, understanding of its composition and function regarding AD pathology is limited.

METHODS

Shallow-shotgun metagenomic data was used to analyze fecal microbiome from participants enrolled in the Wisconsin Microbiome in Alzheimer's Risk Study, leveraging clinical data and cerebrospinal fluid (CSF) biomarkers. Differential abundance and ordinary least squares regression analyses were performed to find differentially abundant gut microbiome features and their associations with CSF biomarkers of AD and related pathologies.

RESULTS

Gut microbiome composition and function differed between people with AD and cognitively unimpaired individuals. The compositional difference was replicated in an independent cohort. Differentially abundant gut microbiome features were associated with CSF biomarkers of AD and related pathologies.

DISCUSSION

These findings enhance our understanding of alterations in gut microbial composition and function in AD, and suggest that gut microbes and their pathways are linked to AD pathology.

Sea anemone (Anthozoa, Actiniaria) diversity in Mo'orea (French Polynesia)

Sea anemones (Order Actiniaria) are a diverse group of marine invertebrates ubiquitous across marine ecosystems. Despite their wide distribution and success, a knowledge gap persists in our understanding of their diversity within tropical systems, owed to sampling bias of larger and more charismatic species overshadowing cryptic lineages. This study aims to delineate the sea anemone diversity in Mo'orea (French Polynesia) with the use of a dataset from the Mo'orea Biocode's "BioBlitz" initiative, which prioritized the sampling of more cryptic and understudied taxa. Implementing a target enrichment approach, we integrate 71 newly sequenced samples into an expansive phylogenetic framework and contextualize Mo'orea's diversity within global distribution patterns of sea anemones. Our analysis corroborates the presence of several previously documented sea anemones in French Polynesia and identifies for the first time the occurrence of members of genera Andvakia and Aiptasiomorpha. This research unveils the diverse sea anemone ecosystem in Mo'orea, spotlighting the area's ecological significance and emphasizing the need for continued exploration. Our methodology, encompassing a broad BLAST search coupled with phylogenetic analysis, proved to be a practical and effective approach for overcoming the limitations posed by the lack of comprehensive sequence data for sea anemones. We discuss the merits and limitations of current molecular methodologies and stress the importance of further research into lesser-studied marine organisms like sea anemones. Our work sets a precedent for future phylogenetic studies stemming from BioBlitz endeavors.

Phylogenomics and biogeography of the small carpenter bees (Apidae: Xylocopinae: Ceratina)

Small carpenter bees in the genus Ceratina are behaviourally diverse, species-rich, and cosmopolitan, with over 370 species and a range including all continents except Antarctica. Here, we present the first comprehensive phylogeny of the genus based on ultraconserved element (UCE) phylogenomic data, covering a total of 185 ingroup specimens representing 22 of the 25 current subgenera. Our results support most recognized subgenera as natural groups, but we also highlight several groups in need of taxonomic revision - particularly the nominate subgenus Ceratina sensu stricto - and several clades that likely need to be described as new subgenera. In addition to phylogeny, we explore the evolutionary history of Ceratina through divergence time estimation and biogeographic reconstruction. Our findings suggest that Ceratinini split from its sister tribe Allodapini about 72 million years ago. The common ancestor of Ceratina emerged in the Afrotropical realm approximately 42 million years ago, near the Middle Eocene Climatic Optimum. Multiple subsequent dispersal events led to the present cosmopolitan distribution of Ceratina, with the majority of transitions occurring between the Afrotropics, Indomalaya, and the Palearctic. Additional movements also led to the arrival of Ceratina in Madagascar, Australasia, and a single colonization of the Americas. Dispersal events were asymmetrical overall, with temperate regions primarily acting as destinations for migrations from tropical source regions.

Testing ultraconserved elements (UCEs) for phylogenetic inference across bivalves (Mollusca: Bivalvia)

Bivalves constitute an important resource for fisheries and as cultural objects. Bivalve phylogenetics has had a long tradition using both morphological and molecular characters, and genomic resources are available for a good number of commercially important species. However, relationships among bivalve families have been unstable and major conflicting results exist between mitogenomics and results based on Sanger-based amplicon sequencing or phylotranscriptomics. Here we design and test an ultraconserved elements probe set for the class Bivalvia with the aim to use hundreds of loci without the need to sequence full genomes or transcriptomes, which are expensive and complex to analyze, and to open bivalve phylogenetics to museum specimens. Our probe set successfully captured 1,513 UCEs for a total of 263,800 bp with an average length of 174.59 ± 3.44 per UCE (ranging from 28 to 842 bp). Phylogenetic testing of this UCE probe set across Bivalvia and within the family Donacidae using different data matrices and methods for phylogenetic inference shows promising results at multiple taxonomic levels. In addition, our probe set was able to capture large numbers of UCEs for museum specimens collected before 1900 and from DNAs properly stored, of which many museums and laboratories are well stocked. Overall, this constitutes a novel and useful resource for bivalve phylogenetics.

Selection and Gene Duplication Associated With High-Elevation Diversification in Pristimantis, the Largest Terrestrial Vertebrate Genus

The genus Pristimantis diversified in the tropical Andes mountains and is the most speciose genus of terrestrial vertebrates. Pristimantis are notable among frogs in that they thrive at high elevations (>2,000 m) and are direct developers without a tadpole stage. Despite their ecological significance, little is known about the genetic and physiological traits enabling their success. We conducted transcriptomic analysis on seven Pristimantis species sampled across elevations in the Ecuadorean Andes to explore three hypotheses for their success: (i) unique genes are under selection relative to all other frogs, (ii) common selection occurs across all direct developers, or (iii) common selection occurs across all high-elevation frog clades. Comparative analysis with 34 frog species revealed unique positive selection in Pristimantis genes related to aerobic respiration, hemostasis, signaling, cellular transportation of proteins and ions, and immunity. Additionally, we detected positive selection across all direct developers for genes associated with oxygenase activity and metal ion binding. While many genes under selection in Pristimantis were not positively selected in other high-elevation frog species, we identified some shared genes and pathways linked to lipid metabolism, innate immunity, and cellular redox processes. We observed more positive selection in duplicated- versus single-copy genes, while relaxed purifying selection was prevalent in single-copy genes. Notably, copy number of an innate immunity complement gene was positively correlated with Pristimantis species elevation. Our findings contribute novel insights into the genetic basis of adaptation in Pristimantis and provide a foundation for future studies on the evolutionary mechanisms leading to direct development and coping with high elevations.

Functional trait mismatch between native and introduced bee pollinators servicing a global fruit crop

BACKGROUND

Understanding connections between biodiversity and ecosystem services can be enhanced by shifting focus from species richness to functional trait-based approaches, that when paired with comparative phylogenetic methods can provide even deeper insights. We investigated the functional ecology and phylogenetic diversity of pollination services provided by hymenopteran insects visiting apple flowers in orchards surrounded by either 'natural' or 'disturbed' landscapes in New South Wales, Australia. We assessed whether morphological and behavioural traits (hairiness, body size, glossa length, pollen load purity, and probability of loose pollen) exhibited non-random phylogenetic patterns. Then, explored whether bees, the primary pollinators in this system, filled unique or overlapping functional entities (FEs). For each landscape, we calculated phylogenetic diversity and used FEs to assess functional richness, evenness, and diversion.

RESULTS

A phylogenomic matrix based on ultraconserved elements (UCEs; 1,382,620 bp from 1,969 loci) was used to infer a fully-resolved and well-supported maximum likelihood phylogeny for 48 hymenopteran morphospecies. There was no significant difference in species richness between landscape categories. Pollinator communities at natural sites had higher phylogenetic complexity (X = 2.37) and functional divergence (x̄ = 0.74 ± 0.02 s.e.) than disturbed sites (X = 1.65 and x̄ = 0.6 ± 0.01 s.e.). Hairiness showed significant phylogenetic clustering (K = 0.94), whereas body size, glossa length, and loose pollen showed weaker non-random phylogenetic patterns (K between 0.3-0.5). Pollen load purity showed no association with phylogeny. The assemblage of 17 bee morphospecies comprised nine FEs: eight FEs consisted of native bees with three containing 65% of all native bee taxa. The introduced honey bee (Apis mellifera) occupied a unique FE, likely due to its different evolutionary history. Both landscape types supported six FEs each with three overlapping: two native bee FEs and the honey bee FE.

CONCLUSIONS

Bee hairiness was the only functional trait to exhibit demonstrable phylogenetic signal. Despite differences in species richness, and functional and phylogenetic diversity between orchard landscape types, both maintained equal bee FE numbers. While no native bee taxon was analogous to the honey bee FE, four native bee FEs shared the same hairiness level as honey bees. Health threats to honey bee populations in Australia will likely disrupt pollination services to apple, and other pollination-dependent food crops, given the low level of functional redundancy within the investigated pollinator assemblages.

UCE-based phylogenomics of the lepidopteran endoparasitoid wasp subfamily Rogadinae (Hymenoptera: Braconidae) unveils a new Neotropical tribe

During the past two decades, the phylogenetic relationships and higher-level classification of the subfamily Rogadinae have received relevant contributions based on Sanger, mitogenome and genome-wide nuclear DNA sequence data. These studies have helped to update the circumscription and tribal classification of this subfamily, with six tribes currently recognised (Aleiodini, Betylobraconini, Clinocentrini, Rogadini, Stiropiini and Yeliconini). The tribal relationships within Rogadinae, however, are yet to be fully resolved, including the status of tribe Facitorini, previously regarded as betylobraconine, with respect to the members of Yeliconini. We conducted a phylogenomic analysis among the tribes of Rogadinae based on genomic ultraconserved element (UCE) data and extensive taxon sampling including three undescribed genera of uncertain tribal placement. Our almost fully supported estimate of phylogeny confirmed the basal position of Rogadini within the subfamily and a Facitorini clade (Yeliconini+Aleiodini) that led us to propose the former group as a valid rogadine tribe (Facitorini stat. res.). Stiropiini, however, was recovered for the first time as sister to the remaining rogadine tribes except Rogadini, and Clinocentrini as sister to a clade with Betylobraconini+the three undescribed genera. The relationships recovered and morphological examination of the material included led us to place the latter three new genera and recently described genus Gondwanocentrus within a new rogadine tribe, Gondwanocentrini Shimbori & Zaldívar-Riverón trib. nov. We described these genera (Ghibli Shimbori & Zaldívar-Riverón gen. nov., Racionais Shimbori & Zaldívar-Riverón gen. nov. and Soraya Shimbori gen. nov.) with two or three new species each (G. miyazakii Shimbori & Zaldívar-Riverón sp. nov., G. totoro Shimbori & Zaldívar-Riverón sp. nov., R. brunus Shimbori & Zaldívar-Riverón sp. nov., R. kaelejay Shimbori & Zaldívar-Riverón sp. nov., R. superstes Shimbori & Zaldívar-Riverón sp. nov., S. alencarae Shimbori sp. nov. and S. venus Shimbori & Zaldívar-Riverón sp. nov.). A new species of Facitorini, Jannya pasargadae Gadelha & Shimbori sp. nov., is also described. Our newly proposed classification expands the number of tribes and genera within Rogadinae to 8 and 66 respectively. ZooBank: urn:lsid:zoobank.org:pub:51951C78-069A-4D8B-B5F0-7EBD4D9D21CE.

Daphnia sp. (Branchiopoda: Cladocera) Mitochondrial Genome Gene Rearrangement and Phylogenetic Position Within Branchiopoda

In high-altitude (4500 m) freshwater lakes, Daphnia is the apex species and the dominant zooplankton. It frequently dwells in the same lake as the Gammarid. Branchiopoda, a class of Arthropoda, Crustacea, is a relatively primitive group in the subphylum Crustacea, which originated in the Cambrian period of the Paleozoic. The complete mitogenome sequence of Daphnia sp. (Branchiopoda: Cladocera) was sequenced and annotated in this study and deposited in GenBank. The sequence structure of this species was studied by comparing the original sequences with BLAST. In addition, we have also researched the mechanisms of their mitochondrial gene rearrangement by establishing a model. We have used the Bayesian inference [BI] and maximum likelihood [ML] methods to proceed with phylogenetic analysis inference, which generates identical phylogenetic topology that reveals the phylogenetic state of Daphnia. The complete mitogenome of Daphnia sp. shows that it was 15,254 bp in length and included two control regions (CRs) and 37 genes (13 protein-coding genes, 22 tRNAs and two ribosomal RNAs [16S and 12S]). In addition to tRNA-Ser (GCT), other tRNAs have a typical cloverleaf secondary structure. Meanwhile, the mitogenome of Daphnia sp. was clearly rearranged when compared to the mitogenome of typical Daphnia. In a word, we report a newly sequenced mitogenome of Daphnia sp. with a unique rearrangement phenomenon. These results will be helpful for further phylogenetic research and provide a foundation for future studies on the characteristics of the mitochondrial gene arrangement process in Daphnia.

Phylogenomics of the geometrid tribe Palyadini (Lepidoptera: Geometridae) reveals contrasting patterns of phylogenetic signal in wing colour characters

Next generation sequencing techniques currently represent a practical and efficient way to infer robust evolutionary hypotheses. Palyadini is a small Neotropical tribe of geometrid moths composed of six genera that feature strikingly colourful wings. Here, we investigated patterns of evolution and amount of phylogenetic signal contained in various colour characters featured in the wings of members of this tribe by (i) inferring a robust phylogenetic hypothesis using ultraconserved elements (UCEs), and afterwards, (ii) mapping the morphological characters onto the molecular topology under a parsimonious ancestral character optimization. Our matrix, obtained with 60% completeness, includes 754 UCE loci and 73 taxa (64 ingroup, nine outgroup). Maximum likelihood and parsimony generated largely identical topologies with strongly supported nodes, except for one node inside the genus Opisthoxia. According to our topology, most wing colour characters are reconstructed as homoplastic, particularly at the tribe level, but five of the seven provide evidence supporting common ancestry at the genus level. Our results emphasize, once again, that no character system is infallible, and that more research is necessary to take our understanding of the evolution of wing colour in moths to a level comparable with the knowledge we have for butterflies.

Insights from the first phylogenomic analysis of flat wasps (Hymenoptera, Bethylidae) reveal two new subfamilies

Despite significant advances in alpha level taxonomy in the past few decades, the higher-level phylogeny of flat wasps (Hymenoptera, Bethylidae) remains poorly explored. Herein we provide the first phylogenomic assessment of the family based on data from ultraconserved elements for 96 species in 61 genera of the family, with material from 29 countries and all biogeographic regions. Dataset cutoffs including ultraconserved element loci recovered for 50 and 70% of terminals resulted in matrices with 1513 and 451 loci, which were analysed in both parsimony and maximum likelihood frameworks. We also provide the first analyses of divergence dating for the family based on the calibration of 12 nodes. All analyses recovered the Bethylidae as a monophyletic group and estimated the origin of the family at around 143 Ma. The results suggest that all extant subfamilies had already diversified by the Late Cretaceous. All topologies suggest that Glenosema and Chilepyris form a clade separate from other Scleroderminae; owing to the morphological distinctiveness of each genus, we propose that they are accommodated in two new subfamilies, Glenoseminae subf. nov. and Chilepyrinae subf. nov. The monotypic sclerodermine genus Galodoxa was consistently recovered within Epyrinae and is transferred to the latter subfamily.

Dietary Acculturation Is Associated With Altered Gut Microbiome, Circulating Metabolites, and Cardiovascular Disease Risk in US Hispanics and Latinos: Results From HCHS/SOL

BACKGROUND

Dietary acculturation, or adoption of dominant culture diet by migrant groups, influences human health. We aimed to examine dietary acculturation and its relationships with cardiovascular disease (CVD), gut microbiota, and blood metabolites among US Hispanic and Latino adults.

METHODS

In the HCHS/SOL (Hispanic Community Health Study/Study of Latinos), US exposure was defined by years in the United States (50 states and Washington, DC) and US nativity. A dietary acculturation pattern was derived from 14 172 participants with two 24-hour dietary recalls at baseline (2008-2011) using least absolute shrinkage and selection operator regression, with food groups as predictors of US exposure. We evaluated associations of dietary acculturation with incident CVD across ≈7 years of follow-up (n=211/14 172 cases/total) and gut microbiota (n=2349; visit 2, 2014 to 2017). Serum metabolites associated with both dietary acculturation-related gut microbiota (n=694) and incident CVD (n=108/5256 cases/total) were used as proxy measures to assess the association of diet-related gut microbiome with incident CVD.

RESULTS

We identified an empirical US-oriented dietary acculturation score that increased with US exposure. Higher dietary acculturation score was associated with higher risk of incident CVD (hazard ratio per SD, 1.33 [95% CI, 1.13-1.57]), adjusted for sociodemographic, lifestyle, and clinical factors. Sixty-nine microbial species (17 enriched from diverse species, 52 depleted mainly from fiber-utilizing Clostridia and Prevotella species) were associated with dietary acculturation, driven by lower intakes of whole grains, beans, and fruits and higher intakes of refined grains. Twenty-five metabolites, involved predominantly in fatty acid and glycerophospholipid metabolism (eg, branched-chain 14:0 dicarboxylic acid** and glycerophosphoethanolamine), were associated with both diet acculturation-related gut microbiota and incident CVD. Proxy association analysis based on these metabolites suggested a positive relationship between diet acculturation-related microbiome and risk of CVD (r=0.70, P<0.001).

CONCLUSIONS

Among US Hispanic and Latino adults, greater dietary acculturation was associated with elevated CVD risk, possibly through alterations in gut microbiota and related metabolites. Diet and microbiota-targeted interventions may offer opportunities to mitigate CVD burdens of dietary acculturation.

Zoonotic Cycle of American Trypanosomiasis in an Endemic Region of the Argentine Chaco, Factors That Influenced a Paradigm Shift

Trypanosoma cruzi, the causative agent of Chagas disease (American trypanosomiasis), is a highly complex zoonosis that is present throughout South America, Central America, and Mexico. The transmission of this disease is influenced by various factors, including human activities like deforestation and land use changes, which may have altered the natural transmission cycles and their connection to the environment. In this study conducted in the Argentine Chaco region, we examined the transmission dynamics of T. cruzi by collecting blood samples from wild and domestic animals, as well as triatomine bugs from human dwellings, across five sites of varying anthropic intervention. Samples were analyzed for T. cruzi infection via qPCR, and we additionally examined triatomines for bloodmeal analysis via NGS amplicon sequencing. Our analysis revealed a 15.3% infection rate among 20 wild species (n = 123) and no T. cruzi presence in 9 species of domestic animals (n = 1359) or collected triatomines via qPCR. Additionally, we found chicken (34.28%), human (21.59%), and goat (19.36%) as the predominant bloodmeal sources across all sites. These findings suggest that anthropic intervention and other variables analyzed may have directly impacted the spillover dynamics of T. cruzi's sylvatic cycle and potentially reduced its prevalence in human habitats.

A haplotype-resolved chromosome-level assembly and annotation of European hazelnut (C. avellana cv. Jefferson) provides insight into mechanisms of eastern filbert blight resistance

European hazelnut (Corylus avellana L.) is an important tree nut crop. Hazelnut production in North America is currently limited in scalability due to Anisogramma anomala, a fungal pathogen that causes Eastern Filbert Blight (EFB) disease in hazelnut. Successful deployment of EFB resistant cultivars has been limited to the state of Oregon, where the breeding program at Oregon State University (OSU) has released cultivars with a dominant allele at a single resistance locus identified by classical breeding, linkage mapping, and molecular markers. C. avellana cultivar "Jefferson" is resistant to the predominant EFB biotype in Oregon and has been selected by the OSU breeding program as a model for hazelnut genetic and genomic research. Here, we present a near complete, haplotype-resolved chromosome-level hazelnut genome assembly for "Jefferson". This new assembly is a significant improvement over a previously published genome draft. Analysis of genomic regions linked to EFB resistance and self-incompatibility confirmed haplotype splitting and identified new gene candidates that are essential for downstream molecular marker development, thereby facilitating breeding efforts.

Synergistic innovations enabled the radiation of anglerfishes in the deep open ocean

Major ecological transitions are thought to fuel diversification, but whether they are contingent on the evolution of certain traits called key innovations1 is unclear. Key innovations are routinely invoked to explain how lineages rapidly exploit new ecological opportunities.1,2,3 However, investigations of key innovations often focus on single traits rather than considering trait combinations that collectively produce effects of interest.4 Here, we investigate the evolution of synergistic trait interactions in anglerfishes, which include one of the most species-rich vertebrate clades in the bathypelagic, or "midnight," zone of the deep sea: Ceratioidea.5 Ceratioids are the only vertebrates that possess sexual parasitism, wherein males temporarily attach or permanently fuse to females to mate.6,7 We show that the rapid transition of ancestrally benthic anglerfishes into pelagic habitats occurred during a period of major global warming 50-35 million years ago.8,9 This transition coincided with the origins of sexual parasitism, which is thought to increase the probability of successful reproduction once a mate is found in the midnight zone, Earth's largest habitat.5,6,7 Our reconstruction of the evolutionary history of anglerfishes and the loss of immune genes support that permanently fusing clades have convergently degenerated their adaptive immunity. We find that degenerate adaptive immune genes and sexual body size dimorphism, both variably present in anglerfishes outside the ceratioid radiation, likely promoted their transition into the bathypelagic zone. These results show how traits from separate physiological, morphological, and reproductive systems can interact synergistically to drive major transitions and subsequent diversification in novel environments.

WHOLE GENOME TARGETED ENRICHMENT AND SEQUENCING OF HUMAN-INFECTING CRYPTOSPORIDIUM spp

Cryptosporidium spp. are protozoan parasites that cause severe illness in vulnerable human populations. Obtaining pure Cryptosporidium DNA from clinical and environmental samples is challenging because the oocysts shed in contaminated feces are limited in quantity, difficult to purify efficiently, may derive from multiple species, and yield limited DNA (<40 fg/oocyst). Here, we develop and validate a set of 100,000 RNA baits (CryptoCap_100k) based on six human-infecting Cryptosporidium spp. (C. cuniculus, C. hominis, C. meleagridis, C. parvum, C. tyzzeri, and C. viatorum) to enrich Cryptosporidium spp. DNA from a wide array of samples. We demonstrate that CryptoCap_100k increases the percentage of reads mapping to target Cryptosporidium references in a wide variety of scenarios, increasing the depth and breadth of genome coverage, facilitating increased accuracy of detecting and analyzing species within a given sample, while simultaneously decreasing costs, thereby opening new opportunities to understand the complex biology of these important pathogens.

Gut Microbiome Multi-Omics and Cognitive Function in the Hispanic Community Health Study/Study of Latinos- Investigation of Neurocognitive Aging

INTRODUCTION

We conducted a study within the Hispanic Community Health Study/Study of Latinos- Investigation of Neurocognitive Aging (HCHS/SOL-INCA) cohort to examine the association between gut microbiome and cognitive function.

METHODS

We analyzed the fecal metagenomes of 2,471 HCHS/SOL-INCA participants to, cross-sectionally, identify microbial taxonomic and functional features associated with global cognitive function. Omnibus (PERMANOVA) and feature-wise analyses (MaAsLin2) were conducted to identify microbiome-cognition associations, and specific microbial species and pathways (Kyoto Encyclopedia of Genes and Genomes (KEGG modules) associated with cognition.

RESULTS

Eubacterium species( E. siraeum and E. eligens ), were associated with better cognition. Several KEGG modules, most strongly Ornithine, Serine biosynthesis and Urea Cycle, were associated with worse cognition.

DISCUSSION

In a large Hispanic/Latino cohort, we identified several microbial taxa and KEGG pathways associated with cognition.

Mitogenome architecture supports the non-monophyly of the cosmopolitan parasitoid wasp subfamily Doryctinae (Hymenoptera: Braconidae) recovered by nuclear and mitochondrial phylogenomics

Mitochondrial DNA gene organisation is an important source of phylogenetic information for various metazoan taxa at different evolutionary timescales, though this has not been broadly tested for all insect groups nor within a phylogenetic context. The cosmopolitan subfamily Doryctinae is a highly diverse group of braconid wasps mainly represented by ectoparasitoids of xylophagous beetle larvae. Previous molecular studies based on Sanger and genome-wide (ultraconserved elements, UCE; and mitochondrial genomes) sequence data have recovered a non-monophyletic Doryctinae, though the relationships involved have always been weakly supported. We characterised doryctine mitogenomes and conducted separate phylogenetic analyses based on mitogenome and UCE sequence data of ~100 representative doryctine genera to assess the monophyly and higher-level classification of the subfamily. We identified rearrangements of mitochondrial transfer RNAs (tRNAs) that support a non-monophyletic Doryctinae consisting of two separate non-related clades with strong geographic structure ('New World' and 'Old World' clades). This geographic structure was also consistently supported by the phylogenetic analyses preformed with mitogenome and UCE sequence data. These results highlight the utility of the mitogenome gene rearrangements as a potential source of phylogenetic information at different evolutionary timescales.

Peripheral neuronal activation shapes the microbiome and alters gut physiology

The gastrointestinal (GI) tract is innervated by intrinsic neurons of the enteric nervous system (ENS) and extrinsic neurons of the central nervous system and peripheral ganglia. The GI tract also harbors a diverse microbiome, but interactions between the ENS and the microbiome remain poorly understood. Here, we activate choline acetyltransferase (ChAT)-expressing or tyrosine hydroxylase (TH)-expressing gut-associated neurons in mice to determine effects on intestinal microbial communities and their metabolites as well as on host physiology. The resulting multi-omics datasets support broad roles for discrete peripheral neuronal subtypes in shaping microbiome structure, including modulating bile acid profiles and fungal colonization. Physiologically, activation of either ChAT+ or TH+ neurons increases fecal output, while only ChAT+ activation results in increased colonic contractility and diarrhea-like fluid secretion. These findings suggest that specific subsets of peripherally activated neurons differentially regulate the gut microbiome and GI physiology in mice without involvement of signals from the brain.

Gut microbiome is associated with recurrence-free survival in patients with resected Stage IIIB-D or Stage IV melanoma treated with immune checkpoint inhibitors

The gut microbiome (GMB) has been associated with outcomes of immune checkpoint blockade therapy in melanoma, but there is limited consensus on the specific taxa involved, particularly across different geographic regions. We analyzed pre-treatment stool samples from 674 melanoma patients participating in a phase-III trial of adjuvant nivolumab plus ipilimumab versus nivolumab, across three continents and five regions. Longitudinal analysis revealed that GMB was largely unchanged following treatment, offering promise for lasting GMB-based interventions. In region-specific and cross-region meta-analyses, we identified pre-treatment taxonomic markers associated with recurrence, including Eubacterium, Ruminococcus, Firmicutes, and Clostridium. Recurrence prediction by these markers was best achieved across regions by matching participants on GMB compositional similarity between the intra-regional discovery and external validation sets. AUCs for prediction ranged from 0.83-0.94 (depending on the initial discovery region) for patients closely matched on GMB composition (e.g., JSD ≤0.11). This evidence indicates that taxonomic markers for prediction of recurrence are generalizable across regions, for individuals of similar GMB composition.

Does urbanisation lead to parallel demographic shifts across the world in a cosmopolitan plant?

Urbanisation is occurring globally, leading to dramatic environmental changes that are altering the ecology and evolution of species. In particular, the expansion of human infrastructure and the loss and fragmentation of natural habitats in cities is predicted to increase genetic drift and reduce gene flow by reducing the size and connectivity of populations. Alternatively, the 'urban facilitation model' suggests that some species will have greater gene flow into and within cities leading to higher diversity and lower differentiation in urban populations. These alternative hypotheses have not been contrasted across multiple cities. Here, we used the genomic data from the GLobal Urban Evolution project (GLUE), to study the effects of urbanisation on non-adaptive evolutionary processes of white clover (Trifolium repens) at a global scale. We found that white clover populations presented high genetic diversity and no evidence of reduced Ne linked to urbanisation. On the contrary, we found that urban populations were less likely to experience a recent decrease in effective population size than rural ones. In addition, we found little genetic structure among populations both globally and between urban and rural populations, which showed extensive gene flow between habitats. Interestingly, white clover displayed overall higher gene flow within urban areas than within rural habitats. Our study provides the largest comprehensive test of the demographic effects of urbanisation. Our results contrast with the common perception that heavily altered and fragmented urban environments will reduce the effective population size and genetic diversity of populations and contribute to their isolation.

Under the hood: Phylogenomics of hooded tick spiders (Arachnida, Ricinulei) uncovers discordance between morphology and molecules

Ricinulei or hooded tick-spiders are a cryptic and ancient group of arachnids. The order consists of around 100 highly endemic extant species restricted to the Afrotropics and the Neotropics along with 22 fossil species. Their antiquity and low vagility make them an excellent group with which to interrogate biogeographic questions. To date, only four molecular analyses have been conducted on the group and they failed to resolve the relationships of the main lineages and even recovering the non-monophyly of the three genera. These studies were limited to a few Sanger loci or phylogenomic analyses with at most seven ingroup samples. To increase phylogenetic resolution in this little-understood and poorly studied group, we present the most comprehensive phylogenomic study of Ricinulei to date leveraging the Arachnida ultra-conserved element probe set. With a data set of 473 loci across 96 ingroup samples, analyses resolved a monophyletic Neotropical clade consisting of four main lineages. Two of them correspond to the current genera Cryptocellus and Pseudocellus while topology testing revealed one lineage to likely be a phylogenetic reconstruction artefact. The fourth lineage, restricted to Northwestern, Andean South America, is consistent with the Cryptocellus magnus group, likely corresponding to the historical genus Heteroricinoides. Since we did not sample the type species for this old genus, we do not formally re-erect Heteroricinoides but our data suggest the need for a thorough morphological re-examination of Neotropical Ricinulei.

Integration of polygenic and gut metagenomic risk prediction for common diseases

Multiomics has shown promise in noninvasive risk profiling and early detection of various common diseases. In the present study, in a prospective population-based cohort with ~18 years of e-health record follow-up, we investigated the incremental and combined value of genomic and gut metagenomic risk assessment compared with conventional risk factors for predicting incident coronary artery disease (CAD), type 2 diabetes (T2D), Alzheimer disease and prostate cancer. We found that polygenic risk scores (PRSs) improved prediction over conventional risk factors for all diseases. Gut microbiome scores improved predictive capacity over baseline age for CAD, T2D and prostate cancer. Integrated risk models of PRSs, gut microbiome scores and conventional risk factors achieved the highest predictive performance for all diseases studied compared with models based on conventional risk factors alone. The present study demonstrates that integrated PRSs and gut metagenomic risk models improve the predictive value over conventional risk factors for common chronic diseases.

Lessons from assembling UCEs: A comparison of common methods and the case of Clavinomia (Halictidae)

Sequence data assembly is a foundational step in high-throughput sequencing, with untold consequences for downstream analyses. Despite this, few studies have interrogated the many methods for assembling phylogenomic UCE data for their comparative efficacy, or for how outputs may be impacted. We study this by comparing the most commonly used assembly methods for UCEs in the under-studied bee lineage Nomiinae and a representative sampling of relatives. Data for 63 UCE-only and 75 mixed taxa were assembled with five methods, including ABySS, HybPiper, SPAdes, Trinity and Velvet, and then benchmarked for their relative performance in terms of locus capture parameters and phylogenetic reconstruction. Unexpectedly, Trinity and Velvet trailed the other methods in terms of locus capture and DNA matrix density, whereas SPAdes performed favourably in most assessed metrics. In comparison with SPAdes, the guided-assembly approach HybPiper generally recovered the highest quality loci but in lower numbers. Based on our results, we formally move Clavinomia to Dieunomiini and render Epinomia once more a subgenus of Dieunomia. We strongly advise that future studies more closely examine the influence of assembly approach on their results, or, minimally, use better-performing assembly methods such as SPAdes or HybPiper. In this way, we can move forward with phylogenomic studies in a more standardized, comparable manner.

WHOLE GENOME TARGETED ENRICHMENT AND SEQUENCING OF HUMAN-INFECTING CRYPTOSPORIDIUM spp

Cryptosporidium spp. are protozoan parasites that cause severe illness in vulnerable human populations. Obtaining pure Cryptosporidium DNA from clinical and environmental samples is challenging because the oocysts shed in contaminated feces are limited in quantity, difficult to purify efficiently, may derive from multiple species, and yield limited DNA (<40 fg/oocyst). Here, we develop and validate a set of 100,000 RNA baits (CryptoCap_100k) based on six human-infecting Cryptosporidium spp. ( C. cuniculus , C. hominis , C. meleagridis , C. parvum , C. tyzzeri , and C. viatorum ) to enrich Cryptosporidium spp. DNA from a wide array of samples. We demonstrate that CryptoCap_100k increases the percentage of reads mapping to target Cryptosporidium references in a wide variety of scenarios, increasing the depth and breadth of genome coverage, facilitating increased accuracy of detecting and analyzing species within a given sample, while simultaneously decreasing costs, thereby opening new opportunities to understand the complex biology of these important pathogens.

An island 'endemic' born out of hybridization between introduced lineages

Humans have profoundly impacted the distribution of plant and animal species over thousands of years. The most direct example of these effects is human-mediated movement of individuals, either through translocation of individuals within their range or through the introduction of species to new habitats. While human involvement may be suspected in species with obvious range disjunctions, it can be difficult to detect natural versus human-mediated dispersal events for populations at the edge of a species' range, and this uncertainty muddles how we understand the evolutionary history of populations and broad biogeographical patterns. Studies combining genetic data with archaeological, linguistic and historical evidence have confirmed prehistoric examples of human-mediated dispersal; however, it is unclear whether these methods can disentangle recent dispersal events, such as species translocated by European colonizers during the past 500 years. We use genomic DNA from historical museum specimens and historical records to evaluate three hypotheses regarding the timing and origin of Northern Bobwhites (Colinus virginianus) in Cuba, whose status as an endemic or introduced population has long been debated. We discovered that bobwhites from southern Mexico arrived in Cuba between the 12th and 16th centuries, followed by the subsequent introduction of bobwhites from the southeastern USA to Cuba between the 18th and 20th centuries. These dates suggest the introduction of bobwhites to Cuba was human-mediated and concomitant with Spanish colonial shipping routes between Veracruz, Mexico and Havana, Cuba during this period. Our results identify endemic Cuban bobwhites as a genetically distinct population born of hybridization between divergent, introduced lineages.

A conserved interdomain microbial network underpins cadaver decomposition despite environmental variables

Microbial breakdown of organic matter is one of the most important processes on Earth, yet the controls of decomposition are poorly understood. Here we track 36 terrestrial human cadavers in three locations and show that a phylogenetically distinct, interdomain microbial network assembles during decomposition despite selection effects of location, climate and season. We generated a metagenome-assembled genome library from cadaver-associated soils and integrated it with metabolomics data to identify links between taxonomy and function. This universal network of microbial decomposers is characterized by cross-feeding to metabolize labile decomposition products. The key bacterial and fungal decomposers are rare across non-decomposition environments and appear unique to the breakdown of terrestrial decaying flesh, including humans, swine, mice and cattle, with insects as likely important vectors for dispersal. The observed lockstep of microbial interactions further underlies a robust microbial forensic tool with the potential to aid predictions of the time since death.

Microgeographic population structuring in a genus of California trapdoor spiders and discovery of an enigmatic new species (Euctenizidae: Promyrmekiaphila korematsui sp. nov.)

The recognition and delineation of cryptic species remains a perplexing problem in systematics, evolution, and species delimitation. Once recognized as such, cryptic species complexes provide fertile ground for studying genetic divergence within the context of phenotypic and ecological divergence (or lack thereof). Herein we document the discovery of a new cryptic species of trapdoor spider, Promyrmekiaphila korematsui sp. nov. Using subgenomic data obtained via target enrichment, we document the phylogeography of the California endemic genus Promyrmekiaphila and its constituent species, which also includes P. clathrata and P. winnemem. Based on these data we show a pattern of strong geographic structuring among populations but cannot entirely discount recent gene flow among populations that are parapatric, particularly for deeply diverged lineages within P. clathrata. The genetic data, in addition to revealing a new undescribed species, also allude to a pattern of potential phenotypic differentiation where species likely come into close contact. Alternatively, phenotypic cohesion among genetically divergent P. clathrata lineages suggests that some level of gene flow is ongoing or occurred in the recent past. Despite considerable field collection efforts over many years, additional sampling in potential zones of contact for both species and lineages is needed to completely resolve the dynamics of divergence in Promyrmekiaphila at the population-species interface.

Quantitative and qualitative mutational impact of ionizing radiation on normal cells

The comprehensive genomic impact of ionizing radiation (IR), a carcinogen, on healthy somatic cells remains unclear. Using large-scale whole-genome sequencing (WGS) of clones expanded from irradiated murine and human single cells, we revealed that IR induces a characteristic spectrum of short insertions or deletions (indels) and structural variations (SVs), including balanced inversions, translocations, composite SVs (deletion-insertion, deletion-inversion, and deletion-translocation composites), and complex genomic rearrangements (CGRs), including chromoplexy, chromothripsis, and SV by breakage-fusion-bridge cycles. Our findings suggest that 1 Gy IR exposure causes an average of 2.33 mutational events per Gb genome, comprising 2.15 indels, 0.17 SVs, and 0.01 CGRs, despite a high level of inter-cellular stochasticity. The mutational burden was dependent on total irradiation dose, regardless of dose rate or cell type. The findings were further validated in IR-induced secondary cancers and single cells without clonalization. Overall, our study highlights a comprehensive and clear picture of IR effects on normal mammalian genomes.

A target enrichment approach for enhanced recovery of Synchytrium endobioticum nuclear genome sequences

Potato wart disease is caused by the obligate fungal pathogen Synchytrium endobioticum. DNA extraction from compost, purified spores and crude wart tissue derived from tuber galls of infected potatoes often results in low S. endobioticum DNA concentration or highly contaminated with DNA coming from other microorganisms and the potato host. Therefore, Illumina sequencing of these samples generally results in suboptimal recovery of the nuclear genome sequences of S. endobioticum. A hybridization-based target enrichment protocol was developed to strongly enhance the recovery of S. endobioticum DNA while off-target organisms DNA remains uncaptured. The design strategy involved creating a set of 180,000 molecular baits targeting both gene and non-gene regions of S. endobioticum. The baits were applied to whole genome amplified DNA samples of various S. endobioticum pathotypes (races) in compost, from purified spores and crude wart tissue samples. This was followed by Illumina sequencing and bioinformatic analyses. Compared to non-enriched samples, target enriched samples: 1) showed a significant increase in the proportion of sequenced bases mapped to the S. endobioticum nuclear genome, especially for crude wart tissue samples; 2) yielded sequencing data with higher and better nuclear genome coverage; 3) biased genome assembly towards S. endobioticum sequences, yielding smaller assembly sizes but higher representation of putative S. endobioticum contigs; 4) showed an increase in the number of S. endobioticum genes detected in the genome assemblies. Our hybridization-based target enrichment protocol offers a valuable tool for enhancing genome sequencing and NGS-based molecular detection of S. endobioticum, especially in difficult samples.

Long-term gut colonization with ESBL-producing Escherichia coli in participants without known risk factors from the southeastern United States

We evaluated gut carriage of extended spectrum beta lactamase producing Enterobacteriaceae (ESBL-E) in southeastern U.S. residents without recent in-patient healthcare exposure. Study enrollment was January 2021-February 2022 in Athens, Georgia, U.S. and included a diverse population of 505 adults plus 50 child participants (age 0-5). Based on culture-based screening of stool samples, 4.5% of 555 participants carried ESBL-Es. This is slightly higher than reported in studies conducted 2012-2015, which found carriage rates of 2.5-3.9% in healthy U.S. residents. All ESBL-E confirmed isolates (n=25) were identified as Escherichia coli. Isolates belonged to 11 sequence types, with 48% classified as ST131. Ninety six percent of ESBL-E isolates carried a blaCTX-M gene. Isolated ESBL-Es frequently carried virulence genes as well as multiple classes of antibiotic resistance genes. Long-term colonization was common, with 64% of ESBL-E positive participants testing positive when rescreened three months later. One participant yielded isolates belonging to two different E. coli sequence types that carried blaCTX-M-1 genes on near-identical plasmids, suggesting intra-gut plasmid transfer. Isolation of E. coli on media without antibiotics revealed that ESBL-E. coli typically made up a minor fraction of the overall gut E. coli population, although in some cases they were the dominant strain. ESBL-E carriage was not associated with a significantly different stool microbiome composition. However, some microbial taxa were differentially abundant in ESBL-E carriers. Together, these results suggest that a small subpopulation of US residents are long-term, asymptomatic carriers of ESBL-Es, and may serve as an important reservoir for community spread of these ESBL genes.

Young evolutionary origins of dioecy in the genus Asparagus

PREMISE

Dioecy (separate sexes) has independently evolved numerous times across the angiosperm phylogeny and is recently derived in many lineages. However, our understanding is limited regarding the evolutionary mechanisms that drive the origins of dioecy in plants. The recent and repeated evolution of dioecy across angiosperms offers an opportunity to make strong inferences about the ecological, developmental, and molecular factors influencing the evolution of dioecy, and thus sex chromosomes. The genus Asparagus (Asparagaceae) is an emerging model taxon for studying dioecy and sex chromosome evolution, yet estimates for the age and origin of dioecy in the genus are lacking.

METHODS

We use plastome sequences and fossil time calibrations in phylogenetic analyses to investigate the age and origin of dioecy in the genus Asparagus. We also review the diversity of sexual systems present across the genus to address contradicting reports in the literature.

RESULTS

We estimate that dioecy evolved once or twice approximately 2.78-3.78 million years ago in Asparagus, of which roughly 27% of the species are dioecious and the remaining are hermaphroditic with monoclinous flowers.

CONCLUSIONS

Our findings support previous work implicating a young age and the possibility of two origins of dioecy in Asparagus, which appear to be associated with rapid radiations and range expansion out of Africa. Lastly, we speculate that paleoclimatic oscillations throughout northern Africa may have helped set the stage for the origin(s) of dioecy in Asparagus approximately 2.78-3.78 million years ago.

Role of Gut Microbiota in Statin-Associated New-Onset Diabetes-A Cross-Sectional and Prospective Analysis of the FINRISK 2002 Cohort

BACKGROUND

Dyslipidemia is treated effectively with statins, but treatment has the potential to induce new-onset type-2 diabetes. Gut microbiota may contribute to this outcome variability. We assessed the associations of gut microbiota diversity and composition with statins. Bacterial associations with statin-associated new-onset type-2 diabetes (T2D) risk were also prospectively evaluated.

METHODS

We examined shallow-shotgun-sequenced fecal samples from 5755 individuals in the FINRISK-2002 population cohort with a 17+-year-long register-based follow-up. Alpha-diversity was quantified using Shannon index and beta-diversity with Aitchison distance. Species-specific differential abundances were analyzed using general multivariate regression. Prospective associations were assessed with Cox regression. Applicable results were validated using gradient boosting.

RESULTS

Statin use associated with differing taxonomic composition (R2, 0.02%; q=0.02) and 13 differentially abundant species in fully adjusted models (MaAsLin; q<0.05). The strongest positive association was with Clostridium sartagoforme (β=0.37; SE=0.13; q=0.02) and the strongest negative association with Bacteroides cellulosilyticus (β=-0.31; SE=0.11; q=0.02). Twenty-five microbial features had significant associations with incident T2D in statin users, of which only Bacteroides vulgatus (HR, 1.286 [1.136-1.457]; q=0.03) was consistent regardless of model adjustment. Finally, higher statin-associated T2D risk was seen with [Ruminococcus] torques (ΔHRstatins, +0.11; q=0.03), Blautia obeum (ΔHRstatins, +0.06; q=0.01), Blautia sp. KLE 1732 (ΔHRstatins, +0.05; q=0.01), and beta-diversity principal component 1 (ΔHRstatin, +0.07; q=0.03) but only when adjusting for demographic covariates.

CONCLUSIONS

Statin users have compositionally differing microbiotas from nonusers. The human gut microbiota is associated with incident T2D risk in statin users and possibly has additive effects on statin-associated new-onset T2D risk.

Sex Hormones, the Stool Microbiome, and Subclinical Atherosclerosis in Women With and Without HIV

CONTEXT

Cardioprotective roles of endogenous estrogens may be particularly important in women with HIV, who have reduced estrogen exposure and elevated cardiovascular disease risk. The gut microbiome metabolically interacts with sex hormones, but little is known regarding possible impact on cardiovascular risk.

OBJECTIVE

To analyze potential interplay of sex hormones and gut microbiome in cardiovascular risk.

METHODS

Among 197 postmenopausal women in the Women's Interagency HIV Study, we measured 15 sex hormones in serum and assessed the gut microbiome in stool. Presence of carotid artery plaque was determined (B-mode ultrasound) in a subset (n = 134). We examined associations of (i) sex hormones and stool microbiome, (ii) sex hormones and plaque, and (iii) sex hormone-related stool microbiota and plaque, adjusting for potential confounders.

RESULTS

Participant median age was 58 years and the majority were living with HIV (81%). Sex hormones (estrogens, androgens, and adrenal precursors) were associated with stool microbiome diversity and specific species, similarly in women with and without HIV. Estrogens were associated with higher diversity, higher abundance of species from Alistipes, Collinsella, Erysipelotrichia, and Clostridia, and higher abundance of microbial β-glucuronidase and aryl-sulfatase orthologs, which are involved in hormone metabolism. Several hormones were associated with lower odds of carotid artery plaque, including dihydrotestosterone, 3α-diol-17G, estradiol, and estrone. Exploratory mediation analysis suggested that estrone-related species, particularly from Collinsella, may mediate the protective association of estrone with plaque.

CONCLUSION

Serum sex hormones are significant predictors of stool microbiome diversity and composition. The gut microbiome may play a role in estrogen-related cardiovascular protection.

A comprehensive dataset on spatiotemporal variation of microbial plankton communities in the Baltic Sea

The Baltic Sea is one of the largest brackish water environments on earth and is characterised by pronounced physicochemical gradients and seasonal dynamics. Although the Baltic Sea has a long history of microscopy-based plankton monitoring, DNA-based metabarcoding has so far mainly been limited to individual transect cruises or time-series of single stations. Here we report a dataset covering spatiotemporal variation in prokaryotic and eukaryotic microbial communities and physicochemical parameters. Within 13-months between January 2019 and February 2020, 341 water samples were collected at 22 stations during monthly cruises along the salinity gradient. Both salinity and seasonality are strongly reflected in the data. Since the dataset was generated with both metabarcoding and microscopy-based methods, it provides unique opportunities for both technical and ecological analyses, and is a valuable biodiversity reference for future studies, in the prospect of climate change.

DNA metabarcoding reveals spatial and temporal variation of fish eye fluke communities in lake ecosystems

Eye flukes (Diplostomidae) are diverse and abundant trematode parasites that form multi-species communities in fish with negative effects on host fitness and survival. However, the environmental factors and host-related characteristics that determine species diversity, composition, and coexistence in such communities remain poorly understood. Here, we developed a cost-effective cox1 region-specific DNA metabarcoding approach to characterize parasitic diplostomid communities in two common fish species (Eurasian perch and common roach) collected from seven temperate lakes in Estonia. We found considerable inter- and intra-lake, as well as inter-host species, variation in diplostomid communities. Sympatric host species characterization revealed that parasite communities were typically more diverse in roach than perch. Additionally, we detected five positive and two negative diplostomid species associations in roach, whereas only a single negative association was observed in perch. These results indicate that diplostomid communities in temperate lakes are complex and dynamic systems exhibiting both spatial and temporal heterogeneity. They are influenced by various environmental factors and by host-parasite and inter-parasite interactions. We expect that the described methodology facilitates ecological and biodiversity research of diplostomid parasites. It is also adaptable to other parasite groups where it could serve to improve current understanding of diversity, distribution, and interspecies interactions of other understudied taxa.

Detection of skunk adenovirus-1 in domestic ferrets (Mustela putorius)

We describe gross pathology, histopathology, ultrastructure, and molecular characterization of skunk adenovirus 1 (SkAdV-1) in 3 ferrets (Mustela putorius) submitted to the Infectious Diseases Laboratory, Zoo and Exotic Animal Pathology Service at the University of Georgia. Case 1 came from a cohort with a 6-week history of multiple ferrets with respiratory disease and mortality. Cases 2 and 3 were submitted 19 days later; they died 3 days after arrival at the vendor's facility. Histopathological evaluation found evidence of suppurative bronchopneumonia, with intralesional gram-positive cocci in case 1. The others had evidence of subacute multifocal histiocytic and lymphoplasmacytic interstitial pneumonia. Viral isolation, polymerase chain reaction (PCR), and DNA sequencing identified SkAdV-1 in lung tissue. In situ hybridization confirmed adenovirus within the lung lesion, and adenovirus particles were visible under an electron microscope. This is the seventh species in which SkAdV-1 has been identified. Cross-species transmission events have important implications for wildlife management and multispecies facilities.

The menopause-related gut microbiome: associations with metabolomics, inflammatory protein markers, and cardiometabolic health in women with HIV

OBJECTIVE

This study aimed to identify menopause-related gut microbial features, as well as their related metabolites and inflammatory protein markers, and link with cardiometabolic risk factors in women with and without HIV.

METHODS

In the Women's Interagency HIV Study, we performed shotgun metagenomic sequencing on 696 stool samples from 446 participants (67% women with HIV), and quantified plasma metabolomics and serum proteomics in a subset (~86%). We examined the associations of menopause (postmenopausal vs premenopausal) with gut microbial features in a cross-sectional repeated-measures design and further evaluated those features in relation to metabolites, proteins, and cardiometabolic risk factors.

RESULTS

Different overall gut microbial composition was observed by menopausal status in women with HIV only. We identified a range of gut microbial features that differed between postmenopausal and premenopausal women with HIV (but none in women without HIV), including abundance of 32 species and functional potentials involving 24 enzymatic reactions and lower β-glucuronidase bacterial gene ortholog. Specifically, highly abundant species Faecalibacterium prausnitzii , Bacteroides species CAG:98 , and Bifidobacterium adolescentis were depleted in postmenopausal versus premenopausal women with HIV. Menopause-depleted species (mainly Clostridia ) in women with HIV were positively associated with several glycerophospholipids, while negatively associated with imidazolepropionic acid and fibroblast growth factor 21. Mediation analysis suggested that menopause may decrease plasma phosphatidylcholine plasmalogen C36:1 and C36:2 levels via reducing abundance of species F. prausnitzii and Acetanaerobacterium elongatum in women with HIV. Furthermore, waist-to-hip ratio was associated with menopause-related microbes, metabolites, and fibroblast growth factor 21 in women with HIV.

CONCLUSIONS

Menopause was associated with a differential gut microbiome in women with HIV, related to metabolite and protein profiles that potentially contribute to elevated cardiometabolic risk.

Primate-specific ZNF808 is essential for pancreatic development in humans

Identifying genes linked to extreme phenotypes in humans has the potential to highlight biological processes not shared with all other mammals. Here, we report the identification of homozygous loss-of-function variants in the primate-specific gene ZNF808 as a cause of pancreatic agenesis. ZNF808 is a member of the KRAB zinc finger protein family, a large and rapidly evolving group of epigenetic silencers which target transposable elements. We show that loss of ZNF808 in vitro results in aberrant activation of regulatory potential contained in the primate-specific transposable elements it represses during early pancreas development. This leads to inappropriate specification of cell fate with induction of genes associated with liver identity. Our results highlight the essential role of ZNF808 in pancreatic development in humans and the contribution of primate-specific regions of the human genome to congenital developmental disease.

Genetic variants associated with hantavirus infection in a reservoir host are related to regulation of inflammation and immune surveillance

The immunogenetics of wildlife populations influence the epidemiology and evolutionary dynamic of the host-pathogen system. Profiling immune gene diversity present in wildlife may be especially important for those species that, while not at risk of disease or extinction themselves, are host to diseases that are a threat to humans, other wildlife, or livestock. Hantaviruses (genus: Orthohantavirus) are globally distributed zoonotic RNA viruses with pathogenic strains carried by a diverse group of rodent hosts. The marsh rice rat (Oryzomys palustris) is the reservoir host of Orthohantavirus bayoui, a hantavirus that causes fatal cases of hantavirus cardiopulmonary syndrome in humans. We performed a genome wide association study (GWAS) using the rice rat "immunome" (i.e., all exons related to the immune response) to identify genetic variants associated with infection status in wild-caught rice rats naturally infected with their endemic strain of hantavirus. First, we created an annotated reference genome using 10× Chromium Linked Reads sequencing technology. This reference genome was used to create custom baits which were then used to target enrich prepared rice rat libraries (n = 128) and isolate their immunomes prior to sequencing. Top SNPs in the association test were present in four genes (Socs5, Eprs, Mrc1, and Il1f8) which have not been previously implicated in hantavirus infections. However, these genes correspond with other loci or pathways with established importance in hantavirus susceptibility or infection tolerance in reservoir hosts: the JAK/STAT, MHC, and NFκB. These results serve as informative markers for future exploration and highlight the importance of immune pathways that repeatedly emerge across hantavirus systems. Our work aids in creating cross-species comparisons for better understanding mechanisms of genetic susceptibility and host-pathogen coevolution in hantavirus systems.

Widespread sympatry in a species-rich clade of marine fishes (Carangoidei)

A universal paradigm describing patterns of speciation across the tree of life has been debated for decades. In marine organisms, inferring patterns of speciation using contemporary and historical patterns of biogeography is challenging due to the deficiency of species-level phylogenies and information on species' distributions, as well as conflicting relationships between species' dispersal, range size and co-occurrence. Most research on global patterns of marine fish speciation and biogeography has focused on coral reef or pelagic species. Carangoidei is an ecologically important clade of marine fishes that use coral reef and pelagic environments. We used sequence capture of 1314 ultraconserved elements (UCEs) from 154 taxa to generate a time-calibrated phylogeny of Carangoidei and its parent clade, Carangiformes. Age-range correlation analyses of the geographical distributions and divergence times of sister species pairs reveal widespread sympatry, with 73% of sister species pairs exhibiting sympatric geographical distributions, regardless of node age. Most species pairs coexist across large portions of their ranges. We also observe greater disparity in body length and maximum depth between sympatric relative to allopatric sister species. These and other ecological or behavioural attributes probably facilitate sympatry among the most closely related carangoids.

Lack of host specialization despite selective host use in brood parasitic cuckoo catfish

Host-parasite dynamics involve coevolutionary arms races, which may lead to host specialization and ensuing diversification. Our general understanding of the evolution of host specialization in brood parasites is compromised by a restricted focus on bird and insect lineages. The cuckoo catfish (Synodontis multipunctatus) is an obligate parasite of parental care of mouthbrooding cichlids in Lake Tanganyika. Given the ecological and taxonomic diversity of mouthbrooding cichlids in the lake, we hypothesized the existence of sympatric host-specific lineages in the cuckoo catfish. In a sample of 779 broods from 20 cichlid species, we found four species parasitized by cuckoo catfish (with prevalence of parasitism of 2%-18%). All parasitized cichlids were from the tribe Tropheini, maternal mouthbrooders that spawn over a substrate (rather than in open water). Phylogenetic analysis based on genomic (ddRAD sequencing) and mitochondrial (Dloop) data from cuckoo catfish embryos showed an absence of host-specific lineages. This was corroborated by analyses of genetic structure and co-ancestry matrix. Within host species, parasitism was not associated with any individual characteristic we recorded (parent size, water depth), but was costly as parasitized parents carried smaller clutches of their own offspring. We conclude that the cuckoo catfish is an intermediate generalist and discuss costs, benefits and constraints of host specialization in this species and brood parasites in general.

Bamboozled! Resolving deep evolutionary nodes within the phylogeny of bamboo corals (Octocorallia: Scleralcyonacea: Keratoisididae)

Keratoisididae is a globally distributed, and exclusively deep-sea, family of octocorals that contains species and genera that are polyphyletic. An alphanumeric system, based on a three-gene-region phylogeny, is widely used to describe the biodiversity within this family. That phylogeny identified 12 major groups although it did not have enough signal to explore the relationships among groups. Using increased phylogenomic resolution generated from Ultraconserved Elements and exons (i.e. conserved elements), we aim to resolve deeper nodes within the family and investigate the relationships among those predefined groups. In total, 109 libraries of conserved elements were generated from individuals representing both the genetic and morphological diversity of our keratoisidids. In addition, the conserved element data of 12 individuals from previous studies were included. Our taxon sampling included 11 of the 12 keratoisidid groups. We present two phylogenies, constructed from a 75% (231 loci) and 50% (1729 loci) taxon occupancy matrix respectively, using both Maximum Likelihood and Multiple Species Coalescence methods. These trees were congruent at deep nodes. As expected, S1 keratoisidids were recovered as a well-supported sister clade to the rest of the bamboo corals. S1 corals do not share the same mitochondrial gene arrangement found in other members of Keratoisididae. All other bamboo corals were recovered within two major clades. Clade I comprises individuals assigned to alphanumeric groups B1, C1, D1&D2, F1, H1, I4, and J3 while Clade II contains representatives from A1, I1, and M1. By combining genomics with already published morphological data, we provide evidence that group H1 is not monophyletic, and that the division between other groups - D1 and D2, and A1 and M1 - needs to be reconsidered. Overall, there is a lack of robust morphological markers within Keratoisididae, but subtle characters such as sclerite microstructure and ornamentation seem to be shared within groups and warrant further investigation as taxonomically diagnostic characters.

Rampant chloroplast capture in Sarracenia revealed by plastome phylogeny

Introgression can produce novel genetic variation in organisms that hybridize. Sympatric species pairs in the carnivorous plant genus Sarracenia L. frequently hybridize, and all known hybrids are fertile. Despite being a desirable system for studying the evolutionary consequences of hybridization, the extent to which introgression occurs in the genus is limited to a few species in only two field sites. Previous phylogenomic analysis of Sarracenia estimated a highly resolved species tree from 199 nuclear genes, but revealed a plastid genome that is highly discordant with the species tree. Such cytonuclear discordance could be caused by chloroplast introgression (i.e. chloroplast capture) or incomplete lineage sorting (ILS). To better understand the extent to which introgression is occurring in Sarracenia, the chloroplast capture and ILS hypotheses were formally evaluated. Plastomes were assembled de-novo from sequencing reads generated from 17 individuals in addition to reads obtained from the previous study. Assemblies of 14 whole plastomes were generated and annotated, and the remaining fragmented assemblies were scaffolded to these whole-plastome assemblies. Coding sequence from 79 homologous genes were aligned and concatenated for maximum-likelihood phylogeny estimation. The plastome tree is extremely discordant with the published species tree. Plastome trees were simulated under the coalescent and tree distance from the species tree was calculated to generate a null distribution of discordance that is expected under ILS alone. A t-test rejected the null hypothesis that ILS could cause the level of discordance seen in the plastome tree, suggesting that chloroplast capture must be invoked to explain the discordance. Due to the extreme level of discordance in the plastome tree, it is likely that chloroplast capture has been common in the evolutionary history of Sarracenia.

Genotyped functional screening of soluble Fab clones enables in-depth analysis of mutation effects

Monoclonal antibodies (mAbs) and their fragments are widely used in therapeutics, diagnostics and basic research. Although display methods such as phage display offer high-throughput, affinities of individual antibodies need to be accurately measured in soluble format. We have developed a screening platform capable of providing genotyped functional data from a total of 9216 soluble, individual antigen binding fragment (Fab) clones by employing next-generation sequencing (NGS) with hierarchical indexing. Full-length, paired variable domain sequences (VL-VH) are linked to functional screening data, enabling in-depth analysis of mutation effects. The platform was applied to four phage display-selected scFv/Fab screening projects and one site-saturation VH affinity maturation project. Genotyped functional screening simultaneously enabled the identification of affinity improving mutations in the VH domain of Fab 49A3 recognizing Dengue virus non-structural protein 1 (NS1) serotype 2 and informed on VH residue positions which cannot be changed from wild-type without decreasing the affinity. Genotype-based identification revealed to us the extent of intraclonal signal variance inherent to single point screening data, a phenomenon often overlooked in the field. Moreover, genotyped screening eliminated the redundant selection of identical genotypes for further study and provided a new analysis tool to evaluate the success of phage display selections and remaining clonal diversity in the screened repertoires.

An elevational shift facilitated the Mesoamerican diversification of Azure-hooded Jays (Cyanolyca cucullata) during the Great American Biotic Interchange

The Great American Biotic Interchange (GABI) was a key biogeographic event in the history of the Americas. The rising of the Panamanian land bridge ended the isolation of South America and ushered in a period of dispersal, mass extinction, and new community assemblages, which sparked competition, adaptation, and speciation. Diversification across many bird groups, and the elevational zonation of others, ties back to events triggered by the GABI. But the exact timing of these events is still being revealed, with recent studies suggesting a much earlier time window for faunal exchange, perhaps as early as 20 million years ago (Mya). Using a time-calibrated phylogenetic tree, we show that the jay genus Cyanolyca is emblematic of bird dispersal trends, with an early, pre-land bridge dispersal from Mesoamerica to South America 6.3-7.3 Mya, followed by a back-colonization of C. cucullata to Mesoamerica 2.3-4.8 Mya, likely after the land bridge was complete. As Cyanolyca species came into contact in Mesoamerica, they avoided competition due to a prior shift to lower elevation in the ancestor of C. cucullata. This shift allowed C. cucullata to integrate itself into the Mesoamerican highland avifauna, which our time-calibrated phylogeny suggests was already populated by higher-elevation, congeneric dwarf-jays (C. argentigula, C. pumilo, C. mirabilis, and C. nanus). The outcome of these events and fortuitous elevational zonation was that C. cucullata could continue colonizing new highland areas farther north during the Pleistocene. Resultingly, four C. cucullata lineages became isolated in allopatric, highland regions from Panama to Mexico, diverging in genetics, morphology, plumage, and vocalizations. At least two of these lineages are best described as species (C. mitrata and C. cucullata). Continued study will further document the influence of the GABI and help clarify how dispersal and vicariance shaped modern-day species assemblages in the Americas.

Library preparation method and DNA source influence endogenous DNA recovery from 100-year-old avian museum specimens

Museum specimens collected prior to cryogenic tissue storage are increasingly being used as genetic resources, and though high-throughput sequencing is becoming more cost-efficient, whole genome sequencing (WGS) of historical DNA (hDNA) remains inefficient and costly due to its short fragment sizes and high loads of exogenous DNA, among other factors. It is also unclear how sequencing efficiency is influenced by DNA sources. We aimed to identify the most efficient method and DNA source for collecting WGS data from avian museum specimens. We analyzed low-coverage WGS from 60 DNA libraries prepared from four American Robin (Turdus migratorius) and four Abyssinian Thrush (Turdus abyssinicus) specimens collected in the 1920s. We compared DNA source (toepad versus incision-line skin clip) and three library preparation methods: (1) double-stranded DNA (dsDNA), single tube (KAPA); (2) single-stranded DNA (ssDNA), multi-tube (IDT); and (3) ssDNA, single tube (Claret Bioscience). We found that the ssDNA, multi-tube method resulted in significantly greater endogenous DNA content, average read length, and sequencing efficiency than the other tested methods. We also tested whether a predigestion step reduced exogenous DNA in libraries from one specimen per species and found promising results that warrant further study. The ~10% increase in average sequencing efficiency of the best-performing method over a commonly implemented dsDNA library preparation method has the potential to significantly increase WGS coverage of hDNA from bird specimens. Future work should evaluate the threshold for specimen age at which these results hold and how the combination of library preparation method and DNA source influence WGS in other taxa.

Phylogenomic analyses using a new 1013-gene Vitaceae bait-set support major groups of North American Vitis

A set of newly designed Vitaceae baits targeting 1013 genes was employed to explore phylogenetic relationships among North American Vitis. Eurasian Vitis taxa including Vitis vinifera were found to be nested within North American Vitis subgenus Vitis. North American Vitis subgenus Vitis can be placed into nine main groups: the Monticola group, the Occidentales group, the Californica group, the Vinifera group (introduced from Eurasia), the Mustangensis group, the Palmata group, the Aestivalis group, the Labrusca group, and the Cinerea group. Strong cytonuclear discordances were detected in North American Vitis, with many species non-monophyletic in the plastid phylogeny, while monophyletic in the nuclear phylogeny. The phylogenomic analyses support recognizing four distinct species in the Vitis cinerea complex in North America: V. cinerea, V. baileyana, V. berlandieri, and V. simpsonii. Such treatment will better serve the conservation of wild Vitis diversity in North America. Yet the evolutionary history of Vitis is highly complex, with the concordance analyses indicating conflicting signals across the phylogeny. Cytonuclear discordances and Analyses using the Species Networks applying Quartets (SNaQ) method support extensive hybridizations in North American Vitis. The results further indicate that plastid genomes alone are insufficient for resolving the evolutionary history of plant groups that have undergone rampant hybridization, like the case in North American Vitis. Nuclear gene data are essential for species delimitation, identification and reconstructing evolutionary relationships; therefore, they are imperative for plant phylogenomic studies.

Phylogeography of the blacklegged tick (Ixodes scapularis) throughout the USA identifies candidate loci for differences in vectorial capacity

The blacklegged tick (Ixodes scapularis (Journal of the Academy of Natural Sciences of Philadelphia, 1821, 2, 59)) is a vector of Borrelia burgdorferi sensu stricto (s.s.) (International Journal of Systematic Bacteriology, 1984, 34, 496), the causative bacterial agent of Lyme disease, part of a slow-moving epidemic of Lyme borreliosis spreading across the northern hemisphere. Well-known geographical differences in the vectorial capacity of these ticks are associated with genetic variation. Despite the need for detailed genetic information in this disease system, previous phylogeographical studies of these ticks have been restricted to relatively few populations or few genetic loci. Here we present the most comprehensive phylogeographical study of genome-wide markers in I. scapularis, conducted by using 3RAD (triple-enzyme restriction-site associated sequencing) and surveying 353 ticks from 33 counties throughout the species' range. We found limited genetic variation among populations from the Northeast and Upper Midwest, where Lyme disease is most common, and higher genetic variation among populations from the South. We identify five spatially associated genetic clusters of I. scapularis. In regions where Lyme disease is increasing in frequency, the I. scapularis populations genetically group with ticks from historically highly Lyme-endemic regions. Finally, we identify 10 variable DNA sites that contribute the most to population differentiation. These variable sites cluster on one of the chromosome-scale scaffolds for I. scapularis and are within identified genes. Our findings illuminate the need for additional research to identify loci causing variation in the vectorial capacity of I. scapularis and where additional tick sampling would be most valuable to further understand disease trends caused by pathogens transmitted by I. scapularis.

Range-wide population genomics of common seadragons shows secondary contact over a former barrier and insights on illegal capture

BACKGROUND

Common seadragons (Phyllopteryx taeniolatus, Syngnathidae) are an emblem of the diverse endemic fauna of Australia's southern rocky reefs, the newly recognized "Great Southern Reef." A lack of assessments spanning this global biodiversity hotspot in its entirety is currently hampering an understanding of the factors that have contributed to its diversity. The common seadragon has a wide range across Australia's entire temperate south and includes a geogenetic break over a former land bridge, which has called its status as a single species into question. As a popular aquarium display that sells for high prices, common seadragons are also vulnerable to illegal capture.

RESULTS

Here, we provide range-wide nuclear sequences (986 variable Ultraconserved Elements) for 198 individuals and mitochondrial genomes for 140 individuals to assess species status, identify genetic units and their diversity, and trace the source of two poached individuals. Using published data of the other two seadragon species, we found that lineages of common seadragons have diverged relatively recently (< 0.63 Ma). Within common seadragons, we found pronounced genetic structure, falling into three major groups in the western, central, and eastern parts of the range. While populations across the Bassian Isthmus were divergent, there is also evidence for secondary contact since the passage opened. We found a strong cline of genetic diversity from the range center tapering symmetrically towards the range peripheries. Based on their genetic similarities, the poached individuals were inferred to have originated from around Albany in southwestern Australia.

CONCLUSIONS

We conclude that common seadragons constitute a single species with strong geographic structure but coherence through gene flow. The low genetic diversity on the east and west coasts is concerning given that these areas are projected to face fast climate change. Our results suggest that in addition to their life history, geological events and demographic expansions have all played a role in shaping populations in the temperate south. These insights are an important step towards understanding the historical determinants of the diversity of species endemic to the Great Southern Reef.

Gut microbiota, circulating inflammatory markers and metabolites, and carotid artery atherosclerosis in HIV infection

BACKGROUND

Alterations in gut microbiota have been implicated in HIV infection and cardiovascular disease. However, how gut microbial alterations relate to host inflammation and metabolite profiles, and their relationships with atherosclerosis, have not been well-studied, especially in the context of HIV infection. Here, we examined associations of gut microbial species and functional components measured by shotgun metagenomics with carotid artery plaque assessed by B-mode carotid artery ultrasound in 320 women with or at high risk of HIV (65% HIV +) from the Women's Interagency HIV Study. We further integrated plaque-associated microbial features with serum proteomics (74 inflammatory markers measured by the proximity extension assay) and plasma metabolomics (378 metabolites measured by liquid chromatography tandem mass spectrometry) in relation to carotid artery plaque in up to 433 women.

RESULTS

Fusobacterium nucleatum, a potentially pathogenic bacteria, was positively associated with carotid artery plaque, while five microbial species (Roseburia hominis, Roseburia inulinivorans, Johnsonella ignava, Odoribacter splanchnicus, Clostridium saccharolyticum) were inversely associated with plaque. Results were consistent between women with and without HIV. Fusobacterium nucleatum was positively associated with several serum proteomic inflammatory markers (e.g., CXCL9), and the other plaque-related species were inversely associated with proteomic inflammatory markers (e.g., CX3CL1). These microbial-associated proteomic inflammatory markers were also positively associated with plaque. Associations between bacterial species (especially Fusobacterium nucleatum) and plaque were attenuated after further adjustment for proteomic inflammatory markers. Plaque-associated species were correlated with several plasma metabolites, including the microbial metabolite imidazole-propionate (ImP), which was positively associated with plaque and several pro-inflammatory markers. Further analysis identified additional bacterial species and bacterial hutH gene (encoding enzyme histidine ammonia-lyase in ImP production) associated with plasma ImP levels. A gut microbiota score based on these ImP-associated species was positively associated with plaque and several pro-inflammatory markers.

CONCLUSION

Among women living with or at risk of HIV, we identified several gut bacterial species and a microbial metabolite ImP associated with carotid artery atherosclerosis, which might be related to host immune activation and inflammation. Video Abstract.