Ocular manifestations of CHARGE syndrome in a pediatric cohort with genotype/phenotype analysis.

CHARGE syndrome is a rare multi-system condition associated with CHD7 variants. However, ocular manifestations and particularly ophthalmic genotype-phenotype associations, are not well-studied. This study evaluated ocular manifestations and genotype-phenotype associations in pediatric patients with CHARGE syndrome. A retrospective chart review included pediatric patients under 20 years-old with clinical diagnosis of CHARGE syndrome and documented ophthalmic examination. Demographics, genetic testing, and ocular findings were collected. Comprehensive literature review enhanced the genotype-phenotype analysis. Forty-two patients (20 male) underwent eye examination at an average age of 9.45 ± 6.52 years-old. Thirty-nine (93%) had ophthalmic manifestations in at least one eye. Optic nerve/chorioretinal colobomas were most common (38 patients), followed by microphthalmia (13), cataract (6), and iris colobomas (4). Extraocular findings included strabismus (32 patients), nasolacrimal duct obstructions (11, 5 with punctal agenesis), and cranial nerve VII palsy (10). Genotype-phenotype analyses (27 patients) showed variability in ocular phenotypes without association to location or variant types. Splicing (10 patients) and frameshift (10) variants were most prevalent. Patients with CHARGE syndrome may present with a myriad of ophthalmic manifestations. There is limited data regarding genotype-phenotype correlations and additional studies are needed.

Inherited cancer predisposing mutations in patients with therapy-related myeloid neoplasms.

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.

PIK3R5 genetic predictors of hypertension induced by VEGF-pathway inhibitors.

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).

The genetics of breast cancer risk in the post-genome era: thoughts on study design to move past BRCA and towards clinical relevance.

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.

Enrichment of inflammatory bowel disease and colorectal cancer risk variants in colon expression quantitative trait loci.

BACKGROUND: Genome-wide association studies (GWAS) have identified single nucleotide polymorphisms (SNPs) associated with diseases of the colon including inflammatory bowel diseases (IBD) and colorectal cancer (CRC). However, the functional role of many of these SNPs is largely unknown and tissue-specific resources are lacking. Expression quantitative trait loci (eQTL) mapping identifies target genes of disease-associated SNPs. This study provides a comprehensive eQTL map of distal colonic samples obtained from 40 healthy African Americans and demonstrates their relevance for GWAS of colonic diseases. RESULTS: 8.4 million imputed SNPs were tested for their associations with 16,252 expression probes representing 12,363 unique genes. 1,941 significant cis-eQTL, corresponding to 122 independent signals, were identified at a false discovery rate (FDR) of 0.01. Overall, among colon cis-eQTL, there was significant enrichment for GWAS variants for IBD (Crohn's disease [CD] and ulcerative colitis [UC]) and CRC as well as type 2 diabetes and body mass index. ERAP2, ADCY3, INPP5E, UBA7, SFMBT1, NXPE1 and REXO2 were identified as target genes for IBD-associated variants. The CRC-associated eQTL rs3802842 was associated with the expression of C11orf93 (COLCA2). Enrichment of colon eQTL near transcription start sites and for active histone marks was demonstrated, and eQTL with high population differentiation were identified. CONCLUSIONS: Through the comprehensive study of eQTL in the human colon, this study identified novel target genes for IBD- and CRC-associated genetic variants. Moreover, bioinformatic characterization of colon eQTL provides a tissue-specific tool to improve understanding of biological differences in diseases between different ethnic groups.

Joint analysis is more efficient than replication-based analysis for two-stage genome-wide association studies.

Genome-wide association is a promising approach to identify common genetic variants that predispose to human disease. Because of the high cost of genotyping hundreds of thousands of markers on thousands of subjects, genome-wide association studies often follow a staged design in which a proportion (pi(samples)) of the available samples are genotyped on a large number of markers in stage 1, and a proportion (pi(samples)) of these markers are later followed up by genotyping them on the remaining samples in stage 2. The standard strategy for analyzing such two-stage data is to view stage 2 as a replication study and focus on findings that reach statistical significance when stage 2 data are considered alone. We demonstrate that the alternative strategy of jointly analyzing the data from both stages almost always results in increased power to detect genetic association, despite the need to use more stringent significance levels, even when effect sizes differ between the two stages. We recommend joint analysis for all two-stage genome-wide association studies, especially when a relatively large proportion of the samples are genotyped in stage 1 (pi(samples) >or= 0.30), and a relatively large proportion of markers are selected for follow-up in stage 2 (pi(markers) >or= 0.01).

Discovery and functional assessment of gene variants in the vascular endothelial growth factor pathway.

Angiogenesis is a host-mediated mechanism in disease pathophysiology. The vascular endothelial growth factor (VEGF) pathway is a major determinant of angiogenesis, and a comprehensive annotation of the functional variation in this pathway is essential to understand the genetic basis of angiogenesis-related diseases. We assessed the allelic heterogeneity of gene expression, population specificity of cis expression quantitative trait loci (eQTLs), and eQTL function in luciferase assays in CEU and Yoruba people of Ibadan, Nigeria (YRI) HapMap lymphoblastoid cell lines in 23 resequenced genes. Among 356 cis-eQTLs, 155 and 174 were unique to CEU and YRI, respectively, and 27 were shared between CEU and YRI. Two cis-eQTLs provided mechanistic evidence for two genome-wide association study findings. Five eQTLs were tested for function in luciferase assays and the effect of two KRAS variants was concordant with the eQTL effect. Two eQTLs found in each of PRKCE, PIK3C2A, and MAP2K6 could predict 44%, 37%, and 45% of the variance in gene expression, respectively. This is the first analysis focusing on the pattern of functional genetic variation of the VEGF pathway genes in CEU and YRI populations and providing mechanistic evidence for genetic association studies of diseases for which angiogenesis plays a pathophysiologic role.

Shared and independent colorectal cancer risk alleles in TGFß-related genes in African and European Americans.

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.

A genome-wide meta-analysis of nodular sclerosing Hodgkin lymphoma identifies risk loci at 6p21.32.

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 (r(2) = 0.257), and the remaining pairs of SNPs were not correlated (r(2) < 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.

A Case of Non-Syndromic Congenital Cataracts Caused by a Novel MAF Variant in the C-Terminal DNA-Binding Domain-Case Report and Literature Review.

The MAF gene encodes a transcription factor in which pathogenic variants have been associated with both isolated and syndromic congenital cataracts. We aim to review the MAF variants in the C-terminal DNA-binding domain associated with non-syndromic congenital cataracts and describe a patient with a novel, disease-causing de novo missense variant. Published reports of C-terminal MAF variants and their associated congenital cataracts and ophthalmic findings were reviewed. The patient we present and his biological parents had genetic testing via a targeted gene panel followed by trio-based whole exome sequencing. A 4-year-old patient with a history of bilateral nuclear and cortical cataracts was found to have a novel, likely pathogenic de novo variant in MAF, NM_005360.5:c.922A>G (p.Lys308Glu). No syndromic findings or anterior segment abnormalities were identified. We report the novel missense variant, c.922A>G (p.Lys308Glu), in the C-terminal DNA-binding domain of MAF classified as likely pathogenic and associated with non-syndromic bilateral congenital cataracts.

The limits of genome-wide methods for pharmacogenomic testing.

OBJECTIVE: The goal of pharmacogenomics is the translation of genomic discoveries into individualized patient care. Recent advances in the means to survey human genetic variation are fundamentally transforming our understanding of the genetic basis of interindividual variation in therapeutic response. The goal of this study was to systematically evaluate high-throughput genotyping technologies for their ability to assay variation in pharmacogenetically important genes (pharmacogenes). These platforms are either being proposed for or are already being widely used for clinical implementation; therefore, knowledge of coverage of pharmacogenes on these platforms would serve to better evaluate current or proposed pharmacogenetic association studies. METHOD: Among the genes included in our study are drug-metabolizing enzymes, transporters, receptors, and drug targets, of interest to the entire pharmacogenetic community. We considered absolute and linkage disequilibrium (LD)-informed coverage, minor allele frequency spectrum, and functional annotation for a Caucasian population. We also examined the effect of LD, effect size, and cohort size on the power to detect single nucleotide polymorphism associations. RESULTS: In our analysis of 253 pharmacogenes, we found that no platform showed more than 85% coverage of these genes (after accounting for LD). Furthermore, the lack of coverage showed a marked increase at minor allele frequencies of less than 20%. Even after accounting for LD, only 30% of the missense polymorphisms (which are enriched for low-frequency alleles) were covered by HapMap, with still lower coverage on the other platforms. CONCLUSION: We have conducted the first systematic evaluation of the Axiom Genomic Database, Omni 2.5 M, and the Drug Metabolizing Enzymes and Transporters chip. This study is the first to utilize the 1000 Genomes Project to present a comprehensive evaluative framework. Our results provide a much-needed assessment of microarray-based genotyping and next-generation sequencing technologies' ability to survey fully the variation in genes of particular interest to the pharmacogenetics community. Our findings demonstrate the limitations of genome-wide methods and the challenges of implementing pharmacogenomic tests into the clinical context.

Gene-expression-guided selection of candidate loci and molecular phenotype analyses enhance genetic discovery in systemic lupus erythematosus.

Systemic lupus erythematosus (SLE) is a highly heterogeneous autoimmune disorder characterized by differences in autoantibody profiles, serum cytokines, and clinical manifestations. We have previously conducted a case-case genome-wide association study (GWAS) of SLE patients to detect associations with autoantibody profile and serum interferon alpha (IFN-α). In this study, we used public gene expression data sets to rationally select additional single nucleotide polymorphisms (SNPs) for validation. The top 200 GWAS SNPs were searched in a database which compares genome-wide expression data to genome-wide SNP genotype data in HapMap cell lines. SNPs were chosen for validation if they were associated with differential expression of 15 or more genes at a significance of P < 9 × 10(-5). This resulted in 11 SNPs which were genotyped in 453 SLE patients and 418 matched controls. Three SNPs were associated with SLE-associated autoantibodies, and one of these SNPs was also associated with serum IFN-α (P < 4.5 × 10(-3) for all). One additional SNP was associated exclusively with serum IFN-α. Case-control analysis was insensitive to these molecular subphenotype associations. This study illustrates the use of gene expression data to rationally select candidate loci in autoimmune disease, and the utility of stratification by molecular phenotypes in the discovery of additional genetic associations in SLE.

Genetic changes associated with relapse in favorable histology Wilms tumor: A Children's Oncology Group AREN03B2 study.

Over the last decade, sequencing of primary tumors has clarified the genetic underpinnings of Wilms tumor but has not affected therapy, outcome, or toxicity. We now sharpen our focus on relapse samples from the umbrella AREN03B2 study. We show that over 40% of relapse samples contain mutations in SIX1 or genes of the MYCN network, drivers of progenitor proliferation. Not previously seen in large studies of primary Wilms tumors, DIS3 and TERT are now identified as recurrently mutated. The analysis of primary-relapse tumor pairs suggests that 11p15 loss of heterozygosity (and other copy number changes) and mutations in WT1 and MLLT1 typically occur early, but mutations in SIX1, MYCN, and WTX are late developments in some individuals. Most strikingly, 75% of relapse samples had gain of 1q, providing strong conceptual support for studying circulating tumor DNA in clinical trials to better detect 1q gain earlier and monitor response.

Genetic heterogeneity in colorectal cancer associations between African and European americans.

BACKGROUND & AIMS: Genome-wide association studies of colorectal cancer (CRC) have identified risk variants in 10 genomic regions. None of these studies included African Americans, who have the highest incidence and mortality from CRC in the United States. For the 10 genomic regions, we performed an association study of Americans of African and European descent. METHODS: We genotyped 22 single nucleotide polymorphisms (SNPs) in DNA samples from 1194 patients with CRC (795 African Americans and 399 European Americans) and 1352 controls (985 African Americans and 367 European Americans). At chromosome 8q24.21 region 3, we analyzed 6 SNPs from 1000 African American cases and 1393 controls. Association testing was done using multivariate logistic regression controlling for ancestry, age, and sex. RESULTS: Among African Americans, the SNP rs6983267 at 8q24.21 was not associated with CRC (odds ratio, 1.18; P = .12); instead, the 8q24.21 SNP rs7014346 (odds ratio, 1.15; P = .03) was associated with CRC in this population. At 15q13.3, rs10318 was associated with CRC in both populations. At 11q23.1, rs3802842 was significantly associated with rectal cancer risk only among African Americans (odds ratio, 1.34; P = .01); this observation was made in previous studies. Among European Americans, SNPs at 8q24.21, 11q23.1, and 16q22.1 were significantly associated with CRC, and the odds ratios were of the same magnitude and direction for all SNPs tested, consistent with previously published studies. In contrast, in African Americans, the opposite allele of rs10795668 at 10p14 was associated with colorectal cancer (odds ratio, 1.35; P = .04), and altogether the odds ratios were in the opposite direction for 9 of the 22 SNPs tested. CONCLUSIONS: There is genetic heterogeneity in CRC associations in Americans of African versus European descent.

Fine-mapping alleles for body weight in LG/J × SM/J F2 and F(34) advanced intercross lines.

The present study measured variation in body weight using a combined analysis in an F(2) intercross and an F(34) advanced intercross line (AIL). Both crosses were derived from inbred LG/J and SM/J mice, which were selected for large and small body size prior to inbreeding. Body weight was measured at 62 (± 5) days of age. Using an integrated GWAS and forward model selection approach, we identified 11 significant QTLs that affected body weight on ten different chromosomes. With these results we developed a full model that explained over 18% of the phenotypic variance. The median 1.5-LOD support interval was 5.55 Mb, which is a significant improvement over most prior body weight QTLs. We identified nonsynonymous coding SNPs between LG/J and SM/J mice in order to further narrow the list of candidate genes. Three of the genes with nonsynonymous coding SNPs (Rad23b, Stk33, and Anks1b) have been associated with adiposity, waist circumference, and body mass index in human GWAS, thus providing evidence that these genes may underlie our QTLs. Our results demonstrate that a relatively small number of loci contribute significantly to the phenotypic variance in body weight, which is in marked contrast to the situation in humans. This difference is likely to be the result of strong selective pressure and the simplified genetic architecture, both of which are important advantages of our system.

Genome-wide association study to identify novel loci associated with therapy-related myeloid leukemia susceptibility.

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 (P(enrich) = .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.

QTLRel: an R package for genome-wide association studies in which relatedness is a concern.

BACKGROUND: Existing software for quantitative trait mapping is either not able to model polygenic variation or does not allow incorporation of more than one genetic variance component. Improperly modeling the genetic relatedness among subjects can result in excessive false positives. We have developed an R package, QTLRel, to enable more flexible modeling of genetic relatedness as well as covariates and non-genetic variance components. RESULTS: We have successfully used the package to analyze many datasets, including F34 body weight data that contains 688 individuals genotyped at 3105 SNP markers and identified 11 QTL. It took 295 seconds to estimate variance components and 70 seconds to perform the genome scan on an Linux machine equipped with a 2.40GHz Intel(R) Core(TM)2 Quad CPU. CONCLUSIONS: QTLRel provides a toolkit for genome-wide association studies that is capable of calculating genetic incidence matrices from pedigrees, estimating variance components, performing genome scans, incorporating interactive covariates and genetic and non-genetic variance components, as well as other functionalities such as multiple-QTL mapping and genome-wide epistasis.

Functional Common and Rare ERBB2 Germline Variants Cooperate in Familial and Sporadic Cancer Susceptibility.

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.