Dataset Information

Canonical, stable, general mapping using context schemes.

ABSTRACT:

Motivation

Sequence mapping is the cornerstone of modern genomics. However, most existing sequence mapping algorithms are insufficiently general.

Results

We introduce context schemes: a method that allows the unambiguous recognition of a reference base in a query sequence by testing the query for substrings from an algorithmically defined set. Context schemes only map when there is a unique best mapping, and define this criterion uniformly for all reference bases. Mappings under context schemes can also be made stable, so that extension of the query string (e.g. by increasing read length) will not alter the mapping of previously mapped positions. Context schemes are general in several senses. They natively support the detection of arbitrary complex, novel rearrangements relative to the reference. They can scale over orders of magnitude in query sequence length. Finally, they are trivially extensible to more complex reference structures, such as graphs, that incorporate additional variation. We demonstrate empirically the existence of high-performance context schemes, and present efficient context scheme mapping algorithms.

Availability and implementation

The software test framework created for this study is available from https://registry.hub.docker.com/u/adamnovak/sequence-graphs/.

Contact

anovak@soe.ucsc.edu

Supplementary information

Supplementary data are available at Bioinformatics online.

SUBMITTER: Novak AM

PROVIDER: S-EPMC4757953 | biostudies-literature |

REPOSITORIES: biostudies-literature

ACCESS DATA

Similar Datasets

Project description:BackgroundMyocardial T1 mapping shows promise for assessment of cardiomyopathies. Most myocardial T1 mapping techniques, such as modified Look-Locker inversion recovery (MOLLI), generate one T1 map per breath-held acquisition (9-17 heartbeats), which prolongs multislice protocols and may be unsuitable for patients with breath-holding difficulties.PurposeTo develop and characterize novel shortened inversion recovery based T1 mapping schemes of 2-5 heartbeats.Study typeProspective.Population/phantomNumerical simulations, agarose/NiCl2 phantom, 16 healthy volunteers, and 24 patients.Field strength/sequence1.5T/MOLLI.AssessmentAll shortened T1 mapping schemes were characterized and compared with a conventional MOLLI scheme (5-(3)-3) in terms of accuracy, precision, spatial variability, and repeatability.Statistical testsKruskal-Wallis, Wilcoxon rank sum tests, analysis of variance, Student's t-tests, Bland-Altman analysis, and Pearson correlation analysis.ResultsAll shortened schemes provided limited T1 time variations (≤2% for T1 times ≤1200 msec) and limited penalty of precision (by a factor of ~1.4-1.5) when compared with MOLLI in numerical simulations. In phantom, differences between all schemes in terms of accuracy, spatial variability, and repeatability did not reach statistical significance (P > 0.71). In healthy volunteers, there were no statistically significant differences between all schemes in terms of native T1 times and repeatability for myocardium (P = 0.21 and P = 0.87, respectively) and blood (P = 0.79 and P = 0.41, respectively). All shortened schemes led to a limited increase of spatial variability for native myocardial T1 mapping with respect to MOLLI (by a factor of 1.2) (P < 0.0001). In both healthy volunteers and patients, the two-heartbeat scheme and MOLLI led to highly linearly correlated T1 times (correlation coefficients ≥0.83).Data conclusionThe proposed two-heartbeat T1 mapping scheme yields a 5-fold acceleration compared with MOLLI, with highly linearly correlated T1 times, no significant difference of repeatability, and limited spatial variability penalty at 1.5T. This approach may enable myocardial T1 mapping in patients with severe breath-holding difficulties and reduce the examination time of multislice protocols.Level of evidence1 Technical Efficacy Stage: 3 J. Magn. Reson. Imaging 2019;50:641-654.

Project description:Simultaneous mapping of multiple behavioral domains into brain networks remains a major challenge. Here, we shed some light on this problem by employing a combination of machine learning, structural and functional brain networks at different spatial resolutions (also known as scales), together with performance scores across multiple neurobehavioral domains, including sensation, motor skills, and cognition. Provided by the Human Connectome Project, we make use of three cohorts: 640 participants for model training, 160 subjects for validation, and 200 subjects for model performance testing thus enhancing prediction generalization. Our modeling consists of two main stages, namely dimensionality reduction in brain network features at multiple scales, followed by canonical correlation analysis, which determines an optimal linear combination of connectivity features to predict multiple behavioral performance scores. To assess the differences in the predictive power of each modality, we separately applied three different strategies: structural unimodal, functional unimodal, and multimodal, that is, structural in combination with functional features of the brain network. Our results show that the multimodal association outperforms any of the unimodal analyses. Then, to answer which human brain structures were most involved in predicting multiple behavioral scores, we simulated different synthetic scenarios in which in each case we completely deleted a brain structure or a complete resting state network, and recalculated performance in its absence. In deletions, we found critical structures to affect performance when predicting single behavioral domains, but this occurred in a lesser manner for prediction of multi-domain behavior. Overall, our results confirm that although there are synergistic contributions between brain structure and function that enhance behavioral prediction, brain networks may also be mutually redundant in predicting multidomain behavior, such that even after deletion of a structure, the connectivity of the others can compensate for its lack in predicting behavior.

Project description:A core component of the M-NM-1-proteobacterial general stress response is the extracytoplasmic function (ECF) sigma factor EcfG, exclusively present in this taxonomic class. Half of the completed M-NM-1-proteobacterial genome sequences contain two or more copies of genes encoding M-OM-^CEcfG-like sigma factors, with the primary copy typically located adjacent to genes coding for a cognate anti-sigma factor (NepR) and two-component response regulator (PhyR). So far, the widespread occurrence of additional, non-canonical M-OM-^CEcfG copies has not satisfactorily been explained. This work explores the hierarchical relation between Rhizobium etli M-OM-^CEcfG1 and M-OM-^CEcfG2, canonical and non-canonical M-OM-^CEcfG proteins, respectively. Contrary to reports in other species, we find that M-OM-^CEcfG1 and M-OM-^CEcfG2 act in parallel, as nodes of a complex regulatory network, rather than in series, as elements of a linear regulatory cascade. We demonstrate that both sigma factors control unique yet also shared target genes, corroborating phenotypic evidence. M-OM-^CEcfG1 drives expression of rpoH2, explaining the increased heat sensitivity of an ecfG1 mutant, while katG is under control of M-OM-^CEcfG2, accounting for reduced oxidative stress resistance of an ecfG2 mutant. We also identify non-coding RNA genes as novel M-OM-^CEcfG targets. We propose a modified model for general stress response regulation in R. etli, in which M-OM-^CEcfG1 and M-OM-^CEcfG2 function largely independently. Based on a phylogenetic analysis and considering the prevalence of M-NM-1-proteobacterial genomes with multiple M-OM-^CEcfG copies, this model may also be applicable to numerous other species. In this study, we investigate to which degree R. etli EcfG1 and EcfG2 regulation is interdependent. Furthermore, we demonstrate that both sigma factors control distinct regulons and pinpoint unique EcfG2 target genes. We also identify non-coding RNAs (ncRNAs) as novel targets of EcfG1 and EcfG2 and show that expression of at least one of these ncRNAs is under direct EcfG control. The presence of ncRNAs in the regulon of an EcfG type sigma factor has not been described previously and adds an extra level of complexity to the regulation of the general stress response in Alpha-proteobacteria. Moreover, considering the widespread existence of Alpha-proteobacterial genomes with multiple EcfG copies, the here presented results not only lead to a better understanding of the general stress response in R. etli, but also contribute to a conceptual paradigm for general stress response regulation by multiple EcfG proteins in Alpha-proteobacteria. Comparative transcriptome analyses were carried out on four samples: the parental strain (WT) and mutants in either sigma factor gene ecfG1, sigma factor gene ecfG2, and both genes combined.

Project description:BackgroundChildhood obesity is high on the policy agenda of wealthier nations, and many interventions have been developed to address it. This work describes an overview of schemes for obese and overweight children and young people in England, and the 'mapping' approach we used.MethodsOur search strategy, inclusion criteria and coding frame had to be suitable for describing a potentially large number of schemes within a short timeframe. Data were collected from key informants, scheme publicity and reports, and via a web-survey. To be included, schemes had to be based in England, follow a structured programme lasting at least two weeks, promote healthy weight, and be delivered exclusively to overweight and/or obese children and young people (age range 4-18). Data were entered into a coding frame recording similar information for each scheme, including any underpinning research evidence, evaluation or monitoring reports. Priority questions were identified in consultation with colleagues from the Department of Health and the Cross Government Obesity Unit.ResultsFifty-one schemes were identified. Some operated in multiple areas, and by using estimates of the number of schemes provided by multi-site scheme leads, we found that between 314 and 375 local programmes were running at any time. Uncertainty is largely due to the largest scheme provider undergoing rapid expansion at the time of the mapping exercise and therefore able to provide only an estimate of the number of programmes running. Many schemes were similar in their approach, had been recently established and were following NICE guidelines on interventions to promote healthy weight. Rigorous evaluation was rare.ConclusionsOur methods enabled us to produce a rapid overview of service activity across a wide geographic area and a range of organisations and sectors. In order to develop the evidence base for childhood obesity interventions, rigorous evaluation of these schemes is required. This overview can serve as a starting point for evaluations of interventions to address obesity. More generally, a rapid and systematic approach of this type is transferable to other types of service activity in health and social care, and may be a tool to inform public health planning.

Project description:A core component of the α-proteobacterial general stress response is the extracytoplasmic function (ECF) sigma factor EcfG, exclusively present in this taxonomic class. Half of the completed α-proteobacterial genome sequences contain two or more copies of genes encoding σEcfG-like sigma factors, with the primary copy typically located adjacent to genes coding for a cognate anti-sigma factor (NepR) and two-component response regulator (PhyR). So far, the widespread occurrence of additional, non-canonical σEcfG copies has not satisfactorily been explained. This work explores the hierarchical relation between Rhizobium etli σEcfG1 and σEcfG2, canonical and non-canonical σEcfG proteins, respectively. Contrary to reports in other species, we find that σEcfG1 and σEcfG2 act in parallel, as nodes of a complex regulatory network, rather than in series, as elements of a linear regulatory cascade. We demonstrate that both sigma factors control unique yet also shared target genes, corroborating phenotypic evidence. σEcfG1 drives expression of rpoH2, explaining the increased heat sensitivity of an ecfG1 mutant, while katG is under control of σEcfG2, accounting for reduced oxidative stress resistance of an ecfG2 mutant. We also identify non-coding RNA genes as novel σEcfG targets. We propose a modified model for general stress response regulation in R. etli, in which σEcfG1 and σEcfG2 function largely independently. Based on a phylogenetic analysis and considering the prevalence of α-proteobacterial genomes with multiple σEcfG copies, this model may also be applicable to numerous other species. In this study, we investigate to which degree R. etli EcfG1 and EcfG2 regulation is interdependent. Furthermore, we demonstrate that both sigma factors control distinct regulons and pinpoint unique EcfG2 target genes. We also identify non-coding RNAs (ncRNAs) as novel targets of EcfG1 and EcfG2 and show that expression of at least one of these ncRNAs is under direct EcfG control. The presence of ncRNAs in the regulon of an EcfG type sigma factor has not been described previously and adds an extra level of complexity to the regulation of the general stress response in Alpha-proteobacteria. Moreover, considering the widespread existence of Alpha-proteobacterial genomes with multiple EcfG copies, the here presented results not only lead to a better understanding of the general stress response in R. etli, but also contribute to a conceptual paradigm for general stress response regulation by multiple EcfG proteins in Alpha-proteobacteria.