Dataset Information

DISSEQT-DIStribution-based modeling of SEQuence space Time dynamics.

ABSTRACT: Rapidly evolving microbes are a challenge to model because of the volatile, complex, and dynamic nature of their populations. We developed the DISSEQT pipeline (DIStribution-based SEQuence space Time dynamics) for analyzing, visualizing, and predicting the evolution of heterogeneous biological populations in multidimensional genetic space, suited for population-based modeling of deep sequencing and high-throughput data. The pipeline is openly available on GitHub (https://github.com/rasmushenningsson/DISSEQT.jl, accessed 23 June 2019) and Synapse (https://www.synapse.org/#!Synapse: syn11425758, accessed 23 June 2019), covering the entire workflow from read alignment to visualization of results. Our pipeline is centered around robust dimension and model reduction algorithms for analysis of genotypic data with additional capabilities for including phenotypic features to explore dynamic genotype-phenotype maps. We illustrate its utility and capacity with examples from evolving RNA virus populations, which present one of the highest degrees of genetic heterogeneity within a given population found in nature. Using our pipeline, we empirically reconstruct the evolutionary trajectories of evolving populations in sequence space and genotype-phenotype fitness landscapes. We show that while sequence space is vastly multidimensional, the relevant genetic space of evolving microbial populations is of intrinsically low dimension. In addition, evolutionary trajectories of these populations can be faithfully monitored to identify the key minority genotypes contributing most to evolution. Finally, we show that empirical fitness landscapes, when reconstructed to include minority variants, can predict phenotype from genotype with high accuracy.

SUBMITTER: Henningsson R

PROVIDER: S-EPMC6680062 | biostudies-literature | 2019 Jul

REPOSITORIES: biostudies-literature

ACCESS DATA

Publications

DISSEQT-DIStribution-based modeling of SEQuence space Time dynamics.

Henningsson R R Moratorio G G Bordería A V AV Vignuzzi M M Fontes M M

Virus evolution 20190701 2

Rapidly evolving microbes are a challenge to model because of the volatile, complex, and dynamic nature of their populations. We developed the DISSEQT pipeline (DIStribution-based SEQuence space Time dynamics) for analyzing, visualizing, and predicting the evolution of heterogeneous biological populations in multidimensional genetic space, suited for population-based modeling of deep sequencing and high-throughput data. The pipeline is openly available on GitHub (https://github.com/rasmushenning ...[more]

PMID: 31392032

Similar Datasets

Project description:The wheat curl mite, Aceria tosichella Keifer, one of the most destructive arthropod pests of bread wheat worldwide, inflicts significant annual reductions in grain yields. Moreover, A. tosichella is the only vector for several economically important wheat viruses in the Americas, Australia and Europe. To date, mite-resistant wheat genotypes have proven to be one of the most effective methods of controlling the A. tosichella-virus complex. Thus, it is important to elucidate A. tosichella population genetic structure, in order to better predict improved mite and virus management. Two genetically distinct A. tosichella lineages occur as pests of wheat in Australia, Europe, North America, South America and the Middle East. These lineages are known as type 1 and type 2 in Australia and North America and in Europe and South America as MT-8 and MT-1, respectively. Type 1 and type 2 mites in Australia and North America are delineated by internal transcribed spacer 1 region (ITS1) and cytochrome oxidase I region (COI) sequence differences. In North America, two A. tosichella genotypes known as biotypes are recognized by their response to the Cmc3 mite resistance gene in wheat. Aceria tosichella biotype 1 is susceptible to Cmc3 and biotype 2 is virulent to Cmc3. In this study, ITS1 and COI sequence differences in 25 different populations of A. tosichella of known biotype 1 or biotype 2 composition were characterized for ITS1 and COI sequence differences and used to model spatio-temporal dynamics based on biotype prevalence. Results showed that the proportion of biotype 1 and 2 varies both spatially and temporally. Greater ranges of cropland and grassland within 5000m of the sample site, as well as higher mean monthly precipitation during the month prior to sampling appear to reduce the probability of occurrence of biotype 1 and increase the probability of occurrence of biotype 2. The results suggest that spatio-temporal modeling can effectively improve A. tosichella management. Continual integration of additional current and future precipitation and ground cover data into the existing model will further improve the accuracy of predicting the occurrence of A. tosichella in annual wheat crops, allowing producers to make informed decisions about the selection of varieties with different A. tosichella resistance genes.

Project description:Coastal ecosystems experience a wide range of stressors including wave forces, storm surge, sea-level rise, and anthropogenic modification and are thus vulnerable to erosion. Urban coastal ecosystems are especially important due to the large populations these limited ecosystems serve. However, few studies have addressed the issue of urban coastal vulnerability at the landscape scale with spatial data that are finely resolved. The purpose of this study was to model and map coastal vulnerability and the role of natural habitats in reducing vulnerability in Jamaica Bay, New York, in terms of nine coastal vulnerability metrics (relief, wave exposure, geomorphology, natural habitats, exposure, exposure with no habitat, habitat role, erodible shoreline, and surge) under past (1609), current (2015), and future (2080) scenarios using InVEST 3.2.0. We analyzed vulnerability results both spatially and across all time periods, by stakeholder (ownership) and by distance to damage from Hurricane Sandy. We found significant differences in vulnerability metrics between past, current and future scenarios for all nine metrics except relief and wave exposure. The marsh islands in the center of the bay are currently vulnerable. In the future, these islands will likely be inundated, placing additional areas of the shoreline increasingly at risk. Significant differences in vulnerability exist between stakeholders; the Breezy Point Cooperative and Gateway National Recreation Area had the largest erodible shoreline segments. Significant correlations exist for all vulnerability (exposure/surge) and storm damage combinations except for exposure and distance to artificial debris. Coastal protective features, ranging from storm surge barriers and levees to natural features (e.g. wetlands), have been promoted to decrease future flood risk to communities in coastal areas around the world. Our methods of combining coastal vulnerability results with additional data and across multiple time periods have considerable potential to provide valuable predictions that resource managers can effectively use to identify areas for restoration and protection.

Project description:A new method has been developed to compute the probability that each amino acid in a protein sequence is in a particular secondary structural element. Each of these probabilities is computed using the entire sequence and a set of predefined structural class models. This set of structural classes is patterned after Jane Richardson's taxonomy for the domains of globular proteins. For each structural class considered, a mathematical model is constructed to represent constraints on the pattern of secondary structural elements characteristic of that class. These are stochastic models having discrete state spaces (referred to as hidden Markov models by researchers in signal processing and automatic speech recognition). Each model is a mathematical generator of amino acid sequences; the sequence under consideration is modeled as having been generated by one model in the set of candidates. The probability that each model generated the given sequence is computed using a filtering algorithm. The protein is then classified as belonging to the structural class having the most probable model. The secondary structure of the sequence is then analyzed using a "smoothing" algorithm that is optimal for that structural class model. For each residue position in the sequence, the smoother computes the probability that the residue is contained within each of the defined secondary structural elements of the model. This method has two important advantages: (1) the probability of each residue being in each of the modeled secondary structural elements is computed using the totality of the amino acid sequence, and (2) these probabilities are consistent with prior knowledge of realizable domain folds as encoded in each model. As an example of the method's utility, we present its application to flavodoxin, a prototypical alpha/beta protein having a central beta-sheet, and to thioredoxin, which belongs to a similar structural class but shares no significant sequence similarity.

Dataset Information

DISSEQT-DIStribution-based modeling of SEQuence space Time dynamics.

Publications

DISSEQT-DIStribution-based modeling of SEQuence space Time dynamics.

Similar Datasets

OmicsDI is part of the ELIXIR infrastructure

Tweets