From SNPs to pathways: integration of functional effect of sequence variations on models of cell signalling pathways

Background: Single nucleotide polymorphisms (SNPs) are the most frequent type of sequence variation between individuals, and represent a promising tool for finding genetic determinants of complex diseases and understanding the differences in drug response. In this regard, it is of particular interest to study the effect of non-synonymous SNPs in the context of biological networks such as cell signalling pathways. UniProt provides curated information about the functional and phenotypic effects of sequence variation, including SNPs, as well as on mutations of protein sequences. However, no strategy has been developed to integrate this information with biological networks, with the ultimate goal of studying the impact of the functional effect of SNPs in the structure and dynamics of biological networks. Results: First, we identified the different challenges posed by the integration of the phenotypic effect of sequence variants and mutations with biological networks. Second, we developed a strategy for the combination of data extracted from public resources, such as UniProt, NCBI dbSNP, Reactome and BioModels. We generated attribute files containing phenotypic and genotypic annotations to the nodes of biological networks, which can be imported into network visualization tools such as Cytoscape. These resources allow the mapping and visualization of mutations and natural variations of human proteins and their phenotypic effect on biological networks (e.g. signalling pathways, protein-protein interaction networks, dynamic models). Finally, an example on the use of the sequence variation data in the dynamics of a network model is presented. Conclusion: In this paper we present a general strategy for the integration of pathway and sequence variation data for visualization, analysis and modelling purposes, including the study of the functional impact of protein sequence variations on the dynamics of signalling pathways. This is of particular interest when the SNP or mutation is known to be associated to disease. We expect that this approach will help in the study of the functional impact of disease-associated SNPs on the behaviour of cell signalling pathways, which ultimately will lead to a better understanding of the mechanisms underlying complex diseases.

Sequence variations in the von Hippel-Lindau tumor suppressor gene in patients with intracranial aneurysms

The rupture of intracranial aneurysms leads to subarachnoid hemorrhage, which is often associated with poor outcome. Preventive treatment of unruptured intracranial aneurysms is possible and recommended. However, the lack of candidate genes precludes identifying patients at risk by genetic analyses. We observed intracranial aneurysms in 2 patients with von Hippel-Lindau (VHL) disease and the known disease-causing mutation c.292T > C (p.Tyr98His) in the VHL tumor suppressor gene. This study investigates whether the VHL gene is a possible candidate gene for aneurysm formation.

8p23 beta-defensin copy number determination by single-locus pseudogene-based paralog ratio tests risk bias due to low-frequency sequence variations.

BACKGROUND The copy number variation (CNV) in beta-defensin genes (DEFB) on human chromosome 8p23 has been proposed to contribute to the phenotypic differences in inflammatory diseases. However, determination of exact DEFB CN is a major challenge in association studies. Quantitative real-time PCR (qPCR), paralog ratio tests (PRT) and multiplex ligation-dependent probe amplification (MLPA) have been extensively used to determine DEFB CN in different laboratories, but inter-method inconsistencies were observed frequently. In this study we asked which one is superior among the three methods for DEFB CN determination. RESULTS We developed a clustering approach for MLPA and PRT to statistically correlate data from a single experiment. Then we compared qPCR, a newly designed PRT and MLPA for DEFB CN determination in 285 DNA samples. We found MLPA had the best convergence and clustering results of the raw data and the highest call rate. In addition, the concordance rates between MLPA or PRT and qPCR (32.12% and 37.99%, respectively) were unacceptably low with underestimated CN by qPCR. Concordance rate between MLPA and PRT (90.52%) was high but PRT systematically underestimated CN by one in a subset of samples. In these samples a sequence variant which caused complete PCR dropout of the respective DEFB cluster copies was found in one primer binding site of one of the targeted paralogous pseudogenes. CONCLUSION MLPA is superior to PRT and even more to qPCR for DEFB CN determination. Although the applied PRT provides in most cases reliable results, such a test is particularly sensitive to low-frequency sequence variations preferably accumulating in loci like pseudogenes which are most likely not under selective pressure. In the light of the superior performance of multiplex assays, the drawbacks of such single PRTs could be overcome by combining more test markers.

Sequence variation of the alpha-lactalbumin gene in Holstein and Nellore cows

The alpha-lactalbumin is a subunit of lactose-synthase, an enzyme responsible for lactose production, a disaccharide that influences milk production. Sequence variations of bovine alpha-lactalbumin have been associated with differences in milk yield. This study aimed to analyze allelic frequency differences at position-1689 (g. AG) and+15 (g. AG) of the alpha-lactalbumin gene in Holstein (Bos taurus) and Nellore (Bos indicus) cows. Blood samples were analyzed from 34 Holstein, 104 Nellore, and 99 Dairy Nellore cows using PCR-RFLP. The different RFLP patterns were sequenced and a novel sequence variation on nucleotide-46 was identified. An adenine at this position was designated as the A allele and a guanine was designated B allele. The frequencies of alleles A-1689, A-46, and A+15 differed between Holstein and both Nellore breeds. The results show that differences in alpha-lactalbumin allelic variants in the 5`-flanking and the 5`-UTR region might be associated with differences in milk production between Holstein cows and cows from Nellore breeds. However, the lack of difference between Nellore and Dairy Nellore suggests that other sequence variantions that regulate milk production might be responsible for the selection of Dairy Nellore cows with superior milk production.

OSIRISv1.2: a named entity recognition system for sequence variants of genes in biomedical literature

Background: Single Nucleotide Polymorphisms, among other type of sequence variants, constitute key elements in genetic epidemiology and pharmacogenomics. While sequence data about genetic variation is found at databases such as dbSNP, clues about the functional and phenotypic consequences of the variations are generally found in biomedical literature. The identification of the relevant documents and the extraction of the information from them are hampered by the large size of literature databases and the lack of widely accepted standard notation for biomedical entities. Thus, automatic systems for the identification of citations of allelic variants of genes in biomedical texts are required. Results: Our group has previously reported the development of OSIRIS, a system aimed at the retrieval of literature about allelic variants of genes http://ibi.imim.es/osirisform.html. Here we describe the development of a new version of OSIRIS (OSIRISv1.2, http://ibi.imim.es/OSIRISv1.2.html webcite) which incorporates a new entity recognition module and is built on top of a local mirror of the MEDLINE collection and HgenetInfoDB: a database that collects data on human gene sequence variations. The new entity recognition module is based on a pattern-based search algorithm for the identification of variation terms in the texts and their mapping to dbSNP identifiers. The performance of OSIRISv1.2 was evaluated on a manually annotated corpus, resulting in 99% precision, 82% recall, and an F-score of 0.89. As an example, the application of the system for collecting literature citations for the allelic variants of genes related to the diseases intracranial aneurysm and breast cancer is presented. Conclusion: OSIRISv1.2 can be used to link literature references to dbSNP database entries with high accuracy, and therefore is suitable for collecting current knowledge on gene sequence variations and supporting the functional annotation of variation databases. The application of OSIRISv1.2 in combination with controlled vocabularies like MeSH provides a way to identify associations of biomedical interest, such as those that relate SNPs with diseases.

Genetic variations of EBV-LMP1 from nasopharyngeal carcinoma biopsies: potential loss of T cell epitopes

To find Epstein-Barr virus (EBV) strains with genetic variations of EBV latent membrane protein 1 (EBV-LMP1) from nasopharyngeal carcinoma (NPC), the full-length DNA of LMP1 genes from 21 NPC biopsies obtained in Hunan province in southern China was amplified and sequenced. Our sequences were compared to those previously reported by the Clustal V method. Results showed that all 21 sequences displayed two amino acid changes most frequently in LMP1 of CD4+ T cell epitopes at codons 144 (F®I, 21/21) and 212 (G®S, 19/21) or (G®N, 2/21). We also show that type A EBV strain is prevalent in the cases of NPC from Hunan province with a 30-bp 18/21 deletion, and we highlight that this deletion resulted in loss of one of the CD4+ T cell-restricted epitopes. The other 3 sequences without this deletion all had a change at codon 344 (G®D). Furthermore, in the major epitope sequence of CD8+ T cells restricted by HLA-A2, all 21 sequences showed changes at codons 126 (L®F) and 129 (M®I). Our study discovered that one of the 21 sequence variations harbored a new change at codon 131 (W®C), and 5/21 specimens showed another novel change at codon 115 (G®A) in the major epitope sequence of CD8+ T cells restricted by HLA-A2. Our study suggests that these sequence variations of NPC-derived LMP1 may lead to a potential escape from host cell immune recognition, protecting latent EBV infection and causing an increase in tumorigenicity.

PPARGC1A sequence variation and cardiovascular risk-factor levels: a study of the main genetic effects and gene x environment interactions in children from the European Youth Heart Study.

AIMS/HYPOTHESIS: The PPARGC1A gene coactivates multiple nuclear transcription factors involved in cellular energy metabolism and vascular stasis. In the present study, we genotyped 35 tagging polymorphisms to capture all common PPARGC1A nucleotide sequence variations and tested for association with metabolic and cardiovascular traits in 2,101 Danish and Estonian boys and girls from the European Youth Heart Study, a multicentre school-based cross-sectional cohort study. METHODS: Fasting plasma glucose concentrations, anthropometric variables and blood pressure were measured. Habitual physical activity and aerobic fitness were objectively assessed using uniaxial accelerometry and a maximal aerobic exercise stress test on a bicycle ergometer, respectively. RESULTS: In adjusted models, nominally significant associations were observed for BMI (rs10018239, p = 0.039), waist circumference (rs7656250, p = 0.012; rs8192678 [Gly482Ser], p = 0.015; rs3755863, p = 0.02; rs10018239, beta = -0.01 cm per minor allele copy, p = 0.043), systolic blood pressure (rs2970869, p = 0.018) and fasting glucose concentrations (rs11724368, p = 0.045). Stronger associations were observed for aerobic fitness (rs7656250, p = 0.005; rs13117172, p = 0.008) and fasting glucose concentrations (rs7657071, p = 0.002). None remained significant after correcting for the number of statistical comparisons. We proceeded by testing for gene x physical activity interactions for the polymorphisms that showed nominal evidence of association in the main effect models. None of these tests was statistically significant. CONCLUSIONS/INTERPRETATION: Variants at PPARGC1A may influence several metabolic traits in this European paediatric cohort. However, variation at PPARGC1A is unlikely to have a major impact on cardiovascular or metabolic health in these children.

Human MAMLD1 Gene Variations Seem Not Sufficient to Explain a 46,XY DSD Phenotype

MAMLD1 is thought to cause disordered sex development in 46,XY patients. But its role is controversial because some MAMLD1 variants are also detected in normal individuals, several MAMLD1 mutations have wild-type activity in functional tests, and the male Mamld1-knockout mouse has normal genitalia and reproduction. Our aim was to search for MAMLD1 variations in 108 46,XY patients with disordered sex development, and to test them functionally. We detected MAMDL1 variations and compared SNP frequencies in controls and patients. We tested MAMLD1 transcriptional activity on promoters involved in sex development and assessed the effect of MAMLD1 on androgen production. MAMLD1 expression in normal steroid-producing tissues and mutant MAMLD1 protein expression were also assessed. Nine MAMLD1 mutations (7 novel) were characterized. In vitro, most MAMLD1 variants acted similarly to wild type. Only the L210X mutation showed loss of function in all tests. We detected no effect of wild-type or MAMLD1 variants on CYP17A1 enzyme activity in our cell experiments, and Western blots revealed no significant differences for MAMLD1 protein expression. MAMLD1 was expressed in human adult testes and adrenals. In conclusion, our data support the notion that MAMLD1 sequence variations may not suffice to explain the phenotype in carriers and that MAMLD1 may also have a role in adult life.

Genotype-phenotype correlation in Brazillian Rett syndrome patients

BACKGROUND: Rett syndrome (RS) is a severe neurodevelopmental X-linked dominant disorder caused by mutations in the MECP2 gene. PURPOSE: To search for point mutations on the MECP2 gene and to establish a correlation between the main point mutations found and the phenotype. METHOD: Clinical evaluation of 105 patients, following a standard protocol. Detection of point mutations on the MECP2 gene was performed on peripheral blood DNA by sequencing the coding region of the gene. RESULTS: Classical RS was seen in 68% of the patients. Pathogenic point mutations were found in 64.1% of all patients and in 70.42% of those with the classical phenotype. Four new sequence variations were found, and their nature suggests patogenicity. Genotype-phenotype correlations were performed. CONCLUSION: Detailed clinical descriptions and identification of the underlying genetic alterations of this Brazilian RS population add to our knowledge of genotype/phenotype correlations, guiding the implementation of mutation searching programs.

Disruption of Urate Transport in Familial Renal Glucosuria and Report on SGLT2 Expression in Normal and Pathological Kidney

Familial renal glucosuria (FRG) is a rare co -dominantly inherited benign phenotype characterized by the presence of glucose in the urine. It is caused by mutations in the SLC5A2 gene that encodes SGLT2, a Na+ -glucose co -transporter. The purpose of our current work was twofold: to characterize the molecular and phenotype findings of an FRG cohort and, in addition, to detail the SGLT2 expression in the adult human kidney. The phenotype of FRG pedigrees was evaluated using direct sequencing for the identification of sequence variations in the SLC5A2 gene. The expression of SGLT2 in the adult human kidney was studied by immunofluorescence on kidney biopsy specimens. In the absence of renal biopsies from FRG individuals, and in order to evaluate the potential disruption of SGLT2 expression in a glucosuric nephropathy, we have selected cases of nucleoside analogues induced proximal tubular toxicity. We identified six novel SLC5A2 mutations in six FRG pedigrees and described the occurrence of hyperuricosuria associated with hypouricaemia in the two probands with the most severe phenotypes. Histopathological studies proved that SGLT2 is localized to the brush -border of the proximal tubular epithelia cell and that this normal pattern was found to be disrupted in cases of nucleoside analogues induced tubulopathy. We present six novel SLC5A2 mutations, further contributing to the allelic heterogeneity in FRG, and identified hyperuricosuria and hypouricaemia as part of the FRG phenotype. SGLT2 is localized to the brush -border of the proximal tubule in the adult human normal kidney, and aberrant expression of the co -transporter may underlie the glucosuria seen with the use of nucleoside analogues.

Mutation screening of the glutamate cysteine ligase modifier (GCLM) gene in patients with schizophrenia.

BACKGROUND: Experimental evidences show that glutathione and its rate-limiting synthesizing enzyme, the glutamate-cysteine ligase (GCL), are involved in the pathogenesis of schizophrenia. Furthermore, genetic association has been previously reported between two single nucleotide polymorphisms lying in noncoding regions of glutamate cysteine ligase modifier (GCLM) gene, which specifies for the modifier subunit of GCL and schizophrenia. OBJECTIVE: We wanted to investigate the presence of GCLM true functional mutations, likely in linkage disequilibrium with the previously identified single nucleotide polymorphism alleles, in the same set of cases that allowed the detection of the original association signal. METHODS: We screened all the coding regions of GCLM and their intronic flanking vicinities in 353 patients with schizophrenia by direct DNA sequencing. RESULTS: Ten sequence variations were identified, five of which were not previously described. None of these DNA changes was within the GCLM coding sequence and in-silico analysis failed to indicate functional impairment induced by these variations. Furthermore, screening of normal controls and downstream statistical analyses revealed no significant relationship of any of these DNA variants with schizophrenia. CONCLUSION: It is unlikely that functional mutations in the GCLM gene could play a major role in genetic predisposition to schizophrenia and further studies will be required to assess its etiological function in the disease.

Prevalence of human papillomavirus type 16 variants in the Federal District, Central Brazil

We report the prevalence of human papillomavirus type 16 (HPV-16) variants in women with cervical lesions from the Federal District, Central Brazil. We analyzed 34 HPV-16 samples, identifying the sequence variations of E6 and L1 genes and correlating variant frequency with disease status. The most prevalent HPV-16 variant was the European (50%), followed by Asian-American (41.2%), African-1 (5.9%), and African-2 (2.9%). European and non-European variants appeared in equal frequencies among the cytological types of lesions - atypical squamous or glandular cells of undetermined significance, cytological alterations suggesting HPV infection, cervical intraepithelial neoplasias, squamous cell carcinoma, and adenocarcinoma.

Genetic variability in the 5' UTR and NS5A regions of hepatitis C virus RNA isolated from non-responding and responding patients with chronic HCV genotype 1 infection

Sequence variation among different hepatitis C virus (HCV) isolates has adaptive significance and reflects the modes and intensities of selection mechanisms operating on the virus. In this work, we sought to investigate using classical population genetics parameters, the genetic variability of HCV genotype 1 using the 5' UTR and NS5A regions from treatment non-responding and responding groups of patients. Both regions showed low genetic varia-bility and the 5' UTR showed neutral deviation. No differences were observed in the nonsynonymous/synonymous nucleotide substitution ratio among groups for NS5A. The analysis of molecular variance test of the 5' UTR region showed an 11.94% variation among groups. Phylogenetic analysis showed no correlation between sequence variations and therapeutic responses.

Allelic variation in the VMD2 gene in best disease and age-related macular degeneration.

PURPOSE: To assess the allelic variation of the VMD2 gene in patients with Best disease and age-related macular degeneration (AMD). METHODS: Three hundred twenty-one AMD patients, 192 ethnically similar control subjects, 39 unrelated probands with familial Best disease, and 57 unrelated probands with the ophthalmoscopic findings of Best disease but no family history were screened for sequence variations in the VMD2 gene by single-strand conformation polymorphism (SSCP) analysis. Amplimers showing a bandshift were reamplified and sequenced bidirectionally. In addition, the coding regions of the VMD2 gene were completely sequenced in six probands with familial Best disease who showed no SSCP shift. RESULTS: Forty different probable or possible disease-causing mutations were found in one or more Best disease or AMD patients. Twenty-nine of these variations are novel. Of the 39 probands with familial Best disease, mutations were detected in all 39 (33 by SSCP and 6 by DNA sequencing). SSCP screening of the 57 probands with a clinical diagnosis of Best disease but no family history revealed 16 with mutations. Mutations were found in 5 of 321 AMD patients (1.5%), a fraction that was not significantly greater than in control individuals (0/192, 0%). CONCLUSIONS: Patients with the clinical diagnosis of Best disease are significantly more likely to have a mutation in the VMD2 gene if they also have a positive family history. These findings suggest that a small fraction of patients with the clinical diagnosis of AMD may actually have a late-onset variant of Best disease, whereas at the same time, a considerable fraction of isolated patients with the ophthalmoscopic features of Best disease are probably affected with some other macular disease.

An integrated genetic and functional analysis of the role of type II transmembrane serine proteases (TMPRSSs) in hearing loss.

Building on our discovery that mutations in the transmembrane serine protease, TMPRSS3, cause nonsyndromic deafness, we have investigated the contribution of other TMPRSS family members to the auditory function. To identify which of the 16 known TMPRSS genes had a strong likelihood of involvement in hearing function, three types of biological evidence were examined: 1) expression in inner ear tissues; 2) location in a genomic interval that contains a yet unidentified gene for deafness; and 3) evaluation of hearing status of any available Tmprss knockout mouse strains. This analysis demonstrated that, besides TMPRSS3, another TMPRSS gene was essential for hearing and, indeed, mice deficient for Hepsin (Hpn) also known as Tmprss1 exhibited profound hearing loss. In addition, TMPRSS2, TMPRSS5, and CORIN, also named TMPRSS10, showed strong likelihood of involvement based on their inner ear expression and mapping position within deafness loci PKSR7, DFNB24, and DFNB25, respectively. These four TMPRSS genes were then screened for mutations in affected members of the DFNB24 and DFNB25 deafness families, and in a cohort of 362 sporadic deaf cases. This large mutation screen revealed numerous novel sequence variations including three potential pathogenic mutations in the TMPRSS5 gene. The mutant forms of TMPRSS5 showed reduced or absent proteolytic activity. Subsequently, TMPRSS genes with evidence of involvement in deafness were further characterized, and their sites of expression were determined. Tmprss1, 3, and 5 proteins were detected in spiral ganglion neurons. Tmprss3 was also present in the organ of Corti. TMPRSS1 and 3 proteins appeared stably anchored to the endoplasmic reticulum membranes, whereas TMPRSS5 was also detected at the plasma membrane. Collectively, these results provide evidence that TMPRSS1 and TMPRSS3 play and TMPRSS5 may play important and specific roles in hearing.