Deletion 17q12 is a recurrent copy number variant that confers high risk of autism and schizophrenia

Autism spectrum disorders (ASD) and schizophrenia are neurodevelopmental disorders for which recent evidence indicates an important etiologic role for rare copy number variants (CNVs) and suggests common genetic mechanisms. We performed cytogenomic array analysis in a discovery sample of patients with neurodevelopmental disorders referred for clinical testing. We detected a recurrent 1.4 Mb deletion at 17q12, which harbors HNF1B, the gene responsible for renal cysts and diabetes syndrome (RCAD), in 18/15,749 patients, including several with ASD, but 0/4,519 controls. We identified additional shared phenotypic features among nine patients available for clinical assessment, including macrocephaly, characteristic facial features, renal anomalies, and neurocognitive impairments. In a large follow-up sample, the same deletion was identified in 2/1,182 ASD/neurocognitive impairment and in 4/6,340 schizophrenia patients, but in 0/47,929 controls (corrected p = 7.37 × 10⁻⁵). These data demonstrate that deletion 17q12 is a recurrent, pathogenic CNV that confers a very high risk for ASD and schizophrenia and show that one or more of the 15 genes in the deleted interval is dosage sensitive and essential for normal brain development and function. In addition, the phenotypic features of patients with this CNV are consistent with a contiguous gene syndrome that extends beyond RCAD, which is caused by HNF1B mutations only.     

Genome scan for Tourette disorder in affected-sibling-pair and multigenerational families

Tourette disorder (TD) is a neuropsychiatric disorder with a complex mode of inheritance and is characterized by multiple waxing and waning motor and phonic tics. This article reports the results of the largest genetic linkage study yet undertaken for TD. The sample analyzed includes 238 nuclear families yielding 304 "independent" sibling pairs and 18 separate multigenerational families, for a total of 2,040 individuals. A whole-genome screen with the use of 390 microsatellite markers was completed. Analyses were completed using two diagnostic classifications: (1) only individuals with TD were included as affected and (2) individuals with either TD or chronic-tic (CT) disorder were included as affected. Strong evidence of linkage was observed for a region on chromosome 2p (-log P = 4.42, P = 3.8 x 10(-5) in the analyses that included individuals with TD or CT disorder as affected. Results in several other regions also provide moderate evidence (-log P >2.0) of additional susceptibility loci for TD.     

Genetic diversity and subdivision of 57 European and Middle-Eastern sheep breeds

The population structure and genetic diversity of 57 European and Middle Eastern marginal and cosmopolitan sheep breeds from 15 countries were analysed by typing 31 microsatellite markers. Mean unbiased expected heterozygosities ranged from 0.63 in British Exmoor Horn to 0.77 in Albanian Ruda. South-eastern European and Middle-Eastern sheep breeds were significantly more variable than northwestern and western European breeds. An overall heterozygote deficiency (f) across all loci was observed (P < 0.001), while genetic differentiation (theta) was 5.7%. Principal component analysis and Bayesian model-based clustering indicate a south-east to north-west cline, but also revealed distinct groups of Middle-Eastern fat-tailed sheep, south-eastern European sheep and north-western/western European sheep. Within the last group, two less-distinct clusters comprised the Merino-type and Alpine breeds respectively. The incomplete demarcations of most clusters probably reflects cross-breeding and/or upgrading.     

Improving Phenotypic Prediction by Combining Genetic and Epigenetic Associations

We tested whether DNA-methylation profiles account for inter-individual variation in body mass index (BMI) and height and whether they predict these phenotypes over and above genetic factors. Genetic predictors were derived from published summary results from the largest genome-wide association studies on BMI (n ∼ 350,000) and height (n ∼ 250,000) to date. We derived methylation predictors by estimating probe-trait effects in discovery samples and tested them in external samples. Methylation profiles associated with BMI in older individuals from the Lothian Birth Cohorts (LBCs, n = 1,366) explained 4.9% of the variation in BMI in Dutch adults from the LifeLines DEEP study (n = 750) but did not account for any BMI variation in adolescents from the Brisbane Systems Genetic Study (BSGS, n = 403). Methylation profiles based on the Dutch sample explained 4.9% and 3.6% of the variation in BMI in the LBCs and BSGS, respectively. Methylation profiles predicted BMI independently of genetic profiles in an additive manner: 7%, 8%, and 14% of variance of BMI in the LBCs were explained by the methylation predictor, the genetic predictor, and a model containing both, respectively. The corresponding percentages for LifeLines DEEP were 5%, 9%, and 13%, respectively, suggesting that the methylation profiles represent environmental effects. The differential effects of the BMI methylation profiles by age support previous observations of age modulation of genetic contributions. In contrast, methylation profiles accounted for almost no variation in height, consistent with a mainly genetic contribution to inter-individual variation. The BMI results suggest that combining genetic and epigenetic information might have greater utility for complex-trait prediction.     

Recessive Osteogenesis Imperfecta Caused by Missense Mutations in SPARC

Secreted protein, acidic, cysteine-rich (SPARC) is a glycoprotein that binds to collagen type I and other proteins in the extracellular matrix. Using whole-exome sequencing to identify the molecular defect in two unrelated girls with severe bone fragility and a clinical diagnosis of osteogenesis imperfecta type IV, we identified two homozygous variants in SPARC (GenBank: {"type":"entrez-nucleotide","attrs":{"text":"NM_003118.3","term_id":"365777426","term_text":"NM_003118.3"}}NM_003118.3; c.497G>A [p.Arg166His] in individual 1; c.787G>A [p.Glu263Lys] in individual 2). Published modeling and site-directed mutagenesis studies had previously shown that the residues substituted by these mutations form an intramolecular salt bridge in SPARC and are essential for the binding of SPARC to collagen type I. The amount of SPARC secreted by skin fibroblasts was reduced in individual 1 but appeared normal in individual 2. The migration of collagen type I alpha chains produced by these fibroblasts was mildly delayed on SDS-PAGE gel, suggesting some overmodification of collagen during triple helical formation. Pulse-chase experiments showed that collagen type I secretion was mildly delayed in skin fibroblasts from both individuals. Analysis of an iliac bone sample from individual 2 showed that trabecular bone was hypermineralized on the material level. In conclusion, these observations show that homozygous mutations in SPARC can give rise to severe bone fragility in humans.     

Meta-analysis of Correlated Traits via Summary Statistics from GWASs with an Application in Hypertension

Genome-wide association studies (GWASs) have identified many genetic variants underlying complex traits. Many detected genetic loci harbor variants that associate with multiple—even distinct—traits. Most current analysis approaches focus on single traits, even though the final results from multiple traits are evaluated together. Such approaches miss the opportunity to systemically integrate the phenome-wide data available for genetic association analysis. In this study, we propose a general approach that can integrate association evidence from summary statistics of multiple traits, either correlated, independent, continuous, or binary traits, which might come from the same or different studies. We allow for trait heterogeneity effects. Population structure and cryptic relatedness can also be controlled. Our simulations suggest that the proposed method has improved statistical power over single-trait analysis in most of the cases we studied. We applied our method to the Continental Origins and Genetic Epidemiology Network (COGENT) African ancestry samples for three blood pressure traits and identified four loci (CHIC2, HOXA-EVX1, IGFBP1/IGFBP3, and CDH17; p < 5.0 × 10⁻⁸) associated with hypertension-related traits that were missed by a single-trait analysis in the original report. Six additional loci with suggestive association evidence (p < 5.0 × 10⁻⁷) were also observed, including CACNA1D and WNT3. Our study strongly suggests that analyzing multiple phenotypes can improve statistical power and that such analysis can be executed with the summary statistics from GWASs. Our method also provides a way to study a cross phenotype (CP) association by using summary statistics from GWASs of multiple phenotypes.     

RTTN Mutations Cause Primary Microcephaly and Primordial Dwarfism in Humans

Primary microcephaly is a developmental brain anomaly that results from defective proliferation of neuroprogenitors in the germinal periventricular zone. More than a dozen genes are known to be mutated in autosomal-recessive primary microcephaly in isolation or in association with a more generalized growth deficiency (microcephalic primordial dwarfism), but the genetic heterogeneity is probably more extensive. In a research protocol involving autozygome mapping and exome sequencing, we recruited a multiplex consanguineous family who is affected by severe microcephalic primordial dwarfism and tested negative on clinical exome sequencing. Two candidate autozygous intervals were identified, and the second round of exome sequencing revealed a single intronic variant therein (c.2885+8A>G [p.Ser963^∗] in RTTN exon 23). RT-PCR confirmed that this change creates a cryptic splice donor and thus causes retention of the intervening 7 bp of the intron and leads to premature truncation. On the basis of this finding, we reanalyzed the exome file of a second consanguineous family affected by a similar phenotype and identified another homozygous change in RTTN as the likely causal mutation. Combined linkage analysis of the two families confirmed that RTTN maps to the only significant linkage peak. Finally, through international collaboration, a Canadian multiplex family affected by microcephalic primordial dwarfism and biallelic mutation of RTTN was identified. Our results expand the phenotype of RTTN-related disorders, hitherto limited to polymicrogyria, to include microcephalic primordial dwarfism with a complex brain phenotype involving simplified gyration.     

Genome-wide Analysis of Body Proportion Classifies Height-Associated Variants by Mechanism of Action and Implicates Genes Important for Skeletal Development

Human height is a composite measurement, reflecting the sum of leg, spine, and head lengths. Many common variants influence total height, but the effects of these or other variants on the components of height (body proportion) remain largely unknown. We studied sitting height ratio (SHR), the ratio of sitting height to total height, to identify such effects in 3,545 African Americans and 21,590 individuals of European ancestry. We found that SHR is heritable: 26% and 39% of the total variance of SHR can be explained by common variants in European and African Americans, respectively, and global European admixture is negatively correlated with SHR in African Americans (r² ≈ 0.03). Six regions reached genome-wide significance (p < 5 × 10⁻⁸) for association with SHR and overlapped biological candidate genes, including TBX2 and IGFBP3. We found that 130 of 670 height-associated variants are nominally associated (p < 0.05) with SHR, more than expected by chance (p = 5 × 10⁻⁴⁰). At these 130 loci, the height-increasing alleles are associated with either a decrease (71 loci) or increase (59 loci) in SHR, suggesting that different height loci disproportionally affect either leg length or spine/head length. Pathway analyses via DEPICT revealed that height loci affecting SHR, and especially those affecting leg length, show enrichment of different biological pathways (e.g., bone/cartilage/growth plate pathways) than do loci with no effect on SHR (e.g., embryonic development). These results highlight the value of using a pair of related but orthogonal phenotypes, in this case SHR with height, as a prism to dissect the biology underlying genetic associations in polygenic traits and diseases.     

Hi-Plex targeted sequencing is effective using DNA derived from archival dried blood spots

Many genetic epidemiology resources have collected dried blood spots (predominantly as Guthrie Cards) as an economical and efficient means of archiving sources of DNA, conferring great value to genetic screening methods that are compatible with this medium. We applied Hi-Plex to screen the breast cancer predisposition gene PALB2 in 93 Guthrie Card-derived DNA specimens previously characterized for PALB2 genetic variants via DNA derived from lymphoblastoid cell lines, whole blood, and buffy coat. Of the 93 archival Guthrie Card-derived DNAs, 92 (99%) were processed successfully and sequenced using approximately half of a MiSeq run. From these 92 DNAs, all 59 known variants were detected and no false-positive variant calls were yielded. Fully 98.13% of amplicons (5417/5520) were represented within 15-fold of the median coverage (2786 reads), and 99.98% of amplicons (5519/5520) were represented at a depth of 10 read-pairs or greater. With Hi-Plex, we show for the first time that a High-Plex amplicon-based massively parallel sequencing (MPS) system can be applied effectively to DNA prepared from dried blood spot archival specimens and, as such, can dramatically increase the scopes of both method and resource.