Sequence diversity studies of rat brain RNA: effects of environmental complexity on rat brain RNA diversity.
Explore the source record for details and available documents.
SEARCH · PubMed Health
Explore indexed PubMed citations for clinical trials, systematic reviews and public health research. Read source abstracts and follow each citation to its original PubMed record.
Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.
Explore the source record for details and available documents.
The most frequently overlooked aspect of conservation of natural resources and economic development is human cultural diversity. However, conservation and development of natural resources are basically human-oriented endeavors and all conservation and developmental efforts ought to start with a clear understanding of the varied needs of the people. In addition, cultural diversity is a natural resource that ought to be protected along with all the more commonly recognized resources of the ecosystems of which humans form an integral part. Cultural diversity in a large measure is an ecological phenomenon because such diversity includes variations in the ways different people perceive and utilize the environments in which they live. Thus, cultural diversity, in large measure, equates with ecological diversity. It has been well established that a high degree of ecological diversity (including taxonomic, niche, biogeochemical and other measures of diversity) is a necessary attribute of humid tropical ecosystems if such ecosystems are to remain viable over long periods of time. The current land-use trends in Central America and Panama are leading toward ever larger areas being devoted to monocultural use with sharply reduced ecological diversity that poses great dangers for the near and long terms. In addition, the resultant removal of people from rural areas results in growing social, economic and political problems, that are not being successfully met by developing nations. Not only is there a growing wastage of human resources, a weakening of the social structures, and an increasing and dangerous dependence upon monocultures oriented toward export markets, but the ecological diversity of the previous existing land-use systems are being lost. This kind of diversity is probably no less valuable to the ecological health of a nation's agriculture, forestry, and general resource utilization than is the genetic diversity of "primitive" crop plant varieties which biologists now recognize and increasingly seek to preserve as "modern" crop plant varieties become ever more simple genetically and hence ever more vulnerable to disease and other perturbations.
This study provides the first metagenomic assessment of microbial diversity from the tea rhizosphere of the Kangra valley. Tea rhizosphere soil samples were collected from 4 locations (Dharamshala, Baijnath, Palampur, and Joginder Nagar) of the Kangra valley. DNA extracts of rhizosphere samples were analysed for bacterial and Archaeal diversity using amplicon sequencing (V3-V4) region of the 16S rRNA gene and Fungal diversity using ITS1 and ITS2 regions. Baijnath and Palampur samples showed the highest bacterial richness, while Dharamshala and Palampur had the highest fungal richness. Proteobacteria was a dominant phylum in all the rhizosphere samples, followed by Firmicutes, Actinobacteria, Acidobacteria, and Bacteroidetes. A total of 11 fungal phyla were identified among all the locations, with abundance of Ascomycota and Basidiomycota. For the Archaea domain, uncultured archaeon and Aeropyrum camini were the most common found among all the locations. A small fraction (< 0.5%) of Bacillus and Pseudomonas species were observed among all the locations. Alpha and beta diversity indices displayed notable differences within and between microbial diversities. Soil factors were variably associated with microbial diversity, with nitrogen positively aligned with fungal diversity, while EC and K were associated with Archaeal diversity. Soil pH and OM% showed moderate associations with bacterial diversity. These findings provided valuable and comprehensive insights into tea rhizosphere microbial ecology and could be used to better understand microbial functions and their role in plant health.
Gout, a painful inflammatory arthritis, is characterized by hyperuricemia and monosodium urate crystal deposition, with growing evidence linking its pathogenesis to gut microbiome dysbiosis. However, traditional diversity metrics fail to capture the complex spatial organization of microbial communities. This study addresses this gap by applying the novel metagenomic Diversity-Area Relationship (m-DAR) model to investigate scaling laws in the gout microbiome-quantifying how metagenomic diversity changes with the number of individuals sampled. Our analysis of gut microbiomes from gout patients and healthy controls revealed fundamental ecological disruptions. We found that gout microbiomes exhibited significantly altered scaling patterns: they showed greater inter-individual dissimilarity (higher z-values) at the level of rare genes (q = 0), but weaker scaling of dominant genes (q = 1-3) compared to healthy controls. Crucially, the maximal accrual diversity (MAD) was substantially lower in gout patients, indicating a severely constrained potential for total microbial gene diversity. Furthermore, profiling of metagenomic functional gene clusters (MFGCs) uncovered widespread functional perturbations, including increased diversity scaling for carbohydrate-active enzymes (CAZy) but decreased scaling in essential metabolic pathways (KEGG, KO). These results demonstrate that the gout gut microbiome is defined by a loss of ecological structure, featuring reduced homogeneity in dominant taxa, expanded rare biosphere variation, and an overall collapsed diversity capacity. This work introduces an ecological framework for characterizing dysbiosis in gout that complements traditional diversity metrics and may inform the development of microbiome-based therapeutic strategies. Further research is needed to translate these ecological patterns into clinical applications.
Phylogenetic diversity (PD) plays an important role in biodiversity, conservation, and evolutionary studies by measuring the diversity of a set of taxa based on their phylogenetic relationships. In phylogenetic trees, a subset of k taxa with maximum PD can be found by a simple and efficient greedy algorithm. However, this algorithmic tractability is lost when considering phylogenetic networks, which incorporate reticulate evolutionary events such as hybridization and horizontal gene transfer. To address this challenge, we introduce PaNDA (Phylogenetic Network Diversity Algorithms), the first software package and interactive graphical user-interface for exploring, visualizing, and maximizing diversity in phylogenetic networks. PaNDA includes a novel algorithm to find a subset of k taxa with maximum diversity, running in polynomial time for networks of bounded scanwidth, a measure of tree-likeness of a network that grows slower than the well-known level measure. This algorithm considers the variant of PD on networks in which the branch lengths of all paths from the root to the selected taxa contribute towards their diversity. We demonstrate the scalability of this algorithm on simulated networks, successfully analyzing level-15 networks with up to 200 taxa in seconds. We also provide a proof-of-concept analysis using a phylogenetic network on Xiphophorus species, illustrating how the tool can support diversity studies based on real genomic data. The software is easily installable and freely available at https://github.com/nholtgrefe/panda. Additionally, we extend the definition of PD to semi-directed phylogenetic networks, which are mixed graphs increasingly used in phylogenetic analysis to model uncertainty of the root location. We prove that finding a subset of k taxa with maximum diversity remains NP-hard on semi-directed networks, but do present a polynomial-time algorithm for networks with bounded level.
Biodiversity has experienced tremendous shifts in community, species, and genetic diversity during the Anthropocene. Understanding temporal diversity shifts is especially critical in biodiversity hotspots, i.e., regions that are exceptionally biodiverse and threatened. Here, we use museomics and temporal genomics approaches to quantify temporal shifts in genomic diversity in an assemblage of eight generalist highland bird species from the Ethiopian Highlands (part of the Eastern Afromontane Biodiversity Hotspot). With genomic data from contemporary and historical samples, we demonstrate an assemblage-wide trend of increased genomic diversity through time, potentially due to improved habitat connectivity within highland regions. Genomic diversity shifts in these generalist species contrast with general trends of genomic diversity declines in specialist or imperiled species. In addition to genetic diversity shifts, we found an assemblage-wide trend of decreased realized mutational load, indicative of overall trends for potentially deleterious variation to be masked or selectively purged. Across this avian assemblage, we also show that shifts in population genomic structure are idiosyncratic, with species-specific trends. These results are in contrast with other charismatic and imperiled African taxa that have largely shown strong increases in population genetic structure over the recent past. This study highlights that not all taxa respond the same to environmental change, and generalists, in some cases, may even respond positively. Future comparative conservation genomics assessments on species groups or assemblages with varied natural history characteristics would help us better understand how diverse taxa respond to anthropogenic landscape changes.
UNLABELLED: Vaccines targeting Neisseria gonorrhoeae are needed to reduce disease burden and help address the problem of antimicrobial resistance, with an understanding of relationships between gonococcal genetics and molecules influencing diversity, infection, and the immune response essential for developing effective vaccine formulations. Whole-genome sequence data can be used to investigate these relationships among thousands of gonococcal isolates, allowing the study of antigenic diversity on a population scale. Such analyses typically examine antigenic diversity occurring in complete protein sequences, generating mean diversity indices and phylogenetic analyses that can inform on vaccine potential; however, to detect and measure the immune responses elicited, epitope characterization within an antigen helps guide vaccine formulations, with epitopes commonly located in surface-exposed regions of a protein. Here, we analyzed the genetic diversity of the major gonococcal antigen, PorB, in WGS from 22,227 N. gonorrhoeae isolates. We characterized the diversity of all eight surface-exposed outer membrane loops, or variable regions (VRs), and generated a PorB VR subtyping scheme to facilitate the global and temporal detection of circulating PorB subtypes. These analyses identified the presence of dominant VR combinations that persisted over time, indicative of (i) epistatic interactions between VRs and (ii) positive selection. Strain-specific, anti-PorB IgG responses directed toward distinct VR subtypes were detected in sera obtained from participants vaccinated with 4CMenB. The deconstruction of PorB into each surface-exposed loop provides a powerful approach for evaluating vaccine candidates: the methods used here allow immunodominant regions to be detected, which is invaluable for further vaccine investigations. IMPORTANCE: In the context of rising global gonorrhea cases, the development of vaccines becomes a priority; however, N. gonorrhoeae antigenic diversity and its ability to evade the immune system complicate vaccine development. This study characterizes the genetic diversity of the outer membrane protein, PorB, a key component of the outer membrane and a major gonococcal antigen. Using genomics and machine-learning techniques, this research identified dominant PorB variants that drive the immune response, proposing potential vaccine candidates and improving our understanding of the evolutionary forces maintaining genome structure and biological fitness. Understanding these processes is crucial for designing vaccines that effectively target N. gonorrhoeae and combat the spread of multidrug-resistant gonococci.
Indigenous chickens play a critical role in food security and climate resilience in smallholder systems, yet their genomic diversity and adaptive potential remain insufficiently characterised. This study employed low-pass whole-genome sequencing (LP-WGS; 0.2-1.99×) to investigate genomic diversity, population structure, inbreeding and candidate environment-associated genomic variation in 33 chickens from highland, midland, and lowland agroecologies in the Tigray region of northern Ethiopia. After imputation and stringent filtering, 23.4 million high-confidence SNPs were retained, including ~ 17% novel variants, indicating substantial uncharacterised genetic diversity in these populations. SNP density (13.8 ± 8.6 SNPs/kb) was comparable to values reported from high-coverage Ethiopian chicken datasets, demonstrating the suitability of LP-WGS for population genomics in resource-limited settings. Marked differences in genomic diversity were observed among ecotypes: midland chickens showed the highest nucleotide diversity (π = 0.00267), followed by lowland (π = 0.00233), whereas highland chickens showed the lowest diversity (π = 0.00203) and elevated genomic inbreeding (FROH and FHOM ≈ 0.18). Population structure analyses revealed clear genetic separation among ecotypes. PCA (13.91% variation explained) distinguished lowland chickens along PC1 and separated highland from midland along PC2, while ADMIXTURE and FST patterns supported three major ancestral genomic backgrounds. Functional annotation of private missense variants uncovered distinct adaptive signatures reflecting the contrasting agroecological conditions. Highland chickens showed enrichment of candidate genes potentially involved in physiological processes relevant to high-altitude environments, including cold response, angiogenesis, cardiovascular regulation and metabolic homeostasis (eg., PARP1, ACOX2, ITGB3, EDNRB, SOX8, and SOX10). Midland chickens exhibited candidate signals of selection in genes with known roles in innate antiviral immunity, bacterial defence and inflammatory regulation (eg., BAK1, CLSTN1, CYSLTR1, CYSLTR2, CXCR7, GIPR, DSCAM, GDAP1, TLR3, TLR4, TLR7, IFIH1, ADORA1, EPHB1, and TMPRSS2). Lowland chickens displayed candidate variants associated with heat-stress response, DNA damage repair, oxidative balance and cardiovascular support under extreme temperatures (e.g., MLH1, BDKRB1, GPR19, FLT1, CCL18, TGM2, and RAMP3). Overall, the results indicate substantial genomic differentiation among ecotypes and suggest candidate environment-associated genetic divergence across Tigray's diverse agroecological zones. These populations may represent important reservoirs of adaptive genetic variation for climate-resilient poultry breeding, warranting further functional validation and conservation-oriented management.
The gut microbiota is increasingly recognized as a target for obesity management; however, whether baseline gut microbial diversity conditions responsiveness to microbiota-targeted interventions remains unclear. We aimed to investigate whether baseline gut microbial diversity is associated with responsiveness to a paraprobiotic derived from Lactiplantibacillus plantarum LRCC5282 (LP5282-P) in overweight adults. In a 12-week, randomized, double-blind, placebo-controlled, multicenter trial of 120 overweight adults, LP5282-P produced no significant between-group differences in any clinical outcome across the overall per-protocol population. However, in the low-diversity subgroup, LP5282-P was associated with significant reductions in body weight, body mass index, and circulating leptin levels. These clinical changes were accompanied by compositional shifts in the gut microbiota, including higher relative abundances of Christensenellaceae, Faecalibacterium, and Alistipes. Fecal metabolite profiles showed elevated acetate and butyrate concentrations and altered bile acid composition. Within the low-diversity subgroup, changes in the relative abundances of Akkermansia and Eubacterium were inversely correlated with changes in body weight, body fat mass, and leptin levels. In contrast, the high-diversity subgroup exhibited no consistent response across the outcome domains examined. Overall, baseline gut microbial diversity was associated with differential responsiveness to LP5282-P, supporting its potential use as a stratification variable in future microbiota-targeted intervention trials. Further studies integrating direct measures of microbial activity and host response are warranted to elucidate the biological pathways underlying this diversity-dependent responsiveness. Trial registration: Clinical Research Information Service (CRIS), KCT0008119.
Genetic diversity within species is the basis for evolutionary adaptive capacity and has recently been included as a target for protection in the United Nations' Global Biodiversity Framework (GBF). However, we lack large-scale mathematical frameworks to quantify how much genetic diversity has already been lost, let alone to predict future losses under 21st century conservation scenarios. To fill this gap, we developed an area-based spatio-temporal predictive framework of genetic diversity calibrated with population-scale genomic data of 29 plant and animal species. To estimate present genetic diversity loss with our framework, we used species' habitat area and population sizes losses reported in the Living Planet Index, the Red List, and new GBF indicators across 13,808 species for the last 5 decades. Applying our evolutionary framework across these species, we estimate genetic diversity loss lags behind population and habitat area declines, with an estimated current 13-22% π genetic diversity loss. However, we forecast future genetic diversity losses will reach 41-76% even if populations are not further contracted. These results highlight that safeguarding existing habitats is insufficient to maintain the genetic health of species and relying solely on continuous genetic monitoring underestimates lagging long term impacts.
Drosophila melanogaster Down Syndrome cell adhesion molecule 1 (Dscam1) gene encodes 38,016 diverse cell surface receptor proteins via alternative splicing, which have both nervous and immune functions. However, it remains elusive why organisms have evolved such an astonishing diversity of isoforms. Here, we show that fitness and immunity properties have driven the modern evolution of Dscam1 isoform diversity. We assess multiple aspects of fly fitness in deletion mutants harboring exon 4, 6, or 9 clusters, respectively, reducing ectodomain isoform diversity stepwise from 18,612 to 396. All fitness-related traits generally improved as the potential number of isoforms increased; however, the magnitude of the changes varied remarkably in a variable cluster-specific manner. Correlation analysis revealed that fitness-related traits were much more sensitive to reductions in Dscam1 diversity compared to canonical neuronal self/non-self discrimination. We conclude that the role of Dscam1 isoforms in canonical neuronal self-avoidance and self/non-self discrimination is mediated by a small fraction of all isoforms (<1/10), whereas a separate role essential for other developmental contexts and resistances, likely in fitness and immunity, requires almost full isoform diversity. Thus, fitness and immunity properties, rather than canonical neuronal functions, are the dominant drivers during the modern diversification of the Dscam1 isoform. Our findings suggest that Dscam1 diversity is closely linked to adaptation and species diversification in arthropods.
Since biliary excretion of metabolites might determine the pattern of intestinal neoplasms induced by azoxymethane, the number and distribution of tumors were studied in rats after pancreaticobiliary diversion to the mid small bowel. Pancreaticobiliary diversion was performed either immediately before the first of 16 weekly injections of azoxymethane or 10 days after the last. Seven months after pancreaticobiliary diversion, persistent ileal hyperplasia was manifested by higher levels of mucosal RNA and DNA compared with controls (34--102%: P less than 0.001), while there was little residual adaptation in the colon. Qualitative and quantitative analysis of fecal bile acids 6--26 wk after pancreaticobiliary diversion showed few differences. Pancreaticobiliary diversion doubled the incidence of colonic tumors, whether operation preceded (P less than 0.005) or followed (P less than 0.02) the course of azoxymethane. Suture-line tumors were common in the small bowel, particularly in the transposed duodenal stump. Despite intense ileal hyperplasia as a consequence of pancreaticobiliary diversion, the ileum remained resistant to chemical carcinogenesis. The potentiation of colonic neoplasms by pancreaticobiliary diversion probably depends on the stimulation of colonic mucosal proliferation.
Azadirachta indica (A. indica) A. Juss., commonly known as Neem, is a valuable multipurpose tree with profound medicinal properties and socioeconomic importance, widely recognized since ancient Ayurvedic times. Despite its prominence, knowledge about its genetic diversity within the metropolitan area of Ahmedabad is limited. This study marks the first in-depth exploration of the genetic diversity and population structure of A. indica in Ahmedabad. The authenticity of the species was validated through DNA barcoding, and a Geographical Information System (GIS) was used to collect the samples. A total of 35 A. indica accessions were analyzed using five Inter Simple Sequence Repeat (ISSR) primers. Genetic diversity and population structure were evaluated using Inter Simple Sequence Repeat (ISSR) markers through polymorphism assessment, clustering, ordination, and Bayesian population structure analyses. ISSRs revealed a high level of polymorphism (75.66%), indicating substantial genetic variability among accessions. An analysis of genetic diversity indices revealed low to moderate diversity (Hs = 0.14, Ht = 0.217, I = 0.217). Analysis of Molecular Variance (AMOVA) analysis depicted 81% variation within the population and 19% among the population. Low to moderate genetic differentiation (Gst = 0.319) and moderate gene flow (Nm = 1.06) indicated that urban development has not hindered gene flow among populations. Mantel's test revealed a weak but significant correlation between genetic and geographic distances, suggesting limited isolation by distance. The estimated ΔK using STRUCTURE exhibited two subpopulations, representing two gene pools for A. indica accessions (K = 2). Collectively, these patterns indicate that urbanization has not severely disrupted genetic connectivity in A. indica, reflecting its resilience and adaptive potential in a metropolitan environment. These findings provide pivotal knowledge for further understanding the genetic diversity and population structure of A. indica in one of the fastest-growing cities in India, which can be utilized for new breeding programmes, sustainable development and future conservation strategies around the globe.
Genetic diversity loss due to anthropogenic factors is occurring rapidly on a global scale, putting many species at risk of extirpation and extinction. Different management strategies have been developed to slow this loss; however, it is often unknown whether these strategies reach their intended goals. In this study, we evaluate population structure and changes in nucleotide diversity (π) in isolated populations of Bonneville cutthroat trout (Oncorhynchus clarkii utah) from the Snake Range (Nevada, USA). Starting in the 1990s, three of these populations were used to reestablish populations in the Snake Range because many of the historic populations were extirpated. Some populations were stocked using a single-source and others were stocked using multiple-sources. Using low-coverage whole-genome sequencing coupled with historic samples (2003-2010) and contemporary samples (2019-2022), we find that single-source populations lost nucleotide diversity while mixed-source populations maintained nucleotide diversity. Further, source populations used to restore populations throughout the Snake Range lost the most nucleotide diversity over the time span evaluated. Our findings provide insight into how small, isolated populations can be managed to maintain genetic diversity.
BACKGROUND: Long-term exposure to indoor fungal bioaerosols is a recognised risk factor for respiratory illness, particularly in damp and poorly ventilated housing. However, the diversity and seasonal variability of these fungal communities are poorly understood. As part of the West London Healthy Home and Environment Study (WellHome), this study aimed to characterise the composition, diversity, and temporal dynamics of indoor fungal bioaerosols in urban UK homes, as compared with outdoor air, to inform future exposure baselines and policy development. METHODS: In this prospective, community-based observational study, 118 households were recruited across West London, UK, via community networks and partner organisations, prioritising families with children aged 5-17 years with asthma or allergies, from diverse socioeconomic backgrounds. Sampling occurred between Oct 3, 2022, and June 14, 2024. Participant data were collected via questionnaires completed by household members, capturing demographics, building characteristics, and respiratory health. Passive-air samplers were used in living rooms for 28 days during two seasonal campaigns, with concurrent outdoor sampling at four fixed community sites. Fungal bioaerosols were identified by ITS2 amplicon sequencing and quantified using broad-range quantitative PCR targeting the 18S rRNA gene. Diversity indexes and temporal dynamics were analysed using ecological statistics and generalised additive models. FINDINGS: 118 households were enrolled, comprising 504 residents (263 women, 237 men, and four not reported). Among 504 participants who self-identified, the largest groups comprised individuals identifying as Black African (n=47), Somali (n=46), White British (n=42), and African (n=38), with additional representation from mixed race ethnic backgrounds (n=29), Black British (n=27), White (n=22), and Black Caribbean (n=18), alongside several other ethnicities each represented at lower frequencies. Of 118 households, 104 completed both seasonal campaigns and 14 completed one, yielding 262 air samples (222 indoor and 40 outdoor). DNA was successfully recovered from all samples, identifying 2027 fungal genera. Indoor environments showed significantly higher richness (mean 646 vs 495 amplicon sequence variants; p<0·0001) and Shannon diversity (4·21 vs 3·53; p<0·0001) than outdoors. Community composition differed markedly (permutational multivariate ANOVA p<0·0001), with Penicillium, Aspergillus, and Wallemia enriched indoors. Indoor fungal communities presented stronger seasonal cycling (R2=0·203) than outdoor communities (R2=0·012). Fungal burden across all homes had a median 11 043 genomic equivalence (GE); IQR 4598-20 579 GE. The highest levels were observed in homes with visible mould; one household showed elevated Aspergillus exposure linked to repeated asthma hospitalisations in a sensitised resident. INTERPRETATION: Indoor fungal bioaerosols are more diverse and dynamic than outdoor communities in urban UK homes. These findings establish foundational exposure data and highlight the need for incorporating fungal bioaerosol monitoring into public health policy to mitigate mould-related health risks. FUNDING: UK Research and Innovation (UKRI) Strategic Priorities Fund (SPF) Clean Air Programme.
Cyanobacterial blooms are an immense global problem that release complex mixtures of poorly characterized biologically active cyanopeptides into freshwater. In this study, metabolomics and genomics were used to assess the diversity and concentrations of cyanopeptides during a dense Microcystis bloom during the late summer of 2023 in Lake Champlain, a large transboundary lake situated between Canada and the United States. Despite the relatively low genetic diversity of the bloom determined by 16S rRNA metabarcoding, 151 cyanopeptides were detected by non-targeted metabolomics. This represents the most recorded cyanopeptides from a single lake plankton bloom event to date. Fifty-two cyanopeptides were previously reported and 99 represent putative new structures. Standards from the microcystin, cyanopeptolin, microginin, and anabaenopeptin groups were used to either quantify or approximate respective cyanopeptide concentrations over the sampling period. Cyanopeptolins were the most diverse (n = 68) cyanopeptides and the second most abundant, reaching 12,892 μg/L. Microginins were the second most diverse (n = 24) and reached the highest concentrations (18,262 μg/L). Anabaenopeptins were the third most diverse (n = 17) cyanopeptides, reaching 4,818 μg/L. Only 8 microcystins were detected, reaching 4,935 μg/L, where MC-LR was the dominant congener. Target cyanopeptide biosynthesis genes for microcystins (mcyE), cyanopeptolins (mcnC), anabaenopeptins (apnD), microviridins (mdnC), and aeruginosins (aerA) were also quantified using digital droplet PCR (ddPCR). The gene copy numbers for mcyE, mcnC, and apnD were highly correlated with their corresponding cyanopeptide concentrations. Overall, the studied Microcystis bloom produced a very diverse cyanopeptide mixture with high cyanopeptide concentrations including non-microcystin groups.
Enterococci appear to have originated in the guts of early terrestrializing arthropods and invertebrates over 425 million years ago-hosts that are now highly diverse and widespread in nature today. Yet most knowledge of the genus comes from human infection-associated lineages with genomes swollen by the recent accretion of foreign DNA conveyed by mobile elements. Because invertebrates dominate terrestrial animal diversity and biomass, they would be predicted to constitute a major but little-explored reservoir of enterococcal diversity. We therefore systematically examined Enterococcus association and species diversification in invertebrate hosts of the comparatively natural, isolated, but well-characterized environment of the Azorean island of Terceira. Over 100 invertebrate specimens were examined for associated enterococci, which were taxonomically classified by whole-genome sequencing. Supporting the existence of a large pool of uncharacterized enterococci and Enterococcus-adapted genes, 40% (eight of 20) of the Enterococcus species identified were either undescribed, including four candidate new species described here, or very recently discovered. In contrast, control isolates from vertebrates were exclusively of known species typical of sampling elsewhere, discounting geographic isolation as a main driver of the novelty observed. Further, because of the abundance of E. casseliflavus and E. flavescens in this collection, we obtained the resolution necessary to quantify the divergence and decipher the drivers of speciation in the controversial division between these naturally vancomycin-resistant species. These findings provide robust support for the existence of a large pool of new species and unexplored adaptive traits in invertebrate-associated enterococci-diverse environmental survival traits optimized for expression in an enterococcal background, and well positioned for transmission into human-associated enterococcal strains.IMPORTANCEEnterococci are auxotrophic gut-associated bacteria that co-evolved with their terrestrial hosts over many eons. In the last 75 years-the "antibiotic era"-E. faecalis and E. faecium gained genes for antibiotic resistance and enhanced virulence, emerging as leading causes of multidrug-resistant infection. Little is known about the source of those genes or the pathway by which they entered human-associated strains. A recent global survey suggested a potentially large repository of uncharacterized genetic diversity in the enterococci of invertebrates. We directly tested this prospect by examining enterococci of invertebrate hosts in a largely natural and pastoral environment. Our findings provide clear evidence that invertebrates naturally harbor vast unexplored enterococcal diversity. Moreover, associations are likely driven by intrinsic host selection factors rather than geographic isolation. This expands our knowledge of Enterococcus biodiversity, including the identification of four novel species, identifying a vast reservoir of enterococcal genes available to species that colonize and infect humans.
Retrons are bacterial genetic elements encoding a specialized reverse transcriptase (RT) that synthesizes multicopy single-stranded DNA and are increasingly recognized as components of bacterial anti-phage defense systems. However, their diversity and ecological distribution across large-scale genomic resources remain poorly characterized. Here, we surveyed retron RTs across the SPIRE representative metagenome collection, a non-redundant, species-level data set spanning diverse microbial habitats. Using a curated panel of type-specific hidden Markov models, we identified retrons representing all canonical types together with additional divergent lineages. Retron distribution showed strong taxonomic and ecological structuring, with some groups restricted to specific bacterial phyla, whereas others were broadly distributed across environmental categories. Systematic novelty assessment identified two candidate type XI-like lineages, TXI_C2like and TXI_noncan_h, characterized by protease-independent architectures and distinct accessory modules associated with WYL- and DnaB_C-containing proteins, respectively. De novo covariance-based analyses further identified candidate msr/msd-like non-coding RNA structures in both lineages, supporting conservation of the canonical RT-ncRNA organizational framework despite extensive sequence divergence. Together, these findings expand the known diversity of retron systems and identify type XI-like retrons as a dynamic and previously underexplored evolutionary group.IMPORTANCERetrons are bacterial genetic elements that are increasingly exploited as programmable tools for genome editing, molecular recording, and biosensing in addition to their natural role in anti-phage defense. Despite this growing biotechnological interest, the true diversity of retrons across the bacterial world has remained largely unmapped. By mining a resource of over 100,000 processed microbial metagenomes, we uncovered thousands of retron sequences spanning known types as well as previously unrecognized lineages and found that their distribution is strongly shaped by both bacterial taxonomy and ecological niche. Among these, we identified two candidate new lineages related to type XI retrons that lack the protease domain typical of this group but instead carry distinct accessory proteins, expanding the known architectural diversity of these systems. These findings broaden the catalog of retron diversity available for functional characterization and biotechnological engineering and provide a framework for prioritizing candidate lineages for future experimental validation.