Search

Supplementary Tables 1-7

Figure S1

We investigated a Spanish and Catalan family in which multiple cancer types tracked across three generations, but for which no genetic etiology had been identified. Whole-exome sequencing of germline DNA from multiple affected family members was performed to identify candidate variants to explain this occurrence of familial cancer. We discovered in all cancer-affected family members a single rare heterozygous germline variant (I654V, rs1801201) in ERBB2/HER2, which is located in a transmembrane glycine zipper motif critical for ERBB2-mediated signaling and in complete linkage disequilibrium (D' = 1) with a common polymorphism (I655V, rs1136201) previously reported in some populations as associated with cancer risk. Because multiple cancer types occurred in this family, we tested both the I654V and the I655V variants for association with cancer across multiple tumor types in 6,371 cases of Northern European ancestry drawn from The Cancer Genome Atlas and 6,647 controls, and found that the rare variant (I654V) was significantly associated with an increased risk for cancer (OR = 1.40; P = 0.021; 95% confidence interval (CI), 1.05–1.89). Functional assays performed in HEK 293T cells revealed that both the I655V single mutant (SM) and the I654V;I655V double mutant (DM) stabilized ERBB2 protein and activated ERBB2 signaling, with the DM activating ERBB2 significantly more than the SM alone. Thus, our results suggest a model whereby heritable genetic variation in the transmembrane domain activating ERBB2 signaling is associated with both sporadic and familial cancer risk, with increased ERBB2 stabilization and activation associated with increased cancer risk.Prevention Relevance:By performing whole-exome sequencing on germline DNA from multiple cancer-affected individuals belonging to a family in which multiple cancer types track across three generations, we identified and then characterized functional common and rare variation in ERBB2 associated with both sporadic and familial cancer. Our results suggest that heritable variation activating ERBB2 signaling is associated with risk for multiple cancer types, with increases in signaling correlated with increases in risk, and modified by ancestry or family history.

Molecular markers of sorafenib efficacy in patients with metastatic renal cell carcinoma (mRCC) are not available. The purpose of this study was to discover genetic markers of survival in patients with mRCC treated with sorafenib. Germline variants from 56 genes were genotyped in 295 patients with mRCC. Variant-overall survival (OS) associations were tested in multivariate regression models. Mechanistic studies were conducted to validate clinical associations. VEGFA rs1885657, ITGAV rs3816375, and WWOX rs8047917 (sorafenib arm), and FLT4 rs307826 and VEGFA rs3024987 (sorafenib and placebo arms combined) were associated with shorter OS. FLT4 rs307826 increased VEGFR-3 phosphorylation, membrane trafficking, and receptor activation. VEGFA rs1885657 and rs58159269 increased transcriptional activity of the constructs containing these variants in endothelial and RCC cell lines, and VEGFA rs58159269 increased endothelial cell proliferation and tube formation. FLT4 rs307826 and VEGFA rs58159269 led to reduced sorafenib cytotoxicity. Genetic variation in VEGFA and FLT4 could affect survival in sorafenib-treated patients with mRCC. These markers should be examined in additional malignancies treated with sorafenib and in other angiogenesis inhibitors used in mRCC.Significance:Clinical and mechanistic data identify germline genetic variants in VEGFA and FLT4 as markers of survival in patients with metastatic renal cell carcinoma.

Genotyped genes.

Expression of WT and T494A VEGFR-3 in HUVECs (A) and replication of phosphorylation experiments testing the effect of FLT4 rs307826 (A>G, T494A) on VEGFR-3 phosphorylation (B).

Supplementary Tables 1-6. Supplementary Table 1. List of genes and variants previously reported as associated with NPC. Supplementary Table 2. Proportion of exome covered at different read depths for each sample sequenced. Supplementary Table 3. Variant subtype for each samples sequenced. Supplementary Table 4. Rare deleterious variants shared in the germline of all three NPC-affected family members. Supplementary Table 6. EBV insertion events identified in the exome of each NPC-affected family member.

SNP selection criteria, imputation and bioinformatic analyses

Supplementary Table 5: Somatic Mutations Identified in Individual I-1.

No biomarkers are available to predict patients at risk of developing hypertension induced by VEGF-pathway inhibitors. This study aimed to identify predictive biomarkers of hypertension induced by these drugs using a discovery-replication approach. The discovery set included 140 sorafenib-treated patients (TARGET study) genotyped for 973 SNPs in 56 genes. The most statistically significant SNPs associated with grade ≥2 hypertension were tested for association with grade ≥2 hypertension in the replication set of a GWAS of 1039 bevacizumab-treated patients from four clinical trials (CALGB/Alliance). In the discovery set, rs444904 (G > A) in PIK3R5 was associated with an increased risk of sorafenib-induced hypertension (p = 0.006, OR = 3.88 95% CI 1.54–9.81). In the replication set, rs427554 (G > A) in PIK3R5 (in complete linkage disequilibrium with rs444904) was associated with an increased risk of bevacizumab-induced hypertension (p = 0.008, OR = 1.39, 95% CI 1.09–1.78). This study identified a predictive marker of drug-induced hypertension that should be evaluated for other VEGF-pathway inhibitors. ClinicalTrials.gov Identifier:NCT00073307 (TARGET).

Some patients with therapy-related myeloid neoplasms (t-MN) may have unsuspected inherited cancer predisposition syndrome (CPS). We propose a set of clinical criteria to identify t-MN patients with high risk of CPS (HR-CPS). Among 225 t-MN patients with an antecedent non-myeloid malignancy, our clinical criteria identified 52 (23%) HR-CPS patients. Germline whole-exome sequencing identified pathogenic or likely pathogenic variants in 10 of 27 HR-CPS patients compared to 0 of 9 low-risk CPS patients (37% vs. 0%, p = 0.04). These simple clinical criteria identify t-MN patients most likely to benefit from genetic testing for inherited CPS.

Acute myeloid leukemia (AML) is a hematological malignancy with an undefined heritable risk. Here we perform a meta-analysis of three genome-wide association studies, with replication in a fourth study, incorporating a total of 4018 AML cases and 10488 controls. We identify a genome-wide significant risk locus for AML at 11q13.2 (rs4930561; P = 2.15 × 10−8; KMT5B). We also identify a genome-wide significant risk locus for the cytogenetically normal AML sub-group (N = 1287) at 6p21.32 (rs3916765; P = 1.51 × 10−10; HLA). Our results inform on AML etiology and identify putative functional genes operating in histone methylation (KMT5B) and immune function (HLA)., Genome wide association studies in cancer are used to understand the heritable genetic contribution to disease risk. Here, the authors perform a genome wide association study in European patients with acute myeloid leukemia and identify loci associated with risk of developing the disease.

Acute myeloid leukemia (AML) is a hematological malignancy with an undefined heritable risk. Here we performed a meta-analysis of three genome-wide association studies, with replication in a fourth study, incorporating a total of 4018 AML cases and 10488 controls. We identified a genome-wide significant risk locus for AML at 11q13.2 (rs4930561; P = 2.15 × 10-8; KMT5B). We also identified a genome-wide significant risk locus for the cytogenetically normal AML sub-group (N=1287) at 6p21.32 (rs3916765; P = 1.51 × 10-10; HLA). Our results inform on AML etiology by identifying putative functional genes operating in histone methylation (KMT5B) and immune function (HLA).

We determined differential gene expression in response to high glucose in lymphoblastoid cell lines derived from matched individuals with type 1 diabetes with and without retinopathy. Those genes exhibiting the largest difference in glucose response were assessed for association with diabetic retinopathy in a genome-wide association study meta-analysis. Expression quantitative trait loci (eQTLs) of the glucose response genes were tested for association with diabetic retinopathy. We detected an enrichment of the eQTLs from the glucose response genes among small association p-values and identified folliculin (FLCN) as a susceptibility gene for diabetic retinopathy. Expression of FLCN in response to glucose was greater in individuals with diabetic retinopathy. Independent cohorts of individuals with diabetes revealed an association of FLCN eQTLs with diabetic retinopathy. Mendelian randomization confirmed a direct positive effect of increased FLCN expression on retinopathy. Integrating genetic association with gene expression implicated FLCN as a disease gene for diabetic retinopathy., eLife digest One of the side effects of diabetes is loss of vision from diabetic retinopathy, which is caused by injury to the light sensing tissue in the eye, the retina. Almost all individuals with diabetes develop diabetic retinopathy to some extent, and it is the leading cause of irreversible vision loss in working-age adults in the United States. How long a person has been living with diabetes, the extent of increased blood sugars and genetics all contribute to the risk and severity of diabetic retinopathy. Unfortunately, virtually no genes associated with diabetic retinopathy have yet been identified. When a gene is activated, it produces messenger molecules known as mRNA that are used by cells as instructions to produce proteins. The analysis of mRNA molecules, as well as genes themselves, can reveal the role of certain genes in disease. The studies of all genes and their associated mRNAs are respectively called genomics and transcriptomics. Genomics reveals what genes are present, while transcriptomics shows how active genes are in different cells. Skol et al. developed methods to study genomics and transcriptomics together to help discover genes that cause diabetic retinopathy. Genes involved in how cells respond to high blood sugar were first identified using cells grown in the lab. By comparing the activity of these genes in people with and without retinopathy the study identified genes associated with an increased risk of retinopathy in diabetes. In people with retinopathy, the activity of the folliculin gene (FLCN) increased more in response to high blood sugar. This was further verified with independent groups of people and using computer models to estimate the effect of different versions of the folliculin gene. The methods used here could be applied to understand complex genetics in other diseases. The results provide new understanding of the effects of diabetes. They may also help in the development of new treatments for diabetic retinopathy, which are likely to improve on the current approach of using laser surgery or injections into the eye.

The Genotype-Tissue Expression (GTEx) project dissects how genetic variation affects gene expression and splicing. Some human genetic variants affect the amount of RNA produced and the splicing of gene transcripts, crucial steps in development and maintaining a healthy individual. However, some of these changes only occur in a small number of tissues within the body. The Genotype-Tissue Expression (GTEx) project has been expanded over time, and, looking at the final data in version 8, Aguet et al. present a deep characterization of genetic associations and gene expression and splicing in 838 individuals over 49 tissues (see the Perspective by Wilson). This large study was able to characterize the details underlying many aspects of gene expression and provides a resource with which to better understand the fundamental molecular mechanisms of how genetic variants affect gene regulation and complex traits in humans. Science, this issue p. 1318; see also p. 1298 The Genotype-Tissue Expression (GTEx) project was established to characterize genetic effects on the transcriptome across human tissues and to link these regulatory mechanisms to trait and disease associations. Here, we present analyses of the version 8 data, examining 15,201 RNA-sequencing samples from 49 tissues of 838 postmortem donors. We comprehensively characterize genetic associations for gene expression and splicing in cis and trans, showing that regulatory associations are found for almost all genes, and describe the underlying molecular mechanisms and their contribution to allelic heterogeneity and pleiotropy of complex traits. Leveraging the large diversity of tissues, we provide insights into the tissue specificity of genetic effects and show that cell type composition is a key factor in understanding gene regulatory mechanisms in human tissues. We thank the donors and their families for their generous gifts of organ donation for transplantation and tissue donations for the GTEx research project; the Genomics Platform at the Broad Institute for data generation; J. Struewing for support and leadership of the GTEx project; M. Khan and C. Stolte for the illustrations in Fig. 1; and R. Do, D. Jordan, and M. Verbanck for providing GWAS pleiotropy scores. Funding: This work was supported by the Common Fund of the Office of the Director, U.S. National Institutes of Health, and by NCI, NHGRI, NHLBI, NIDA, NIMH, NIA, NIAID, and NINDS through NIH contracts HHSN261200800001E (Leidos Prime contract with NCI: A.M.S., D.E.T., N.V.R., J.A.M., L.S., M.E.B., L.Q., T.K., D.B., K.R., and A.U.), 10XS170 (NDRI: W.F.L., J.A.T., G.K., A.M., S.S., R.H., G.Wa., M.J., M.Wa., L.E.B., C.J., J.W., B.R., M.Hu., K.M., L.A.S., H.M.G., M.Mo., and L.K.B.), 10XS171 (Roswell Park Cancer Institute: B.A.F., M.T.M., E.K., B.M.G., K.D.R., and J.B.), 10X172 (Science Care Inc.), 12ST1039 (IDOX), 10ST1035 (Van Andel Institute: S.D.J., D.C.R., and D.R.V.), HHSN268201000029C (Broad Institute: F.A., G.G., K.G.A., A.V.S., X.Li., E.T., S.G., A.G., S.A., K.H.H., D.T.N., K.H., S.R.M., and J.L.N.), 5U41HG009494 (F.A., G.G., and K.G.A.), and through NIH grants R01 DA006227-17 (University of Miami Brain Bank: D.C.M. and D.A.D.), Supplement to University of Miami grant DA006227 (D.C.M. and D.A.D.), R01 MH090941 (University of Geneva), R01 MH090951 and R01 MH090937 (University of Chicago), R01 MH090936 (University of North Carolina–Chapel Hill), R01MH101814 (M.M.-A., V.W., S.B.M., R.G., E.T.D., D.G.-M., and A.V.), U01HG007593 (S.B.M.), R01MH101822 (C.D.B.), U01HG007598 (M.O. and B.E.S.), U01MH104393 (A.P.F.), extension H002371 to 5U41HG002371 (W.J.K.), as well as other funding sources: R01MH106842 (T.L., P.M., E.F., and P.J.H.), R01HL142028 (T.L., Si.Ka., and P.J.H.), R01GM122924 (T.L. and S.E.C.), R01MH107666 (H.K.I.), P30DK020595 (H.K.I.), UM1HG008901 (T.L.), R01GM124486 (T.L.), R01HG010067 (Y.Pa.), R01HG002585 (G.Wa. and M.St.), Gordon and Betty Moore Foundation GBMF 4559 (G.Wa. and M.St.), 1K99HG009916-01 (S.E.C.), R01HG006855 (Se.Ka. and R.E.H.), BIO2015-70777-P, Ministerio de Economia y Competitividad and FEDER funds (M.M.-A., V.W., R.G., and D.G.-M.), la Caixa Foundation ID 100010434 under agreement LCF/BQ/SO15/52260001 (D.G.-M.), NIH CTSA grant UL1TR002550-01 (P.M.), Marie-Skłodowska Curie fellowship H2020 Grant 706636 (S.K.-H.), R35HG010718 (E.R.G.), FPU15/03635, Ministerio de Educación, Cultura y Deporte (M.M.-A.),R01MH109905, 1R01HG010480 (A.Ba.), Searle Scholar Program (A.Ba.), R01HG008150 (S.B.M.), 5T32HG000044-22, NHGRI Institutional Training Grant in Genome Science (N.R.G.), EU IMI program (UE7-DIRECT-115317-1) (E.T.D. and A.V.), FNS funded project RNA1 (31003A_149984) (E.T.D. and A.V.), DK110919 (F.H.), F32HG009987 (F.H.), Massachusetts Lions Eye Research Fund Grant (A.R.H.), Wellcome grant WT108749/Z/15/Z (P.F.), and European Molecular Biology Laboratory (P.F. and D.Z.). Peer Reviewed "Article signat per 1 autors/es del BSC membres del THE GTEX CONSORTIUM: Marta Mele Messeguer"

Many complex human phenotypes exhibit sex-differentiated characteristics. However, the molecular mechanisms underlying these differences remain largely unknown. We generated a catalog of sex differences in gene expression and in the genetic regulation of gene expression across 44 human tissue sources surveyed by the Genotype-Tissue Expression project (GTEx, v8 release). We demonstrate that sex influences gene expression levels and cellular composition of tissue samples across the human body. A total of 37% of all genes exhibit sex-biased expression in at least one tissue. We identify cis expression quantitative trait loci (eQTLs) with sex-differentiated effects and characterize their cellular origin. By integrating sex-biased eQTLs with genome-wide association study data, we identify 58 gene-trait associations that are driven by genetic regulation of gene expression in a single sex. These findings provide an extensive characterization of sex differences in the human transcriptome and its genetic regulation. INTRODUCTION: Many complex human pheno-types, including diseases, exhibit sex-differentiated characteristics. These sex differences have been variously attributed to hormones, sex chromosomes, genotype × sex effects, differences in behavior, and differences in environmental exposures; however, their mechanisms and underlying biology remain largely unknown. The Genotype-Tissue Expression (GTEx) project provides an opportunity to investigate the prevalence and genetic mechanisms of sex differences in the human transcriptome by surveying many tissues that have not previously been characterized in this manner. RATIONALE: To characterize sex differences in the human transcriptome and its regulation, and to discover how sex and genetics interact to influence complex traits and disease, we generated a catalog of sex differences in gene expression and its genetic regulation across 44 human tissue sources surveyed by the GTEx project (v8 data release), analyzing 16,245 RNA-sequencing samples and genotypes of 838 adult individuals. We report sex differences in gene expression levels, tissue cell type composition, and cis expression quantitative trait loci (cis-eQTLs). To assess their impact, we integrated these results with gene function, transcription factor binding annotation, and genome-wide association study (GWAS) summary statistics of 87 GWASs. RESULTS: Sex effects on gene expression are ubiquitous (13,294 sex-biased genes across all tissues). However, these effects are small and largely tissue-specific. Genes with sex-differentiated expression are not primarily driven by tissue-specific gene expression and are involved in a diverse set of biological functions, such as drug and hormone response, embryonic development and tissue morphogenesis, fertilization, sexual reproduction and spermatogenesis, fat metabolism, cancer, and immune response. Whereas X-linked genes with higher expression in females suggest candidates for escape from X-chromosome inactivation, sex-biased expression of autosomal genes suggests hormone-related transcription factor regulation and a role for additional transcription factors, as well as sex-differentiated distribution of epigenetic marks, particularly histone H3 Lys(27) trimethylation (H3K27me3). Sex differences in the genetic regulation of gene expression are much less common (369 sex-biased eQTLs across all tissues) and are highly tissue-specific. We identified 58 gene-trait associations driven by genetic regulation of gene expression in a single sex. These include loci where sex-differentiated cell type abundances mediate genotype-phenotype associations, as well as loci where sex may play a more direct role in the underlying molecular mechanism of the association. For example, we identified a female-specific eQTL in liver for the hexokinase HKDC1 that influences glucose metabolism in pregnant females, which is subsequently reflected in the birth weight of the offspring. CONCLUSION: By integrating sex-aware analyses of GTEx data with gene function and transcription factor binding annotations, we describe tissue-specific and tissue-shared drivers and mechanisms contributing to sex differences in the human transcriptome and eQTLs. We discovered multiple sex-differentiated genetic effects on gene expression that colocalize with complex trait genetic associations, thereby facilitating the mechanistic interpretation of GWAS signals. Because the causative tissue is unknown for many phenotypes, analysis of the diverse GTEx tissue collection can serve as a powerful resource for investigations into the basis of sex-biased traits. This work provides an extensive characterization of sex differences in the human transcriptome and its genetic regulation. ■

The goal of the study was to identify genes whose aberrant expression can contribute to diabetic retinopathy. We determined differential gene expression in response to high glucose in lymphoblastoid cell lines derived from matched individuals with type 1 diabetes (T1D) with and without retinopathy. Those genes exhibiting the largest difference in glucose response between individuals with diabetes with and without retinopathy were assessed for association to diabetic retinopathy utilizing genotype data from a genome-wide association study meta-analysis. All genetic variants associated with gene expression (expression Quantitative Trait Loci, eQTLs) of the glucose response genes were tested for association with diabetic retinopathy. We detected an enrichment of the eQTLs from the glucose response genes among small association p-values and identified folliculin (FLCN) as a susceptibility gene for diabetic retinopathy. We show that expression of FLCN in response to glucose was greater in individuals with diabetic retinopathy compared to individuals with diabetes without retinopathy. Three large, independent cohorts of individuals with diabetes revealed an association of FLCN eQTLs to diabetic retinopathy. Mendelian randomization further confirmed a direct positive effect of increased FLCN expression on retinopathy in individuals with diabetes. Together, our studies integrating genetic association and gene expression implicate FLCN as a disease gene for diabetic retinopathy.

We investigated a Spanish and Catalan family in which multiple cancer types tracked across three generations, but for which no genetic etiology had been identified. Whole exome sequencing of germline DNA from multiple affected family members was performed to identify candidate variants to explain this occurrence of familial cancer. We discovered in all cancer-affected family members a single rare heterozygous germline variant (I654V, rs1801201) in ERBB2/HER2 that is located in a transmembrane glycine zipper motif critical for ERBB2-mediated signaling and in complete linkage disequilibrium (D’=1) with a common polymorphism (I655V, rs1136201) previously reported in some populations as associated with cancer risk. Because multiple cancer types occurred in this family, we tested both the I654V and the I655V variants for association with cancer across multiple tumor types in 6,371 cases of Northern European ancestry drawn from The Cancer Genome Atlas and 6,647 controls, and found that the rare variant (I654V) was significantly associated with an increased risk for cancer (OR=1.40, p=0.021, 95% CI=1.05–1.89). Functional assays performed in HEK 293T cells revealed that both the I655V single mutant (SM) and the I654V;I655V double mutant (DM) stabilized ERBB2 protein and activated ERBB2 signaling, with the DM activating ERBB2 significantly more than the SM alone. Thus, our results suggest a model whereby heritable genetic variation in the transmembrane domain activating ERBB2 signaling is associated with both sporadic and familial cancer risk, with increased ERBB2 stabilization and activation associated with increased cancer risk.

Allele expression (AE) analysis robustly measures cis-regulatory effects. Here, we present and demonstrate the utility of a vast AE resource generated from the GTEx v8 release, containing 15,253 samples spanning 54 human tissues for a total of 431 million measurements of AE at the SNP level and 153 million measurements at the haplotype level. In addition, we develop an extension of our tool phASER that allows effect sizes of cis-regulatory variants to be estimated using haplotype-level AE data. This AE resource is the largest to date, and we are able to make haplotype-level data publicly available. We anticipate that the availability of this resource will enable future studies of regulatory variation across human tissues.

Background We previously demonstrated a gene-by-prenatal-environment interaction whereby the monoamine oxidase A gene (MAOA) modified the impact of prenatal tobacco exposure (PTE) on adolescent disruptive behavior (DB), with the MAOA risk genotype varying by sex. We extend this work by examining whether this mechanism is evident with another common adversity, prenatal stress exposure (PSE), and whether sex differences are present earlier in development in closer proximity to exposure. Methods Participants were 281 mothers and their 285 children derived from a prenatal cohort with in-depth prospective measures of PSE and PTE. We assessed DB at age 5 via dimensional developmentally-sensitive measurement. Analyses were stratified by sex based on prior evidence for sex differences. Results Concurrent stress exposure predicted DB in children (β = 0.310, p = 0.001), while main effects of prenatal exposures were seen only in boys. We found a three-way interaction of MAOA × PSE × sex on DB (β = 0.813, p = 0.022). Boys with MAOA-H had more DB as a function of PSE, controlling for PTE (β = 0.774, p = 0.015), and as a function of PTE, controlling for PSE (β = 0.362, p = 0.037). Boys with MAOA-L did not show this susceptibility. MAOA did not interact with PSE (β = − 0.133, p = 0.561) nor PTE (β = − 0.144; p = 0.505) in predicting DB in girls. Examination of gene-environment correlation (rGE) showed a correlation between paternal MAOA-L and daughters' concurrent stress exposure (r = − 0.240, p = 0.013). Discussion Findings underscore complex mechanisms linking genetic susceptibility and early adverse exposures. Replication in larger cohorts followed from the pregnancy through adolescence is suggested to elucidate mechanisms that appear to have varying developmental expression.

Molecular markers of sorafenib efficacy in patients with metastatic renal cell carcinoma (mRCC) are not available. The purpose of this study was to discover genetic markers of survival in patients with mRCC treated with sorafenib. Germline variants from 56 genes were genotyped in 295 patients with mRCC. Variant-overall survival (OS) associations were tested in multivariate regression models. Mechanistic studies were conducted to validate clinical associations. VEGFA rs1885657, ITGAV rs3816375, and WWOX rs8047917 (sorafenib arm), and FLT4 rs307826 and VEGFA rs3024987 (sorafenib and placebo arms combined) were associated with shorter OS. FLT4 rs307826 increased VEGFR-3 phosphorylation, membrane trafficking, and receptor activation. VEGFA rs1885657 and rs58159269 increased transcriptional activity of the constructs containing these variants in endothelial and RCC cell lines, and VEGFA rs58159269 increased endothelial cell proliferation and tube formation. FLT4 rs307826 and VEGFA rs58159269 led to reduced sorafenib cytotoxicity. Genetic variation in VEGFA and FLT4 could affect survival in sorafenib-treated patients with mRCC. These markers should be examined in additional malignancies treated with sorafenib and in other angiogenesis inhibitors used in mRCC. Significance: Clinical and mechanistic data identify germline genetic variants in VEGFA and FLT4 as markers of survival in patients with metastatic renal cell carcinoma.

Long non-coding RNA (lncRNA) genes have well-established and important impacts on molecular and cellular functions. However, among the thousands of lncRNA genes, it is still a major challenge to identify the subset with disease or trait relevance. To systematically characterize these lncRNA genes, we used Genotype Tissue Expression (GTEx) project v8 genetic and multi-tissue transcriptomic data to profile the expression, genetic regulation, cellular contexts, and trait associations of 14,100 lncRNA genes across 49 tissues for 101 distinct complex genetic traits. Using these approaches, we identified 1,432 lncRNA gene-trait associations, 800 of which were not explained by stronger effects of neighboring protein-coding genes. This included associations between lncRNA quantitative trait loci and inflammatory bowel disease, type 1 and type 2 diabetes, and coronary artery disease, as well as rare variant associations to body mass index.

Determining protein levels in each tissue and how they compare with RNA levels is important for understanding human biology and disease as well as regulatory processes that control protein levels. We quantified the relative protein levels from 12,627 genes across 32 normal human tissue types prepared by the GTEx project. Known and new tissue specific or enriched proteins (5,499) were identified and compared to transcriptome data. Many ubiquitous transcripts are found to encode highly tissue specific proteins. Discordance in the sites of RNA expression and protein detection also revealed potential sites of synthesis and action of protein signaling molecules. Overall, these results provide an extraordinary resource, and demonstrate that understanding protein levels can provide insights into metabolism, regulation, secretome, and human diseases.SummaryQuantitative proteome study of 32 human tissues and integrated analysis with transcriptome data revealed that understanding protein levels could provide in-depth knowledge to post transcriptional or translational regulations, human metabolism, secretome, and diseases.

Molecular markers of sorafenib efficacy in metastatic renal cell carcinoma (mRCC) patients are not available. The purpose of this study was to discover genetic markers of survival in mRCC patients treated with sorafenib. Germline variants from 56 genes were genotyped in 295 mRCC patients. Variant-overall survival (OS) associations were tested in multivariate regression models. Mechanistic studies were conducted to validate clinical associations. VEGFA rs1885657, ITGAV rs3816375, and WWOX rs8047917 (sorafenib arm), and FLT4 rs307826 and VEGFA rs3024987 (sorafenib and placebo arms combined) were associated with shorter OS. FLT4 rs307826 increased VEGFR-3 phosphorylation, membrane trafficking and receptor activation. VEGFA rs1885657 and rs58159269 increased transcriptional activity of the constructs containing these variants in endothelial and RCC cell lines, and VEGFA rs58159269 increased endothelial cell proliferation and tube formation. FLT4 rs307826 and VEGFA rs58159269 led to reduced sorafenib cytotoxicity. Genetic variation in VEGFA and FLT4 could affect survival in sorafenib-treated mRCC patients. These markers should be examined in additional malignancies treated with sorafenib, and in other angiogenesis inhibitors used in mRCC.

Gene co-expression networks capture biologically important patterns in gene expression data, enabling functional analyses of genes, discovery of biomarkers, and interpretation of regulatory genetic variants. Most network analyses to date have been limited to assessing correlation between total gene expression levels in a single or small sets of tissues. Here, we have reconstructed networks that capture a much more complete set of regulatory relationships, specifically including regulation of relative isoform abundance and splicing, and tissue-specific connections unique to each of a diverse set of tissues. Using the Genotype-Tissue Expression (GTEx) project v6 RNA-sequencing data across 44 tissues in 449 individuals, we evaluated shared and tissue-specific network relationships. First, we developed a framework called Transcriptome Wide Networks (TWNs) for combining total expression and relative isoform levels into a single sparse network, capturing the complex interplay between the regulation of splicing and transcription. We built TWNs for sixteen tissues, and found that hubs with isoform node neighbors in these networks were strongly enriched for splicing and RNA binding genes, demonstrating their utility in unraveling regulation of splicing in the human transcriptome, and providing a set of candidate shared and tissue-specific regulatory hub genes. Next, we used a Bayesian biclustering model that identifies network edges between genes with co-expression in a single tissue to reconstruct tissue-specific networks (TSNs) for 27 distinct GTEx tissues and for four subsets of related tissues. Using both TWNs and TSNs, we characterized gene co-expression patterns shared across tissues. Finally, we found genetic variants associated with multiple neighboring nodes in our networks, supporting the estimated network structures and identifying 33 genetic variants with distant regulatory impact on transcription and splicing. Our networks provide an improved understanding of the complex relationships between genes in the human transcriptome, including tissue-specificity of gene co-expression, regulation of splicing, and the coordinated impact of genetic variation on transcription.

The incidence and mortality rates of B-cell acute lymphoblastic leukemia (B-ALL) differ by age and sex. To determine if inherited genetic susceptibility contributes to these differences we performed 2 genome-wide association studies (GWAS) by age, sex, and subtype and subsequent meta-analyses. The GWAS included 446 B-ALL cases, and 3027 healthy unrelated blood and marrow transplant (BMT) donors as controls from the Determining the Influence of Susceptibility Conveying Variants Related to One-Year Mortality after BMT (DISCOVeRY-BMT) study. We identified 1 novel variant, rs189434316, significantly associated with odds of normal cytogenetic B-ALL (odds ratio from meta-analysis [ORmeta] = 3.7; 95% confidence interval [CI], 2.5, 6.2; P value from meta-analysis [Pmeta] = 6.0 × 10-9). The previously reported pediatric B-ALL GWAS variant, rs11980379 (IKZF1), replicated in B-ALL pediatric patients (ORmeta = 2.3; 95% CI, 1.5, 3.7; Pmeta = 1.0 × 10-9), with evidence of heterogeneity (P = .02) between males and females. Sex differences in single-nucleotide polymorphism effect were seen in those >15 years (OR = 1.7; 95% CI, 1.4, 2.2, PMales = 6.38 × 10-6/OR = 1.1; 95% CI, 0.8, 1.5; PFemales = .6) but not ≤15 years (OR = 2.3; 95% CI, 1.4, 3.8; PMales = .0007/OR = 1.9; 95% CI, 1.2, 3.2; PFemales = .007). The latter association replicated in independent pediatric B-ALL cohorts. A previously identified adolescent and young-adult onset ALL-associated variant in GATA3 is associated with B-ALL risk in those >40 years. Our findings provide more evidence of the influence of genetics on B-ALL age of onset and we have shown the first evidence that IKZF1 associations with B-ALL may be sex and age specific.

PURPOSE: Although there are considerable data on mechanisms of radiation-induced apoptosis in vitro and in animal models, little is known about functional variation in these pathways in humans. We sought to develop a tractable system to evaluate this. MATERIALS AND METHODS: Peripheral blood mononuclear cells were isolated from 90 healthy volunteers, divided into two aliquots, one irradiated with a 5 Gy dose and the other sham-treated (0 Gy), and assessed for damage-induced apoptosis after 24 hours. To investigate reproducibility, 10 individuals spanning the entire radiation-induced apoptotic range were tested three times each, with 3-6 months between replicates. RESULTS: We observed surprising heterogeneity in apoptosis among individuals, ranging from 21-62%. Biological replicates from a single individual, however, were completely concordant, suggesting the variability observed across individuals is not the result of stochastic or short-term effects. We found significantly higher radiation-induced apoptosis in males than in females (Mean: 41.0% vs. 30.7%; p < 3.5 × 10(-7)). Moreover, advancing age was associated with decreasing radiation-induced apoptosis in males (p = 0.01) but not females (p = 0.82). CONCLUSIONS: Our results provide evidence that the function of cellular pathways crucial for stress-induced apoptosis varies by sex and could decline with age in humans.

Genome-wide association studies (GWAS) in colorectal cancer (CRC) identified five regions near transforming growth factor β-related genes BMP4, GREM1, CDH1, SMAD7 and RPHN2. The true risk alleles remain to be identified in these regions, and their role in CRC risk in non-European populations has been understudied. Our previous work noted significant genetic heterogeneity between African Americans (AAs) and European Americans (EAs) for single nucleotide polymorphisms (SNPs) identified in GWAS. We hypothesized that associations may not have been replicated in AAs due to differential or independent genetic structures. In order to test this hypothesis, we genotyped 195 tagging SNPs across these five gene regions in 1194 CRC cases (795 AAs and 399 EAs) and 1352 controls (985 AAs and 367 EAs). Imputation was performed, and association testing of genotyped and imputed SNPs included ancestry, age and sex as covariates. In two of the five genes originally associated with CRC, we found evidence for association in AAs including rs1862748 in CDH1 (OR(Add) = 0.82, P = 0.02) and in GREM1 the SNPs rs10318 (OR(Rec) = 60.1, P = 0.01), rs11632715 (OR(Rec) = 2.36; P = 0.004) and rs12902616 (OR(Rec) = 1.28, P = 0.005), the latter which is in linkage disequilibrium with the previously identified SNP rs4779584. Testing more broadly for associations in these gene regions in AAs, we noted three statistically significant association peaks in GREM1 and RHPN2 that were not identified in EAs. We conclude that some CRC risk alleles are shared between EAs and AAs and others are population specific.

To investigate genetic predispositions for MYCN-amplified neuroblastoma, we performed a meta-analysis of three genome-wide association studies totaling 615 MYCN-amplified high-risk neuroblastoma cases and 1869 MYCN-nonamplified non-high-risk neuroblastoma cases as controls using a fixed-effects model with inverse variance weighting. All statistical tests were two-sided. We identified a novel locus at 3p21.31 indexed by the single nucleotide polymorphism (SNP) rs80059929 (odds ratio [OR] = 2.95, 95% confidence interval [CI] = 2.17 to 4.02, Pmeta = 6.47 × 10-12) associated with MYCN-amplified neuroblastoma, which was replicated in 127 MYCN-amplified cases and 254 non-high-risk controls (OR = 2.30, 95% CI = 1.12 to 4.69, Preplication = .02). To confirm this signal is exclusive to MYCN-amplified tumors, we performed a second meta-analysis comparing 728 MYCN-nonamplified high-risk patients to identical controls. rs80059929 was not statistically significant in MYCN-nonamplified high-risk patients (OR = 1.24, 95% CI = 0.90 to 1.71, Pmeta = .19). SNP rs80059929 is within intron 16 in the KIF15 gene. Additionally, the previously reported LMO1 neuroblastoma risk locus was statistically significant only in patients with MYCN-nonamplified high-risk tumors (OR = 0.63, 95% CI = 0.53 to 0.75, Pmeta = 1.51 × 10-8; Pmeta = .95). Our results indicate that common genetic variation predisposes to different neuroblastoma genotypes, including the likelihood of somatic MYCN-amplification.

More than 12 % of women will be diagnosed with breast cancer in their lifetime. Although there have been tremendous advances in elucidating genetic risk factors underlying both familial and sporadic breast cancer, much of the genetic contribution to breast cancer etiology remains unknown. The discovery of BRCA1 and BRCA2 over 20 years ago remains the seminal event in the field and has paved the way for the discovery of other high-penetrance susceptibility genes by linkage analysis. The advent of genome-wide association studies made possible the next wave of discoveries, in which over 80 low-penetrance and moderate-penetrance variants were identified. Although these studies were highly successful at discovering variants associated with both familial and sporadic breast cancer, the variants identified to date explain only 50 % of the heritability of breast cancer. In this review, we look back at the investigative strategies that have led to our current understanding of breast cancer genetics, consider the challenges of performing association studies in heterogeneous complex diseases such as breast cancer, and look ahead toward the types of study designs that may lead to the identification of the genetic variation accounting for the remaining missing heritability.

Paediatric B-cell precursor acute lymphoblastic leukaemia (BCP-ALL) is the most common cancer of childhood, yet little is known about BCP-ALL predisposition. In this study, in 2,187 cases of European ancestry and 5,543 controls, we discover and replicate a locus indexed by rs77728904 at 9p21.3 associated with BCP-ALL susceptibility (Pcombined=3.32 × 10−15, OR=1.72) and independent from rs3731217, the previously reported ALL-associated variant in this region. Of correlated SNPs tagged by this locus, only rs662463 is significant in African Americans, suggesting it is a plausible causative variant. Functional analysis shows that rs662463 is a cis-eQTL for CDKN2B, with the risk allele associated with lower expression, and suggests that rs662463 influences BCP-ALL risk by regulating CDKN2B expression through CEBPB signalling. Functional analysis of rs3731217 suggests it is associated with BCP-ALL by acting within a splicing regulatory element determining CDKN2A exon 3 usage (P=0.01). These findings provide new insights into the critical role of the CDKN2 locus in BCP-ALL aetiology., A risk variant located at 9p21.3 is associated with cancer risk in pediatric B-cell precursor acute lymphoblastic leukaemia. Here, the authors show that this variant affects the gene expression of the tumour suppressor gene Cdkn2b.

Nodular sclerosing Hodgkin lymphoma (NSHL) is a distinct, highly heritable Hodgkin lymphoma subtype. We undertook a genome-wide meta-analysis of 393 European-origin adolescent/young adult NSHL patients and 3315 controls using the Illumina Human610-Quad Beadchip and Affymetrix Genome-Wide Human SNP Array 6.0. We identified 3 single nucleotide polymorphisms (SNPs) on chromosome 6p21.32 that were significantly associated with NSHL risk: rs9268542 (P = 5.35 × 10−10), rs204999 (P = 1.44 × 10−9), and rs2858870 (P = 1.69 × 10−8). We also confirmed a previously reported association in the same region, rs6903608 (P = 3.52 × 10−10). rs204999 and rs2858870 were weakly correlated (r2 = 0.257), and the remaining pairs of SNPs were not correlated (r2 < 0.1). In an independent set of 113 NSHL cases and 214 controls, 2 SNPs were significantly associated with NSHL and a third showed a comparable odds ratio (OR). These SNPs are found on 2 haplotypes associated with NSHL risk (rs204999-rs9268528-rs9268542-rs6903608-rs2858870; AGGCT, OR = 1.7, P = 1.71 × 10−6; GAATC, OR = 0.4, P = 1.16 × 10−4). All individuals with the GAATC haplotype also carried the HLA class II DRB1*0701 allele. In a separate analysis, the DRB1*0701 allele was associated with a decreased risk of NSHL (OR = 0.5, 95% confidence interval = 0.4, 0.7). These data support the importance of the HLA class II region in NSHL etiology.

Purpose: To discover genetic determinants of overall survival (OS) in metastatic renal cell carcinoma (mRCC) patients treated with sorafenib. Patients and Methods: Information for 11,117 germline DNA variants from 56 genes was obtained, using both Illumina GoldenGate genotyping and imputation, from 295 patients originally enrolled in the phase III TARGET trial (155 sorafenib and 140 placebo patients). OS was this genotyping study's primary endpoint, and associations between variants and OS were tested using multivariate Cox regression. Associations were considered significant if P Results: Five intragenic variants significantly associated with OS. Among patients in the sorafenib arm, VEGFA rs1885657 (HR=17.3, 95% CI 5.7-52.7; P=1.4x10-4), ITGAV rs3816375 (HR=5.9, 95% CI 2.1-16.4; P=4.9x10-4), and WWOX rs8047917 (HR=4.1, 95% CI 1.9-8.8; P=3.3x10-4) associated with shorter OS. Among patients from both arms, FLT4 rs307826 (HR=13.8, 95% CI 3.0-62.6; P=1.2x10-4) and VEGFA rs3024987 (HR=3.0, 95% CI 1.7-5.4; P=8.3x10-5) associated with shorter OS. No significant associations were discovered in the placebo arm. Bioinformatic evidence was used to prioritize rs1885657, rs58159269 (in complete linkage disequilibrium with rs1885657), and rs307826. rs1885657 and rs58159269 increased luciferase activity in human endothelial (TIME, LPEC) and mRCC (Caki-1) cells, and rs58159269 increased cell proliferation/reduced sorafenib cytotoxicity, and reduced endothelial tube formation in TIME cells. rs307826 increased VEGFR-3 phosphorylation, and led to post-translational VEGFR-3 processing differences in HUVECs. Conclusions: Regulatory VEGFA and FLT4 variants associated with OS in mRCC patients induce cellular alterations that affect angiogenesis and sorafenib activity. These are novel and potentially impactful clinical biomarkers for mRCC and sorafenib, and their relevance to other multikinase inhibitors used in mRCC should be tested. Citation Format: Daniel James Crona, Andrew D. Skol, Veli-Matti Leppänen, Dylan M. Glubb, Amy S. Etheridge, Eleanor Hilliard, Carol Peña, Yuri K. Peterson, Nancy Klauber-DeMore, Kari K. Alitalo, Federico Innocenti. Clinical biomarkers of survival in renal cell carcinoma patients treated with sorafenib [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2018; 2018 Apr 14-18; Chicago, IL. Philadelphia (PA): AACR; Cancer Res 2018;78(13 Suppl):Abstract nr 1625.

Model organisms offer many advantages for the genetic analysis of complex traits. However, identification of specific genes is often hampered by a lack of recombination between the genomes of inbred progenitors. Recently, genome-wide association studies (GWAS) in humans have offered gene-level mapping resolution that is possible because of the large number of accumulated recombinations among unrelated human subjects. To obtain analogous improvements in mapping resolution in mice, we used a 34th generation advanced intercross line (AIL) derived from two inbred strains (SM/J and LG/J). We used simulations to show that familial relationships among subjects must be accounted for when analyzing these data; we then used a mixed model that included polygenic effects to address this problem in our own analysis. Using a combination of F2 and AIL mice derived from the same inbred progenitors, we identified genome-wide significant, subcentimorgan loci that were associated with methamphetamine sensitivity, (e.g., chromosome 18; LOD = 10.5) and non-drug-induced locomotor activity (e.g., chromosome 8; LOD = 18.9). The 2-LOD support interval for the former locus contains no known genes while the latter contains only one gene (Csmd1). This approach is broadly applicable in terms of phenotypes and model organisms and allows GWAS to be performed in multigenerational crosses between and among inbred strains where familial relatedness is often unavoidable.

Therapy-related acute myeloid leukemia (t-AML) is a rare but fatal complication of cytotoxic therapy. Whereas sporadic cancer results from interactions between complex exposures and low-penetrance alleles, t-AML results from an acute exposure to a limited number of potent genotoxins. Consequently, we hypothesized that the effect sizes of variants associated with t-AML would be greater than in sporadic cancer, and, therefore, that these variants could be detected even in a modest-sized cohort. To test this, we undertook an association study in 80 cases and 150 controls using Affymetrix Mapping 10K arrays. Even at nominal significance thresholds, we found a significant excess of associations over chance; for example, although 6 associations were expected at P less than .001, we found 15 (Penrich = .002). To replicate our findings, we genotyped the 10 most significantly associated single nucleotide polymorphisms (SNPs) in an independent t-AML cohort (n = 70) and obtained evidence of association with t-AML for 3 SNPs in the subset of patients with loss of chromosomes 5 or 7 or both, acquired abnormalities associated with prior exposure to alkylator chemotherapy. Thus, we conclude that the effect of genetic factors contributing to cancer risk is potentiated and more readily discernable in t-AML compared with sporadic cancer.

Genetic susceptibility to antisocial behavior may increase fetal sensitivity to prenatal exposure to cigarette smoke. Testing putative gene x exposure mechanisms requires precise measurement of exposure and outcomes. We tested whether a functional polymorphism in the gene encoding the enzyme monoamine oxidase A (MAOA) interacts with exposure to predict pathways to adolescent antisocial behavior. We assessed both clinical and information-processing outcomes. One hundred seventy-six adolescents and their mothers participated in a follow-up of a pregnancy cohort with well-characterized exposure. A sex-specific pattern of gene x exposure interaction was detected. Exposed boys with the low-activity MAOA 5' uVNTR (untranslated region variable number of tandem repeats) genotype were at increased risk for conduct disorder (CD) symptoms. In contrast, exposed girls with the high-activity MAOA uVNTR genotype were at increased risk for both CD symptoms and hostile attribution bias on a face-processing task. There was no evidence of a gene-environment correlation (rGE). Findings suggest that the MAOA uVNTR genotype, prenatal exposure to cigarettes and sex interact to predict antisocial behavior and related information-processing patterns. Future research to replicate and extend these findings should focus on elucidating how gene x exposure interactions may shape behavior through associated changes in brain function.

More than 120 published reports have described associations between single nucleotide polymorphisms (SNPs) and type 2 diabetes. However, multiple studies of the same variant have often been discordant. From a literature search, we identified previously reported type 2 diabetes–associated SNPs. We initially genotyped 134 SNPs on 786 index case subjects from type 2 diabetes families and 617 control subjects with normal glucose tolerance from Finland and excluded from analysis 20 SNPs in strong linkage disequilibrium (r2 > 0.8) with another typed SNP. Of the 114 SNPs examined, we followed up the 20 most significant SNPs (P < 0.10) on an additional 384 case subjects and 366 control subjects from a population-based study in Finland. In the combined data, we replicated association (P < 0.05) for 12 SNPs: PPARG Pro12Ala and His447, KCNJ11 Glu23Lys and rs5210, TNF −857, SLC2A2 Ile110Thr, HNF1A/TCF1 rs2701175 and GE117881_360, PCK1 −232, NEUROD1 Thr45Ala, IL6 −598, and ENPP1 Lys121Gln. The replication of 12 SNPs of 114 tested was significantly greater than expected by chance under the null hypothesis of no association (P = 0.012). We observed that SNPs from genes that had three or more previous reports of association were significantly more likely to be replicated in our sample (P = 0.03), although we also replicated 4 of 58 SNPs from genes that had only one previous report of association.

BACKGROUND Rare variants (

Background: Little is known about genetic factors associated with nasopharyngeal carcinoma (NPC). To gain insight into NPC etiology, we performed whole exome sequencing on germline and tumor DNA from three closely related family members with NPC. Methods: The family was ascertained through the Pediatric Familial Cancer Clinic at The University of Chicago (Chicago, IL). The diagnosis of NPC was confirmed pathologically for each individual. For each sample sequenced, 97.3% of the exome was covered at 5×, with an average depth of 44×. Candidate germline and somatic variants associated with NPC were identified and prioritized using a custom pipeline. Results: We discovered 72 rare deleterious germline variants in 56 genes shared by all three individuals. Of these, only three are in previously identified NPC-associated genes, all of which are located within MLL3, a gene known to be somatically altered in NPC. One variant introduces an early stop codon in MLL3, which predicts complete loss-of-function. Tumor DNA analysis revealed somatic mutations and Epstein–Barr virus (EBV) integration events; none, however, were shared among all three individuals. Conclusions: These data suggest that inherited mutations in MLL3 may have predisposed these three individuals from a single family to develop NPC, and may cooperate with individually acquired somatic mutations or EBV integration events in NPC etiology. Impact: Our finding is the first instance of a plausible candidate high penetrance inherited mutation predisposing to NPC. Cancer Epidemiol Biomarkers Prev; 24(8); 1222–8. ©2015 AACR.

Wei and Hemmings [2000: Nat Genet 25:376-377], using 80 British parent-offspring trios, identified a number of NOTCH4 variants and haplotypes that showed statistically significant evidence of association to schizophrenia. Specifically, the 10 repeat allele of a (CTG)(n) marker and the 8 repeat allele of a (TAA)(n) marker demonstrated excess transmission to affected individuals; SNP21 and haplotypes SNP2-(CTG)(n) and SNP12-SNP2-(CTG)(n) also showed significant associations. In an attempt to replicate these findings, we tested for linkage and association between the same five markers used by Wei and Hemmings in 166 families collected from a multi-center study conducted by the Department of Veterans Affairs (DVA) Cooperative Study Program (CSP). The families include 392 affected subjects (schizophrenia or schizoaffective disorder, depressed) and 216 affected sibling pairs. The families represent a mix of European Americans (n = 62, 37%), African Americans (n = 60, 36%), and racially mixed or other races (n = 44, 27%). We identified moderate evidence for linkage in the pooled race sample (LOD = 1.25) and found excess transmission of the 8 (P = 0.06) and 13 (P = 0.04) repeat alleles of the (TAA)(n) marker to African American schizophrenic subjects. The 8 and 13 repeat alleles were previously identified to be positively associated with schizophrenia by Wei and Hemmings [2000: Nat Genet 25:376-377] and Sklar et al. [2001: Nat Genet 28:126-128], respectively.

Gene-mapping studies routinely rely on checking for Mendelian transmission of marker alleles in a pedigree, as a means of screening for genotyping errors and mutations, with the implicit assumption that, if a pedigree is consistent with Mendel's laws of inheritance, then there are no genotyping errors. However, the occurrence of inheritance inconsistencies alone is an inadequate measure of the number of genotyping errors, since the rate of occurrence depends on the number and relationships of genotyped pedigree members, the type of errors, and the distribution of marker-allele frequencies. In this article, we calculate the expected probability of detection of a genotyping error or mutation as an inheritance inconsistency in nuclear-family data, as a function of both the number of genotyped parents and offspring and the marker-allele frequency distribution. Through computer simulation, we explore the sensitivity of our analytic calculations to the underlying error model. Under a random-allele-error model, we find that detection rates are 51%-77% for multiallelic markers and 13%-75% for biallelic markers; detection rates are generally lower when the error occurs in a parent than in an offspring, unless a large number of offspring are genotyped. Errors are especially difficult to detect for biallelic markers with equally frequent alleles, even when both parents are genotyped; in this case, the maximum detection rate is 34% for four-person nuclear families. Error detection in families in which parents are not genotyped is limited, even with multiallelic markers. Given these results, we recommend that additional error checking (e.g., on the basis of multipoint analysis) be performed, beyond routine checking for Mendelian consistency. Furthermore, our results permit assessment of the plausibility of an observed number of inheritance inconsistencies for a family, allowing the detection of likely pedigree-rather than genotyping-errors in the early stages of a genome scan. Such early assessments are valuable in either the targeting of families for resampling or discontinued genotyping.

Melanoma is the deadliest form of skin cancer and presents a significant health care burden in many countries. In addition to ultraviolet radiation in sunlight, the main causal factor for melanoma, genetic factors also play an important role in melanoma susceptibility. Although genome-wide association studies have identified many single nucleotide polymorphisms associated with melanoma, little is known about the proportion of disease risk attributable to these loci and their distribution throughout the genome. Here, we investigated the genetic architecture of melanoma in 1,888 cases and 990 controls of European non-Hispanic ancestry. We estimated the overall narrow-sense heritability of melanoma to be 0.18 (P < 0.03), indicating that genetics contributes significantly to the risk of sporadically-occurring melanoma. We then demonstrated that only a small proportion of this risk is attributable to known risk variants, suggesting that much remains unknown of the role of genetics in melanoma. To investigate further the genetic architecture of melanoma, we partitioned the heritability by chromosome, minor allele frequency, and functional annotations. We showed that common genetic variation contributes significantly to melanoma risk, with a risk model defined by a handful of genomic regions rather than many risk loci distributed throughout the genome. We also demonstrated that variants affecting gene expression in skin account for a significant proportion of the heritability, and are enriched among melanoma risk loci. Finally, by incorporating skin color into our analyses, we observed both a shift in significance for melanoma-associated loci and an enrichment of expression quantitative trait loci among melanoma susceptibility variants. These findings suggest that skin color may be an important modifier of melanoma risk. We speculate that incorporating skin color and other non-genetic factors into genetic studies may allow for an improved understanding of melanoma susceptibility and guide future investigations to identify melanoma risk genes.

Previous studies have reported genetic linkage evidence for a schizophrenia gene on chromosome 15q. Here, chromosome 15 was examined by genetic linkage analysis using 166 schizophrenia families, each with two or more affected subjects. The families, assembled from multiple centers by the Department of Veterans Affairs Cooperative Study Program, consisted of 392 sampled affected subjects and 216 affected sibling pairs. By DSM-III-R criteria, 360 subjects (91.8%) had a diagnosis of schizophrenia and 32 (8.2%) were classified as schizo-affective disorder, depressed. Participating families had diverse ethnic backgrounds. The largest single group were northern European American families (n = 62, 37%), but a substantial proportion was African American kindreds (n = 60, 36%). The chromosome 15 markers tested were spaced at intervals of approximately 10 cM over the entire chromosome and 2-5 cM for the region surrounding the alpha-7 nicotinic cholinergic receptor subunit gene (CHRNA7). These markers were genotyped and the data analyzed using semiparametric affecteds-only linkage analysis. In the European American families, there was a maximum Z-score of 1.65 between markers D15S165 and D15S1010. These markers are within 1 cM from CHRNA-7, the site previously implicated in schizophrenia. However, there was no evidence for linkage to this region in the African America kindreds.