New Tmc1 Deafness Mutations Impact Mechanotransduction in Auditory Hair Cells.

Transmembrane channel-like protein isoform 1 (TMC1) is a major component of the mechano-electrical transducer (MET) channel in cochlear hair cells and is subject to numerous mutations causing deafness. We report a new dominant human deafness mutation, TMC1 p.T422K, and have characterized the homologous mouse mutant, Tmc1 p.T416K, which caused deafness and outer hair cell (OHC) loss by the fourth postnatal week. MET channels showed decreased Ca2+ permeability and resting open probability, but no change in single-channel conductance or expression. Three adjacent deafness mutations are TMC1 p.L416R, p.G417R, and p.M418K, the last homologous to the mouse Beethoven that exhibits similar channel effects. All substitute a positive for a neutral residue, which could produce charge screening in the channel pore or influence binding of an accessory subunit. Channel properties were compared in mice of both sexes between dominant (Tmc1 p.T416K, Tmc1 p.D569N) and recessive (Tmc1 p.W554L, Tmc1 p.D528N) mutations of residues near the putative pore of the channel. Tmc1 p.W554L and p.D569N exhibit reduced maximum current with no effect on single-channel conductance, implying a smaller number of channels transported to the stereociliary tips; this may stem from impaired TMC1 binding to LHFPL5. Tmc1 p.D528N, located in the pore's narrowest region, uniquely caused large reductions in MET channel conductance and block by dihydrostreptomycin (DHS). For Tmc1 p.T416K and Tmc1 p.D528N, transduction loss occurred between P15 and P20. We propose two mechanisms linking channel mutations and deafness: decreased Ca2+ permeability, common to all mutants, and decreased resting open probability in low Ca2+, confined to dominant mutations.SIGNIFICANCE STATEMENT Transmembrane channel-like protein isoform 1 (TMC1) is thought to be a major component of the mechanotransducer channel in auditory hair cells, but the protein organization and channel structure are still uncertain. We made four mouse lines harboring Tmc1 point mutations that alter channel properties, causing hair cell degeneration and deafness. These include a mouse homolog of a new human deafness mutation pT416K that decreased channel Ca2+ permeability by introducing a positively-charged amino acid in the putative pore. All mutations are consistent with the channel structure predicted from modeling, but only one, p.D528N near the external face of the pore, substantially reduced channel conductance and Ca2+ permeability and virtually abolished block by dihydrostreptomycin (DHS), strongly endorsing its siting within the pore.

Characterization of a possible founder synonymous variant in TECTA in multiple individuals with autosomal recessive hearing loss.

Synonymous variants have been shown to alter the correct splicing of pre-mRNAs and generate disease-causing transcripts. These variants are not an uncommon etiology of genetic disease; however, they are frequently overlooked during genetic testing in the absence of functional and clinical data. Here, we describe the occurrence of a synonymous variant [NM_005422.4 (TECTA):c.327C>T, p.(Gly109=)] in seven individuals with hearing loss from six unrelated families. The variant is not located near exonic/intronic boundaries but is predicted to impact splicing by activating a cryptic splicing donor site in exon 4 of TECTA. In vitro minigene assays show that the variant disrupts the reading frame of the canonical transcript, which is predicted to cause a premature termination codon 48 amino acids downstream of the variant, leading to nonsense-mediated decay. The variant is present in population databases, predominantly in Latinos of African ancestry, but is rare in other ethnic groups. Our findings suggest that this synonymous variant is likely pathogenic for TECTA-associated autosomal recessive hearing loss and seems to have arisen as a founder variant in this specific Latino subpopulation. This study demonstrates that synonymous variants need careful splicing assessment and support from additional testing methodologies to determine their clinical impact.

Lung function, airway and peripheral basophils and eosinophils are associated with molecular pharmacogenomic endotypes of steroid response in severe asthma.

INTRODUCTION: Asthma is a complex disease with heterogeneous expression/severity. There is growing interest in defining asthma endotypes consistently associated with different responses to therapy, focusing on type 2 inflammation (Th2) as a key pathological mechanism. Current asthma endotypes are defined primarily by clinical/laboratory criteria. Each endotype is likely characterised by distinct molecular mechanisms that identify optimal therapies. METHODS: We applied unsupervised (without a priori clinical criteria) principal component analysis on sputum airway cells RNA-sequencing transcriptomic data from 19 asthmatics from the Severe Asthma Research Program at baseline and 6-8 weeks follow-up after a 40 mg dose of intramuscular corticosteroids. We investigated principal components PC1, PC3 for association with 55 clinical variables. RESULTS: PC3 was associated with baseline Th2 clinical features including blood (rank-sum p=0.0082) and airway (rank-sum p=0.0024) eosinophilia, FEV1 change (Kendall tau-b R=-0.333 (-0.592 to -0.012)) and follow-up FEV1 albuterol response (Kendall tau-b R=0.392 (0.079 to 0.634)). PC1 with blood basophlia (rank-sum p=0.0191). The top 5% genes contributing to PC1, PC3 were enriched for distinct immune system/inflammation ontologies suggesting distinct subject-specific clusters of transcriptomic response to corticosteroids. PC3 association with FEV1 change was reproduced in silico in a comparable independent 14-subject (baseline, 8 weeks after daily inhaled corticosteroids (ICS)) airway epithelial cells microRNAome dataset. CONCLUSIONS: Transcriptomic PCs from this unsupervised methodology define molecular pharmacogenomic endotypes that may yield novel biology underlying different subject-specific responses to corticosteroid therapy in asthma, and optimal personalised asthma care. Top contributing genes to these PCs may suggest new therapeutic targets.

Molecular characterization of pathogenic OTOA gene conversions in hearing loss patients.

Bi-allelic loss-of-function variants of OTOA are a well-known cause of moderate-to-severe hearing loss. Whereas non-allelic homologous recombination-mediated deletions of the gene are well known, gene conversions to pseudogene OTOAP1 have been reported in the literature but never fully described nor their pathogenicity assessed. Here, we report two unrelated patients with moderate hearing-loss, who were compound heterozygotes for a converted allele and a deletion of OTOA. The conversions were initially detected through sequencing depths anomalies at the OTOA locus after exome sequencing, then confirmed with long range polymerase chain reactions. Both conversions lead to loss-of-function by introducing a premature stop codon in exon 22 (p.Glu787*). Using genomic alignments and long read nanopore sequencing, we found that the two probands carry stretches of converted DNA of widely different lengths (at least 9 kbp and around 900 bp, respectively).

Disease-specific ACMG/AMP guidelines improve sequence variant interpretation for hearing loss.

PURPOSE: The ClinGen Variant Curation Expert Panels (VCEPs) provide disease-specific rules for accurate variant interpretation. Using the hearing loss-specific American College of Medical Genetics and Genomics/Association for Molecular Pathology (ACMG/AMP) guidelines, the Hearing Loss VCEP (HL VCEP) illustrates the utility of expert specifications in variant interpretation. METHODS: A total of 157 variants across nine HL genes, previously submitted to ClinVar, were curated by the HL VCEP. The curation process involved collecting published and unpublished data for each variant by biocurators, followed by bimonthly meetings of an expert curation subgroup that reviewed all evidence and applied the HL-specific ACMG/AMP guidelines to reach a final classification. RESULTS: Before expert curation, 75% (117/157) of variants had single or multiple variants of uncertain significance (VUS) submissions (17/157) or had conflicting interpretations in ClinVar (100/157). After applying the HL-specific ACMG/AMP guidelines, 24% (4/17) of VUS and 69% (69/100) of discordant variants were resolved into benign (B), likely benign (LB), likely pathogenic (LP), or pathogenic (P). Overall, 70% (109/157) variants had unambiguous classifications (B, LB, LP, P). We quantify the contribution of the HL-specified ACMG/AMP codes to variant classification. CONCLUSION: Expert specification and application of the HL-specific ACMG/AMP guidelines effectively resolved discordant interpretations in ClinVar. This study highlights the utility of ClinGen VCEPs in supporting more consistent clinical variant interpretation.

The eMERGE genotype set of 83,717 subjects imputed to ~40 million variants genome wide and association with the herpes zoster medical record phenotype.

The Electronic Medical Records and Genomics (eMERGE) network is a network of medical centers with electronic medical records linked to existing biorepository samples for genomic discovery and genomic medicine research. The network sought to unify the genetic results from 78 Illumina and Affymetrix genotype array batches from 12 contributing medical centers for joint association analysis of 83,717 human participants. In this report, we describe the imputation of eMERGE results and methods to create the unified imputed merged set of genome-wide variant genotype data. We imputed the data using the Michigan Imputation Server, which provides a missing single-nucleotide variant genotype imputation service using the minimac3 imputation algorithm with the Haplotype Reference Consortium genotype reference set. We describe the quality control and filtering steps used in the generation of this data set and suggest generalizable quality thresholds for imputation and phenotype association studies. To test the merged imputed genotype set, we replicated a previously reported chromosome 6 HLA-B herpes zoster (shingles) association and discovered a novel zoster-associated loci in an epigenetic binding site near the terminus of chromosome 3 (3p29).

ClinGen expert clinical validity curation of 164 hearing loss gene-disease pairs.

PURPOSE: Proper interpretation of genomic variants is critical to successful medical decision making based on genetic testing results. A fundamental prerequisite to accurate variant interpretation is the clear understanding of the clinical validity of gene-disease relationships. The Clinical Genome Resource (ClinGen) has developed a semiquantitative framework to assign clinical validity to gene-disease relationships. METHODS: The ClinGen Hearing Loss Gene Curation Expert Panel (HL GCEP) uses this framework to perform evidence-based curations of genes present on testing panels from 17 clinical laboratories in the Genetic Testing Registry. The HL GCEP curated and reviewed 142 genes and 164 gene-disease pairs, including 105 nonsyndromic and 59 syndromic forms of hearing loss. RESULTS: The final outcome included 82 Definitive (50%), 12 Strong (7%), 25 Moderate (15%), 32 Limited (20%), 10 Disputed (6%), and 3 Refuted (2%) classifications. The summary of each curation is date stamped with the HL GCEP approval, is live, and will be kept up-to-date on the ClinGen website ( https://search.clinicalgenome.org/kb/gene-validity ). CONCLUSION: This gene curation approach serves to optimize the clinical sensitivity of genetic testing while reducing the rate of uncertain or ambiguous test results caused by the interrogation of genes with insufficient evidence of a disease link.

Correction: ClinGen expert clinical validity curation of 164 hearing loss gene-disease pairs.

An amendment to this paper has been published and can be accessed via a link at the top of the paper.

Consensus interpretation of the p.Met34Thr and p.Val37Ile variants in GJB2 by the ClinGen Hearing Loss Expert Panel.

PURPOSE: Pathogenic variants in GJB2 are the most common cause of autosomal recessive sensorineural hearing loss. The classification of c.101T>C/p.Met34Thr and c.109G>A/p.Val37Ile in GJB2 are controversial. Therefore, an expert consensus is required for the interpretation of these two variants. METHODS: The ClinGen Hearing Loss Expert Panel collected published data and shared unpublished information from contributing laboratories and clinics regarding the two variants. Functional, computational, allelic, and segregation data were also obtained. Case-control statistical analyses were performed. RESULTS: The panel reviewed the synthesized information, and classified the p.Met34Thr and p.Val37Ile variants utilizing professional variant interpretation guidelines and professional judgment. We found that p.Met34Thr and p.Val37Ile are significantly overrepresented in hearing loss patients, compared with population controls. Individuals homozygous or compound heterozygous for p.Met34Thr or p.Val37Ile typically manifest mild to moderate hearing loss. Several other types of evidence also support pathogenic roles for these two variants. CONCLUSION: Resolving controversies in variant classification requires coordinated effort among a panel of international multi-institutional experts to share data, standardize classification guidelines, review evidence, and reach a consensus. We concluded that p.Met34Thr and p.Val37Ile variants in GJB2 are pathogenic for autosomal recessive nonsyndromic hearing loss with variable expressivity and incomplete penetrance.

Expert specification of the ACMG/AMP variant interpretation guidelines for genetic hearing loss.

Due to the high genetic heterogeneity of hearing loss (HL), current clinical testing includes sequencing large numbers of genes, which often yields a significant number of novel variants. Therefore, the standardization of variant interpretation is crucial to provide consistent and accurate diagnoses. The Hearing Loss Variant Curation Expert Panel was created within the Clinical Genome Resource to provide expert guidance for standardized genomic interpretation in the context of HL. As one of its major tasks, our Expert Panel has adapted the American College of Medical Genetics and Genomics/Association for Molecular Pathology (ACMG/AMP) guidelines for the interpretation of sequence variants in HL genes. Here, we provide a comprehensive illustration of the newly specified ACMG/AMP HL rules. Three rules remained unchanged, four rules were removed, and the remaining 21 rules were specified. These rules were further validated and refined using a pilot set of 51 variants assessed by curators and disease experts. Of the 51 variants evaluated in the pilot, 37% (19/51) changed category based upon application of the expert panel specified rules and/or aggregation of evidence across laboratories. These HL-specific ACMG/AMP rules will help standardize variant interpretation, ultimately leading to better care for individuals with HL.

Recurrent variants in OTOF are significant contributors to prelingual nonsydromic hearing loss in Saudi patients.

PurposeHearing loss is more prevalent in the Saudi Arabian population than in other populations; however, the full range of genetic etiologies in this population is unknown. We report the genetic findings from 33 Saudi hearing-loss probands of tribal ancestry, with predominantly prelingual severe to profound hearing loss.MethodsTesting was performed over the course of 2012-2016, and involved initial GJB2 sequence and GJB6-D13S1830 deletion screening, with negative cases being reflexed to a next-generation sequencing panel with 70, 71, or 87 hearing-loss genes.ResultsA "positive" result was reached in 63% of probands, with two recurrent OTOF variants (p.Glu57* and p.Arg1792His) accountable for a third of all "positive" cases. The next most common cause was pathogenic variants in MYO7A and SLC26A4, each responsible for three "positive" cases. Interestingly, only one "positive" diagnosis had a DFNB1-related cause, due to a homozygous GJB6-D13S1830 deletion, and no sequence variants in GJB2 were detected.ConclusionOur findings implicate OTOF as a potential major contributor to hearing loss in the Saudi population, while highlighting the low contribution of GJB2, thus offering important considerations for clinical testing strategies for Saudi patients. Further screening of Saudi patients is needed to characterize the genetic spectrum in this population.

Allele-Specific Droplet Digital PCR Combined with a Next-Generation Sequencing-Based Algorithm for Diagnostic Copy Number Analysis in Genes with High Homology: Proof of Concept Using Stereocilin.

BACKGROUND: Copy number variants (CNVs) can substantially contribute to the pathogenic variant spectrum in several disease genes. The detection of this type of variant is complicated in genes with high homology to other genomic sequences, yet such genomics regions are more likely to lead to CNVs, making it critical to address detection in these settings. METHODS: We developed a copy number analysis approach for high homology genes/regions that consisted of next-generation sequencing (NGS)-based dosage analysis accompanied by allele-specific droplet digital PCR (ddPCR) confirmatory testing. We applied this approach to copy number analysis in STRC, a gene with 98.9% homology to a nonfunctional pseudogene, pSTRC, and characterized its accuracy in detecting different copy number states by use of known samples. RESULTS: Using a cohort of 517 patients with hearing loss, we prospectively demonstrated the clinical utility of the approach, which contributed 30 of the 122 total positives (6%) to the diagnostic yield, increasing the overall yield from 17.6% to 23.6%. Positive STRC genotypes included homozygous (n = 15) or compound heterozygous (n = 8) deletions, or heterozygous deletions in trans with pathogenic sequence variants (n = 7). Finally, this approach limited ddPCR testing to cases with NGS copy number findings, thus markedly reducing the number of costly and laborious, albeit specific, ddPCR tests. CONCLUSIONS: NGS-based CNV detection followed by allele-specific ddPCR confirmatory testing is a reliable and affordable approach for copy number analysis in medically relevant genes with homology issues.

Mitochondrial gene sequence variants in children with severe malaria anaemia with or without lactic acidosis: a case control study.

BACKGROUND: Evolutionary pressure by Plasmodium falciparum malaria is known to have favoured a large number of human gene adaptations, but there is surprisingly little investigation of the effect of malaria on human mitochondrial sequence variation. Plasmodium falciparum infection can cause severe malaria anaemia (SMA) with insufficient tissue oxygenation, lactic acidosis and death. Despite equal degrees of severe anaemia, some individuals develop lactic acidosis while others do not. A case-control study design was used to investigate whether differences in host mitochondrial gene sequences were associated with lactic acidosis in SMA. Full mitochondrial sequences were obtained from 36 subjects with SMA complicated by lactic acidosis and 37 subjects with SMA without lactic acidosis. The two groups were matched for age, sex, and degree of anaemia. RESULTS: Compared with the reference sequence, a median of 60 nucleotide variants per individual (interquartile range 4-91) was found, with an average frequency of 3.97 variants per 1000 nucleotides. The frequency and distribution of non-synonymous DNA variants in genes associated with oxidative phosphorylation were not statistically different between the two groups. Non-synonymous variants predicted to have the most disruptive effect on proteins responsible for oxidative phosphorylation were present at a similar frequency in both groups. CONCLUSIONS: Lactic acidosis in SMA does not appear to be consistently associated with the high prevalence of any mitochondrial gene variant.

HIV Integration Site Analysis of Cellular Models of HIV Latency with a Probe-Enriched Next-Generation Sequencing Assay.

UNLABELLED: Antiretroviral therapy (ART) is successful in the suppression of HIV but cannot target and eradicate the latent proviral reservoir. The location of retroviral integration into the human genome is thought to play a role in the clonal expansion of infected cells and HIV persistence. We developed a high-throughput targeted sequence capture assay that uses a pool of HIV-specific probes to enrich Illumina libraries prior to deep sequencing. Using an expanded clonal population of ACH-2 cells, we demonstrate that this sequence capture assay has an extremely low false-positive rate. This assay assessed four cellular models commonly used to study HIV latency and latency-reversing agents: ACH-2 cells, J-Lat cells, the Bcl-2-transduced primary CD4(+)model, and the cultured TCM(central memory) CD4(+)model. HIV integration site characteristics and genes were compared between these cellular models and to previously reported patient data sets. Across these cellular models, there were significant differences in integration site characteristics, including orientation relative to that of the host gene, the proportion of clonally expanded sites, and the proportion located within genic regions and exons. Despite a greater diversity of minority integration sites than expected in ACH-2 cells, their integration site characteristics consistently differed from those of the other models and from the patient samples. Gene ontology analysis of highly represented genes from the patient samples found little overlap with HIV-containing genes from the cell lines. These findings show that integration site differences exist among the commonly used cellular models of HIV latency and in comparison to integration sites found in patient samples. IMPORTANCE: Despite the success of ART, currently there is no successful therapy to eradicate integrated proviruses. Cellular models of HIV latency are used to test the efficacy of latency-reversing agents, but it is unclear how well these models reflect HIV integration into the human genome in vivo We have developed a novel probe-based sequence enrichment assay to sequence and analyze integrated HIV. We compared HIV integration site characteristics between four cellular models and to previously described patient data sets. Significant differences were detected in the distribution of HIV integration sites between cellular models of HIV latency and compared to data sets from patient samples. The results from this study have implications for how well these cellular models of HIV infection truly reflect HIV integration in vivo and their applicability in drug discovery for novel latency-reversing agents.

Using large sequencing data sets to refine intragenic disease regions and prioritize clinical variant interpretation.

PURPOSE: Classification of novel variants is a major challenge facing the widespread adoption of comprehensive clinical genomic sequencing and the field of personalized medicine in general. This is largely because most novel variants do not have functional, genetic, or population data to support their clinical classification. METHODS: To improve variant interpretation, we leveraged the Exome Aggregation Consortium (ExAC) data set (N = ~60,000) as well as 7,000 clinically curated variants in 132 genes identified in more than 11,000 probands clinically tested for cardiomyopathies, rasopathies, hearing loss, or connective tissue disorders to perform a systematic evaluation of domain level disease associations. RESULTS: We statistically identify regions that are most sensitive to functional variation in the general population and also most commonly impacted in symptomatic individuals. Our data show that a significant number of exons and domains in genes strongly associated with disease can be defined as disease-sensitive or disease-tolerant, leading to potential reclassification of at least 26% (450 out of 1,742) of variants of uncertain clinical significance in the 132 genes. CONCLUSION: This approach leverages domain functional annotation and associated disease in each gene to prioritize candidate disease variants, increasing the sensitivity and specificity of novel variant assessment within these genes.Genet Med advance online publication 22 September 2016.

Targeted Droplet-Digital PCR as a Tool for Novel Deletion Discovery at the DFNB1 Locus.

Pathogenic variants at the DFNB1 locus encompassing the GJB2 and GJB6 genes account for 50% of autosomal-recessive, congenital nonsyndromic hearing loss in the United States. Most cases are caused by sequence variants within the GJB2 gene, but a significant number of DFNB1 patients carry a large deletion (GJB6-D13S1830) in trans with a GJB2 variant. This deletion lies upstream of GJB2 and was shown to reduce GJB2 expression by disrupting unidentified regulatory elements. First-tier genetic testing for hearing loss includes GJB2 sequence and GJB6-D13S1830 deletion analysis; however, several other deletions in this locus, each with distinct breakpoints, have been reported in DFNB1 patients and are missed by current panels. Here, we report the development of a targeted droplet digital polymerase chain reaction-based assay for comprehensive copy-number analysis at the DFNB1 locus that detects all deletions reported to date. This assay increased detection rates in a multiethnic cohort of 87 hearing loss patients with only one identified pathogenic GJB2 variant. We identify two deletions, one of which is novel, in two patients (2/87 or 2.3%), suggesting that other pathogenic deletions at the DFNB1 locus may be missed. Mapping the assayed DFNB1 deletions also revealed a ∼ 95 kb critical region, which may harbor the GJB2 regulatory element(s).

Improving hearing loss gene testing: a systematic review of gene evidence toward more efficient next-generation sequencing-based diagnostic testing and interpretation.

PURPOSE: With next generation sequencing technology improvement and cost reductions, it has become technically feasible to sequence a large number of genes in one diagnostic test. This is especially relevant for diseases with large genetic and/or phenotypic heterogeneity, such as hearing loss. However, variant interpretation remains the major bottleneck. This is further exacerbated by the lack in the clinical genetics community of consensus criteria for defining the evidence necessary to include genes on targeted disease panels or in genomic reports, and the consequent risk of reporting variants in genes with no relevance to disease. METHODS: We describe a systematic evidence-based approach for assessing gene-disease associations and for curating relevant genes for different disease aspects, including mode of inheritance, phenotypic severity, and mutation spectrum. RESULTS: By applying this approach to clinically available hearing loss gene panels with a total of 163 genes, we show that a significant number (45%) of genes lack sufficient evidence of association with disease and thus are expected to increase uncertainty and patient anxiety, in addition to intensifying the interpretation burden. Information about all curated genes is summarized. Our retrospective analysis of 539 hearing loss cases tested by our previous OtoGenomeV2 panel demonstrates the impact of including genes with weak disease association in laboratory wet-bench and interpretation processes. CONCLUSION: Our study is, to our knowledge, the first to highlight the urgent need for defining the clinical validity of gene-disease relationships for more efficient and accurate clinical testing and reporting.Genet Med 18 6, 545-553.

VisCap: inference and visualization of germ-line copy-number variants from targeted clinical sequencing data.

PURPOSE: To develop and validate VisCap, a software program targeted to clinical laboratories for inference and visualization of germ-line copy-number variants (CNVs) from targeted next-generation sequencing data. METHODS: VisCap calculates the fraction of overall sequence coverage assigned to genomic intervals and computes log2 ratios of these values to the median of reference samples profiled using the same test configuration. Candidate CNVs are called when log2 ratios exceed user-defined thresholds. RESULTS: We optimized VisCap using 14 cases with known CNVs, followed by prospective analysis of 1,104 cases referred for diagnostic DNA sequencing. To verify calls in the prospective cohort, we used droplet digital polymerase chain reaction (PCR) to confirm 10/27 candidate CNVs and 72/72 copy-neutral genomic regions scored by VisCap. We also used a genome-wide bead array to confirm the absence of CNV calls across panels applied to 10 cases. To improve specificity, we instituted a visual scoring system that enabled experienced reviewers to differentiate true-positive from false-positive calls with minimal impact on laboratory workflow. CONCLUSIONS: VisCap is a sensitive method for inferring CNVs from targeted sequence data from targeted gene panels. Visual scoring of data underlying CNV calls is a critical step to reduce false-positive calls for follow-up testing.Genet Med 18 7, 712-719.Genetics in Medicine (2016); 18 7, 712-719. doi:10.1038/gim.2015.156.

Rare germline disorders implicate long non-coding RNAs disrupted by chromosomal structural rearrangements.

In recent years, there has been increased focus on exploring the role the non-protein-coding genome plays in Mendelian disorders. One class of particular interest is long non-coding RNAs (lncRNAs), which has recently been implicated in the regulation of diverse molecular processes. However, because lncRNAs do not encode protein, there is uncertainty regarding what constitutes a pathogenic lncRNA variant, and thus annotating such elements is challenging. The Developmental Genome Anatomy Project (DGAP) and similar projects recruit individuals with apparently balanced chromosomal abnormalities (BCAs) that disrupt or dysregulate genes in order to annotate the human genome. We hypothesized that rearrangements disrupting lncRNAs could be the underlying genetic etiology for the phenotypes of a subset of these individuals. Thus, we assessed 279 cases with BCAs and selected 191 cases with simple BCAs (breakpoints at only two genomic locations) for further analysis of lncRNA disruptions. From these, we identified 66 cases in which the chromosomal rearrangements directly disrupt lncRNAs. Strikingly, the lncRNAs MEF2C-AS1 and ENSG00000257522 are each disrupted in two unrelated cases. Furthermore, in 30 cases, no genes of any other class aside from lncRNAs are directly disrupted, consistent with the hypothesis that lncRNA disruptions could underly the phenotypes of these individuals. To showcase the power of this genomic approach for annotating lncRNAs, here we focus on clinical reports and genetic analysis of two individuals with BCAs and additionally highlight six individuals with likely developmental etiologies due to lncRNA disruptions.

SOX12 and NRSN2 are candidate genes for 20p13 subtelomeric deletions associated with developmental delay.

20p13 telomeric/subtelomeric deletions are clinically significant but are currently under-investigated. So far only five molecularly delineated cases have been reported in literature and no candidate genes have been sufficiently implicated. Here, we present six new deletion cases identified by chromosomal microarray analysis (CMA). We also review 32 cases combined from literature and databases. We found that most 20p13 deletion patients exhibit significant developmental delay. Dysmorphic features are common but a consistent pattern was not recognized. Reduced cognitive ability was frequent. Based on pathogenic deletions delineated in this study, we mapped the smallest overlapping region and identified two nervous system expressing genes (SOX12 and NRSN2) as candidate genes that may be involved in the developmental defects in 20p13 microdeletion.