Search

PURPOSE: Mendelian disease genomic research has undergone a massive transformation over the past decade. With increasing availability of exome and genome sequencing, the role of Mendelian research has expanded beyond data collection, sequencing, and analysis to worldwide data sharing and collaboration. METHODS: Over the past 10 years, the National Institutes of Health–supported Centers for Mendelian Genomics (CMGs) have played a major role in this research and clinical evolution. RESULTS: We highlight the cumulative gene discoveries facilitated by the program, biomedical research leveraged by the approach, and the larger impact on the research community. Beyond generating a list of gene-phenotype relationships and participating in widespread data sharing, the CMGs have created resources, tools, and training for the larger community to foster understanding of genes and genome variation. The CMGs have participated in a wide range of data sharing activities, including deposition of all eligible CMG data into the Analysis, Visualization, and Informatics Lab-space (AnVIL), sharing candidate genes through the Matchmaker Exchange and the CMG website, and sharing variants in Genotypes to Mendelian Phenotypes (Geno2MP) and VariantMatcher. CONCLUSION: The work is far from complete; strengthening communication between research and clinical realms, continued development and sharing of knowledge and tools, and improving access to richly characterized data sets are all required to diagnose the remaining molecularly undiagnosed patients.

Background and AimsThe microbiome has long been suspected of a role in colorectal cancer (CRC) tumorigenesis. The mutational signature SBS88 mechanistically links CRC development with the strain ofEscherichia coliharboring thepksisland that produces the genotoxin colibactin, but the genomic, pathological and survival characteristics associated with SBS88-positive tumors are unknown.MethodsSBS88-positive CRCs were identified from targeted sequencing data from 5,292 CRCs from 17 studies and tested for their association with clinico-pathological features, oncogenic pathways, genomic characteristics and survival.ResultsIn total, 7.5% (398/5,292) of the CRCs were SBS88-positive, of which 98.7% (392/398) were microsatellite stable/microsatellite instability low (MSS/MSI-L), compared with 80% (3916/4894) of SBS88 negative tumors (p=1.5×10-28). Analysis of MSS/MSI-L CRCs demonstrated that SBS88 positive CRCs were associated with the distal colon (OR=1.84, 95% CI=1.40-2.42, p=1×10-5) and rectum (OR=1.90, 95% CI=1.44-2.51, p=6×10-6) tumor sites compared with the proximal colon. The top seven recurrent somatic mutations associated with SBS88-positive CRCs demonstrated mutational contexts associated with colibactin-induced DNA damage, the strongest of which was theAPC:c.835-8A>G mutation (OR=65.5, 95%CI=39.0-110.0, p=3×10-80). Large copy number alterations (CNAs) including CNA loss on 14q and gains on 13q, 16q and 20p were significantly enriched in SBS88- positive CRCs. SBS88-positive CRCs were associated with better CRC-specific survival (p=0.007; hazard ratio of 0.69, 95% CI=0.52-0.90) when stratified by age, sex, study, and by stage.ConclusionSBS88-positivity, a biomarker of colibactin-induced DNA damage, can identify a novel subtype of CRC characterized by recurrent somatic mutations, copy number alterations and better survival. These findings provide new insights for treatment and prevention strategies for this subtype of CRC.

Little is known regarding the potential relationship between clonal hematopoiesis (CH) of indeterminate potential (CHIP), which is the expansion of hematopoietic stem cells with somatic mutations, and risk of prostate cancer, the fifth leading cause of cancer death of men worldwide. We evaluated the association of age-related CHIP with overall and aggressive prostate cancer risk in two large whole-exome sequencing studies of 75 047 European ancestry men, including 7663 prostate cancer cases, 2770 of which had aggressive disease, and 3266 men carrying CHIP variants. We found that CHIP, defined by over 50 CHIP genes individually and in aggregate, was not significantly associated with overall (aggregate HR = 0.93, 95% CI = 0.76–1.13, P = 0.46) or aggressive (aggregate OR = 1.14, 95% CI = 0.92–1.41, P = 0.22) prostate cancer risk. CHIP was weakly associated with genetic risk of overall prostate cancer, measured using a polygenic risk score (OR = 1.05 per unit increase, 95% CI = 1.01–1.10, P = 0.01). CHIP was not significantly associated with carrying pathogenic/likely pathogenic/deleterious variants in DNA repair genes, which have previously been found to be associated with aggressive prostate cancer. While findings from this study suggest that CHIP is likely not a risk factor for prostate cancer, it will be important to investigate other types of CH in association with prostate cancer risk.

Here we describe MyGene2, Geno2MP, VariantMatcher, and Franklin; databases that provide variant-level information and phenotypic features to researchers, clinicians, healthcare providers and patients. Following the footsteps of the Matchmaker Exchange project that connects exome, genome, and phenotype databases at the gene level, these databases have as one goal to facilitate connection to one another using Data Connect, a standard for discovery and search of biomedical data from the Global Alliance for Genomics and Health (GA4GH).

GeneMatcher (genematcher.org) is a tool designed to connect individuals with an interest in the same gene. Now used around the world to create collaborations and generate the evidence needed to support novel disease gene identification, GeneMatcher is a founding member of the Matchmaker Exchange (MME; matchmakerexchange.org) and strongest possible advocate for global data sharing including those in resource-limited environments. As of October 1, 2021, there are 12,531 submitters from 94 countries who have submitted 58,134 submissions with 13,498 unique genes in the database. Among these genes, 8970 (64%) have matched at least once and the total number of matches is 378,806, growing by about 10,000 per month. GeneMatcher submitters increase by 80-120 each month and submissions grow by800 per month, while unique genes and gene matches continue to grow steadily at rate of about 80 per month. The number of genes without a match peaked at 4371 in February of 2019 and despite the increase in the number of new submissions, the number of unique genes without a match continues to slowly decline, currently standing at 4,016. All submissions in GeneMatcher are available for matching across the MME.

TET3 encodes an essential dioxygenase involved in epigenetic regulation through DNA demethylation. TET3 deficiency, or Beck-Fahrner syndrome (BEFAHRS; MIM: 618798), is a recently described neurodevelopmental disorder of the DNA demethylation machinery with a nonspecific phenotype resembling other chromatin-modifying disorders, but inconsistent variant types and inheritance patterns pose diagnostic challenges. Given TET3’s direct role in regulating 5-methylcytosine and recent identification of syndrome-specific DNA methylation profiles, we analyzed genome-wide DNA methylation in whole blood of TET3-deficient individuals and identified an episignature that distinguishes affected and unaffected individuals and those with mono-allelic and bi-allelic pathogenic variants. Validation and testing of the episignature correctly categorized known TET3 variants and determined pathogenicity of variants of uncertain significance. Clinical utility was demonstrated when the episignature alone identified an affected individual from over 1000 undiagnosed cases and was confirmed upon distinguishing TET3-deficient individuals from those with 46 other disorders. The TET3-deficient signature - and the signature resulting from activating mutations in DNMT1 which normally opposes TET3 - are characterized by hypermethylation, which for BEFAHRS involves CpG sites that may be biologically relevant. This work expands the role of epi-phenotyping in molecular diagnosis and reveals genome-wide DNA methylation profiling as a quantitative, functional readout for characterization of this new biochemical category of disease.

Carriers of germline biallelic pathogenic variants in the MUTYH gene have a high risk of colorectal cancer. We test 5649 colorectal cancers to evaluate the discriminatory potential of a tumor mutational signature specific to MUTYH for identifying biallelic carriers and classifying variants of uncertain clinical significance (VUS). Using a tumor and matched germline targeted multi-gene panel approach, our classifier identifies all biallelic MUTYH carriers and all known non-carriers in an independent test set of 3019 colorectal cancers (accuracy = 100% (95% confidence interval 99.87–100%)). All monoallelic MUTYH carriers are classified with the non-MUTYH carriers. The classifier provides evidence for a pathogenic classification for two VUS and a benign classification for five VUS. Somatic hotspot mutations KRAS p.G12C and PIK3CA p.Q546K are associated with colorectal cancers from biallelic MUTYH carriers compared with non-carriers (p = 2 × 10−23 and p = 6 × 10−11, respectively). Here, we demonstrate the potential application of mutational signatures to tumor sequencing workflows to improve the identification of biallelic MUTYH carriers.

Background With the advent of whole exome (ES) and genome sequencing (GS) as tools for disease gene discovery, rare variant filtering, prioritization and data sharing have become essential components of the search for disease genes and variants potentially contributing to disease phenotypes. The computational storage, data manipulation, and bioinformatic interpretation of thousands to millions of variants identified in ES and GS, respectively, is a challenging task. To aid in that endeavor, we constructed PhenoDB, GeneMatcher and VariantMatcher. Results PhenoDB is an accessible, freely available, web-based platform that allows users to store, share, analyze and interpret their patients’ phenotypes and variants from ES/GS data. GeneMatcher is accessible to all stakeholders as a web-based tool developed to connect individuals (researchers, clinicians, health care providers and patients) around the globe with interest in the same gene(s), variant(s) or phenotype(s). Finally, VariantMatcher was developed to enable public sharing of variant-level data and phenotypic information from individuals sequenced as part of multiple disease gene discovery projects. Here we provide updates on PhenoDB and GeneMatcher applications and implementation and introduce VariantMatcher. Conclusion Each of these tools has facilitated worldwide data sharing and data analysis and improved our ability to connect genes to phenotypic traits. Further development of these platforms will expand variant analysis, interpretation, novel disease-gene discovery and facilitate functional annotation of the human genome for clinical genomics implementation and the precision medicine initiative.

The All of Us Research Program (AoURP, ‘the program’) is an initiative, sponsored by the National Institutes of Health (NIH), that aims to enroll one million people (or more) across the United States. Through repeated engagement of participants, a research resource is being created to enable a variety of future observational and interventional studies. The program has also committed to genomic data generation and returning important health-related information to participants. To do so, whole genome sequencing (WGS), variant calling processes, data interpretation, and return-of-results procedures had to be created and receive an Investigational Device Exemption (IDE) from the United States Food and Drug Administration (FDA). The performance of the entire workflow was assessed through the largest known cross-center, WGS-based, validation activity that was refined iteratively through interactions with the FDA over many months. The accuracy and precision of the WGS process as a device for the return of certain health-related genomic results was determined to be sufficient, and an IDE was granted. We present here both the process of navigating the IDE application process with the FDA and the results of the validation study as a guide to future projects which may need to follow a similar path. Future supplements to the IDE will be submitted to support additional variant classes, sample types, and any expansion to the reportable regions.

Rare pathogenic variants in DNA repair genes have been found to influence risk of aggressive prostate cancer. We conducted a large case-only exome sequencing study to further understand the role of rare coding variation in aggressive prostate cancer in a study of 9,185 aggressive (prostate cancer death, metastatic disease, T4, or both T3 and Gleason≥8) and 8,361 non-aggressive cases (T1/T2 and Gleason≤6) of European ancestry from 19 international studies. Stage 1 samples (N=5,545) had whole exome-sequencing, and stage 2 samples (N=12,001) had targeted exome sequencing for 1,459 genes selected based on stage 1 results and previous evidence. Logistic regression models were used to evaluate gene-based tests and the aggregate effect of multiple genes to investigate whether carrying pathogenic/likely pathogenic/deleterious (P/LP/D) variants (18,759 identified) was associated with risk of aggressive prostate cancer, prostate cancer death (N=6,033), or metastatic disease (N=1,730) compared to non-aggressive disease. Gene-based tests were meta-analyzed across stages 1 and 2. BRCA2, ATM, and NBN had the most statistically significant gene-based results: BRCA2 P/LP/D variant carriers had 4.3-fold higher odds of aggressive disease (95% CI=3.15-5.86, P=4x10-20), 4.7-fold higher odds of prostate cancer death (95% CI=3.41-6.59, P=2x10-20), and 5.7-fold higher odds of metastatic disease (95% CI=3.71-8.76, P=2x10-15); ATM P/LP/D variant carriers had 2.2-fold higher odds of aggressive disease (95% CI=1.58-2.99, P=2x10-6), 2.2-fold higher odds of prostate cancer death (95% CI=1.52-3.05, P=2x10-5), and 3.0-fold higher odds of metastatic disease (95% CI=1.93-4.61, P=9x10-7); and NBN P/LP/D variant carriers had 5.9-fold higher odds of metastatic disease (95% CI=2.56-13.84, P=3x10-5). Among potentially novel genes with strong but not exome-wide significant statistical evidence were MMP19, involved in reproduction and metastasis, with carriers having 2.8-fold higher odds of prostate cancer death (95% CI=1.53-5.05, P=8x10-4); PKD2L2, involved in fertility, with carriers having 3.5-fold higher odds of prostate cancer death (95% CI=1.76-7.04, P=5x10-4); and SMPD1, involved in converting sphingomyelin to ceramide, with carriers having 5.3-fold higher odds of metastatic disease (95% CI=1.85-14.98, P=0.002). At least one P/LP/D variant within a subset of 24 previously curated candidate prostate cancer DNA repair genes was carried by 12.8% of aggressive cases (OR=1.48, 95% CI=1.34-1.64, P=3x10-14), 12.6% of cases who died due to prostate cancer (OR=1.47, 95% CI=1.31-1.65, P=3x10-11), and 15.1% of metastatic cases (OR=2.16, 95% CI=1.57-2.16, P=5x10-14) compared to 9.4% of non-aggressive cases. These findings support the importance of rare genetic variation in aggressive prostate cancer risk and may have important implications for prostate cancer risk stratification and screening. Citation Format: Burcu F. Darst, Ed Saunders, Tokhir Dadaev, Xin Sheng, Peggy Wan, Loreall Pooler, Lucy Y. Xia, Stephen Chanock, Sonja I. Berndt, Susan M. Gapstur, Victoria Stevens, Demetrius Albanes, Stephanie J. Weinstein, Vincent Gnanapragasam, Graham G. Giles, Tu Nguyen-Dumont, Roger L. Milne, Mark M. Pomerantz, Julie A. Schmidt, Ruth C. Travis, Timothy J. Key, Konrad H. Stopsack, Lorelei A. Mucci, William J. Catalona, Beth Marosy, Kurt N. Hetrick, Kimberly F. Doheny, Robert J. MacInnis, Melissa C. Southey, Rosalind A. Eeles, Fredrik Wiklund, Zsofia Kote-Jarai, David V. Conti, Christopher A. Haiman. Multi-stage exome sequencing study of 17,546 aggressive and non-aggressive prostate cancer cases [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2022; 2022 Apr 8-13. Philadelphia (PA): AACR; Cancer Res 2022;82(12_Suppl):Abstract nr 688.

Background There is an urgent need to identify factors specifically associated with aggressive prostate cancer (PCa) risk. We investigated whether rare pathogenic, likely pathogenic, or deleterious (P/LP/D) germline variants in DNA repair genes are associated with aggressive PCa risk in a case-case study of aggressive vs nonaggressive disease. Methods Participants were 5545 European-ancestry men, including 2775 nonaggressive and 2770 aggressive PCa cases, which included 467 metastatic cases (16.9%). Samples were assembled from 12 international studies and germline sequenced together. Rare (minor allele frequency Results BRCA2 and PALB2 had the most statistically significant gene-based associations, with 2.5% of aggressive and 0.8% of nonaggressive cases carrying P/LP/D BRCA2 alleles (odds ratio [OR] = 3.19, 95% confidence interval [CI] = 1.94 to 5.25, P = 8.58 × 10-7) and 0.65% of aggressive and 0.11% of nonaggressive cases carrying P/LP/D PALB2 alleles (OR = 6.31, 95% CI = 1.83 to 21.68, P = 4.79 × 10-4). ATM had a nominal association, with 1.6% of aggressive and 0.8% of nonaggressive cases carrying P/LP/D ATM alleles (OR = 1.88, 95% CI = 1.10 to 3.22, P = .02). In aggregate, P/LP/D alleles within 24 literature-curated candidate PCa DNA repair genes were more common in aggressive than nonaggressive cases (carrier frequencies = 14.2% vs 10.6%, respectively; P = 5.56 × 10-5). However, this difference was non-statistically significant (P = .18) on excluding BRCA2, PALB2, and ATM. Among these 24 genes, P/LP/D carriers had a 1.06-year younger diagnosis age (95% CI = -1.65 to 0.48, P = 3.71 × 10-4). Conclusions Risk conveyed by DNA repair genes is largely driven by rare P/LP/D alleles within BRCA2, PALB2, and ATM. These findings support the importance of these genes in both screening and disease management considerations.

To further dissect the genetic architecture of colorectal cancer (CRC), we performed whole-genome sequencing of 1,439 cases and 720 controls, imputed discovered sequence variants and Haplotype Reference Consortium panel variants into genome-wide association study data, and tested for association in 34,869 cases and 29,051 controls. Findings were followed up in an additional 23,262 cases and 38,296 controls. We discovered a strongly protective 0.3% frequency variant signal at CHD1. In a combined meta-analysis of 125,478 individuals, we identified 40 new independent signals at P, Reporting Summary. Further information on experimental design is available in the Life Sciences Reporting Summary linked to this article.

Importance Advanced age-related macular degeneration (AMD) is a leading cause of blindness in Western countries. Causal, modifiable risk factors need to be identified to develop preventive measures for advanced AMD. Objective To assess whether smoking, alcohol consumption, blood pressure, body mass index, and glycemic traits are associated with increased risk of advanced AMD. Design, Setting, Participants This study used 2-sample mendelian randomization. Genetic instruments composed of variants associated with risk factors at genome-wide significance (P

Kabuki syndrome is a monogenic disorder caused by loss of function variants in either of two genes encoding histone-modifying enzymes. We performed targeted sequencing in a cohort of 27 probands with a clinical diagnosis of Kabuki syndrome. Of these, 12 had causative variants in the two known Kabuki syndrome genes. In 2, we identified presumptive loss of function de novo variants in KMT2A (missense and splice site variants), a gene that encodes another histone modifying enzyme previously exclusively associated with Wiedermann-Steiner syndrome. Although Kabuki syndrome is a disorder of histone modification, we also find alterations in DNA methylation among individuals with a Kabuki syndrome diagnosis relative to matched normal controls, regardless of whether they carry a variant in KMT2A or KMT2D or not. Furthermore, we observed characteristic global abnormalities of DNA methylation that distinguished patients with a loss of function variant in KMT2D or missense or splice site variants in either KMT2D or KMT2A from normal controls. Our results provide new insights into the relationship of genotype to epigenotype and phenotype and indicate cross-talk between histone and DNA methylation machineries exposed by inborn errors of the epigenetic apparatus.

There has been extensive debate about both the necessity of orthogonal confirmation of next-generation sequencing (NGS) results in Clinical Laboratory Improvement Amendments-approved laboratories and return of research NGS results to participants enrolled in research studies. In eMERGE-PGx, subjects underwent research NGS using PGRNseq and orthogonal targeted genotyping in clinical laboratories, which prompted a comparison of genotyping results between platforms. Concordance (percentage agreement) was reported for 4077 samples tested across nine combinations of research and clinical laboratories. Retesting was possible on a subset of 1792 samples, and local laboratory directors determined sources of genotype discrepancy. Research NGS and orthogonal clinical genotyping had an overall per sample concordance rate of 0.972 and per variant concordance rate of 0.997. Genotype discrepancies attributed to research NGS were because of sample switching (preanalytical errors), whereas the majority of genotype discrepancies (92.3%) attributed to clinical genotyping were because of allele dropout as a result of rare variants interfering with primer hybridization (analytical errors). These results highlight the analytical quality of clinically significant pharmacogenetic variants derived from NGS and reveal important areas for research and clinical laboratories to address with quality management programs.

The breast cancer risk variants identified in genome-wide association studies explain only a small fraction of the familial relative risk, and the genes responsible for these associations remain largely unknown. To identify novel risk loci and likely causal genes, we performed a transcriptome-wide association study evaluating associations of genetically predicted gene expression with breast cancer risk in 122,977 cases and 105,974 controls of European ancestry. We used data from the Genotype-Tissue Expression Project to establish genetic models to predict gene expression in breast tissue and evaluated model performance using data from The Cancer Genome Atlas. Of the 8,597 genes evaluated, significant associations were identified for 48 at a Bonferroni-corrected threshold of P −6, including 14 genes at loci not yet reported for breast cancer. We silenced 13 genes and showed an effect for 11 on cell proliferation and/or colony-forming efficiency. Our study provides new insights into breast cancer genetics and biology.

Identifying genes and variants contributing to rare disease phenotypes and Mendelian conditions informs biology and medicine, yet potential phenotypic consequences for variation of >75% of the ~20,000 annotated genes in the human genome are lacking. Technical advances to assess rare variation genome-wide, particularly exome sequencing (ES), enabled establishment in the United States of the National Institutes of Health (NIH)-supported Centers for Mendelian Genomics (CMGs) and have facilitated collaborative studies resulting in novel “disease gene” discoveries. Pedigree-based genomic studies and rare variant analyses in families with suspected Mendelian conditions have led to the elucidation of hundreds of novel disease genes and highlighted the impact of de novo mutational events, somatic variation underlying nononcologic traits, incompletely penetrant alleles, phenotypes with high locus heterogeneity, and multilocus pathogenic variation. Herein, we highlight CMG collaborative discoveries that have contributed to understanding both rare and common diseases and discuss opportunities for future discovery in single-locus Mendelian disorder genomics. Phenotypic annotation of all human genes; development of bioinformatic tools and analytic methods; exploration of non-Mendelian modes of inheritance including reduced penetrance, multilocus variation, and oligogenic inheritance; construction of allelic series at a locus; enhanced data sharing worldwide; and integration with clinical genomics are explored. Realizing the full contribution of rare disease research to functional annotation of the human genome, and further illuminating human biology and health, will lay the foundation for the Precision Medicine Initiative.

Recent technological advancements have permitted high-throughput measurement of the human genome, epigenome, metabolome, transcriptome, and proteome at the population level. We hypothesized that subsets of genes identified from omic studies might have closely related biological functions and thus might interact directly at the network level. Therefore, we conducted an integrative analysis of multi-omic datasets of non-small cell lung cancer (NSCLC) to search for association patterns beyond the genome and transcriptome. A large, complex, and robust gene network containing well-known lung cancer-related genes, including EGFR and TERT, was identified from combined gene lists for lung adenocarcinoma. Members of the hypoxia-inducible factor (HIF) gene family were at the center of this network. Subsequent sequencing of network hub genes within a subset of samples from the Transdisciplinary Research in Cancer of the Lung-International Lung Cancer Consortium (TRICL-ILCCO) consortium revealed a SNP (rs12614710) in EPAS1 associated with NSCLC that reached genome-wide significance (OR = 1.50; 95% CI: 1.31–1.72; p = 7.75 × 10−9). Using imputed data, we found that this SNP remained significant in the entire TRICL-ILCCO consortium (p = .03). Additional functional studies are warranted to better understand interrelationships among genetic polymorphisms, DNA methylation status, and EPAS1 expression., • Previous Evidence • As a major participant of Transdisciplinary Research in Cancer of the Lung-International Lung Cancer Consortium (TRICL-ILCCO), we continued monitoring the progress of genomic and other omics studies of lung cancer through regular communications within the consortium and literature search using the PubMed database. The TRICL-ILCCO is one of the largest international consortium of lung cancer which includes many major ongoing lung cancer case-control and cohort studies with the aim of sharing comparable data. Although GWAS successfully defined many lung cancer-associated genomic loci, the integrative studies of multi-dimensional high throughput “-omics” measurements from tumor tissues and corresponding blood specimens are limited. • Added value of this study • We identified a HIFs-EGFR-HDAC4-TERT network associated with lung adenocarcinoma, and subsequent sequencing of network hub genes identified a new locus within EPAS1 that is associated with lung cancer risk. This locus is in hub gene EPAS1, which is a key member of the HIF family involved in every aspect of cancer development and progression. • Implications of all the available evidence • We developed a network building approach for the integrative analysis of multi-omic datasets. The integration of multi-dimensional high throughput “-omics” measurements from tumor tissues and corresponding blood specimens, together with new systems strategies for diagnostics, enables the identification of cancer biomarkers that will facilitate pre-symptomatic diagnosis, stratification of disease, assessment of disease progression, evaluation of patient response to therapy, and identification of recurrences.

The common nonsynonymous variant rs16969968 in the α5 nicotinic receptor subunit gene (CHRNA5) is the strongest genetic risk factor for nicotine dependence in European Americans and contributes to risk in African Americans. To comprehensively examine whether other CHRNA5 coding variation influences nicotine dependence risk, we performed targeted sequencing on 1582 nicotine-dependent cases (Fagerstrom Test for Nicotine Dependence score⩾4) and 1238 non-dependent controls, with independent replication of common and low frequency variants using 12 studies with exome chip data. Nicotine dependence was examined using logistic regression with individual common variants (minor allele frequency (MAF)⩾0.05), aggregate low frequency variants (0.05>MAF⩾0.005) and aggregate rare variants (MAF

Discovering the genetic basis of a Mendelian phenotype establishes a causal link between genotype and phenotype, making possible carrier and population screening and direct diagnosis. Such discoveries also contribute to our knowledge of gene function, gene regulation, development, and biological mechanisms that can be used for developing new therapeutics. As of February 2015, 2,937 genes underlying 4,163 Mendelian phenotypes have been discovered, but the genes underlying ∼50% (i.e., 3,152) of all known Mendelian phenotypes are still unknown, and many more Mendelian conditions have yet to be recognized. This is a formidable gap in biomedical knowledge. Accordingly, in December 2011, the NIH established the Centers for Mendelian Genomics (CMGs) to provide the collaborative framework and infrastructure necessary for undertaking large-scale whole-exome sequencing and discovery of the genetic variants responsible for Mendelian phenotypes. In partnership with 529 investigators from 261 institutions in 36 countries, the CMGs assessed 18,863 samples from 8,838 families representing 579 known and 470 novel Mendelian phenotypes as of January 2015. This collaborative effort has identified 956 genes, including 375 not previously associated with human health, that underlie a Mendelian phenotype. These results provide insight into study design and analytical strategies, identify novel mechanisms of disease, and reveal the extensive clinical variability of Mendelian phenotypes. Discovering the gene underlying every Mendelian phenotype will require tackling challenges such as worldwide ascertainment and phenotypic characterization of families affected by Mendelian conditions, improvement in sequencing and analytical techniques, and pervasive sharing of phenotypic and genomic data among researchers, clinicians, and families.

Bioinformatics approaches to examine gene-gene models provide a means to discover interactions between multiple genes that underlie complex disease. Extensive computational demands and adjusting for multiple testing make uncovering genetic interactions a challenge. Here, we address these issues using our knowledge-driven filtering method, Biofilter, to identify putative single nucleotide polymorphism (SNP) interaction models for cataract susceptibility, thereby reducing the number of models for analysis. Models were evaluated in 3,377 European Americans (1,185 controls, 2,192 cases) from the Marshfield Clinic, a study site of the Electronic Medical Records and Genomics (eMERGE) Network, using logistic regression. All statistically significant models from the Marshfield Clinic were then evaluated in an independent dataset of 4,311 individuals (742 controls, 3,569 cases), using independent samples from additional study sites in the eMERGE Network: Mayo Clinic, Group Health/University of Washington, Vanderbilt University Medical Center, and Geisinger Health System. Eighty-three SNP-SNP models replicated in the independent dataset at likelihood ratio test P < 0.05. Among the most significant replicating models was rs12597188 (intron of CDH1)–rs11564445 (intron of CTNNB1). These genes are known to be involved in processes that include: cell-to-cell adhesion signaling, cell-cell junction organization, and cell-cell communication. Further Biofilter analysis of all replicating models revealed a number of common functions among the genes harboring the 83 replicating SNP-SNP models, which included signal transduction and PI3K-Akt signaling pathway. These findings demonstrate the utility of Biofilter as a biology-driven method, applicable for any genome-wide association study dataset.

Summary/AbstractGenome-wide association studies (GWAS) have laid the foundation for investigations into the biology of complex traits, drug development, and clinical guidelines. However, the dominance of European-ancestry populations in GWAS creates a biased view of the role of human variation in disease, and hinders the equitable translation of genetic associations into clinical and public health applications. The Population Architecture using Genomics and Epidemiology (PAGE) study conducted a GWAS of 26 clinical and behavioral phenotypes in 49,839 non-European individuals. Using strategies designed for analysis of multi-ethnic and admixed populations, we confirm 574 GWAS catalog variants across these traits, and find 38 secondary signals in known loci and 27 novel loci. Our data shows strong evidence of effect-size heterogeneity across ancestries for published GWAS associations, substantial benefits for fine-mapping using diverse cohorts, and insights into clinical implications. We strongly advocate for continued, large genome-wide efforts in diverse populations to reduce health disparities.

Genome-wide association studies (GWAS) have laid the foundation for investigations into the biology of complex traits, drug development and clinical guidelines. However, the majority of discovery efforts are based on data from populations of European ancestry1–3. In light of the differential genetic architecture that is known to exist between populations, bias in representation can exacerbate existing disease and healthcare disparities. Critical variants may be missed if they have a low frequency or are completely absent in European populations, especially as the field shifts its attention towards rare variants, which are more likely to be population-specific4–10. Additionally, effect sizes and their derived risk prediction scores derived in one population may not accurately extrapolate to other populations11,12. Here we demonstrate the value of diverse, multi-ethnic participants in large-scale genomic studies. The Population Architecture using Genomics and Epidemiology (PAGE) study conducted a GWAS of 26 clinical and behavioural phenotypes in 49,839 non-European individuals. Using strategies tailored for analysis of multi-ethnic and admixed populations, we describe a framework for analysing diverse populations, identify 27 novel loci and 38 secondary signals at known loci, as well as replicate 1,444 GWAS catalogue associations across these traits. Our data show evidence of effect-size heterogeneity across ancestries for published GWAS associations, substantial benefits for fine-mapping using diverse cohorts and insights into clinical implications. In the United States—where minority populations have a disproportionately higher burden of chronic conditions13—the lack of representation of diverse populations in genetic research will result in inequitable access to precision medicine for those with the highest burden of disease. We strongly advocate for continued, large genome-wide efforts in diverse populations to maximize genetic discovery and reduce health disparities. © 2019, The Author(s), under exclusive licence to Springer Nature Limited.

Genotyping arrays have been widely adopted as an efficient means to interrogate variation across the human genome. Genetic variants may be observed either directly, via genotyping, or indirectly, through linkage disequilibrium with a genotyped variant. The total proportion of genomic variation captured by an array, either directly or indirectly, is referred to as “genomic coverage.” Here we use genotype imputation and Phase 3 of the 1000 Genomes Project to assess genomic coverage of several modern genotyping arrays. We find that in general, coverage increases with increasing array density. However, arrays designed to cover specific populations may yield better coverage in those populations compared to denser arrays not tailored to the given population. Ultimately, array choice involves trade-offs between cost, density, and coverage, and our work helps inform investigators weighing these choices and trade-offs.

Author(s): Taylor, Kimberly E; Wong, Quenna; Levine, David M; McHugh, Caitlin; Laurie, Cathy; Doheny, Kimberly; Lam, Mi Y; Baer, Alan N; Challacombe, Stephen; Lanfranchi, Hector; Schiodt, Morten; Srinivasan, M; Umehara, Hisanori; Vivino, Frederick B; Zhao, Yan; Shiboski, Stephen C; Daniels, Troy E; Greenspan, John S; Shiboski, Caroline H; Criswell, Lindsey A | Abstract: ObjectiveThe Sjogren's International Collaborative Clinical Alliance (SICCA) is an international data registry and biorepository derived from a multisite observational study of participants in whom genotyping was performed on the Omni2.5M platform and who had undergone deep phenotyping using common protocol-directed methods. The aim of this study was to examine the genetic etiology of Sjogren's syndrome (SS) across ancestry and disease subsets.MethodsWe performed genome-wide association study analyses using SICCA subjects and external controls obtained from dbGaP data sets, one using all participants (1,405 cases, 1,622 SICCA controls, and 3,125 external controls), one using European participants (585, 966, and 580, respectively), and one using Asian participants (460, 224, and 901, respectively) with ancestry adjustments via principal components analyses. We also investigated whether subphenotype distributions differ by ethnicity, and whether this contributes to the heterogeneity of genetic associations.ResultsWe observed significant associations in established regions of the major histocompatibility complex (MHC), IRF5, and STAT4 (P = 3 × 10-42 , P = 3 × 10-14 , and P = 9 × 10-10 , respectively), and several novel suggestive regions (those with 2 or more associations at P l 1 × 10-5 ). Two regions have been previously implicated in autoimmune disease: KLRG1 (P = 6 × 10-7 [Asian cluster]) and SH2D2A (P = 2 × 10-6 [all participants]). We observed striking differences between the associations in Europeans and Asians, with high heterogeneity especially in the MHC; representative single-nucleotide polymorphisms from established and suggestive regions had highly significant differences in the allele frequencies in the study populations. We showed that SSA/SSB autoantibody production and the labial salivary gland focus score criteria were associated with the first worldwide principal component, indicative of higher non-European ancestry (P = 4 × 10-15 and P = 4 × 10-5 , respectively), but that subphenotype differences did not explain most of the ancestry differences in genetic associations.ConclusionGenetic associations with SS differ markedly according to ancestry; however, this is not explained by differences in subphenotypes.

The OCAC OncoArray genotyping project: et al., To identify common alleles associated with different histotypes of epithelial ovarian cancer (EOC), we pooled data from multiple genome-wide genotyping projects totaling 25,509 EOC cases and 40,941 controls. We identified nine new susceptibility loci for different EOC histotypes: six for serous EOC histotypes (3q28, 4q32.3, 8q21.11, 10q24.33, 18q11.2 and 22q12.1), two for mucinous EOC (3q22.3 and 9q31.1) and one for endometrioid EOC (5q12.3). We then performed meta-analysis on the results for high-grade serous ovarian cancer with the results from analysis of 31,448 BRCA1 and BRCA2 mutation carriers, including 3,887 mutation carriers with EOC. This identified three additional susceptibility loci at 2q13, 8q24.1 and 12q24.31. Integrated analyses of genes and regulatory biofeatures at each locus predicted candidate susceptibility genes, including OBFC1, a new candidate susceptibility gene for low-grade and borderline serous EOC., The OCAC OncoArray genotyping project was funded through grants from the US National Institutes of Health (CA1X01HG007491-01 (C.I.A.), U19-CA148112 (T.A.S.), R01-CA149429 (C.M.P.) and R01-CA058598 (M.T.G.)); Canadian Institutes of Health Research (MOP-86727 (L.E.K.)); and the Ovarian Cancer Research Fund (A.B.). Funding for the CIMBA OncoArray genotyping was provided by the Government of Canada through Genome Canada and the Canadian Institutes of Health Research, the Ministère de l’Économie, de la Science et de l’Innovation du Québec through Génome Québec, the Quebec Breast Cancer Foundation for the PERSPECTIVE project, the US National Institutes of Health (CA1X01HG007491-01 (C.I.A.)), the Odense University Hospital Research Foundation (M.T.), the National R&D Program for Cancer Control, Ministry of Health and Welfare, Republic of Korea (1420190 (S.K.P.)), the Italian Association for Cancer Research (IG16933 (L.O.)) and German Cancer Aid (110837 (R.K.S.).

Lateral meningocele syndrome (LMS, OMIM%130720), also known as Lehman syndrome, is a very rare skeletal disorder with facial anomalies, hypotonia and meningocele-related neurologic dysfunction. The characteristic lateral meningoceles represent the severe end of the dural ectasia spectrum and are typically most severe in the lower spine. Facial features of LMS include hypertelorism and telecanthus, high arched eyebrows, ptosis, midfacial hypoplasia, micrognathia, high and narrow palate, low-set ears and a hypotonic appearance. Hyperextensibility, hernias and scoliosis reflect a connective tissue abnormality, and aortic dilation, a high-pitched nasal voice, wormian bones and osteolysis may be present. Lateral meningocele syndrome has phenotypic overlap with Hajdu–Cheney syndrome. We performed exome resequencing in five unrelated individuals with LMS and identified heterozygous truncating NOTCH3 mutations. In an additional unrelated individual Sanger sequencing revealed a deleterious variant in the same exon 33. In total, five novel de novo NOTCH3 mutations were identified in six unrelated patients. One had a 26 bp deletion (c.6461_6486del, p.G2154fsTer78), two carried the same single base pair insertion (c.6692_93insC, p.P2231fsTer11), and three individuals had a nonsense point mutation at c.6247A > T (pK2083*), c.6663C > G (p.Y2221*) or c.6732C > A, (p.Y2244*). All mutations cluster into the last coding exon, resulting in premature termination of the protein and truncation of the negative regulatory proline-glutamate-serine-threonine rich PEST domain. Our results suggest that mutant mRNA products escape nonsense mediated decay. The truncated NOTCH3 may cause gain-of-function through decreased clearance of the active intracellular product, resembling NOTCH2 mutations in the clinically related Hajdu–Cheney syndrome and contrasting the NOTCH3 missense mutations causing CADASIL. © 2014 Wiley Periodicals, Inc.

A dozen genes/regions have been confirmed as genetic risk factors for oral clefts in human association and linkage studies, and animal models argue even more genes may be involved. Genomic sequencing studies should identify specific causal variants and may reveal additional genes as influencing risk to oral clefts, which have a complex and heterogeneous etiology. We conducted a whole exome sequencing (WES) study to search for potentially causal variants using affected relatives drawn from multiplex cleft families. Two or three affected second, third, and higher degree relatives from 55 multiplex families were sequenced. We examined rare single nucleotide variants (SNVs) shared by affected relatives in 348 recognized candidate genes. Exact probabilities that affected relatives would share these rare variants were calculated, given pedigree structures, and corrected for the number of variants tested. Five novel and potentially damaging SNVs shared by affected distant relatives were found and confirmed by Sanger sequencing. One damaging SNV in CDH1, shared by three affected second cousins from a single family, attained statistical significance (P = 0.02 after correcting for multiple tests). Family-based designs such as the one used in this WES study offer important advantages for identifying genes likely to be causing complex and heterogeneous disorders.

We describe here the design and initial implementation of the eMERGE-PGx project. eMERGE-PGx, a partnership of the Electronic Medical Records and Genomics Network and the Pharmacogenomics Research Network, has three objectives: (i) to deploy PGRNseq, a next-generation sequencing platform assessing sequence variation in 84 proposed pharmacogenes, in nearly 9,000 patients likely to be prescribed drugs of interest in a 1- to 3-year time frame across several clinical sites; (ii) to integrate well-established clinically validated pharmacogenetic genotypes into the electronic health record with associated clinical decision support and to assess process and clinical outcomes of implementation; and (iii) to develop a repository of pharmacogenetic variants of unknown significance linked to a repository of electronic health record-based clinical phenotype data for ongoing pharmacogenomics discovery. We describe site-specific project implementation and anticipated products, including genetic variant and phenotype data repositories, novel variant association studies, clinical decision support modules, clinical and process outcomes, approaches to managing incidental findings, and patient and clinician education methods.

Autosomal recessive osteogenesis imperfecta (OI) accounts for 10% of all OI cases, and, currently, mutations in 10 genes (CRTAP, LEPRE1, PPIB, SERPINH1, FKBP10, SERPINF1, SP7, BMP1, TMEM38B, and WNT1) are known to be responsible for this form of the disease. PEDF is a secreted glycoprotein of the serpin superfamily that maintains bone homeostasis and regulates osteoid mineralization, and it is encoded by SERPINF1, currently associated with OI type VI (MIM 172860). Here, we report a consanguineous Brazilian family in which multiple individuals from at least 4 generations are affected with a severe form of OI, and we also report an unrelated individual from the same small city in Brazil with a similar but more severe phenotype. In both families the same homozygous SERPINF1 19-bp deletion was identified which is not known in the literature yet. We described intra- and interfamilial clinical and radiological phenotypic variability of OI type VI caused by the same homozygous SERPINF1 19-bp deletion and suggest a founder effect. Furthermore, the SERPINF1 genotypes/phenotypes reported so far in the literature are reviewed.

Microarray single-nucleotide polymorphism genotyping, combined with imputation of untyped variants, has been widely adopted as an efficient means to interrogate variation across the human genome. “Genomic coverage” is the total proportion of genomic variation captured by an array, either by direct observation or through an indirect means such as linkage disequilibrium or imputation. We have performed imputation-based genomic coverage assessments of eight current genotyping arrays that assay from ~0.3 to ~5 million variants. Coverage was determined separately in each of the four continental ancestry groups in the 1000 Genomes Project phase 1 release. We used the subset of 1000 Genomes variants present on each array to impute the remaining variants and assessed coverage based on correlation between imputed and observed allelic dosages. More than 75% of common variants (minor allele frequency > 0.05) are covered by all arrays in all groups except for African ancestry, and up to ~90% in all ancestries for the highest density arrays. In contrast, less than 40% of less common variants (0.01 < minor allele frequency < 0.05) are covered by low density arrays in all ancestries and 50–80% in high density arrays, depending on ancestry. We also calculated genome-wide power to detect variant-trait association in a case-control design, across varying sample sizes, effect sizes, and minor allele frequency ranges, and compare these array-based power estimates with a hypothetical array that would type all variants in 1000 Genomes. These imputation-based genomic coverage and power analyses are intended as a practical guide to researchers planning genetic studies.

This unit describes a technique for generating exome-enriched sequencing libraries using DNA extracted from formalin-fixed paraffin-embedded (FFPE) samples. Utilizing commercially available kits, we present a low-input FFPE workflow starting with 50 ng of DNA. This procedure includes a repair step to address damage caused by FFPE preservation that improves sequence quality. Subsequently, libraries undergo an in-solution-targeted selection for exons, followed by sequencing using the Illumina next-generation short-read sequencing platform. © 2017 by John Wiley & Sons, Inc.

Breast cancer risk is influenced by rare coding variants in susceptibility genes, such as BRCA1, and many common, mostly non-coding variants. However, much of the genetic contribution to breast cancer risk remains unknown. Here we report the results of a genome-wide association study of breast cancer in 122,977 cases and 105,974 controls of European ancestry and 14,068 cases and 13,104 controls of East Asian ancestry. We identified 65 new loci that are associated with overall breast cancer risk at P

Background: Common cancers develop through a multistep process often including inherited susceptibility. Collaboration among multiple institutions, and funding from multiple sources, has allowed the development of an inexpensive genotyping microarray, the OncoArray. The array includes a genome-wide backbone, comprising 230,000 SNPs tagging most common genetic variants, together with dense mapping of known susceptibility regions, rare variants from sequencing experiments, pharmacogenetic markers and cancer related traits. Methods: The OncoArray can be genotyped using a novel technology developed by Illumina to facilitate efficient genotyping. The consortium developed standard approaches for selecting SNPs for study, for quality control of markers and for ancestry analysis. The array was genotyped at selected sites and with prespecified replicate samples to permit evaluation of genotyping accuracy among centers and by ethnic background. Results: The OncoArray consortium genotyped 447,705 samples. A total of 494,763 SNPs passed quality control steps with a sample success rate of 97% of the samples. Participating sites performed ancestry analysis using a common set of markers and a scoring algorithm based on principal components analysis. Conclusions: Results from these analyses will enable researchers to identify new susceptibility loci, perform fine mapping of new or known loci associated with either single or multiple cancers, assess the degree of overlap in cancer causation and pleiotropic effects of loci that have been identified for disease-specific risk, and jointly model genetic, environmental and lifestyle related exposures. Impact: Ongoing analyses will shed light on etiology and risk assessment for many types of cancer. TRICL (Transdisciplinary Research for Cancer of Lung) and International Lung Cancer Consortium (ILCCO): National Institute of Health U19 CA148127-01 (C.I. Amos, J. Byun, Y. Li, X. Xiao, J. L. Forman, A. Risch, H. Bickeböller, M. Waldenberger, I. Brüske, J. Manz, P. Brennan, R. Hung, H. Shen, Y. Shi, A. Kamal, C. I. Amos, J.K. Field), Canadian Cancer Society Research Institute (no. 020214, R. Hung). DRIVE (Discovery, Biology, and Risk of Inherited Variants in Breast Cancer): National Institute of Health U19 CA148065 (D.J. Hazelett, D. F. Easton, S. Lindström, P. Kraft, J. Dennis, A. Dunning, K. Michailidou, L. Fachal, S. Benlloch, J. Cunningham, K. Butterbach, K. Offit, R. Schmutzler, L. Ottini, D. Vincent). CORECT (ColoRectal Transdisciplinary Study): National Institute of Health U19 CA148107; R01 CA81488, P30 CA014089 (S. Gruber, U. Peters, G. Casey). ELLIPSE (ELLIPSE, Elucidating Loci in Prostate Cancer Susceptibility): This work was support by U19 CA148537(F. Schumacher, S.A. Gayther, S.L. Schmit, C.K. Edlund, D.J. Hazelett, G.A. Coetzee, C. Haiman, S. Demetriades, D. Van Den Berg). FOCI (Transdisciplinary Cancer Genetic Association and Interacting Studies): National Institutes of Health U19 CA148112- 01 (T.A. Sellers, C. Phelan), R01-CA122443(B. Hicks) ) P50-CA116201, CA192393 (F. Couch), P50-CA136393(J. Cunningham, F. Couch), OCAC - CA-149429 (C. Phelan). , P30-CA15083 (B. Hicks), Cancer Research UK (C490/A8339(A. Antoniou, P. Pharoah, A. Dunning), C490/A16561(P. Pharoah, A. Dunning), C490/A10119 (A. Dunning, P. Pharoah), C490/A10124 (A. Dunning, P. Pharoah)). ASTERISK: a Hospital Clinical Research Program (PHRC) and supported by the Regional Council of Pays de la Loire, the Groupement des Entreprises Françaises dans la Lutte contre le Cancer (GEFLUC), the Association Anne de Bretagne Génétique and the Ligue Régionale Contre le Cancer (LRCC) (F. Bacot). DACHS: German Research Council (Deutsche Downloaded from cebp.aacrjournals.org on October 31, 2016. © 2016 American Association for Cancer Research. Author manuscripts have been peer reviewed and accepted for publication but have not yet been edited. Author Manuscript Published OnlineFirst on October 3, 2016; DOI: 10.1158/1055-9965.EPI-16-0106 Forschungsgemeinschaft, BR 1704/6-1, BR 1704/6-3, BR 1704/6-4, and CH 117/1-1, C.K. Edlund), and the German Federal Ministry of Education and Research (01KH0404 and 01ER0814). Prostate Lung Colon Ovarian Screening trial: Intramural Research Program of the Division of Cancer Epidemiology and Genetics and supported by contracts from the Division of Cancer Prevention, National Cancer Institute, NIH, DHHS. The GC-HBOC (German Consortium of Hereditary Breast and Ovarian Cancer) is supported by the German Cancer Aid (grant no 110837, coordinator: R. Schmutzler, Cologne). The PERSPECTIVE (J. Simard) project was supported by the Government of Canada through Genome Canada and the Canadian Institutes of Health Research (MOP-86727, L.E. Kelemen), the Ministère de l’Économie, Innovation et Exportation du Québec through Genome Québec, and thee Quebec Breast Cancer Foundation, including additional funding from the Canadian Breast Cancer Foundation and the National Cancer Institute USA through Breast Cancer Family Registry Cohort (#1UM1CA164920-01A1). European Union Framework Programme for Research and Innovation (MSCA-IF-2014-EF-656144)(L. Fachal). The Breast Cancer Research Foundation, P30 CA08748, Andrew Sabin Foundation (K. Offit). CIDR genotyping for the Oncoarray was conducted under contract 268201200008I to K.Doheny and through grant 1X01HG007491- 01 to C.I. Amos. This is the author accepted manuscript. The final version is available from American Association for Cancer Research via http://dx.doi.org/10.1158/1055-9965.EPI-16-0106

Karen Mohlke, Markku Laakso, Michael Boehnke and colleagues report the first application of the Illumina HumanExome Beadchip array, examining association with insulin and glycemic traits in 8,229 nondiabetic Finnish males from the population-based Metabolic Syndrome in Men (METSIM) study. They identify low-frequency coding variants at both known and newly associated loci with insulin processing and secretion.

DNA sample contamination is a serious problem in DNA sequencing studies and may result in systematic genotype misclassification and false positive associations. Although methods exist to detect and filter out cross-species contamination, few methods to detect within-species sample contamination are available. In this paper, we describe methods to identify within-species DNA sample contamination based on (1) a combination of sequencing reads and array-based genotype data, (2) sequence reads alone, and (3) array-based genotype data alone. Analysis of sequencing reads allows contamination detection after sequence data is generated but prior to variant calling; analysis of array-based genotype data allows contamination detection prior to generation of costly sequence data. Through a combination of analysis of in silico and experimentally contaminated samples, we show that our methods can reliably detect and estimate levels of contamination as low as 1%. We evaluate the impact of DNA contamination on genotype accuracy and propose effective strategies to screen for and prevent DNA contamination in sequencing studies.

When smokers were separated by their nicotinic receptor gene variants, those with the low-risk genotype responded equally well to pharmacological treatments, including both nicotine replacement and bupropion, and nonpharmacological therapies. Those with the high-risk genotype, as identified by DNA sequencing, responded only to pharmacological treatments. Clinicians advising patients on smoking cessation can suspect genetic risk on the basis of early onset of heavy smoking and direct those smokers specifically to pharmacological treatments.

Aims Nicotine dependence is a highly heritable disorder associated with severe medical morbidity and mortality. Recent meta-analyses have found novel genetic loci associated with cigarettes per day (CPD), a proxy for nicotine dependence. The aim of this paper is to evaluate the importance of phenotype definition (i.e. CPD versus Fagerstrom Test for Cigarette Dependence (FTCD) score as a measure of nicotine dependence) on genome-wide association studies of nicotine dependence. Design Genome-wide association study. Setting Community sample. Participants A total of 3365 subjects who had smoked at least one cigarette were selected from the Study of Addiction: Genetics and Environment (SAGE). Of the participants, 2267 were European Americans, 999 were African Americans. Measurements Nicotine dependence defined by FTCD score ≥4, CPD. Findings The genetic locus most strongly associated with nicotine dependence was rs1451240 on chromosome 8 in the region of CHRNB3 [odds ratio (OR) = 0.65, P = 2.4 × 10−8]. This association was further strengthened in a meta-analysis with a previously published data set (combined P = 6.7 × 10−16, total n = 4200). When CPD was used as an alternate phenotype, the association no longer reached genome-wide significance (β = −0.08, P = 0.0004). Conclusions Daily cigarette consumption and the Fagerstrom Test for Cigarette Dependence show different associations with polymorphisms in genetic loci.

Clonal mosaicism for large chromosomal anomalies (duplications, deletions and uniparental disomy) was detected using SNP microarray data from over 50,000 subjects recruited for genome-wide association studies. This detection method requires a relatively high frequency of cells (>5–10%) with the same abnormal karyotype (presumably of clonal origin) in the presence of normal cells. The frequency of detectable clonal mosaicism in peripheral blood is low (

Genome-wide association studies (GWAS) are a useful approach in the study of the genetic components of complex phenotypes. Aside from large cohorts, GWAS have generally been limited to the study of one or a few diseases or traits. The emergence of biobanks linked to electronic medical records (EMRs) allows the efficient re-use of genetic data to yield meaningful genotype-phenotype associations for multiple phenotypes or traits. Phase I of the electronic MEdical Records and GEnomics (eMERGE-I) Network is a National Human Genome Research Institute (NHGRI)-supported consortium composed of five sites to perform various genetic association studies using DNA repositories and EMR systems. Each eMERGE site has developed EMR-based algorithms to comprise a core set of fourteen phenotypes for extraction of study samples from each site’s DNA repository. Each eMERGE site selected samples for a specific phenotype, and these samples were genotyped at either the Broad Institute or at the Center for Inherited Disease Research (CIDR) using the Illumina Infinium BeadChip technology. In all, approximately 17,000 samples from across the five sites were genotyped. A unified quality control (QC) pipeline was developed by the eMERGE Genomics Working Group and used to ensure thorough cleaning of the data. This process includes examination of sample quality, marker quality, and various batch effects. Upon completion of the genotyping and QC analyses for each site’s primary study, the eMERGE Coordinating Center merged the datasets from all five sites. This larger merged dataset re-entered the established eMERGE QC pipeline. Based on lessons learned during the process, additional analyses and QC checkpoints were added to the pipeline to ensure proper merging. Here we explore the challenges associated with combining datasets from different genotyping centers and describe the expansion to the eMERGE QC pipeline for merged datasets. These additional steps will be useful as the eMERGE project expands to include additional sites in eMERGE-II and also serve as a starting point for investigators merging multiple genotype data sets accessible through the National Center for Biotechnology Information (NCBI) in the database of Genotypes and Phenotypes (dbGaP). Our experience demonstrates that merging multiple datasets after additional QC can be an efficient use of genotype data despite new challenges that appear in the process.

Ischemic stroke (IS) is among the leading causes of death in Western countries. There is a significant genetic component to IS susceptibility, especially among young adults. To date, research to identify genetic loci predisposing to stroke has met only with limited success. We performed a genome-wide association (GWA) analysis of early-onset IS to identify potential stroke susceptibility loci. The GWA analysis was conducted by genotyping 1 million SNPs in a biracial population of 889 IS cases and 927 controls, ages 15–49 years. Genotypes were imputed using the HapMap3 reference panel to provide 1.4 million SNPs for analysis. Logistic regression models adjusting for age, recruitment stages, and population structure were used to determine the association of IS with individual SNPs. Although no single SNP reached genome-wide significance (P < 5 × 10−8), we identified two SNPs in chromosome 2q23.3, rs2304556 (in FMNL2; P = 1.2 × 10−7) and rs1986743 (in ARL6IP6; P = 2.7 × 10−7), strongly associated with early-onset stroke. These data suggest that a novel locus on human chromosome 2q23.3 may be associated with IS susceptibility among young adults.

White blood cell count (WBC) is unique among identified inflammatory predictors of chronic disease in that it is routinely measured in asymptomatic patients in the course of routine patient care. We led a genome-wide association analysis to identify variants associated with WBC levels in 13,923 subjects in the electronic Medical Records and Genomics (eMERGE) Network. We identified two regions of interest that were each unique to subjects of genetically determined ancestry to the African continent (AA) or to the European continent (EA). WBC varies among different ancestry groups. Despite being ancestry specific, these regions were identifiable in the combined analysis. In AA subjects, the region surrounding the Duffy antigen/chemokine receptor gene (DARC) on 1q21 exhibited significant association (p value = 6.71e–55). These results validate the previously reported association between WBC and of the regulatory variant rs2814778 in the promoter region, which causes the Duffy negative phenotype (Fy−/−). A second missense variant (rs12075) is responsible for the two principal antigens, Fya and Fyb of the Duffy blood group system. The two variants, consisting of four alleles, act in concert to produce five antigens and subsequent phenotypes. We were able to identify the marginal and novel interaction effects of these two variants on WBC. In the EA subjects, we identified significantly associated SNPs tagging three separate genes in the 17q21 region: (1) GSDMA, (2) MED24, and (3) PSMD3. Variants in this region have been reported to be associated with WBC, neutrophil count, and inflammatory diseases including asthma and Crohn’s disease.