Structural diversity and African origin of the 17q21.31 inversion polymorphism

The 17q21.31 inversion polymorphism exists either as direct (H1) or inverted (H2) haplotypes with differential predispositions to disease and selection. We investigated its genetic diversity in 2,700 individuals, with an emphasis on African populations. We characterize eight structural haplotypes due to complex rearrangements that vary in size from 1.08-1.49 Mb and provide evidence for a 30-kb H1-H2 double recombination event. We show that recurrent partial duplications of the KANSL1 gene have occurred on both the H1 and H2 haplotypes and have risen to high frequency in European populations. We identify a likely ancestral H2 haplotype (H2') lacking these duplications that is enriched among African hunter-gatherer groups yet essentially absent from West African populations. Whereas H1 and H2 segmental duplications arose independently and before human migration out of Africa, they have reached high frequencies recently among Europeans, either because of extraordinary genetic drift or selective sweeps.     

A common inversion under selection in Europeans

A refined physical map of chromosome 17q21.31 uncovered a 900-kb inversion polymorphism. Chromosomes with the inverted segment in different orientations represent two distinct lineages, H1 and H2, that have diverged for as much as 3 million years and show no evidence of having recombined. The H2 lineage is rare in Africans, almost absent in East Asians but found at a frequency of 20% in Europeans, in whom the haplotype structure is indicative of a history of positive selection. Here we show that the H2 lineage is undergoing positive selection in the Icelandic population, such that carrier females have more children and have higher recombination rates than noncarriers.     

Mutations in KANSL1 cause the 17q21.31 microdeletion syndrome phenotype

The chromosome 17q21.31 deletion syndrome is a genomic disorder characterized by highly distinctive facial features, moderate-to-severe intellectual disability, hypotonia and friendly behavior. Here, we show that de novo loss-of-function mutations in KANSL1 (also called KIAA1267) cause a full del(17q21.31) phenotype in two unrelated individuals that lack deletion at 17q21.31. These findings indicate that 17q21.31 deletion syndrome is a monogenic disorder caused by haploinsufficiency of KANSL1.     

Mutations in the chromatin modifier gene KANSL1 cause the 17q21.31 microdeletion syndrome

We show that haploinsufficiency of KANSL1 is sufficient to cause the 17q21.31 microdeletion syndrome, a multisystem disorder characterized by intellectual disability, hypotonia and distinctive facial features. The KANSL1 protein is an evolutionarily conserved regulator of the chromatin modifier KAT8, which influences gene expression through histone H4 lysine 16 (H4K16) acetylation. RNA sequencing studies in cell lines derived from affected individuals and the presence of learning deficits in Drosophila melanogaster mutants suggest a role for KANSL1 in neuronal processes.     

Discovery of previously unidentified genomic disorders from the duplication architecture of the human genome

Genomic disorders are characterized by the presence of flanking segmental duplications that predispose these regions to recurrent rearrangement. Based on the duplication architecture of the genome, we investigated 130 regions that we hypothesized as candidates for previously undescribed genomic disorders. We tested 290 individuals with mental retardation by BAC array comparative genomic hybridization and identified 16 pathogenic rearrangements, including de novo microdeletions of 17q21.31 found in four individuals. Using oligonucleotide arrays, we refined the breakpoints of this microdeletion, defining a 478-kb critical region containing six genes that were deleted in all four individuals. We mapped the breakpoints of this deletion and of four other pathogenic rearrangements in 1q21.1, 15q13, 15q24 and 17q12 to flanking segmental duplications, suggesting that these are also sites of recurrent rearrangement. In common with the 17q21.31 deletion, these breakpoint regions are sites of copy number polymorphism in controls, indicating that these may be inherently unstable genomic regions.     

A new chromosome 17q21.31 microdeletion syndrome associated with a common inversion polymorphism

Submicroscopic genomic copy number changes have been identified only recently as an important cause of mental retardation. We describe the detection of three interstitial, overlapping 17q21.31 microdeletions in a cohort of 1,200 mentally retarded individuals associated with a clearly recognizable clinical phenotype of mental retardation, hypotonia and a characteristic face. The deletions encompass the MAPT and CRHR1 genes and are associated with a common inversion polymorphism.     

Microdeletion encompassing MAPT at chromosome 17q21.3 is associated with developmental delay and learning disability

Recently, the application of array-based comparative genomic hybridization (array CGH) has improved rates of detection of chromosomal imbalances in individuals with mental retardation and dysmorphic features. Here, we describe three individuals with learning disability and a heterozygous deletion at chromosome 17q21.3, detected in each case by array CGH. FISH analysis demonstrated that the deletions occurred as de novo events in each individual and were between 500 kb and 650 kb in size. A recently described 900-kb inversion that suppresses recombination between ancestral H1 and H2 haplotypes encompasses the deletion. We show that, in each trio, the parent of origin of the deleted chromosome 17 carries at least one H2 chromosome. This region of 17q21.3 shows complex genomic architecture with well-described low-copy repeats (LCRs). The orientation of LCRs flanking the deleted segment in inversion heterozygotes is likely to facilitate the generation of this microdeletion by means of non-allelic homologous recombination.     

A copy number variation morbidity map of developmental delay

To understand the genetic heterogeneity underlying developmental delay, we compared copy number variants (CNVs) in 15,767 children with intellectual disability and various congenital defects (cases) to CNVs in 8,329 unaffected adult controls. We estimate that ～14.2% of disease in these children is caused by CNVs >400 kb. We observed a greater enrichment of CNVs in individuals with craniofacial anomalies and cardiovascular defects compared to those with epilepsy or autism. We identified 59 pathogenic CNVs, including 14 new or previously weakly supported candidates, refined the critical interval for several genomic disorders, such as the 17q21.31 microdeletion syndrome, and identified 940 candidate dosage-sensitive genes. We also developed methods to opportunistically discover small, disruptive CNVs within the large and growing diagnostic array datasets. This evolving CNV morbidity map, combined with exome and genome sequencing, will be critical for deciphering the genetic basis of developmental delay, intellectual disability and autism spectrum disorders.     

Genome-wide association study identifies six new loci influencing pulse pressure and mean arterial pressure

Numerous genetic loci have been associated with systolic blood pressure (SBP) and diastolic blood pressure (DBP) in Europeans. We now report genome-wide association studies of pulse pressure (PP) and mean arterial pressure (MAP). In discovery (N = 74,064) and follow-up studies (N = 48,607), we identified at genome-wide significance (P = 2.7 × 10(-8) to P = 2.3 × 10(-13)) four new PP loci (at 4q12 near CHIC2, 7q22.3 near PIK3CG, 8q24.12 in NOV and 11q24.3 near ADAMTS8), two new MAP loci (3p21.31 in MAP4 and 10q25.3 near ADRB1) and one locus associated with both of these traits (2q24.3 near FIGN) that has also recently been associated with SBP in east Asians. For three of the new PP loci, the estimated effect for SBP was opposite of that for DBP, in contrast to the majority of common SBP- and DBP-associated variants, which show concordant effects on both traits. These findings suggest new genetic pathways underlying blood pressure variation, some of which may differentially influence SBP and DBP.     

Inverted genomic segments and complex triplication rearrangements are mediated by inverted repeats in the human genome

We identified complex genomic rearrangements consisting of intermixed duplications and triplications of genomic segments at the MECP2 and PLP1 loci. These complex rearrangements were characterized by a triplicated segment embedded within a duplication in 11 unrelated subjects. Notably, only two breakpoint junctions were generated during each rearrangement formation. All the complex rearrangement products share a common genomic organization, duplication-inverted triplication-duplication (DUP-TRP/INV-DUP), in which the triplicated segment is inverted and located between directly oriented duplicated genomic segments. We provide evidence that the DUP-TRP/INV-DUP structures are mediated by inverted repeats that can be separated by >300 kb, a genomic architecture that apparently leads to susceptibility to such complex rearrangements. A similar inverted repeat-mediated mechanism may underlie structural variation in many other regions of the human genome. We propose a mechanism that involves both homology-driven events, via inverted repeats, and microhomologous or nonhomologous events.     

Limb-girdle muscular dystrophy type 2G is caused by mutations in the gene encoding the sarcomeric protein telethonin

Autosomal recessive limb-girdle muscular dystrophies (AR LGMDs) are a genetically heterogeneous group of disorders that affect mainly the proximal musculature. There are eight genetically distinct forms of AR LGMD, LGMD 2A-H (refs 2-10), and the genetic lesions underlying these forms, except for LGMD 2G and 2H, have been identified. LGMD 2A and LGMD 2B are caused by mutations in the genes encoding calpain 3 (ref. 11) and dysferlin, respectively, and are usually associated with a mild phenotype. Mutations in the genes encoding gamma-(ref. 14), alpha-(ref. 5), beta-(refs 6,7) and delta (ref. 15)-sarcoglycans are responsible for LGMD 2C to 2F, respectively. Sarcoglycans, together with sarcospan, dystroglycans, syntrophins and dystrobrevin, constitute the dystrophin-glycoprotein complex (DGC). Patients with LGMD 2C-F predominantly have a severe clinical course. The LGMD 2G locus maps to a 3-cM interval in 17q11-12 in two Brazilian families with a relatively mild form of AR LGMD (ref. 9). To positionally clone the LGMD 2G gene, we constructed a physical map of the 17q11-12 region and refined its localization to an interval of 1.2 Mb. The gene encoding telethonin, a sarcomeric protein, lies within this candidate region. We have found that mutations in the telethonin gene cause LGMD 2G, identifying a new molecular mechanism for AR LGMD.     

Identification of large ancient duplications associated with human gene deserts

We identified 15 regions of >1 Mb in the human genome composed of large ancient local duplications corresponding to gene deserts. We detected these intrachromosomal duplications in mouse and dog but not in chicken; they present as patches of similarity as low as 60%. These findings suggest that some human gene deserts originated from duplications of segments lacking genes in a mammalian common ancestor.     

Segregation at three loci explains familial and population risk in Hirschsprung disease

Hirschsprung disease (HSCR), the most common hereditary cause of intestinal obstruction, shows considerable variation and complex inheritance. Coding sequence mutations in RET, GDNF, EDNRB, EDN3 and SOX10 lead to long-segment (L-HSCR) and syndromic HSCR but fail to explain the transmission of the much more common short-segment form (S-HSCR). We conducted a genome scan in families with S-HSCR and identified susceptibility loci at 3p21, 10q11 and 19q12 that seem to be necessary and sufficient to explain recurrence risk and population incidence. The gene at 10q11 is probably RET, supporting its crucial role in all forms of HSCR; however, coding sequence mutations are present in only 40% of linked families, suggesting the importance of noncoding variation. Here we show oligogenic inheritance of S-HSCR, the 3p21 and 19q12 loci as RET-dependent modifiers, and a parent-of-origin effect at RET. This study demonstrates by a complete genetic dissection why the inheritance pattern of S-HSCR is nonmendelian.     

A 1.5 million-base pair inversion polymorphism in families with Williams-Beuren syndrome.

Williams-Beuren syndrome (WBS) is most often caused by hemizygous deletion of a 1.5-Mb interval encompassing at least 17 genes at 7q11.23 (refs. 1,2). As with many other haploinsufficiency diseases, the mechanism underlying the WBS deletion is thought to be unequal meiotic recombination, probably mediated by the highly homologous DNA that flanks the commonly deleted region. Here, we report the use of interphase fluorescence in situ hybridization (FISH) and pulsed-field gel electrophoresis (PFGE) to identify a genomic polymorphism in families with WBS, consisting of an inversion of the WBS region. We have observed that the inversion is hemizygous in 3 of 11 (27%) atypical affected individuals who show a subset of the WBS phenotypic spectrum but do not carry the typical WBS microdeletion. Two of these individuals also have a parent who carries the inversion. In addition, in 4 of 12 (33%) families with a proband carrying the WBS deletion, we observed the inversion exclusively in the parent transmitting the disease-related chromosome. These results suggest the presence of a newly identified genomic variant within the population that may be associated with the disease. It may result in predisposition to primarily WBS-causing microdeletions, but may also cause translocations and inversions.     

Positional cloning of a novel gene influencing asthma from chromosome 2q14

Asthma is a common disease in children and young adults. Four separate reports have linked asthma and related phenotypes to an ill-defined interval between 2q14 and 2q32 (refs. 1-4), and two mouse genome screens have linked bronchial hyper-responsiveness to the region homologous to 2q14 (refs. 5,6). We found and replicated association between asthma and the D2S308 microsatellite, 800 kb distal to the IL1 cluster on 2q14. We sequenced the surrounding region and constructed a comprehensive, high-density, single-nucleotide polymorphism (SNP) linkage disequilibrium (LD) map. SNP association was limited to the initial exons of a solitary gene of 3.6 kb (DPP10), which extends over 1 Mb of genomic DNA. DPP10 encodes a homolog of dipeptidyl peptidases (DPPs) that cleave terminal dipeptides from cytokines and chemokines, and it presents a potential new target for asthma therapy.     

A common CFH haplotype, with deletion of CFHR1 and CFHR3, is associated with lower risk of age-related macular degeneration

Age-related macular degeneration (AMD; OMIM #603075) is the most frequent cause of visual impairment in the elderly population, with severe disease affecting nearly 10% of individuals of European descent over the age of 75 years. It is a complex disease in which genetic and environmental factors contribute to susceptibility. Complement factor H (CFH) has recently been identified as a major AMD susceptibility gene, and the Y402H polymorphism has been proposed as the likely causative factor. We genotyped polymorphisms spanning the cluster of CFH and five CFH-related genes on chromosome 1q23 in 173 individuals with severe neovascular AMD and 170 elderly controls with no signs of AMD. Detailed analysis showed a common haplotype associated with decreased risk of AMD that was present on 20% of chromosomes of controls and 8% of chromosomes of individuals with AMD. We found that this haplotype carried a deletion of CFHR1 and CFHR3, and the proteins encoded by these genes were absent in serum of homozygotes. The protective effect of the deletion haplotype cannot be attributed to linkage disequilibrium with Y402H and was replicated in an independent sample.     

Complex SNP-related sequence variation in segmental genome duplications

There is uncertainty about the true nature of predicted single-nucleotide polymorphisms (SNPs) in segmental duplications (duplicons) and whether these markers genuinely exist at increased density as indicated in public databases. We explored these issues by genotyping 157 predicted SNPs in duplicons and control regions in normal diploid genomes and fully homozygous complete hydatidiform moles. Our data identified many true SNPs in duplicon regions and few paralogous sequence variants. Twenty-eight percent of the polymorphic duplicon sequences we tested involved multisite variation, a new type of polymorphism representing the sum of the signals from many individual duplicon copies that vary in sequence content due to duplication, deletion or gene conversion. Multisite variations can masquerade as normal SNPs when genotyped. Given that duplicons comprise at least 5% of the genome and many are yet to be annotated in the genome draft, effective strategies to identify multisite variation must be established and deployed.     

Discovery and genotyping of genome structural polymorphism by sequencing on a population scale

Accurate and complete analysis of genome variation in large populations will be required to understand the role of genome variation in complex disease. We present an analytical framework for characterizing genome deletion polymorphism in populations using sequence data that are distributed across hundreds or thousands of genomes. Our approach uses population-level concepts to reinterpret the technical features of sequence data that often reflect structural variation. In the 1000 Genomes Project pilot, this approach identified deletion polymorphism across 168 genomes (sequenced at 4 × average coverage) with sensitivity and specificity unmatched by other algorithms. We also describe a way to determine the allelic state or genotype of each deletion polymorphism in each genome; the 1000 Genomes Project used this approach to type 13,826 deletion polymorphisms (48-995,664 bp) at high accuracy in populations. These methods offer a way to relate genome structural polymorphism to complex disease in populations.