Single Nucleotide Sequence Variants Research Articles

BlastFrost: fast querying of 100,000s of bacterial genomes in Bifrost graphs

BlastFrost is a highly efficient method for querying 100,000s of genome assemblies, building on Bifrost, a dynamic data structure for compacted and colored de Bruijn graphs. BlastFrost queries a Bifrost data structure for sequences of interest and extracts local subgraphs, enabling the identification of the presence or absence of individual genes or single nucleotide sequence variants. We show two examples using Salmonella genomes: finding within minutes the presence of genes in the SPI-2 pathogenicity island in a collection of 926 genomes and identifying single nucleotide polymorphisms associated with fluoroquinolone resistance in three genes among 190,209 genomes. BlastFrost is available at https://github.com/nluhmann/BlastFrost/tree/master/data.

Influence of oxytocin receptor single nucleotide sequence variants on contractility of human myometrium: an in vitro functional study

BackgroundOxytocin receptor (OXTR) gene variants have been shown to affect the prevalence of preterm birth, mode of delivery and oxytocin (OXT) requirements for labor induction and augmentation. We hypothesized that this might be associated with different myometrium responses to oxytocin. Our aim was to investigate the influence of a selection of eight OXTR gene single nucleotide variants on oxytocin-induced stimulation of human myometrium contractility in vitro.MethodsHuman myometrium biopsies were collected during elective cesarean sections at term, if patients had given informed consent. Myometrial strips were submerged under tension in an organ bath and allowed to contract; the remaining material was stored at − 80 °C for further determination of relevant genetics and mRNA level. The area under the curve (AUC) of all contractions taking place in the absence of OXT and of those occurring upon OXT addition (for 30 min each) was measured. OXT stimulation, defined as the ratio between AUC measurements after OXT addition and those in the absence of OXT was calculated for each strip. TaqMan™ Assays were used to detect the allele distribution of the eight OXTR variants and to determine the relative amounts of OXTR-mRNA in the samples. For each variant, oxytocin stimulation of contractility was compared between samples homozygous for the reference allele (reference group) and samples with at least one variant allele (variant group) by linear regression.ResultsSixty samples were included in the present study. For rs1042778, rs11706648, rs4686301, rs53576, rs237895, and rs237902, OXT stimulation was similar in the reference and in the variant groups. However, the values of OXT stimulation differed significantly between the reference and the variant groups for rs4686302 (3.1 vs. 4.1 times; p = 0.022) and rs237888 (3.2 vs. 5.5 times; p = 0.001). No significant differences between the levels of OXTR-mRNA in the various reference and corresponding variant groups were detected.ConclusionsPatients with variant alleles of rs237888 and/or rs4686302 may be more sensitive to oxytocin stimulation, explaining why these sequence variants have been associated with lower cesarean section prevalence and premature birth, respectively.

Genetic influences on the human oral microbiome

BackgroundThe human oral microbiome is formed early in development. Its composition is influenced by environmental factors including diet, substance use, oral health, and overall health and disease. The influence of human genes on the composition and stability of the oral microbiome is still poorly understood. We studied both environmental and genetic characteristics on the oral microbiome in a large twin sample as well as in a large cohort of unrelated individuals. We identify several significantly heritable features of the oral microbiome. The heritability persists in twins even when their cohabitation changes. The heritability of these traits correlates with the cumulative genetic contributions of over half a million single nucleotide sequence variants measured in a different population of unrelated individuals. Comparison of same-sex and opposite sex cotwins showed no significant differences. We show that two new loci on chromosomes 7 and 12 are associated with the most heritable traits.ResultsAn analysis of 752 twin pairs from the Colorado Twin Registry, shows that the beta-diversity of monozygotic twins is significantly lower than for dizygotic or unrelated individuals. This is independent of cohabitation status. Intraclass correlation coefficients of nearly all taxa examined were higher for MZ than DZ twin pairs. A comparison of individuals sampled over 2-7 years confirmed previous reports that the oral microbiome remains relatively more stable in individuals over that time than to unrelated people. Twin modeling shows that a number of microbiome phenotypes were more than 50% heritable consistent with the hypothesis that human genes influence microbial populations. To identify loci that could influence microbiome phenotypes, we carried out an unbiased GWAS analysis which identified one locus on chromosome 7 near the gene IMMPL2 that reached genome-wide significance after correcting for multiple testing. Another locus on chromosome 12 near the non-coding RNA gene INHBA-AS1 achieved genome-wide significance when analyzed using KGG4 that sums SNP significance across coding genes.DiscussionUsing multiple methods, we have demonstrated that some aspects of the human oral microbiome are heritable and that with a relatively small sample we were able to identify two previously unidentified loci that may be involved.

Napier grass stunt disease prevalence, incidence, severity and genetic variability of the associated phytoplasma in Uganda

The prevalence, incidence and severity of Napier grass stunt disease (NGSD) caused by phytoplasma on Pennisetum purpureum, the main fodder for livestock under intensive and semi-intensive management systems in Uganda were determined following a field survey carried out in 17 districts. A total of 298 Napier grass fields were visited and NGSD status visually assessed and 1192 samples collected for identification and confirmation of the phytoplasma by polymerase chain reaction (PCR) assays using universal primers P1/P6 nested with R16F2n/R16R2n and, tuf primers 890/340 nested with 835 and 400. From these, 221 PCR products were sequenced and sequences aligned. Napier grass stunt disease is widely spread at an epidemic proportion, with the districts at different risk levels. The most affected districts are in central, East and North parts of the country while those in the west are least affected. Sequence alignments and Blast searches showed that the phytoplasma causing NGSD in Uganda belonged to the phytoplasma group 16SrXI, with single nucleotide sequence variants in a few districts. Therefore, there is a need for development of an area wide NGSD management strategy to contain the disease.

An integrated model of the transcriptome of HER2-positive breast cancer.

Our goal in these analyses was to use genomic features from a test set of primary breast tumors to build an integrated transcriptome landscape model that makes relevant hypothetical predictions about the biological and/or clinical behavior of HER2-positive breast cancer. We interrogated RNA-Seq data from benign breast lesions, ER+, triple negative, and HER2-positive tumors to identify 685 differentially expressed genes, 102 alternatively spliced genes, and 303 genes that expressed single nucleotide sequence variants (eSNVs) that were associated with the HER2-positive tumors in our survey panel. These features were integrated into a transcriptome landscape model that identified 12 highly interconnected genomic modules, each of which represents a cellular processes pathway that appears to define the genomic architecture of the HER2-positive tumors in our test set. The generality of the model was confirmed by the observation that several key pathways were enriched in HER2-positive TCGA breast tumors. The ability of this model to make relevant predictions about the biology of breast cancer cells was established by the observation that integrin signaling was linked to lapatinib sensitivity in vitro and strongly associated with risk of relapse in the NCCTG N9831 adjuvant trastuzumab clinical trial dataset. Additional modules from the HER2 transcriptome model, including ubiquitin-mediated proteolysis, TGF-beta signaling, RHO-family GTPase signaling, and M-phase progression, were linked to response to lapatinib and paclitaxel in vitro and/or risk of relapse in the N9831 dataset. These data indicate that an integrated transcriptome landscape model derived from a test set of HER2-positive breast tumors has potential for predicting outcome and for identifying novel potential therapeutic strategies for this breast cancer subtype.

High-Throughput Detection of Actionable Genomic Alterations in Clinical Tumor Samples by Targeted, Massively Parallel Sequencing

Knowledge of "actionable" somatic genomic alterations present in each tumor (e.g., point mutations, small insertions/deletions, and copy-number alterations that direct therapeutic options) should facilitate individualized approaches to cancer treatment. However, clinical implementation of systematic genomic profiling has rarely been achieved beyond limited numbers of oncogene point mutations. To address this challenge, we utilized a targeted, massively parallel sequencing approach to detect tumor genomic alterations in formalin-fixed, paraffin-embedded (FFPE) tumor samples. Nearly 400-fold mean sequence coverage was achieved, and single-nucleotide sequence variants, small insertions/deletions, and chromosomal copynumber alterations were detected simultaneously with high accuracy compared with other methods in clinical use. Putatively actionable genomic alterations, including those that predict sensitivity or resistance to established and experimental therapies, were detected in each tumor sample tested. Thus, targeted deep sequencing of clinical tumor material may enable mutation-driven clinical trials and, ultimately, "personalized" cancer treatment. Despite the rapid proliferation of targeted therapeutic agents, systematic methods to profile clinically relevant tumor genomic alterations remain underdeveloped. We describe a sequencingbased approach to identifying genomic alterations in FFPE tumor samples. These studies affirm the feasibility and clinical utility of targeted sequencing in the oncology arena and provide a foundation for genomics-based stratification of cancer patients.

Inferring relative proportions of DNA variants from sequencing electropherograms

Determination of the relative copy number of single-nucleotide sequence variants (SNVs) within a DNA sample is a frequent experimental goal. Various methods can be applied to this problem, although hybridization-based approaches tend to suffer from high-setup cost and poor adaptability, while others (such as pyrosequencing) may not be accessible to all laboratories. The potential to extract relative copy number information from standard dye-terminator electropherograms has been little explored, yet this technology is cheap and widely accessible. Since several biologically important loci have paralogous copies that interfere with genotyping, and which may also display copy number variation (CNV), there are many situations in which determination of the relative copy number of SNVs is desirable. We have developed a desktop application, QSVanalyzer, which allows high-throughput quantification of the proportions of DNA sequences containing SNVs. In reconstruction experiments, QSVanalyzer accurately estimated the known relative proportions of SNVs. By analyzing a large panel of genomic DNA samples, we demonstrate the ability of the software to analyze not only common biallelic SNVs, but also SNVs within a locus at which gene conversion between four genomic paralogs operates, and within another that is subject to CNV. QSVanalyzer is freely available at http://dna.leeds.ac.uk/qsv/. It requires the Microsoft .NET framework version 2.0, which can be installed on all Microsoft operating systems from Windows 98 onwards. msjimc@leeds.ac.uk Supplementary data are available at Bioinformatics online.

Non-synonymous sequence variants within the oxygen-dependent degradation domain of the HIF1Agene are not associated with pre-eclampsia in the Finnish population

BackgroundReduced placental perfusion predisposes to the maternal syndrome pre-eclampsia characterized by systemically reduced perfusion. Considerable data support the role of angiogenic factors in the development of the maternal syndrome. Hypoxia-inducible factor (HIF-1) mediates the cellular responses to hypoxia e.g. by promoting angiogenesis.MethodsHere we studied whether two single nucleotide sequence variants, c.1744 C>T that changes residue 582 of HIF-1α from proline to serine (P582S) and c.1762 G>A that changes residue 588 of HIF-1α from alanine to threonine (A588T) in the exon 12 of the HIF1A gene, are associated with pre-eclampsia. We studied 108 women with pre-eclampsia in their first pregnancy, and 101 controls with normotensive pregnancies. Pre-eclampsia was defined as a blood pressure level of at least 140/90 mmHg in a woman who was normotensive before 20 weeks of gestation, and proteinuria at least of 0.3 g per 24-hour urine collection. The patients and controls were genotyped for variations in the exon 12 of HIF1A gene by sequencingResultsThe frequencies of the c.1744 C>T and c.1762G>A sequence variants were not significantly different between women with pre-eclamptic first pregnancies and women with normotensive pregnancies. In addition, two synonymous variants (c.1740G>A and c.1800A>T) were detected at comparable levels in the two groups. All variants were identified in the heterozygous form.ConclusionThe sequence variants in the exon 12 of the HIF1A gene were not associated with pre-eclampsia in the Finnish population.

Genetic variation in human P450 oxidoreductase

Catalysis by all 50 Type II (microsomal) P450 enzymes, including steroidogenic P450c17, P450c21, and P450aro and hepatic drug-metabolizing enzymes requires electron donation from P450 oxidoreductase (POR). POR knockout mice are embryonic lethal, but human POR mutations cause a complex disorder of steroidogenesis. Disorders of hepatic drug metabolism in human POR deficiency have not yet been described. To understand the potential contribution of POR to pharmacogenetics, we sequenced the POR gene in 842 normal persons from 4 ethnic groups. We detected 140 single nucleotide sequence variants of which 43 were in >1% of alleles, including 15 missense mutants; this brings the total of known POR missense mutants to 35. A503V was found on 28% of alleles, varying from 19% in African Americans to 37% in Chinese Americans. We expressed all 35 missense mutants in E. coli and assayed their activities to: oxidize NADPH, reduce cytochrome c, support the 17alpha-hydroxylase and 17,20 lyase activities of bacterially expressed human P450c17, and support the metabolism of fluorogenic EOMCC by bacterially expressed human CYP1A2 and CYP2C19. These data show that there are great differences in the activities of some POR mutants depending on the electron recipient assayed; for example, Q153R causes severely impaired steroid biosynthesis in human patients and in vitro, but is a gain-of-function mutant with CYP1A2 and 2C19. A503V reduces both activities of P450c17 in half, but had no effect on CYP1A2 or 2C19. POR variants are a previously unappreciated source of genetic variation in patterns of steroid synthesis and drug metabolism.

Sequence analysis of 21 genes located in the Kartagener syndrome linkage region on chromosome 15q

Primary ciliary dyskinesia (PCD) is a rare genetic disorder, which shows extensive genetic heterogeneity and is mostly inherited in an autosomal recessive fashion. There are four genes with a proven pathogenetic role in PCD. DNAH5 and DNAI1 are involved in 28 and 10% of PCD cases, respectively, while two other genes, DNAH11 and TXNDC3, have been identified as causal in one PCD family each. We have previously identified a 3.5 cM (2.82 Mb) region on chromosome 15q linked to Kartagener syndrome (KS), a subtype of PCD characterized by the randomization of body organ positioning. We have now refined the KS candidate region to a 1.8 Mb segment containing 18 known genes. The coding regions of these genes and three neighboring genes were subjected to sequence analysis in seven KS probands, and we were able to identify 60 single nucleotide sequence variants, 35 of which resided in mRNA coding sequences. However, none of the variations alone could explain the occurrence of the disease in these patients.

Human dopamine transporter gene: coding region conservation among normal, Tourette's disorder, alcohol dependence and attention-deficit hyperactivity disorder populations.

The dopamine transporter (DAT) provides major regulation of the synaptic levels of dopamine and is a principal target of psychostimulant drugs. Associations between DAT gene polymorphisms and human disorders with possible links to dopaminergic neurotransmission, including attention-deficit/hyperactivity disorder (ADHD) and consequences of cocaine and alcohol administration, have been reported. We now report approximately 60000 bp of genomic sequence containing the entire DAT gene. This sequence was used to amplify each of the 15 DAT gene exons and several introns and analyze these amplification products by single-stranded sequence conformation (SSCP) and/or direct sequencing. These results define silent allelic single nucleotide sequence variants in DAT gene exons 2, 6, 9 and 15. Rare conservative mutations are identified in amino acids encoded by DAT exons 2 and 8. Analyses of the common nucleotide variants and the previously reported VNTR in the non-coding region of exon 15 define the pattern of linkage disequilibrium across the DAT locus. These comprehensive analyses, however, fail to identify any common protein coding DAT sequence variant in more than 150 unrelated individuals free of neuropsychiatric disease, 109 individuals meeting City of Hope criteria for Tourette's syndrome, 64 individuals with DSM-IV diagnoses of ethanol dependence, or 15 individuals with ADHD. These data are consistent with substantial evolutionary conservation of the DAT protein sequence. They suggest that gene variants that alter levels of DAT expression provide the best current candidate mechanism for reported associations between DAT gene markers, ADHD and other more tentatively associated neuropsychiatric disorders.