- Current Issue
Volume 10, Issue 3, 2024
- Research Articles
-
- Functional Genomics and Microbe–Niche Interactions
-
-
Phylometagenomics of cycad coralloid roots reveals shared symbiotic signals
Cycads are known to host symbiotic cyanobacteria, including Nostocales species, as well as other sympatric bacterial taxa within their specialized coralloid roots. Yet, it is unknown if these bacteria share a phylogenetic origin and/or common genomic functions that allow them to engage in facultative symbiosis with cycad roots. To address this, we obtained metagenomic sequences from 39 coralloid roots sampled from diverse cycad species and origins in Australia and Mexico. Culture-independent shotgun metagenomic sequencing was used to validate sub-community co-cultures as an efficient approach for functional and taxonomic analysis. Our metanalysis shows a host-independent microbiome core consisting of seven bacterial orders with high species diversity within the identified taxa. Moreover, we recovered 43 cyanobacterial metagenome-assembled genomes, and in addition to Nostoc spp., symbiotic cyanobacteria of the genus Aulosira were identified for the first time. Using this robust dataset, we used phylometagenomic analysis to reveal three monophyletic cyanobiont clades, two host-generalist and one cycad-specific that includes Aulosira spp. Although the symbiotic clades have independently arisen, they are enriched in certain functional genes, such as those related to secondary metabolism. Furthermore, the taxonomic composition of associated sympatric bacterial taxa remained constant. Our research quadruples the number of cycad cyanobiont genomes and provides a robust framework to decipher cyanobacterial symbioses, with the potential of improving our understanding of symbiotic communities. This study lays a solid foundation to harness cyanobionts for agriculture and bioprospection, and assist in conservation of critically endangered cycads.
-
- Pathogens and Epidemiology
-
-
Genomic analysis of clinical Aeromonas isolates reveals genetic diversity but little evidence of genetic determinants for diarrhoeal disease
Aeromonas spp. are associated with a number of infectious syndromes in humans including gastroenteritis and dysentery. Our understanding of the genetic diversity, population structure, virulence determinants and antimicrobial resistance of the genus has been limited by a lack of sequenced genomes linked to metadata. We performed a comprehensive analysis of the whole genome sequences of 447 Aeromonas isolates from children in Karachi, Pakistan, with moderate-to-severe diarrhoea (MSD) and from matched controls without diarrhoea that were collected as part of the Global Enteric Multicenter Study (GEMS). Human-associated Aeromonas isolates exhibited high species diversity and extensive antimicrobial and virulence gene content. Aeromonas caviae, A. dhankensis, A. veronii and A. enteropelogenes were all significantly associated with MSD in at least one cohort group. The maf2 and lafT genes that encode components of polar and lateral flagella, respectively, exhibited a weak association with isolates originating from cases of gastroenteritis.
-
-
-
A metagenomic investigation of phytoplasma diversity in Australian vegetable growing regions
In this study, metagenomic sequence data was used to investigate the phytoplasma taxonomic diversity in vegetable-growing regions across Australia. Metagenomic sequencing was performed on 195 phytoplasma-positive samples, originating either from historic collections (n=46) or during collection efforts between January 2015 and June 2022 (n=149). The sampled hosts were classified as crop (n=155), weed (n=24), ornamental (n=7), native plant (n=6), and insect (n=3) species. Most samples came from Queensland (n=78), followed by Western Australia (n=46), the Northern Territory (n=32), New South Wales (n=17), and Victoria (n=10). Of the 195 draft phytoplasma genomes, 178 met our genome criteria for comparison using an average nucleotide identity approach. Ten distinct phytoplasma species were identified and could be classified within the 16SrII, 16SrXII (PCR only), 16SrXXV, and 16SrXXXVIII phytoplasma groups, which have all previously been recorded in Australia. The most commonly detected phytoplasma taxa in this study were species and subspecies classified within the 16SrII group (n=153), followed by strains within the 16SrXXXVIII group (‘Ca. Phytoplasma stylosanthis’; n=6). Several geographic- and host-range expansions were reported, as well as mixed phytoplasma infections of 16SrII taxa and ‘Ca. Phytoplasma stylosanthis’. Additionally, six previously unrecorded 16SrII taxa were identified, including five putative subspecies of ‘Ca. Phytoplasma australasiaticum’ and a new putative 16SrII species. PCR and sequencing of the 16S rRNA gene was a suitable triage tool for preliminary phytoplasma detection. Metagenomic sequencing, however, allowed for higher-resolution identification of the phytoplasmas, including mixed infections, than was afforded by only direct Sanger sequencing of the 16S rRNA gene. Since the metagenomic approach theoretically obtains sequences of all organisms in a sample, this approach was useful to confirm the host family, genus, and/or species. In addition to improving our understanding of the phytoplasma species that affect crop production in Australia, the study also significantly expands the genomic sequence data available in public sequence repositories to contribute to phytoplasma molecular epidemiology studies, revision of taxonomy, and improved diagnostics.
-
- Evolution and Responses to Interventions
-
-
Streptococcus pneumoniae serotype 3 population structure in the era of conjugate vaccines, 2001–2018
Background. Despite use of highly effective conjugate vaccines, invasive pneumococcal disease (IPD) remains a leading cause of morbidity and mortality and disproportionately affects Indigenous populations. Although included in the 13-valent pneumococcal conjugate vaccine (PCV13), which was introduced in 2010, serotype 3 continues to cause disease among Indigenous communities in the Southwest USA. In the Navajo Nation, serotype 3 IPD incidence increased among adults (3.8/100 000 in 2001–2009 and 6.2/100 000 in 2011–2019); in children the disease persisted although the rates dropped from 5.8/100 000 to 2.3/100 000.
Methods. We analysed the genomic epidemiology of serotype 3 isolates collected from 129 adults and 63 children with pneumococcal carriage (n=61) or IPD (n=131) from 2001 to 2018 of the Navajo Nation. Using whole-genome sequencing data, we determined clade membership and assessed changes in serotype 3 population structure over time.
Results. The serotype 3 population structure was characterized by three dominant subpopulations: clade II (n=90, 46.9 %) and clade Iα (n=59, 30.7 %), which fall into Clonal Complex (CC) 180, and a non-CC180 clade (n=43, 22.4 %). The proportion of clade II-associated IPD cases increased significantly from 2001 to 2010 to 2011–2018 among adults (23.1–71.8 %; P<0.001) but not in children (27.3–33.3 %; P=0.84). Over the same period, the proportion of clade II-associated carriage increased; this was statistically significant among children (23.3–52.6 %; P=0.04) but not adults (0–50.0 %, P=0.08).
Conclusions. In this setting with persistent serotype 3 IPD and carriage, clade II has increased since 2010. Genomic changes may be contributing to the observed trends in serotype 3 carriage and disease over time.
-
-
-
Genomic insights into local-scale evolution of ocular Chlamydia trachomatis strains within and between individuals in Gambian trachoma-endemic villages
Trachoma, a neglected tropical disease caused by Chlamydia trachomatis (Ct) serovars A–C, is the leading infectious cause of blindness worldwide. Africa bears the highest burden, accounting for over 86 % of global trachoma cases. We investigated Ct serovar A (SvA) and B (SvB) whole genome sequences prior to the induction of mass antibiotic drug administration in The Gambia. Here, we explore the factors contributing to Ct strain diversification and the implications for Ct evolution within the context of ocular infection. A cohort study in 2002–2003 collected ocular swabs across nine Gambian villages during a 6 month follow-up study. To explore the genetic diversity of Ct within and between individuals, we conducted whole-genome sequencing (WGS) on a limited number (n=43) of Ct-positive samples with an omcB load ≥10 from four villages. WGS was performed using target enrichment with SureSelect and Illumina paired-end sequencing. Out of 43 WGS samples, 41 provided sufficient quality for further analysis. ompA analysis revealed that 11 samples had highest identity to ompA from strain A/HAR13 (NC_007429) and 30 had highest identity to ompA from strain B/Jali20 (NC_012686). While SvB genome sequences formed two distinct village-driven subclades, the heterogeneity of SvA sequences led to the formation of many individual branches within the Gambian SvA subclade. Comparing the Gambian SvA and SvB sequences with their reference strains, Ct A/HAR13 and Ct B/Jali20, indicated an single nucleotide polymorphism accumulation rate of 2.4×10−5 per site per year for the Gambian SvA and 1.3×10−5 per site per year for SvB variants (P<0.0001). Variant calling resulted in a total of 1371 single nucleotide variants (SNVs) with a frequency >25 % in SvA sequences, and 438 SNVs in SvB sequences. Of note, in SvA variants, highest evolutionary pressure was recorded on genes responsible for host cell modulation and intracellular survival mechanisms, whereas in SvB variants this pressure was mainly on genes essential for DNA replication/repair mechanisms and protein synthesis. A comparison of the sequences between observed separate infection events (4–20 weeks between infections) suggested that the majority of the variations accumulated in genes responsible for host–pathogen interaction such as CTA_0166 (phospholipase D-like protein), CTA_0498 (TarP) and CTA_0948 (deubiquitinase). This comparison of Ct SvA and SvB variants within a trachoma endemic population focused on their local evolutionary adaptation. We found a different variation accumulation pattern in the Gambian SvA chromosomal genes compared with SvB, hinting at the potential of Ct serovar-specific variation in diversification and evolutionary fitness. These findings may have implications for optimizing trachoma control and prevention strategies.
-