Whole-Exome Sequencing Efficiently Detects Rare Mutations in Autosomal Recessive Nonsyndromic Hearing Loss

Identification of the pathogenic mutations underlying autosomal recessive nonsyndromic hearing loss (ARNSHL) is difficult, since causative mutations in 39 different genes have so far been reported. After excluding mutations in the most common ARNSHL gene, GJB2, via Sanger sequencing, we performed whole-exome sequencing (WES) in 30 individuals from 20 unrelated multiplex consanguineous families with ARNSHL. Agilent SureSelect Human All Exon 50 Mb kits and an Illumina Hiseq2000 instrument were used. An average of 93%, 84% and 73% of bases were covered to 1X, 10X and 20X within the ARNSHL-related coding RefSeq exons, respectively. Uncovered regions with WES included those that are not targeted by the exome capture kit and regions with high GC content. Twelve homozygous mutations in known deafness genes, of which eight are novel, were identified in 12 families: MYO15A-p.Q1425X, -p.S1481P, -p.A1551D; LOXHD1-p.R1494X, -p.E955X; GIPC3-p.H170N; ILDR1-p.Q274X; MYO7A-p.G2163S; TECTA-p.Y1737C; TMC1-p.S530X; TMPRSS3-p.F13Lfs*10; TRIOBP-p.R785Sfs*50. Each mutation was within a homozygous run documented via WES. Sanger sequencing confirmed co-segregation of the mutation with deafness in each family. Four rare heterozygous variants, predicted to be pathogenic, in known deafness genes were detected in 12 families where homozygous causative variants were already identified. Six heterozygous variants that had similar characteristics to those abovementioned variants were present in 15 ethnically-matched individuals with normal hearing. Our results show that rare causative mutations in known ARNSHL genes can be reliably identified via WES. The excess of heterozygous variants should be considered during search for causative mutations in ARNSHL genes, especially in small-sized families.     

Mutations in C12orf57 Cause a Syndromic Form of Colobomatous Microphthalmia

Microphthalmia is an important developmental eye disorder. Although mutations in several genes have been linked to this condition, they only account for a minority of cases. We performed autozygome analysis and exome sequencing on a multiplex consanguineous family in which colobomatous microphthalmia is associated with profound global developmental delay, intractable seizures, and corpus callosum abnormalities, and we identified a homozygous truncating mutation in C12orf57 [c.1A>G; p.Met1?]. In a simplex case with a similar phenotype, we identified compound heterozygosity for the same mutation and another missense mutation [c.152T>A; p.Leu51Gln]. Little is known about C12orf57 but we show that it is expressed in several mouse tissues, including the eye and brain. Our data strongly implicate mutations in C12orf57 in the pathogenesis of a clinically distinct autosomal-recessive syndromic form of colobomatous microphthalmia.     

Whole-Exome Sequencing Identifies LRIT3 Mutations as a Cause of Autosomal-Recessive Complete Congenital Stationary Night Blindness

Congenital stationary night blindness (CSNB) is a clinically and genetically heterogeneous retinal disorder. Two forms can be distinguished clinically: complete CSNB (cCSNB) and incomplete CSNB. Individuals with cCSNB have visual impairment under low-light conditions and show a characteristic electroretinogram (ERG). The b-wave amplitude is severely reduced in the dark-adapted state of the ERG, representing abnormal function of ON bipolar cells. Furthermore, individuals with cCSNB can show other ocular features such as nystagmus, myopia, and strabismus and can have reduced visual acuity and abnormalities of the cone ERG waveform. The mode of inheritance of this form can be X-linked or autosomal recessive, and the dysfunction of four genes (NYX, GRM6, TRPM1, and GPR179) has been described so far. Whole-exome sequencing in one simplex cCSNB case lacking mutations in the known genes led to the identification of a missense mutation (c.983G>A [p.Cys328Tyr]) and a nonsense mutation (c.1318C>T [p.Arg440^∗]) in LRIT3, encoding leucine-rich-repeat (LRR), immunoglobulin-like, and transmembrane-domain 3 (LRIT3). Subsequent Sanger sequencing of 89 individuals with CSNB identified another cCSNB case harboring a nonsense mutation (c.1151C>G [p.Ser384^∗]) and a deletion predicted to lead to a premature stop codon (c.1538_1539del [p.Ser513Cysfs^∗59]) in the same gene. Human LRIT3 antibody staining revealed in the outer plexiform layer of the human retina a punctate-labeling pattern resembling the dendritic tips of bipolar cells; similar patterns have been observed for other proteins implicated in cCSNB. The exact role of this LRR protein in cCSNB remains to be elucidated.     

Targeted Capture and Next-Generation Sequencing Identifies C9orf75, Encoding Taperin, as the Mutated Gene in Nonsyndromic Deafness DFNB79

Targeted genome capture combined with next-generation sequencing was used to analyze 2.9 Mb of the DFNB79 interval on chromosome 9q34.3, which includes 108 candidate genes. Genomic DNA from an affected member of a consanguineous family segregating recessive, nonsyndromic hearing loss was used to make a library of fragments covering the DFNB79 linkage interval defined by genetic analyses of four pedigrees. Homozygosity for eight previously unreported variants in transcribed sequences was detected by evaluating a library of 402,554 sequencing reads and was later confirmed by Sanger sequencing. Of these variants, six were determined to be polymorphisms in the Pakistani population, and one was in a noncoding gene that was subsequently excluded genetically from the DFNB79 linkage interval. The remaining variant was a nonsense mutation in a predicted gene, C9orf75, renamed TPRN. Evaluation of the other three DFNB79-linked families identified three additional frameshift mutations, for a total of four truncating alleles of this gene. Although TPRN is expressed in many tissues, immunolocalization of the protein product in the mouse cochlea shows prominent expression in the taper region of hair cell stereocilia. Consequently, we named the protein taperin.     

Whole Exome Sequencing Identifies Novel Recurrently Mutated Genes in Patients with Splenic Marginal Zone Lymphoma

The pathogenesis of splenic marginal zone lymphoma (SMZL) remains largely unknown. Recent high-throughput sequencing studies have identified recurrent mutations in key pathways, most notably NOTCH2 mutations in >25% of patients. These studies are based on small, heterogeneous discovery cohorts, and therefore only captured a fraction of the lesions present in the SMZL genome. To identify further novel pathogenic mutations within related biochemical pathways, we applied whole exome sequencing (WES) and copy number (CN) analysis to a biologically and clinically homogeneous cohort of seven SMZL patients with 7q abnormalities and IGHV1-2*04 gene usage. We identified 173 somatic non-silent variants, affecting 160 distinct genes. In additional to providing independent validation of the presence of mutation in several previously reported genes (NOTCH2, TNFAIP3, MAP3K14, MLL2 and SPEN), our study defined eight additional recurrently mutated genes in SMZL; these genes are CREBBP, CBFA2T3, AMOTL1, FAT4, FBXO11, PLA2G4D, TRRAP and USH2A. By integrating our WES and CN data we identified three mutated putative candidate genes targeted by 7q deletions (CUL1, EZH2 and FLNC), with FLNC positioned within the well-characterized 7q minimally deleted region. Taken together, this work expands the reported directory of recurrently mutated cancer genes in this disease, thereby expanding our understanding of SMZL pathogenesis. Ultimately, this work will help to establish a stratified approach to care including the possibility of targeted therapy.     

Whole-Exome Sequencing Identifies a Novel Genotype-Phenotype Correlation in the Entactin Domain of the Known Deafness Gene TECTA

Postlingual progressive hearing loss, affecting primarily the high frequencies, is the clinical finding in most cases of autosomal dominant nonsyndromic hearing loss (ADNSHL). The molecular genetic etiology of ADNSHL is extremely heterogeneous. We applied whole-exome sequencing to reveal the genetic etiology of high-frequency hearing loss in a mid-sized Korean family without any prior linkage data. Whole-exome sequencing of four family members (two affected and two unaffected), together with our filtering strategy based on comprehensive bioinformatics analyses, identified 21 potential pathogenic candidates. Sanger validation of an additional five family members excluded 20 variants, leaving only one novel variant, TECTA c.710C>T (p.T237I), as the strongest candidate. This variant resides in the entactin (ENT) domain and co-segregated perfectly with non-progressive high-frequency hearing loss in the family. It was absent among 700 ethnically matched control chromosomes, and the T237 residue is conserved among species, which supports its pathogenicity. Interestingly, this finding contrasted with a previously proposed genotype-phenotype correlation in which variants of the ENT domain of TECTA were associated with mid-frequency hearing loss. Based upon what we observed, we propose a novel “genotype to phenotype” correlation in the ENT domain of TECTA. Our results shed light on another important application of whole-exome sequencing: the establishment of a novel genotype-phenotype in the molecular genetic diagnosis of autosomal dominant hearing loss.     

Whole Exome Sequencing Identifies Recessive PKHD1 Mutations in a Chinese Twin Family with Caroli Disease

Background

Mutations in PKHD1 cause autosomal recessive Caroli disease, which is a rare congenital disorder involving cystic dilatation of the intrahepatic bile ducts. However, the mutational spectrum of PKHD1 and the phenotype-genotype correlations have not yet been fully established.

Methods

Whole exome sequencing (WES) was performed on one twin sample with Caroli disease from a Chinese family from Shandong province. Routine Sanger sequencing was used to validate the WES and to carry out segregation studies. We also described the PKHD1 mutation associated with the genotype-phenotype of this twin.

Results

A combination of WES and Sanger sequencing revealed the genetic defect to be a novel compound heterozygous genotype in PKHD1, including the missense mutation c.2507 T>C, predicted to cause a valine to alanine substitution at codon 836 (c.2507T>C, p.Val836Ala), and the nonsense mutation c.2341C>T, which is predicted to result in an arginine to stop codon at codon 781 (c.2341C>T, p.Arg781*). This compound heterozygous genotype co-segregates with the Caroli disease-affected pedigree members, but is absent in 200 normal chromosomes.

Conclusions

Our findings indicate exome sequencing can be useful in the diagnosis of Caroli disease patients and associate a compound heterozygous genotype in PKHD1 with Caroli disease, which further increases our understanding of the mutation spectrum of PKHD1 in association with Caroli disease.     

Both Rare and De Novo Copy Number Variants Are Prevalent in Agenesis of the Corpus Callosum but Not in Cerebellar Hypoplasia or Polymicrogyria

Agenesis of the corpus callosum (ACC), cerebellar hypoplasia (CBLH), and polymicrogyria (PMG) are severe congenital brain malformations with largely undiscovered causes. We conducted a large-scale chromosomal copy number variation (CNV) discovery effort in 255 ACC, 220 CBLH, and 147 PMG patients, and 2,349 controls. Compared to controls, significantly more ACC, but unexpectedly not CBLH or PMG patients, had rare genic CNVs over one megabase (p = 1.48×10⁻³; odds ratio [OR] = 3.19; 95% confidence interval [CI] = 1.89–5.39). Rare genic CNVs were those that impacted at least one gene in less than 1% of the combined population of patients and controls. Compared to controls, significantly more ACC but not CBLH or PMG patients had rare CNVs impacting over 20 genes (p = 0.01; OR = 2.95; 95% CI = 1.69–5.18). Independent qPCR confirmation showed that 9.4% of ACC patients had de novo CNVs. These, in comparison to inherited CNVs, preferentially overlapped de novo CNVs previously observed in patients with autism spectrum disorders (p = 3.06×10⁻⁴; OR = 7.55; 95% CI = 2.40–23.72). Interestingly, numerous reports have shown a reduced corpus callosum area in autistic patients, and diminished social and executive function in many ACC patients. We also confirmed and refined previously known CNVs, including significantly narrowing the 8p23.1-p11.1 duplication present in 2% of our current ACC cohort. We found six novel CNVs, each in a single patient, that are likely deleterious: deletions of 1p31.3-p31.1, 1q31.2-q31.3, 5q23.1, and 15q11.2-q13.1; and duplications of 2q11.2-q13 and 11p14.3-p14.2. One ACC patient with microcephaly had a paternally inherited deletion of 16p13.11 that included NDE1. Exome sequencing identified a recessive maternally inherited nonsense mutation in the non-deleted allele of NDE1, revealing the complexity of ACC genetics. This is the first systematic study of CNVs in congenital brain malformations, and shows a much higher prevalence of large gene-rich CNVs in ACC than in CBLH and PMG.     

Morphometric Changes of the Corpus Callosum in Congenital Blindness

We examined the effects of visual deprivation at birth on the development of the corpus callosum in a large group of congenitally blind individuals. We acquired high-resolution T1-weighted MRI scans in 28 congenitally blind and 28 normal sighted subjects matched for age and gender. There was no overall group effect of visual deprivation on the total surface area of the corpus callosum. However, subdividing the corpus callosum into five subdivisions revealed significant regional changes in its three most posterior parts. Compared to the sighted controls, congenitally blind individuals showed a 12% reduction in the splenium, and a 20% increase in the isthmus and the posterior part of the body. A shape analysis further revealed that the bending angle of the corpus callosum was more convex in congenitally blind compared to the sighted control subjects. The observed morphometric changes in the corpus callosum are in line with the well-described cross-modal functional and structural neuroplastic changes in congenital blindness.     

In-Cell Intrabody Selection from a Diverse Human Library Identifies C12orf4 Protein as a New Player in Rodent Mast Cell Degranulation

The high specificity of antibodies for their antigen allows a fine discrimination of target conformations and post-translational modifications, making antibodies the first choice tool to interrogate the proteome. We describe here an approach based on a large-scale intracellular expression and selection of antibody fragments in eukaryotic cells, so-called intrabodies, and the subsequent identification of their natural target within living cell. Starting from a phenotypic trait, this integrated system allows the identification of new therapeutic targets together with their companion inhibitory intrabody. We applied this system in a model of allergy and inflammation. We first cloned a large and highly diverse intrabody library both in a plasmid and a retroviral eukaryotic expression vector. After transfection in the RBL-2H3 rat basophilic leukemia cell line, we performed seven rounds of selection to isolate cells displaying a defect in FcεRI-induced degranulation. We used high throughput sequencing to identify intrabody sequences enriched during the course of selection. Only one intrabody was common to both plasmid and retroviral selections, and was used to capture and identify its target from cell extracts. Mass spectrometry analysis identified protein RGD1311164 (C12orf4), with no previously described function. Our data demonstrate that RGD1311164 is a cytoplasmic protein implicated in the early signaling events following FcεRI-induced cell activation. This work illustrates the strength of the intrabody-based in-cell selection, which allowed the identification of a new player in mast cell activation together with its specific inhibitor intrabody.     

Detecting Identity by Descent and Homozygosity Mapping in Whole-Exome Sequencing Data

The detection of genetic segments of Identical by Descent (IBD) in Genome-Wide Association Studies has proven successful in pinpointing genetic relatedness between reportedly unrelated individuals and leveraging such regions to shortlist candidate genes. These techniques depend on high-density genotyping arrays and their effectiveness in diverse sequence data is largely unknown. Due to decreasing costs and increasing effectiveness of high throughput techniques for whole-exome sequencing, an influx of exome sequencing data has become available. Studies using exomes and IBD-detection methods within known pedigrees have shown that IBD can be useful in finding hidden genetic candidates where known relatives are available. We set out to examine the viability of using IBD-detection in whole exome sequencing data in population-wide studies. In doing so, we extend GERMLINE, a method to detect IBD from exome sequencing data by finding small slices of matching alleles between pairs of individuals and extending them into full IBD segments. This algorithm allows for efficient population-wide detection in dense data. We apply this algorithm to a cohort of Crohn''s Disease cases where whole-exome and GWAS array data is available. We confirm that GWAS-based detected segments are highly accurate and predictive of underlying shared variation. Where segments inferred from GWAS are expected to be of high accuracy, we compare exome-based detection accuracy of multiple detection strategies. We find detection accuracy to be prohibitively low in all assessments, both in terms of segment sensitivity and specificity. Even after isolating relatively long segments beyond 10cM, exome-based detection continued to offer poor specificity/sensitivity tradeoffs. We hypothesize that the variable coverage and platform biases of exome capture account for this decreased accuracy and look toward whole genome sequencing data as a higher quality source for detecting population-wide IBD.     

Genetics of Microstructure of the Corpus Callosum in Older Adults

The current study sought to examine the relative influence of genetic and environmental factors on corpus callosum (CC) microstructure in a community sample of older adult twins. Analyses were undertaken in 284 healthy older twins (66% female; 79 MZ and 63 DZ pairs) from the Older Australian Twins Study. The average age of the sample was 69.82 (SD = 4.76) years. Brain imaging scans were collected and DTI measures were estimated for the whole CC as well as its five subregions. Parcellation of the CC was performed using Analyze. In addition, white matter lesion (WMLs) burden was estimated. Heritability and genetic correlation analyses were undertaken using the SOLAR software package. Age, sex, scanner, handedness and blood pressure were considered as covariates. Heritability (h²) analysis for the DTI metrics of whole CC, indicated significant h² for fractional anisotropy (FA) (h² = 0.56; p = 2.89×10⁻¹⁰), mean diffusivity (MD) (h² = 0.52; p = 0.30×10⁻⁶), radial diffusivity (RD) (h² = 0.49; p = 0.2×10⁻⁶) and axial diffusivity (AD) (h² = 0.37; p = 8.15×10⁻⁵). We also performed bivariate genetic correlation analyses between (i) whole CC DTI measures and (ii) whole CC DTI measures with total brain WML burden. Across the DTI measures for the whole CC, MD and RD shared 84% of the common genetic variance, followed by MD- AD (77%), FA - RD (52%), RD - AD (37%) and FA – MD (11%). For total WMLs, significant genetic correlations indicated that there was 19% shared common genetic variance with whole CC MD, followed by CC RD (17%), CC AD (16%) and CC FA (5%). Our findings suggest that the CC microstructure is under moderate genetic control. There was also evidence of shared genetic factors between the CC DTI measures. In contrast, there was less shared genetic variance between WMLs and the CC DTI metrics, suggesting fewer common genetic variants.     

Exome Sequencing Identifies a Founder Frameshift Mutation in an Alternative Exon of USH1C as the Cause of Autosomal Recessive Retinitis Pigmentosa with Late-Onset Hearing Loss

We used a combined approach of homozygosity mapping and whole exome sequencing (WES) to search for the genetic cause of autosomal recessive retinitis pigmentosa (arRP) in families of Yemenite Jewish origin. Homozygosity mapping of two arRP Yemenite Jewish families revealed a few homozygous regions. A subsequent WES analysis of the two index cases revealed a shared homozygous novel nucleotide deletion (c.1220delG) leading to a frameshift (p.Gly407Glufs*56) in an alternative exon (#15) of USH1C. Screening of additional Yemenite Jewish patients revealed a total of 16 homozygous RP patients (with a carrier frequency of 0.008 in controls). Funduscopic and electroretinography findings were within the spectrum of typical RP. While other USH1C mutations usually cause Usher type I (including RP, vestibular dysfunction and congenital deafness), audiometric screening of 10 patients who are homozygous for c.1220delG revealed that patients under 40 years of age had normal hearing while older patients showed mild to severe high tone sensorineural hearing loss. This is the first report of a mutation in a known USH1 gene that causes late onset rather than congenital sensorineural hearing loss. The c.1220delG mutation of USH1C accounts for 23% of RP among Yemenite Jewish patients in our cohort.     

Abnormality of the Corpus Callosum in Coalmine Gas Explosion-Related Posttraumatic Stress Disorder

Abnormal corpus callosum (CC) has been reported in childhood trauma-related posttraumatic stress disorder (PTSD); however, the nature of white matter (WM) integrity alterations in the CC of young adult-onset PTSD patients is unknown. In this study, 14 victims of a coal mine gas explosion with PTSD and 23 matched coal miners without experiencing the coal mine explosion were enrolled. The differences in fractional anisotropy (FA) within 7 sub-regions of the CC were compared between the two groups. Compared to the controls, PTSD coal miners exhibited significantly reduced FA values in the anterior sub-regions of the CC (P < 0.05, Bonferroni-corrected), which mainly interconnect the bilateral frontal cortices. Our findings indicated that the anterior part of the CC was more severely impaired than the posterior part in young adult-onset PTSD, which suggested the patterns of CC impairment may depend on the developmental stage of the structure when the PTSD occurs.     

Correlation between Corpus Callosum Sub-Segmental Area and Cognitive Processes in School-Age Children

We assessed the relationship between structural characteristics (area) and microstructure (apparent diffusion coefficient; ADC) of the corpus callosum (CC) in 57 healthy children aged 7.0 to 9.1 years, with diverse cognitive and academic abilities as well as executive functions evaluated with a neuropsychological battery for children. The CC was manually delineated and sub-segmented into six regions, and their ADC and area were measured. There were no significant differences between genders in the callosal region area or in ADC. The CC area and ADC, mainly of anterior regions, correlated with different cognitive abilities for each gender. Our results suggest that the relationship between cognitive abilities and CC characteristics is different between girls and boys and between the anterior and posterior regions of the CC. Furthermore, these findings strenghten the idea that regardless of the different interhemispheric connectivity schemes per gender, the results of cognitive tasks are very similar for girls and boys throughout childhood.     

Characterization of Glial Populations in the Aging and Remyelinating Mouse Corpus Callosum

Cells in the white matter of the adult brain have a characteristic distribution pattern in which several cells are contiguously connected to each other, making a linear array (LA) resembling pearls-on-a-string parallel to the axon axis. We have been interested in how this pattern of cell distribution changes during aging and remyelination after demyelination. In the present study, with a multiplex staining method, semi-quantitative analysis of the localization of oligodendrocyte lineage cells (oligodendrocyte progenitors, premyelinating oligodendrocytes, and mature oligodendrocytes), astrocytes, and microglia in 8-week-old (young adult) and 32-week-old (aged) corpus callosum showed that young adult cells still include immature oligodendrocytes and that LAs contain a higher proportion of microglia than isolated cells. In aged mice, premyelinating oligodendrocytes were decreased, but microglia continued to be present in the LAs. These results suggest that the presence of microglia is important for the characteristic cell localization pattern of LAs. In a cuprizone-induced demyelination model, we observed re-formation of LAs after completion of cuprizone treatment, concurrent with remyelination. These re-formed LAs again contained more microglia than the isolated cells. This finding supports the hypothesis that microglia contribute to the formation and maintenance of LAs. In addition, regardless of the distribution of cells (LAs or isolated cells), astrocytes were found to be more abundant than in the normal corpus callosum at 24 weeks after cuprizone treatment when remyelination is completed. This suggests that astrocytes are involved in maintaining the functions of remyelinated white matter.