Multisite verification of the accuracy of a multi-gene next generation sequencing panel for detection of mutations and copy number alterations in solid tumours.

Molecular variants including single nucleotide variants (SNVs), copy number variants (CNVs) and fusions can be detected in the clinical setting using deep targeted sequencing. These assays support low limits of detection using little genomic input material. They are gaining in popularity in clinical laboratories, where sample volumes are limited, and low variant allele fractions may be present. However, data on reproducibility between laboratories is limited. Using a ring study, we evaluated the performance of 7 Ontario laboratories using targeted sequencing panels. All laboratories analysed a series of control and clinical samples for SNVs/CNVs and gene fusions. High concordance was observed across laboratories for measured CNVs and SNVs. Over 97% of SNV calls in clinical samples were detected by all laboratories. Whilst only a single CNV was detected in the clinical samples tested, all laboratories were able to reproducibly report both the variant and copy number. Concordance for information derived from RNA was lower than observed for DNA, due largely to decreased quality metrics associated with the RNA components of the assay, suggesting that the RNA portions of comprehensive NGS assays may be more vulnerable to variations in approach and workflow. Overall the results of this study support the use of the OFA for targeted sequencing for testing of clinical samples and suggest specific internal quality metrics that can be reliable indicators of assay failure. While we believe this evidence can be interpreted to support deep targeted sequencing in general, additional studies should be performed to confirm this.

Whole-exome sequencing reveals a rapid change in the frequency of rare functional variants in a founding population of humans.

Whole-exome or gene targeted resequencing in hundreds to thousands of individuals has shown that the majority of genetic variants are at low frequency in human populations. Rare variants are enriched for functional mutations and are expected to explain an important fraction of the genetic etiology of human disease, therefore having a potential medical interest. In this work, we analyze the whole-exome sequences of French-Canadian individuals, a founder population with a unique demographic history that includes an original population bottleneck less than 20 generations ago, followed by a demographic explosion, and the whole exomes of French individuals sampled from France. We show that in less than 20 generations of genetic isolation from the French population, the genetic pool of French-Canadians shows reduced levels of diversity, higher homozygosity, and an excess of rare variants with low variant sharing with Europeans. Furthermore, the French-Canadian population contains a larger proportion of putatively damaging functional variants, which could partially explain the increased incidence of genetic disease in the province. Our results highlight the impact of population demography on genetic fitness and the contribution of rare variants to the human genetic variation landscape, emphasizing the need for deep cataloguing of genetic variants by resequencing worldwide human populations in order to truly assess disease risk.

Rare allelic forms of PRDM9 associated with childhood leukemogenesis.

One of the most rapidly evolving genes in humans, PRDM9, is a key determinant of the distribution of meiotic recombination events. Mutations in this meiotic-specific gene have previously been associated with male infertility in humans and recent studies suggest that PRDM9 may be involved in pathological genomic rearrangements. In studying genomes from families with children affected by B-cell precursor acute lymphoblastic leukemia (B-ALL), we characterized meiotic recombination patterns within a family with two siblings having hyperdiploid childhood B-ALL and observed unusual localization of maternal recombination events. The mother of the family carries a rare PRDM9 allele, potentially explaining the unusual patterns found. From exomes sequenced in 44 additional parents of children affected with B-ALL, we discovered a substantial and significant excess of rare allelic forms of PRDM9. The rare PRDM9 alleles are transmitted to the affected children in half the cases; nonetheless there remains a significant excess of rare alleles among patients relative to controls. We successfully replicated this latter observation in an independent cohort of 50 children with B-ALL, where we found an excess of rare PRDM9 alleles in aneuploid and infant B-ALL patients. PRDM9 variability in humans is thought to influence genomic instability, and these data support a potential role for PRDM9 variation in risk of acquiring aneuploidies or genomic rearrangements associated with childhood leukemogenesis.

Identification of functional DNA variants in the constitutive promoter region of MDM2.

Although mutations in the oncoprotein murine double minute 2 (MDM2) are rare, MDM2 gene overexpression has been observed in several human tumors. Given that even modest changes in MDM2 levels might influence the p53 tumor suppressor signaling pathway, we postulated that sequence variation in the promoter region of MDM2 could lead to disregulated expression and variation in gene dosage. Two promoters have been reported for MDM2; an internal promoter (P2), which is located near the end of intron 1 and is p53-responsive, and an upstream constitutive promoter (P1), which is p53-independent. Both promoter regions contain DNA variants that could influence the expression levels of MDM2, including the well-studied single nucleotide polymorphism (SNP) SNP309, which is located in the promoter P2; i.e., upstream of exon 2. In this report, we screened the promoter P1 for DNA variants and assessed the functional impact of the corresponding SNPs. Using the dbSNP database and genotyping validation in individuals of European descent, we identified three common SNPs (-1494 G > A; indel 40 bp; and -182 C > G). Three major promoter haplotypes were inferred by using these three promoter SNPs together with rs2279744 (SNP309). Following subcloning into a gene reporter system, we found that two of the haplotypes significantly influenced MDM2 promoter activity in a haplotype-specific manner. Site-directed mutagenesis experiments indicated that the 40 bp insertion/deletion variation is causing the observed allelic promoter activity. This study suggests that part of the variability in the MDM2 expression levels could be explained by allelic p53-independent P1 promoter activity.

Functional analysis of promoter variants in KU70 and their role in cancer susceptibility.

KU70 is involved in the DNA double-strand break repair pathway, which plays a critical role in maintaining genomic stability and preventing cancer. Genetic variation within the KU70 gene has been shown to be associated with increased risk of several types of cancers including breast cancer. Here, we used gene reporter and gel shift assays combined with site-directed mutagenesis to characterize genetic variation within the KU70 proximal promoter region and investigate the putative functional role of regulatory variation and altered KU70 expression in modulating an individual's susceptibility to disease. We show that the variant rs2267437C>G significantly influences KU70 transcriptional activity in an allele- specific manner and alters DNA-protein complex formation in breast cancer cell lines. Our data provide a possible molecular explanation for the associations observed between the KU70 regulatory variant rs2267437 and breast cancer risk.

ALG: automated genotype calling of Luminex assays.

Single nucleotide polymorphisms (SNPs) are the most commonly used polymorphic markers in genetics studies. Among the different platforms for SNP genotyping, Luminex is one of the less exploited mainly due to the lack of a robust (semi-automated and replicable) freely available genotype calling software. Here we describe a clustering algorithm that provides automated SNP calls for Luminex genotyping assays. We genotyped 3 SNPs in a cohort of 330 childhood leukemia patients, 200 parents of patient and 325 healthy individuals and used the Automated Luminex Genotyping (ALG) algorithm for SNP calling. ALG genotypes were called twice to test for reproducibility and were compared to sequencing data to test for accuracy. Globally, this analysis demonstrates the accuracy (99.6%) of the method, its reproducibility (99.8%) and the low level of no genotyping calls (3.4%). The high efficiency of the method proves that ALG is a suitable alternative to the current commercial software. ALG is semi-automated, and provides numerical measures of confidence for each SNP called, as well as an effective graphical plot. Moreover ALG can be used either through a graphical user interface, requiring no specific informatics knowledge, or through command line with access to the open source code. The ALG software has been implemented in R and is freely available for non-commercial use either at http://alg.sourceforge.net or by request to mathieu.bourgey@umontreal.ca.

Allele-specific chromatin remodeling in the ZPBP2/GSDMB/ORMDL3 locus associated with the risk of asthma and autoimmune disease.

Common SNPs in the chromosome 17q12-q21 region alter the risk for asthma, type 1 diabetes, primary biliary cirrhosis, and Crohn disease. Previous reports by us and others have linked the disease-associated genetic variants with changes in expression of GSDMB and ORMDL3 transcripts in human lymphoblastoid cell lines (LCLs). The variants also alter regulation of other transcripts, and this domain-wide cis-regulatory effect suggests a mechanism involving long-range chromatin interactions. Here, we further dissect the disease-linked haplotype and identify putative causal DNA variants via a combination of genetic and functional analyses. First, high-throughput resequencing of the region and genotyping of potential candidate variants were performed. Next, additional mapping of allelic expression differences in Yoruba HapMap LCLs allowed us to fine-map the basis of the cis-regulatory differences to a handful of candidate functional variants. Functional assays identified allele-specific differences in nucleosome distribution, an allele-specific association with the insulator protein CTCF, as well as a weak promoter activity for rs12936231. Overall, this study shows a common disease allele linked to changes in CTCF binding and nucleosome occupancy leading to altered domain-wide cis-regulation. Finally, a strong association between asthma and cis-regulatory haplotypes was observed in three independent family-based cohorts (p = 1.78 x 10(-8)). This study demonstrates the requirement of multiple parallel allele-specific tools for the investigation of noncoding disease variants and functional fine-mapping of human disease-associated haplotypes.

Functional impact of sequence variation in the promoter region of TGFB1.

Pathological deregulation of the transforming growth factor, beta 1 (TGFB1) pathway has been implicated in the development of several major diseases, including cancers. Regulatory variation in the TGFB1 gene may lead to altered TGFB1 expression and activity, and thus, modulate an individual's susceptibility to disease. Here, we performed a study of the functional relevance of cis-acting regulatory variation in the proximal promoter region of the TGFB1 gene. In a previous study, 9 promoter polymorphisms were identified in the 2kb region upstream of the transcription start site and 9 distinct promoter haplotypes were inferred from a panel of individuals from 5 distinct continental population groups. Following experimental validation, we found that the 2 major haplotypes significantly influenced TGFB1 transcriptional activity in an allele-specific manner and that 3 of the SNPs (-1886G>A, -509C>T and -1550DEL/AGG) altered DNA-protein complex formation. Though the biological relevance of these findings remains to be verified, our study suggests that polymorphisms in the TGFB1 promoter could indeed influence gene expression and potentially contribute to the pathogenesis of TGFB1 related diseases.

In vivo footprinting analysis of the Glypican 3 (GPC3) promoter region in neuroblastoma cells.

Glypican 3 (GPC3) is an X-linked gene that has its peak expression during development and is down-regulated in all studied tissues after birth. We have shown that GPC3 was expressed in neuroblastoma and Wilms' tumor. To understand the mechanisms regulating the transcription of this gene in neuroblastoma cells, we have focused our study on the identification of putative transcription factors binding the promoter. In this report we performed in vivo dimethylsulfate, UV type C irradiation and DNaseI footprinting analyses coupled with ligation-mediated PCR on nearly 1000 bp of promoter in two neuroblastoma cell lines, SJNB-7 (expressing GPC3) and SK-N-FI (not expressing GPC3). Nucleosome signature footprints were observed in the most distal part of the studied region in both cell lines. We detected eight large differentially protected regions, suggesting the presence of binding proteins in both cell lines but more DNA-protein interactions in GPC3-expressing cells. Sp1 was previously shown to be able to bind some of these regions. Here by combining electromobility shift assays and chromatin immunoprecipitations we showed that the transcription factor NFY was part of the DNA-protein complex found in footprinted regions upstream of the described minimal promoter. These studies performed on chromatin in situ suggest that NFY and yet unknown cell type-specific factors may play an important role in the regulation of GPC3.

Promoter SNPs in G1/S checkpoint regulators and their impact on the susceptibility to childhood leukemia.

Mutations leading to the alteration of cell-cycle checkpoint functions are a common feature of most cancers. Because of the highly regulated nature of the cell cycle, it seems likely that variation in gene dosage of key components due to functional regulatory polymorphisms could play an important role in cancer development. Here we provide evidence of the involvement of promoter single-nucleotide polymorphisms (pSNPs) in the cyclin-dependent-kinase inhibitor genes CDKN2A, CDKN2B, CDKN1A, and CDKN1B in the etiology of childhood pre-B acute lymphoblastic leukemia (ALL). A case-control study, conducted in 240 patients with pre-B ALL and 277 healthy controls, combined with a family-based analysis using 135 parental trios, all of French-Canadian origin, were used to evaluate single-site genotypic as well as multilocus haplotypic associations for a total of 10 pSNPs. Using both study designs, we showed evidence of association between variants CDKN2A -222A, CDKN2B -593A, and CDKN1B -1608A, and an increased risk of ALL. These findings suggest that variable expression levels of cell-cycle inhibitor genes CDKN2A, CDKN2B, and CDKN1B due to regulatory polymorphisms could indeed influence the risk of childhood pre-B ALL and contribute to carcinogenesis.

Relationship between cholesteryl ester transfer protein and LDL heterogeneity in familial hypercholesterolemia.

Small, dense LDL particles have been associated with an increased risk of coronary artery disease, and cholesteryl ester transfer protein (CETP) has been suggested to play a role in LDL particle remodeling. We examined the relationship between LDL heterogeneity and plasma CETP mass in familial hypercholesterolemia (FH). LDL particles were characterized by polyacrylamide gradient gel electrophoresis in a total of 259 FH heterozygotes and 208 nonFH controls. CETP mass was measured by enzyme-linked immunosorbent assay in a subgroup of 240 participants, which included 120 FH patients matched with 120 controls. As compared with controls, FH subjects had an 11% higher CETP mass. Moreover, LDL-peak particle diameter (LDL-PPD) was significantly smaller in FH heterozygotes than in controls (258.1 +/- 4.8 vs. 259.2 +/- 4.1 A; P = 0.01) after adjustment for covariates. There was also an inverse relationship between LDL-PPD and CETP mass (R = -0.15; P = 0.02), and this relationship was abolished by adjustment for the FH/control status, indicating that LDL-PPD changes in FH are mediated, at least in part, by an increase in plasma CETP mass concentrations. These results suggest that increased plasma CETP mass concentrations could lead to significant LDL particle remodeling in FH heterozygotes and could contribute to the pathogenesis of atherosclerosis.

Effects of atorvastatin on electrophoretic characteristics of LDL particles among subjects with heterozygous familial hypercholesterolemia.

The effects of the HMG CoA reductase inhibitor atorvastatin on electrophoretic characteristics of LDL particles were evaluated in 46 patients (28 males and 18 females) with heterozygous familial hypercholesterolemia (FH) aged 20-61 carrying either a negative or a defective LDL receptor gene mutation. Following a 6 week drug-free baseline period, FH heterozygotes were treated with atorvastatin (median dose: 20 mg/day, range 10-80 mg/day)) for 6 months to maintain their plasma LDL-cholesterol concentrations between 4.0 and 5.0 mmol/l. Atorvastatin treatment significantly reduced plasma total cholesterol, LDL-cholesterol and triglyceride levels and increased plasma HDL-cholesterol. Furthermore, atorvastatin treatment significantly increased LDL peak particle diameter (LDL-PPD) by 0.5% (from 255.0+/-6.2 to 256.4+/-5.5 A, P=0.004) and reduced the absolute concentration of cholesterol among small (<255 A) and large (>260 A) LDL particles by 35% (P<0.001). Changes in LDL-PPD and plasma triglyceride levels were inversely correlated (R=-0.34; P=0.02). Stepwise multiple linear regression analyses showed that 41.6% of the variation in the LDL-PPD response to atorvastatin was attributable to the initial LDL-PPD (14.4%, P=0.003), the apo E polymorphism (12.4%, P=0.02), the nature of the LDL receptor gene mutation (9.6%, P=0.01) and change in triglyceride levels (5.2%, P=0.04). Moreover, the reduction in the cholesterol content of LDL <255 A was directly correlated with the daily dosage of atorvastatin (P=0.05). Results of the present study showed that atorvastatin alters significantly LDL heterogeneity in patients at high risk of coronary heart disease (CHD) such as FH heterozygotes. These results also suggest that genetic and metabolic factors may be important determinants of atorvastatin-induced changes of LDL particle size and distribution among FH heterozygotes.