In vivo modeling of metastatic human high-grade serous ovarian cancer in mice.

Metastasis is responsible for 90% of human cancer mortality, yet it remains a challenge to model human cancer metastasis in vivo. Here we describe mouse models of high-grade serous ovarian cancer, also known as high-grade serous carcinoma (HGSC), the most common and deadliest human ovarian cancer type. Mice genetically engineered to harbor Dicer1 and Pten inactivation and mutant p53 robustly replicate the peritoneal metastases of human HGSC with complete penetrance. Arising from the fallopian tube, tumors spread to the ovary and metastasize throughout the pelvic and peritoneal cavities, invariably inducing hemorrhagic ascites. Widespread and abundant peritoneal metastases ultimately cause mouse deaths (100%). Besides the phenotypic and histopathological similarities, mouse HGSCs also display marked chromosomal instability, impaired DNA repair, and chemosensitivity. Faithfully recapitulating the clinical metastases as well as molecular and genomic features of human HGSC, this murine model will be valuable for elucidating the mechanisms underlying the development and progression of metastatic ovarian cancer and also for evaluating potential therapies.

CNVs cause autosomal recessive genetic diseases with or without involvement of SNV/indels.

PURPOSE: Improved resolution of molecular diagnostic technologies enabled detection of smaller sized exonic level copy-number variants (CNVs). The contribution of CNVs to autosomal recessive (AR) conditions may be better recognized using a large clinical cohort. METHODS: We retrospectively investigated the CNVs' contribution to AR conditions in cases subjected to chromosomal microarray analysis (CMA, N = ~70,000) and/or clinical exome sequencing (ES, N = ~12,000) at Baylor Genetics; most had pediatric onset neurodevelopmental disorders. RESULTS: CNVs contributed to biallelic variations in 87 cases, including 81 singletons and three affected sibling pairs. Seventy cases had CNVs affecting both alleles, and 17 had a CNV and a single-nucleotide variant (SNV)/indel in trans. In total, 94.3% of AR-CNVs affected one gene; among these 41.4% were single-exon and 35.0% were multiexon partial-gene events. Sixty-nine percent of homozygous AR-CNVs were embedded in homozygous genomic intervals. Five cases had large deletions unmasking an SNV/indel on the intact allele for a recessive condition, resulting in multiple molecular diagnoses. CONCLUSIONS: AR-CNVs are often smaller in size, transmitted through generations, and underrecognized due to limitations in clinical CNV detection methods. Our findings from a large clinical cohort emphasized integrated CNV and SNV/indel analyses for precise clinical and molecular diagnosis especially in the context of genomic disorders.

Recurrent De Novo and Biallelic Variation of ATAD3A, Encoding a Mitochondrial Membrane Protein, Results in Distinct Neurological Syndromes.

ATPase family AAA-domain containing protein 3A (ATAD3A) is a nuclear-encoded mitochondrial membrane protein implicated in mitochondrial dynamics, nucleoid organization, protein translation, cell growth, and cholesterol metabolism. We identified a recurrent de novo ATAD3A c.1582C>T (p.Arg528Trp) variant by whole-exome sequencing (WES) in five unrelated individuals with a core phenotype of global developmental delay, hypotonia, optic atrophy, axonal neuropathy, and hypertrophic cardiomyopathy. We also describe two families with biallelic variants in ATAD3A, including a homozygous variant in two siblings, and biallelic ATAD3A deletions mediated by nonallelic homologous recombination (NAHR) between ATAD3A and gene family members ATAD3B and ATAD3C. Tissue-specific overexpression of borR534W, the Drosophila mutation homologous to the human c.1582C>T (p.Arg528Trp) variant, resulted in a dramatic decrease in mitochondrial content, aberrant mitochondrial morphology, and increased autophagy. Homozygous null bor larvae showed a significant decrease of mitochondria, while overexpression of borWT resulted in larger, elongated mitochondria. Finally, fibroblasts of an affected individual exhibited increased mitophagy. We conclude that the p.Arg528Trp variant functions through a dominant-negative mechanism that results in small mitochondria that trigger mitophagy, resulting in a reduction in mitochondrial content. ATAD3A variation represents an additional link between mitochondrial dynamics and recognizable neurological syndromes, as seen with MFN2, OPA1, DNM1L, and STAT2 mutations.

Recurrent Muscle Weakness with Rhabdomyolysis, Metabolic Crises, and Cardiac Arrhythmia Due to Bi-allelic TANGO2 Mutations.

The underlying genetic etiology of rhabdomyolysis remains elusive in a significant fraction of individuals presenting with recurrent metabolic crises and muscle weakness. Using exome sequencing, we identified bi-allelic mutations in TANGO2 encoding transport and Golgi organization 2 homolog (Drosophila) in 12 subjects with episodic rhabdomyolysis, hypoglycemia, hyperammonemia, and susceptibility to life-threatening cardiac tachyarrhythmias. A recurrent homozygous c.460G>A (p.Gly154Arg) mutation was found in four unrelated individuals of Hispanic/Latino origin, and a homozygous ∼34 kb deletion affecting exons 3-9 was observed in two families of European ancestry. One individual of mixed Hispanic/European descent was found to be compound heterozygous for c.460G>A (p.Gly154Arg) and the deletion of exons 3-9. Additionally, a homozygous exons 4-6 deletion was identified in a consanguineous Middle Eastern Arab family. No homozygotes have been reported for these changes in control databases. Fibroblasts derived from a subject with the recurrent c.460G>A (p.Gly154Arg) mutation showed evidence of increased endoplasmic reticulum stress and a reduction in Golgi volume density in comparison to control. Our results show that the c.460G>A (p.Gly154Arg) mutation and the exons 3-9 heterozygous deletion in TANGO2 are recurrent pathogenic alleles present in the Latino/Hispanic and European populations, respectively, causing considerable morbidity in the homozygotes in these populations.

Atlas-CNV: a validated approach to call single-exon CNVs in the eMERGESeq gene panel.

PURPOSE: To provide a validated method to confidently identify exon-containing copy-number variants (CNVs), with a low false discovery rate (FDR), in targeted sequencing data from a clinical laboratory with particular focus on single-exon CNVs. METHODS: DNA sequence coverage data are normalized within each sample and subsequently exonic CNVs are identified in a batch of samples, when the target log2 ratio of the sample to the batch median exceeds defined thresholds. The quality of exonic CNV calls is assessed by C-scores (Z-like scores) using thresholds derived from gold standard samples and simulation studies. We integrate an ExonQC threshold to lower FDR and compare performance with alternate software (VisCap). RESULTS: Thirteen CNVs were used as a truth set to validate Atlas-CNV and compared with VisCap. We demonstrated FDR reduction in validation, simulation, and 10,926 eMERGESeq samples without sensitivity loss. Sixty-four multiexon and 29 single-exon CNVs with high C-scores were assessed by Multiplex Ligation-dependent Probe Amplification (MLPA). CONCLUSION: Atlas-CNV is validated as a method to identify exonic CNVs in targeted sequencing data generated in the clinical laboratory. The ExonQC and C-score assignment can reduce FDR (identification of targets with high variance) and improve calling accuracy of single-exon CNVs respectively. We propose guidelines and criteria to identify high confidence single-exon CNVs.

Homozygous and hemizygous CNV detection from exome sequencing data in a Mendelian disease cohort.

We developed an algorithm, HMZDelFinder, that uses whole exome sequencing (WES) data to identify rare and intragenic homozygous and hemizygous (HMZ) deletions that may represent complete loss-of-function of the indicated gene. HMZDelFinder was applied to 4866 samples in the Baylor-Hopkins Center for Mendelian Genomics (BHCMG) cohort and detected 773 HMZ deletion calls (567 homozygous or 206 hemizygous) with an estimated sensitivity of 86.5% (82% for single-exonic and 88% for multi-exonic calls) and precision of 78% (53% single-exonic and 96% for multi-exonic calls). Out of 773 HMZDelFinder-detected deletion calls, 82 were subjected to array comparative genomic hybridization (aCGH) and/or breakpoint PCR and 64 were confirmed. These include 18 single-exon deletions out of which 8 were exclusively detected by HMZDelFinder and not by any of seven other CNV detection tools examined. Further investigation of the 64 validated deletion calls revealed at least 15 pathogenic HMZ deletions. Of those, 7 accounted for 17-50% of pathogenic CNVs in different disease cohorts where 7.1-11% of the molecular diagnosis solved rate was attributed to CNVs. In summary, we present an algorithm to detect rare, intragenic, single-exon deletion CNVs using WES data; this tool can be useful for disease gene discovery efforts and clinical WES analyses.

CRIPT exonic deletion and a novel missense mutation in a female with short stature, dysmorphic features, microcephaly, and pigmentary abnormalities.

Mutations in CRIPT encoding cysteine-rich PDZ domain-binding protein are rare, and to date have been reported in only two patients with autosomal recessive primordial dwarfism and distinctive facies. Here, we describe a female with biallelic mutations in CRIPT presenting with postnatal growth retardation, global developmental delay, and dysmorphic features including frontal bossing, high forehead, and sparse hair and eyebrows. Additional clinical features included high myopia, admixed hyper- and hypopigmented macules primarily on the face, arms, and legs, and syndactyly of 4-5 toes bilaterally. Using whole exome sequencing (WES) and chromosomal microarray analysis (CMA), we detected a c.8G>A (p.C3Y) missense variant in exon 1 of the CRIPT gene inherited from the mother and a 1,331 bp deletion encompassing exon 1, inherited from the father. The c.8G>A (p.C3Y) missense variant in CRIPT was apparently homozygous in the proband due to the exon 1 deletion. Our findings illustrate the clinical utility of combining WES with copy number variant (CNV) analysis to provide a molecular diagnosis to patients with rare Mendelian disorders. Our findings also illustrate the clinical spectrum of CRIPT related mutations. © 2016 Wiley Periodicals, Inc.

Genome-wide association scan identifies a colorectal cancer susceptibility locus on chromosome 8q24.

Using a multistage genetic association approach comprising 7,480 affected individuals and 7,779 controls, we identified markers in chromosomal region 8q24 associated with colorectal cancer. In stage 1, we genotyped 99,632 SNPs in 1,257 affected individuals and 1,336 controls from Ontario. In stages 2-4, we performed serial replication studies using 4,024 affected individuals and 4,042 controls from Seattle, Newfoundland and Scotland. We identified one locus on chromosome 8q24 and another on 9p24 having combined odds ratios (OR) for stages 1-4 of 1.18 (trend; P = 1.41 x 10(-8)) and 1.14 (trend; P = 1.32 x 10(-5)), respectively. Additional analyses in 2,199 affected individuals and 2,401 controls from France and Europe supported the association at the 8q24 locus (OR = 1.16, trend; 95% confidence interval (c.i.): 1.07-1.26; P = 5.05 x 10(-4)). A summary across all seven studies at the 8q24 locus was highly significant (OR = 1.17, c.i.: 1.12-1.23; P = 3.16 x 10(-11)). This locus has also been implicated in prostate cancer.

Association analysis using next-generation sequence data from publicly available control groups: the robust variance score statistic.

MOTIVATION: Sufficiently powered case-control studies with next-generation sequence (NGS) data remain prohibitively expensive for many investigators. If feasible, a more efficient strategy would be to include publicly available sequenced controls. However, these studies can be confounded by differences in sequencing platform; alignment, single nucleotide polymorphism and variant calling algorithms; read depth; and selection thresholds. Assuming one can match cases and controls on the basis of ethnicity and other potential confounding factors, and one has access to the aligned reads in both groups, we investigate the effect of systematic differences in read depth and selection threshold when comparing allele frequencies between cases and controls. We propose a novel likelihood-based method, the robust variance score (RVS), that substitutes genotype calls by their expected values given observed sequence data. RESULTS: We show theoretically that the RVS eliminates read depth bias in the estimation of minor allele frequency. We also demonstrate that, using simulated and real NGS data, the RVS method controls Type I error and has comparable power to the 'gold standard' analysis with the true underlying genotypes for both common and rare variants. AVAILABILITY AND IMPLEMENTATION: An RVS R script and instructions can be found at strug.research.sickkids.ca, and at https://github.com/strug-lab/RVS. CONTACT: lisa.strug@utoronto.ca SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Unraveling the complex genetic model for cystic fibrosis: pleiotropic effects of modifier genes on early cystic fibrosis-related morbidities.

The existence of pleiotropy in disorders with multi-organ involvement can suggest therapeutic targets that could ameliorate overall disease severity. Here we assessed pleiotropy of modifier genes in cystic fibrosis (CF). CF, caused by mutations in the cystic fibrosis transmembrane conductance regulator (CFTR) gene, affects the lungs, liver, pancreas and intestines. However, modifier genes contribute to variable disease severity across affected organs, even in individuals with the same CFTR genotype. We sought to determine whether SLC26A9, SLC9A3 and SLC6A14, that contribute to meconium ileus in CF, are pleiotropic for other early-affecting CF co-morbidities. In the Canadian CF population, we assessed evidence for pleiotropic effects on (1) pediatric lung disease severity (n = 815), (2) age at first acquisition of Pseudomonas aeruginosa (P. aeruginosa) (n = 730), and (3) prenatal pancreatic damage measured by immunoreactive trypsinogen (n = 126). A multiple-phenotype analytic strategy assessed evidence for pleiotropy in the presence of phenotypic correlation. We required the same alleles to be associated with detrimental effects. SLC26A9 was pleiotropic for meconium ileus and pancreatic damage (p = 0.002 at rs7512462), SLC9A3 for meconium ileus and lung disease (p = 1.5 × 10(-6) at rs17563161), and SLC6A14 for meconium ileus and both lung disease and age at first P. aeruginosa infection (p = 0.0002 and p = 0.006 at rs3788766, respectively). The meconium ileus risk alleles in SLC26A9, SLC9A3 and SLC6A14 are pleiotropic, increasing risk for other early CF co-morbidities. Furthermore, co-morbidities affecting the same organ tended to associate with the same genes. The existence of pleiotropy within this single disorder suggests that complementary therapeutic strategies to augment solute transport will benefit multiple CF-associated tissues.

Molecular findings among patients referred for clinical whole-exome sequencing.

IMPORTANCE: Clinical whole-exome sequencing is increasingly used for diagnostic evaluation of patients with suspected genetic disorders. OBJECTIVE: To perform clinical whole-exome sequencing and report (1) the rate of molecular diagnosis among phenotypic groups, (2) the spectrum of genetic alterations contributing to disease, and (3) the prevalence of medically actionable incidental findings such as FBN1 mutations causing Marfan syndrome. DESIGN, SETTING, AND PATIENTS: Observational study of 2000 consecutive patients with clinical whole-exome sequencing analyzed between June 2012 and August 2014. Whole-exome sequencing tests were performed at a clinical genetics laboratory in the United States. Results were reported by clinical molecular geneticists certified by the American Board of Medical Genetics and Genomics. Tests were ordered by the patient's physician. The patients were primarily pediatric (1756 [88%]; mean age, 6 years; 888 females [44%], 1101 males [55%], and 11 fetuses [1% gender unknown]), demonstrating diverse clinical manifestations most often including nervous system dysfunction such as developmental delay. MAIN OUTCOMES AND MEASURES: Whole-exome sequencing diagnosis rate overall and by phenotypic category, mode of inheritance, spectrum of genetic events, and reporting of incidental findings. RESULTS: A molecular diagnosis was reported for 504 patients (25.2%) with 58% of the diagnostic mutations not previously reported. Molecular diagnosis rates for each phenotypic category were 143/526 (27.2%; 95% CI, 23.5%-31.2%) for the neurological group, 282/1147 (24.6%; 95% CI, 22.1%-27.2%) for the neurological plus other organ systems group, 30/83 (36.1%; 95% CI, 26.1%-47.5%) for the specific neurological group, and 49/244 (20.1%; 95% CI, 15.6%-25.8%) for the nonneurological group. The Mendelian disease patterns of the 527 molecular diagnoses included 280 (53.1%) autosomal dominant, 181 (34.3%) autosomal recessive (including 5 with uniparental disomy), 65 (12.3%) X-linked, and 1 (0.2%) mitochondrial. Of 504 patients with a molecular diagnosis, 23 (4.6%) had blended phenotypes resulting from 2 single gene defects. About 30% of the positive cases harbored mutations in disease genes reported since 2011. There were 95 medically actionable incidental findings in genes unrelated to the phenotype but with immediate implications for management in 92 patients (4.6%), including 59 patients (3%) with mutations in genes recommended for reporting by the American College of Medical Genetics and Genomics. CONCLUSIONS AND RELEVANCE: Whole-exome sequencing provided a potential molecular diagnosis for 25% of a large cohort of patients referred for evaluation of suspected genetic conditions, including detection of rare genetic events and new mutations contributing to disease. The yield of whole-exome sequencing may offer advantages over traditional molecular diagnostic approaches in certain patients.

Copy number variant and runs of homozygosity detection by microarrays enabled more precise molecular diagnoses in 11,020 clinical exome cases.

BACKGROUND: Exome sequencing (ES) has been successfully applied in clinical detection of single nucleotide variants (SNVs) and small indels. However, identification of copy number variants (CNVs) using ES data remains challenging. The purpose of this study is to understand the contribution of CNVs and copy neutral runs of homozygosity (ROH) in molecular diagnosis of patients referred for ES. METHODS: In a cohort of 11,020 consecutive ES patients, an Illumina SNP array analysis interrogating mostly coding SNPs was performed as a quality control (QC) measurement and for CNV/ROH detection. Among these patients, clinical chromosomal microarray analysis (CMA) was performed at Baylor Genetics (BG) on 3229 patients, either before, concurrently, or after ES. We retrospectively analyzed the findings from CMA and the QC array. RESULTS: The QC array can detect ~ 70% of pathogenic/likely pathogenic CNVs (PCNVs) detectable by CMA. Out of the 11,020 ES cases, the QC array identified PCNVs in 327 patients and uniparental disomy (UPD) disorder-related ROH in 10 patients. The overall PCNV/UPD detection rate was 5.9% in the 3229 ES patients who also had CMA at BG; PCNV/UPD detection rate was higher in concurrent ES and CMA than in ES with prior CMA (7.2% vs 4.6%). The PCNVs/UPD contributed to the molecular diagnoses in 17.4% (189/1089) of molecularly diagnosed ES cases with CMA and were estimated to contribute in 10.6% of all molecularly diagnosed ES cases. Dual diagnoses with both PCNVs and SNVs were detected in 38 patients. PCNVs affecting single recessive disorder genes in a compound heterozygous state with SNVs were detected in 4 patients, and homozygous deletions (mostly exonic deletions) were detected in 17 patients. A higher PCNV detection rate was observed for patients with syndromic phenotypes and/or cardiovascular abnormalities. CONCLUSIONS: Our clinical genomics study demonstrates that detection of PCNV/UPD through the QC array or CMA increases ES diagnostic rate, provides more precise molecular diagnosis for dominant as well as recessive traits, and enables more complete genetic diagnoses in patients with dual or multiple molecular diagnoses. Concurrent ES and CMA using an array with exonic coverage for disease genes enables most effective detection of both CNVs and SNVs and therefore is recommended especially in time-sensitive clinical situations.

Lessons learned from additional research analyses of unsolved clinical exome cases.

BACKGROUND: Given the rarity of most single-gene Mendelian disorders, concerted efforts of data exchange between clinical and scientific communities are critical to optimize molecular diagnosis and novel disease gene discovery. METHODS: We designed and implemented protocols for the study of cases for which a plausible molecular diagnosis was not achieved in a clinical genomics diagnostic laboratory (i.e. unsolved clinical exomes). Such cases were recruited to a research laboratory for further analyses, in order to potentially: (1) accelerate novel disease gene discovery; (2) increase the molecular diagnostic yield of whole exome sequencing (WES); and (3) gain insight into the genetic mechanisms of disease. Pilot project data included 74 families, consisting mostly of parent-offspring trios. Analyses performed on a research basis employed both WES from additional family members and complementary bioinformatics approaches and protocols. RESULTS: Analysis of all possible modes of Mendelian inheritance, focusing on both single nucleotide variants (SNV) and copy number variant (CNV) alleles, yielded a likely contributory variant in 36% (27/74) of cases. If one includes candidate genes with variants identified within a single family, a potential contributory variant was identified in a total of ~51% (38/74) of cases enrolled in this pilot study. The molecular diagnosis was achieved in 30/63 trios (47.6%). Besides this, the analysis workflow yielded evidence for pathogenic variants in disease-associated genes in 4/6 singleton cases (66.6%), 1/1 multiplex family involving three affected siblings, and 3/4 (75%) quartet families. Both the analytical pipeline and the collaborative efforts between the diagnostic and research laboratories provided insights that allowed recent disease gene discoveries (PURA, TANGO2, EMC1, GNB5, ATAD3A, and MIPEP) and increased the number of novel genes, defined in this study as genes identified in more than one family (DHX30 and EBF3). CONCLUSION: An efficient genomics pipeline in which clinical sequencing in a diagnostic laboratory is followed by the detailed reanalysis of unsolved cases in a research environment, supplemented with WES data from additional family members, and subject to adjuvant bioinformatics analyses including relaxed variant filtering parameters in informatics pipelines, can enhance the molecular diagnostic yield and provide mechanistic insights into Mendelian disorders. Implementing these approaches requires collaborative clinical molecular diagnostic and research efforts.

Use of Exome Sequencing for Infants in Intensive Care Units: Ascertainment of Severe Single-Gene Disorders and Effect on Medical Management.

Importance: While congenital malformations and genetic diseases are a leading cause of early infant death, to our knowledge, the contribution of single-gene disorders in this group is undetermined. Objective: To determine the diagnostic yield and use of clinical exome sequencing in critically ill infants. Design, Setting, and Participants: Clinical exome sequencing was performed for 278 unrelated infants within the first 100 days of life who were admitted to Texas Children's Hospital in Houston, Texas, during a 5-year period between December 2011 and January 2017. Exome sequencing types included proband exome, trio exome, and critical trio exome, a rapid genomic assay for seriously ill infants. Main Outcomes and Measures: Indications for testing, diagnostic yield of clinical exome sequencing, turnaround time, molecular findings, patient age at diagnosis, and effect on medical management among a group of critically ill infants who were suspected to have genetic disorders. Results: The mean (SEM) age for infants participating in the study was 28.5 (1.7) days; of these, the mean (SEM) age was 29.0 (2.2) days for infants undergoing proband exome sequencing, 31.5 (3.9) days for trio exome, and 22.7 (3.9) days for critical trio exome. Clinical indications for exome sequencing included a range of medical concerns. Overall, a molecular diagnosis was achieved in 102 infants (36.7%) by clinical exome sequencing, with relatively low yield for cardiovascular abnormalities. The diagnosis affected medical management for 53 infants (52.0%) and had a substantial effect on informed redirection of care, initiation of new subspecialist care, medication/dietary modifications, and furthering life-saving procedures in select patients. Critical trio exome sequencing revealed a molecular diagnosis in 32 of 63 infants (50.8%) at a mean (SEM) of 33.1 (5.6) days of life with a mean (SEM) turnaround time of 13.0 (0.4) days. Clinical care was altered by the diagnosis in 23 of 32 patients (71.9%). The diagnostic yield, patient age at diagnosis, and medical effect in the group that underwent critical trio exome sequencing were significantly different compared with the group who underwent regular exome testing. For deceased infants (n = 81), genetic disorders were molecularly diagnosed in 39 (48.1%) by exome sequencing, with implications for recurrence risk counseling. Conclusions and Relevance: Exome sequencing is a powerful tool for the diagnostic evaluation of critically ill infants with suspected monogenic disorders in the neonatal and pediatric intensive care units and its use has a notable effect on clinical decision making.

The genetics of reading disability in an often excluded sample: novel loci suggested for reading disability in Rolandic epilepsy.

BACKGROUND: Reading disability (RD) is a common neurodevelopmental disorder with genetic basis established in families segregating "pure" dyslexia. RD commonly occurs in neurodevelopmental disorders including Rolandic Epilepsy (RE), a complex genetic disorder. We performed genomewide linkage analysis of RD in RE families, testing the hypotheses that RD in RE families is genetically heterogenenous to pure dyslexia, and shares genetic influences with other sub-phenotypes of RE. METHODS: We initially performed genome-wide linkage analysis using 1000 STR markers in 38 US families ascertained through a RE proband; most of these families were multiplex for RD. We analyzed the data by two-point and multipoint parametric LOD score methods. We then confirmed the linkage evidence in a second US dataset of 20 RE families. We also resequenced the SEMA3C gene at the 7q21 linkage locus in members of one multiplex RE/RD pedigree and the DISC1 gene in affected pedigrees at the 1q42 locus. RESULTS: In the discovery dataset there was suggestive evidence of linkage for RD to chromosome 7q21 (two-point LOD score 3.05, multipoint LOD 3.08) and at 1q42 (two-point LOD 2.87, multipoint LOD 3.03). Much of the linkage evidence at 7q21 derived from families of French-Canadian origin, whereas the linkage evidence at 1q42 was well distributed across all the families. There was little evidence for linkage at known dyslexia loci. Combining the discovery and confirmation datasets increased the evidence at 1q42 (two-point LOD = 3.49, multipoint HLOD = 4.70), but decreased evidence at 7q21 (two-point LOD = 2.28, multipoint HLOD  = 1.81), possibly because the replication sample did not have French Canadian representation. DISCUSSION: Reading disability in rolandic epilepsy has a genetic basis and may be influenced by loci at 1q42 and, in some populations, at 7q21; there is little evidence of a role for known DYX loci discovered in "pure" dyslexia pedigrees. 1q42 and 7q21 are candidate novel dyslexia loci.

Multiple apical plasma membrane constituents are associated with susceptibility to meconium ileus in individuals with cystic fibrosis.

Variants associated with meconium ileus in cystic fibrosis were identified in 3,763 affected individuals by genome-wide association study (GWAS). Five SNPs at two loci near SLC6A14 at Xq23-24 (minimum P = 1.28 × 10(-12) at rs3788766) and SLC26A9 at 1q32.1 (minimum P = 9.88 × 10(-9) at rs4077468) accounted for ~5% of phenotypic variability and were replicated in an independent sample of affected individuals (n = 2,372; P = 0.001 and 0.0001, respectively). By incorporating the knowledge that disease-causing mutations in CFTR alter electrolyte and fluid flux across surface epithelium into a hypothesis-driven GWAS (GWAS-HD), we identified associations with the same SNPs in SLC6A14 and SLC26A9 and established evidence for the involvement of SNPs in a third solute carrier gene, SLC9A3. In addition, GWAS-HD provided evidence of association between meconium ileus and multiple genes encoding constituents of the apical plasma membrane where CFTR resides (P = 0.0002; testing of 155 apical membrane genes jointly and in replication, P = 0.022). These findings suggest that modulating activities of apical membrane constituents could complement current therapeutic paradigms for cystic fibrosis.

Assessment of methylation level prediction accuracy in methyl-DNA immunoprecipitation and sodium bisulfite based microarray platforms.

In this study, we verified the accuracy of two array methods--methylated DNA immunoprecipitation coupled with CpG island microarrays (MeDIP-CGI-arrays) and sodium bisulfite conversion based microarrays (BC-arrays)--in predicting regional methylation levels as measured by pyrosequencing of bisulfite converted DNA (BC-pyrosequencing). To test the accuracy of these methods we used the Agilent Human CpG island and the Illumina HumanMethylation27 microarrays respectively, and compared microarray outputs to the data from targeted BC-pyrosequencing assays from several genomic regions of corresponding samples. We observed relatively high correlation with BC-pyrosequencing data for both array platforms, R = 0.87 for BC-Array and R = 0.79 for MeDIP-CGI array. However, MeDIP-CGI array were less reliable in predicting intermediate levels of DNA methylation. Several bioinformatics strategies, to ameliorate the performance of the MeDIP-CGI-Arrays did not improve the correlation with BC-pyrosequencing data. The high scalability, low cost and simpler analysis of BC-arrays, together with the recent extended coverage may make them a more versatile methylation analysis tool.

A pure likelihood approach to the analysis of genetic association data: an alternative to Bayesian and frequentist analysis.

Investigators performing genetic association studies grapple with how to measure strength of association evidence, choose sample size, and adjust for multiple testing. We apply the evidential paradigm (EP) to genetic association studies, highlighting its strengths. The EP uses likelihood ratios (LRs), as opposed to P-values or Bayes' factors, to measure strength of association evidence. We derive EP methodology to estimate sample size, adjust for multiple testing, and provide informative graphics for drawing inferences, as illustrated with a Rolandic Epilepsy (RE) fine-mapping study. We focus on controlling the probability of observing weak evidence for or against association (W) rather than type I errors (M). For example, for LR> or =32 representing strong evidence, at one locus with n=200 cases, n=200 controls, W=0.134, whereas M=0.005. For n=300 cases and controls, W=0.039 and M=0.004. These calculations are based on detecting an OR=1.5. Despite the common misconception, one is not tied to this planning value for analysis; rather one calculates the likelihood at all possible values to assess evidence for association. We provide methodology to adjust for multiple tests across m loci, which adjusts M and W for m. We do so for (a) single-stage designs, (b) two-stage designs, and (c) simultaneously controlling family-wise error rate (FWER) and W. Method (c) chooses larger sample sizes than (a) or (b), whereas (b) has smaller bounds on the FWER than (a). The EP, using our innovative graphical display, identifies important SNPs in elongator protein complex 4 (ELP4) associated with RE that may not have been identified using standard approaches.

Specific variants in the MLH1 gene region may drive DNA methylation, loss of protein expression, and MSI-H colorectal cancer.

BACKGROUND: We previously identified an association between a mismatch repair gene, MLH1, promoter SNP (rs1800734) and microsatellite unstable (MSI-H) colorectal cancers (CRCs) in two samples. The current study expanded on this finding as we explored the genetic basis of DNA methylation in this region of chromosome 3. We hypothesized that specific polymorphisms in the MLH1 gene region predispose it to DNA methylation, resulting in the loss of MLH1 gene expression, mismatch-repair function, and consequently to genome-wide microsatellite instability. METHODOLOGY/PRINCIPAL FINDINGS: We first tested our hypothesis in one sample from Ontario (901 cases, 1,097 controls) and replicated major findings in two additional samples from Newfoundland and Labrador (479 cases, 336 controls) and from Seattle (591 cases, 629 controls). Logistic regression was used to test for association between SNPs in the region of MLH1 and CRC, MSI-H CRC, MLH1 gene expression in CRC, and DNA methylation in CRC. The association between rs1800734 and MSI-H CRCs, previously reported in Ontario and Newfoundland, was replicated in the Seattle sample. Two additional SNPs, in strong linkage disequilibrium with rs1800734, showed strong associations with MLH1 promoter methylation, loss of MLH1 protein, and MSI-H CRC in all three samples. The logistic regression model of MSI-H CRC that included MLH1-promoter-methylation status and MLH1 immunohistochemistry status fit most parsimoniously in all three samples combined. When rs1800734 was added to this model, its effect was not statistically significant (P-value  = 0.72 vs. 2.3×10(-4) when the SNP was examined alone). CONCLUSIONS/SIGNIFICANCE: The observed association of rs1800734 with MSI-H CRC occurs through its effect on the MLH1 promoter methylation, MLH1 IHC deficiency, or both.