Recommendations for risk allele evidence curation, classification, and reporting from the ClinGen Low Penetrance/Risk Allele Working Group.

PURPOSE: Genetic variants at the low end of the penetrance spectrum have historically been challenging to interpret because their high population frequencies exceed the disease prevalence of the associated condition, leading to a lack of clear segregation between the variant and disease. There is currently substantial variation in the classification of these variants, and no formal classification framework has been widely adopted. The Clinical Genome Resource Low Penetrance/Risk Allele Working Group was formed to address these challenges and promote harmonization within the clinical community. METHODS: The work presented here is the product of internal and community Likert-scaled surveys in combination with expert consensus within the Working Group. RESULTS: We formally recognize risk alleles and low-penetrance variants as distinct variant classes from those causing highly penetrant disease that require special considerations regarding their clinical classification and reporting. First, we provide a preferred terminology for these variants. Second, we focus on risk alleles and detail considerations for reviewing relevant studies and present a framework for the classification these variants. Finally, we discuss considerations for clinical reporting of risk alleles. CONCLUSION: These recommendations support harmonized interpretation, classification, and reporting of variants at the low end of the penetrance spectrum.

Detecting regions of homozygosity improves the diagnosis of pathogenic variants and uniparental disomy in pediatric patients.

Chromosomal microarray analysis using single nucleotide polymorphism probes can detect regions of homozygosity (ROH). This confers a potential utility in revealing autosomal recessive (AR) diseases and uniparental disomy (UPD). Results of genetic testing among pediatric patients from 2015 to 2019 were evaluated. Diagnostic findings with detected ROH from large consecutive case series in the literature were reviewed. Of 2050 pediatric patients, 65 (3%) had one or more ROH and 31 (53%) had follow-up whole exome sequencing (WES) and methylation studies. Seven homozygous variants were detected and four of them from three patients (9.6%) were within the detected ROH and classified as pathogenic or likely pathogenic variants for AR diseases. One patient (3%) had segmental UPD15q for a diagnosis of Prader-Willi syndrome. Additive diagnostic yield from ROH reporting was at least 0.2% (4/2050) of pediatric patients. These results were consistent with findings from several large case series reported in the literature. Detecting ROH had an estimated baseline predictive value of 10% for AR diseases and 3% for UPD. Consanguinity revealed by multiple ROH was a strong predictor for AR diseases. These results provide evidence for genetic counseling and recommendation of follow-up WES and methylation studies for pediatric patients reported with ROH.

Common genetic variants on 5p14.1 associate with autism spectrum disorders.

Autism spectrum disorders (ASDs) represent a group of childhood neurodevelopmental and neuropsychiatric disorders characterized by deficits in verbal communication, impairment of social interaction, and restricted and repetitive patterns of interests and behaviour. To identify common genetic risk factors underlying ASDs, here we present the results of genome-wide association studies on a cohort of 780 families (3,101 subjects) with affected children, and a second cohort of 1,204 affected subjects and 6,491 control subjects, all of whom were of European ancestry. Six single nucleotide polymorphisms between cadherin 10 (CDH10) and cadherin 9 (CDH9)-two genes encoding neuronal cell-adhesion molecules-revealed strong association signals, with the most significant SNP being rs4307059 (P = 3.4 x 10(-8), odds ratio = 1.19). These signals were replicated in two independent cohorts, with combined P values ranging from 7.4 x 10(-8) to 2.1 x 10(-10). Our results implicate neuronal cell-adhesion molecules in the pathogenesis of ASDs, and represent, to our knowledge, the first demonstration of genome-wide significant association of common variants with susceptibility to ASDs.

Evaluation of copy number variations reveals novel candidate genes in autism spectrum disorder-associated pathways.

Autism spectrum disorders (ASDs) are highly heritable, yet relatively few associated genetic loci have been replicated. Copy number variations (CNVs) have been implicated in autism; however, the majority of loci contribute to <1% of the disease population. Therefore, independent studies are important to refine associated CNV regions and discover novel susceptibility genes. In this study, a genome-wide SNP array was utilized for CNV detection by two distinct algorithms in a European ancestry case-control data set. We identify a significantly higher burden in the number and size of deletions, and disrupting more genes in ASD cases. Moreover, 18 deletions larger than 1 Mb were detected exclusively in cases, implicating novel regions at 2q22.1, 3p26.3, 4q12 and 14q23. Case-specific CNVs provided further evidence for pathways previously implicated in ASDs, revealing new candidate genes within the GABAergic signaling and neural development pathways. These include DBI, an allosteric binder of GABA receptors, GABARAPL1, the GABA receptor-associated protein, and SLC6A11, a postsynaptic GABA transporter. We also identified CNVs in COBL, deletions of which cause defects in neuronal cytoskeleton morphogenesis in model vertebrates, and DNER, a neuron-specific Notch ligand required for cerebellar development. Moreover, we found evidence of genetic overlap between ASDs and other neurodevelopmental and neuropsychiatric diseases. These genes include glutamate receptors (GRID1, GRIK2 and GRIK4), synaptic regulators (NRXN3, SLC6A8 and SYN3), transcription factor (ZNF804A) and RNA-binding protein FMR1. Taken together, these CNVs may be a few of the missing pieces of ASD heritability and lead to discovering novel etiological mechanisms.

Maternal Uniparental Isodisomy of Chromosome 2 Leading to Homozygous Variants in SPR and ZNF142 : A Case Report and Review of the UPD2 Literature.

We report a 4-year-old girl with neurodevelopmental abnormalities who has maternal uniparental isodisomy of chromosome 2 leading to homozygosity for a likely pathogenic variant in SPR , and a variant of uncertain significance in ZNF142 . Biallelic pathogenic variants in SPR lead to sepiapterin reductase deficiency (SRD), a dopa-responsive dystonia. Pathogenic variants in ZNF142 are associated with an autosomal recessive neurodevelopmental disorder characterized by impaired speech and hyperkinetic movements, which has significant clinical overlap with SRD. Our patient showed dramatic improvement in motor skills after treatment with levodopa. We also reviewed 67 published reports of uniparental disomy of chromosome 2 (UPD2) associated with various clinical outcomes. These include autosomal recessive disorders associated with loci on chromosome 2, infants with UPD2 whose gestations were associated with confined placental mosaicism for trisomy 2 leading to intrauterine growth restriction with good postnatal catchup growth, and normal phenotypes in children and adults with an incidental finding of either maternal or paternal UPD2. These latter reports provide support for the conclusion that genes located on chromosome 2 are not subject to imprinting. We also explore the mechanisms giving rise to UPD2.

KIF6 Trp719Arg Genetic Variant Increases Risk for Thoracic Aortic Dissection.

BACKGROUND: KIF6 (kinesin family member 6), a protein coded by the KIF6 gene, serves an important intracellular function to transport organelles along microtubules. In a pilot study, we found that a common KIF6 Trp719Arg variant increased the propensity of thoracic aortic aneurysms (TAA) to suffer dissection (AD). The present study aims for a definite investigation of the predictive ability of KIF6 719Arg vis à vis AD. Confirmatory findings would enhance natural history prediction in TAA. METHODS: 1108 subjects (899 aneurysm and 209 dissection patients) had KIF6 719Arg variant status determined. RESULTS: The 719Arg variant in the KIF6 gene correlated strongly with occurrence of AD. Specifically, KIF6 719Arg positivity (homozygous or heterozygous) was substantially more common in dissectors (69.8%) than non-dissectors (58.5%) (p = 0.003). Odds ratios (OR) for suffering aortic dissection ranged from 1.77 to 1.94 for Arg carriers in various dissection categories. These high OR associations were noted for both ascending and descending aneurysms and for homozygous and heterozygous Arg variant patients. The rate of aortic dissection over time was significantly higher for carriers of the Arg allele (p = 0.004). Additionally, Arg allele carriers were more likely to reach the combined endpoint of dissection or death (p = 0.03). CONCLUSIONS: We demonstrate the marked adverse impact of the 719Arg variant of the KIF6 gene on the likelihood that a TAA patient will suffer aortic dissection. Clinical assessment of the variant status of this molecularly important gene may provide a valuable "non-size" criterion to enhance surgical decision making above and beyond the currently used metric of aortic size (diameter).

Genetic analysis of biological pathway data through genomic randomization.

Genome Wide Association Studies (GWAS) are a standard approach for large-scale common variation characterization and for identification of single loci predisposing to disease. However, due to issues of moderate sample sizes and particularly multiple testing correction, many variants of smaller effect size are not detected within a single allele analysis framework. Thus, small main effects and potential epistatic effects are not consistently observed in GWAS using standard analytical approaches that consider only single SNP alleles. Here, we propose unique methodology that aggregates variants of interest (for example, genes in a biological pathway) using GWAS results. Multiple testing and type I error concerns are minimized using empirical genomic randomization to estimate significance. Randomization corrects for common pathway-based analysis biases, such as SNP coverage and density, linkage disequilibrium, gene size and pathway size. Pathway Analysis by Randomization Incorporating Structure (PARIS) applies this randomization and in doing so directly accounts for linkage disequilibrium effects. PARIS is independent of association analysis method and is thus applicable to GWAS datasets of all study designs. Using the KEGG database as an example, we apply PARIS to the publicly available Autism Genetic Resource Exchange GWAS dataset, revealing pathways with a significant enrichment of positive association results.

Microduplications in an autism multiplex family narrow the region of susceptibility for developmental disorders on 15q24 and implicate 7p21.

Copy number variations (CNVs) play a crucial role in the intricate genetics of autism spectrum disorders. A region on chromosome 15q24 vulnerable to both deletions and duplications has been previously implicated in a range of phenotypes including autism, Asperger's syndrome, delayed development, and mild to severe mental retardation. Prior studies have delineated a minimal critical region of approximately 1.33 Mb. In this study, a multiplex autism family was evaluated for CNVs using genotyping data from the Illumina 1 M BeadChip and analyzed with the PennCNV algorithm. Variants were then identified that co-segregate with autism features in this family. Here, we report autistic first cousins who carry two microduplications concordant with disease. Both duplications were inherited maternally and found to be identical by descent. The first is an approximately 10,000 base pair microduplication within the minimal region on 15q24 that falls across a single gene, ubiquitin-like 7. This is the smallest duplication in the region to result in a neuropsychiatric disorder, potentially narrowing the critical region for susceptibility to developmental and autism spectrum disorders. The second is a novel, 352 kb tandem duplication on 7p21 that replicates part of the neurexophilin 1 and islet cell autoantigen 1 genes. The breakpoint junction falls within the intronic regions of these genes and demonstrates a microhomology of four base pairs. Each of these microduplications may contribute to the complex etiology of autism spectrum disorders.

Variation in WNT genes is associated with non-syndromic cleft lip with or without cleft palate.

Non-syndromic cleft lip with or without cleft palate (NSCLP) is a common birth defect. Genetic and environmental factors have been causally implicated and studies have begun to delineate genetic contributions. The Wnt genes are involved in regulating mid-face development and upper lip fusion and are therefore strong candidates for an etiological role in NSCLP. Furthermore, the clf1 region in A/WyN clefting susceptible mice contains the Wnt3 and Wnt9B genes. To assess the role of the Wnt family of genes in NSCLP, we interrogated seven Wnt genes (Wnt3, Wnt3A, Wnt5A, Wnt7A, Wnt8A, Wnt9B and Wnt11) in our well-defined NSCLP dataset. Thirty-eight single nucleotide polymorphisms were genotyped in 132 multiplex NSCLP families and 354 simplex parent-child trios. In the entire dataset, single-nucleotide polymorphisms (SNPs) in three genes, Wnt3A (P = 0.006), Wnt 5A (P = 0.002) and Wnt11 (P = 0.0001) were significantly associated with NSCLP after correction for multiple testing. When stratified by ethnicity, the strongest associations were found for SNPs in Wnt3A (P = 0.0007), Wnt11 (P = 0.0012) and Wnt8A (P = 0.0013). Multiple haplotypes in Wnt genes were associated with NSCLP, and gene-gene interactions were observed between Wnt3A and both Wnt3 and Wnt5A (P = 0.004 and P = 0.039, respectively). This data suggests that alteration in Wnt gene function may perturb formation and/or fusion of the facial processes and predispose to NSCLP.

A genome-wide association study of autism reveals a common novel risk locus at 5p14.1.

Although autism is one of the most heritable neuropsychiatric disorders, its underlying genetic architecture has largely eluded description. To comprehensively examine the hypothesis that common variation is important in autism, we performed a genome-wide association study (GWAS) using a discovery dataset of 438 autistic Caucasian families and the Illumina Human 1M beadchip. 96 single nucleotide polymorphisms (SNPs) demonstrated strong association with autism risk (p-value < 0.0001). The validation of the top 96 SNPs was performed using an independent dataset of 487 Caucasian autism families genotyped on the 550K Illumina BeadChip. A novel region on chromosome 5p14.1 showed significance in both the discovery and validation datasets. Joint analysis of all SNPs in this region identified 8 SNPs having improved p-values (3.24E-04 to 3.40E-06) than in either dataset alone. Our findings demonstrate that in addition to multiple rare variations, part of the complex genetic architecture of autism involves common variation.

Determinants of skeletal age deviation in a cross-sectional study.

OBJECTIVE: Skeletal age deviation (SAD) is associated with bone mass and fracture risk in children, but factors determining this are unknown. The aim of this population-based cross-sectional study was to describe the factors associated with SAD. METHODS: A convenience sample of 640 male and female children aged 7-17 yr was studied. All were assessed for body composition (dual-energy x-ray absorptiometry), diet, strength, dexterity, habitual physical activity, sunlight exposure, smoking, and medication use. Skeletal age was assigned using the Tanner-Whitehouse-2 method. RESULTS: Subjects with a SAD greater than the 75th percentile had significantly higher height, weight, and Tanner stage compared with all other subjects. Bone-free lean mass, fat mass, and grip strength were positively associated with SAD. In multivariate analysis, ever smoking and use of inhaled corticosteroids were negatively associated with SAD, whereas milk drinking was positively associated with SAD. There was no significant association between sunlight exposure, television watching, light, or strenuous exercise and SAD. CONCLUSIONS: The results of this study should be regarded as hypothesis generating but are biologically plausible and suggest that body composition, strength, diet, ever smoking, and inhaled corticosteroid use may be determinants of bone maturity relative to age and thus affect fracture risk in children. However, more studies are necessary to explore other determinants of SAD such as genetic and perinatal factors and whether SAD influences peak bone mass.

Bone density interpretation and relevance in Caucasian children aged 9-17 years of age: insights from a population-based fracture study.

The interpretation of bone density measurement in children is difficult due to a number of factors including rapid change in body size and uncertain clinical significance of bone density in children. This study asked two questions. (1) Is there a preferred bone density measurement site or type for fracture risk in children? (2) What is the best way to interpret bone density in children? This population-based case control study included 321 upper limb fracture cases and 321 class- and sex- matched randomly selected controls. Bone density at the hip, spine, and total body (including the arm) was measured by a Hologic QDR2000 densitometer (Waltham, MA) and examined as bone area (BA), bone mineral content (BMC), bone mineral density (BMD), bone mineral apparent density (BMAD), and BMC/lean mass (BMCLM). The only dual-energy X-ray absorptiometry (DXA) variables that were consistently associated with fracture risk in both boys and girls were spine BMD and BMAD for total upper limb fractures, and spine and hip BMAD for wrist and forearm fractures. No significant associations were observed for BA and BMCLM and inconsistent associations for BMC and other BMD sites. Five-yr fracture risk varied from 15-24% depending on site and gender in a child with a Z-score of -3. In the controls, all DXA variables were associated with age, height, and weight, but the weakest associations were with BMAD. In conclusion, in this study the spine BMAD had the strongest and most consistent association with upper limb fracture risk in children. The associations with age and body size imply that age specific Z-scores will be the most convenient for interpretation of DXA measures in children. Five-yr wrist and forearm fracture risk has potential as a clinical endpoint of immediate relevance.

Synovial haemangioma of the knee: a case report.

Synovial haemangioma is a rare but important cause of knee symptoms, which, when undiagnosed, can lead to significant morbidity. Diagnosis is frequently difficult and delayed. We report on a case of synovial haemangioma, which demonstrates the difficulties inherent in diagnosis and the morbidity associated with diagnostic delay, in a young woman. Magnetic resonance imaging (MRI) is a useful tool for diagnosis, but detection on MRI can also be problematic, as shown by this case, demonstrating the need for greater awareness of this condition by both clinicians and radiologists. Arthroscopy is important in both the diagnosis and treatment of these lesions.

Skeletal age deviation assessed by the Tanner-Whitehouse 2 method is associated with bone mass and fracture risk in children.

The aim of this population-based case-control study was to describe the association among skeletal age deviation (SAD), bone density, and upper limb fracture risk in male and female children aged 9-16 years. A total of 321 fracture cases and 321 randomly selected individually matched controls were studied. Skeletal age was assessed by standard atlas. Bone mineral density (BMD) was measured by dual-energy X-ray absorptiometry (DXA) and metacarpal index (MI). There were no significant differences in mean skeletal age or chronological age between fracture cases and controls. However, SAD was associated with total, hand, and female fracture risk (all P<0.05). The fracture associations became nonsignificant after adjustment for BMD and MI in all subgroups with the exception of hand fractures (OR, 0.67/year; 95% CI, 0.47-0.96). SAD was also positively associated with BMD at all sites (r=0.33-0.35, all P<0.05) and MI (r=0.20, P<0.05). The strength of association reduced but remained significant at most sites after adjustment for body size, maturity, age, and sex. In conclusion, SAD is positively associated with measures of bone strength and negatively associated with upper limb fracture risk (especially those of the hand) in children. SAD is simple to measure and gives additional information regarding bone health and fracture risk in children.

Extension of multifactor dimensionality reduction for identifying multilocus effects in the GAW14 simulated data.

The multifactor dimensionality reduction (MDR) is a model-free approach that can identify gene x gene or gene x environment effects in a case-control study. Here we explore several modifications of the MDR method. We extended MDR to provide model selection without crossvalidation, and use a chi-square statistic as an alternative to prediction error (PE). We also modified the permutation test to provide different levels of stringency. The extended MDR (EMDR) includes three permutation tests (fixed, non-fixed, and omnibus) to obtain p-values of multilocus models. The goal of this study was to compare the different approaches implemented in the EMDR method and evaluate the ability to identify genetic effects in the Genetic Analysis Workshop 14 simulated data. We used three replicates from the simulated family data, generating matched pairs from family triads. The results showed: 1) chi-square and PE statistics give nearly consistent results; 2) results of EMDR without cross-validation matched that of EMDR with 10-fold cross-validation; 3) the fixed permutation test reports false-positive results in data from loci unrelated to the disease, but the non-fixed and omnibus permutation tests perform well in preventing false positives, with the omnibus test being the most conservative. We conclude that the non-cross-validation test can provide accurate results with the advantage of high efficiency compared to 10-cross-validation, and the non-fixed permutation test provides a good compromise between power and false-positive rate.

Television, computer, and video viewing; physical activity; and upper limb fracture risk in children: a population-based case control study.

UNLABELLED: The effect of physical activity on upper limb fractures was examined in this population-based case control study with 321 age- and gender-matched pairs. Sports participation increased fracture risk in boys and decreased risk in girls. Television viewing had a deleterious dose response association with wrist and forearm fractures while light physical activity was protective. INTRODUCTION: The aim of this population-based case control study was to examine the association between television, computer, and video viewing; types and levels of physical activity; and upper limb fractures in children 9-16 years of age. MATERIALS AND METHODS: A total of 321 fracture cases and 321 randomly selected individually matched controls were studied. Television, computer, and video viewing and types and levels of physical activity were determined by interview-administered questionnaire. Bone strength was assessed by DXA and metacarpal morphometry. RESULTS: In general, sports participation increased total upper limb fracture risk in boys and decreased risk in girls. Gender-specific risk estimates were significantly different for total, contact, noncontact, and high-risk sports participation as well as four individual sports (soccer, cricket, surfing, and swimming). In multivariate analysis, time spent television, computer, and video viewing in both sexes was positively associated with wrist and forearm fracture risk (OR 1.6/category, 95% CI: 1.1-2.2), whereas days involved in light physical activity participation decreased fracture risk (OR 0.8/category, 95% CI: 0.7-1.0). Sports participation increased hand (OR 1.5/sport, 95% CI: 1.1-2.0) and upper arm (OR 29.8/sport, 95% CI: 1.7-535) fracture risk in boys only and decreased wrist and forearm fracture risk in girls only (OR 0.5/sport, 95% CI: 0.3-0.9). Adjustment for bone density and metacarpal morphometry did not alter these associations. CONCLUSION: There is gender discordance with regard to sports participation and fracture risk in children, which may reflect different approaches to sport. Importantly, television, computer, and video viewing has a dose-dependent association with wrist and forearm fractures, whereas light physical activity is protective. The mechanism is unclear but may involve bone-independent factors, or less likely, changes in bone quality not detected by DXA or metacarpal morphometry.

The association between bone mineral density, metacarpal morphometry, and upper limb fractures in children: a population-based case-control study.

The aim of this population-based case-control study was to examine the association between bone mass and upper limb fractures in children aged 9-16 yr. Areal bone mineral density and bone mineral apparent density (BMAD) were measured by both dual energy absorptiometry (DXA) and metacarpal index (MI) by hand radiograph. A total of 321 fracture cases and 321 randomly selected individually matched controls were studied. For all fractures, cases had lower DXA measures at all sites (1.1-3.3%; all P < 0.05). A larger reduction was observed for those with wrist and forearm fractures (1.2-4.5%; all P < 0.05, except total body BMAD) but not other upper limb fractures (hand, -1.6 to +1.2%; upper arm: 0.9-4.8%; all P > 0.05). For metacarpal measures, cases had a thinner cortical width and lower MI for wrist and forearm fractures only. In multivariate modeling, both spine BMAD (odds ratio, 1.4/SD reduction) and MI (odds ratio, 1.5/SD reduction) remained statistically significant predictors of wrist and forearm fractures. In conclusion, both DXA measures and MI are independently associated with wrist and forearm but not other upper limb fractures. The magnitude of this association is somewhat weaker than in adults but suggests that optimizing age-appropriate bone mass will lessen the risk of fracture in children.

A de novo 10.79 Mb interstitial deletion at 2q13q14.2 involving PAX8 causing hypothyroidism and mullerian agenesis: a novel case report and literature review.

Reports of interstitial deletions involving proximal long arm of chromosome 2 are limited. Based on early chromosomal analysis studies, the phenotypic consequence of deletions at the ancestral chromosome fusion site at chromosome 2q13q14.1 remains unclear. A recurrent 1.71 Mb deletion at 2q13 has recently been proposed as a new genomic disorder, associated with an increased risk of intellectual disability and craniofacial dysmorphism. Herein, we report the case of a 12 year-old girl with unique clinical features including global developmental delay, mullerian agenesis, and hypothyroidism associated with a normal size and position of the thyroid gland, as well as negative thyroid antibodies. Microarray-based comparative genomic hybridization study revealed a de novo 10.79 Mb deletion at 2q13q14.2 (111,548,932-122,336,492), which involves more than 88 UCSC genes, 38 of which are OMIM genes, 7 of which are disease-causing and 3 of which (including GLI2, IL1B and PAX8) show a dominant inheritance pattern.. Interestingly, PAX8 (chr2:113,973,574-114,036,498), a member of the paired-box gene family, is essential for the formation of thyroxine-producing follicular cells. Autosomal dominant transmission of congenital thyroid hypoplasia due to loss-of-function mutation of PAX8 suggests a possible haploinsufficiency effect. Additionally, PAX8 is also expressed in the tissue primordia that form both the mullerian duct derivatives and the upper urinary tracts. A recent study has associated a novel PAX8 mutation with a severe form of hypothyroidism and abnormalities in the urogenital tract. Taken together, the unique clinical manifestation seen in this patient could be attributed to the heterozygous deletion of PAX8 gene. A prospective investigation is merited to fully evaluate the pathogenic effect of the interstitial deletion of 2q13q14.2.

Evidence of novel fine-scale structural variation at autism spectrum disorder candidate loci.

BACKGROUND: Autism spectrum disorders (ASD) represent a group of neurodevelopmental disorders characterized by a core set of social-communicative and behavioral impairments. Gamma-aminobutyric acid (GABA) is the major inhibitory neurotransmitter in the brain, acting primarily via the GABA receptors (GABR). Multiple lines of evidence, including altered GABA and GABA receptor expression in autistic patients, indicate that the GABAergic system may be involved in the etiology of autism. METHODS: As copy number variations (CNVs), particularly rare and de novo CNVs, have now been implicated in ASD risk, we examined the GABA receptors and genes in related pathways for structural variation that may be associated with autism. We further extended our candidate gene set to include 19 genes and regions that had either been directly implicated in the autism literature or were directly related (via function or ancestry) to these primary candidates. For the high resolution CNV screen we employed custom-designed 244 k comparative genomic hybridization (CGH) arrays. Collectively, our probes spanned a total of 11 Mb of GABA-related and additional candidate regions with a density of approximately one probe every 200 nucleotides, allowing a theoretical resolution for detection of CNVs of approximately 1 kb or greater on average. One hundred and sixty-eight autism cases and 149 control individuals were screened for structural variants. Prioritized CNV events were confirmed using quantitative PCR, and confirmed loci were evaluated on an additional set of 170 cases and 170 control individuals that were not included in the original discovery set. Loci that remained interesting were subsequently screened via quantitative PCR on an additional set of 755 cases and 1,809 unaffected family members. RESULTS: Results include rare deletions in autistic individuals at JAKMIP1, NRXN1, Neuroligin4Y, OXTR, and ABAT. Common insertion/deletion polymorphisms were detected at several loci, including GABBR2 and NRXN3. Overall, statistically significant enrichment in affected vs. unaffected individuals was observed for NRXN1 deletions. CONCLUSIONS: These results provide additional support for the role of rare structural variation in ASD.

The expanding role of MBD genes in autism: identification of a MECP2 duplication and novel alterations in MBD5, MBD6, and SETDB1.

The methyl-CpG-binding domain (MBD) gene family was first linked to autism over a decade ago when Rett syndrome, which falls under the umbrella of autism spectrum disorders (ASDs), was revealed to be predominantly caused by MECP2 mutations. Since that time, MECP2 alterations have been recognized in idiopathic ASD patients by us and others. Individuals with deletions across the MBD5 gene also present with ASDs, impaired speech, intellectual difficulties, repetitive behaviors, and epilepsy. These findings suggest that further investigations of the MBD gene family may reveal additional associations related to autism. We now describe the first study evaluating individuals with ASD for rare variants in four autosomal MBD family members, MBD5, MBD6, SETDB1, and SETDB2, and expand our initial screening in the MECP2 gene. Each gene was sequenced over all coding exons and evaluated for copy number variations in 287 patients with ASD and an equal number of ethnically matched control individuals. We identified 186 alterations through sequencing, approximately half of which were novel (96 variants, 51.6%). We identified 17 ASD specific, nonsynonymous variants, four of which were concordant in multiplex families: MBD5 Tyr1269Cys, MBD6 Arg883Trp, MECP2 Thr240Ser, and SETDB1 Pro1067del. Furthermore, a complex duplication spanning of the MECP2 gene was identified in two brothers who presented with developmental delay and intellectual disability. From our studies, we provide the first examples of autistic patients carrying potentially detrimental alterations in MBD6 and SETDB1, thereby demonstrating that the MBD gene family potentially plays a significant role in rare and private genetic causes of autism.