Mendelian randomization and mediation analysis of leukocyte telomere length and risk of lung and head and neck cancers

Background: Evidence from observational studies of telomere length (TL) has been conflicting regarding its direction of association with cancer risk. We investigated the causal relevance of TL for lung and head and neck cancers using Mendelian Randomization (MR) and mediation analyses. Methods: We developed a novel genetic instrument for TL in chromosome 5p15.33, using variants identified through deep-sequencing, that were genotyped in 2051 cancer-free subjects. Next, we conducted an MR analysis of lung (16 396 cases, 13 013 controls) and head and neck cancer (4415 cases, 5013 controls) using eight genetic instruments for TL. Lastly, the 5p15.33 instrument and distinct 5p15.33 lung cancer risk loci were evaluated using two-sample mediation analysis, to quantify their direct and indirect, telomere-mediated, effects. Results: The multi-allelic 5p15.33 instrument explained 1.49-2.00% of TL variation in our data (p = 2.6 × 10-9). The MR analysis estimated that a 1000 base-pair increase in TL increases risk of lung cancer [odds ratio (OR) = 1.41, 95% confidence interval (CI): 1.20-1.65] and lung adenocarcinoma (OR = 1.92, 95% CI: 1.51-2.22), but not squamous lung carcinoma (OR = 1.04, 95% CI: 0.83-1.29) or head and neck cancers (OR = 0.90, 95% CI: 0.70-1.05). Mediation analysis of the 5p15.33 instrument indicated an absence of direct effects on lung cancer risk (OR = 1.00, 95% CI: 0.95-1.04). Analysis of distinct 5p15.33 susceptibility variants estimated that TL mediates up to 40% of the observed associations with lung cancer risk. Conclusions: Our findings support a causal role for long telomeres in lung cancer aetiology, particularly for adenocarcinoma, and demonstrate that telomere maintenance partially mediates the lung cancer susceptibility conferred by 5p15.33 loci

Association of apolipoprotein E polymorphism in late-onset Alzheimer's disease and vascular dementia in Brazilians

The genetic basis for dementias is complex. A common polymorphism in the apolipoprotein E (APOE) gene is considered to be the major risk factor in families with sporadic and late-onset Alzheimer's disease as well as in the general population. The distribution of alleles and genotypes of the APOE gene in late-onset Alzheimer's disease (N = 68), other late-life dementias (N = 39), and in cognitively normal controls (N = 58) was determined, as also was the risk for Alzheimer's disease associated with the epsilon4 allele. Peripheral blood samples were obtained from a total of 165 individuals living in Brazil aged 65-82 years. Genomic DNA was amplified by the polymerase chain reaction and the products were digested with HhaI restriction enzyme. APOE epsilon2 frequency was considerably lower in the Alzheimer's disease group (1%), and the epsilon3 allele and epsilon3/epsilon3 genotype frequencies were higher in the controls (84 and 72%, respectively) as were the epsilon4 allele and epsilon3/epsilon4 genotype frequencies in Alzheimer's disease (25 and 41%, respectively). The higher frequency of the epsilon4 allele in Alzheimer's disease confirmed its role as a risk factor, while epsilon2 provided a weak protection against development of the disease. However, in view of the unexpectedly low frequency of the epsilon4 allele, additional analyses in a more varied Brazilian sample are needed to clarify the real contribution of apolipoprotein E to the development of Alzheimer's disease in this population

High frequencies of plexiform neurofibromas, mental retardation, learning difficulties, and scoliosis in Brazilian patients with neurofibromatosis type 1

A clinical study of Brazilian patients with neurofibromatosis type 1 (NF1) was performed in a multidisciplinary Neurofibromatosis Program called CEPAN (Center of Research and Service in Neurofibromatosis). Among 55 patients (60% females, 40% males) who met the NIH criteria for the diagnosis of NF1, 98% had more than six café-au-lait patches, 94.5% had axillary freckling, 45% had inguinal freckling, and 87.5% had Lisch nodules. Cutaneous neurofibromas were observed in 96%, and 40% presented plexiform neurofibromas. A positive family history of NF1 was found in 60%, and mental retardation occurred in 35%. Some degree of scoliosis was noted in 49%, 51% had macrocephaly, 40% had short stature, 76% had learning difficulties, and 2% had optic gliomas. Unexpectedly high frequencies of plexiform neurofibromas, mental retardation, learning difficulties, and scoliosis were observed, probably reflecting the detailed clinical analysis methods adopted by the Neurofibromatosis Program. These same patients were screened for mutations in the GAP-related domain/GRD (exons 20-27a) by single-strand conformation polymorphism. Four different mutations (Q1189X, 3525-3526delAA, E1356G, c.4111-1G>A) and four polymorphisms (c.3315-27G>A, V1146I, V1317A, c.4514+11C>G) were identified. These data were recently published

High frequencies of plexiform neurofibromas, mental retardation, learning difficulties, and scoliosis in Brazilian patients with neurofibromatosis type 1

A clinical study of Brazilian patients with neurofibromatosis type 1 (NF1) was performed in a multidisciplinary Neurofibromatosis Program called CEPAN (Center of Research and Service in Neurofibromatosis). Among 55 patients (60% females, 40% males) who met the NIH criteria for the diagnosis of NF1, 98% had more than six café-au-lait patches, 94.5% had axillary freckling, 45% had inguinal freckling, and 87.5% had Lisch nodules. Cutaneous neurofibromas were observed in 96%, and 40% presented plexiform neurofibromas. A positive family history of NF1 was found in 60%, and mental retardation occurred in 35%. Some degree of scoliosis was noted in 49%, 51% had macrocephaly, 40% had short stature, 76% had learning difficulties, and 2% had optic gliomas. Unexpectedly high frequencies of plexiform neurofibromas, mental retardation, learning difficulties, and scoliosis were observed, probably reflecting the detailed clinical analysis methods adopted by the Neurofibromatosis Program. These same patients were screened for mutations in the GAP-related domain/GRD (exons 20-27a) by single-strand conformation polymorphism. Four different mutations (Q1189X, 3525-3526delAA, E1356G, c.4111-1G>A) and four polymorphisms (c.3315-27G>A, V1146I, V1317A, c.4514+11C>G) were identified. These data were recently published

Investigating the effect of sexual behaviour on oropharyngeal cancer risk: a methodological assessment of Mendelian randomization

Background Human papilloma virus infection is known to influence oropharyngeal cancer (OPC) risk, likely via sexual transmission. However, sexual behaviour has been correlated with other risk factors including smoking and alcohol, meaning independent effects are difficult to establish. We aimed to evaluate the causal effect of sexual behaviour on the risk of OPC using Mendelian randomization (MR). Methods Genetic variants robustly associated with age at first sex (AFS) and the number of sexual partners (NSP) were used to perform both univariable and multivariable MR analyses with summary data on 2641 OPC cases and 6585 controls, obtained from the largest available genome-wide association studies (GWAS). Given the potential for genetic pleiotropy, we performed a number of sensitivity analyses: (i) MR methods to account for horizontal pleiotropy, (ii) MR of sexual behaviours on positive (cervical cancer and seropositivity for Chlamydia trachomatis) and negative control outcomes (lung and oral cancer), (iii) Causal Analysis Using Summary Effect estimates (CAUSE), to account for correlated and uncorrelated horizontal pleiotropic effects, (iv) multivariable MR analysis to account for the effects of smoking, alcohol, risk tolerance and educational attainment. Results In univariable MR, we found evidence supportive of an effect of both later AFS (IVW OR = 0.4, 95%CI (0.3, 0.7), per standard deviation (SD), p = < 0.001) and increasing NSP (IVW OR = 2.2, 95%CI (1.3, 3.8) per SD, p = < 0.001) on OPC risk. These effects were largely robust to sensitivity analyses accounting for horizontal pleiotropy. However, negative control analysis suggested potential violation of the core MR assumptions and subsequent CAUSE analysis implicated pleiotropy of the genetic instruments used to proxy sexual behaviours. Finally, there was some attenuation of the univariable MR results in the multivariable models (AFS IVW OR = 0.7, 95%CI (0.4, 1.2), p = 0.21; NSP IVW OR = 0.9, 95%CI (0.5 1.7), p = 0.76). Conclusions Despite using genetic variants strongly related sexual behaviour traits in large-scale GWAS, we found evidence for correlated pleiotropy. This emphasizes a need for multivariable approaches and the triangulation of evidence when performing MR of complex behavioural traits

Identification And Complete Sequencing Of Novel Human Transcripts Through The Use Of Mouse Orthologs And Testis Cdna Sequences

The correct identification of all human genes, and their derived transcripts, has not yet been achieved, and it remains one of the major aims of the worldwide genomics community. Computational programs suggest the existence of 30,000 to 40,000 human genes. However, definitive gene identification can only be achieved by experimental approaches. We used two distinct methodologies, one based on the alignment of mouse orthologous sequences to the human genome, and another based on the construction of a high-quality human testis cDNA library, in an attempt to identify new human transcripts within the human genome sequence. We generated 47 complete human transcript sequences, comprising 27 unannotated and 20 annotated sequences. Eight of these transcripts are variants of previously known genes. These transcripts were characterized according to size, number of exons, and chromosomal localization, and a search for protein domains was undertaken based on their putative open reading frames. In silico expression analysis suggests that some of these transcripts are expressed at low levels and in a restricted set of tissues.34493511Adams, M.D., Kelley, J.M., Gocayne, J.D., Dubnick, M., Polymeropoulos, M.H., Xiao, H., Merril, C.R., Moreno, R.F., Complementary DNA sequencing: Expressed sequence tags and human genome project (1991) Science, 252, pp. 1651-1656Blanco, E., Parra, G., Guigó, R., Finding genes (2002) Current Protocols in Bioinformatics, , (Baxevanis, A., ed.). John Wiley & Sons Ltd., New York (in press)Bonaldo, M.F., Lennon, G., Soares, M.B., Normalization and subtraction: Two approaches to facilitate gene discovery (1996) Genome Res., 6, pp. 791-806Boon, K., Osorio, E.C., Greenhut, S.F., Schaefer, C.F., Shoemaker, J., Polyak, K., Morin, P.J., Riggins, G.J., An anatomy of normal and malignant gene expression (2002) Proc. Natl. Acad. Sci. USA, 99, pp. 11287-11292Brentani, H., Caballero, O.L., Camargo, A.A., Da Silva, A.M., Da Silva Jr., W.A., Dias Neto, E., Grivet, M., Luna, A.M., The generation and utilization of a cancer-oriented representation of the human transcriptome by using expressed sequence tags (2003) Proc. Natl. Acad. Sci. USA, 100, pp. 13418-13423Burge, C., Karlin, S., Prediction of complete gene structures in human genomic DNA (1997) J. Mol. Biol., 268, pp. 78-94Burset, M., Guigó, R., Evaluation of gene structure prediction programs (1996) Genomics, 34, pp. 353-367Camargo, A.A., Samaia, H.P., Dias-Neto, E., Simao, D.F., Migotto, I.A., Briones, M.R., Costa, F.F., Simpson, A.J., The contribution of 700,000 ORF sequence tags to the definition of the human transcriptome (2001) Proc. Natl. Acad. Sci. USA, 98, pp. 12103-12108Carninci, P., Shibata, Y., Hayatsu, N., Sugahara, Y., Shibata, K., Itoh, M., Konno, H., Hayashizaki, Y., Normalization and subtraction of cap-trapper-selected cDNA molecules to prepare full-length cDNA libraries for rapid discovery of new genes (2000) Genome Res., 10, pp. 1617-1630Carninci, P., Shibata, Y., Hayatsu, N., Itoh, M., Shiraki, T., Hirozane, T., Watahiki, A., Hayashizaki, Y., Balanced-size and long-size cloning of full-length, cap-trapped cDNAs into vectors of the novel lambda-FLC family allows enhanced gene discovery rate and functional analysis (2001) Genomics, 77, pp. 79-90Chirgwin, J.M., Przybyla, A.E., McDonald, R.J., Rutter, W.J., Isolation of biologically active ribonucleic acid from sources enriched in ribonuclease (1979) Biochemistry, 18, pp. 5294-5299De Souza, S.J., Camargo, A.A., Briones, M.R., Costa, F.F., Nagai, M.A., Verjovski-Almeida, S., Zago, M.A., Simpson, A.J., Identification of human chromosome 22 transcribed sequences with ORF expressed sequence tags (2000) Proc. Natl. Acad. Sci. USA, 97, pp. 12690-12693Dias-Neto, E., Correa, R.G., Verjovski-Almeida, S., Briones, M.R., Nagai, M.A., Da Silva Jr., W., Zago, M.A., Simpson, A.J., Shotgun sequencing of the human transcriptome with ORF expressed sequence tags (2000) Proc. Natl. Acad. Sci. USA, 97, pp. 3491-3496Guigó, R., Knudsen, S., Drake, N., Smith, T., Prediction of gene structure (1992) J. Mol. Biol., 226, pp. 141-157Henke, W., Herdel, K., Jung, K., Schoor, D., Lorning, S.A., Betaine improves the PCR amplification of GC-rich DNA sequences (1997) Nucleic Acids Res., 25, pp. 3957-3958Hickox, D.M., Gibbs, G., Morrison, J.R., Sebire, K., Edgar, K., Keah, H.H., Alter, K., O'Bryan, M.K., Identification of a novel testis-specific member of the phosphatidylethanolamine binding protein family, pebp-2 (2002) Biol. Reprod., 67, pp. 917-927Houlgatte, R., Mariage-Samson, R., Duprat, S., Tessier, A., Bentolila, S., Lamy, B., Auffray, C., The Genexpress Index: A resource for gene discovery and the genic map of the human genome (1995) Genome Res., 5, pp. 272-304Huang, X., Madan, A., CAP3: A DNA sequence assembly program (1999) Genome Res., 9, pp. 868-877Kent, W.J., Sugnet, C.W., Furey, T.S., Roskin, K.M., Pringle, T.H., Zahler, A.M., Haussler, D., The Human Genome Browser at UCSC (2002) Genome Res., 12, pp. 996-1006Kikuno, R., Nagase, T., Waki, M., Ohara, O., HUGE: A database for human large proteins identified in the Kazusa cDNA sequencing project (2002) Nucleic Acids Res., 30, pp. 166-168Lander, E.S., Linton, L.M., Birren, B., Nusbaum, C., Zody, M.C., Baldwin, J., Devon, K., Uberbacher, E., Initial sequencing and analysis of the human genome (2001) Nature, 409, pp. 860-921Liang, F., Holt, I., Pertea, G., Karamycheva, S., Salzberg, S.L., Quackenbush, J., Gene index analysis of the human genome estimates approximately 120,000 genes (2000) Nature, 25, pp. 239-240Makalowski, W., Zhang, J., Boguski, M.S., Comparative analysis of 1196 orthologous mouse and human full-length mRNA and protein sequences (1996) Genome Res., 6, pp. 846-857Modrek, B., Lee, C., A genomic view of alternative splicing (2002) Nat. Genet., 30, pp. 13-19Nakajima, D., Okazaki, N., Yamakawa, H., Kikuno, R., Ohara, O., Nagase, T., Construction of expression-ready cDNA clones for KIAA genes: Manual curation of 330 KIAA cDNA clones (2002) DNA Res., 9, pp. 99-106Rogic, S., Mackworth, A.K., Ouellette, F.B., Evaluation of gene-finding programs on mammalian sequences (2001) Genome Res., 11, pp. 817-832Sakharkar, M.K., Kangueane, P., Genome SEGE: A database for 'intronless' genes in eukaryotic genomes (2004) BMC Bioinformatics, 5, p. 67Sakharkar, M.K., Kangueane, P., Petrov, D.A., Kolaskar, A.S., Subbiah, S., SEGE: A database on 'intron less/single exonic' genes from eukaryotes (2002) Bioinformatics, 18, pp. 1266-1267Sambrook, J., Fritsch, E., Maniatis, T., (1989) Molecular Cloning, , Cold Spring Harbor Laboratory Press, Cold Spring Harbor, NY, USAShabalina, S.A., Ogurtsov, A.Y., Rogozin, I.B., Koonin, E.V., Lipman, D.J., Comparative analysis of orthologous eukaryotic mRNAs: Potential hidden functional signals (2004) Nucleic Acids Res., 32, pp. 1774-1782Shendure, J., Church, G.M., Computational discovery of sense-antisense transcription in the human and mouse genomes (2002) Genome Biol., 3, pp. 1-14Sogayar, M.C., Camargo, A.A., Bettoni, F., Carraro, D.M., Pires, L.C., Parmigiani, R.B., Ferreira, E.N., Zanette, D.L., A transcript finishing initiative for closing gaps in the human transcriptome (2004) Genome Res., 14, pp. 1413-1423Solovyev, V., Statistical approaches in eukaryotic gene prediction (2001) Handbook of Statistical Genetics, pp. 83-127. , (Balding, D.J., Bishop, M. and Cannings, C., eds.). John Wiley & Sons Ltd., UKSorek, R., Safer, H.M., A novel algorithm for computational identification of contaminated EST libraries (2003) Nucleic Acids Res., 31, pp. 1067-1074Strausberg, R.L., Feingold, E.A., Klausner, R.D., Collins, F.S., The mammalian gene collection (1999) Science, 286, pp. 455-457Strausberg, R.L., Feingold, E.A., Grouse, L.H., Derge, J.G., Klausner, R.D., Collins, F.S., Wagner, L., Marra, M.A., Generation and initial analysis of more than 15,000 full-length human and mouse cDNA sequences (2002) Proc. Natl. Acad. Sci. USA, 99, pp. 16899-16903Genome sequence of the nematode C. elegans: A platform for investigating biology (1998) Science, 282, pp. 2012-2018Velculescu, V.E., Zhang, L., Vogelstein, B., Kinzler, K.W., Serial analysis of gene expression (1995) Science, 270, pp. 484-487Venter, J.C., Adams, M.D., Myers, E.W., Li, P.W., Mural, R.J., Sutton, G.G., Smith, H.O., Xiao, C., The Sequence of the Human Genome (2001) Science, 291, pp. 1304-1351Vettore, A.L., Da Silva, F.R., Kemper, E.L., Arruda, P., The libraries that made SUCEST (2001) Gen. Mol. Biol., 24, pp. 1-7Warrington, J.A., Nair, A., Mahadevappa, M., Tsyganskaya, M., Comparison of human adult and fetal expression and identification of 535 housekeeping/maintenance genes (2000) Physiol. Genomics, 2, pp. 143-147Wiemann, S., Weil, B., Wellenreuther, R., Gassenhuber, J., Glassl, S., Ansorge, W., Bocher, M., Poustka, A., Toward a catalog of human genes and proteins: Sequencing and analysis of 500 novel complete protein coding human cDNAs (2001) Genome Res., 11, pp. 422-435Yao, J., Chiba, T., Sakai, J., Hirose, K., Yamamoto, M., Hada, A., Kuramoto, K., Mori, M., Mouse testis transcriptome revealed using serial analysis of gene expression (2004) Mamm. Genome, 15, pp. 433-451Yelin, R., Dahary, D., Sorek, R., Levanon, E.Y., Goldstein, O., Shoshan, A., Diber, A., Rotman, G., Widespread occurrence of antisense transcription in the human genome (2003) Nat. Biotechnol., 21, pp. 379-38