With recent rapid advances in genomic technologies, precise delineation of structural chromosome rearrangements at the nucleotide level is becoming increasingly feasible. In this era of "next-generation cytogenetics" (i.e., an integration of traditional cytogenetic techniques and next-generation sequencing), a consensus nomenclature is essential for accurate communication and data sharing. Currently, nomenclature for describing the sequencing data of these aberrations is lacking. Herein, we present a system called Next-Gen Cytogenetic Nomenclature, which is concordant with the International System for Human Cytogenetic Nomenclature (2013). This system starts with the alignment of rearrangement sequences by BLAT or BLAST (alignment tools) and arrives at a concise and detailed description of chromosomal changes. To facilitate usage and implementation of this nomenclature, we are developing a program designated BLA(S)T Output Sequence Tool of Nomenclature (BOSToN), a demonstrative version of which is accessible online. A standardized characterization of structural chromosomal rearrangements is essential both for research analyses and for application in the clinical setting.

Persons with neurodevelopmental disorders or autism spectrum disorder (ASD) often harbor chromosomal microdeletions, yet the individual genetic contributors within these regions have not been systematically evaluated. We established a consortium of clinical diagnostic and research laboratories to accumulate a large cohort with genetic alterations of chromosomal region 2q23.1 and acquired 65 subjects with microdeletion or translocation. We sequenced translocation breakpoints; aligned microdeletions to determine the critical region; assessed effects on mRNA expression; and examined medical records, photos, and clinical evaluations. We identified a single gene, methyl-CpG-binding domain 5 (MBD5), as the only locus that defined the critical region. Partial or complete deletion of MBD5 was associated with haploinsufficiency of mRNA expression, intellectual disability, epilepsy, and autistic features. Fourteen alterations, including partial deletions of noncoding regions not typically captured or considered pathogenic by current diagnostic screening, disrupted MBD5 alone. Expression profiles and clinical characteristics were largely indistinguishable between MBD5-specific alteration and deletion of the entire 2q23.1 interval. No copy-number alterations of MBD5 were observed in 7878 controls, suggesting MBD5 alterations are highly penetrant. We surveyed MBD5 coding variations among 747 ASD subjects compared to 2043 non-ASD subjects analyzed by whole-exome sequencing and detected an association with a highly conserved methyl-CpG-binding domain missense variant, p.79Gly>Glu (c.236G>A) (p = 0.012). These results suggest that genetic alterations of MBD5 cause features of 2q23.1 microdeletion syndrome and that this epigenetic regulator significantly contributes to ASD risk, warranting further consideration in research and clinical diagnostic screening and highlighting the importance of chromatin remodeling in the etiology of these complex disorders.

We report a female patient with a de novo balanced translocation, 46,X,t(X;2)(p11.2;q37)dn, who exhibits seizures, capillary abnormality, developmental delay, infantile hypotonia, and obesity. The 2q37 breakpoint observed in association with the seizure phenotype is of particular interest, because it lies near loci implicated in epilepsy in humans and mice. Fluorescence in situ hybridization mapping of the translocation breakpoints showed that no known genes are disrupted at Xp11.2, whereas diacylglycerol kinase delta (DGKD) is disrupted at 2q37. Expression studies in Drosophila and mouse suggest that DGKD is involved in central nervous system development and function. Electroencephalographic assessment of Dgkd mutant mice revealed abnormal epileptic discharges and electrographic seizures in three of six homozygotes. These findings implicate DGKD disruption by the t(X;2)(p11.2;q37)dn in the observed phenotype and support a more general role for DGKD in the etiology of seizures.

Structural variation, composed of balanced and unbalanced genomic rearrangements, is an important contributor to human genetic diversity with prominent roles in somatic and congenital disease. At the nucleotide level, structural variants (SVs) have been shown to frequently harbor additional breakpoints and copy-number imbalances, a complexity predicted to emerge wholly as a single-cell division event. Chromothripsis, chromoplexy, and chromoanasynthesis, collectively referred to as chromoanagenesis, are three major mechanisms that explain the occurrence of complex germline and somatic SVs. While chromothripsis and chromoplexy have been shown to be key signatures of cancer, chromoanagenesis has been detected in numerous cases of developmental disease and phenotypically normal individuals. Such observations advocate for a deeper study of the polymorphic and pathogenic properties of complex germline SVs, many of which go undetected by traditional clinical molecular and cytogenetic methods. This review focuses on congenital chromoanagenesis, mechanisms leading to occurrence of these complex rearrangements, and their impact on chromosome organization and genome function. We highlight future applications of routine screening of complex and balanced SVs in the clinic, as these represent a potential and often neglected genetic disease source, a true "iceberg under water."

Apparently balanced chromosomal rearrangements in individuals with major congenital anomalies represent natural experiments of gene disruption and dysregulation. These individuals can be studied to identify novel genes critical in human development and to annotate further the function of known genes. Identification and characterization of these genes is the goal of the Developmental Genome Anatomy Project (DGAP). DGAP is a multidisciplinary effort that leverages the recent advances resulting from the Human Genome Project to increase our understanding of birth defects and the process of human development. Clinically significant phenotypes of individuals enrolled in DGAP are varied and, in most cases, involve multiple organ systems. Study of these individuals' chromosomal rearrangements has resulted in the mapping of 77 breakpoints from 40 chromosomal rearrangements by FISH with BACs and fosmids, array CGH, Southern-blot hybridization, MLPA, RT-PCR, and suppression PCR. Eighteen chromosomal breakpoints have been cloned and sequenced. Unsuspected genomic imbalances and cryptic rearrangements were detected, but less frequently than has been reported previously. Chromosomal rearrangements, both balanced and unbalanced, in individuals with multiple congenital anomalies continue to be a valuable resource for gene discovery and annotation.

Hearing loss is one of the top contributors to years lived with disability and is a risk factor for dementia. Molecular evidence on the cellular origins of hearing loss in humans is growing. Here, we performed a genome-wide association meta-analysis of clinically diagnosed and self-reported hearing impairment on 723,266 individuals and identified 48 significant loci, 10 of which are novel. A large proportion of associations comprised missense variants, half of which lie within known familial hearing loss loci. We used single-cell RNA-sequencing data from mouse cochlea and brain and mapped common-variant genomic results to spindle, root, and basal cells from the stria vascularis, a structure in the cochlea necessary for normal hearing. Our findings indicate the importance of the stria vascularis in the mechanism of hearing impairment, providing future paths for developing targets for therapeutic intervention in hearing loss.

It has long been appreciated that genetic analysis of fetal or trophoblast cells in maternal blood could revolutionize prenatal diagnosis. We implemented a protocol for single circulating trophoblast (SCT) testing using positive selection by magnetic-activated cell sorting and single-cell low-coverage whole-genome sequencing to detect fetal aneuploidies and copy-number variants (CNVs) at ∼1 Mb resolution. In 95 validation cases, we identified on average 0.20 putative trophoblasts/mL, of which 55% were of high quality and scorable for both aneuploidy and CNVs. We emphasize the importance of analyzing individual cells because some cells are apoptotic, in S-phase, or otherwise of poor quality. When two or more high-quality trophoblast cells were available for singleton pregnancies, there was complete concordance between all trophoblasts unless there was evidence of confined placental mosaicism. SCT results were highly concordant with available clinical data from chorionic villus sampling (CVS) or amniocentesis procedures. Although determining the exact sensitivity and specificity will require more data, this study further supports the potential for SCT testing to become a diagnostic prenatal test.

Genetic mutations responsible for oblique facial clefts (ObFC), a unique class of facial malformations, are largely unknown. We show that loss-of-function mutations in SPECC1L are pathogenic for this human developmental disorder and that SPECC1L is a critical organizer of vertebrate facial morphogenesis. During murine embryogenesis, Specc1l is expressed in cell populations of the developing facial primordial, which proliferate and fuse to form the face. In zebrafish, knockdown of a SPECC1L homolog produces a faceless phenotype with loss of jaw and facial structures, and knockdown in Drosophila phenocopies mutants in the integrin signaling pathway that exhibit cell-migration and -adhesion defects. Furthermore, in mammalian cells, SPECC1L colocalizes with both tubulin and actin, and its deficiency results in defective actin-cytoskeleton reorganization, as well as abnormal cell adhesion and migration. Collectively, these data demonstrate that SPECC1L functions in actin-cytoskeleton reorganization and is required for proper facial morphogenesis.

Uterine leiomyomata (UL), the most prevalent pelvic tumors in women of reproductive age, pose a major public health problem given their high frequency, associated morbidities, and most common indication for hysterectomies. A genetic component to UL predisposition is supported by analyses of ethnic predisposition, twin studies, and familial aggregation. A genome-wide SNP linkage panel was genotyped and analyzed in 261 white UL-affected sister-pair families from the Finding Genes for Fibroids study. Two significant linkage regions were detected in 10p11 (LOD = 4.15) and 3p21 (LOD = 3.73), and five additional linkage regions were identified with LOD scores > 2.00 in 2q37, 5p13, 11p15, 12q14, and 17q25. Genome-wide association studies were performed in two independent cohorts of white women, and a meta-analysis was conducted. One SNP (rs4247357) was identified with a p value (p = 3.05 × 10(-8)) that reached genome-wide significance (odds ratio = 1.299). The candidate SNP is under a linkage peak and in a block of linkage disequilibrium in 17q25.3, which spans fatty acid synthase (FASN), coiled-coil-domain-containing 57 (CCDC57), and solute-carrier family 16, member 3 (SLC16A3). By tissue microarray immunohistochemistry, we found elevated (3-fold) FAS levels in UL-affected tissue compared to matched myometrial tissue. FAS transcripts and/or protein levels are upregulated in various neoplasms and implicated in tumor cell survival. FASN represents the initial UL risk allele identified in white women by a genome-wide, unbiased approach and opens a path to management and potential therapeutic intervention.

Congenital anomalies of the kidney and urinary tract (CAKUT) include vesicoureteral reflux (VUR). VUR is a complex, genetically heterogeneous developmental disorder characterized by the retrograde flow of urine from the bladder into the ureter and is associated with reflux nephropathy, the cause of 15% of end-stage renal disease in children and young adults. We investigated a man with a de novo translocation, 46,X,t(Y;3)(p11;p12)dn, who exhibits multiple congenital abnormalities, including severe bilateral VUR with ureterovesical junction defects. This translocation disrupts ROBO2, which encodes a transmembrane receptor for SLIT ligand, and produces dominant-negative ROBO2 proteins that abrogate SLIT-ROBO signaling in vitro. In addition, we identified two novel ROBO2 intracellular missense variants that segregate with CAKUT and VUR in two unrelated families. Adult heterozygous and mosaic mutant mice with reduced Robo2 gene dosage also exhibit striking CAKUT-VUR phenotypes. Collectively, these results implicate the SLIT-ROBO signaling pathway in the pathogenesis of a subset of human VUR.

Autism is a neurodevelopmental disorder of complex etiology in which genetic factors play a major role. We have implicated the neurexin 1 (NRXN1) gene in two independent subjects who display an autism spectrum disorder (ASD) in association with a balanced chromosomal abnormality involving 2p16.3. In the first, with karyotype 46,XX,ins(16;2)(q22.1;p16.1p16.3)pat, NRXN1 is directly disrupted within intron 5. Importantly, the father possesses the same chromosomal abnormality in the absence of ASD, indicating that the interruption of alpha-NRXN1 is not fully penetrant and must interact with other factors to produce ASD. The breakpoint in the second subject, with 46,XY,t(1;2)(q31.3;p16.3)dn, occurs approximately 750 kb 5' to NRXN1 within a 2.6 Mb genomic segment that harbors no currently annotated genes. A scan of the NRXN1 coding sequence in a cohort of ASD subjects, relative to non-ASD controls, revealed that amino acid alterations in neurexin 1 are not present at high frequency in ASD. However, a number of rare sequence variants in the coding region, including two missense changes in conserved residues of the alpha-neurexin 1 leader sequence and of an epidermal growth factor (EGF)-like domain, respectively, suggest that even subtle changes in NRXN1 might contribute to susceptibility to ASD.

In this exciting era of "next-gen cytogenetics," integrating genomic sequencing into the prenatal diagnostic setting is possible within an actionable time frame and can provide precise delineation of balanced chromosomal rearrangements at the nucleotide level. Given the increased risk of congenital abnormalities in newborns with de novo balanced chromosomal rearrangements, comprehensive interpretation of breakpoints could substantially improve prediction of phenotypic outcomes and support perinatal medical care. Herein, we present and evaluate sequencing results of balanced chromosomal rearrangements in ten prenatal subjects with respect to the location of regulatory chromatin domains (topologically associated domains [TADs]). The genomic material from all subjects was interpreted to be "normal" by microarray analyses, and their rearrangements would not have been detected by cell-free DNA (cfDNA) screening. The findings of our systematic approach correlate with phenotypes of both pregnancies with untoward outcomes (5/10) and with healthy newborns (3/10). Two pregnancies, one with a chromosomal aberration predicted to be of unknown clinical significance and another one predicted to be likely benign, were terminated prior to phenotype-genotype correlation (2/10). We demonstrate that the clinical interpretation of structural rearrangements should not be limited to interruption, deletion, or duplication of specific genes and should also incorporate regulatory domains of the human genome with critical ramifications for the control of gene expression. As detailed in this study, our molecular approach to both detecting and interpreting the breakpoints of structural rearrangements yields unparalleled information in comparison to other commonly used first-tier diagnostic methods, such as non-invasive cfDNA screening and microarray analysis, to provide improved genetic counseling for phenotypic outcome in the prenatal setting.

Unexplained infertility affects 2%-3% of reproductive-aged couples. One approach to identifying genes involved in infertility is to study subjects with this clinical phenotype and a de novo balanced chromosomal aberration (BCA). While BCAs may reduce fertility by production of unbalanced gametes, a chromosomal rearrangement may also disrupt or dysregulate genes important in fertility. One such subject, DGAP230, has severe oligozoospermia and 46,XY,t(20;22)(q13.3;q11.2). We identified exclusive overexpression of SYCP2 from the der(20) allele that is hypothesized to result from enhancer adoption. Modeling the dysregulation in budding yeast resulted in disrupted structural integrity of the synaptonemal complex, a common cause of defective spermatogenesis in mammals. Exome sequencing of infertile males revealed three heterozygous SYCP2 frameshift variants in additional subjects with cryptozoospermia and azoospermia. In sum, this investigation illustrates the power of precision cytogenetics for annotation of the infertile genome, suggests that these mechanisms should be considered as an alternative etiology to that of segregation of unbalanced gametes in infertile men harboring a BCA, and provides evidence of SYCP2-mediated male infertility in humans.

Although somatic mutations in a number of genes have been associated with development of human tumors, such as lipomas, relatively few examples exist of germline mutations in these genes. Here we describe an 8-year-old boy who has a de novo pericentric inversion of chromosome 12, with breakpoints at p11.22 and q14.3, and a phenotype including extreme somatic overgrowth, advanced endochondral bone and dental ages, a cerebellar tumor, and multiple lipomas. His chromosomal inversion was found to truncate HMGA2, a gene that encodes an architectural factor involved in the etiology of many benign mesenchymal tumors and that maps to the 12q14.3 breakpoint. Similar truncations of murine Hmga2 in transgenic mice result in somatic overgrowth and, in particular, increased abundance of fat and lipomas, features strikingly similar to those observed in the child. This represents the first report of a constitutional rearrangement affecting HMGA2 and demonstrates the role of this gene in human growth and development. Systematic genetic analysis and clinical studies of this child may offer unique insights into the role of HMGA2 in adipogenesis, osteogenesis, and general growth control.

Potocki-Shaffer syndrome (PSS) is a contiguous gene disorder due to the interstitial deletion of band p11.2 of chromosome 11 and is characterized by multiple exostoses, parietal foramina, intellectual disability (ID), and craniofacial anomalies (CFAs). Despite the identification of individual genes responsible for multiple exostoses and parietal foramina in PSS, the identity of the gene(s) associated with the ID and CFA phenotypes has remained elusive. Through characterization of independent subjects with balanced translocations and supportive comparative deletion mapping of PSS subjects, we have uncovered evidence that the ID and CFA phenotypes are both caused by haploinsufficiency of a single gene, PHF21A, at 11p11.2. PHF21A encodes a plant homeodomain finger protein whose murine and zebrafish orthologs are both expressed in a manner consistent with a function in neurofacial and craniofacial development, and suppression of the latter led to both craniofacial abnormalities and neuronal apoptosis. Along with lysine-specific demethylase 1 (LSD1), PHF21A, also known as BHC80, is a component of the BRAF-histone deacetylase complex that represses target-gene transcription. In lymphoblastoid cell lines from two translocation subjects in whom PHF21A was directly disrupted by the respective breakpoints, we observed derepression of the neuronal gene SCN3A and reduced LSD1 occupancy at the SCN3A promoter, supporting a direct functional consequence of PHF21A haploinsufficiency on transcriptional regulation. Our finding that disruption of PHF21A by translocations in the PSS region is associated with ID adds to the growing list of ID-associated genes that emphasize the critical role of transcriptional regulation and chromatin remodeling in normal brain development and cognitive function.

Large intergenic noncoding (linc) RNAs represent a newly described class of ribonucleic acid whose importance in human disease remains undefined. We identified a severely developmentally delayed 16-year-old female with karyotype 46,XX,t(2;11)(p25.1;p15.1)dn in the absence of clinically significant copy number variants (CNVs). DNA capture followed by next-generation sequencing of the translocation breakpoints revealed disruption of a single noncoding gene on chromosome 2, LINC00299, whose RNA product is expressed in all tissues measured, but most abundantly in brain. Among a series of additional, unrelated subjects referred for clinical diagnostic testing who showed CNV affecting this locus, we identified four with exon-crossing deletions in association with neurodevelopmental abnormalities. No disruption of the LINC00299 coding sequence was seen in almost 14,000 control subjects. Together, these subjects with disruption of LINC00299 implicate this particular noncoding RNA in brain development and raise the possibility that, as a class, abnormalities of lincRNAs may play a significant role in human developmental disorders.

The contribution of balanced chromosomal rearrangements to complex disorders remains unclear because they are not detected routinely by genome-wide microarrays and clinical localization is imprecise. Failure to consider these events bypasses a potentially powerful complement to single nucleotide polymorphism and copy-number association approaches to complex disorders, where much of the heritability remains unexplained. To capitalize on this genetic resource, we have applied optimized sequencing and analysis strategies to test whether these potentially high-impact variants can be mapped at reasonable cost and throughput. By using a whole-genome multiplexing strategy, rearrangement breakpoints could be delineated at a fraction of the cost of standard sequencing. For rearrangements already mapped regionally by karyotyping and fluorescence in situ hybridization, a targeted approach enabled capture and sequencing of multiple breakpoints simultaneously. Importantly, this strategy permitted capture and unique alignment of up to 97% of repeat-masked sequences in the targeted regions. Genome-wide analyses estimate that only 3.7% of bases should be routinely omitted from genomic DNA capture experiments. Illustrating the power of these approaches, the rearrangement breakpoints were rapidly defined to base pair resolution and revealed unexpected sequence complexity, such as co-occurrence of inversion and translocation as an underlying feature of karyotypically balanced alterations. These findings have implications ranging from genome annotation to de novo assemblies and could enable sequencing screens for structural variations at a cost comparable to that of microarrays in standard clinical practice.

Interpretation of variants of uncertain significance, especially chromosomal rearrangements in non-coding regions of the human genome, remains one of the biggest challenges in modern molecular diagnosis. To improve our understanding and interpretation of such variants, we used high-resolution three-dimensional chromosomal structural data and transcriptional regulatory information to predict position effects and their association with pathogenic phenotypes in 17 subjects with apparently balanced chromosomal abnormalities. We found that the rearrangements predict disruption of long-range chromatin interactions between several enhancers and genes whose annotated clinical features are strongly associated with the subjects' phenotypes. We confirm gene-expression changes for a couple of candidate genes to exemplify the utility of our analysis of position effect. These results highlight the important interplay between chromosomal structure and disease and demonstrate the need to utilize chromatin conformational data for the prediction of position effects in the clinical interpretation of non-coding chromosomal rearrangements.

Genomic rearrangements involving AUTS2 (7q11.22) are associated with autism and intellectual disability (ID), although evidence for causality is limited. By combining the results of diagnostic testing of 49,684 individuals, we identified 24 microdeletions that affect at least one exon of AUTS2, as well as one translocation and one inversion each with a breakpoint within the AUTS2 locus. Comparison of 17 well-characterized individuals enabled identification of a variable syndromic phenotype including ID, autism, short stature, microcephaly, cerebral palsy, and facial dysmorphisms. The dysmorphic features were more pronounced in persons with 3'AUTS2 deletions. This part of the gene is shown to encode a C-terminal isoform (with an alternative transcription start site) expressed in the human brain. Consistent with our genetic data, suppression of auts2 in zebrafish embryos caused microcephaly that could be rescued by either the full-length or the C-terminal isoform of AUTS2. Our observations demonstrate a causal role of AUTS2 in neurocognitive disorders, establish a hitherto unappreciated syndromic phenotype at this locus, and show how transcriptional complexity can underpin human pathology. The zebrafish model provides a valuable tool for investigating the etiology of AUTS2 syndrome and facilitating gene-function analysis in the future.

Age-related hearing impairment (ARHI) is the most common sensory impairment in the aging population; a third of individuals are affected by disabling hearing loss by the age of 65. It causes social isolation and depression and has recently been identified as a risk factor for dementia. The genetic risk factors and underlying pathology of ARHI are largely unknown, meaning that targets for new therapies remain elusive, yet heritability estimates range between 35% and 55%. We performed genome-wide association studies (GWASs) for two self-reported hearing phenotypes, using more than 250,000 UK Biobank (UKBB) volunteers aged between 40 and 69 years. Forty-four independent genome-wide significant loci (p < 5E-08) were identified, considerably increasing the number of established trait loci. Thirty-four loci are novel associations with hearing loss of any form, and only one of the ten known hearing loci has a previously reported association with an ARHI-related trait. Gene sets from these loci are enriched in auditory processes such as synaptic activities, nervous system processes, inner ear morphology, and cognition, while genetic correlation analysis revealed strong positive correlations with multiple personality and psychological traits for the first time. Immunohistochemistry for protein localization in adult mouse cochlea implicate metabolic, sensory, and neuronal functions for NID2, CLRN2, and ARHGEF28. These results provide insight into the genetic landscape underlying ARHI, opening up novel therapeutic targets for further investigation. In a wider context, our study also highlights the viability of using self-report phenotypes for genetic discovery in very large samples when deep phenotyping is unavailable.