Whole genomic analysis reveals atypical non-homologous off-target large structural variants induced by CRISPR-Cas9-mediated genome editing.

CRISPR-Cas9 genome editing has promising therapeutic potential for genetic diseases and cancers, but safety could be a concern. Here we use whole genomic analysis by 10x linked-read sequencing and optical genome mapping to interrogate the genome integrity after editing and in comparison to four parental cell lines. In addition to the previously reported large structural variants at on-target sites, we identify heretofore unexpected large chromosomal deletions (91.2 and 136 Kb) at atypical non-homologous off-target sites without sequence similarity to the sgRNA in two edited lines. The observed large structural variants induced by CRISPR-Cas9 editing in dividing cells may result in pathogenic consequences and thus limit the usefulness of the CRISPR-Cas9 editing system for disease modeling and gene therapy. In this work, our whole genomic analysis may provide a valuable strategy to ensure genome integrity after genomic editing to minimize the risk of unintended effects in research and clinical applications.

Whole-exome sequencing reveals damaging gene variants associated with hypoalphalipoproteinemia.

Low levels of high density lipoprotein-cholesterol (HDL-C) are associated with an elevated risk of arteriosclerotic coronary heart disease. Heritability of HDL-C levels is high. In this research discovery study, we used whole-exome sequencing to identify damaging gene variants that may play significant roles in determining HDL-C levels. We studied 204 individuals with a mean HDL-C level of 27.8 ± 6.4 mg/dl (range: 4-36 mg/dl). Data were analyzed by statistical gene burden testing and by filtering against candidate gene lists. We found 120 occurrences of probably damaging variants (116 heterozygous; four homozygous) among 45 of 104 recognized HDL candidate genes. Those with the highest prevalence of damaging variants were ABCA1 (n = 20), STAB1 (n = 9), OSBPL1A (n = 8), CPS1 (n = 8), CD36 (n = 7), LRP1 (n = 6), ABCA8 (n = 6), GOT2 (n = 5), AMPD3 (n = 5), WWOX (n = 4), and IRS1 (n = 4). Binomial analysis for damaging missense or loss-of-function variants identified the ABCA1 and LDLR genes at genome-wide significance. In conclusion, whole-exome sequencing of individuals with low HDL-C showed the burden of damaging rare variants in the ABCA1 and LDLR genes is particularly high and revealed numerous occurrences in HDL candidate genes, including many genes identified in genome-wide association study reports. Many of these genes are involved in cancer biology, which accords with epidemiologic findings of the association of HDL deficiency with increased risk of cancer, thus presenting a new area of interest in HDL genomics.

Integrated genomic analyses of cutaneous T-cell lymphomas reveal the molecular bases for disease heterogeneity.

Cutaneous T-cell lymphomas (CTCLs) are a clinically heterogeneous collection of lymphomas of the skin-homing T cell. To identify molecular drivers of disease phenotypes, we assembled representative samples of CTCLs from patients with diverse disease subtypes and stages. Via DNA/RNA-sequencing, immunophenotyping, and ex vivo functional assays, we identified the landscape of putative driver genes, elucidated genetic relationships between CTCLs across disease stages, and inferred molecular subtypes in patients with stage-matched leukemic disease. Collectively, our analysis identified 86 putative driver genes, including 19 genes not previously implicated in this disease. Two mutations have never been described in any cancer. Functionally, multiple mutations augment T-cell receptor-dependent proliferation, highlighting the importance of this pathway in lymphomagenesis. To identify putative genetic causes of disease heterogeneity, we examined the distribution of driver genes across clinical cohorts. There are broad similarities across disease stages. Many driver genes are shared by mycosis fungoides (MF) and Sezary syndrome (SS). However, there are significantly more structural variants in leukemic disease, leading to highly recurrent deletions of putative tumor suppressors that are uncommon in early-stage skin-centered MF. For example, TP53 is deleted in 7% and 87% of MF and SS, respectively. In both human and mouse samples, PD1 mutations drive aggressive behavior. PD1 wild-type lymphomas show features of T-cell exhaustion. PD1 deletions are sufficient to reverse the exhaustion phenotype, promote a FOXM1-driven transcriptional signature, and predict significantly worse survival. Collectively, our findings clarify CTCL genetics and provide novel insights into pathways that drive diverse disease phenotypes.

Genetic profiles of 103,106 individuals in the Taiwan Biobank provide insights into the health and history of Han Chinese.

Personalized medical care focuses on prediction of disease risk and response to medications. To build the risk models, access to both large-scale genomic resources and human genetic studies is required. The Taiwan Biobank (TWB) has generated high-coverage, whole-genome sequencing data from 1492 individuals and genome-wide SNP data from 103,106 individuals of Han Chinese ancestry using custom SNP arrays. Principal components analysis of the genotyping data showed that the full range of Han Chinese genetic variation was found in the cohort. The arrays also include thousands of known functional variants, allowing for simultaneous ascertainment of Mendelian disease-causing mutations and variants that affect drug metabolism. We found that 21.2% of the population are mutation carriers of autosomal recessive diseases, 3.1% have mutations in cancer-predisposing genes, and 87.3% carry variants that affect drug response. We highlight how TWB data provide insight into both population history and disease burden, while showing how widespread genetic testing can be used to improve clinical care.

A Large-Scale Association Study Detects Novel Rare Variants, Risk Genes, Functional Elements, and Polygenic Architecture of Prostate Cancer Susceptibility.

To identify rare variants associated with prostate cancer susceptibility and better characterize the mechanisms and cumulative disease risk associated with common risk variants, we conducted an integrated study of prostate cancer genetic etiology in two cohorts using custom genotyping microarrays, large imputation reference panels, and functional annotation approaches. Specifically, 11,984 men (6,196 prostate cancer cases and 5,788 controls) of European ancestry from Northern California Kaiser Permanente were genotyped and meta-analyzed with 196,269 men of European ancestry (7,917 prostate cancer cases and 188,352 controls) from the UK Biobank. Three novel loci, including two rare variants (European ancestry minor allele frequency < 0.01, at 3p21.31 and 8p12), were significant genome wide in a meta-analysis. Gene-based rare variant tests implicated a known prostate cancer gene (HOXB13), as well as a novel candidate gene (ILDR1), which encodes a receptor highly expressed in prostate tissue and is related to the B7/CD28 family of T-cell immune checkpoint markers. Haplotypic patterns of long-range linkage disequilibrium were observed for rare genetic variants at HOXB13 and other loci, reflecting their evolutionary history. In addition, a polygenic risk score (PRS) of 188 prostate cancer variants was strongly associated with risk (90th vs. 40th-60th percentile OR = 2.62, P = 2.55 × 10-191). Many of the 188 variants exhibited functional signatures of gene expression regulation or transcription factor binding, including a 6-fold difference in log-probability of androgen receptor binding at the variant rs2680708 (17q22). Rare variant and PRS associations, with concomitant functional interpretation of risk mechanisms, can help clarify the full genetic architecture of prostate cancer and other complex traits. SIGNIFICANCE: This study maps the biological relationships between diverse risk factors for prostate cancer, integrating different functional datasets to interpret and model genome-wide data from over 200,000 men with and without prostate cancer.See related commentary by Lachance, p. 1637.

Three patients with homozygous familial hypercholesterolemia: Genomic sequencing and kindred analysis.

BACKGROUND: Homozygous Familial Hypercholesterolemia (HoFH) is an inherited recessive condition associated with extremely high levels of low-density lipoprotein (LDL) cholesterol in affected individuals. It is usually caused by homozygous or compound heterozygous functional mutations in the LDL receptor (LDLR). A number of mutations causing FH have been reported in literature and such genetic heterogeneity presents great challenges for disease diagnosis. OBJECTIVE: We aim to determine the likely genetic defects responsible for three cases of pediatric HoFH in two kindreds. METHODS: We applied whole exome sequencing (WES) on the two probands to determine the likely functional variants among candidate FH genes. We additionally applied 10x Genomics (10xG) Linked-Reads whole genome sequencing (WGS) on one of the kindreds to identify potentially deleterious structural variants (SVs) underlying HoFH. A PCR-based screening assay was also established to detect the LDLR structural variant in a cohort of 641 patients with elevated LDL. RESULTS: In the Caucasian kindred, the FH homozygosity can be attributed to two compound heterozygous LDLR damaging variants, an exon 12 p.G592E missense mutation and a novel 3kb exon 1 deletion. By analyzing the 10xG phased data, we ascertained that this deletion allele was most likely to have originated from a Russian ancestor. In the Mexican kindred, the strikingly elevated LDL cholesterol level can be attributed to a homozygous frameshift LDLR variant p.E113fs. CONCLUSIONS: While the application of WES can provide a cost-effective way of identifying the genetic causes of FH, it often lacks sensitivity for detecting structural variants. Our finding of the LDLR exon 1 deletion highlights the broader utility of Linked-Read WGS in detecting SVs in the clinical setting, especially when HoFH patients remain undiagnosed after WES.

Targeted Genomic Profiling of Acral Melanoma.

BACKGROUND: Acral melanoma is a rare type of melanoma that affects world populations irrespective of skin color and has worse survival than other cutaneous melanomas. It has relatively few single nucleotide mutations without the UV signature of cutaneous melanomas, but instead has a genetic landscape characterized by structural rearrangements and amplifications. BRAF mutations are less common than in other cutaneous melanomas, and knowledge about alternative therapeutic targets is incomplete. METHODS: To identify alternative therapeutic targets, we performed targeted deep-sequencing on 122 acral melanomas. We confirmed the loss of the tumor suppressors p16 and NF1 by immunohistochemistry in select cases. RESULTS: In addition to BRAF (21.3%), NRAS (27.9%), and KIT (11.5%) mutations, we identified a broad array of MAPK pathway activating alterations, including fusions of BRAF (2.5%), NTRK3 (2.5%), ALK (0.8%), and PRKCA (0.8%), which can be targeted by available inhibitors. Inactivation of NF1 occurred in 18 cases (14.8%). Inactivation of the NF1 cooperating factor SPRED1 occurred in eight cases (6.6%) as an alternative mechanism of disrupting the negative regulation of RAS. Amplifications recurrently affected narrow loci containing PAK1 and GAB2 (n = 27, 22.1%), CDK4 (n = 27, 22.1%), CCND1 (n = 24, 19.7%), EP300 (n = 20, 16.4%), YAP1 (n = 15, 12.3%), MDM2 (n = 13, 10.7%), and TERT (n = 13, 10.7%) providing additional and possibly complementary therapeutic targets. Acral melanomas with BRAFV600E mutations harbored fewer genomic amplifications and were more common in patients with European ancestry. CONCLUSION: Our findings support a new, molecularly based subclassification of acral melanoma with potential therapeutic implications: BRAFV600E mutant acral melanomas with characteristics similar to nonacral melanomas that could benefit from BRAF inhibitor therapy, and non-BRAFV600E mutant acral melanomas. Acral melanomas without BRAFV600E mutations harbor a broad array of therapeutically relevant alterations. Expanded molecular profiling would increase the detection of potentially targetable alterations for this subtype of acral melanoma.

The Clinical Sequencing Evidence-Generating Research Consortium: Integrating Genomic Sequencing in Diverse and Medically Underserved Populations.

The Clinical Sequencing Evidence-Generating Research (CSER) consortium, now in its second funding cycle, is investigating the effectiveness of integrating genomic (exome or genome) sequencing into the clinical care of diverse and medically underserved individuals in a variety of healthcare settings and disease states. The consortium comprises a coordinating center, six funded extramural clinical projects, and an ongoing National Human Genome Research Institute (NHGRI) intramural project. Collectively, these projects aim to enroll and sequence over 6,100 participants in four years. At least 60% of participants will be of non-European ancestry or from underserved settings, with the goal of diversifying the populations that are providing an evidence base for genomic medicine. Five of the six clinical projects are enrolling pediatric patients with various phenotypes. One of these five projects is also enrolling couples whose fetus has a structural anomaly, and the sixth project is enrolling adults at risk for hereditary cancer. The ongoing NHGRI intramural project has enrolled primarily healthy adults. Goals of the consortium include assessing the clinical utility of genomic sequencing, exploring medical follow up and cascade testing of relatives, and evaluating patient-provider-laboratory level interactions that influence the use of this technology. The findings from the CSER consortium will offer patients, healthcare systems, and policymakers a clearer understanding of the opportunities and challenges of providing genomic medicine in diverse populations and settings, and contribute evidence toward developing best practices for the delivery of clinically useful and cost-effective genomic sequencing in diverse healthcare settings.

Genome-wide association study of prostate-specific antigen levels identifies novel loci independent of prostate cancer.

Prostate-specific antigen (PSA) levels have been used for detection and surveillance of prostate cancer (PCa). However, factors other than PCa-such as genetics-can impact PSA. Here we present findings from a genome-wide association study (GWAS) of PSA in 28,503 Kaiser Permanente whites and 17,428 men from replication cohorts. We detect 40 genome-wide significant (P<5 × 10-8) single-nucleotide polymorphisms (SNPs): 19 novel, 15 previously identified for PSA (14 of which were also PCa-associated), and 6 previously identified for PCa only. Further analysis incorporating PCa cases suggests that at least half of the 40 SNPs are PSA-associated independent of PCa. The 40 SNPs explain 9.5% of PSA variation in non-Hispanic whites, and the remaining GWAS SNPs explain an additional 31.7%; this percentage is higher in younger men, supporting the genetic basis of PSA levels. These findings provide important information about genetic markers for PSA that may improve PCa screening, thereby reducing over-diagnosis and over-treatment.

Arrhythmogenic Right Ventricular Cardiomyopathy Caused by a Novel Frameshift Mutation.

Arrhythmogenic right ventricular cardiomyopathy is a rare cardiomyopathy that might be asymptomatic or symptomatic, causing palpations or syncope, and might lead to sudden cardiac death. It is recommended that physical exertion be reduced. It is also recommended that those with syncope and ventricular tachycardia/ventricular fibrillation have an implantable cardioverter-defibrillator placed. ß-Blockers, antiarrhythmic drugs, and radiofrequency ablation should be used to control the ventricular arrhythmia burden in arrhythmogenic right ventricular cardiomyopathy.

A large multiethnic genome-wide association study of prostate cancer identifies novel risk variants and substantial ethnic differences.

UNLABELLED: A genome-wide association study (GWAS) of prostate cancer in Kaiser Permanente health plan members (7,783 cases, 38,595 controls; 80.3% non-Hispanic white, 4.9% African-American, 7.0% East Asian, and 7.8% Latino) revealed a new independent risk indel rs4646284 at the previously identified locus 6q25.3 that replicated in PEGASUS (N = 7,539) and the Multiethnic Cohort (N = 4,679) with an overall P = 1.0 × 10(-19) (OR, 1.18). Across the 6q25.3 locus, rs4646284 exhibited the strongest association with expression of SLC22A1 (P = 1.3 × 10(-23)) and SLC22A3 (P = 3.2 × 10(-52)). At the known 19q13.33 locus, rs2659124 (P = 1.3 × 10(-13); OR, 1.18) nominally replicated in PEGASUS. A risk score of 105 known risk SNPs was strongly associated with prostate cancer (P < 1.0 × 10(-8)). Comparing the highest to lowest risk score deciles, the OR was 6.22 for non-Hispanic whites, 5.82 for Latinos, 3.77 for African-Americans, and 3.38 for East Asians. In non-Hispanic whites, the 105 risk SNPs explained approximately 7.6% of disease heritability. The entire GWAS array explained approximately 33.4% of heritability, with a 4.3-fold enrichment within DNaseI hypersensitivity sites (P = 0.004). SIGNIFICANCE: Taken together, our findings of independent risk variants, ethnic variation in existing SNP replication, and remaining unexplained heritability have important implications for further clarifying the genetic risk of prostate cancer. Our findings also suggest that there may be much promise in evaluating understudied variation, such as indels and ethnically diverse populations.

The genetics of splicing in neuroblastoma.

UNLABELLED: Regulation of mRNA splicing, a critical and tightly regulated cellular function, underlies the majority of proteomic diversity and is frequently disrupted in disease. Using an integrative genomics approach, we combined both genomic data and exon-level transcriptome data in two somatic tissues (cerebella and peripheral ganglia) from a transgenic mouse model of neuroblastoma, a tumor that arises from the peripheral neural crest. Here, we describe splicing quantitative trait loci associated with differential splicing across the genome that we use to identify genes with previously unknown functions within the splicing pathway and to define de novo intronic splicing motifs that influence splicing from hundreds of bases away. Our results show that these splicing motifs represent sites for functional recurrent mutations and highlight novel candidate genes in human cancers, including childhood neuroblastoma. SIGNIFICANCE: Somatic mutations with predictable downstream effects are largely relegated to coding regions, which comprise less than 2% of the human genome. Using an unbiased in vivo analysis of a mouse model of neuroblastoma, we have identified intronic splicing motifs that translate into sites for recurrent somatic mutations in human cancers.

Targeted deep sequencing reveals no definitive evidence for somatic mosaicism in atrial fibrillation.

BACKGROUND: Studies of ≤15 atrial fibrillation (AF) patients have identified atrial-specific mutations within connexin genes, suggesting that somatic mutations may account for sporadic cases of the arrhythmia. We sought to identify atrial somatic mutations among patients with and without AF using targeted deep next-generation sequencing of 560 genes, including genetic culprits implicated in AF, the Mendelian cardiomyopathies and channelopathies, and all ion channels within the genome. METHODS AND RESULTS: Targeted gene capture and next-generation sequencing were performed on DNA from lymphocytes and left atrial appendages of 34 patients (25 with AF). Twenty AF patients had undergone cardiac surgery exclusively for pulmonary vein isolation and 17 had no structural heart disease. Sequence alignment and variant calling were performed for each atrial-lymphocyte pair using the Burrows-Wheeler Aligner, the Genome Analysis Toolkit, and MuTect packages. Next-generation sequencing yielded a median 265-fold coverage depth (interquartile range, 64-369). Comparison of the 3 million base pairs from each atrial-lymphocyte pair revealed a single potential somatic missense mutation in 3 AF patients and 2 in a single control (12 versus 11%; P=1). All potential discordant variants had low allelic fractions (range, 2.3%-7.3%) and none were detected with conventional sequencing. CONCLUSIONS: Using high-depth next-generation sequencing and state-of-the art somatic mutation calling approaches, no pathogenic atrial somatic mutations could be confirmed among 25 AF patients in a comprehensive cardiac arrhythmia genetic panel. These findings indicate that atrial-specific mutations are rare and that somatic mosaicism is unlikely to exert a prominent role in AF pathogenesis.

Novel mutations in PXDN cause microphthalmia and anterior segment dysgenesis.

We used exome sequencing to study a non-consanguineous family with two children who had anterior segment dysgenesis, sclerocornea, microphthalmia, hypotonia and developmental delays. Sanger sequencing verified two Peroxidasin (PXDN) mutations in both sibs--a maternally inherited, nonsense mutation, c.1021C>T predicting p.(Arg341*), and a paternally inherited, 23-basepair deletion causing a frameshift and premature protein truncation, c.2375_2397del23, predicting p.(Leu792Hisfs*67). We re-examined exome data from 20 other patients with structural eye defects and identified two additional PXDN mutations in a sporadic male with bilateral microphthalmia, cataracts and anterior segment dysgenesis--a maternally inherited, frameshift mutation, c.1192delT, predicting p.(Tyr398Thrfs*40) and a paternally inherited, missense substitution that was predicted to be deleterious, c.947 A>C, predicting p.(Gln316Pro). Mutations in PXDN were previously reported in three families with congenital cataracts, microcornea, sclerocornea and developmental glaucoma. The gene is expressed in corneal epithelium and is secreted into the extracellular matrix. Defective peroxidasin has been shown to impair sulfilimine bond formation in collagen IV, a constituent of the basement membrane, implying that the eye defects result because of loss of basement membrane integrity in the developing eye. Our finding of a broader phenotype than previously appreciated for PXDN mutations is typical for exome-sequencing studies, which have proven to be highly effective for mutation detection in patients with atypical presentations. We conclude that PXDN sequencing should be considered in microphthalmia with anterior segment dysgenesis.

Expression quantitative trait loci and receptor pharmacology implicate Arg1 and the GABA-A receptor as therapeutic targets in neuroblastoma.

The development of targeted therapeutics for neuroblastoma, the third most common tumor in children, has been limited by a poor understanding of growth signaling mechanisms unique to the peripheral nerve precursors from which tumors arise. In this study, we combined genetics with gene-expression analysis in the peripheral sympathetic nervous system to implicate arginase 1 and GABA signaling in tumor formation in vivo. In human neuroblastoma cells, either blockade of ARG1 or benzodiazepine-mediated activation of GABA-A receptors induced apoptosis and inhibited mitogenic signaling through AKT and MAPK. These results suggest that ARG1 and GABA influence both neural development and neuroblastoma and that benzodiazepines in clinical use may have potential applications for neuroblastoma therapy.

Effects of the absence of apolipoprotein e on lipoproteins, neurocognitive function, and retinal function.

IMPORTANCE: The identification of a patient with a rare form of severe dysbetalipoproteinemia allowed the study of the consequences of total absence of apolipoprotein E (apoE). OBJECTIVES: To discover the molecular basis of this rare disorder and to determine the effects of complete absence of apoE on neurocognitive and visual function and on lipoprotein metabolism. DESIGN, SETTING, AND PARTICIPANTS: Whole-exome sequencing was performed on the patient's DNA. He underwent detailed neurological and visual function testing and lipoprotein analysis. Lipoprotein analysis was also performed in the Cardiovascular Research Institute, University of California, San Francisco, on blood samples from the proband's mother, wife, 2 daughters, and normolipidemic control participants. MAIN OUTCOME MEASURES: Whole-exome sequencing, lipoprotein analysis, and neurocognitive function. RESULTS: The patient was homozygous for an ablative APOE frameshift mutation (c.291del, p.E97fs). No other mutations likely to contribute to the phenotype were discovered, with the possible exception of two, in ABCC2 (p.I670T) and LIPC (p.G137R). Despite complete absence of apoE, he had normal vision, exhibited normal cognitive, neurological, and retinal function, had normal findings on brain magnetic resonance imaging, and had normal cerebrospinal fluid levels of ß-amyloid and tau proteins. He had no significant symptoms of cardiovascular disease except a suggestion of myocardial ischemia on treadmill testing and mild atherosclerosis noted on carotid ultrasonography. He had exceptionally high cholesterol content (760 mg/dL; to convert to millimoles per liter, multiply by 0.0259) and a high cholesterol to triglycerides ratio (1.52) in very low-density lipoproteins with elevated levels of small-diameter high-density lipoproteins, including high levels of prebeta-1 high-density lipoprotein. Intermediate-density lipoproteins, low-density lipoproteins, and very low-density lipoproteins contained elevated apoA-I and apoA-IV levels. The patient's apoC-III and apoC-IV levels were decreased in very low-density lipoproteins. Electron microscopy revealed large lamellar particles having electron-opaque cores attached to electron-lucent zones in intermediate-density and low-density lipoproteins. Low-density lipoprotein particle diameters were distributed bimodally. CONCLUSIONS AND RELEVANCE: Despite a profound effect on lipoprotein metabolism, detailed neurocognitive and retinal studies failed to demonstrate any defects. This suggests that functions of apoE in the brain and eye are not essential or that redundant mechanisms exist whereby its role can be fulfilled. Targeted knockdown of apoE in the central nervous system might be a therapeutic modality in neurodegenerative disorders.

Organic cation transporter variation and response to smoking cessation therapies.

INTRODUCTION: We evaluated chr6q25.3 organic cation transporter gene (SLC22A1, SLC22A2, SLC22A3) variation and response to smoking cessation therapies. The corresponding proteins are low-affinity transporters of choline, acetylcholine and monoamines, and smoking cessation pharmacotherapies expressed in multiple tissues. METHODS: We selected 7 common polymorphisms for mega-regression analysis. We assessed additive model association of polymorphisms with 7-day point prevalence abstinence overall and by assigned pharmacotherapy at end of treatment and at 6 months among European-ancestry participants of 7 randomized controlled trials adjusted for demographic, population genetic, and trial covariates. RESULTS: Initial results were obtained in 6 trials with 1,839 participants. Nominally statistically significant associations of 2 SLC22A2 polymorphisms were observed: (1) with rs316019 at 6 months, overall ([c.808T>G; p.Ser270Ala], OR = 1.306, 95% CI = 1.034-1.649, p = .025), and among those randomized to nicotine replacement therapy (NRT) (OR = 1.784, 95% CI = 1.072-2.970, p = .026); and (2) with rs316006 (c.1502-529A>T) among those randomized to varenicline (OR = 1.420, 95% CI = 1.038-1.944, p = .028, OR = 1.362, 95% CI = 1.001-1.853, p = .04) at end of treatment and 6 months. Individuals randomized to NRT from a seventh trial were genotyped for rs316019; rs316019 was associated with a nominally statistically significant effect on abstinence overall at 6 months among 2,233 participants (OR = 1.249, 95% CI = 1.007-1.550, p = .043). CONCLUSIONS: The functional OCT2 Ser270Ala polymorphism is nominally statistically significantly associated with abstinence among European-ancestry treatment-seeking smokers after adjustments for pharmacotherapy, demographics, population genetics, and without adjustment for multiple testing of 7 SNPs. Replication of these preliminary findings in additional randomized controlled trials of smoking cessation therapies and from multiple continental populations would describe another pharmacogenetic role for SLC22A2/OCT2.

Finished sequence and assembly of the DUF1220-rich 1q21 region using a haploid human genome.

BACKGROUND: Although the reference human genome sequence was declared finished in 2003, some regions of the genome remain incomplete due to their complex architecture. One such region, 1q21.1-q21.2, is of increasing interest due to its relevance to human disease and evolution. Elucidation of the exact variants behind these associations has been hampered by the repetitive nature of the region and its incomplete assembly. This region also contains 238 of the 270 human DUF1220 protein domains, which are implicated in human brain evolution and neurodevelopment. Additionally, examinations of this protein domain have been challenging due to the incomplete 1q21 build. To address these problems, a single-haplotype hydatidiform mole BAC library (CHORI-17) was used to produce the first complete sequence of the 1q21.1-q21.2 region. RESULTS: We found and addressed several inaccuracies in the GRCh37sequence of the 1q21 region on large and small scales, including genomic rearrangements and inversions, and incorrect gene copy number estimates and assemblies. The DUF1220-encoding NBPF genes required the most corrections, with 3 genes removed, 2 genes reassigned to the 1p11.2 region, 8 genes requiring assembly corrections for DUF1220 domains (~91 DUF1220 domains were misassigned), and multiple instances of nucleotide changes that reassigned the domain to a different DUF1220 subtype. These corrections resulted in an overall increase in DUF1220 copy number, yielding a haploid total of 289 copies. Approximately 20 of these new DUF1220 copies were the result of a segmental duplication from 1q21.2 to 1p11.2 that included two NBPF genes. Interestingly, this duplication may have been the catalyst for the evolutionarily important human lineage-specific chromosome 1 pericentric inversion. CONCLUSIONS: Through the hydatidiform mole genome sequencing effort, the 1q21.1-q21.2 region is complete and misassemblies involving inter- and intra-region duplications have been resolved. The availability of this single haploid sequence path will aid in the investigation of many genetic diseases linked to 1q21, including several associated with DUF1220 copy number variations. Finally, the corrected sequence identified a recent segmental duplication that added 20 additional DUF1220 copies to the human genome, and may have facilitated the chromosome 1 pericentric inversion that is among the most notable human-specific genomic landmarks.

Application of next generation sequencing to CEPH cell lines to discover variants associated with FDA approved chemotherapeutics.

BACKGROUND: The goal of this study was to perform candidate gene association with cytotoxicity of chemotherapeutics in cell line models through resequencing and discovery of rare and low frequency variants along with common variations. Here, an association study of cytotoxicity response to 30 FDA approved drugs was conducted and we applied next generation targeted sequencing technology to discover variants from 103 candidate genes in 95 lymphoblastoid cell lines from 14 CEPH pedigrees. In this article, we called variants across 95 cell lines and performed association analysis for cytotoxic response using the Family Based Association Testing method and software. RESULTS: We called 2281 variable SNP genotypes across the 103 genes for these cell lines and identified three genes of significant association within this marker set. Specifically, ATP-binding cassette, sub-family C, member 5 (ABCC5), metallothionein 1A (MT1A) and NAD(P)H dehydrogenase quinone1 (NQO1) were significantly associated with oxaliplatin drug response. The significant SNP on NQO1 (rs1800566) has been linked with poor survival rates in patients with non-small cell lung cancer treated with cisplatin (which belongs to the same class of drugs as oxaliplatin). A SNP (rs1846692) near the 5' region of MT1A was associated with arsenic trioxide. CONCLUSIONS: The results from this study are promising and this serves as a proof-of-principle demonstration of the use of sequencing data in the cytotoxicity models of human cell lines. With increased sample sizes, such studies will be a fast and powerful way to associate common and rare variants with drug response; while overcoming the cost and time limitations to recruit cohorts for association study.

A genome-wide investigation of copy number variation in patients with sporadic brain arteriovenous malformation.

BACKGROUND: Brain arteriovenous malformations (BAVM) are clusters of abnormal blood vessels, with shunting of blood from the arterial to venous circulation and a high risk of rupture and intracranial hemorrhage. Most BAVMs are sporadic, but also occur in patients with Hereditary Hemorrhagic Telangiectasia, a Mendelian disorder caused by mutations in genes in the transforming growth factor beta (TGFß) signaling pathway. METHODS: To investigate whether copy number variations (CNVs) contribute to risk of sporadic BAVM, we performed a genome-wide association study in 371 sporadic BAVM cases and 563 healthy controls, all Caucasian. Cases and controls were genotyped using the Affymetrix 6.0 array. CNVs were called using the PennCNV and Birdsuite algorithms and analyzed via segment-based and gene-based approaches. Common and rare CNVs were evaluated for association with BAVM. RESULTS: A CNV region on 1p36.13, containing the neuroblastoma breakpoint family, member 1 gene (NBPF1), was significantly enriched with duplications in BAVM cases compared to controls (P = 2.2×10(-9)); NBPF1 was also significantly associated with BAVM in gene-based analysis using both PennCNV and Birdsuite. We experimentally validated the 1p36.13 duplication; however, the association did not replicate in an independent cohort of 184 sporadic BAVM cases and 182 controls (OR = 0.81, P = 0.8). Rare CNV analysis did not identify genes significantly associated with BAVM. CONCLUSION: We did not identify common CNVs associated with sporadic BAVM that replicated in an independent cohort. Replication in larger cohorts is required to elucidate the possible role of common or rare CNVs in BAVM pathogenesis.