Disruptive CHD8 mutations define a subtype of autism early in development.

Autism spectrum disorder (ASD) is a heterogeneous disease in which efforts to define subtypes behaviorally have met with limited success. Hypothesizing that genetically based subtype identification may prove more productive, we resequenced the ASD-associated gene CHD8 in 3,730 children with developmental delay or ASD. We identified a total of 15 independent mutations; no truncating events were identified in 8,792 controls, including 2,289 unaffected siblings. In addition to a high likelihood of an ASD diagnosis among patients bearing CHD8 mutations, characteristics enriched in this group included macrocephaly, distinct faces, and gastrointestinal complaints. chd8 disruption in zebrafish recapitulates features of the human phenotype, including increased head size as a result of expansion of the forebrain/midbrain and impairment of gastrointestinal motility due to a reduction in postmitotic enteric neurons. Our findings indicate that CHD8 disruptions define a distinct ASD subtype and reveal unexpected comorbidities between brain development and enteric innervation.

Population screening shows risk of inherited cancer and familial hypercholesterolemia in Oregon.

The Healthy Oregon Project (HOP) is a statewide effort that aims to build a large research repository and influence the health of Oregonians through providing no-cost genetic screening to participants for a next-generation sequencing 32-gene panel comprising genes related to inherited cancers and familial hypercholesterolemia. This type of unbiased population screening can detect at-risk individuals who may otherwise be missed by conventional medical approaches. However, challenges exist for this type of high-throughput testing in an academic setting, including developing a low-cost high-efficiency test and scaling up the clinical laboratory for processing large numbers of samples. Modifications to our academic clinical laboratory including efficient test design, robotics, and a streamlined analysis approach increased our ability to test more than 1,000 samples per month for HOP using only one dedicated HOP laboratory technologist. Additionally, enrollment using a HIPAA-compliant smartphone app and sample collection using mouthwash increased efficiency and reduced cost. Here, we present our experience three years into HOP and discuss the lessons learned, including our successes, challenges, opportunities, and future directions, as well as the genetic screening results for the first 13,670 participants tested. Overall, we have identified 730 pathogenic/likely pathogenic variants in 710 participants in 24 of the 32 genes on the panel. The carrier rate for pathogenic/likely pathogenic variants in the inherited cancer genes on the panel for an unselected population was 5.0% and for familial hypercholesterolemia was 0.3%. Our laboratory experience described here may provide a useful model for population screening projects in other states.

Shared and Distinct Functional Effects of Patient-Specific Tbr1 Mutations on Cortical Development.

T-Box Brain Transcription Factor 1 (TBR1) plays essential roles in brain development, mediating neuronal migration, fate specification, and axon tract formation. While heterozygous loss-of-function and missense TBR1 mutations are associated with neurodevelopmental conditions, the effects of these heterogeneous mutations on brain development have yet to be fully explored. We characterized multiple mouse lines carrying Tbr1 mutations differing by type and exonic location, including the previously generated Tbr1 exon 2-3 knock-out (KO) line, and we analyzed male and female mice at neonatal and adult stages. The frameshift patient mutation A136PfsX80 (A136fs) caused reduced TBR1 protein in cortex similar to Tbr1 KO, while the missense patient mutation K228E caused significant TBR1 upregulation. Analysis of cortical layer formation found similar defects between KO and A136fs homozygotes in their CUX1+ and CTIP2+ layer positions, while K228E homozygosity produced layering defects distinct from these mutants. Meanwhile, the examination of cortical apoptosis found extensive cell death in KO homozygotes but limited cell death in A136fs or K228E homozygotes. Despite their discordant cortical phenotypes, these Tbr1 mutations produced several congruent phenotypes, including anterior commissure reduction in heterozygotes, which was previously observed in humans with TBR1 mutations. These results indicate that patient-specific Tbr1 mutant mice will be valuable translational models for pinpointing shared and distinct etiologies among patients with TBR1-related developmental conditions.SIGNIFICANCE STATEMENT Mutations of the TBR1 gene increase the likelihood of neurodevelopmental conditions such as intellectual disability and autism. Therefore, the study of TBR1 can offer insights into the biological mechanisms underlying these conditions, which affect millions worldwide. To improve the modeling of TBR1-related conditions over current Tbr1 knock-out mice, we created mouse lines carrying Tbr1 mutations identical to those found in human patients. Mice with one mutant Tbr1 copy show reduced amygdalar connections regardless of mutation type, suggesting a core biomarker for TBR1-related disorders. In mice with two mutant Tbr1 copies, brain phenotypes diverge by mutation type, suggesting differences in Tbr1 gene functionality in different patients. These mouse models will serve as valuable tools for understanding genotype-phenotype relationships among patients with neurodevelopmental conditions.

An Integrated Deep-Mutational-Scanning Approach Provides Clinical Insights on PTEN Genotype-Phenotype Relationships.

Germline variation in PTEN results in variable clinical presentations, including benign and malignant neoplasia and neurodevelopmental disorders. Despite decades of research, it remains unclear how the PTEN genotype is related to clinical outcomes. In this study, we combined two recent deep mutational scanning (DMS) datasets probing the effects of single amino acid variation on enzyme activity and steady-state cellular abundance with a large, well-curated clinical cohort of PTEN-variant carriers. We sought to connect variant-specific molecular phenotypes to the clinical outcomes of individuals with PTEN variants. We found that DMS data partially explain quantitative clinical traits, including head circumference and Cleveland Clinic (CC) score, which is a semiquantitative surrogate of disease burden. We built logistic regression models that use DMS and CADD scores to separate clinical PTEN variation from gnomAD control-only variation with high accuracy. By using a survival-like analysis, we identified molecular phenotype groups with differential risk of early cancer onset as well as lifetime risk of cancer. Finally, we identified classes of DMS-defined variants with significantly different risk levels for classical hamartoma-related features (odds ratio [OR] range of 4.1-102.9). In stark contrast, the risk for developing autism or developmental delay does not significantly change across variant classes (OR range of 5.4-12.4). Together, these findings highlight the potential impact of combining DMS datasets with rich clinical data and provide new insights that might guide personalized clinical decisions for PTEN-variant carriers.

A Saturation Mutagenesis Approach to Understanding PTEN Lipid Phosphatase Activity and Genotype-Phenotype Relationships.

Phosphatase and tensin homolog (PTEN) is a tumor suppressor frequently mutated in diverse cancers. Germline PTEN mutations are also associated with a range of clinical outcomes, including PTEN hamartoma tumor syndrome (PHTS) and autism spectrum disorder (ASD). To empower new insights into PTEN function and clinically relevant genotype-phenotype relationships, we systematically evaluated the effect of PTEN mutations on lipid phosphatase activity in vivo. Using a massively parallel approach that leverages an artificial humanized yeast model, we derived high-confidence estimates of functional impact for 7,244 single amino acid PTEN variants (86% of possible). We identified 2,273 mutations with reduced cellular lipid phosphatase activity, which includes 1,789 missense mutations. These data recapitulated known functional findings but also uncovered new insights into PTEN protein structure, biochemistry, and mutation tolerance. Several residues in the catalytic pocket showed surprising mutational tolerance. We identified that the solvent exposure of wild-type residues is a critical determinant of mutational tolerance. Further, we created a comprehensive functional map by leveraging correlations between amino acid substitutions to impute functional scores for all variants, including those not present in the assay. Variant functional scores can reliably discriminate likely pathogenic from benign alleles. Further, 32% of ClinVar unclassified missense variants are phosphatase deficient in our assay, supporting their reclassification. ASD-associated mutations generally had less severe fitness scores relative to PHTS-associated mutations (p = 7.16 × 10-5) and a higher fraction of hypomorphic mutations, arguing for continued genotype-phenotype studies in larger clinical datasets that can further leverage these rich functional data.

Autism and attention-deficit/hyperactivity disorders and symptoms in children with neurofibromatosis type 1.

AIM: To evaluate if autism symptoms and diagnoses are more common in children with neurofibromatosis type 1 (NF1) than in typically developing children, to which levels, and to determine if co-occurring attention-deficit/hyperactivity disorder (ADHD) symptomatology accounts for this increase. METHOD: We searched hospital electronic medical records (EMR) for International Classification of Diseases, 10th Revision NF1 and co-occurring diagnoses codes. We recruited a subsample of 45 children (mean age 9y 2mo; SD 2y 7mo; range 5-12y; 22 males, 23 females) and collected parental reports of autism symptomatology, adaptive behavior, and behavioral problems that were compared to those of 360 age- and sex-matched controls from the Simons Simplex Collection (SSC) with autism spectrum disorder (ASD; SSC-ASD) or typically developing (SSC-TD). RESULTS: The EMR search identified 968 children with NF1; 8.8% had ADHD and 2.1% had ASD co-occurring diagnoses. In the subsample, the mean autism scale score for participants with NF1 was below cut-off for significant autism symptoms. Participants with NF1 had significantly more autism and behavioral symptoms than SSC-TD participants, and significantly less than SSC-ASD participants, with one exception: ADHD symptom levels were similar to those of SSC-ASD participants. In analyses that controlled for internalizing, ADHD, and communication scores, the difference in autism symptom levels between participants with NF1 and typically developing controls disappeared almost entirely. INTERPRETATION: Our results do not support an association between NF1 and autism, both at the symptom and disorder levels. WHAT THIS PAPER ADDS: Diagnoses of attention-deficit/hyperactivity disorder (ADHD) were more common in children with neurofibromatosis type 1 (NF1) than in the general child population. Diagnoses of autism spectrum disorder were no more common in children with NF1 than in the general child population. Increases in autism symptoms did not reach clinically significant thresholds. Co-occurring ADHD symptoms accounted for increased autism questionnaire scores. Adaptive behavior in participants with NF1 showed normal socialization but lower communication proficiency.

Exonic Mosaic Mutations Contribute Risk for Autism Spectrum Disorder.

Genetic risk factors for autism spectrum disorder (ASD) have yet to be fully elucidated. Postzygotic mosaic mutations (PMMs) have been implicated in several neurodevelopmental disorders and overgrowth syndromes. By leveraging whole-exome sequencing data on a large family-based ASD cohort, the Simons Simplex Collection, we systematically evaluated the potential role of PMMs in autism risk. Initial re-evaluation of published single-nucleotide variant (SNV) de novo mutations showed evidence consistent with putative PMMs for 11% of mutations. We developed a robust and sensitive SNV PMM calling approach integrating complementary callers, logistic regression modeling, and additional heuristics. In our high-confidence call set, we identified 470 PMMs in children, increasing the proportion of mosaic SNVs to 22%. Probands have a significant burden of synonymous PMMs and these mutations are enriched for computationally predicted impacts on splicing. Evidence of increased missense PMM burden was not seen in the full cohort. However, missense burden signal increased in subcohorts of families where probands lacked nonsynonymous germline mutations, especially in genes intolerant to mutations. Parental mosaic mutations that were transmitted account for 6.8% of the presumed de novo mutations in the children. PMMs were identified in previously implicated high-confidence neurodevelopmental disorder risk genes, such as CHD2, CTNNB1, SCN2A, and SYNGAP1, as well as candidate risk genes with predicted functions in chromatin remodeling or neurodevelopment, including ACTL6B, BAZ2B, COL5A3, SSRP1, and UNC79. We estimate that PMMs potentially contribute risk to 3%-4% of simplex ASD case subjects and future studies of PMMs in ASD and related disorders are warranted.

The contribution of de novo coding mutations to autism spectrum disorder.

Whole exome sequencing has proven to be a powerful tool for understanding the genetic architecture of human disease. Here we apply it to more than 2,500 simplex families, each having a child with an autistic spectrum disorder. By comparing affected to unaffected siblings, we show that 13% of de novo missense mutations and 43% of de novo likely gene-disrupting (LGD) mutations contribute to 12% and 9% of diagnoses, respectively. Including copy number variants, coding de novo mutations contribute to about 30% of all simplex and 45% of female diagnoses. Almost all LGD mutations occur opposite wild-type alleles. LGD targets in affected females significantly overlap the targets in males of lower intelligence quotient (IQ), but neither overlaps significantly with targets in males of higher IQ. We estimate that LGD mutation in about 400 genes can contribute to the joint class of affected females and males of lower IQ, with an overlapping and similar number of genes vulnerable to contributory missense mutation. LGD targets in the joint class overlap with published targets for intellectual disability and schizophrenia, and are enriched for chromatin modifiers, FMRP-associated genes and embryonically expressed genes. Most of the significance for the latter comes from affected females.

Sporadic autism exomes reveal a highly interconnected protein network of de novo mutations.

It is well established that autism spectrum disorders (ASD) have a strong genetic component; however, for at least 70% of cases, the underlying genetic cause is unknown. Under the hypothesis that de novo mutations underlie a substantial fraction of the risk for developing ASD in families with no previous history of ASD or related phenotypes--so-called sporadic or simplex families--we sequenced all coding regions of the genome (the exome) for parent-child trios exhibiting sporadic ASD, including 189 new trios and 20 that were previously reported. Additionally, we also sequenced the exomes of 50 unaffected siblings corresponding to these new (n = 31) and previously reported trios (n = 19), for a total of 677 individual exomes from 209 families. Here we show that de novo point mutations are overwhelmingly paternal in origin (4:1 bias) and positively correlated with paternal age, consistent with the modest increased risk for children of older fathers to develop ASD. Moreover, 39% (49 of 126) of the most severe or disruptive de novo mutations map to a highly interconnected ß-catenin/chromatin remodelling protein network ranked significantly for autism candidate genes. In proband exomes, recurrent protein-altering mutations were observed in two genes: CHD8 and NTNG1. Mutation screening of six candidate genes in 1,703 ASD probands identified additional de novo, protein-altering mutations in GRIN2B, LAMC3 and SCN1A. Combined with copy number variant (CNV) data, these results indicate extreme locus heterogeneity but also provide a target for future discovery, diagnostics and therapeutics.

Current Perspectives in Autism Spectrum Disorder: From Genes to Therapy.

Autism spectrum disorder (ASD) is a constellation of neurodevelopmental presentations with high heritability and both phenotypic and genetic heterogeneity. To date, mutations in hundreds of genes have been associated to varying degrees with increased ASD risk. A better understanding of the functions of these genes and whether they fit together in functional groups or impact similar neuronal circuits is needed to develop rational treatment strategies. We will review current areas of emphasis in ASD research, starting from human genetics and exploring how mouse models of human mutations have helped identify specific molecular pathways (protein synthesis and degradation, chromatin remodeling, intracellular signaling), which are linked to alterations in circuit function and cognitive/social behavior. We will conclude by discussing how we can leverage the findings on molecular and cellular alterations found in ASD to develop therapies for neurodevelopmental disorders.

Rare-variant extensions of the transmission disequilibrium test: application to autism exome sequence data.

Many population-based rare-variant (RV) association tests, which aggregate variants across a region, have been developed to analyze sequence data. A drawback of analyzing population-based data is that it is difficult to adequately control for population substructure and admixture, and spurious associations can occur. For RVs, this problem can be substantial, because the spectrum of rare variation can differ greatly between populations. A solution is to analyze parent-child trio data, by using the transmission disequilibrium test (TDT), which is robust to population substructure and admixture. We extended the TDT to test for RV associations using four commonly used methods. We demonstrate that for all RV-TDT methods, using proper analysis strategies, type I error is well-controlled even when there are high levels of population substructure or admixture. For trio data, unlike for population-based data, RV allele-counting association methods will lead to inflated type I errors. However type I errors can be properly controlled by obtaining p values empirically through haplotype permutation. The power of the RV-TDT methods was evaluated and compared to the analysis of case-control data with a number of genetic and disease models. The RV-TDT was also used to analyze exome data from 199 Simons Simplex Collection autism trios and an association was observed with variants in ABCA7. Given the problem of adequately controlling for population substructure and admixture in RV association studies and the growing number of sequence-based trio studies, the RV-TDT is extremely beneficial to elucidate the involvement of RVs in the etiology of complex traits.

Mutations in CSPP1 cause primary cilia abnormalities and Joubert syndrome with or without Jeune asphyxiating thoracic dystrophy.

Joubert syndrome (JBTS) is a recessive ciliopathy in which a subset of affected individuals also have the skeletal dysplasia Jeune asphyxiating thoracic dystrophy (JATD). Here, we have identified biallelic truncating CSPP1 (centrosome and spindle pole associated protein 1) mutations in 19 JBTS-affected individuals, four of whom also have features of JATD. CSPP1 mutations explain ∼5% of JBTS in our cohort, and despite truncating mutations in all affected individuals, the range of phenotypic severity is broad. Morpholino knockdown of cspp1 in zebrafish caused phenotypes reported in other zebrafish models of JBTS (curved body shape, pronephric cysts, and cerebellar abnormalities) and reduced ciliary localization of Arl13b, further supporting loss of CSPP1 function as a cause of JBTS. Fibroblasts from affected individuals with CSPP1 mutations showed reduced numbers of primary cilia and/or short primary cilia, as well as reduced axonemal localization of ciliary proteins ARL13B and adenylyl cyclase III. In summary, CSPP1 mutations are a major cause of the Joubert-Jeune phenotype in humans; however, the mechanism by which these mutations lead to both JBTS and JATD remains unknown.

Large-scale targeted sequencing comparison highlights extreme genetic heterogeneity in nephronophthisis-related ciliopathies.

BACKGROUND: The term nephronophthisis-related ciliopathies (NPHP-RC) describes a group of rare autosomal-recessive cystic kidney diseases, characterised by broad genetic and clinical heterogeneity. NPHP-RC is frequently associated with extrarenal manifestations and accounts for the majority of genetically caused chronic kidney disease (CKD) during childhood and adolescence. Generation of a molecular diagnosis has been impaired by this broad genetic heterogeneity. However, recently developed high-throughput exon sequencing techniques represent powerful and efficient tools to screen large cohorts for dozens of causative genes. METHODS: Therefore, we performed massively multiplexed targeted sequencing using the modified molecular inversion probe strategy (MIPs) in an international cohort of 384 patients diagnosed with NPHP-RC. RESULTS: As a result, we established the molecular diagnoses in 81/384 unrelated individuals (21.1%). We detected 127 likely disease-causing mutations in 18 of 34 evaluated NPHP-RC genes, 22 of which were novel. We further compared a subgroup of current findings to the results of a previous study in which we used an array-based microfluidic PCR technology in the same cohort. While 78 likely disease-causing mutations were previously detected by the array-based microfluidic PCR, the MIPs approach identified 94 likely pathogenic mutations. Compared with the previous approach, MIPs redetected 66 out of 78 variants and 28 previously unidentified variants, for a total of 94 variants. CONCLUSIONS: In summary, we demonstrate that the modified MIPs technology is a useful approach to screen large cohorts for a multitude of established NPHP genes in order to identify the underlying molecular cause. Combined application of two independent library preparation and sequencing techniques, however, may still be indicated for Mendelian diseases with extensive genetic heterogeneity in order to further increase diagnostic sensitivity.

Transmission disequilibrium of small CNVs in simplex autism.

We searched for disruptive, genic rare copy-number variants (CNVs) among 411 families affected by sporadic autism spectrum disorder (ASD) from the Simons Simplex Collection by using available exome sequence data and CoNIFER (Copy Number Inference from Exome Reads). Compared to high-density SNP microarrays, our approach yielded ∼2× more smaller genic rare CNVs. We found that affected probands inherited more CNVs than did their siblings (453 versus 394, p = 0.004; odds ratio [OR] = 1.19) and that the probands' CNVs affected more genes (921 versus 726, p = 0.02; OR = 1.30). These smaller CNVs (median size 18 kb) were transmitted preferentially from the mother (136 maternal versus 100 paternal, p = 0.02), although this bias occurred irrespective of affected status. The excess burden of inherited CNVs among probands was driven primarily by sibling pairs with discordant social-behavior phenotypes (p < 0.0002, measured by Social Responsiveness Scale [SRS] score), which contrasts with families where the phenotypes were more closely matched or less extreme (p > 0.5). Finally, we found enrichment of brain-expressed genes unique to probands, especially in the SRS-discordant group (p = 0.0035). In a combined model, our inherited CNVs, de novo CNVs, and de novo single-nucleotide variants all independently contributed to the risk of autism (p < 0.05). Taken together, these results suggest that small transmitted rare CNVs play a role in the etiology of simplex autism. Importantly, the small size of these variants aids in the identification of specific genes as additional risk factors associated with ASD.

Single molecule molecular inversion probes for targeted, high-accuracy detection of low-frequency variation.

The detection and quantification of genetic heterogeneity in populations of cells is fundamentally important to diverse fields, ranging from microbial evolution to human cancer genetics. However, despite the cost and throughput advances associated with massively parallel sequencing, it remains challenging to reliably detect mutations that are present at a low relative abundance in a given DNA sample. Here we describe smMIP, an assay that combines single molecule tagging with multiplex targeted capture to enable practical and highly sensitive detection of low-frequency or subclonal variation. To demonstrate the potential of the method, we simultaneously resequenced 33 clinically informative cancer genes in eight cell line and 45 clinical cancer samples. Single molecule tagging facilitated extremely accurate consensus calling, with an estimated per-base error rate of 8.4 × 10(-6) in cell lines and 2.6 × 10(-5) in clinical specimens. False-positive mutations in the single molecule consensus base-calls exhibited patterns predominantly consistent with DNA damage, including 8-oxo-guanine and spontaneous deamination of cytosine. Based on mixing experiments with cell line samples, sensitivity for mutations above 1% frequency was 83% with no false positives. At clinically informative sites, we identified seven low-frequency point mutations (0.2%-4.7%), including BRAF p.V600E (melanoma, 0.2% alternate allele frequency), KRAS p.G12V (lung, 0.6%), JAK2 p.V617F (melanoma, colon, two lung, 0.3%-1.4%), and NRAS p.Q61R (colon, 4.7%). We anticipate that smMIP will be broadly adoptable as a practical and effective method for accurately detecting low-frequency mutations in both research and clinical settings.

Rapid and accurate large-scale genotyping of duplicated genes and discovery of interlocus gene conversions.

Over 900 genes have been annotated within duplicated regions of the human genome, yet their functions and potential roles in disease remain largely unknown. One major obstacle has been the inability to accurately and comprehensively assay genetic variation for these genes in a high-throughput manner. We developed a sequencing-based method for rapid and high-throughput genotyping of duplicated genes using molecular inversion probes designed to target unique paralogous sequence variants. We applied this method to genotype all members of two gene families, SRGAP2 and RH, among a diversity panel of 1,056 humans. The approach could accurately distinguish copy number in paralogs having up to ∼99.6% sequence identity, identify small gene-disruptive deletions, detect single-nucleotide variants, define breakpoints of unequal crossover and discover regions of interlocus gene conversion. The ability to rapidly and accurately genotype multiple gene families in thousands of individuals at low cost enables the development of genome-wide gene conversion maps and 'unlocks' many previously inaccessible duplicated genes for association with human traits.

PI3K/AKT pathway mutations cause a spectrum of brain malformations from megalencephaly to focal cortical dysplasia.

Malformations of cortical development containing dysplastic neuronal and glial elements, including hemimegalencephaly and focal cortical dysplasia, are common causes of intractable paediatric epilepsy. In this study we performed multiplex targeted sequencing of 10 genes in the PI3K/AKT pathway on brain tissue from 33 children who underwent surgical resection of dysplastic cortex for the treatment of intractable epilepsy. Sequencing results were correlated with clinical, imaging, pathological and immunohistological phenotypes. We identified mosaic activating mutations in PIK3CA and AKT3 in this cohort, including cancer-associated hotspot PIK3CA mutations in dysplastic megalencephaly, hemimegalencephaly, and focal cortical dysplasia type IIa. In addition, a germline PTEN mutation was identified in a male with hemimegalencephaly but no peripheral manifestations of the PTEN hamartoma tumour syndrome. A spectrum of clinical, imaging and pathological abnormalities was found in this cohort. While patients with more severe brain imaging abnormalities and systemic manifestations were more likely to have detected mutations, routine histopathological studies did not predict mutation status. In addition, elevated levels of phosphorylated S6 ribosomal protein were identified in both neurons and astrocytes of all hemimegalencephaly and focal cortical dysplasia type II specimens, regardless of the presence or absence of detected PI3K/AKT pathway mutations. In contrast, expression patterns of the T308 and S473 phosphorylated forms of AKT and in vitro AKT kinase activities discriminated between mutation-positive dysplasia cortex, mutation-negative dysplasia cortex, and non-dysplasia epilepsy cortex. Our findings identify PI3K/AKT pathway mutations as an important cause of epileptogenic brain malformations and establish megalencephaly, hemimegalencephaly, and focal cortical dysplasia as part of a single pathogenic spectrum.

Genome-wide association study and admixture mapping reveal new loci associated with total IgE levels in Latinos.

BACKGROUND: IgE is a key mediator of allergic inflammation, and its levels are frequently increased in patients with allergic disorders. OBJECTIVE: We sought to identify genetic variants associated with IgE levels in Latinos. METHODS: We performed a genome-wide association study and admixture mapping of total IgE levels in 3334 Latinos from the Genes-environments & Admixture in Latino Americans (GALA II) study. Replication was evaluated in 454 Latinos, 1564 European Americans, and 3187 African Americans from independent studies. RESULTS: We confirmed associations of 6 genes identified by means of previous genome-wide association studies and identified a novel genome-wide significant association of a polymorphism in the zinc finger protein 365 gene (ZNF365) with total IgE levels (rs200076616, P = 2.3 × 10(-8)). We next identified 4 admixture mapping peaks (6p21.32-p22.1, 13p22-31, 14q23.2, and 22q13.1) at which local African, European, and/or Native American ancestry was significantly associated with IgE levels. The most significant peak was 6p21.32-p22.1, where Native American ancestry was associated with lower IgE levels (P = 4.95 × 10(-8)). All but 22q13.1 were replicated in an independent sample of Latinos, and 2 of the peaks were replicated in African Americans (6p21.32-p22.1 and 14q23.2). Fine mapping of 6p21.32-p22.1 identified 6 genome-wide significant single nucleotide polymorphisms in Latinos, 2 of which replicated in European Americans. Another single nucleotide polymorphism was peak-wide significant within 14q23.2 in African Americans (rs1741099, P = 3.7 × 10(-6)) and replicated in non-African American samples (P = .011). CONCLUSION: We confirmed genetic associations at 6 genes and identified novel associations within ZNF365, HLA-DQA1, and 14q23.2. Our results highlight the importance of studying diverse multiethnic populations to uncover novel loci associated with total IgE levels.

Copy number variation detection and genotyping from exome sequence data.

While exome sequencing is readily amenable to single-nucleotide variant discovery, the sparse and nonuniform nature of the exome capture reaction has hindered exome-based detection and characterization of genic copy number variation. We developed a novel method using singular value decomposition (SVD) normalization to discover rare genic copy number variants (CNVs) as well as genotype copy number polymorphic (CNP) loci with high sensitivity and specificity from exome sequencing data. We estimate the precision of our algorithm using 122 trios (366 exomes) and show that this method can be used to reliably predict (94% overall precision) both de novo and inherited rare CNVs involving three or more consecutive exons. We demonstrate that exome-based genotyping of CNPs strongly correlates with whole-genome data (median r(2) = 0.91), especially for loci with fewer than eight copies, and can estimate the absolute copy number of multi-allelic genes with high accuracy (78% call level). The resulting user-friendly computational pipeline, CoNIFER (copy number inference from exome reads), can reliably be used to discover disruptive genic CNVs missed by standard approaches and should have broad application in human genetic studies of disease.