Mind the gap: the relevance of the genome reference to resolve rare and pathogenic inversions.

Both long-read genome sequencing (lrGS) and the recently published Telomere to Telomere (T2T) reference genome provide increased coverage and resolution across repetitive regions promising heightened structural variant detection and improved mapping. Inversions (INV), intrachromosomal segments which are rotated 180° and inserted back into the same chromosome, are a class of structural variants particularly challenging to detect due to their copy-number neutral state and association with repetitive regions. Inversions represent about 1/20 of all balanced structural chromosome aberrations and can lead to disease by gene disruption or altering regulatory regions of dosage sensitive genes in cis . Here we remapped the genome data from six individuals carrying unsolved cytogenetically detected inversions. An INV6 and INV10 were resolved using GRCh38 and T2T-CHM13. Finally, an INV9 required optical genome mapping, de novo assembly of lrGS data and T2T-CHM13. This inversion disrupted intron 25 of EHMT1, confirming a diagnosis of Kleefstra syndrome 1 (MIM#610253). These three inversions, only mappable in specific references, prompted us to investigate the presence and population frequencies of differential reference regions (DRRs) between T2T-CHM13, GRCh37, GRCh38, the chimpanzee and bonobo, and hundreds of megabases of DRRs were identified. Our results emphasize the significance of the chosen reference genome and the added benefits of lrGS and optical genome mapping in solving rearrangements in challenging regions of the genome. This is particularly important for inversions and may impact clinical diagnostics.

Long-read sequencing and optical mapping generates near T2T assemblies that resolves a centromeric translocation.

Long-read genome sequencing (lrGS) is a promising method in genetic diagnostics. Here we investigate the potential of lrGS to detect a disease-associated chromosomal translocation between 17p13 and the 19 centromere. We constructed two sets of phased and non-phased de novo assemblies; (i) based on lrGS only and (ii) hybrid assemblies combining lrGS with optical mapping using lrGS reads with a median coverage of 34X. Variant calling detected both structural variants (SVs) and small variants and the accuracy of the small variant calling was compared with those called with short-read genome sequencing (srGS). The de novo and hybrid assemblies had high quality and contiguity with N50 of 62.85 Mb, enabling a near telomere to telomere assembly with less than a 100 contigs per haplotype. Notably, we successfully identified the centromeric breakpoint of the translocation. A concordance of 92% was observed when comparing small variant calling between srGS and lrGS. In summary, our findings underscore the remarkable potential of lrGS as a comprehensive and accurate solution for the analysis of SVs and small variants. Thus, lrGS could replace a large battery of genetic tests that were used for the diagnosis of a single symptomatic translocation carrier, highlighting the potential of lrGS in the realm of digital karyotyping.

Rare coding variants in NOX4 link high ROS levels to psoriatic arthritis mutilans.

Psoriatic arthritis mutilans (PAM) is the rarest and most severe form of psoriatic arthritis, characterized by erosions of the small joints and osteolysis leading to joint disruption. Despite its severity, the underlying mechanisms are unknown, and no susceptibility genes have hitherto been identified. We aimed to investigate the genetic basis of PAM by performing massive parallel sequencing in sixty-one patients from the PAM Nordic cohort. We found rare variants in the NADPH oxidase 4 (NOX4) in four patients. In silico predictions show that the identified variants are potentially damaging. NOXs are the only enzymes producing reactive oxygen species (ROS). NOX4 is specifically involved in the differentiation of osteoclasts, the cells implicated in bone resorption. Functional follow-up studies using cell culture, zebrafish models, and measurement of ROS in patients uncovered that these NOX4 variants increase ROS levels both in vitro and in vivo. We propose NOX4 as the first candidate susceptibility gene for PAM. Our study links high levels of ROS caused by NOX4 variants to the development of PAM, offering a potential therapeutic target.

Feasibility to use whole-genome sequencing as a sole diagnostic method to detect genomic aberrations in pediatric B-cell acute lymphoblastic leukemia.

Introduction: The suitability of whole-genome sequencing (WGS) as the sole method to detect clinically relevant genomic aberrations in B-cell acute lymphoblastic leukemia (ALL) was investigated with the aim of replacing current diagnostic methods. Methods: For this purpose, we assessed the analytical performance of 150 bp paired-end WGS (90x leukemia/30x germline). A set of 88 retrospective B-cell ALL samples were selected to represent established ALL subgroups as well as ALL lacking stratifying markers by standard-of-care (SoC), so-called B-other ALL. Results: Both the analysis of paired leukemia/germline (L/N)(n=64) as well as leukemia-only (L-only)(n=88) detected all types of aberrations mandatory in the current ALLTogether trial protocol, i.e., aneuploidies, structural variants, and focal copy-number aberrations. Moreover, comparison to SoC revealed 100% concordance and that all patients had been assigned to the correct genetic subgroup using both approaches. Notably, WGS could allocate 35 out of 39 B-other ALL samples to one of the emerging genetic subgroups considered in the most recent classifications of ALL. We further investigated the impact of high (90x; n=58) vs low (30x; n=30) coverage on the diagnostic yield and observed an equally perfect concordance with SoC; low coverage detected all relevant lesions. Discussion: The filtration of the WGS findings with a short list of genes recurrently rearranged in ALL was instrumental to extract the clinically relevant information efficiently. Nonetheless, the detection of DUX4 rearrangements required an additional customized analysis, due to multiple copies of this gene embedded in the highly repetitive D4Z4 region. We conclude that the diagnostic performance of WGS as the standalone method was remarkable and allowed detection of all clinically relevant genomic events in the diagnostic setting of B-cell ALL.

Insights into cellular behavior and micromolecular communication in urothelial micrografts.

Autologous micrografting is a technique currently applied within skin wound healing, however, the potential use for surgical correction of other organs with epithelial lining, including the urinary bladder, remains largely unexplored. Currently, little is known about the micrograft expansion potential and the micromolecular events that occur in micrografted urothelial cells. In this study, we aimed to evaluate the proliferative potential of different porcine urothelial micrograft sizes in vitro, and, furthermore, to explore how urothelial micrografts communicate and which microcellular events are triggered. We demonstrated that increased tissue fragmentation subsequently potentiated the yield of proliferative cells and the cellular expansion potential, which confirms, that the micrografting principles of skin epithelium also apply to uroepithelium. Furthermore, we targeted the expression of the extracellular signal-regulated kinase (ERK) pathway and demonstrated that ERK activation occurred predominately at the micrograft borders and that ERK inhibition led to decreased urothelial migration and proliferation. Finally, we successfully isolated extracellular vesicles from the micrograft culture medium and evaluated their contents and relevance within various enriched biological processes. Our findings substantiate the potential of applying urothelial micrografting in future tissue-engineering models for reconstructive urological surgery, and, furthermore, highlights certain mechanisms as potential targets for future wound healing treatments.

Long-read whole-genome analysis of human single cells.

Long-read sequencing has dramatically increased our understanding of human genome variation. Here, we demonstrate that long-read technology can give new insights into the genomic architecture of individual cells. Clonally expanded CD8+ T-cells from a human donor were subjected to droplet-based multiple displacement amplification (dMDA) to generate long molecules with reduced bias. PacBio sequencing generated up to 40% genome coverage per single-cell, enabling detection of single nucleotide variants (SNVs), structural variants (SVs), and tandem repeats, also in regions inaccessible by short reads. 28 somatic SNVs were detected, including one case of mitochondrial heteroplasmy. 5473 high-confidence SVs/cell were discovered, a sixteen-fold increase compared to Illumina-based results from clonally related cells. Single-cell de novo assembly generated a genome size of up to 598 Mb and 1762 (12.8%) complete gene models. In summary, our work shows the promise of long-read sequencing toward characterization of the full spectrum of genetic variation in single cells.

Case report: Extending the spectrum of clinical and molecular findings in FOXC1 haploinsufficiency syndrome.

FOXC1 is a ubiquitously expressed forkhead transcription factor that plays a critical role during early development. Germline pathogenic variants in FOXC1 are associated with anterior segment dysgenesis and Axenfeld-Rieger syndrome (ARS, #602482), an autosomal dominant condition with ophthalmologic anterior segment abnormalities, high risk for glaucoma and extraocular findings including distinctive facial features, as well as dental, skeletal, audiologic, and cardiac anomalies. De Hauwere syndrome is an ultrarare condition previously associated with 6p microdeletions and characterized by anterior segment dysgenesis, joint instability, short stature, hydrocephalus, and skeletal abnormalities. Here, we report clinical findings of two unrelated adult females with FOXC1 haploinsufficiency who have ARS and skeletal abnormalities. Final molecular diagnoses of both patients were achieved using genome sequencing. Patient 1 had a complex rearrangement involving a 4.9 kB deletion including FOXC1 coding region (Hg19; chr6:1,609,721-1,614,709), as well as a 7 MB inversion (Hg19; chr6:1,614,710-8,676,899) and a second deletion of 7.1 kb (Hg19; chr6:8,676,900-8,684,071). Patient 2 had a heterozygous single nucleotide deletion, resulting in a frameshift and a premature stop codon in FOXC1 (NM_001453.3): c.467del, p.(Pro156Argfs*25). Both individuals had moderate short stature, skeletal abnormalities, anterior segment dysgenesis, glaucoma, joint laxity, pes planovalgus, dental anomalies, hydrocephalus, distinctive facial features, and normal intelligence. Skeletal surveys revealed dolichospondyly, epiphyseal hypoplasia of femoral and humeral heads, dolichocephaly with frontal bossin gand gracile long bones. We conclude that haploinsufficiency of FOXC1 causes ARS and a broad spectrum of symptoms with variable expressivity that at its most severe end also includes a phenotype overlapping with De Hauwere syndrome.

Transposable element insertions in 1000 Swedish individuals.

The majority of rare diseases are genetic, and regardless of advanced high-throughput genomics-based investigations, 60% of patients remain undiagnosed. A major factor limiting our ability to identify disease-causing alterations is a poor understanding of the morbid and normal human genome. A major genomic contributor of which function and distribution remain largely unstudied are the transposable elements (TE), which constitute 50% of our genome. Here we aim to resolve this knowledge gap and increase the diagnostic yield of rare disease patients investigated with clinical genome sequencing. To this end we characterized TE insertions in 1000 Swedish individuals from the SweGen dataset and 2504 individuals from the 1000 Genomes Project (1KGP), creating seven population-specific TE insertion databases. Of note, 66% of TE insertions in SweGen were present at >1% in the 1KGP databases, proving that most insertions are common across populations. Focusing on the rare TE insertions, we show that even though ~0.7% of those insertions affect protein coding genes, they rarely affect known disease casing genes (<0.1%). Finally, we applied a TE insertion identification workflow on two clinical cases where disease causing TE insertions were suspected and could verify the presence of pathogenic TE insertions in both. Altogether we demonstrate the importance of TE insertion detection and highlight possible clinical implications in rare disease diagnostics.

Discovery of non-reference processed pseudogenes in the Swedish population.

The vast majority of the human genome is non-coding. There is a diversity of non-coding features, some of which have functional importance. Although the non-coding regions constitute the majority of the genome, they remain understudied, and for a long time, these regions have been referred to as junk DNA. Pseudogenes are one of these features. A pseudogene is a non-functional copy of a protein-coding gene. Pseudogenes may arise through a variety of genetic mechanisms. Processed pseudogenes are formed through reverse transcription of mRNA by LINE elements, after which the cDNA is integrated into the genome. Processed pseudogenes are known to be variable across populations; however, the variability and distribution remains unknown. Herein, we apply a custom-designed processed pseudogene pipeline on the whole genome sequencing data of 3,500 individuals; 2,500 individuals from the thousand genomes dataset, as well as 1,000 Swedish individuals. Through these analyses, we discover over 3,000 pseudogenes missing from the GRCh38 reference. Utilising our pipeline, we position 74% of the detected processed pseudogenes-allowing for analyses of formation. Notably, we find that common structural variant callers, such as Delly, classify the processed pseudogenes as deletion events, which are later predicted to be truncating variants. By compiling lists of non-reference processed pseudogenes and their frequencies, we find a great variability of pseudogenes; indicating that non-reference processed pseudogenes may be useful for DNA testing and as population-specific markers. In summary, our findings highlight a great diversity of processed pseudogenes, that processed pseudogenes are actively formed in the human genome; and that our pipeline may be used to reduce false positive structural variation caused by the misalignment and subsequent misclassification of non-reference processed pseudogenes.

Genome sequencing is a sensitive first-line test to diagnose individuals with intellectual disability.

PURPOSE: Individuals with intellectual disability (ID) and/or neurodevelopment disorders (NDDs) are currently investigated with several different approaches in clinical genetic diagnostics. METHODS: We compared the results from 3 diagnostic pipelines in patients with ID/NDD: genome sequencing (GS) first (N = 100), GS as a secondary test (N = 129), or chromosomal microarray (CMA) with or without FMR1 analysis (N = 421). RESULTS: The diagnostic yield was 35% (GS-first), 26% (GS as a secondary test), and 11% (CMA/FMR1). Notably, the age of diagnosis was delayed by 1 year when GS was performed as a secondary test and the cost per diagnosed individual was 36% lower with GS first than with CMA/FMR1. Furthermore, 91% of those with a negative result after CMA/FMR1 analysis (338 individuals) have not yet been referred for additional genetic testing and remain undiagnosed. CONCLUSION: Our findings strongly suggest that genome analysis outperforms other testing strategies and should replace traditional CMA and FMR1 analysis as a first-line genetic test in individuals with ID/NDD. GS is a sensitive, time- and cost-effective method that results in a confirmed molecular diagnosis in 35% of all referred patients.

Mosaic Deletions of Known Genes Explain Skeletal Dysplasias With High and Low Bone Mass.

Mosaicism, a state in which an individual has two or more genetically distinct populations of cells in the body, can be difficult to detect because of either mild or atypical clinical presentation and limitations in the commonly used detection methods. Knowledge of the role of mosaicism is limited in many skeletal disorders, including osteopathia striata with cranial sclerosis (OSCS) and cleidocranial dysplasia (CCD). We used whole-genome sequencing (WGS) with coverage >40× to identify the genetic causes of disease in two clinically diagnosed patients. In a female patient with OSCS, we identified a mosaic 7-nucleotide frameshift deletion in exon 2 of AMER1, NM_152424.4:c.855_861del:p.(His285Glnfs*7), affecting 8.3% of the WGS reads. In a male patient with CCD, approximately 34% of the WGS reads harbored a 3710-basepair mosaic deletion, NC_000006.11:g.45514471_45518181del, starting in intron 8 of RUNX2 and terminating in the 3' untranslated region. Droplet digital polymerase chain reaction was used to validate these deletions and quantify the absolute level of mosaicism in each patient. Although constitutional variants in AMER1 and RUNX2 are a known cause of OSCS and CCD, respectively, the mosaic changes here reported have not been described previously. Our study indicates that mosaicism should be considered in unsolved cases of skeletal dysplasia and should be investigated with comprehensive and sensitive detection methods. © 2022 The Authors. JBMR Plus published by Wiley Periodicals LLC on behalf of American Society for Bone and Mineral Research.

Multi-Omic Investigations of a 17-19 Translocation Links MINK1 Disruption to Autism, Epilepsy and Osteoporosis.

Balanced structural variants, such as reciprocal translocations, are sometimes hard to detect with sequencing, especially when the breakpoints are located in repetitive or insufficiently mapped regions of the genome. In such cases, long-range information is required to resolve the rearrangement, identify disrupted genes and, in symptomatic carriers, pinpoint the disease-causing mechanisms. Here, we report an individual with autism, epilepsy and osteoporosis and a de novo balanced reciprocal translocation: t(17;19) (p13;p11). The genomic DNA was analyzed by short-, linked- and long-read genome sequencing, as well as optical mapping. Transcriptional consequences were assessed by transcriptome sequencing of patient-specific neuroepithelial stem cells derived from induced pluripotent stem cells (iPSC). The translocation breakpoints were only detected by long-read sequencing, the first on 17p13, located between exon 1 and exon 2 of MINK1 (Misshapen-like kinase 1), and the second in the chromosome 19 centromere. Functional validation in induced neural cells showed that MINK1 expression was reduced by >50% in the patient's cells compared to healthy control cells. Furthermore, pathway analysis revealed an enrichment of changed neural pathways in the patient's cells. Altogether, our multi-omics experiments highlight MINK1 as a candidate monogenic disease gene and show the advantages of long-read genome sequencing in capturing centromeric translocations.

Multi-omics analysis reveals multiple mechanisms causing Prader-Willi like syndrome in a family with a X;15 translocation.

Prader-Willi syndrome (PWS; MIM# 176270) is a neurodevelopmental disorder caused by the loss of expression of paternally imprinted genes within the PWS region located on 15q11.2. It is usually caused by either maternal uniparental disomy of chromosome 15 (UPD15) or 15q11.2 recurrent deletion(s). Here, we report a healthy carrier of a balanced X;15 translocation and her two daughters, both with the karyotype 45,X,der(X)t(X;15)(p22;q11.2),-15. Both daughters display symptoms consistent with haploinsufficiency of the SHOX gene and PWS. We explored the architecture of the derivative chromosomes and investigated effects on gene expression in patient-derived neural cells. First, a multiplex ligation-dependent probe amplification methylation assay was used to determine the methylation status of the PWS-region revealing maternal UPD15 in daughter 2, explaining her clinical symptoms. Next, short read whole genome sequencing and 10X genomics linked read sequencing was used to pinpoint the exact breakpoints of the translocation. Finally, we performed transcriptome sequencing on neuroepithelial stem cells from the mother and from daughter 1 and observed biallelic expression of genes in the PWS region (including SNRPN) in daughter 1. In summary, our multi-omics analysis highlights two different PWS mechanisms in one family and provide an example of how structural variation can affect imprinting through long-range interactions.

Linked-read whole-genome sequencing resolves common and private structural variants in multiple myeloma.

Multiple myeloma (MM) is an incurable and aggressive plasma cell malignancy characterized by a complex karyotype with multiple structural variants (SVs) and copy-number variations (CNVs). Linked-read whole-genome sequencing (lrWGS) allows for refined detection and reconstruction of SVs by providing long-range genetic information from standard short-read sequencing. This makes lrWGS an attractive solution for capturing the full genomic complexity of MM. Here we show that high-quality lrWGS data can be generated from low numbers of cells subjected to fluorescence-activated cell sorting (FACS) without DNA purification. Using this protocol, we analyzed MM cells after FACS from 37 patients with MM using lrWGS. We found high concordance between lrWGS and fluorescence in situ hybridization (FISH) for the detection of recurrent translocations and CNVs. Outside of the regions investigated by FISH, we identified >150 additional SVs and CNVs across the cohort. Analysis of the lrWGS data allowed for resolution of the structure of diverse SVs affecting the MYC and t(11;14) loci, causing the duplication of genes and gene regulatory elements. In addition, we identified private SVs causing the dysregulation of genes recurrently involved in translocations with the IGH locus and show that these can alter the molecular classification of MM. Overall, we conclude that lrWGS allows for the detection of aberrations critical for MM prognostics and provides a feasible route for providing comprehensive genetics. Implementing lrWGS could provide more accurate clinical prognostics, facilitate genomic medicine initiatives, and greatly improve the stratification of patients included in clinical trials.

A somatic UBA2 variant preceded ETV6-RUNX1 in the concordant BCP-ALL of monozygotic twins.

Genetic analysis of leukemic clones in monozygotic twins with concordant acute lymphoblastic leukemia (ALL) has proved a unique opportunity to gain insight into the molecular phylogenetics of leukemogenesis. Using whole-genome sequencing, we characterized constitutional and somatic single nucleotide variants/insertion-deletions (indels) and structural variants in a monozygotic twin pair with concordant ETV6-RUNX1+ B-cell precursor ALL (BCP-ALL). In addition, digital PCR (dPCR) was applied to evaluate the presence of and quantify selected somatic variants at birth, diagnosis, and remission. A shared somatic complex rearrangement involving chromosomes 11, 12, and 21 with identical fusion sequences in leukemias of both twins offered direct proof of a common clonal origin. The ETV6-RUNX1 fusion detected at diagnosis was found to originate from this complex rearrangement. A shared somatic frameshift deletion in UBA2 was also identified in diagnostic samples. In addition, each leukemia independently acquired analogous deletions of 3 genes recurrently targeted in BCP-ALLs (ETV6, ATF7IP, and RAG1/RAG2), providing evidence of a convergent clonal evolution only explained by a strong concurrent selective pressure. Quantification of the UBA2 deletion by dPCR surprisingly indicated it persisted in remission. This, for the first time to our knowledge, provided evidence of a UBA2 variant preceding the well-established initiating event ETV6-RUNX1. Further, we suggest the UBA2 deletion exerted a leukemia predisposing effect and that its essential role in Small Ubiquitin-like Modifier (SUMO) attachment (SUMOylation), regulating nearly all physiological and pathological cellular processes such as DNA-repair by nonhomologous end joining, may hold a mechanistic explanation for the predisposition.

Genomic profile - a possible diagnostic and prognostic marker in upper tract urothelial carcinoma.

OBJECTIVES: To investigate gene alterations as diagnostic and prognostic markers in upper tract urothelial carcinoma (UTUC). PATIENTS AND METHODS: Patients with UTUC who underwent nephroureterectomy between 2005 and 2012 were followed until November 2020. DNA was extracted from paraffin-embedded tumour tissue. Next-generation sequencing using a 388-gene panel was performed. First a blinded analysis using principal component analysis and hierarchical clustering was used to search for patterns of mutations. Then a comparative analysis using analysis of variance (ANOVA) was used to search for mutations enriched in groups of various grades, stages, and survival. In addition, careful manual annotation was used to identify pathogenic mutations over-represented in tumours of high grade/stage and/or poor survival. RESULTS: A total of 39 patients were included. All tumour stages and grades were represented in the cohort. The median follow-up was 10.6 years. In all, 11 patients died from UTUC during the follow-up. Tumour mutational burden showed a statistically significant correlation with stage, grade, and stage + grade. Grade 1, Grade 2, and Grade 3 tumours had different mutational patterns. Patients who died from UTUC had pathogenic mutations in specific genes e.g. tumour protein p53 (TP53) and HRas proto-oncogene, GTPase (HRAS). Patients with Ta Grade 1 tumours with a known pathogenic fibroblast growth factor receptor 3 (FGFR3) mutation did not die from UTUC. CONCLUSION: The genetic analysis was highly concordant with histopathological features and added prognostic information in some cases. Thus, results from genomic profiling may contribute to the choice of treatment and follow-up regimens in the future.

A database on differentially expressed microRNAs during rodent bladder healing.

Urinary bladder wound healing relies on multiple biological events that are finely tuned in a spatial-temporal manner. MicroRNAs are small non-coding RNA molecules with regulatory functions. We hypothesized that microRNAs are important molecules in the coordination of normal urinary bladder wound healing. We aimed at identifying microRNAs expressed during bladder wound healing using Affymetrix global array for microRNA profiling of the rodent urinary bladder during healing of a surgically created wound. Results were validated in the rat bladders by real-time PCR (RT-PCR) using three of the differentially expressed (DE) microRNAs. The model was thereafter validated in human cells, by measuring the expression of eight of the DE microRNAs upon in vitro wound-healing assays in primary urothelial cells. Our results indicated that 508 (40%) of all rodent microRNAs were expressed in the urinary bladder during wound healing. Thirteen of these microRNAs (1%) were DE (false discovery rate (FDR) < 0.05, P < 0.05, |logfold|> 0.25) in wounded compared to non-wounded bladders. Bioinformatic analyses helped us to identify target molecules for the DE microRNAs, and biological pathways involved in tissue repair. All data are made available in an open-access database for other researchers to explore.

Chromoanagenesis Event Underlies a de novo Pericentric and Multiple Paracentric Inversions in a Single Chromosome Causing Coffin-Siris Syndrome.

Chromoanagenesis is a descriptive term that encompasses classes of catastrophic mutagenic processes that generate localized and complex chromosome rearrangements in both somatic and germline genomes. Herein, we describe a 5-year-old female presenting with a constellation of clinical features consistent with a clinical diagnosis of Coffin-Siris syndrome 1 (CSS1). Initial G-banded karyotyping detected a 90-Mb pericentric and a 47-Mb paracentric inversion on a single chromosome. Subsequent analysis of short-read whole-genome sequencing data and genomic optical mapping revealed additional inversions, all clustered on chromosome 6, one of them disrupting ARID1B for which haploinsufficiency leads to the CSS1 disease trait (MIM:135900). The aggregate structural variant data show that the resolved, the resolved derivative chromosome architecture presents four de novo inversions, one pericentric and three paracentric, involving six breakpoint junctions in what appears to be a shuffling of genomic material on this chromosome. Each junction was resolved to nucleotide-level resolution with mutational signatures suggestive of non-homologous end joining. The disruption of the gene ARID1B is shown to occur between the fourth and fifth exon of the canonical transcript with subsequent qPCR studies confirming a decrease in ARID1B expression in the patient versus healthy controls. Deciphering the underlying genomic architecture of chromosomal rearrangements and complex structural variants may require multiple technologies and can be critical to elucidating the molecular etiology of a patient's clinical phenotype or resolving unsolved Mendelian disease cases.

Single-cell multimodal analysis in a case with reduced penetrance of Progranulin-Frontotemporal Dementia.

We identified an autosomal dominant progranulin mutation carrier without symptoms of dementia in her lifetime (Reduced Penetrance Mutation Carrier, RedPenMC). This resistance to develop expected pathology presents a unique opportunity to interrogate neurodegenerative mechanisms. We performed multimodal single-nuclei analyses of post-mortem frontal cortex from RedPenMC, including transcriptomics and global levels of chromatin marks. RedPenMC had an increased ratio of GRN-expressing microglia, higher levels of activating histone mark H3k4me3 in microglia and lower levels of the repressive chromatin marks H3k9me1 and H3k9me3 in the frontal cortex than her affected mutation carrier son and evidence of higher protein levels of progranulin in both plasma and brain homogenates. Although the study is limited to one case, the results support that restoring brain progranulin levels may be sufficient to escape neurodegeneration and FTD. In addition to previously identified modifier genes, it is possible that epigenetic marks may contribute to the increased progranulin expression in cases of reduced penetrance. These findings may stimulate similar follow-up studies and new therapeutic approaches.

High diagnostic yield in skeletal ciliopathies using massively parallel genome sequencing, structural variant screening and RNA analyses.

Skeletal ciliopathies are a heterogenous group of disorders with overlapping clinical and radiographic features including bone dysplasia and internal abnormalities. To date, pathogenic variants in at least 30 genes, coding for different structural cilia proteins, are reported to cause skeletal ciliopathies. Here, we summarize genetic and phenotypic features of 34 affected individuals from 29 families with skeletal ciliopathies. Molecular diagnostic testing was performed using massively parallel sequencing (MPS) in combination with copy number variant (CNV) analyses and in silico filtering for variants in known skeletal ciliopathy genes. We identified biallelic disease-causing variants in seven genes: DYNC2H1, KIAA0753, WDR19, C2CD3, TTC21B, EVC, and EVC2. Four variants located in non-canonical splice sites of DYNC2H1, EVC, and KIAA0753 led to aberrant splicing that was shown by sequencing of cDNA. Furthermore, CNV analyses showed an intragenic deletion of DYNC2H1 in one individual and a 6.7 Mb de novo deletion on chromosome 1q24q25 in another. In five unsolved cases, MPS was performed in family setting. In one proband we identified a de novo variant in PRKACA and in another we found a homozygous intragenic deletion of IFT74, removing the first coding exon and leading to expression of a shorter message predicted to result in loss of 40 amino acids at the N-terminus. These findings establish IFT74 as a new skeletal ciliopathy gene. In conclusion, combined single nucleotide variant, CNV and cDNA analyses lead to a high yield of genetic diagnoses (90%) in a cohort of patients with skeletal ciliopathies.