Mapping recurrent mosaic copy number variation in human neurons.

When somatic cells acquire complex karyotypes, they often are removed by the immune system. Mutant somatic cells that evade immune surveillance can lead to cancer. Neurons with complex karyotypes arise during neurotypical brain development, but neurons are almost never the origin of brain cancers. Instead, somatic mutations in neurons can bring about neurodevelopmental disorders, and contribute to the polygenic landscape of neuropsychiatric and neurodegenerative disease. A subset of human neurons harbors idiosyncratic copy number variants (CNVs, "CNV neurons"), but previous analyses of CNV neurons are limited by relatively small sample sizes. Here, we develop an allele-based validation approach, SCOVAL, to corroborate or reject read-depth based CNV calls in single human neurons. We apply this approach to 2,125 frontal cortical neurons from a neurotypical human brain. SCOVAL identifies 226 CNV neurons, which include a subclass of 65 CNV neurons with highly aberrant karyotypes containing whole or substantial losses on multiple chromosomes. Moreover, we find that CNV location appears to be nonrandom. Recurrent regions of neuronal genome rearrangement contain fewer, but longer, genes.

Genomic data resources of the Brain Somatic Mosaicism Network for neuropsychiatric diseases.

Somatic mosaicism is defined as an occurrence of two or more populations of cells having genomic sequences differing at given loci in an individual who is derived from a single zygote. It is a characteristic of multicellular organisms that plays a crucial role in normal development and disease. To study the nature and extent of somatic mosaicism in autism spectrum disorder, bipolar disorder, focal cortical dysplasia, schizophrenia, and Tourette syndrome, a multi-institutional consortium called the Brain Somatic Mosaicism Network (BSMN) was formed through the National Institute of Mental Health (NIMH). In addition to genomic data of affected and neurotypical brains, the BSMN also developed and validated a best practices somatic single nucleotide variant calling workflow through the analysis of reference brain tissue. These resources, which include >400 terabytes of data from 1087 subjects, are now available to the research community via the NIMH Data Archive (NDA) and are described here.

Mapping the Complex Genetic Landscape of Human Neurons.

When somatic cells acquire complex karyotypes, they are removed by the immune system. Mutant somatic cells that evade immune surveillance can lead to cancer. Neurons with complex karyotypes arise during neurotypical brain development, but neurons are almost never the origin of brain cancers. Instead, somatic mutations in neurons can bring about neurodevelopmental disorders, and contribute to the polygenic landscape of neuropsychiatric and neurodegenerative disease. A subset of human neurons harbors idiosyncratic copy number variants (CNVs, "CNV neurons"), but previous analyses of CNV neurons have been limited by relatively small sample sizes. Here, we developed an allele-based validation approach, SCOVAL, to corroborate or reject read-depth based CNV calls in single human neurons. We applied this approach to 2,125 frontal cortical neurons from a neurotypical human brain. This approach identified 226 CNV neurons, as well as a class of CNV neurons with complex karyotypes containing whole or substantial losses on multiple chromosomes. Moreover, we found that CNV location appears to be nonrandom. Recurrent regions of neuronal genome rearrangement contained fewer, but longer, genes.

Comprehensive identification of somatic nucleotide variants in human brain tissue.

BACKGROUND: Post-zygotic mutations incurred during DNA replication, DNA repair, and other cellular processes lead to somatic mosaicism. Somatic mosaicism is an established cause of various diseases, including cancers. However, detecting mosaic variants in DNA from non-cancerous somatic tissues poses significant challenges, particularly if the variants only are present in a small fraction of cells. RESULTS: Here, the Brain Somatic Mosaicism Network conducts a coordinated, multi-institutional study to examine the ability of existing methods to detect simulated somatic single-nucleotide variants (SNVs) in DNA mixing experiments, generate multiple replicates of whole-genome sequencing data from the dorsolateral prefrontal cortex, other brain regions, dura mater, and dural fibroblasts of a single neurotypical individual, devise strategies to discover somatic SNVs, and apply various approaches to validate somatic SNVs. These efforts lead to the identification of 43 bona fide somatic SNVs that range in variant allele fractions from ~ 0.005 to ~ 0.28. Guided by these results, we devise best practices for calling mosaic SNVs from 250× whole-genome sequencing data in the accessible portion of the human genome that achieve 90% specificity and sensitivity. Finally, we demonstrate that analysis of multiple bulk DNA samples from a single individual allows the reconstruction of early developmental cell lineage trees. CONCLUSIONS: This study provides a unified set of best practices to detect somatic SNVs in non-cancerous tissues. The data and methods are freely available to the scientific community and should serve as a guide to assess the contributions of somatic SNVs to neuropsychiatric diseases.

Identification and characterization of occult human-specific LINE-1 insertions using long-read sequencing technology.

Long Interspersed Element-1 (LINE-1) retrotransposition contributes to inter- and intra-individual genetic variation and occasionally can lead to human genetic disorders. Various strategies have been developed to identify human-specific LINE-1 (L1Hs) insertions from short-read whole genome sequencing (WGS) data; however, they have limitations in detecting insertions in complex repetitive genomic regions. Here, we developed a computational tool (PALMER) and used it to identify 203 non-reference L1Hs insertions in the NA12878 benchmark genome. Using PacBio long-read sequencing data, we identified L1Hs insertions that were absent in previous short-read studies (90/203). Approximately 81% (73/90) of the L1Hs insertions reside within endogenous LINE-1 sequences in the reference assembly and the analysis of unique breakpoint junction sequences revealed 63% (57/90) of these L1Hs insertions could be genotyped in 1000 Genomes Project sequences. Moreover, we observed that amplification biases encountered in single-cell WGS experiments led to a wide variation in L1Hs insertion detection rates between four individual NA12878 cells; under-amplification limited detection to 32% (65/203) of insertions, whereas over-amplification increased false positive calls. In sum, these data indicate that L1Hs insertions are often missed using standard short-read sequencing approaches and long-read sequencing approaches can significantly improve the detection of L1Hs insertions present in individual genomes.

Y-Chromosome Structural Diversity in the Bonobo and Chimpanzee Lineages.

The male-specific regions of primate Y-chromosomes (MSY) are enriched for multi-copy genes highly expressed in the testis. These genes are located in large repetitive sequences arranged as palindromes, inverted-, and tandem repeats termed amplicons. In humans, these genes have critical roles in male fertility and are essential for the production of sperm. The structure of human and chimpanzee amplicon sequences show remarkable difference relative to the remainder of the genome, a difference that may be the result of intense selective pressure on male fertility. Four subspecies of common chimpanzees have undergone extended periods of isolation and appear to be in the early process of subspeciation. A recent study found amplicons enriched for testis-expressed genes on the primate X-chromosome the target of hard selective sweeps, and male-fertility genes on the Y-chromosome may also be the targets of selection. However, little is understood about Y-chromosome amplicon diversity within and across chimpanzee populations. Here, we analyze nine common chimpanzee (representing three subspecies: Pan troglodytes schweinfurthii, Pan troglodytes ellioti, and Pan troglodytes verus) and two bonobo (Pan paniscus) male whole-genome sequences to assess Y ampliconic copy-number diversity across the Pan genus. We observe that the copy number of Y chromosome amplicons is variable among chimpanzees and bonobos, and identify several lineage-specific patterns, including variable copy number of azoospermia candidates RBMY and DAZ We detect recurrent switchpoints of copy-number change along the ampliconic tracts across chimpanzee populations, which may be the result of localized genome instability or selective forces.

Resolving complex structural genomic rearrangements using a randomized approach.

Complex chromosomal rearrangements are structural genomic alterations involving multiple instances of deletions, duplications, inversions, or translocations that co-occur either on the same chromosome or represent different overlapping events on homologous chromosomes. We present SVelter, an algorithm that identifies regions of the genome suspected to harbor a complex event and then resolves the structure by iteratively rearranging the local genome structure, in a randomized fashion, with each structure scored against characteristics of the observed sequencing data. SVelter is able to accurately reconstruct complex chromosomal rearrangements when compared to well-characterized genomes that have been deeply sequenced with both short and long reads.

The genomic landscape of polymorphic human nuclear mitochondrial insertions.

The transfer of mitochondrial genetic material into the nuclear genomes of eukaryotes is a well-established phenomenon that has been previously limited to the study of static reference genomes. The recent advancement of high throughput sequencing has enabled an expanded exploration into the diversity of polymorphic nuclear mitochondrial insertions (NumtS) within human populations. We have developed an approach to discover and genotype novel Numt insertions using whole genome, paired-end sequencing data. We have applied this method to a thousand individuals in 20 populations from the 1000 Genomes Project and other datasets and identified 141 new sites of Numt insertions, extending our current knowledge of existing NumtS by almost 20%. We find that recent Numt insertions are derived from throughout the mitochondrial genome, including the D-loop, and have integration biases that differ in some respects from previous studies on older, fixed NumtS in the reference genome. We determined the complete inserted sequence for a subset of these events and have identified a number of nearly full-length mitochondrial genome insertions into nuclear chromosomes. We further define their age and origin of insertion and present an analysis of their potential impact to ongoing studies of mitochondrial heteroplasmy and disease.

CHD7 mutations and CHARGE syndrome in semicircular canal dysplasia.

OBJECTIVE: To determine whether patients with semicircular canal dysplasia have mutations in CHD7. BACKGROUND: CHARGE syndrome is a nonrandom clustering of congenital anomalies, including ocular coloboma, heart defects, choanal atresia or stenosis, retarded growth and development, genital hypoplasia, and inner and outer ear anomalies including deafness. Semicircular canal dysplasia has been included as a major diagnostic criterion for CHARGE syndrome. Mutations in the gene CHD7 on chromosome 8q12.1 are a major cause of CHARGE syndrome, but the extent to which patients with semicircular canal dysplasia have CHD7 mutations is not fully understood. STUDY DESIGN: Cross-sectional analysis of CHD7 in 12 patients with semicircular canal dysplasia and variable clinical features of CHARGE syndrome. RESULTS: We identified 6 CHD7 mutations, 5 of which occurred in patients who fulfilled Verloes' diagnostic criteria for typical CHARGE syndrome, and three of which were previously unreported. Of the 3 remaining CHD7 mutation-positive patients, one had atypical CHARGE by diagnostic criteria. Four MRI records were available, which revealed 2 patients with cochlear nerve aplasia and 1 patient with Chiari 1 malformation. CONCLUSION: These data provide additional evidence that CHD7 mutations are a significant cause of semicircular canal atresia in children with full or partial CHARGE syndrome.

A human homeotic transformation resulting from mutations in PLCB4 and GNAI3 causes auriculocondylar syndrome.

Auriculocondylar syndrome (ACS) is a rare, autosomal-dominant craniofacial malformation syndrome characterized by variable micrognathia, temporomandibular joint ankylosis, cleft palate, and a characteristic "question-mark" ear malformation. Careful phenotypic characterization of severely affected probands in our cohort suggested the presence of a mandibular patterning defect resulting in a maxillary phenotype (i.e., homeotic transformation). We used exome sequencing of five probands and identified two novel (exclusive to the patient and/or family studied) missense mutations in PLCB4 and a shared mutation in GNAI3 in two unrelated probands. In confirmatory studies, three additional novel PLCB4 mutations were found in multigenerational ACS pedigrees. All mutations were confirmed by Sanger sequencing, were not present in more than 10,000 control chromosomes, and resulted in amino-acid substitutions located in highly conserved protein domains. Additionally, protein-structure modeling demonstrated that all ACS substitutions disrupt the catalytic sites of PLCB4 and GNAI3. We suggest that PLCB4 and GNAI3 are core signaling molecules of the endothelin-1-distal-less homeobox 5 and 6 (EDN1-DLX5/DLX6) pathway. Functional studies demonstrated a significant reduction in downstream DLX5 and DLX6 expression in ACS cases in assays using cultured osteoblasts from probands and controls. These results support the role of the previously implicated EDN1-DLX5/6 pathway in regulating mandibular specification in other species, which, when disrupted, results in a maxillary phenotype. This work defines the molecular basis of ACS as a homeotic transformation (mandible to maxilla) in humans.

Identification of p.A684V missense mutation in the WFS1 gene as a frequent cause of autosomal dominant optic atrophy and hearing impairment.

Optic atrophy (OA) and sensorineural hearing loss (SNHL) are key abnormalities in several syndromes, including the recessively inherited Wolfram syndrome, caused by mutations in WFS1. In contrast, the association of autosomal dominant OA and SNHL without other phenotypic abnormalities is rare, and almost exclusively attributed to mutations in the Optic Atrophy-1 gene (OPA1), most commonly the p.R445H mutation. We present eight probands and their families from the US, Sweden, and UK with OA and SNHL, whom we analyzed for mutations in OPA1 and WFS1. Among these families, we found three heterozygous missense mutations in WFS1 segregating with OA and SNHL: p.A684V (six families), and two novel mutations, p.G780S and p.D797Y, all involving evolutionarily conserved amino acids and absent from 298 control chromosomes. Importantly, none of these families harbored the OPA1 p.R445H mutation. No mitochondrial DNA deletions were detected in muscle from one p.A684V patient analyzed. Finally, wolframin p.A684V mutant ectopically expressed in HEK cells showed reduced protein levels compared to wild-type wolframin, strongly indicating that the mutation is disease-causing. Our data support OA and SNHL as a phenotype caused by dominant mutations in WFS1 in these additional eight families. Importantly, our data provide the first evidence that a single, recurrent mutation in WFS1, p.A684V, may be a common cause of ADOA and SNHL, similar to the role played by the p.R445H mutation in OPA1. Our findings suggest that patients who are heterozygous for WFS1 missense mutations should be carefully clinically examined for OA and other manifestations of Wolfram syndrome.

Autosomal dominant progressive sensorineural hearing loss due to a novel mutation in the KCNQ4 gene.

OBJECTIVE: To identify the genetic etiology in a family with autosomal dominant progressive sensorineural hearing loss. DESIGN: Prospective molecular genetic research study. SETTING: Academic genetic research laboratory. PARTICIPANTS: Seventeen members of a family with dominant progressive nonsyndromic sensorineural hearing loss: 9 affected, 6 unaffected, and 2 spouses. INTERVENTIONS: Clinical data from questionnaires, interviews, serial audiograms, and medical records; genetic data from genome-wide linkage analysis and candidate gene mutation analysis. MAIN OUTCOME MEASURES: Symptoms, age at onset, serial audiometric data, and the presence or absence of a deafness-associated mutation. RESULTS: Affected individuals in this family presented with autosomal dominant nonsyndromic high-frequency progressive sensorineural hearing loss, with age at onset ranging from 1 to 21 years. Genome-wide linkage analysis of single-nucleotide polymorphisms yielded evidence of linkage to an 18.9-Mb region on chromosome 1p34-p36, with a multipoint logarithm of odds score of 3.6. This interval contains a known deafness gene, KCNQ4, which underlies DNFA2 deafness. Sequencing of the 14 coding exons and intron-exon junctions of KCNQ4 revealed a novel heterozygous missense mutation, c.859G>C, p.Gly287Arg. The mutation disrupts the highly conserved GYG motif (glycine-tyrosine-glycine) of the phosphate-binding loop, hypothesized to be critical in maintaining pore structure and function. All 274 controls were negative for the mutation. CONCLUSIONS: Autosomal dominant high-frequency hearing loss is genetically heterogeneous, and linkage analysis is an efficient means of identifying the etiology in larger families. Deafness in this family is caused by a novel mutation in KCNQ4.

Congenital cataracts in two siblings with Wolfram syndrome.

BACKGROUND: Wolfram syndrome is characterized by optic atrophy, insulin dependent diabetes mellitus, diabetes insipidus and deafness. There are several other associated conditions reported in the literature, but congenital or early childhood cataracts are not among them. MATERIALS AND METHODS: Observational case series with confirmatory genetic analysis. RESULTS: A pair of siblings, followed over 17 years, who manifest congenital or early childhood cataracts, diabetes insipidus, diabetes mellitus, optic atrophy, and deafness. They are both compound heterozygotes for mutations (V415 deletion and A684V substitution) in the WFS1 gene. Their father has congenital sensorineural hearing loss and developed optic atrophy. He is heterozygous for A684V in WFS1. CONCLUSIONS: Wolfram syndrome should be in the differential diagnosis of genetic syndromes associated with congenital and early childhood cataracts. Here, we report on a mother who is a phenotypically normal carrier of an autosomal recessive Wolfram syndrome gene, and a father who has some of the findings of the syndrome and carries a single mutation that appears to be responsible for his hearing loss and optic atrophy. Their 2 children are compound heterozygotes and manifest the full Wolfram syndrome, in addition to cataracts.

Increased activity of Diaphanous homolog 3 (DIAPH3)/diaphanous causes hearing defects in humans with auditory neuropathy and in Drosophila.

Auditory neuropathy is a rare form of deafness characterized by an absent or abnormal auditory brainstem response with preservation of outer hair cell function. We have identified Diaphanous homolog 3 (DIAPH3) as the gene responsible for autosomal dominant nonsyndromic auditory neuropathy (AUNA1), which we previously mapped to chromosome 13q21-q24. Genotyping of additional family members narrowed the interval to an 11-Mb, 3.28-cM gene-poor region containing only four genes, including DIAPH3. DNA sequencing of DIAPH3 revealed a c.-172G>A, g. 48G>A mutation in a highly conserved region of the 5' UTR. The c.-172G>A mutation occurs within a GC box sequence element and was not found in 379 controls. Using genome-wide expression arrays and quantitative RT-PCR, we demonstrate a 2- to 3-fold overexpression of DIAPH3 mRNA in lymphoblastoid cell lines from affected individuals. Likewise, a significant increase (approximately 1.5-fold) in DIAPH3 protein was found by quantitative immunoblotting of lysates from lymphoblastoid cell lines derived from affected individuals in comparison with controls. In addition, the c.-172G>A mutation is sufficient to drive overexpression of a luciferase reporter. Finally, the expression of a constitutively active form of diaphanous protein in the auditory organ of Drosophila melanogaster recapitulates the phenotype of impaired response to sound. To date, only two genes, the otoferlin gene OTOF and the pejvakin gene PJVK, are known to underlie nonsyndromic auditory neuropathy. Genetic testing for DIAPH3 may be useful for individuals with recessive as well as dominant inheritance of nonsyndromic auditory neuropathy.

Otosclerosis or congenital stapes ankylosis? The diagnostic role of genetic analysis.

HYPOTHESIS: Molecular genetic testing is useful to differentiate otosclerosis from syndromic stapes ankylosis. BACKGROUND: Congenital stapes ankylosis is genetically heterogeneous. Mutations in the NOG gene are known to be associated with a variety of rare stapes ankylosis syndromes including stapes ankylosis with broad thumbs and toes, multiple synostoses syndrome, and proximal symphalangism. These syndromes have overlapping clinical features that may be unrecognized. METHODS: The proband was a 54-year-old woman diagnosed in childhood with bilateral maximal conductive hearing loss. Audiologic, medical, and surgical records were reviewed. Deoxyribonucleic acid (DNA) was obtained from peripheral lymphocytes. DNA sequencing was used to assay for mutations in the NOG gene. RESULTS: Clinical genetics evaluation was most consistent with proximal symphalangism, but features of multiple synostoses syndrome were identified as well. DNA sequencing revealed a heterozygous p.W205C mutation in the NOG gene, not found in 100 controls. CONCLUSION: Evaluation of the patient with stapes ankylosis should include a family history and specific inquiry into features associated with stapes ankylosis syndromes, such as bony anomalies of the spine, hands, and feet. However, a negative family history does not exclude the possibility of a syndrome. Many patients who are thought to have nonsyndromic otosclerosis actually have syndromes caused by mutations in the NOG gene. Identifying a syndrome has implications for surgical management and prognosis.