Search

We introduce the operational genomic unit (OGU) method, a metagenome analysis strategy that directly exploits sequence alignment hits to individual reference genomes as the minimum unit for assessing the diversity of microbial communities and their relevance to environmental factors. This approach is independent of taxonomic classification, granting the possibility of maximal resolution of community composition, and organizes features into an accurate hierarchy using a phylogenomic tree. The outputs are suitable for contemporary analytical protocols for community ecology, differential abundance, and supervised learning while supporting phylogenetic methods, such as UniFrac and phylofactorization, that are seldom applied to shotgun metagenomics despite being prevalent in 16S rRNA gene amplicon studies. As demonstrated in two real-world case studies, the OGU method produces biologically meaningful patterns from microbiome data sets. Such patterns further remain detectable at very low metagenomic sequencing depths. Compared with taxonomic unit-based analyses implemented in currently adopted metagenomics tools, and the analysis of 16S rRNA gene amplicon sequence variants, this method shows superiority in informing biologically relevant insights, including stronger correlation with body environment and host sex on the Human Microbiome Project data set and more accurate prediction of human age by the gut microbiomes of Finnish individuals included in the FINRISK 2002 cohort. We provide Woltka, a bioinformatics tool to implement this method, with full integration with the QIIME 2 package and the Qiita web platform, to facilitate adoption of the OGU method in future metagenomics studies. IMPORTANCE Shotgun metagenomics is a powerful, yet computationally challenging, technique compared to 16S rRNA gene amplicon sequencing for decoding the composition and structure of microbial communities. Current analyses of metagenomic data are primarily based on taxonomic classification, which is limited in feature resolution. To solve these challenges, we introduce operational genomic units (OGUs), which are the individual reference genomes derived from sequence alignment results, without further assigning them taxonomy. The OGU method advances current read-based metagenomics in two dimensions: (i) providing maximal resolution of community composition and (ii) permitting use of phylogeny-aware tools. Our analysis of real-world data sets shows that it is advantageous over currently adopted metagenomic analysis methods and the finest-grained 16S rRNA analysis methods in predicting biological traits. We thus propose the adoption of OGUs as an effective practice in metagenomic studies.

ABSTRACT We introduce the operational genomic unit (OGU) method, a metagenome analysis strategy that directly exploits sequence alignment hits to individual reference genomes as the minimum unit for assessing the diversity of microbial communities and their relevance to environmental factors. This approach is independent of taxonomic classification, granting the possibility of maximal resolution of community composition, and organizes features into an accurate hierarchy using a phylogenomic tree. The outputs are suitable for contemporary analytical protocols for community ecology, differential abundance, and supervised learning while supporting phylogenetic methods, such as UniFrac and phylofactorization, that are seldom applied to shotgun metagenomics despite being prevalent in 16S rRNA gene amplicon studies. As demonstrated in two real-world case studies, the OGU method produces biologically meaningful patterns from microbiome data sets. Such patterns further remain detectable at very low metagenomic sequencing depths. Compared with taxonomic unit-based analyses implemented in currently adopted metagenomics tools, and the analysis of 16S rRNA gene amplicon sequence variants, this method shows superiority in informing biologically relevant insights, including stronger correlation with body environment and host sex on the Human Microbiome Project data set and more accurate prediction of human age by the gut microbiomes of Finnish individuals included in the FINRISK 2002 cohort. We provide Woltka, a bioinformatics tool to implement this method, with full integration with the QIIME 2 package and the Qiita web platform, to facilitate adoption of the OGU method in future metagenomics studies. IMPORTANCE Shotgun metagenomics is a powerful, yet computationally challenging, technique compared to 16S rRNA gene amplicon sequencing for decoding the composition and structure of microbial communities. Current analyses of metagenomic data are primarily based on taxonomic classification, which is limited in feature resolution. To solve these challenges, we introduce operational genomic units (OGUs), which are the individual reference genomes derived from sequence alignment results, without further assigning them taxonomy. The OGU method advances current read-based metagenomics in two dimensions: (i) providing maximal resolution of community composition and (ii) permitting use of phylogeny-aware tools. Our analysis of real-world data sets shows that it is advantageous over currently adopted metagenomic analysis methods and the finest-grained 16S rRNA analysis methods in predicting biological traits. We thus propose the adoption of OGUs as an effective practice in metagenomic studies.

Abstract Background Modern sequencing technologies have generated low-cost microbiome survey datasets, across sample sites, conditions, and treatments, on an unprecedented scale and throughput. These datasets often come with a phylogenetic tree that provides a unique opportunity to examine how shared evolutionary history affects the different patterns in host-associated microbial communities. Results In this paper, we describe an R package, phyloMDA, for phylogeny-aware microbiome data analysis. It includes the Dirichlet-tree multinomial model for multivariate abundance data, tree-guided empirical Bayes estimation of microbial compositions, and tree-based multiscale regression methods with relative abundances as predictors. Conclusion phyloMDA is a versatile and user-friendly tool to analyze microbiome data while incorporating the phylogenetic information and addressing some of the challenges posed by the data.

Abstract Background Advances in DNA sequencing have offered researchers an unprecedented opportunity to better study the variety of species living in and on the human body. However, the analysis of microbiome data is complicated by several challenges. First, the sequencing depth may vary by orders of magnitude across samples. Second, species are rare and the data often contain many zeros. Third, the specimen is a fraction of the microbial ecosystem, and so the data are compositional carrying only relative information. Other characteristics of microbiome data include pronounced over-dispersion in taxon abundances, and the existence of a phylogenetic tree that relates all bacterial species. To address some of these challenges, microbiome analysis workflows often normalize the read counts prior to downstream analysis. However, there are limitations in the current literature on the normalization of microbiome data. Results Under the multinomial distribution for the read counts and a prior for the unknown proportions, we propose an empirical Bayes approach to microbiome data normalization. Using a tree-based extension of the Dirichlet prior, we further extend our method by incorporating the phylogenetic tree into the normalization process. We study the impact of normalization on differential abundance analysis. In the presence of tree structure, we propose a phylogeny-aware detection procedure. Conclusions Extensive simulations and gut microbiome data applications are conducted to demonstrate the superior performance of our empirical Bayes method over other normalization methods, and over commonly-used methods for differential abundance testing. Original R scripts are available at GitHub ( https://github.com/liudoubletian/eBay ).

Tintinnid ciliates build loricae, whose structure, shape, and size still largely represent the basis for taxonomy and classification, although genetic analyses demonstrated their limited utility for inferring evolutionary relationships. The textures of the lorica walls, however, result from the chemical and physical properties of the forming material, which is supposed to be rather conserved in closely related taxa, viz., congeners and confamilial genera. Within a particular texture, small deviations in the chemical composition might affect the wall's stickiness and accordingly its capability to adhere foreign particles, explaining the intertwining of tintinnids with hyaline and agglutinated loricae in phylogenetic inferences. In a comprehensive comparative study, the lorica textures were electron microscopically and morphometrically analyzed in 21 species from 17 genera and more than nine families together with literature data. Most species were investigated for the first time, and the taxa cover a substantial portion of the molecular genealogy. The phylogeny‐aware analysis of the lorica‐related features provides a preliminary hypothesis on lorica evolution. Eventually, this conspectus suggests the dominance of hard lorica walls with an alveolar texture and proposes different modes of lorica formation. [ABSTRACT FROM AUTHOR]

The COVID-19 pandemic caused by SARS-CoV-2 has led to a wide range of clinical presentations, with respiratory symptoms being common. However, emerging evidence suggests that the gastrointestinal (GI) tract is also affected, with angiotensin-converting enzyme 2, a key receptor for SARS-CoV-2, abundantly expressed in the ileum and colon. The virus has been detected in GI tissues and fecal samples, even in cases with negative results of the reverse transcription polymerase chain reaction in the respiratory tract. GI symptoms have been associated with an increased risk of ICU admission and mortality. The gut microbiome, a complex ecosystem of around 40 trillion bacteria, plays a crucial role in immunological and metabolic pathways. Dysbiosis of the gut microbiota, characterized by a loss of beneficial microbes and decreased microbial diversity, has been observed in COVID-19 patients, potentially contributing to disease severity. We conducted a comprehensive gut microbiome study in 204 hospitalized COVID-19 patients using both shallow and deep shotgun sequencing methods. We aimed to track microbiota composition changes induced by hospitalization, link these alterations to clinical procedures (antibiotics administration) and outcomes (ICU referral, survival), and assess the predictive potential of the gut microbiome for COVID-19 prognosis. Shallow shotgun sequencing was evaluated as a cost-effective diagnostic alternative for clinical settings. Our study demonstrated the diverse effects of various combinations of clinical parameters, microbiome profiles, and patient metadata on the precision of outcome prognostication in patients. It indicates that microbiological data possesses greater reliability in forecasting patient outcomes when contrasted with clinical data or metadata. Furthermore, we established that shallow shotgun sequencing presents a viable and costeffective diagnostic alternative to deep sequencing within clinical environments. [ABSTRACT FROM AUTHOR]

Phylogenetic trees provide a framework for organizing evolutionary histories across the tree of life and aid downstream comparative analyses such as metagenomic identification. Methods that rely on single-marker genes such as 16S rRNA have produced trees of limited accuracy with hundreds of thousands of organisms, whereas methods that use genome-wide data are not scalable to large numbers of genomes. We introduce updating trees using divide-and-conquer (uDance), a method that enables updatable genome-wide inference using a divide-and-conquer strategy that refines different parts of the tree independently and can build off of existing trees, with high accuracy and scalability. With uDance, we infer a species tree of roughly 200,000 genomes using 387 marker genes, totaling 42.5 billion amino acid residues. Large, updatable phylogenetic trees are constructed using a divide-and-conquer strategy. [ABSTRACT FROM AUTHOR]

Studies using 16S rRNA and shotgun metagenomics typically yield different results, usually attributed to PCR amplification biases. We introduce Greengenes2, a reference tree that unifies genomic and 16S rRNA databases in a consistent, integrated resource. By inserting sequences into a whole-genome phylogeny, we show that 16S rRNA and shotgun metagenomic data generated from the same samples agree in principal coordinates space, taxonomy and phenotype effect size when analyzed with the same tree. A comprehensive microbial resource reconciles genomic and 16S rRNA data in a single tree. [ABSTRACT FROM AUTHOR]

Motivation Taxonomic classification of short reads and taxonomic profiling of metagenomic samples are well-studied yet challenging problems. The presence of species belonging to groups without close representation in a reference dataset is particularly challenging. While k -mer-based methods have performed well in terms of running time and accuracy, they tend to have reduced accuracy for such novel species. Thus, there is a growing need for methods that combine the scalability of k -mers with increased sensitivity. Results Here, we show that using locality-sensitive hashing (LSH) can increase the sensitivity of the k -mer-based search. Our method, which combines LSH with several heuristics techniques including soft lowest common ancestor labeling and voting, is more accurate than alternatives in both taxonomic classification of individual reads and abundance profiling. Availability and implementation CONSULT-II is implemented in C ++ , and the software, together with reference libraries, is publicly available on GitHub https://github.com/bo1929/CONSULT-II. [ABSTRACT FROM AUTHOR]

The brain–gut axis forms a bidirectional communication system between the gastrointestinal (GI) tract and cognitive brain areas. Disturbances to this system in disease states such as inflammatory bowel disease have consequences for neuronal activity and subsequent cognitive function. The gut–microbiota–brain axis refers to the communication between gut-resident bacteria and the brain. This circuits exists to detect gut microorganisms and relay information to specific areas of the central nervous system (CNS) that in turn, regulate gut physiology. Changes in both the stability and diversity of the gut microbiota have been implicated in several neuronal disorders, including depression, autism spectrum disorder Parkinson's disease, Alzheimer's disease and multiple sclerosis. Correcting this imbalance with medicinal herbs, the metabolic products of dysregulated bacteria and probiotics have shown hope for the treatment of these neuronal disorders. In this review, we focus on recent advances in our understanding of the intricate connections between the gut–microbiota and the brain. We discuss the contribution of gut microbiota to neuronal disorders and the tangible links between diseases of the GI tract with cognitive function and behaviour. In this regard, we focus on irritable bowel syndrome (IBS) given its strong links to brain function and anxiety disorders. This adds to the growing body of evidence supporting targeted therapeutic strategies to modulate the gut microbiota for the treatment of brain/mental-health-related disease. [ABSTRACT FROM AUTHOR]

The use of museum specimens for research in microbial evolutionary ecology remains an under-utilized investigative dimension with important potential. Despite this potential, there remain barriers in methodology and analysis to the wide-spread adoption of museum specimens for such studies. Here, we hypothesized that there would be significant differences in taxonomic prediction and related diversity among sample type (museum or fresh) and sequencing strategy (medium-depth shotgun metagenomic or 16S rRNA gene). We found dramatically higher predicted diversity from shotgun metagenomics when compared to 16S rRNA gene sequencing in museum and fresh samples, with this differential being larger in museum specimens. Broadly confirming these hypotheses, the highest diversity found in fresh samples was with shotgun sequencing using the Rep200 reference inclusive of viruses and microeukaryotes, followed by the WoL reference database. In museum-specimens, community diversity metrics also differed significantly between sequencing strategies, with the alpha-diversity ACE differential being significantly greater than the same comparisons made for fresh specimens. Beta diversity results were more variable, with significance dependent on reference databases used. Taken together, these findings demonstrate important differences in diversity results and prompt important considerations for future experiments and downstream analyses aiming to incorporate microbiome datasets from museum specimens. [ABSTRACT FROM AUTHOR]

Background: Human milk is unquestionably beneficial for preterm infants. We investigated how the transition from tube to oral/breastfeeding impacts the preterm infants' oral and gut microbiome and metabolome., Methods: We analyzed stool, saliva, and milk samples collected from a cohort of preterm infants enrolled in the MAP Study, a prospective observational trial. The microbiome and metabolome of the samples were analyzed from 4 longitudinal sample time points, 2 during tube feeds only and 2 after the initiation of oral/breastfeeding., Results: We enrolled 11 mother-infant dyads (gestational age = 27.9 (23.4-32.2)) and analyzed a total of 39 stool, 44 saliva, and 43 milk samples over 4 timepoints. In saliva samples, there was a shift towards increased Streptococcus and decreased Staphylococcus after oral feeding/breastfeeding initiation (p < 0.05). Milk sample metabolites were strongly influenced by the route of feeding and milk type (p < 0.05) and represented the pathways of Vitamin E metabolism, Vitamin B12 metabolism, and Tryptophan metabolism., Conclusion: Our analysis demonstrated that the milk and preterm infant's saliva microbiome and metabolome changed over the course of the first four to 5 months of life, coinciding with the initiation of oral/breastfeeds., Impact: The microbiome and metabolome is altered in the infant's saliva but not their stool, and in mother's milk when feeds are transitioned from tube to oral/breastfeeding. We assessed the relationship between the gut and oral microbiome/metabolome with the milk microbiome/metabolome over a longitudinal period of time in preterm babies. Metabolites that changed in the infants saliva after the initiation of oral feeds have the potential to be used as biomarkers for disease risk., (© 2024. The Author(s).)

The COVID-19 pandemic, caused by the coronavirus SARS-CoV-2, has generated a renewed interest in the larger family of Coronaviridae, which causes a variety of different respiratory infections in a variety of different hosts. Understanding the mechanisms behind the ability of a family of viruses to spill over into different hosts is an ongoing study. In this work, we studied the relationship between specific amino acid sites and the solvent accessibility of the surface (or spike) protein of different Coronaviridae. Since host specificity hinges on the portion(s) of the protein that interfaces with the host cell membrane, there could be a relationship between information gain in specific amino acid sites and solvent accessibility. We found a connection between sites with high information gain and solvent accessibility within several major subgenera of Coronaviridae. Such a connection could be used to study other lesser-known families of viruses, which is desirable because information gain is much easier to compute when the number of sequences is large, as we show. Finally, we produced a visualization of the sequences within each major subgenus and discussed several regions of interest, as well as focused on some pairs of Coronaviridae hosts of interest. [ABSTRACT FROM AUTHOR]

Periodontal disease is more common in individuals with rheumatoid arthritis (RA) who have detectable anti-citrullinated protein antibodies (ACPAs), implicating oral mucosal inflammation in RA pathogenesis. Here, we performed paired analysis of human and bacterial transcriptomics in longitudinal blood samples from RA patients. We found that patients with RA and periodontal disease experienced repeated oral bacteremias associated with transcriptional signatures of ISG15+HLADRhi and CD48highS100A2pos monocytes, recently identified in inflamed RA synovia and blood of those with RA flares. The oral bacteria observed transiently in blood were broadly citrullinated in the mouth, and their in situ citrullinated epitopes were targeted by extensively somatically hypermutated ACPAs encoded by RA blood plasmablasts. Together, these results suggest that (i) periodontal disease results in repeated breaches of the oral mucosa that release citrullinated oral bacteria into circulation, which (ii) activate inflammatory monocyte subsets that are observed in inflamed RA synovia and blood of RA patients with flares and (iii) activate ACPA B cells, thereby promoting affinity maturation and epitope spreading to citrullinated human antigens. Periodontal disease and rheumatoid arthritis: The incidence of periodontal disease is high in individuals with rheumatoid arthritis (RA) who also have anti-citrullinated protein antibodies (ACPAs), suggesting a link between these two diseases. Brewer et al. carried out a paired analysis of human and bacterial transcriptomes from blood samples collected longitudinally from RA patients with and without periodontal disease. They identified transcriptional signatures within inflammatory monocyte subsets that correlated with repeated oral bacteremias and clinical arthritis flares in patients with RA and periodontal disease. These oral bacteria were broadly citrullinated, and some of these citrullinated epitopes were the targets of ACPA expressed by RA blood plasmablasts that have undergone affinity maturation. These results confirm that periodontal disease can cause breaches in oral mucosa that release citrullinated bacteria into the blood, which activates inflammatory monocytes and ACPA-specific B cells. —CF [ABSTRACT FROM AUTHOR]

Cow's milk is frequently included in the human diet, but the relationship between milk intake and type 2 diabetes (T2D) remains controversial. Here, using data from the Hispanic Community Health Study/Study of Latinos, we show that in both sexes, higher milk intake is associated with lower risk of T2D in lactase non-persistent (LNP) individuals (determined by a variant of the lactase LCT gene, single nucleotide polymorphism rs4988235 ) but not in lactase persistent individuals. We validate this finding in the UK Biobank. Further analyses reveal that among LNP individuals, higher milk intake is associated with alterations in gut microbiota (for example, enriched Bifidobacterium and reduced Prevotella) and circulating metabolites (for example, increased indolepropionate and reduced branched-chain amino acid metabolites). Many of these metabolites are related to the identified milk-associated bacteria and partially mediate the association between milk intake and T2D in LNP individuals. Our study demonstrates a protective association between milk intake and T2D among LNP individuals and a potential involvement of gut microbiota and blood metabolites in this association., (© 2024. The Author(s), under exclusive licence to Springer Nature Limited.)

Background: Crohn's disease (CD), an inflammatory bowel disease (IBD) subtype, results from pathologic interactions between host cells and its resident gut microbes. CD manifests in both isolated disease locations (ileum or colon) or a combination of locations (ileocolonic). To date, a comprehensive understanding of how isolated CD subtypes influence molecular profiles remains outstanding. To address this, we sought to define CD location signatures by leveraging a large cross-sectional feature set captured from the stool of over 200 IBD patients and healthy controls using metaproteomics, shotgun metagenomics, 16S rRNA sequencing, metabolomic profiling, and host genetics paired with clinical endoscopic assessments. Results: Neither metagenomic nor host genetics alone distinguished CD location subtypes. In contrast, ileal and colonic CD were distinguished using mass spectrometry-based methods (metabolomics or metaproteomics) or a combined multi-omic feature set. This multi-omic feature set revealed colonic CD was strongly associated with neutrophil-related proteins. Additionally, colonic CD displayed a disease-severity-related association with Bacteroides vulgatus. Colonic CD and ulcerative colitis profiles harbored strikingly similar feature enrichments compared to ileal CD, including neutrophil-related protein enrichments. Compared to colonic CD, ileal CD profiles displayed increased primary and secondary bile acid levels and concomitant shifts in taxa with noted sensitivities such as Faecalibacterium prausnitzii or affinities for bile acid-rich environments, including Gammaproteobacteria and Blautia sp. Having shown robust molecular and microbial distinctions tied to CD locations, we leveraged these profiles to generate location-specific disease severity biomarkers that surpass the performance of Calprotectin. Conclusions: When compared using multi-omics features, colonic- and ileal-isolated CD subtypes display striking differences that suggest separate location-specific pathologies. Colonic CD's strong similarity to ulcerative colitis, including neutrophil and Bacteroides vulgatus involvement, is also evidence of a shared pathology for colonic-isolated IBD subtypes, while ileal CD maintains a unique, bile acid-driven profile. More broadly, this study demonstrates the power of multi-omics approaches for IBD biomarker discovery and elucidating the underlying biology. E_-RTD2us73nf2wD8LJ9fb Video Abstract [ABSTRACT FROM AUTHOR]

Preterm infants are at a greater risk for the development of asthma and atopic disease, which can lead to lifelong negative health consequences. This may be due, in part, to alterations that occur in the gut microbiome and metabolome during their stay in the Neonatal Intensive Care Unit (NICU). To explore the differential roles of family history (i.e., predisposition due to maternal asthma diagnosis) and hospital-related environmental and clinical factors that alter microbial exposures early in life, we considered a unique cohort of preterm infants born ≤ 34 weeks gestational age from two local level III NICUs, as part of the MAP (Microbiome, Atopic disease, and Prematurity) Study. From MAP participants, we chose a sub-cohort of infants whose mothers had a history of asthma and matched gestational age and sex to infants of mothers without a history of asthma diagnosis (control). We performed a prospective, paired metagenomic and metabolomic analysis of stool and milk feed samples collected at birth, 2 weeks, and 6 weeks postnatal age. Although there were clinical factors associated with shifts in the diversity and composition of stool-associated bacterial communities, maternal asthma diagnosis did not play an observable role in shaping the infant gut microbiome during the study period. There were significant differences, however, in the metabolite profile between the maternal asthma and control groups at 6 weeks postnatal age. The most notable changes occurred in the linoleic acid spectral network, which plays a role in inflammatory and immune pathways, suggesting early metabolomic changes in the gut of preterm infants born to mothers with a history of asthma. Our pilot study suggests that a history of maternal asthma alters a preterm infants' metabolomic pathways in the gut, as early as the first 6 weeks of life. [ABSTRACT FROM AUTHOR]

Autism spectrum disorder (ASD) is a neurodevelopmental disorder characterized by heterogeneous cognitive, behavioral and communication impairments. Disruption of the gut-brain axis (GBA) has been implicated in ASD although with limited reproducibility across studies. In this study, we developed a Bayesian differential ranking algorithm to identify ASD-associated molecular and taxa profiles across 10 cross-sectional microbiome datasets and 15 other datasets, including dietary patterns, metabolomics, cytokine profiles and human brain gene expression profiles. We found a functional architecture along the GBA that correlates with heterogeneity of ASD phenotypes, and it is characterized by ASD-associated amino acid, carbohydrate and lipid profiles predominantly encoded by microbial species in the genera Prevotella, Bifidobacterium, Desulfovibrio and Bacteroides and correlates with brain gene expression changes, restrictive dietary patterns and pro-inflammatory cytokine profiles. The functional architecture revealed in age-matched and sex-matched cohorts is not present in sibling-matched cohorts. We also show a strong association between temporal changes in microbiome composition and ASD phenotypes. In summary, we propose a framework to leverage multi-omic datasets from well-defined cohorts and investigate how the GBA influences ASD., (© 2023. The Author(s).)

Interacting proteins and protein domains coevolve on multiple scales, from their correlated presence across species, to correlations in amino-acid usage. Genomic databases provide rapidly growing data for variability in genomic protein content and in protein sequences, calling for computational predictions of unknown interactions. We first introduce the concept of direct phyletic couplings, based on global statistical models of phylogenetic profiles. They strongly increase the accuracy of predicting pairs of related protein domains beyond simpler correlation-based approaches like phylogenetic profiling (80% vs. 30–50% positives out of the 1000 highest-scoring pairs). Combined with the direct coupling analysis of inter-protein residue-residue coevolution, we provide multi-scale evidence for direct but unknown interaction between protein families. An in-depth discussion shows these to be biologically sensible and directly experimentally testable. Negative phyletic couplings highlight alternative solutions for the same functionality, including documented cases of convergent evolution. Thereby our work proves the strong potential of global statistical modeling approaches to genome-wide coevolutionary analysis, far beyond the established use for individual protein complexes and domain-domain interactions. Interactions between proteins and their domains are at the basis of almost all biological processes. To complement labor intensive and error-prone experimental approaches to the genome-scale characterization of such interactions, we propose a computational approach based upon rapidly growing protein-sequence databases. To maintain interaction in the course of evolution, proteins and their domains are required to coevolve: evolutionary changes in the interaction partners appear correlated across several scales, from correlated presence-absence patterns of proteins across species, up to correlations in the amino-acid usage. Our approach combines these different scales within a common mathematical-statistical inference framework, which is inspired by the so-called direct coupling analysis. It is able to predict currently unknown, but biologically sensible interaction, and to identify cases of convergent evolution leading to alternative solutions for a common biological task. Thereby our work illustrates the potential of global statistical inference for the genome-scale coevolutionary analysis of interacting proteins and protein domains. [ABSTRACT FROM AUTHOR]

Background: Deforestation in the Amazon and the social vulnerability of its settler communities has been associated with increased malaria incidence. The feeding biology of the most important malaria vectors in the region, notably Nyssorhynchus darlingi, compounds efforts to control vectors and reduce transmission of what has become known as "Frontier Malaria". Exploring Anophelinae mosquito diversity is fundamental to understanding the species responsible for transmission and developing appropriate management and intervention strategies for malaria control in the Amazon River basin. Methods: This study describes Anophelinae mosquito diversity from settler communities affected by Frontier Malaria in the states of Acre, Amazonas and Rondônia by analysing COI gene data using cluster and tree-based species delimitation approaches. Results: In total, 270 specimens from collection sites were sequenced and these were combined with 151 reference (GenBank) sequences in the analysis to assist in species identification. Conservative estimates found that the number of species collected at these sites was between 23 (mPTP partition) and 27 (strict ABGD partition) species, up to 13 of which appeared to be new. Nyssorhynchus triannulatus and Nyssorhynchus braziliensis displayed exceptional levels of intraspecific genetic diversity but there was little to no support for putative species complex status. Conclusions: This study demonstrates that Anophelinae mosquito diversity continues to be underestimated in poorly sampled areas where frontier malaria is a major public health concern. The findings will help shape future studies of vector incrimination and transmission dynamics in these areas and support efforts to develop more effective vector control and transmission reduction strategies in settler communities in the Amazon River basin. [ABSTRACT FROM AUTHOR]