Search NASA⌕ Search

SEARCH · Search NASA

Results for “phylogenetics”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 19 records

Phylogenetic structure of specialization: A new approach that integrates partner availability and phylogenetic diversity to quantify biotic specialization in ecological networks

Abstract Biotic specialization holds information about the assembly, evolution, and stability of biological communities. Partner availabilities can play an important role in enabling species interactions, where uneven partner availabilities can bias estimates of biotic specialization when using phylogenetic diversity indices. It is therefore important to account for partner availability when characterizing biotic specialization using phylogenies. We developed an index, phylogenetic structure of specialization (PSS), that avoids bias from uneven partner availabilities by uncoupling the null models for interaction frequency and phylogenetic distance. We incorporate the deviation between observed and random interaction frequencies as weights into the calculation of partner phylogenetic α‐diversity. To calculate the PSS index, we then compare observed partner phylogenetic α‐diversity to a null distribution generated by randomizing phylogenetic distances among the same number of partners. PSS quantifies the phylogenetic structure (i.e., clustered, overdispersed, or random) of the partners of a focal species. We show with simulations that the PSS index is not correlated with network properties, which allows comparisons across multiple systems. We also implemented PSS on empirical networks of host–parasite, avian seed‐dispersal, lichenized fungi–cyanobacteria, and hummingbird pollination interactions. Across these systems, a large proportion of taxa interact with phylogenetically random partners according to PSS, sometimes to a larger extent than detected with an existing method that does not account for partner availability. We also found that many taxa interact with phylogenetically clustered partners, while taxa with overdispersed partners were rare. We argue that species with phylogenetically overdispersed partners have often been misinterpreted as generalists when they should be considered specialists. Our results highlight the important role of randomness in shaping interaction networks, even in highly intimate symbioses, and provide a much‐needed quantitative framework to assess the role that evolutionary history and symbiotic specialization play in shaping patterns of biodiversity. PSS is available as an R package at https://github.com/cjpardodelahoz/pss .

59 BASIC BIOLOGICAL SCIENCES↗

Rare species do not disproportionately contribute to phylogenetic diversity in a subalpine plant community

Abstract Premise Within plant communities, few species are abundant, and most are locally rare. Worldwide, 36% of plant species are exceedingly rare and often face high extinction risk. However, the community phylogenetic impact of the loss of rare plants is largely unknown in many systems. We address this gap by investigating how rare species contribute to phylogenetic diversity, considering multiple metrics of rarity and multiple elevations in a subalpine plant community. Methods We collected abundance data at three sites near the Rocky Mountain Biological Laboratory (Colorado, USA). We calculated each species' range size from public occurrence data. We calculated phylogenetic signal for abundance and range size, compared community phylogenetic metrics weighted by range size and abundance to unweighted metrics, and quantified the change in phylogenetic diversity when removing single species and groups of species ranked by rarity. Results We found phylogenetic signal for abundance, but not range size. There was no difference between rarity‐weighted and ‐unweighted phylogenetic diversity metrics. Finally, phylogenetic diversity did not decline more when we removed single rare species or groups of rare species than when we removed single common species and groups of common species. Conclusions We found that rare species, whether at low abundance or with a small range, do not disproportionately contribute to phylogenetic diversity in our subalpine plant community. These results were consistent across elevations. Instead, rare species might provide phylogenetic redundancy with common species. Deeper understanding of functional differentiation is needed to understand contributions of rare species to this system.

Plant Sciences↗

The ecological assembly of bacterial communities in Antarctic wetlands varies across levels of phylogenetic resolution

Summary As functional traits are conserved at different phylogenetic depths, the ability to detect community assembly processes can be conditional on the phylogenetic resolution; yet most previous work quantifying their influence has focused on a single level of phylogenetic resolution. Here, we have studied the ecological assembly of bacterial communities from an Antarctic wetland complex, applying null models across different levels of phylogenetic resolution (i.e. clustering ASVs into OTUs with decreasing sequence identity thresholds). We found that the relative influence of the community assembly processes varies with phylogenetic resolution. More specifically, selection processes seem to impose stronger influence at finer (100% sequence similarity ASV) than at coarser (99%–97% sequence similarity OTUs) resolution. We identified environmental features related with the ecological processes and propose a conceptual model for the bacterial community assembly in this Antarctic ecosystem. Briefly, eco‐evolutionary processes appear to be leading to different but very closely related ASVs in lotic, lentic and terrestrial environments. In all, this study shows that assessing community assembly processes at different phylogenetic resolutions is key to improve our understanding of microbial ecology. More importantly, a failure to detect selection processes at coarser phylogenetic resolution does not imply the absence of such processes at finer resolutions.

59 BASIC BIOLOGICAL SCIENCES↗

Full ribosomal operon sequencing of anaerobic gut fungi (phylum Neocallimastigomycota ): insights on its markers and phylogenetic resolution

The phylogenetic affiliations of anaerobic gut fungi (Neocallimastigomycota) are typically evaluated using single-gene markers. However, this approach often fails to resolve relationships between closely related lineages. To address this issue and identify alternative markers, we created a curated database comprising the complete ribosomal operon sequences of 156 isolates, representing 20 of the 22 recognized genera and two new genus-level clades. Using long-read sequencing, we obtained ~9 kbp operon sequences and developed a robust analysis pipeline. Incorporating both coding genes and non-coding regions (excluding IGS1) improved phylogenetic resolution. This phylogenetic approach successfully resolved the Cyllamyces and Caecomyces clades (hard-to-distinguish genetically), as well as seven analysed Piromyces species. We also scanned the operon for markers that are suitable for short-read sequencing platforms, with the aim of enhancing biodiversity and phylogenetic studies. Notably, the ETS1 genetic region also enabled the distinction between these lineages, indicating its phylogenetic value within the ribosomal operon. The resulting database is a valuable resource for expanding and strengthening phylogenetic frameworks.

High-throughput sequencing↗

Phylogenetic Diversity and Mycotoxin Potential of Emergent Phytopathogens Within the Fusarium tricinctum Species Complex

Recent studies on multiple continents indicate members of the Fusarium tricinctum species complex (FTSC) are emerging as prevalent pathogens of small-grain cereals, pulses, and other economically important crops. These understudied fusaria produce structurally diverse mycotoxins, among which enniatins (ENNs) and moniliformin (MON) are the most frequent and of greatest concern to food and feed safety. Herein a large survey of fusaria in the Fusarium Research Center and Agricultural Research Service culture collections was undertaken to assess species diversity and mycotoxin potential within the FTSC. A 151-strain collection originating from diverse hosts and substrates from different agroclimatic regions throughout the world was selected from 460 FTSC strains to represent the breadth of FTSC phylogenetic diversity. Evolutionary relationships inferred from a five-locus dataset, using maximum likelihood and parsimony, resolved the 151 strains as 24 phylogenetically distinct species, including nine that are new to science. Of the five genes analyzed, nearly full-length phosphate permease sequences contained the most phylogenetically informative characters, establishing its suitability for species-level phylogenetics within the FTSC. Fifteen of the species produced ENNs, MON, the sphingosine analog 2-amino-14,16-dimethyloctadecan-3-ol (AOD), and the toxic pigment aurofusarin (AUR) on a cracked corn kernel substrate. Interestingly, the five earliest diverging species in the FTSC phylogeny (i.e., F. iranicum, F. flocciferum, F. torulosum, and Fusarium spp. FTSC 8 and 24) failed to produce AOD and MON, but synthesized ENNs and/or AUR. Moreover, our reassessment of nine published phylogenetic studies on the FTSC identified 11 additional novel taxa, suggesting this complex comprises at least 36 species.

Plant Sciences↗

Phylogenetic organization in the assimilation of chemically distinct substrates by soil bacteria

Soils are among the most biodiverse habitats on earth and while the species composition of microbial communities can influence decomposition rates and pathways, the functional significance of many microbial species and phylogenetic groups remains unknown. If bacteria exhibit phylogenetic organization in their function, this could enable ecologically meaningful classification of bacterial clades. We show non-random phylogenetic organization in the rates of relative carbon assimilation for both rapidly mineralized substrates (amino acids and glucose) assimilated by many microbial taxa and slowly mineralized substrates (lipids and cellulose) assimilated by relatively few microbial taxa. When mapped onto bacterial phylogeny using ancestral character estimation this phylogenetic organization enabled the identification of clades involved in the decomposition of specific soil organic matter substrates. Phylogenetic organization in substrate assimilation could provide a basis for predicting the functional attributes of uncharacterized microbial taxa and understanding the significance of microbial community composition for soil organic matter decomposition.

59 BASIC BIOLOGICAL SCIENCES↗

Phylogenetic conservatism of mycoparasitism and its contribution to pathogen antagonism

Abstract Closely related species are expected to have similar functional traits due to shared ancestry and phylogenetic inertia. However, few tests of this hypothesis are available for plant‐associated fungal symbionts. Fungal leaf endophytes occur in all land plants and can protect their host plant from disease by a variety of mechanisms, including by parasitizing pathogens (e.g., mycoparasitism). Here, we tested whether phylogenetic relatedness among species of Cladosporium , a widespread genus that includes mycoparasitic species, predicts the effect of this endophyte on the severity of leaf rust disease. First, we used congruence among different marker sequences (i.e., genealogical concordance phylogenetic species recognition criterion) to delimit species of Cladosporium . Next, in a controlled experiment, we quantified both mycoparasitism and disease modification for the selected Cladosporium species. We identified 17 species of Cladosporium ; all the species reduced rust disease severity in our experiment. Cladosporium phylogeny was a significant predictor of mycoparasitism. However, we did not observe a phylogenetic effect on disease severity overall, indicating that other mechanism/s operating independently of shared ancestry also contributed to endophyte effects on disease severity. Indeed, a second experiment showed that Cladosporium endophyte exudates (no live organism) from divergent species groups equally reduced disease severity. Our results reveal that multiple mechanisms contribute to the protective effects of an endophyte against a plant pathogen, but not all traits underlying these mechanisms are phylogenetically conserved.

Barge, Edward G.↗

Islands Within Islands: Bacterial Phylogenetic Structure and Consortia in Hawaiian Lava Caves and Fumaroles

Lava caves, tubes, and fumaroles in Hawai‘i present a range of volcanic, oligotrophic environments from different lava flows and host unexpectedly high levels of bacterial diversity. These features provide an opportunity to study the ecological drivers that structure bacterial community diversity and assemblies in volcanic ecosystems and compare the older, more stable environments of lava tubes, to the more variable and extreme conditions of younger, geothermally active caves and fumaroles. Using 16S rRNA amplicon-based sequencing methods, we investigated the phylogenetic distinctness and diversity and identified microbial interactions and consortia through co-occurrence networks in 70 samples from lava tubes, geothermal lava caves, and fumaroles on the island of Hawai‘i. Our data illustrate that lava caves and geothermal sites harbor unique microbial communities, with very little overlap between caves or sites. We also found that older lava tubes (500–800 yrs old) hosted greater phylogenetic diversity (Faith's PD) than sites that were either geothermally active or younger (<400 yrs old). Geothermally active sites had a greater number of interactions and complexity than lava tubes. Average phylogenetic distinctness, a measure of the phylogenetic relatedness of a community, was higher than would be expected if communities were structured at random. This suggests that bacterial communities of Hawaiian volcanic environments are phylogenetically over-dispersed and that competitive exclusion is the main driver in structuring these communities. This was supported by network analyses that found that taxa (Class level) co-occurred with more distantly related organisms than close relatives, particularly in geothermal sites. Network “hubs” (taxa of potentially higher ecological importance) were not the most abundant taxa in either geothermal sites or lava tubes and were identified as unknown families or genera of the phyla, Chloroflexi and Acidobacteria. These results highlight the need for further study on the ecological role of microbes in caves through targeted culturing methods, metagenomics, and long-read sequence technologies.

59 BASIC BIOLOGICAL SCIENCES↗

Variation in Sphingomonas traits across habitats and phylogenetic clades

Whether microbes show habitat preferences is a fundamental question in microbial ecology. If different microbial lineages have distinct traits, those lineages may occur more frequently in habitats where their traits are advantageous. Sphingomonas is an ideal bacterial clade in which to investigate how habitat preference relates to traits because these bacteria inhabit diverse environments and hosts. Here we downloaded 440 publicly available Sphingomonas genomes, assigned them to habitats based on isolation source, and examined their phylogenetic relationships. We sought to address whether: (1) there is a relationship between Sphingomonas habitat and phylogeny, and (2) whether there is a phylogenetic correlation between key, genome-based traits and habitat preference. We hypothesized that Sphingomonas strains from similar habitats would cluster together in phylogenetic clades, and key traits that improve fitness in specific environments should correlate with habitat. Genome-based traits were categorized into the Y-A-S trait-based framework for high growth yield, resource acquisition, and stress tolerance. We selected 252 high quality genomes and constructed a phylogenetic tree with 12 well-defined clades based on an alignment of 404 core genes. Sphingomonas strains from the same habitat clustered together within the same clades, and strains within clades shared similar clusters of accessory genes. Additionally, key genome-based trait frequencies varied across habitats. We conclude that Sphingomonas gene content reflects habitat preference. This knowledge of how environment and host relate to phylogeny may also help with future functional predictions about Sphingomonas and facilitate applications in bioremediation.

59 BASIC BIOLOGICAL SCIENCES↗

Phylogenetic analyses and reclassification of the oleaginous marine species Nannochloris sp. “desiccata" (Trebouxiophyceae, Chlorophyta), formerly Chlorella desiccata , supported by a high-quality genome assembly

Microalgae are diverse, with many gaps remaining in phylogenetic and physiological understanding. Thus, studying new microalgae species increases our broader comprehension of biological diversity, and evaluation of new candidates as algal production platforms can lead to improved productivity under a variety of cultivation conditions. Chlorella is a genus of fast-growing species often isolated from freshwater habitats and cultivated as a source of nutritional supplements. However, the use of freshwater increases competition with other freshwater needs. We identified Chlorella desiccata to be worthy of further investigation as a potential algae production strain, due to its isolation from a marine environment and its promising growth and biochemical composition properties. Long-read genomic sequencing was conducted for C. desiccata UTEX 2526, resulting in a high-quality, near chromosome level, diploid genome with an assembly length of 21.55 Mbp in only 18 contigs. We also report complete circular mitochondrial and chloroplast genomes. Phylogenomic and phylogenetic analyses using nuclear, chloroplast, 18S rRNA, and actin sequences revealed that this species clades within strains currently identified as Nannochloris (Trebouxiophyceae, Chlorophyta), leading to its reclassification as Nannochloris sp. “desiccata” UTEX 2526. The mode of cell division for this species is autosporulation, differing from the type species N. bacillaris. As has occurred across multiple microalgae genera, there are repeated examples of Nannochloris species reclassification in the literature. This high-quality genome assembly and phylogenetic analysis of the potential algal production strain Nannochloris sp. “desiccata” UTEX 2526 provides an important reference and useful tool for further studying this region of the phylogenetic tree.

59 BASIC BIOLOGICAL SCIENCES↗

Molecular phylogenetic analysis of the Paguristes tortugae Schmitt, 1933 complex and selected other Paguroidea (Crustacea: Decapoda: Anomura)

Morphological characters, as presently applied to describe members of the Paguristes tortugae Schmitt, 1933 species complex, appear to be of limited value in inferring phylogenetic relationships within the genus, and may have similarly misinformed understanding of relationships between members of this complex and those presently assigned to the related genera Areopaguristes Rahayu & McLaughlin, 2010 and Pseudopaguristes McLaughlin, 2002. Previously undocumented observations of similarities and differences in color patterns among populations additionally suggest genetic divergences within some species, or alternatively seem to support phylogenetic groupings of some species. In the present study, a Maximum Likelihood (ML) phylogenetic analysis was undertaken based on the H3, 12S mtDNA, and 16S mtDNA sequences of 148 individuals, primarily representatives of paguroid species from the western Atlantic. This molecular analysis supported a polyphyletic Diogenidae Ortmann, 1892, although incomplete taxonomic sampling among the genera of Diogenidae limits the utility of this finding for resolving family level relationships. Several hypotheses regarding the evolutionary relationships among hermit crab genera were refuted by the Kishino-Hasegawa (KH). ShimodairaHasegawa (SH) and Approximately Unbiased (AU) tree topology tests, among them the hypothesis that Areopaguristes is monophyletic. A lack of support for the monophyly of Areopaguristes calls into question the phylogenetic validity of gill number for the differentiation of Paguristes, Areopaguristes, and Pseudopaguristes. The study was inconclusive with regard to the relationships among these three genera, but previously unknown diversity within both Paguristes and Areopaguristes was demonstrated. Existence of an undescribed species confounded under the name Paguristes tortugae Schmitt, 1933 was supported by genetics, morphology, and coloration. A second undescribed species with remarkable similarity to Areopaguristes hummi Wass, 1955 was discovered based on genetics and coloration.

12S↗

Weak phylogenetic and habitat effects on root trait variation of 218 Neotropical tree species

Tropical forests harbor a large diversity of closely related tree species that can thrive across habitats. This biodiversity has been found to correspond to large functional diversity in aboveground traits, and likely also relates to belowground trait variation. Globally, root trait (co-)variation is driven by different belowground resource strategies of species, environmental variation, and phylogeny; however, these patterns mostly reflect observations from temperate biomes and remain unconfirmed in tropical trees. We examine phylogenetic and environmental effects on root trait (co-)variation of trees across habitats in an Amazonian rainforest. Roots of 218 tree species from ten dominant families were sampled across three major habitats near Manaus, Brazil. We quantified five morphological and architectural root traits to (i) investigate how they reflected different resource strategies across species, (ii) compare them between families and superorders to test phylogenetic effects, and (iii) compare them between habitats to determine environmental effects on root trait expressions and variability. Root traits discriminated species along a tradeoff between root diameter and root branching and, secondly, due to variation in root tissue density. Our results further show weak phylogenetic effects on tropical tree root variation, for example, families from the same superorder showed large divergence in their root traits, while those from different superorders often overlapped in their root morphology and architecture. Root traits differed significantly between habitats but habitat type had only little effect on overall root trait variation. Our work suggests that the dimensions and drivers that underlie (co-)variation in tropical root traits may differ from global patterns defined by mostly temperate datasets. Due to (a)biotic environmental differences, different root trait dimensions may underlie the belowground functional diversity in (Neo)tropical forests, and we found little evidence for the strong phylogenetic conservatism observed in root traits in temperate biomes. We highlight important avenues for future research on tropical roots in order to determine the degree of, and shifts in functional diversity belowground as communities and environments change in tropical forests.

59 BASIC BIOLOGICAL SCIENCES↗

Novel symmetry-preserving neural network model for phylogenetic inference

Abstract Motivation Scientists world-wide are putting together massive efforts to understand how the biodiversity that we see on Earth evolved from single-cell organisms at the origin of life and this diversification process is represented through the Tree of Life. Low sampling rates and high heterogeneity in the rate of evolution across sites and lineages produce a phenomenon denoted “long branch attraction” (LBA) in which long nonsister lineages are estimated to be sisters regardless of their true evolutionary relationship. LBA has been a pervasive problem in phylogenetic inference affecting different types of methodologies from distance-based to likelihood-based. Results Here, we present a novel neural network model that outperforms standard phylogenetic methods and other neural network implementations under LBA settings. Furthermore, unlike existing neural network models in phylogenetics, our model naturally accounts for the tree isomorphisms via permutation invariant functions which ultimately result in lower memory and allows the seamless extension to larger trees. Availability and implementation We implement our novel theory on an open-source publicly available GitHub repository: https://github.com/crsl4/nn-phylogenetics.

59 BASIC BIOLOGICAL SCIENCES↗

Combining biomarker and virus phylogenetic models improves HIV-1 epidemiological source identification

To identify and stop active HIV transmission chains new epidemiological techniques are needed. Here, we describe the development of a multi-biomarker augmentation to phylogenetic inference of the underlying transmission history in a local population. HIV biomarkers are measurable biological quantities that have some relationship to the amount of time someone has been infected with HIV. To train our model, we used five biomarkers based on real data from serological assays, HIV sequence data, and target cell counts in longitudinally followed, untreated patients with known infection times. The biomarkers were modeled with a mixed effects framework to allow for patient specific variation and general trends, and fit to patient data using Markov Chain Monte Carlo (MCMC) methods. Subsequently, the density of the unobserved infection time conditional on observed biomarkers were obtained by integrating out the random effects from the model fit. This probabilistic information about infection times was incorporated into the likelihood function for the transmission history and phylogenetic tree reconstruction, informed by the HIV sequence data. To critically test our methodology, we developed a coalescent-based simulation framework that generates phylogenies and biomarkers given a specific or general transmission history. Testing on many epidemiological scenarios showed that biomarker augmented phylogenetics can reach 90% accuracy under idealized situations. Under realistic within-host HIV-1 evolution, involving substantial within-host diversification and frequent transmission of multiple lineages, the average accuracy was at about 50% in transmission clusters involving 5–50 hosts. Realistic biomarker data added on average 16 percentage points over using the phylogeny alone. Using more biomarkers improved the performance. Shorter temporal spacing between transmission events and increased transmission heterogeneity reduced reconstruction accuracy, but larger clusters were not harder to get right. More sequence data per infected host also improved accuracy. We show that the method is robust to incomplete sampling and that adding biomarkers improves reconstructions of real HIV-1 transmission histories. The technology presented here could allow for better prevention programs by providing data for locally informed and tailored strategies.

60 APPLIED LIFE SCIENCES↗

Genome sequence, phylogenetic analysis, and structure-based annotation reveal metabolic potential of Chlorella sp. SLA-04

Algae are a broad class of photosynthetic eukaryotes that are phylogenetically and physiologically diverse. Most of the phylogenetic diversity has been inferred from 18S rDNA sequencing since there are only a few complete genomes available in public databases. Here we use ultra-long-read Nanopore sequencing to determine a gapless, telomere-to-telomere complete genome sequence of Chlorella sp. SLA-04, previously described as Chlorella sorokiniana SLA-04. Chlorella sp. SLA-04 is a green alga that grows to high cell density in a wide variety of environments - high and neutral pH, high and low alkalinity, and high and low salinity. SLA-04's ability to grow in high pH and high alkalinity media without external CO 2 supply is favorable for large-scale algal biomass production. Phylogenetic analysis performed using ribosomal DNA and conserved protein sequences consistently reveal that Chlorella sp. SLA-04 forms a distinct lineage from other strains of Chlorella sorokiniana. We complement traditional genome annotation methods with high throughput structural predictions and demonstrate that this approach expands functional prediction of the SLA-04 proteome. Genomic analysis of the SLA-04 genome identifies the genes capable of utilizing TCA cycle intermediates to replenish cytosolic acetyl-CoA pools for lipid production. We also identify a complete metabolic pathway for sphingolipid anabolism that may allow SLA-04 to readily adapt to changing environmental conditions and facilitate robust cultivation in mass production systems. Altogether, this work clarifies the phylogeny of Chlorella sp. SLA-04 within Trebouxiophyceae and demonstrates how structural predictions can be used to improve annotation beyond sequencebased methods.

59 BASIC BIOLOGICAL SCIENCES↗

Phylogenetic distribution and experimental characterization of corrinoid production and dependence in soil bacterial isolates

Abstract Soil microbial communities impact carbon sequestration and release, biogeochemical cycling, and agricultural yields. These global effects rely on metabolic interactions that modulate community composition and function. However, the physicochemical and taxonomic complexity of soil and the scarcity of available isolates for phenotypic testing are significant barriers to studying soil microbial interactions. Corrinoids—the vitamin B12 family of cofactors—are critical for microbial metabolism, yet they are synthesized by only a subset of microbiome members. Here, we evaluated corrinoid production and dependence in soil bacteria as a model to investigate the ecological roles of microorganisms involved in metabolic interactions. We isolated and characterized a taxonomically diverse collection of 161 soil bacteria from a single study site. Most corrinoid-dependent bacteria in the collection prefer B12 over other corrinoids, while all tested producers synthesize B12, indicating metabolic compatibility between producers and dependents in the collection. Furthermore, a subset of producers release B12 at levels sufficient to support dependent isolates in laboratory culture at estimated ratios of up to 1000 dependents per producer. Within our isolate collection, we did not find strong phylogenetic patterns in corrinoid production or dependence. Upon investigating trends in the phylogenetic dispersion of corrinoid metabolism categories across sequenced bacteria from various environments, we found that these traits are conserved in 47 out of 85 genera. Together, these phenotypic and genomic results provide evidence for corrinoid-based metabolic interactions among bacteria and provide a framework for the study of nutrient-sharing ecological interactions in microbial communities.

59 BASIC BIOLOGICAL SCIENCES↗

Identifying impacts of contact tracing on HIV epidemiological inference from phylogenetic data

Abstract Robust sampling methods are foundational to inferences using phylogenies. Yet the impact of using contact tracing, a type of non-uniform sampling used in public health applications such as infectious disease outbreak investigations, has not been investigated in the molecular epidemiology field. To understand how contact tracing influences a recovered phylogeny, we developed a new simulation tool called SEEPS (Sequence Evolution and Epidemiological Process Simulator) that allows for the simulation of contact tracing and the resulting transmission tree, pathogen phylogeny, and corresponding virus genetic sequences. Importantly, SEEPS takes within-host evolution into account when generating pathogen phylogenies and sequences from transmission histories. Using SEEPS, we demonstrate that contact tracing can significantly impact the structure of the resulting tree, as described by popular tree statistics. Contact tracing generates phylogenies that are less balanced than the underlying transmission process, less representative of the larger epidemiological process, and affects the internal/external branch length ratios that characterize specific epidemiological scenarios. We also examined real data from a 2007–2008 Swedish HIV-1 outbreak and the broader 1998–2010 European HIV-1 epidemic to highlight the differences in contact tracing and expected phylogenies. Aided by SEEPS, we show that the data collection of the Swedish outbreak was strongly influenced by contact tracing even after downsampling, while the broader European Union epidemic showed little evidence of universal contact tracing, agreeing with the known epidemiological information about sampling and spread. Overall, our results highlight the importance of including possible non-uniform sampling schemes when examining phylogenetic trees. For that, SEEPS serves as a useful tool to evaluate such impacts, thereby facilitating better phylogenetic inferences of the characteristics of a disease outbreak. SEEPS is available at https://github.com/MolEvolEpid/SEEPS.

Virology↗