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.

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.

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.

Comprehensive structural variation genome map of individuals carrying complex chromosomal rearrangements.

Complex chromosomal rearrangements (CCRs) are rearrangements involving more than two chromosomes or more than two breakpoints. Whole genome sequencing (WGS) allows for outstanding high resolution characterization on the nucleotide level in unique sequences of such rearrangements, but problems remain for mapping breakpoints in repetitive regions of the genome, which are known to be prone to rearrangements. Hence, multiple complementary WGS experiments are sometimes needed to solve the structures of CCRs. We have studied three individuals with CCRs: Case 1 and Case 2 presented with de novo karyotypically balanced, complex interchromosomal rearrangements (46,XX,t(2;8;15)(q35;q24.1;q22) and 46,XY,t(1;10;5)(q32;p12;q31)), and Case 3 presented with a de novo, extremely complex intrachromosomal rearrangement on chromosome 1. Molecular cytogenetic investigation revealed cryptic deletions in the breakpoints of chromosome 2 and 8 in Case 1, and on chromosome 10 in Case 2, explaining their clinical symptoms. In Case 3, 26 breakpoints were identified using WGS, disrupting five known disease genes. All rearrangements were subsequently analyzed using optical maps, linked-read WGS, and short-read WGS. In conclusion, we present a case series of three unique de novo CCRs where we by combining the results from the different technologies fully solved the structure of each rearrangement. The power in combining short-read WGS with long-molecule sequencing or optical mapping in these unique de novo CCRs in a clinical setting is demonstrated.

Chromatin interactions in differentiating keratinocytes reveal novel atopic dermatitis- and psoriasis-associated genes.

BACKGROUND: Hundreds of variants associated with atopic dermatitis (AD) and psoriasis, 2 common inflammatory skin disorders, have previously been discovered through genome-wide association studies (GWASs). The majority of these variants are in noncoding regions, and their target genes remain largely unclear. OBJECTIVE: We sought to understand the effects of these noncoding variants on the development of AD and psoriasis by linking them to the genes that they regulate. METHODS: We constructed genomic 3-dimensional maps of human keratinocytes during differentiation by using targeted chromosome conformation capture (Capture Hi-C) targeting more than 20,000 promoters and 214 GWAS variants and combined these data with transcriptome and epigenomic data sets. We validated our results with reporter assays, clustered regularly interspaced short palindromic repeats activation, and examination of patient gene expression from previous studies. RESULTS: We identified 118 target genes of 82 AD and psoriasis GWAS variants. Differential expression of 58 of the 118 target genes (49%) occurred in either AD or psoriatic lesions, many of which were not previously linked to any skin disease. We highlighted the genes AFG1L, CLINT1, ADO, LINC00302, and RP1-140J1.1 and provided further evidence for their potential roles in AD and psoriasis. CONCLUSIONS: Our work focused on skin barrier pathology through investigation of the interaction profile of GWAS variants during keratinocyte differentiation. We have provided a catalogue of candidate genes that could modulate the risk of AD and psoriasis. Given that only 35% of the target genes are the gene nearest to the known GWAS variants, we expect that our work will contribute to the discovery of novel pathways involved in AD and psoriasis.

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.

Discovery of Novel Sequences in 1,000 Swedish Genomes.

Novel sequences (NSs), not present in the human reference genome, are abundant and remain largely unexplored. Here, we utilize de novo assembly to study NS in 1,000 Swedish individuals first sequenced as part of the SweGen project revealing a total of 46 Mb in 61,044 distinct contigs of sequences not present in GRCh38. The contigs were aligned to recently published catalogs of Icelandic and Pan-African NSs, as well as the chimpanzee genome, revealing a great diversity of shared sequences. Analyzing the positioning of NS across the chimpanzee genome, we find that 2,807 NS align confidently within 143 chimpanzee orthologs of human genes. Aligning the whole genome sequencing data to the chimpanzee genome, we discover ancestral NS common throughout the Swedish population. The NSs were searched for repeats and repeat elements: revealing a majority of repetitive sequence (56%), and enrichment of simple repeats (28%) and satellites (15%). Lastly, we align the unmappable reads of a subset of the thousand genomes data to our collection of NS, as well as the previously published Pan-African NS: revealing that both the Swedish and Pan-African NS are widespread, and that the Swedish NSs are largely a subset of the Pan-African NS. Overall, these results highlight the importance of creating a more diverse reference genome and illustrate that significant amounts of the NS may be of ancestral origin.

Hybrid sequencing resolves two germline ultra-complex chromosomal rearrangements consisting of 137 breakpoint junctions in a single carrier.

Chromoanagenesis is a genomic event responsible for the formation of complex structural chromosomal rearrangements (CCRs). Germline chromoanagenesis is rare and the majority of reported cases are associated with an affected phenotype. Here, we report a healthy female carrying two de novo CCRs involving chromosomes 4, 19, 21 and X and chromosomes 7 and 11, respectively, with a total of 137 breakpoint junctions (BPJs). We characterized the CCRs using a hybrid-sequencing approach, combining short-read sequencing, nanopore sequencing, and optical mapping. The results were validated using multiple cytogenetic methods, including fluorescence in situ hybridization, spectral karyotyping, and Sanger sequencing. We identified 137 BPJs, which to our knowledge is the highest number of reported breakpoint junctions in germline chromoanagenesis. We also performed a statistical assessment of the positioning of the breakpoints, revealing a significant enrichment of BPJ-affecting genes (96 intragenic BPJs, 26 genes, p < 0.0001), indicating that the CCRs formed during active transcription of these genes. In addition, we find that the DNA fragments are unevenly and non-randomly distributed across the derivative chromosomes indicating a multistep process of scattering and re-joining of DNA fragments. In summary, we report a new maximum number of BPJs (137) in germline chromoanagenesis. We also show that a hybrid sequencing approach is necessary for the correct characterization of complex CCRs. Through in-depth statistical assessment, it was found that the CCRs most likely was formed through an event resembling chromoplexy-a catastrophic event caused by erroneous transcription factor binding.

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.

Novel pathogenic genomic variants leading to autosomal dominant and recessive Robinow syndrome.

Robinow syndrome (RS) is a genetically heterogeneous disorder characterized by skeletal dysplasia and a distinctive facial appearance. Previous studies have revealed locus heterogeneity with rare variants in DVL1, DVL3, FZD2, NXN, ROR2, and WNT5A underlying the etiology of RS. The aforementioned "Robinow-associated genes" and their gene products all play a role in the WNT/planar cell polarity signaling pathway. We performed gene-targeted Sanger sequencing, exome sequencing, genome sequencing, and array comparative genomic hybridization on four subjects with a clinical diagnosis of RS who had not had prior DNA testing. Individuals in our cohort were found to carry pathogenic or likely pathogenic variants in three RS related genes: DVL1, ROR2, and NXN. One subject was found to have a nonsense variant (c.817C > T [p.Gln273*]) in NXN in trans with an ~1 Mb telomeric deletion on chromosome 17p containing NXN, which supports our contention that biallelic NXN variant alleles are responsible for a novel autosomal recessive RS locus. These findings provide increased understanding of the role of WNT signaling in skeletal development and maintenance. These data further support the hypothesis that dysregulation of the noncanonical WNT pathway in humans gives rise to RS.

Replicative and non-replicative mechanisms in the formation of clustered CNVs are indicated by whole genome characterization.

Clustered copy number variants (CNVs) as detected by chromosomal microarray analysis (CMA) are often reported as germline chromothripsis. However, such cases might need further investigations by massive parallel whole genome sequencing (WGS) in order to accurately define the underlying complex rearrangement, predict the occurrence mechanisms and identify additional complexities. Here, we utilized WGS to delineate the rearrangement structure of 21 clustered CNV carriers first investigated by CMA and identified a total of 83 breakpoint junctions (BPJs). The rearrangements were further sub-classified depending on the patterns observed: I) Cases with only deletions (n = 8) often had additional structural rearrangements, such as insertions and inversions typical to chromothripsis; II) cases with only duplications (n = 7) or III) combinations of deletions and duplications (n = 6) demonstrated mostly interspersed duplications and BPJs enriched with microhomology. In two cases the rearrangement mutational signatures indicated both a breakage-fusion-bridge cycle process and haltered formation of a ring chromosome. Finally, we observed two cases with Alu- and LINE-mediated rearrangements as well as two unrelated individuals with seemingly identical clustered CNVs on 2p25.3, possibly a rare European founder rearrangement. In conclusion, through detailed characterization of the derivative chromosomes we show that multiple mechanisms are likely involved in the formation of clustered CNVs and add further evidence for chromoanagenesis mechanisms in both "simple" and highly complex chromosomal rearrangements. Finally, WGS characterization adds positional information, important for a correct clinical interpretation and deciphering mechanisms involved in the formation of these rearrangements.

pyCancerSig: subclassifying human cancer with comprehensive single nucleotide, structural and microsatellite mutational signature deconstruction from whole genome sequencing.

BACKGROUND: DNA damage accumulates over the course of cancer development. The often-substantial amount of somatic mutations in cancer poses a challenge to traditional methods to characterize tumors based on driver mutations. However, advances in machine learning technology can take advantage of this substantial amount of data. RESULTS: We developed a command line interface python package, pyCancerSig, to perform sample profiling by integrating single nucleotide variation (SNV), structural variation (SV) and microsatellite instability (MSI) profiles into a unified profile. It also provides a command to decipher underlying cancer processes, employing an unsupervised learning technique, Non-negative Matrix Factorization, and a command to visualize the results. The package accepts common standard file formats (vcf, bam). The program was evaluated using a cohort of breast- and colorectal cancer from The Cancer Genome Atlas project (TCGA). The result showed that by integrating multiple mutations modes, the tool can correctly identify cases with known clear mutational signatures and can strengthen signatures in cases with unclear signal from an SNV-only profile. The software package is available at https://github.com/jessada/pyCancerSig. CONCLUSIONS: pyCancerSig has demonstrated its capability in identifying known and unknown cancer processes, and at the same time, illuminates the association within and between the mutation modes.

Loqusdb: added value of an observations database of local genomic variation.

BACKGROUND: Exome and genome sequencing is becoming the method of choice for rare disease diagnostics. One of the key challenges remaining is distinguishing the disease causing variants from the benign background variation. After analysis and annotation of the sequencing data there are typically thousands of candidate variants requiring further investigation. One of the most effective and least biased ways to reduce this number is to assess the rarity of a variant in any population. Currently, there are a number of reliable sources of information for major population frequencies when considering single nucleotide variants (SNVs) and small insertion and deletions (INDELs), with gnomAD as the most prominent public resource available. However, local variation or frequencies in sub-populations may be underrepresented in these public resources. In contrast, for structural variation (SV), the background frequency in the general population is more or less unknown mostly due to challenges in calling SVs in a consistent way. Keeping track of local variation is one way to overcome these problems and significantly reduce the number of potential disease causing variants retained for manual inspection, both for SNVs and SVs. RESULTS: Here, we present loqusdb, a tool to solve the challenge of keeping track of any type of variant observations from genome sequencing data. Loqusdb was designed to handle a large flow of samples and unlike other solutions, samples can be added continuously to the database without rebuilding it, facilitating improvements and additions. We assessed the added value of a local observations database using 98 samples annotated with information from a background of 888 unrelated individuals. CONCLUSIONS: We show both how powerful SV analysis can be when filtering for population frequencies and how the number of apparently rare SNVs/INDELs can be reduced by adding local population information even after annotating the data with other large frequency databases, such as gnomAD. In conclusion, we show that a local frequency database is an attractive, and a necessary addition to the publicly available databases that facilitate the analysis of exome and genome data in a clinical setting.

Cytogenetically visible inversions are formed by multiple molecular mechanisms.

Cytogenetically detected inversions are generally assumed to be copy number and phenotypically neutral events. While nonallelic homologous recombination is thought to play a major role, recent data suggest the involvement of other molecular mechanisms in inversion formation. Using a combination of short-read whole-genome sequencing (WGS), 10X Genomics Chromium WGS, droplet digital polymerase chain reaction and array comparative genomic hybridization we investigated the genomic structure of 18 large unique cytogenetically detected chromosomal inversions and achieved nucleotide resolution of at least one chromosomal inversion junction for 13/18 (72%). Surprisingly, we observed that seemingly copy number neutral inversions can be accompanied by a copy-number gain of up to 350 kb and local genomic complexities (3/18, 17%). In the resolved inversions, the mutational signatures are consistent with nonhomologous end-joining (8/13, 62%) or microhomology-mediated break-induced replication (5/13, 38%). Our study indicates that short-read 30x coverage WGS can detect a substantial fraction of chromosomal inversions. Moreover, replication-based mechanisms are responsible for approximately 38% of those events leading to a significant proportion of inversions that are actually accompanied by additional copy-number variation potentially contributing to the overall phenotypic presentation of those patients.

Rare variants in dynein heavy chain genes in two individuals with situs inversus and developmental dyslexia: a case report.

BACKGROUND: Developmental dyslexia (DD) is a neurodevelopmental learning disorder with high heritability. A number of candidate susceptibility genes have been identified, some of which are linked to the function of the cilium, an organelle regulating left-right asymmetry development in the embryo. Furthermore, it has been suggested that disrupted left-right asymmetry of the brain may play a role in neurodevelopmental disorders such as DD. However, it is unknown whether there is a common genetic cause to DD and laterality defects or ciliopathies. CASE PRESENTATION: Here, we studied two individuals with co-occurring situs inversus (SI) and DD using whole genome sequencing to identify genetic variants of importance for DD and SI. Individual 1 had primary ciliary dyskinesia (PCD), a rare, autosomal recessive disorder with oto-sino-pulmonary phenotype and SI. We identified two rare nonsynonymous variants in the dynein axonemal heavy chain 5 gene (DNAH5): a previously reported variant c.7502G > C; p.(R2501P), and a novel variant c.12043 T > G; p.(Y4015D). Both variants are predicted to be damaging. Ultrastructural analysis of the cilia revealed a lack of outer dynein arms and normal inner dynein arms. MRI of the brain revealed no significant abnormalities. Individual 2 had non-syndromic SI and DD. In individual 2, one rare variant (c.9110A > G;p.(H3037R)) in the dynein axonemal heavy chain 11 gene (DNAH11), coding for another component of the outer dynein arm, was identified. CONCLUSIONS: We identified the likely genetic cause of SI and PCD in one individual, and a possibly significant heterozygosity in the other, both involving dynein genes. Given the present evidence, it is unclear if the identified variants also predispose to DD and further studies into the association between laterality, ciliopathies and DD are needed.

Whole-genome sequencing reveals complex chromosome rearrangement disrupting NIPBL in infant with Cornelia de Lange syndrome.

Clinical laboratory diagnostic evaluation of the genomes of children with suspected genetic disorders, including chromosomal microarray and exome sequencing, cannot detect copy number neutral genomic rearrangements such as inversions, balanced translocations, and complex chromosomal rearrangements (CCRs). We describe an infant with a clinical diagnosis of Cornelia de Lange syndrome (CdLS) in whom chromosome analysis revealed a de novo complex balanced translocation, 46,XY,t(5;7;6)(q11.2;q32;q13)dn. Subsequent molecular characterization by whole-genome sequencing (WGS) identified 23 breakpoints, delineating segments derived from four chromosomes (5;6;7;21) in ancestral or inverted orientation. One of the breakpoints disrupted a known CdLS gene, NIPBL. Further investigation revealed paternal origin of the CCR allele, clustering of the breakpoint junctions, and molecular repair signatures suggestive of a single catastrophic event. Notably, very short DNA segments (25 and 41 bp) were included in the reassembled chromosomes, lending additional support that the DNA repair machinery can detect and repair such segments. Interestingly, there was an independent paternally derived miniscule complex rearrangement, possibly predisposing to subsequent genomic instability. In conclusion, we report a CCR causing a monogenic Mendelian disorder, urging WGS analysis of similar unsolved cases with suspected Mendelian disorders. Breakpoint analysis allowed for identification of the underlying molecular diagnosis and implicated chromoanagenesis in CCR formation.

Targeted copy number screening highlights an intragenic deletion of WDR63 as the likely cause of human occipital encephalocele and abnormal CNS development in zebrafish.

Congenital malformations affecting the neural tube can present as isolated malformations or occur in association with other developmental abnormalities and syndromes. Using high-resolution copy number screening in 66 fetuses with neural tube defects, we identified six fetuses with likely pathogenic mutations, three aneuploidies (one trisomy 13 and two trisomy 18) and three deletions previously reported in NTDs (one 22q11.2 deletion and two 1p36 deletions) corresponding to 9% of the cohort. In addition, we identified five rare deletions and two duplications of uncertain significance including a rare intragenic heterozygous in-frame WDR63 deletion in a fetus with occipital encephalocele. Whole genome sequencing verified the deletion and excluded known pathogenic variants. The deletion spans exons 14-17 resulting in the expression of a protein missing the third and fourth WD-repeat domains. These findings were supported by CRISPR/Cas9-mediated somatic deletions in zebrafish. Injection of two different sgRNA-pairs targeting relevant intronic regions resulted in a deletion mimicking the human deletion and a concomitant increase of abnormal embryos with body and brain malformations (41%, n = 161 and 62%, n = 224, respectively), including a sac-like brain protrusion (7% and 9%, P < 0.01). Similar results were seen with overexpression of RNA encoding the deleted variant in zebrafish (total abnormal; 46%, n = 255, P < 0.001) compared with the overexpression of an equivalent amount of wild-type RNA (total abnormal; 3%, n = 177). We predict the in-frame WDR63 deletion to result in a dominant negative or gain-of-function form of WDR63. These are the first findings supporting a role for WDR63 in encephalocele formation.

Alu-Alu mediated intragenic duplications in IFT81 and MATN3 are associated with skeletal dysplasias.

Skeletal dysplasias are a diverse group of rare Mendelian disorders with clinical and genetic heterogeneity. Here, we used targeted copy number variant (CNV) screening and identified intragenic exonic duplications, formed through Alu-Alu fusion events, in two individuals with skeletal dysplasia and negative exome sequencing results. First, we detected a homozygous tandem duplication of exon 9 and 10 in IFT81 in a boy with Jeune syndrome, or short-rib thoracic dysplasia (SRTD) (MIM# 208500). Western blot analysis did not detect any wild-type IFT81 protein in fibroblasts from the patient with the IFT81 duplication, but only a shorter isoform of IFT81 that was also present in the normal control samples. Complementary zebrafish studies suggested that loss of full-length IFT81 protein but expression of a shorter form of IFT81 protein affects the phenotype while being compatible with life. Second, a de novo tandem duplication of exons 2 to 5 in MATN3 was identified in a girl with multiple epiphyseal dysplasia (MED) type 5 (MIM# 607078). Our data highlights the importance of detection and careful characterization of intragenic duplication CNVs, presenting them as a novel and very rare genetic mechanism in IFT81-related Jeune syndrome and MATN3-related MED.

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.

Uncommon Variants in FLG2 and NOD2 Are Associated with Atopic Dermatitis in the Ethiopian Population.

Loss-of-function variants in the FLG gene have been identified as the strongest cause of susceptibility to atopic dermatitis (AD) in Europeans and Asians. However, very little is known about the genetic etiology behind AD in African populations, where the prevalence of AD is notably high. We sought to investigate the genetic origins of AD by performing whole-genome sequencing in an Ethiopian family with 12 individuals and several affected in different generations. We identified 2 variants within FLG2 (p.D13Y) and NOD2 (p.A918S) genes cosegregating with AD in the affected individuals. Further genotyping analyses in both Ethiopian and Swedish AD cases and controls revealed a significant association with the FLG2 variant (p.D13Y, P < .0013) only in the Ethiopian cohort. However, the NOD2 variant (p.A918S) did not show any association in our Ethiopian cohort. Instead, 2 previously recognized NOD2 variants (p.A849V, P < .0085 and p.G908R, P < .0036) were significantly associated with AD in our Ethiopian cohort. Our study indicates that the FLG2 and NOD2 genes might be important in the etiology of AD in Ethiopians. Additional genetic and functional studies are needed to confirm the role of these genes and the associated variants into the development of AD.