Detection of chromosomal breakpoints in patients with developmental delay and speech disorders

Identification of a complex intrachromosomal inverted insertion in the long arm of chromosome 9 as a cause of tuberous sclerosis complex in a Korean family

BACKGROUND

Tuberous sclerosis complex (TSC) is an autosomal dominant multisystem disorder, caused by a loss-of-function of either TSC1 or TSC2 gene. However, in 10%-15% TSC patients there is no pathogenic variant identified in either TSC1 or TSC2 genes based on standard clinical testing.

METHODS

In this study, genome sequencing was performed for families with clinical diagnosis of TSC with negative results from TSC1 and TSC2 single-gene tests.

RESULTS

Herein, we report a family presenting a classical TSC phenotype with an unusual, complex structural variant involving the TSC1 gene: an intrachromosomal inverted insertion in the long arm of chromosome 9. We speculate that the inverted 9q33.3q34.13 region was inserted into the q31.2 region with the 3'-end of the breakpoint of the inversion being located within the TSC1 gene, resulting in premature termination of TSC1.

CONCLUSIONS

In this study, we demonstrate the utility of genome sequencing for the identification of complex chromosomal rearrangement. Because the breakpoints are located within the deep intronic/intergenic region, this copy-neutral variant was missed by the TSC1 and TSC2 single-gene tests and contributed to an unknown etiology. Together, this finding suggests that complex structural variants may be underestimated causes for the etiology of TSC.

Hypothermic Protection in Neocortex Is Topographic and Laminar, Seizure Unmitigating, and Partially Rescues Neurons Depleted of RNA Splicing Protein Rbfox3/NeuN in Neonatal Hypoxic-Ischemic Male Piglets

The effects of hypothermia on neonatal encephalopathy may vary topographically and cytopathologically in the neocortex with manifestations potentially influenced by seizures that alter the severity, distribution, and type of neuropathology. We developed a neonatal piglet survival model of hypoxic-ischemic (HI) encephalopathy and hypothermia (HT) with continuous electroencephalography (cEEG) for seizures. Neonatal male piglets received HI-normothermia (NT), HI-HT, sham-NT, or sham-HT treatments. Randomized unmedicated sham and HI piglets underwent cEEG during recovery. Survival was 2-7 days. Normal and pathological neurons were counted in different neocortical areas, identified by cytoarchitecture and connectomics, using hematoxylin and eosin staining and immunohistochemistry for RNA-binding FOX-1 homolog 3 (Rbfox3/NeuN). Seizure burden was determined. HI-NT piglets had a reduced normal/total neuron ratio and increased ischemic-necrotic/total neuron ratio relative to sham-NT and sham-HT piglets with differing severities in the anterior and posterior motor, somatosensory, and frontal cortices. Neocortical neuropathology was attenuated by HT. HT protection was prominent in layer III of the inferior parietal cortex. Rbfox3 immunoreactivity distinguished cortical neurons as: Rbfox3-positive/normal, Rbfox3-positive/ischemic-necrotic, and Rbfox3-depleted. HI piglets had an increased Rbfox3-depleted/total neuron ratio in layers II and III compared to sham-NT piglets. Neuronal Rbfox3 depletion was partly rescued by HT. Seizure burdens in HI-NT and HI-HT piglets were similar. We conclude that the neonatal HI piglet neocortex has: (1) suprasylvian vulnerability to HI and seizures; (2) a limited neuronal cytopathological repertoire in functionally different regions that engages protective mechanisms with HT; (3) higher seizure burden, insensitive to HT, that is correlated with more panlaminar ischemic-necrotic neurons in the somatosensory cortex; and (4) pathological RNA splicing protein nuclear depletion that is sensitive to HT. This work demonstrates that HT protection of the neocortex in neonatal HI is topographic and laminar, seizure unmitigating, and restores neuronal depletion of RNA splicing factor.

17q25.3 copy number changes: association with neurodevelopmental disorders and cardiac malformation

Copy number variants (CNVs) have been identified as common genomic variants that play a significant role in inter-individual variability. Conversely, rare recurrent CNVs have been found to be causal for many disorders with well-established genotype-phenotype relationships. However, the phenotypic implications of rare non-recurrent CNVs remain poorly understood. Herein, we re-investigated 18,542 cases reported from chromosomal microarray at Greenwood Genetic Center from 2010 to 2022 and identified 15 cases with CNVs involving the 17q25.3 region. We report the detailed clinical features of these subjects, and compare with the cases reported in the literature to determine genotype-phenotype correlations for a subset of genes in this region. The CNVs in the 17q25.3 region were found to be rare events, with a prevalence of 0.08% (15/18542) observed in our cohort. The CNVs were dispersed across the entire 17q25.3 region with variable breakpoints and no smallest region of overlap. The subjects presented with a wide range of clinical features, with neurodevelopmental disorders (autism spectrum disorder, intellectual disability, developmental delay) being the most common features (80%), then expressive language disorder (33%), and finally cardiovascular malformations (26%). The association of CNVs involving the critical gene-rich region of 17q25.3 with neurodevelopmental disorders and cardiac malformation, implicates several genes as plausible drivers for these events.

Gene sequencing and result analysis of balanced translocation carriers by third-generation gene sequencing technology

Because the total gene copy number remains constant and all genes are normally expressed, carriers of balanced chromosomal translocations usually have a normal phenotype but are able to produce many different types of gametes during meiosis, and unbalanced gametes lead to increased risks of infertility, recurrent spontaneous abortion, stillbirth, neonatal death or malformations and intellectual abnormalities in offspring. The key to balanced translocations lies in finding the breakpoints, but current genetic testing techniques are all short-read sequencing, with the disadvantage of procedural complexity and imprecision for precisely identifying the breakpoints. The latest third-generation sequencing technology overcomes these drawbacks and uses robust long-read sequencing to accurately and rapidly detect genome-wide information and identify breakpoint locations. In this paper, we performed whole genome long-read sequencing using an Oxford Nanopore sequencer to detect the breakpoints of 4 balanced chromosomal translocation carriers. The results showed that employing about ~ 10× coverage confirmed 6 of the 8 breakpoints, of which, 2 had microdeletions/insertions identified near the breakpoints and 4 had breakpoints that disrupted the normal gene structure and were simultaneously tested for genome-wide structural variation (SV). The results show that whole genome long-read sequencing is an efficient method for pinpointing translocation breakpoints and providing genome-wide information, which is essential for medical genetics and preimplantation genetic testing.

Genome-wide association study of population-standardised cognitive performance phenotypes in a rural South African community

Cognitive function is an indicator for global physical and mental health, and cognitive impairment has been associated with poorer life outcomes and earlier mortality. A standard cognition test, adapted to a rural-dwelling African community, and the Oxford Cognition Screen-Plus were used to capture cognitive performance as five continuous traits (total cognition score, verbal episodic memory, executive function, language, and visuospatial ability) for 2,246 adults in this population of South Africans. A novel common variant, rs73485231, reached genome-wide significance for association with episodic memory using data for ~14 million markers imputed from the H3Africa genotyping array data. Window-based replication of previously implicated variants and regions of interest support the discovery of African-specific associated variants despite the small population size and low allele frequency. This African genome-wide association study identifies suggestive associations with general cognition and domain-specific cognitive pathways and lays the groundwork for further genomic studies on cognition in Africa.

Inversions on human chromosomes

Human chromosome inversions are types of balanced structural variations, making them difficult to analyze. Thanks to PEM (paired-end sequencing and mapping), there has been tremendous progress in studying inversions. Inversions play an important role as an evolutionary factor, contributing to the formation of gonosomes, speciation of chimpanzees and humans, and inv17q21.3 or inv8p23.1 exhibit the features of natural selection. Both inversions have been related to pathogenic phenotype by directly affecting a gene structure (e.g., inv5p15.1q14.1), regulating gene expression (e.g., inv7q21.3q35) and by predisposing to other secondary arrangements (e.g., inv7q11.23). A polymorphism of human inversions is documented by the InvFEST database (a database that stores information about clinical predictions, validations, frequency of inversions, etc.), but only a small fraction of these inversions is validated, and a detailed analysis is complicated by the frequent location of breakpoints within regions of repetitive sequences.

Exploring the Genetic Causality of Discordant Phenotypes in Familial Apparently Balanced Translocation Cases Using Whole Exome Sequencing

Familial apparently balanced translocations (ABTs) are usually not associated with a phenotype; however, rarely, ABTs segregate with discordant phenotypes in family members carrying identical rearrangements. The current study was a follow-up investigation of four familial ABTs, where whole exome sequencing (WES) was implemented as a diagnostic tool to identify the underlying genetic aetiology of the patients' phenotypes. Data were analysed using an in-house bioinformatics pipeline alongside VarSome Clinical. WES findings were validated with Sanger sequencing, while the impact of splicing and missense variants was assessed by reverse-transcription PCR and in silico tools, respectively. Novel candidate variants were identified in three families. In family 1, it was shown that the de novo pathogenic STXBP1 variant (NM_003165.6:c.1110+2T>G) affected splicing and segregated with the patient's phenotype. In family 2, a likely pathogenic TUBA1A variant (NM_006009.4:c.875C>T, NP_006000.2:p.(Thr292Ile)) could explain the patient's symptoms. In family 3, an SCN1A variant of uncertain significance (NM_006920.6:c.5060A>G, NP_008851.3:p.(Glu1687Gly)) required additional evidence to sufficiently support causality. This first report of WES application in familial ABT carriers with discordant phenotypes supported our previous findings describing such rearrangements as coincidental. Thus, WES can be recommended as a complementary test to find the monogenic cause of aberrant phenotypes in familial ABT carriers.

Mild Traumatic Brain Injury-Induced Disruption of the Blood-Brain Barrier Triggers an Atypical Neuronal Response

Mild TBI (mTBI), which affects 75% of TBI survivors or more than 50 million people worldwide each year, can lead to consequences including sleep disturbances, cognitive impairment, mood swings, and post-traumatic epilepsy in a subset of patients. To interrupt the progression of these comorbidities, identifying early pathological events is key. Recent studies have shown that microbleeds, caused by mechanical impact, persist for months after mTBI and are correlated to worse mTBI outcomes. However, the impact of mTBI-induced blood-brain barrier damage on neurons is yet to be revealed. We used a well-characterized mouse model of mTBI that presents with frequent and widespread but size-restricted damage to the blood-brain barrier to assess how neurons respond to exposure of blood-borne factors in this pathological context. We used immunohistochemistry and histology to assess the expression of neuronal proteins in excitatory and inhibitory neurons after mTBI. We observed that the expression of NeuN, Parvalbumin, and CamKII was lost within minutes in areas with blood-brain barrier disruption. Yet, the neurons remained alive and could be detected using a fluorescent Nissl staining even 6 months later. A similar phenotype was observed after exposure of neurons to blood-borne factors due to endothelial cell ablation in the absence of a mechanical impact, suggesting that entrance of blood-borne factors into the brain is sufficient to induce the neuronal atypical response. Changes in postsynaptic spines were observed indicative of functional changes. Thus, this study demonstrates That exposure of neurons to blood-borne factors causes a rapid and sustained loss of neuronal proteins and changes in spine morphology in the absence of neurodegeneration, a finding that is likely relevant to many neuropathologies.

Mitochondrial dysfunction-induced H3K27 hyperacetylation perturbs enhancers in Parkinson's disease

Mitochondrial dysfunction is a major pathophysiological contributor to the progression of Parkinson’s disease (PD); however, whether it contributes to epigenetic dysregulation remains unknown. Here, we show that both chemically and genetically driven mitochondrial dysfunctions share a common mechanism of epigenetic dysregulation. Under both scenarios, lysine 27 acetylation of likely variant H3.3 (H3.3K27ac) increased in dopaminergic neuronal models of PD, thereby opening that region to active enhancer activity via H3K27ac. These vulnerable epigenomic loci represent potential transcription factor motifs for PD pathogenesis. We further confirmed that mitochondrial dysfunction induces H3K27ac in ex vivo and in vivo (MitoPark) neurodegenerative models of PD. Notably, the significantly increased H3K27ac in postmortem PD brains highlights the clinical relevance to the human PD population. Our results reveal an exciting mitochondrial dysfunction-metabolism-H3K27ac-transcriptome axis for PD pathogenesis. Collectively, the mechanistic insights link mitochondrial dysfunction to epigenetic dysregulation in dopaminergic degeneration and offer potential new epigenetic intervention strategies for PD.

Feasibility study of using unbalanced embryos as a reference to distinguish euploid carrier from noncarrier embryos by single nucleotide polymorphism array for reciprocal translocations

OBJECTIVES

To study the feasibility of using unbalanced embryos as a reference in distinguishing euploid carrier and noncarrier embryos by single nucleotide polymorphism (SNP) array-based preimplantation genetic testing (PGT) for reciprocal translocations.

METHODS

After comprehensive chromosome screening (CCS), euploid embryos were identified as normal or carriers using a family member as a reference. Next, unbalanced embryos were used as a reference, and the results were compared with the previous ones. Karyotypes of transferred embryos were validated by prenatal diagnosis.

RESULTS

Of 995 embryos from 110 couples, 288 were found to be euploid. Using a family member as a reference, 142 and 144 embryos were tested to be euploid noncarrier and carrier respectively, and the remaining 2 embryos were undetermined. When unbalanced embryos were selected as references, all the results were consistent with the previous ones. A total of 107 embryos were transferred, resulting in 66 clinical pregnancies. Karyotypes of prenatal diagnosis were all in accordance with the results of tested embryos.

CONCLUSIONS

SNP array-based haplotyping is a rapid and effective way to distinguish between euploid carrier and noncarrier embryos. In case no family member is available as a reference, unbalanced embryos can be used for identification of euploid carrier and noncarrier embryos.

X-linked muscular dystrophy in a Labrador Retriever strain: phenotypic and molecular characterisation

Background

Canine models of Duchenne muscular dystrophy (DMD) are a valuable tool to evaluate potential therapies because they faithfully reproduce the human disease. Several cases of dystrophinopathies have been described in canines, but the Golden Retriever muscular dystrophy (GRMD) model remains the most used in preclinical studies. Here, we report a new spontaneous dystrophinopathy in a Labrador Retriever strain, named Labrador Retriever muscular dystrophy (LRMD).

Methods

A colony of LRMD dogs was established from spontaneous cases. Fourteen LRMD dogs were followed-up and compared to the GRMD standard using several functional tests. The disease causing mutation was studied by several molecular techniques and identified using RNA-sequencing.

Results

The main clinical features of the GRMD disease were found in LRMD dogs; the functional tests provided data roughly overlapping with those measured in GRMD dogs, with similar inter-individual heterogeneity. The LRMD causal mutation was shown to be a 2.2-Mb inversion disrupting the DMD gene within intron 20 and involving the TMEM47 gene. In skeletal muscle, the Dp71 isoform was ectopically expressed, probably as a consequence of the mutation. We found no evidence of polymorphism in either of the two described modifier genes LTBP4 and Jagged1. No differences were found in Pitpna mRNA expression levels that would explain the inter-individual variability.

Conclusions

This study provides a full comparative description of a new spontaneous canine model of dystrophinopathy, found to be phenotypically equivalent to the GRMD model. We report a novel large DNA mutation within the DMD gene and provide evidence that LRMD is a relevant model to pinpoint additional DMD modifier genes.

Location of Balanced Chromosome-Translocation Breakpoints by Long-Read Sequencing on the Oxford Nanopore Platform

Genomic structural variants, including translocations, inversions, insertions, deletions, and duplications, are challenging to be reliably detected by traditional genomic technologies. In particular, balanced translocations and inversions can neither be identified by microarrays since they do not alter chromosome copy numbers, nor by short-read sequencing because of the unmappability of short reads against repetitive genomic regions. The precise localization of breakpoints is vital for exploring genetic causes in patients with balanced translocations or inversions. Long-read sequencing techniques may detect these structural variants in a more direct, efficient, and accurate manner. Here, we performed whole-genome, long-read sequencing using the Oxford Nanopore GridION sequencer to detect breakpoints in six balanced chromosome translocation carriers and one inversion carrier. The results showed that all the breakpoints were consistent with the karyotype results with only ~10× coverage. Polymerase chain reaction (PCR) and Sanger sequencing confirmed 8 out of 14 breakpoints; however, other breakpoint loci were slightly missed since they were either in highly repetitive regions or pericentromeric regions. Some of the breakpoints interrupted normal gene structure, and in other cases, micro-deletions/insertions were found just next to the breakpoints. We also detected haplotypes around the breakpoint regions. Our results suggest that long-read, whole-genome sequencing is an ideal strategy for precisely localizing translocation breakpoints and providing haplotype information, which is essential for medical genetics and preimplantation genetic testing.

Whole genome paired-end sequencing elucidates functional and phenotypic consequences of balanced chromosomal rearrangement in patients with developmental disorders

BACKGROUND

Balanced chromosomal rearrangements associated with abnormal phenotype are rare events, but may be challenging for genetic counselling, since molecular characterisation of breakpoints is not performed routinely. We used next-generation sequencing to characterise breakpoints of balanced chromosomal rearrangements at the molecular level in patients with intellectual disability and/or congenital anomalies.

METHODS

Breakpoints were characterised by a paired-end low depth whole genome sequencing (WGS) strategy and validated by Sanger sequencing. Expression study of disrupted and neighbouring genes was performed by RT-qPCR from blood or lymphoblastoid cell line RNA.

RESULTS

Among the 55 patients included (41 reciprocal translocations, 4 inversions, 2 insertions and 8 complex chromosomal rearrangements), we were able to detect 89% of chromosomal rearrangements (49/55). Molecular signatures at the breakpoints suggested that DNA breaks arose randomly and that there was no major influence of repeated elements. Non-homologous end-joining appeared as the main mechanism of repair (55% of rearrangements). A diagnosis could be established in 22/49 patients (44.8%), 15 by gene disruption (KANSL1, FOXP1, SPRED1, TLK2, MBD5, DMD, AUTS2, MEIS2, MEF2C, NRXN1, NFIX, SYNGAP1, GHR, ZMIZ1) and 7 by position effect (DLX5, MEF2C, BCL11B, SATB2, ZMIZ1). In addition, 16 new candidate genes were identified. Systematic gene expression studies further supported these results. We also showed the contribution of topologically associated domain maps to WGS data interpretation.

CONCLUSION

Paired-end WGS is a valid strategy and may be used for structural variation characterisation in a clinical setting.

Rare Copy Number Variations and Predictors in Children With Intellectual Disability and Epilepsy

Introduction: The concurrence of intellectual disability/global developmental delay and epilepsy (ID/GDD-EP) is very common in the pediatric population. The etiologies for both conditions are complex and largely unknown. The predictors of significant copy number variations (CNVs) are known for the cases with ID/GDD, but unknown for those with exclusive ID/GDD-EP. Importantly, the known predictors are largely from the same ethnic group; hence, they lack replication.

Purpose: We aimed to determine and investigate the diagnostic yield of CNV tests, new causative CNVs, and the independent predictors of significant CNVs in Chinese children with unexplained ID/GDD-EP.

Materials and methods: A total of 100 pediatric patients with unexplained ID/GDD-EP and 1,000 healthy controls were recruited. The American College of Medical Genetics guideline was used to classify the CNVs. Additionally, clinical information was collected and compared between those with significant and non-significant CNVs.

Results: Twenty-eight percent of the patients had significant CNVs, 16% had variants of unknown significance, and 56% had non-significant CNVs. In total, 31 CNVs were identified in 28% (28/100) of cases: 25 pathogenic and 6 likely pathogenic. Eighteen known syndromes were diagnosed in 17 cases. Thirteen rare CNVs (8 novel and 5 reported in literature) were identified, of which three spanned dosage-sensitive genes: 19q13.2 deletion (ATP1A3), Xp11.4-p11.3 deletion (CASK), and 6q25.3-q25.3 deletion (ARID1B). By comparing clinical features in patients with significant CNVs against those with non-significant CNVs, a statistically significant association was found between the presence of significant CNVs and speech and language delay for those aged above 2 years and for those with facial malformations, microcephaly, congenital heart disease, fair skin, eye malformations, and mega cisterna magna. Multivariate logistic regression analysis allowed the identification of two independent significant CNV predictors, which are eye malformations and facial malformations.

Conclusion: Our study supports the performance of CNV tests in pediatric patients with unexplained ID/GDD-EP, as there is high diagnostic yield, which informs genetic counseling. It adds 13 rare CNVs (8 novel), which can be accountable for both conditions. Moreover, congenital eye and facial malformations are clinical markers that can aid clinicians to understand which patients can benefit from the CNV testing and which will not, thus helping patients to avoid unnecessary and expensive tests.

RBFOX3/NeuN is dispensable for visual function

RBFOX3/NeuN is a neuronal splicing regulator involved in neural circuitry balance, as well as neurogenesis and synaptogenesis. Rbfox3 is expressed in neurons; however, in the retina, expression is restricted to cells in the ganglion cell layer and some cells of the inner nuclear layer. Rbfox3 is expressed in a layer-specific manner in the retina, which implies a functional role, however, the role of RBFOX3 in the retina is unknown. Rbfox3 homozygous knockout (Rbfox3-/-) mice exhibit deficits in visual learning; therefore, understanding the role of RBFOX3 in the retina is critical for interpreting behavioral results. We found Rbfox3 expression was developmentally regulated in the retina and specifically expressed in ganglion cells, amacrine cells and horizontal cells of the retina. We demonstrate deletion of Rbfox3 resulted in a reduction in the thickness of the inner plexiform layer of the retina, where synapses are formed. Number of ganglion cells and amacrine cells is normal with loss of Rbfox3. Innervation of retinal ganglion cells into their targeted brain regions is normal in Rbfox3-/- mice. Importantly, Rbfox3-/- mice displayed normal non-image and image forming functions. Taken together, our results suggest RBFOX3 is dispensable for visual function.

Partial monosomy 8p and trisomy 16q in two children with developmental delay detected by array comparative genomic hybridization

Two cases of liveborn unrelated children with developmental delay and overlapping unbalanced translocations der(8)t(8;16)(p23.2;q23.3) and der (8)t(8;16)(p23.1;q23.1), leading to partial monosomy 8p and partial trisomy 16q, are reported in the present study. The first patient was a 10-year-old boy with mild developmental delay and minor congenital anomalies (borderline microcephaly, clinodactyly, hypertelorism, epicanthus, mild systolic murmur and kidney reflux). The second patient was a 3 year-old girl with developmental delay, gross motor milestone delay and dysmorphic features. Array-comparative genomic hybridization analysis revealed that partial chromosome 8p monosomy extended from 8p23.2 to 8pter (4.8 Mb) in Patient 1 and from 8p23.1 to 8pter (9.5 Mb) in Patient 2, and partial chromosome 16 trisomy extended from 16q23.3 to 16qter (5.6 Mb) in Patient 1 and from 16q23.1 to 16qter (11.7 Mb) in Patient 2. The mechanism of appearance of the rearrangement in association with the genes involved and the architecture of the region is discussed.

Personalized genome sequencing coupled with iPSC technology identifies GTDC1 as a gene involved in neurodevelopmental disorders

The cellular and molecular mechanisms underlying neurodevelopmental conditions such as autism spectrum disorders have been studied intensively for decades. The ability to generate patient-specific induced pluripotent stem cells (iPSCs) now offers a novel strategy for modelling human diseases. Recent studies have reported the derivation of iPSCs from patients with neurological disorders. The key challenge remains the demonstration of disease-related phenotypes and the ability to model the disease. Here we report a case study with signs of neurodevelopmental disorders (NDDs) harbouring chromosomal rearrangements that were sequenced using long-insert DNA paired-end tag (DNA-PET) sequencing approach. We identified the disruption of a specific gene, GTDC1. By deriving iPSCs from this patient and differentiating them into neural progenitor cells (NPCs) and neurons we dissected the disease process at the cellular level and observed defects in both NPCs and neuronal cells. We also showed that disruption of GTDC1 expression in wild type human NPCs and neurons showed a similar phenotype as patient's iPSCs. Finally, we utilized a zebrafish model to demonstrate a role for GTDC1 in the development of the central nervous system. Our findings highlight the importance of combining sequencing technologies with the iPSC technology for NDDs modelling that could be applied for personalized medicine.

Genome-wide gene by lead exposure interaction analysis identifies UNC5D as a candidate gene for neurodevelopment

BACKGROUND

Neurodevelopment is a complex process involving both genetic and environmental factors. Prenatal exposure to lead (Pb) has been associated with lower performance on neurodevelopmental tests. Adverse neurodevelopmental outcomes are more frequent and/or more severe when toxic exposures interact with genetic susceptibility.

METHODS

To explore possible loci associated with increased susceptibility to prenatal Pb exposure, we performed a genome-wide gene-environment interaction study (GWIS) in young children from Mexico (n = 390) and Bangladesh (n = 497). Prenatal Pb exposure was estimated by cord blood Pb concentration. Neurodevelopment was assessed using the Bayley Scales of Infant Development.

RESULTS

We identified a locus on chromosome 8, containing UNC5D, and demonstrated evidence of its genome-wide significance with mental composite scores (rs9642758, p _meta = 4.35 × 10^-6). Within this locus, the joint effects of two independent single nucleotide polymorphisms (SNPs, rs9642758 and rs10503970) had a p-value of 4.38 × 10^-9 for mental composite scores. Correlating GWIS results with in vitro transcriptomic profiles identified one common gene, SLC1A5, which is involved in synaptic function, neuronal development, and excitotoxicity. Further analysis revealed interconnected interactions that formed a large network of 52 genes enriched with oxidative stress genes and neurodevelopmental genes.

CONCLUSIONS

Our findings suggest that certain genetic polymorphisms within/near genes relevant to neurodevelopment might modify the toxic effects of Pb exposure via oxidative stress.

Copy number variation in a hospital-based cohort of children with epilepsy

Objective

To evaluate the diagnostic yield of microarray analysis in a hospital‐based cohort of children with seizures and to identify novel candidate genes and susceptibility loci for epilepsy.

Methods

Of all children who presented with their first seizure in the University Medical Center Groningen (January 2000 through May 2013) (n = 1,368), we included 226 (17%) children who underwent microarray analysis before June 2014. All 226 children had a definite diagnosis of epilepsy. All their copy number variants (CNVs) on chromosomes 1–22 and X that contain protein‐coding genes and have a prevalence of <1% in healthy controls were evaluated for their pathogenicity.

Results

Children selected for microarray analysis more often had developmental problems (82% vs. 25%, p < 0.001), facial dysmorphisms (49% vs. 8%, p < 0.001), or behavioral problems (41% vs. 13%, p < 0.001) than children who were not selected. We found known clinically relevant CNVs for epilepsy in 24 of the 226 children (11%). Seventeen of these 24 children had been diagnosed with symptomatic focal epilepsy not otherwise specified (71%) and five with West syndrome (21%). Of these 24 children, many had developmental problems (100%), behavioral problems (54%) or facial dysmorphisms (46%). We further identified five novel CNVs comprising four potential candidate genes for epilepsy: MYT1L, UNC5D, SCN4B, and NRXN3.

Significance

The 11% yield in our hospital‐based cohort underscores the importance of microarray analysis in diagnostic evaluation of children with epilepsy.

Accurate Breakpoint Mapping in Apparently Balanced Translocation Families with Discordant Phenotypes Using Whole Genome Mate-Pair Sequencing

Familial apparently balanced translocations (ABTs) segregating with discordant phenotypes are extremely challenging for interpretation and counseling due to the scarcity of publications and lack of routine techniques for quick investigation. Recently, next generation sequencing has emerged as an efficacious methodology for precise detection of translocation breakpoints. However, studies so far have mainly focused on de novo translocations. The present study focuses specifically on familial cases in order to shed some light to this diagnostic dilemma. Whole-genome mate-pair sequencing (WG-MPS) was applied to map the breakpoints in nine two-way ABT carriers from four families. Translocation breakpoints and patient-specific structural variants were validated by Sanger sequencing and quantitative Real Time PCR, respectively. Identical sequencing patterns and breakpoints were identified in affected and non-affected members carrying the same translocations. PTCD1, ATP5J2-PTCD1, CADPS2, and STPG1 were disrupted by the translocations in three families, rendering them initially as possible disease candidate genes. However, subsequent mutation screening and structural variant analysis did not reveal any pathogenic mutations or unique variants in the affected individuals that could explain the phenotypic differences between carriers of the same translocations. In conclusion, we suggest that NGS-based methods, such as WG-MPS, can be successfully used for detailed mapping of translocation breakpoints, which can also be used in routine clinical investigation of ABT cases. Unlike de novo translocations, no associations were determined here between familial two-way ABTs and the phenotype of the affected members, in which the presence of cryptic imbalances and complex chromosomal rearrangements has been excluded. Future whole-exome or whole-genome sequencing will potentially reveal unidentified mutations in the patients underlying the discordant phenotypes within each family. In addition, larger studies are needed to determine the exact percentage for phenotypic risk in families with ABTs.

The genomic landscape of balanced cytogenetic abnormalities associated with human congenital anomalies

Despite the clinical significance of balanced chromosomal abnormalities (BCAs), their characterization has largely been restricted to cytogenetic resolution. We explored the landscape of BCAs at nucleotide resolution in 273 subjects with a spectrum of congenital anomalies. Whole-genome sequencing revised 93% of karyotypes and demonstrated complexity that was cryptic to karyotyping in 21% of BCAs, highlighting the limitations of conventional cytogenetic approaches. At least 33.9% of BCAs resulted in gene disruption that likely contributed to the developmental phenotype, 5.2% were associated with pathogenic genomic imbalances, and 7.3% disrupted topologically associated domains (TADs) encompassing known syndromic loci. Remarkably, BCA breakpoints in eight subjects altered a single TAD encompassing MEF2C, a known driver of 5q14.3 microdeletion syndrome, resulting in decreased MEF2C expression. We propose that sequence-level resolution dramatically improves prediction of clinical outcomes for balanced rearrangements and provides insight into new pathogenic mechanisms, such as altered regulation due to changes in chromosome topology.

Neuronal Splicing Regulator RBFOX3 (NeuN) Regulates Adult Hippocampal Neurogenesis and Synaptogenesis

Dysfunction of RBFOX3 has been identified in neurodevelopmental disorders such as autism spectrum disorder, cognitive impairments and epilepsy and a causal relationship with these diseases has been previously demonstrated with Rbfox3 homozygous knockout mice. Despite the importance of RBFOX3 during neurodevelopment, the function of RBFOX3 regarding neurogenesis and synaptogenesis remains unclear. To address this critical question, we profiled the developmental expression pattern of Rbfox3 in the brain of wild-type mice and analyzed brain volume, disease-relevant behaviors, neurogenesis, synaptic plasticity, and synaptogenesis in Rbfox3 homozygous knockout mice and their corresponding wild-type counterparts. Here we report that expression of Rbfox3 differs developmentally for distinct brain regions. Moreover, Rbfox3 homozygous knockout mice exhibited cold hyperalgesia and impaired cognitive abilities. Focusing on hippocampal phenotypes, we found Rbfox3 homozygous knockout mice displayed deficits in neurogenesis, which was correlated with cognitive impairments. Furthermore, RBFOX3 regulates the exons of genes with synapse-related function. Synaptic plasticity and density, which are related to cognitive behaviors, were altered in the hippocampal dentate gyrus of Rbfox3 homozygous knockout mice; synaptic plasticity decreased and the density of synapses increased. Taken together, our results demonstrate the important role of RBFOX3 during neural development and maturation. In addition, abnormalities in synaptic structure and function occur in Rbfox3 homozygous knockout mice. Our findings may offer mechanistic explanations for human brain diseases associated with dysfunctional RBFOX3.

Genetic variants in RBFOX3 are associated with sleep latency

Time to fall asleep (sleep latency) is a major determinant of sleep quality. Chronic, long sleep latency is a major characteristic of sleep-onset insomnia and/or delayed sleep phase syndrome. In this study we aimed to discover common polymorphisms that contribute to the genetics of sleep latency. We performed a meta-analysis of genome-wide association studies (GWAS) including 2 572 737 single nucleotide polymorphisms (SNPs) established in seven European cohorts including 4242 individuals. We found a cluster of three highly correlated variants (rs9900428, rs9907432 and rs7211029) in the RNA-binding protein fox-1 homolog 3 gene (RBFOX3) associated with sleep latency (P-values=5.77 × 10(-08), 6.59 × 10(-)(08) and 9.17 × 10(-)(08)). These SNPs were replicated in up to 12 independent populations including 30 377 individuals (P-values=1.5 × 10(-)(02), 7.0 × 10(-)(03) and 2.5 × 10(-)(03); combined meta-analysis P-values=5.5 × 10(-07), 5.4 × 10(-07) and 1.0 × 10(-07)). A functional prediction of RBFOX3 based on co-expression with other genes shows that this gene is predominantly expressed in brain (P-value=1.4 × 10(-316)) and the central nervous system (P-value=7.5 × 10(-)(321)). The predicted function of RBFOX3 based on co-expression analysis with other genes shows that this gene is significantly involved in the release cycle of neurotransmitters including gamma-aminobutyric acid and various monoamines (P-values<2.9 × 10(-11)) that are crucial in triggering the onset of sleep. To conclude, in this first large-scale GWAS of sleep latency we report a novel association of variants in RBFOX3 gene. Further, a functional prediction of RBFOX3 supports the involvement of RBFOX3 with sleep latency.

Organ-specific adaptive signaling pathway activation in metastatic breast cancer cells

Breast cancer metastasizes to bone, visceral organs, and/or brain depending on the subtype, which may involve activation of a host organ-specific signaling network in metastatic cells. To test this possibility, we determined gene expression patterns in MDA-MB-231 cells and its mammary fat pad tumor (TMD-231), lung-metastasis (LMD-231), bone-metastasis (BMD-231), adrenal-metastasis (ADMD-231) and brain-metastasis (231-BR) variants. When gene expression between metastases was compared, 231-BR cells showed the highest gene expression difference followed by ADMD-231, LMD-231, and BMD-231 cells. Neuronal transmembrane proteins SLITRK2, TMEM47, and LYPD1 were specifically overexpressed in 231-BR cells. Pathway-analyses revealed activation of signaling networks that would enable cancer cells to adapt to organs of metastasis such as drug detoxification/oxidative stress response/semaphorin neuronal pathway in 231-BR, Notch/orphan nuclear receptor signals involved in steroidogenesis in ADMD-231, acute phase response in LMD-231, and cytokine/hematopoietic stem cell signaling in BMD-231 cells. Only NF-κB signaling pathway activation was common to all except BMD-231 cells. We confirmed NF-κB activation in 231-BR and in a brain metastatic variant of 4T1 cells (4T1-BR). Dimethylaminoparthenolide inhibited NF-κB activity, LYPD1 expression, and proliferation of 231-BR and 4T1-BR cells. Thus, transcriptome change enabling adaptation to host organs is likely one of the mechanisms associated with organ-specific metastasis and could potentially be targeted therapeutically.

Chronic HIV-1 Tat and HIV reduce Rbfox3/NeuN: evidence for sex-related effects

The NeuN antibody has been widely used to identify and quantify neurons in normal and disease situations based on binding to a nuclear epitope in most types of neurons. This epitope was recently identified as the RNA-binding, feminizing locus on X-3 (Rbfox3), a member of the larger, mammalian Fox1 family of RNA binding proteins. Fox1 proteins recognize a unique UGCAUG mRNA motif and regulate alternative splicing of precursor mRNA to control post-transcriptional events important in neuronal differentiation and central nervous system development. Recent clinical findings show that Rbfox3/NeuN gene dosage is altered in certain human neurodevelopmental disorders, and redistribution has been noted in HIV(+) tissue. We hypothesized that HIV-1 Tat might affect Rbfox3/NeuN expression, and examined this question in vivo using inducible transgenic mice, and in vitro using human mesencephalic-derived neurons. Rbfox3/NeuN expression and localization in HIV+ basal ganglia and hippocampus was also examined. Chronic Tat exposure reduced Rbfox3/NeuN protein levels and increased cytoplasmic localization, similar to the effect of HIV exposure. Cytoplasmic Rbfox3/NeuN signal has occasionally been reported, although the meaning or function of cytoplasmic versus nuclear localization remains speculative. Importantly, Rbfox3/NeuN reductions were more significant in male mice. Although Rbfox3/NeuN-expressing cells were significantly decreased by Tat exposure, stereology showed that Nissl(+) neuron numbers remained normal. Thus, loss of Rbfox3/NeuN may relate more to functional change than to neuron loss. The effects of Tat by itself are highly relevant to HIV(+) individuals maintained on antiretroviral therapy, since Tat is released from infected cells even when viral replication is inhibited.

RBFOX3/NeuN is Required for Hippocampal Circuit Balance and Function

RBFOX3 mutations are linked to epilepsy and cognitive impairments, but the underlying pathophysiology of these disorders is poorly understood. Here we report replication of human symptoms in a mouse model with disrupted Rbfox3. Rbfox3 knockout mice displayed increased seizure susceptibility and decreased anxiety-related behaviors. Focusing on hippocampal phenotypes, we found Rbfox3 knockout mice showed increased expression of plasticity genes Egr4 and Arc, and the synaptic transmission and plasticity were defective in the mutant perforant pathway. The mutant dentate granules cells exhibited an increased frequency, but normal amplitude, of excitatory synaptic events, and this change was associated with an increase in the neurotransmitter release probability and dendritic spine density. Together, our results demonstrate anatomical and functional abnormality in Rbfox3 knockout mice, and may provide mechanistic insights for RBFOX3-related human brain disorders.

Identification of disrupted AUTS2 and EPHA6 genes by array painting in a patient carrying a de novo balanced translocation t(3;7) with intellectual disability and neurodevelopment disorder

Intellectual disability (ID) is a frequent feature but is highly clinically and genetically heterogeneous. The establishment of the precise diagnosis in patients with ID is challenging due to this heterogeneity but crucial for genetic counseling and appropriate care for the patients. Among the etiologies of patients with ID, apparently balanced de novo rearrangements represent 0.6%. Several mechanisms explain the ID in patients with apparently balanced de novo rearrangement. Among them, disruption of a disease gene at the breakpoint, is frequently evoked. In this context, technologies recently developed are used to characterize precisely such chromosomal rearrangements. Here, we report the case of a boy with ID, facial features and autistic behavior who is carrying a de novo balanced reciprocal translocation t(3;7)(q11.2;q11.22)dn. Using microarray analysis, array painting (AP) technology combined with molecular study, we have identified the interruption of the autism susceptibility candidate 2 gene (AUTS2) and EPH receptor A6 gene (EPHA6). We consider that the disruption of AUTS2 explains the phenotype of the patient; the exact role of EPHA6 in human pathology is not well defined. Based on the observation of recurrent germinal and somatic translocations involving AUTS2 and the molecular environment content, we put forward the hypothesis that the likely chromosomal mechanism responsible for the translocation could be due either to replicative stress or to recombination-based mechanisms.

Human Structural Variation: Mechanisms of Chromosome Rearrangements

Chromosome structural variation (SV) is a normal part of variation in the human genome, but some classes of SV can cause neurodevelopmental disorders. Analysis of the DNA sequence at SV breakpoints can reveal mutational mechanisms and risk factors for chromosome rearrangement. Large-scale SV breakpoint studies have become possible recently owing to advances in next-generation sequencing (NGS) including whole-genome sequencing (WGS). These findings have shed light on complex forms of SV such as triplications, inverted duplications, insertional translocations, and chromothripsis. Sequence-level breakpoint data resolve SV structure and determine how genes are disrupted, fused, and/or misregulated by breakpoints. Recent improvements in breakpoint sequencing have also revealed non-allelic homologous recombination (NAHR) between paralogous long interspersed nuclear element (LINE) or human endogenous retrovirus (HERV) repeats as a cause of deletions, duplications, and translocations. This review covers the genomic organization of simple and complex constitutional SVs, as well as the molecular mechanisms of their formation.

Human inversions and their functional consequences

Polymorphic inversions are a type of structural variants that are difficult to analyze owing to their balanced nature and the location of breakpoints within complex repeated regions. So far, only a handful of inversions have been studied in detail in humans and current knowledge about their possible functional effects is still limited. However, inversions have been related to phenotypic changes and adaptation in multiple species. In this review, we summarize the evidences of the functional impact of inversions in the human genome. First, given that inversions have been shown to inhibit recombination in heterokaryotes, chromosomes displaying different orientation are expected to evolve independently and this may lead to distinct gene-expression patterns. Second, inversions have a role as disease-causing mutations both by directly affecting gene structure or regulation in different ways, and by predisposing to other secondary arrangements in the offspring of inversion carriers. Finally, several inversions show signals of being selected during human evolution. These findings illustrate the potential of inversions to have phenotypic consequences also in humans and emphasize the importance of their inclusion in genome-wide association studies.

Escape from X inactivation varies in mouse tissues

X chromosome inactivation (XCI) silences most genes on one X chromosome in female mammals, but some genes escape XCI. To identify escape genes in vivo and to explore molecular mechanisms that regulate this process we analyzed the allele-specific expression and chromatin structure of X-linked genes in mouse tissues and cells with skewed XCI and distinguishable alleles based on single nucleotide polymorphisms. Using a binomial model to assess allelic expression, we demonstrate a continuum between complete silencing and expression from the inactive X (Xi). The validity of the RNA-seq approach was verified using RT-PCR with species-specific primers or Sanger sequencing. Both common escape genes and genes with significant differences in XCI status between tissues were identified. Such genes may be candidates for tissue-specific sex differences. Overall, few genes (3-7%) escape XCI in any of the mouse tissues examined, suggesting stringent silencing and escape controls. In contrast, an in vitro system represented by the embryonic-kidney-derived Patski cell line showed a higher density of escape genes (21%), representing both kidney-specific escape genes and cell-line specific escape genes. Allele-specific RNA polymerase II occupancy and DNase I hypersensitivity at the promoter of genes on the Xi correlated well with levels of escape, consistent with an open chromatin structure at escape genes. Allele-specific CTCF binding on the Xi clustered at escape genes and was denser in brain compared to the Patski cell line, possibly contributing to a more compartmentalized structure of the Xi and fewer escape genes in brain compared to the cell line where larger domains of escape were observed.

Novel Insights into NeuN: from Neuronal Marker to Splicing Regulator

Neuronal nuclei (NeuN) is a well-recognized "marker" that is detected exclusively in post-mitotic neurons and was initially identified through an immunological screen to produce neuron-specific antibodies. Immunostaining evidence indicates that NeuN is distributed in the nuclei of mature neurons in nearly all parts of the vertebrate nervous system. NeuN is highly conserved among species and is stably expressed during specific stages of development. Therefore, NeuN has been considered to be a reliable marker of mature neurons for the past two decades. However, this role has been challenged by recent studies indicating that NeuN staining is variable and even absent during certain diseases and specific physiological states. More importantly, despite the widespread use of the anti-NeuN antibody, the natural identity of the NeuN protein remained elusive for 17 years. NeuN was recently eventually identified as an epitope of Rbfox3, which is a novel member of the Rbfox1 family of splicing factors. This identification might provide a novel perspective on NeuN expression during both physiological and pathological conditions. This review summarizes the current progress on the biochemical identity and biological significance of NeuN and recommends caution when applying NeuN immunoreactivity as a definitive marker of mature neurons in certain diseases and specific physiological states.

Deciphering the pathogenic consequences of chromosomal aberrations in human genetic disease

Chromosomal aberrations include translocations, deletions, duplications, inversions, aneuploidies and complex rearrangements. They underlie genetic disease in roughly 15% of patients with multiple congenital abnormalities and/or mental retardation (MCA/MR). In genetic diagnostics, the pathogenicity of chromosomal aberrations in these patients is typically assessed based on criteria such as phenotypic similarity to other patients with the same or overlapping aberration, absence in healthy individuals, de novo occurrence, and protein coding gene content. However, a thorough understanding of the molecular mechanisms that lead to MCA/MR as a result of chromosome aberrations is often lacking. Chromosome aberrations can affect one or more genes in a complex manner, such as by changing the regulation of gene expression, by disrupting exons, and by creating fusion genes. The precise delineation of breakpoints by whole-genome sequencing enables the construction of local genomic architecture and facilitates the prediction of the molecular determinants of the patient's phenotype. Here, we review current methods for breakpoint identification and their impact on the interpretation of chromosome aberrations in patients with MCA/MR. In addition, we discuss opportunities to dissect disease mechanisms based on large-scale genomic technologies and studies in model organisms.

Precise detection of chromosomal translocation or inversion breakpoints by whole-genome sequencing

Structural variations (SVs), including translocations, inversions, deletions and duplications, are potentially associated with Mendelian diseases and contiguous gene syndromes. Determination of SV-related breakpoints at the nucleotide level is important to reveal the genetic causes for diseases. Whole-genome sequencing (WGS) by next-generation sequencers is expected to determine structural abnormalities more directly and efficiently than conventional methods. In this study, 14 SVs (9 balanced translocations, 1 inversion and 4 microdeletions) in 9 patients were analyzed by WGS with a shallow (5 × ) to moderate read coverage (20 × ). Among 28 breakpoints (as each SV has two breakpoints), 19 SV breakpoints had been determined previously at the nucleotide level by any other methods and 9 were uncharacterized. BreakDancer and Integrative Genomics Viewer determined 20 breakpoints (16 translocation, 2 inversion and 2 deletion breakpoints), but did not detect 8 breakpoints (2 translocation and 6 deletion breakpoints). These data indicate the efficacy of WGS for the precise determination of translocation and inversion breakpoints.