Genomic disorders: structural features of the genome can lead to DNA rearrangements and human disease traits

Genome integrity as a potential index of longevity in Ashkenazi Centenarian's families

The aging process, or senescence, is characterized by age-specific decline in physical and physiological function, and increased frailty and genomic changes, including mutation accumulation. However, the mechanisms through which changes in genomic architecture influence human longevity have remained obscure. Copy number variants (CNVs), an abundant class of genomic variants, offer unique opportunities for understanding age-related genomic changes. Here we report the spectrum of CNVs in a cohort of 670 Ashkenazi Jewish centenarians, their progeny, and unrelated controls. The average ages of these groups were 97.4 ± 2.8, 69.2 ± 9.2, and 66.5 ± 7.0 respectively. For the first time, we compared different size classes of CNVs, from 1 kB to 100 MB in size. Using a high-resolution custom Affymetrix array, targeting 44,639 genomic regions, we identified a total of 12,166, 22,188, and 10,285 CNVs in centenarians, their progeny, and control groups, respectively. Interestingly, the offspring group showed the highest number of unique CNVs, followed by control and centenarians. While both gains and losses were found in all three groups, centenarians showed a significantly higher average number of both total gains and losses relative to their controls (p < 0.0327, 0.0182, respectively). Moreover, centenarians showed a lower total length of genomic material lost, suggesting that they may maintain superior genomic integrity over time. We also observe a significance fold increase of CNVs among the offspring, implying greater genomic integrity and a putative mechanism for longevity preservation. Genomic regions that experienced loss or gains appear to be distributed across many sites in the genome and contain genes involved in DNA transcription, cellular transport, developmental pathways, and metabolic functions. Our findings suggest that the exceptional longevity observed in centenarians may be attributed to the prolonged maintenance of functionally important genes. These genes are intrinsic to specific genomic regions as well as to the overall integrity of the genomic architecture. Additionally, a strong association between longer CNVs and differential gene expression observed in this study supports the notion that genomic integrity could positively influence longevity.

Fc gamma receptors: Their evolution, genomic architecture, genetic variation, and impact on human disease

Fc gamma receptors (FcγRs) are a family of receptors that bind IgG antibodies and interface at the junction of humoral and innate immunity. Precise regulation of receptor expression provides the necessary balance to achieve healthy immune homeostasis by establishing an appropriate immune threshold to limit autoimmunity but respond effectively to infection. The underlying genetics of the FCGR gene family are central to achieving this immune threshold by regulating affinity for IgG, signaling efficacy, and receptor expression. The FCGR gene locus was duplicated during evolution, retaining very high homology and resulting in a genomic region that is technically difficult to study. Here, we review the recent evolution of the gene family in mammals, its complexity and variation through copy number variation and single-nucleotide polymorphism, and impact of these on disease incidence, resolution, and therapeutic antibody efficacy. We also discuss the progress and limitations of current approaches to study the region and emphasize how new genomics technologies will likely resolve much of the current confusion in the field. This will lead to definitive conclusions on the impact of genetic variation within the FCGR gene locus on immune function and disease.

Development of Lasing Silica Microsphere for High-Speed DNA Molecular Detection

Laser and molecular detection techniques that have been used to overcome the limitations of fluorescent DNA labeling have presented new challenges. To address some of these challenges, we developed a DNA laser that uses a solid-state silica microsphere as a ring resonator and a site for DNA-binding reactions, as well as a platform to detect and sequence target DNA molecules. We detected target DNA using laser emission from a DNA-labeling dye and a developed solid-state silica microsphere ring resonator. The microsphere was sensitive; a single base mismatch in the DNA resulted in the absence of an optical signal. As each individual microsphere can be utilized as a parallel DNA analysis chamber, this optical digital detection scheme allows for high-throughput and rapid analysis. More importantly, the solid-state DNA laser is free from deformation, which guarantees stable lasing characteristics, and can be manipulated freely outside the solution. Thus, this promising advanced DNA laser scheme can be implemented on platforms other than optofluidic chips.

GGTyper: genotyping complex structural variants using short-read sequencing data

MOTIVATION

Complex structural variants (SVs) are genomic rearrangements that involve multiple segments of DNA. They contribute to human diversity and have been shown to cause Mendelian disease. Nevertheless, our abilities to analyse complex SVs are very limited. As opposed to deletions and other canonical types of SVs, there are no established tools that have explicitly been designed for analysing complex SVs.

RESULTS

Here, we describe a new computational approach that we specifically designed for genotyping complex SVs in short-read sequenced genomes. Given a variant description, our approach computes genotype-specific probability distributions for observing aligned read pairs with a wide range of properties. Subsequently, these distributions can be used to efficiently determine the most likely genotype for any set of aligned read pairs observed in a sequenced genome. In addition, we use these distributions to compute a genotyping difficulty for a given variant, which predicts the amount of data needed to achieve a reliable call. Careful evaluation confirms that our approach outperforms other genotypers by making reliable genotype predictions across both simulated and real data. On up to 7829 human genomes, we achieve high concordance with population-genetic assumptions and expected inheritance patterns. On simulated data, we show that precision correlates well with our prediction of genotyping difficulty. This together with low memory and time requirements makes our approach well-suited for application in biomedical studies involving small to very large numbers of short-read sequenced genomes.

AVAILABILITY AND IMPLEMENTATION

Source code is available at https://github.com/kehrlab/Complex-SV-Genotyping.

Analysis of parental origin of de novo pathogenic CNVs in patients with intellectual disability

Chromosomal Microarray Analysis (CMA) has increased the comprehension of the mechanisms of copy number variation (CNV) formation, classification of these rearrangements, type of recurrence, and its origin, and has also been a powerful approach to identifying CNVs in individuals with intellectual disability. The aim of this study was to establish the parental origin of de novo pathogenic CNV in a cohort of patients with intellectual disability from the public health system of Goiás-Brazil. CMA was done in 76 trios and we identified 15 de novo pathogenic CNVs in 12 patients with intellectual disability. In a total of 15 de novo pathogenic CNV, 60% were derived from the maternal germline and 40% from the paternal germline. CNV flanked by low copy repeats (LCR) were identified in 46.7% and most of them were of maternal origin. No significant association was observed between paternal age and the mutation rate of de novo CNVs. The presence of high-identity LCRs increases the occurrence of CNV formation mediated by non-allelic homologous recombination and the majority of paternal CNVs are non-recurrent. The mechanism of formation of these CNV may have been by microhomology-mediated break-induced replication or non-homologous end joining.

Modelling DNA damage-repair and beyond

The paper presents a review of mechanistic modelling studies of DNA damage and DNA repair, and consequences to follow in mammalian cell nucleus. We hypothesize DNA deletions are consequences of repair of double strand breaks leading to the modifications of genome that play crucial role in long term development of genetic inheritance and diseases. The aim of the paper is to review formation mechanisms underlying naturally occurring DNA deletions in the human genome and their potential relevance for bridging the gap between induced DNA double strand breaks and deletions in damaged human genome from endogenous and exogenous events. The model of the cell nucleus presented enables simulation of DNA damage at molecular level identifying the spectrum of damage induced in all chromosomal territories and loops. Our mechanistic modelling of DNA repair for double stand breaks (DSB), single strand breaks (SSB) and base damage (BD), shows the complexity of DNA damage is responsible for the longer repair times and the reason for the biphasic feature of mammalian cells repair curves. In the absence of experimentally determined data, the mechanistic model of repair predicts the in vivo rate constants for the proteins involved in the repair of DSB, SSB, and of BD.

NanoRanger enables rapid single-base-pair resolution of genomic disorders

BACKGROUND

Delineating base-resolution breakpoints of complex rearrangements is crucial for an accurate clinical understanding of pathogenic variants and for carrier screening within family networks or the broader population. However, despite advances in genetic testing using short-read sequencing (SRS), this task remains costly and challenging.

METHODS

This study addresses the challenges of resolving missing disease-causing breakpoints in complex genomic disorders with suspected homozygous rearrangements by employing multiple long-read sequencing (LRS) strategies, including a novel and efficient strategy named nanopore-based rapid acquisition of neighboring genomic regions (NanoRanger). NanoRanger does not require large amounts of ultrahigh-molecular-weight DNA and stands out for its ease of use and rapid acquisition of large genomic regions of interest with deep coverage.

FINDINGS

We describe a cohort of 16 familial cases, each harboring homozygous rearrangements that defied breakpoint determination by SRS and optical genome mapping (OGM). NanoRanger identified the breakpoints with single-base-pair resolution, enabling accurate determination of the carrier status of unaffected family members as well as the founder nature of these genomic lesions and their frequency in the local population. The resolved breakpoints revealed that repetitive DNA, gene regulatory elements, and transcription activity contribute to genome instability in these novel recessive rearrangements.

CONCLUSIONS

Our data suggest that NanoRanger greatly improves the success rate of resolving base-resolution breakpoints of complex genomic disorders and expands access to LRS for the benefit of patients with Mendelian disorders.

FUNDING

M.L. is supported by KAUST Baseline Award no. BAS/1/1080-01-01 and KAUST Research Translation Fund Award no. REI/1/4742-01.

Low-pass whole genome sequencing as a cost-effective alternative to chromosomal microarray analysis for low- and middle-income countries

Low-pass whole genome sequencing (LP-WGS) has been applied as alternative method to detect copy number variants (CNVs) in the clinical setting. Compared with chromosomal microarray analysis (CMA), the sequencing-based approach provides a similar resolution of CNV detection at a lower cost. In this study, we assessed the efficiency and reliability of LP-WGS as a more affordable alternative to CMA. A total of 1363 patients with unexplained neurodevelopmental delay/intellectual disability, autism spectrum disorders, and/or multiple congenital anomalies were enrolled. Those patients were referred from 15 nonprofit organizations and university centers located in different states in Brazil. The analysis of LP-WGS at 1x coverage (>50kb) revealed a positive testing result in 22% of the cases (304/1363), in which 219 and 85 correspond to pathogenic/likely pathogenic (P/LP) CNVs and variants of uncertain significance (VUS), respectively. The 16% (219/1363) diagnostic yield observed in our cohort is comparable to the 15%-20% reported for CMA in the literature. The use of commercial software, as demonstrated in this study, simplifies the implementation of the test in clinical settings. Particularly for countries like Brazil, where the cost of CMA presents a substantial barrier to most of the population, LP-WGS emerges as a cost-effective alternative for investigating copy number changes in cytogenetics.

Olfactory receptor genes and chromosome 11 structural aberrations: Players or spectators?

The largest multi-gene family in metazoans is the family of olfactory receptor (OR) genes. Human ORs are organized in clusters over most chromosomes and seem to include >0.1% the human genome. Because 369 out of 856 OR genes are mapped on chromosome 11 (HSA11), we sought to determine whether they mediate structural rearrangements involving this chromosome. To this aim, we analyzed 220 specimens collected during diagnostic procedures involving structural rearrangements of chromosome 11. A total of 222 chromosomal abnormalities were included, consisting of inversions, deletions, translocations, duplications, and one insertion, detected by conventional chromosome analysis and/or fluorescence in situ hybridization (FISH) and array comparative genomic hybridization (array-CGH). We verified by bioinformatics and statistical approaches the occurrence of breakpoints in cytobands with or without OR genes. We found that OR genes are not involved in chromosome 11 reciprocal translocations, suggesting that different DNA motifs and mechanisms based on homology or non-homology recombination can cause chromosome 11 structural alterations. We also considered the proximity between the chromosomal territories of chromosome 11 and its partner chromosomes involved in the translocations by using the deposited Hi-C data concerning the possible occurrence of chromosome interactions. Interestingly, most of the breakpoints are located in regions highly involved in chromosome interactions. Further studies should be carried out to confirm the potential role of chromosome territories' proximity in promoting genome structural variation, so fundamental in our understanding of the molecular basis of medical genetics and evolutionary genetics.

Chromosome Transplantation: Opportunities and Limitations

There are thousands of rare genetic diseases that could be treated with classical gene therapy strategies such as the addition of the defective gene via viral or non-viral delivery or by direct gene editing. However, several genetic defects are too complex for these approaches. These "genomic mutations" include aneuploidies, intra and inter chromosomal rearrangements, large deletions, or inversion and copy number variations. Chromosome transplantation (CT) refers to the precise substitution of an endogenous chromosome with an exogenous one. By the addition of an exogenous chromosome and the concomitant elimination of the endogenous one, every genetic defect, irrespective of its nature, could be resolved. In the current review, we analyze the state of the art of this technique and discuss its possible application to human pathology. CT might not be limited to the treatment of human diseases. By working on sex chromosomes, we showed that female cells can be obtained from male cells, since chromosome-transplanted cells can lose either sex chromosome, giving rise to 46,XY or 46,XX diploid cells, a modification that could be exploited to obtain female gametes from male cells. Moreover, CT could be used in veterinary biology, since entire chromosomes containing an advantageous locus could be transferred to animals of zootechnical interest without altering their specific genetic background and the need for long and complex interbreeding. CT could also be useful to rescue extinct species if only male cells were available. Finally, the generation of "synthetic" cells could be achieved by repeated CT into a recipient cell. CT is an additional tool for genetic modification of mammalian cells.

Extended application of PGT-M strategies for small pathogenic CNVs

PURPOSE

The preimplantation genetic testing for aneuploidy (PGT-A) platform is not currently available for small copy-number variants (CNVs), especially those < 1 Mb. Through strategies used in PGT for monogenic disease (PGT-M), this study intended to perform PGT for families with small pathogenic CNVs.

METHODS

Couples who carried small pathogenic CNVs and underwent PGT at the Reproductive and Genetic Hospital of CITIC-Xiangya (Hunan, China) between November 2019 and April 2023 were included in this study. Haplotype analysis was performed through two platforms (targeted sequencing and whole-genome arrays) to identify the unaffected embryos, which were subjected to transplantation. Prenatal diagnosis using amniotic fluid was performed during 18-20 weeks of pregnancy.

RESULTS

PGT was successfully performed for 20 small CNVs (15 microdeletions and 5 microduplications) in 20 families. These CNVs distributed on chromosomes 1, 2, 6, 7, 13, 15, 16, and X with sizes ranging from 57 to 2120 kb. Three haplotyping-based PGT-M strategies were applied. A total of 89 embryos were identified in 25 PGT cycles for the 20 families. The diagnostic yield was 98.9% (88/89). Nineteen transfers were performed for 17 women, resulting in a 78.9% (15/19) clinical pregnancy rate after each transplantation. Of the nine women who had healthy babies, eight accepted prenatal diagnosis and the results showed no related pathogenic CNVs.

CONCLUSION

Our results show that the extended haplotyping-based PGT-M strategy application for small pathogenic CNVs compensated for the insufficient resolution of PGT-A. These three PGT-M strategies could be applied to couples with small pathogenic CNVs.

Low-pass whole genome sequencing is a reliable and cost-effective approach for copy number variant analysis in the clinical setting

INTRODUCTION

Next generation sequencing technology has greatly reduced the cost and time required for sequencing a genome. An approach that is rapidly being adopted as an alternative method for CNV analysis is the low-pass whole genome sequencing (LP-WGS). Here, we evaluated the performance of LP-WGS to detect copy number variants (CNVs) in clinical cytogenetics.

MATERIALS AND METHODS

DNA samples with known CNVs detected by chromosomal microarray analyses (CMA) were selected for comparison and used as positive controls; our panel included 44 DNA samples (12 prenatal and 32 postnatal), comprising a total of 55 chromosome imbalances. The selected cases were chosen to provide a wide range of clinically relevant CNVs, the vast majority being associated with intellectual disability or recognizable syndromes. The chromosome imbalances ranged in size from 75 kb to 90.3 Mb, including aneuploidies and two cases of mosaicism.

RESULTS

All CNVs were successfully detected by LP-WGS, showing a high level of consistency and robust performance of the sequencing method. Notably, the size of chromosome imbalances detected by CMA and LP-WGS were compatible between the two different platforms, which indicates that the resolution and sensitivity of the LP-WGS approach are at least similar to those provided by CMA.

DISCUSSION

Our data show the potential use of LP-WGS to detect CNVs in clinical diagnosis and confirm the method as an alternative for chromosome imbalances detection. The diagnostic effectiveness and feasibility of LP-WGS, in this technical validation study, were evidenced by a clinically representative dataset of CNVs that allowed a systematic assessment of the detection power and the accuracy of the sequencing approach. Further, since the software used in this study is commercially available, the method can easily be tested and implemented in a routine diagnostic setting.

Scrambling the genome in cancer: causes and consequences of complex chromosome rearrangements

Complex chromosome rearrangements, known as chromoanagenesis, are widespread in cancer. Based on large-scale DNA sequencing of human tumours, the most frequent type of complex chromosome rearrangement is chromothripsis, a massive, localized and clustered rearrangement of one (or a few) chromosomes seemingly acquired in a single event. Chromothripsis can be initiated by mitotic errors that produce a micronucleus encapsulating a single chromosome or chromosomal fragment. Rupture of the unstable micronuclear envelope exposes its chromatin to cytosolic nucleases and induces chromothriptic shattering. Found in up to half of tumours included in pan-cancer genomic analyses, chromothriptic rearrangements can contribute to tumorigenesis through inactivation of tumour suppressor genes, activation of proto-oncogenes, or gene amplification through the production of self-propagating extrachromosomal circular DNAs encoding oncogenes or genes conferring anticancer drug resistance. Here, we discuss what has been learned about the mechanisms that enable these complex genomic rearrangements and their consequences in cancer.

Evaluation of genetic risk of apparently balanced chromosomal rearrangement carriers by breakpoint characterization

PURPOSE

To report genetic characteristics and associated risk of chromosomal breaks due to chromosomal rearrangements in large samples.

METHODS

MicroSeq, a technique that combines chromosome microdissection and next-generation sequencing, was used to identify chromosomal breakpoints. Long-range PCR and Sanger sequencing were used to precisely characterize 100 breakpoints in 50 ABCR carriers.

RESULTS

In addition to the recurrent regions of balanced rearrangement breaks in 8q24.13, 11q11.23, and 22q11.21 that had been documented, we have discovered a 10-Mb region of 12q24.13-q24.3 that could potentially be a sparse region of balanced rearrangement breaks. We found that 898 breakpoints caused gene disruption and a total of 188 breakpoints interrupted genes recorded in OMIM. The percentage of breakpoints that disrupted autosomal dominant genes recorded in OMIM was 25.53% (48/188). Fifty-four of the precisely characterized breakpoints had 1-8-bp microhomologous sequences.

CONCLUSION

Our findings provide a reference for the evaluation of the pathogenicity of mutations in related genes that cause protein truncation in clinical practice. According to the characteristics of breakpoints, non-homologous end joining and microhomology-mediated break-induced replication may be the main mechanism for ABCRs formation.

Copy Number Variations in Pancreatic Cancer: From Biological Significance to Clinical Utility

Pancreatic ductal adenocarcinoma (PDAC) is the most common type of pancreatic cancer, characterized by high tumor heterogeneity and a poor prognosis. Inter- and intra-tumoral heterogeneity in PDAC is a major obstacle to effective PDAC treatment; therefore, it is highly desirable to explore the tumor heterogeneity and underlying mechanisms for the improvement of PDAC prognosis. Gene copy number variations (CNVs) are increasingly recognized as a common and heritable source of inter-individual variation in genomic sequence. In this review, we outline the origin, main characteristics, and pathological aspects of CNVs. We then describe the occurrence of CNVs in PDAC, including those that have been clearly shown to have a pathogenic role, and further highlight some key examples of their involvement in tumor development and progression. The ability to efficiently identify and analyze CNVs in tumor samples is important to support translational research and foster precision oncology, as copy number variants can be utilized to guide clinical decisions. We provide insights into understanding the CNV landscapes and the role of both somatic and germline CNVs in PDAC, which could lead to significant advances in diagnosis, prognosis, and treatment. Although there has been significant progress in this field, understanding the full contribution of CNVs to the genetic basis of PDAC will require further research, with more accurate CNV assays such as single-cell techniques and larger cohorts than have been performed to date.

Delineation of an inverted tandem Xq23-26.3 duplication in a female featuring extremely short stature and mild mental deficiency

BACKGROUND

Partial duplications involving the long arm of the X chromosome are associated with mental retardation, short stature, microcephaly, and a wide range of physical findings. Female carriers usually have no clinical phenotype. Occasionally, they may also have heterogeneous features due to non-random inactivation of the X chromosome.

METHODS

The peripheral blood sample was collected from the patient and subjected to a few genetic testing, including chromosomal karyotyping, Chromosomal microarray analysis (CMA), Optical genome mapping, short tandem repeat (STR) analysis for Determination of parental origin, and X chromosome inactivation (XCI) analysis.

RESULTS

We have identified a de novo Xq23-Xq26.3 duplication in an adult female featuring extremely short stature and mild mental deficiency. Chromosome analysis detected a duplication on Xq23-q26.3 with a size of approximately 20 Mb. The duplication region has encompassed a number of genes, among which ARHGEF6, PHF6, HPRT1 and SLC9A6 are associated with X-linked mental retardation. Further analysis suggested that the duplication has derived from her father, was of the inversion duplication type and involved various degrees of skewed X chromosome inactivation.

CONCLUSION

Correlation with her phenotypes might indicate new mechanisms by which the X chromosome may lead to short stature and mental retardation. Our findings thereby may shed more light on the phenotypic implication of functional disomy of X-chromosome genes.

Tracing cancer evolution and heterogeneity using Hi-C

Chromosomal rearrangements can initiate and drive cancer progression, yet it has been challenging to evaluate their impact, especially in genetically heterogeneous solid cancers. To address this problem we developed HiDENSEC, a new computational framework for analyzing chromatin conformation capture in heterogeneous samples that can infer somatic copy number alterations, characterize large-scale chromosomal rearrangements, and estimate cancer cell fractions. After validating HiDENSEC with in silico and in vitro controls, we used it to characterize chromosome-scale evolution during melanoma progression in formalin-fixed tumor samples from three patients. The resulting comprehensive annotation of the genomic events includes copy number neutral translocations that disrupt tumor suppressor genes such as NF1, whole chromosome arm exchanges that result in loss of CDKN2A, and whole-arm copy-number neutral loss of homozygosity involving PTEN. These findings show that large-scale chromosomal rearrangements occur throughout cancer evolution and that characterizing these events yields insights into drivers of melanoma progression.

CRISPR/Cas9-induced DNA breaks trigger crossover, chromosomal loss, and chromothripsis-like rearrangements

DNA double-stranded breaks (DSBs) generated by the Cas9 nuclease are commonly repaired via nonhomologous end-joining (NHEJ) or homologous recombination (HR). However, little is known about unrepaired DSBs and the type of damage they trigger in plants. We designed an assay that detects loss of heterozygosity (LOH) in somatic cells, enabling the study of a broad range of DSB-induced genomic events. The system relies on a mapped phenotypic marker which produces a light purple color (betalain pigment) in all plant tissues. Plants with sectors lacking the Betalain marker upon DSB induction between the marker and the centromere were tested for LOH events. Using this assay, we detected a tomato (Solanum lycopersicum) flower with a twin yellow and dark purple sector, corresponding to a germinally transmitted somatic crossover event. We also identified instances of small deletions of genomic regions spanning the T-DNA and whole chromosome loss. In addition, we show that major chromosomal rearrangements including loss of large fragments, inversions, and translocations were clearly associated with the CRISPR-induced DSB. Detailed characterization of complex rearrangements by whole-genome sequencing and molecular and cytological analyses supports a model in which a breakage-fusion-bridge cycle followed by chromothripsis-like rearrangements had been induced. Our LOH assay provides a tool for precise breeding via targeted crossover detection. It also uncovers CRISPR-mediated chromothripsis-like events in plants.

Chromosomal deletions on 16p11.2 encompassing SH2B1 are associated with accelerated metabolic disease

New approaches are needed to treat people whose obesity and type 2 diabetes (T2D) are driven by specific mechanisms. We investigate a deletion on chromosome 16p11.2 (breakpoint 2-3 [BP2-3]) encompassing SH2B1, a mediator of leptin and insulin signaling. Phenome-wide association scans in the UK (N = 502,399) and Estonian (N = 208,360) biobanks show that deletion carriers have increased body mass index (BMI; p = 1.3 × 10-10) and increased rates of T2D. Compared with BMI-matched controls, deletion carriers have an earlier onset of T2D, with poorer glycemic control despite higher medication usage. Cystatin C, a biomarker of kidney function, is significantly elevated in deletion carriers, suggesting increased risk of renal impairment. In a Mendelian randomization study, decreased SH2B1 expression increases T2D risk (p = 8.1 × 10-6). We conclude that people with 16p11.2 BP2-3 deletions have early, complex obesity and T2D and may benefit from therapies that enhance leptin and insulin signaling.

Germline CNV Detection through Whole-Exome Sequencing (WES) Data Analysis Enhances Resolution of Rare Genetic Diseases

Whole-Exome Sequencing (WES) has proven valuable in the characterization of underlying genetic defects in most rare diseases (RDs). Copy Number Variants (CNVs) were initially thought to escape detection. Recent technological advances enabled CNV calling from WES data with the use of accurate and highly sensitive bioinformatic tools. Amongst 920 patients referred for WES, 454 unresolved cases were further analysed using the ExomeDepth algorithm. CNVs were called, evaluated and categorized according to ACMG/ClinGen recommendations. Causative CNVs were identified in 40 patients, increasing the diagnostic yield of WES from 50.7% (466/920) to 55% (506/920). Twenty-two CNVs were available for validation and were all confirmed; of these, five were novel. Implementation of the ExomeDepth tool promoted effective identification of phenotype-relevant and/or novel CNVs. Among the advantages of calling CNVs from WES data, characterization of complex genotypes comprising both CNVs and SNVs minimizes cost and time to final diagnosis, while allowing differentiation between true or false homozygosity, as well as compound heterozygosity of variants in AR genes. The use of a specific algorithm for calling CNVs from WES data enables ancillary detection of different types of causative genetic variants, making WES a critical first-tier diagnostic test for patients with RDs.

Transmembrane nuclease NUMEN/ENDOD1 regulates DNA repair pathway choice at the nuclear periphery

Proper repair of DNA damage lesions is essential to maintaining genome integrity and preventing the development of human diseases, including cancer. Increasing evidence suggests the importance of the nuclear envelope in the spatial regulation of DNA repair, although the mechanisms of such regulatory processes remain poorly defined. Through a genome-wide synthetic viability screen for PARP-inhibitor resistance using an inducible CRISPR-Cas9 platform and BRCA1-deficient breast cancer cells, we identified a transmembrane nuclease (renamed NUMEN) that could facilitate compartmentalized and non-homologous end joining-dependent repair of double-stranded DNA breaks at the nuclear periphery. Collectively, our data demonstrate that NUMEN generates short 5' overhangs through its endonuclease and 3'→5' exonuclease activities, promotes the repair of DNA lesions-including heterochromatic lamina-associated domain breaks as well as deprotected telomeres-and functions as a downstream effector of DNA-dependent protein kinase catalytic subunit. These findings underline the role of NUMEN as a key player in DNA repair pathway choice and genome-stability maintenance, and have implications for ongoing research into the development and treatment of genome instability disorders.

Fold-back mechanism originating inv-dup-del rearrangements in chromosomes 13 and 15

Intrachromosomal rearrangements involve a single chromosome and can be formed by several proposed mechanisms. We reported two patients with intrachromosomal duplications and deletions, whose rearrangements and breakpoints were characterized through karyotyping, chromosomal microarray, fluorescence in situ hybridization, whole-genome sequencing, and Sanger sequencing. Inverted duplications associated with terminal deletions, known as inv-dup-del rearrangements, were found in 13q and 15q in these patients. The presence of microhomology at the junction points led to the proposal of the Fold-back mechanism for their formation. The use of different high-resolution techniques allowed for a better characterization of the rearrangements, with Sanger sequencing of the junction points being essential to infer the mechanisms of formation as it revealed microhomologies that were missed by the previous techniques. A karyotype-phenotype correlation was also performed for the characterized rearrangements.

Late-Onset Autosomal Dominant Macular Degeneration Caused by Deletion of the CRX Gene

PURPOSE

To characterize the phenotype observed in a case series with macular disease and determine the cause.

DESIGN

Multicenter case series.

PARTICIPANTS

Six families (7 patients) with sporadic or multiplex macular disease with onset at 20 to 78 years, and 1 patient with age-related macular degeneration.

METHODS

Patients underwent ophthalmic examination; exome, genome, or targeted sequencing; and/or polymerase chain reaction (PCR) amplification of the breakpoint, followed by cloning and Sanger sequencing or direct Sanger sequencing.

MAIN OUTCOME MEASURES

Clinical phenotypes, genomic findings, and a hypothesis explaining the mechanism underlying disease in these patients.

RESULTS

All 8 cases carried the same deletion encompassing the genes TPRX1, CRX, and SULT2A1, which was absent from 382 control individuals screened by breakpoint PCR and 13 096 Clinical Genetics patients with a range of other inherited conditions screened by array comparative genomic hybridization. Microsatellite genotypes showed that these 7 families are not closely related, but genotypes immediately adjacent to the deletion breakpoints suggest they may share a distant common ancestor.

CONCLUSIONS

Previous studies had found that carriers for a single defective CRX allele that was predicted to produce no functional CRX protein had a normal ocular phenotype. Here, we show that CRX whole-gene deletion in fact does cause a dominant late-onset macular disease.

Prenatally diagnosed 16p11.2 copy number variations by SNP Array: A retrospective case series

OBJECTIVE

The 16p11.2 copy number variations (CNVs) are increasingly recognized as one of the most frequent genomic disorders, with a broad spectrum of phenotypes. The fetal phenotype associated with 16p11.2 CNVs is poorly described. The current study presents prenatal series of 16p11.2 CNVs and provides a better understanding of this submicroscopic imbalance in prenatal diagnosis.

METHOD

Retrospective case series were extracted from a single tertiary referral center performing prenatal single nucleotide polymorphism (SNP) array from April 2017 to December 2021. The maternal demographics, indication for amniocentesis, ultrasound findings, SNP array results, inheritance of the CNVs, and pregnancy outcomes were studied.

RESULTS

We indentified 30 fetuses carrying 16p11.2 CNVs, representing 0.35% (30/8578) of prenatal SNP array results. The series included 17 fetuses with a proximal deletion, 7 with a distal deletion, 4 with a proximal duplication, and 2 with a distal duplication. Prenatal ultrasound anomalies were reported in 80% of these cases. The most common presentation was vertebralanomalies (9/30). Other features noted in more than one fetus were increased nuchal translucency/nuchal fold (NT/NF) (5/30), absent/hypoplastic nasal bone (3/30), polyhydramnios (3/30), ventricular septal defect (VSD) (2/30), unilateral mild ventriculomegaly (2/30), fetal growth restriction (FGR) (2/30), right aortic arch (2/30). All the 9 vertebralanomalies were present in fetuses harboring proximal deletion (9/17). Familial transmission was confirmed in 44% of cases (11/25) and termination of pregnancy was requested in 62.1% (18/29) of cases.

CONCLUSION

The 16p11.2 CNVs can have variable prenatal phenotypes and these CNVs are frequently inherited from parents with a milder or normal phenotype. Our results underline that vertebral deformities were frequent in cases of 16p11.2 proximal deletion, and further demonstrate the incomplete penetrance of the CNVs.

Comprehensive Analysis of Copy Number Variations on Glycoside Hydrolase 45 Genes among Different Bursaphelenchus xylophilus Strains

Bursaphelenchus xylophilus is considered the most dangerous quarantine pest in China. It causes enormous economic and ecological losses in many countries from Asia and Europe. The glycoside hydrolase 45 gene family has been demonstrated in early studies to contribute to the cell wall degradation ability of B. xylophilus during its infection. However, the copy number variation (CNV) of the GH45 gene and its association with B. xylophilus pathogenicity were not fully elucidated. In this study, we found that the GH45 gene with two copies is the most predominant type among 259 B. xylophilus strains collected from China and Japan. Additionally, 18 strains are identified as GH45 genes with a single copy, and only two strains are verified to have three copies. Subsequent expression analysis and inoculation test suggest that the copy numbers of the GH45 gene are correlated with gene expression as well as the B. xylophilus pathogenicity. B. xylophilus strains with more copies of the GH45 gene usually exhibit more abundant expression and cause more severe wilt symptoms on pine trees. The aforementioned results indicated the potential regulatory effects of CNV in B. xylophilus and provided novel information to better understand the molecular pathogenesis of this devastating pest.

Reversibility and therapeutic development for neurodevelopmental disorders, insights from genetic animal models

Neurodevelopmental Disorders (NDDs) encompass a broad spectrum of conditions resulting from atypical brain development. Over the past decades, we have had the fortune to witness enormous progress in diagnosis, etiology discovery, modeling, and mechanistic understanding of NDDs from both fundamental and clinical research. Here, we review recent neurobiological advances from experimental models of NDDs. We introduce several examples and highlight breakthroughs in reversal studies of phenotypes using genetically engineered models of NDDs. The in-depth understanding of brain pathophysiology underlying NDDs and evaluations of reversibility in animal models paves the foundation for discovering novel treatment options. We discuss how the expanding property of cutting-edge technologies, such as gene editing and AAV-mediated gene delivery, are leveraged in animal models for the therapeutic development of NDDs. We envision opportunities and challenges toward faithful modeling and fruitful clinical translation.

Genes that are Used Together are More Likely to be Fused Together in Evolution by Mutational Mechanisms: A Bioinformatic Test of the Used-Fused Hypothesis

Cases of parallel or recurrent gene fusions in evolution as well as in genetic disease and cancer are difficult to explain, because unlike point mutations, they can require the repetition of a similar configuration of multiple breakpoints rather than the repetition of a single point mutation. The used-together-fused-together hypothesis holds that genes that are used together repeatedly and persistently in a specific context are more likely to undergo fusion mutation in the course of evolution for mechanistic reasons. This hypothesis offers to explain gene fusion in both evolution and disease under one umbrella. Using bioinformatic data, we tested this hypothesis against alternatives, including that all gene pairs can fuse by random mutation, but among pairs thus fused, those that had interacted previously are more likely to be favored by selection. Results show that across multiple measures of gene interaction, human genes whose orthologs are fused in one or more species are more likely to interact with each other than random pairs of genes of the same genomic distance between pair members; that an overlap exists between genes that fused in the course of evolution in non-human species and genes that undergo fusion in human cancers; and that across six primate species studied, fusions predominate over fissions and exhibit substantial evolutionary parallelism. Together, these results support the used-together-fused-together hypothesis over its alternatives. Multiple implications are discussed, including the relevance of mutational mechanisms to the evolution of genome organization, to the distribution of fitness effects of mutation, to evolutionary parallelism and more.

Supplementary Information

The online version contains supplementary material available at 10.1007/s11692-022-09579-9.

Complex genomic rearrangements: an underestimated cause of rare diseases

Complex genomic rearrangements (CGRs) are known contributors to disease but are often missed during routine genetic screening. Identifying CGRs requires (i) identifying copy number variants (CNVs) concurrently with inversions, (ii) phasing multiple breakpoint junctions incis, as well as (iii) detecting and resolving structural variants (SVs) within repeats. We demonstrate how combining cytogenetics and new sequencing methodologies is being successfully applied to gain insights into the genomic architecture of CGRs. In addition, we review CGR patterns and molecular features revealed by studying constitutional genomic disorders. These data offer invaluable lessons to individuals interested in investigating CGRs, evaluating their clinical relevance and frequency, as well as assessing their impact(s) on rare genetic diseases.

Developmental genomics of limb malformations: Allelic series in association with gene dosage effects contribute to the clinical variability

Genetic heterogeneity, reduced penetrance, and variable expressivity, the latter including asymmetric body axis plane presentations, have all been described in families with congenital limb malformations (CLMs). Interfamilial and intrafamilial heterogeneity highlight the complexity of the underlying genetic pathogenesis of these developmental anomalies. Family-based genomics by exome sequencing (ES) and rare variant analyses combined with whole-genome array-based comparative genomic hybridization were implemented to investigate 18 families with limb birth defects. Eleven of 18 (61%) families revealed explanatory variants, including 7 single-nucleotide variant alleles and 3 copy number variants (CNVs), at previously reported "disease trait associated loci": BHLHA9, GLI3, HOXD cluster, HOXD13, NPR2, and WNT10B. Breakpoint junction analyses for all three CNV alleles revealed mutational signatures consistent with microhomology-mediated break-induced replication, a mechanism facilitated by Alu/Alu-mediated rearrangement. Homozygous duplication of BHLHA9 was observed in one Turkish kindred and represents a novel contributory genetic mechanism to Gollop-Wolfgang Complex (MIM: 228250), where triplication of the locus has been reported in one family from Japan (i.e., 4n = 2n + 2n versus 4n = 3n + 1n allelic configurations). Genes acting on limb patterning are sensitive to a gene dosage effect and are often associated with an allelic series. We extend an allele-specific gene dosage model to potentially assist, in an adjuvant way, interpretations of interconnections among an allelic series, clinical severity, and reduced penetrance of the BHLHA9-related CLM spectrum.

Chromosomal microarray analyses from 5778 patients with neurodevelopmental disorders and congenital anomalies in Brazil

Chromosomal microarray analysis (CMA) has been recommended and practiced routinely since 2010 both in the USA and Europe as the first-tier cytogenetic test for patients with unexplained neurodevelopmental delay/intellectual disability, autism spectrum disorders, and/or multiple congenital anomalies. However, in Brazil, the use of CMA is still limited, due to its high cost and complexity in integrating the results from both the private and public health systems. Although Brazil has one of the world’s largest single-payer public healthcare systems, nearly all patients referred for CMA come from the private sector, resulting in only a small number of CMA studies in Brazilian cohorts. To date, this study is by far the largest Brazilian cohort (n = 5788) studied by CMA and is derived from a joint collaboration formed by the University of São Paulo and three private genetic diagnostic centers to investigate the genetic bases of neurodevelopmental disorders and congenital abnormalities. We identified 2,279 clinically relevant CNVs in 1886 patients, not including the 26 cases of UPD found. Among detected CNVs, the corresponding frequency of each category was 55.6% Pathogenic, 4.4% Likely Pathogenic and 40% VUS. The diagnostic yield, by taking into account Pathogenic, Likely Pathogenic and UPDs, was 19.7%. Since the rational for the classification is mostly based on Mendelian or highly penetrant variants, it was not surprising that a second event was detected in 26% of those cases of predisposition syndromes. Although it is common practice to investigate the inheritance of VUS in most laboratories around the world to determine the inheritance of the variant, our results indicate an extremely low cost–benefit of this approach, and strongly suggest that in cases of a limited budget, investigation of the parents of VUS carriers using CMA should not be prioritized.

A missense mutation in TTC8/BBS8 affecting mRNA splicing in patients with non-syndromic retinitis pigmentosa

Splicing disruption is one type of mutation mechanism for disease-predisposing alleles. To date, less than 30 mutations in TTC8/BBS8 have been reported; however, mutations affecting the splice site are rare. Generally missense mutations are assumed to alter protein function; however, reports have shown that mutations in protein coding exons can disrupt splicing by altering exonic splicing silencer or enhancer motifs. Hence, a missense mutation c.1347G > C (p.Q449H) involving final base of the exon 13 in the TTC8, previously identified by us to be linked with non-syndromic autosomal recessive retinitis pigmentosa (arRP), in an Indian family, that might deleteriously affect splicing has been functionally characterized. RNA was isolated, cDNA prepared and amplified using region-specific primers. PCR products were purified and sequenced bi-directionally by Sanger sequencing. Effect of mutation (c.1347G > C) on mRNA splicing has been predicted using bioinformatics tools. We reported that missense mutation (c.1347G > C) at the last base of exon 13 of TTC8 disrupted the canonical donor splice-site resulting in aberrant RNA splicing. A cryptic donor splice-site got activated 77 bases downstream of the authentic splice donor site in intron 13, resulting in the retention of 77 bases of intron 13, and a frameshift leading to pre-mature termination codon in exon 14 at codon 486. Further, duplication of exon 15 and fusion of its duplicated copy occurred with exon 13. The binding site for SC35 protein, normally involved in splicing, also got disrupted (as predicted by SpliceAid2 software), hence, leading to alternative splicing. Our findings strongly suggest that a missense mutation c.1347G > C in TTC8 disrupted the splice donor site causing retention of 77 bases of intron 13, resulting in a frameshift and subsequently introduced a pre-mature termination codon into exon 14, hence creating an altered mRNA transcript. These findings emphasize the significance of examining missense mutations especially in TTC8, to determine their pathogenic role through alternative splicing. Present findings also reiterate the notion that mutations in the TTC8/BBS8 cause phenotypic heterogeneity and does not always follow Mendelian genetics in this ciliopathy.

Comparative analysis of transient receptor potential channel 5 opposite strand-induced gene expression patterns and protein-protein interactions in triple-negative breast cancer

In patients with triple-negative breast cancer (TNBC), high tumour mutation burden and aberrant oncogene expression profiles are some of the causes of poor prognosis. Therefore, it is necessary to identify aberrantly expressed oncogenes, since they have the potential to serve as therapeutic targets. Transient receptor potential channel 5 opposite strand (TRPC5OS) has been previously shown to function as a novel tumour inducer. However, the underlying mechanism of TRPC5OS function in TNBC remain to be elucidated. Therefore, in the present study TRPC5OS expression was first measured in tissue samples of patients with TNBC and a panel of breast cancer cell lines (ZR-75-1, MDA-MB-453, SK-BR-3, JIMT-1, BT474 and HCC1937) by using qRT-PCR and Western blotting. Subsequently, the possible effects of TRPC5OS on MDA-MB-231 cells proliferation were determined using Cell Counting Kit-8 and 5-Ethynyl-2′-deoxyuridine assays after Lentiviral transfection of MDA-MB-231. In addition, potential interaction partners of TRPC5OS were explored using liquid chromatography-mass spectrometry (LC-MS)/MS. Gene expression patterns following TRPC5OS overexpression were also detected in MDA-MB-231 cells by using High-throughput sequencing. Gene Ontology (GO) and Kyoto Encyclopaedia of Genes and Genomes (KEGG) analysis were then used to systematically verify the potential interactions among the TRPC5OS-regulated genes. The potential relationship between TRPC5OS-interacting proteins and gene expression patterns were studied using Search Tool for the Retrieval of Interacting Genes/Proteins (STRING) analysis. TRPC5OS expression was found to be significantly higher in TNBC tumour tissues and breast cancer cell lines compared with luminal tumour tissues and ZR-75-1. In addition, the overexpression of TRPC5OS significantly increased cell proliferation. High-throughput sequencing results revealed that 5,256 genes exhibited differential expression following TRPC5OS overexpression, including 3,269 upregulated genes and 1,987 downregulated genes. GO analysis results indicated that the functions of these differentially expressed genes were enriched in the categories of ‘cell division’ and ‘cell proliferation’ regulation. KEGG analysis showed that the TRPC5OS-regulated genes were associated with processes of ‘homologous recombination’ and ‘TNF signalling pathways’. Subsequently, 17 TRPC5OS-interacting proteins were found using LC-MS/MS and STRING analysis. The most important protein among interacting proteins was ENO1 which was associated with glycolysis and regulated proliferation of cancer. In summary, data from the present study suggest that TRPC5OS overexpression can increase TNBC cell proliferation and ENO1 may be a potential target protein mediated by TRPC5OS. Therefore, TRPC5OS may serve as a novel therapeutic target for TNBC.

Dynamic Features of Chromosomal Instability during Culture of Induced Pluripotent Stem Cells

Induced pluripotent stem cells (iPSCs) hold great potential for regenerative medicine. By reprogramming a patient′s own cells, immunological rejection can be avoided during transplantation. For expansion and gene editing, iPSCs are grown in artificial culture for extended times. Culture affords potential danger for the accumulation of genetic aberrations. To study these, two induced pluripotent stem (iPS) cell lines were cultured and periodically analyzed using advanced optical mapping to detect and classify chromosome numerical and segmental changes that included deletions, insertions, balanced translocations and inversions. In one of the lines, a population trisomic for chromosome 12 gained dominance over a small number of passages. This appearance and dominance of the culture by chromosome 12 trisomic cells was tracked through intermediate passages by the analysis of chromosome spreads. Mathematical modeling suggested that the proliferation rates of diploid versus trisomic cells could not account for the rapid dominance of the trisomic population. In addition, optical mapping revealed hundreds of structural variations distinct from those generally found within the human population. Many of these structural variants were detected in samples obtained early in the culturing process and were maintained in late passage samples, while others were acquired over the course of culturing.

Mechanisms of structural chromosomal rearrangement formation

Structural chromosomal rearrangements result from different mechanisms of formation, usually related to certain genomic architectural features that may lead to genetic instability. Most of these rearrangements arise from recombination, repair, or replication mechanisms that occur after a double-strand break or the stalling/breakage of a replication fork. Here, we review the mechanisms of formation of structural rearrangements, highlighting their main features and differences. The most important mechanisms of constitutional chromosomal alterations are discussed, including Non-Allelic Homologous Recombination (NAHR), Non-Homologous End-Joining (NHEJ), Fork Stalling and Template Switching (FoSTeS), and Microhomology-Mediated Break-Induced Replication (MMBIR). Their involvement in chromoanagenesis and in the formation of complex chromosomal rearrangements, inverted duplications associated with terminal deletions, and ring chromosomes is also outlined. We reinforce the importance of high-resolution analysis to determine the DNA sequence at, and near, their breakpoints in order to infer the mechanisms of formation of structural rearrangements and to reveal how cells respond to DNA damage and repair broken ends.

SHIMS 3.0: Highly efficient single-haplotype iterative mapping and sequencing using ultra-long nanopore reads

The reference sequence of structurally complex regions can only be obtained through a highly accurate clone-based approach that we call Single-Haplotype Iterative Mapping and Sequencing (SHIMS). In recent years, improvements to SHIMS have reduced the cost and time required by two orders of magnitude, but internally repetitive clones still require extensive manual effort to transform draft assemblies into reference-quality finished sequences. Here we describe SHIMS 3.0, using ultra-long nanopore reads to augment the Illumina data from SHIMS 2.0 assemblies and resolve internally repetitive structures. This greatly minimizes the need for manual finishing of Illumina-based draft assemblies, allowing a small team with no prior finishing experience to sequence challenging targets with high accuracy. This protocol proceeds from clone-picking to finished assemblies in 2 weeks for about $80 (USD) per clone. We recently used this protocol to produce reference sequence of structurally complex palindromes on chimpanzee and rhesus macaque X chromosomes. Our protocol provides access to structurally complex regions that would otherwise be inaccessible from whole-genome shotgun data or require an impractical amount of manual effort to generate an accurate assembly.

16p13.11p11.2 triplication syndrome: a new recognizable genomic disorder characterized by optical genome mapping and whole genome sequencing

Highly identical segmental duplications (SDs) account for over 5% of the human genome and are enriched in the short arm of the chromosome 16. These SDs are susceptibility factors for recurrent chromosomal rearrangements mediated by non-allelic homologous recombination (NAHR). Chromosomal microarray analysis (CMA) has been widely used as the first-tier test for individuals with developmental disabilities and/or congenital anomalies and several genomic disorders involving the 16p-arm have been identified with this technique. However, the resolution of CMA and the limitations of short-reads whole genome sequencing (WGS) technology do not allow the full characterization of the most complex chromosomal rearrangements. Herein, we report on two unrelated patients with a de novo 16p13.11p11.2 triplication associated with a 16p11.2 duplication, detected by CMA. These patients share a similar phenotype including hypotonia, severe neurodevelopmental delay with profound speech impairment, hyperkinetic behavior, conductive hearing loss, and distinctive facial features. Short-reads WGS could not map precisely any of the rearrangement's breakpoints that lie within SDs. We used optical genome mapping (OGM) to determine the relative orientation of the triplicated and duplicated segments as well as the genomic positions of the breakpoints, allowing us to propose a mechanism involving recombination between allelic SDs and a NAHR event. In conclusion, we report a new clinically recognizable genomic disorder. In addition, the mechanism of these complex chromosomal rearrangements involving SDs could be unraveled by OGM.

Biology in balance: human diploid genome integrity, gene dosage, and genomic medicine

The path to completion of the functional annotation of the haploid human genome reference build, exploration of the clan genomics hypothesis, understanding human gene and genome functional biology, and gene genome and organismal evolution, is in reach.

Recurrent inversion polymorphisms in humans associate with genetic instability and genomic disorders

Unlike copy number variants (CNVs), inversions remain an underexplored genetic variation class. By integrating multiple genomic technologies, we discover 729 inversions in 41 human genomes. Approximately 85% of inversions <2 kbp form by twin-priming during L1 retrotransposition; 80% of the larger inversions are balanced and affect twice as many nucleotides as CNVs. Balanced inversions show an excess of common variants, and 72% are flanked by segmental duplications (SDs) or retrotransposons. Since flanking repeats promote non-allelic homologous recombination, we developed complementary approaches to identify recurrent inversion formation. We describe 40 recurrent inversions encompassing 0.6% of the genome, showing inversion rates up to 2.7 × 10-4 per locus per generation. Recurrent inversions exhibit a sex-chromosomal bias and co-localize with genomic disorder critical regions. We propose that inversion recurrence results in an elevated number of heterozygous carriers and structural SD diversity, which increases mutability in the population and predisposes specific haplotypes to disease-causing CNVs.

Characterization of Chromosomal Breakpoints in 12 Cases with 8p Rearrangements Defines a Continuum of Fragility of the Region

Improvements in microarray-based comparative genomic hybridization technology have allowed for high-resolution detection of genome wide copy number alterations, leading to a better definition of rearrangements and supporting the study of pathogenesis mechanisms. In this study, we focused our attention on chromosome 8p. We report 12 cases of 8p rearrangements, analyzed by molecular karyotype, evidencing a continuum of fragility that involves the entire short arm. The breakpoints seem more concentrated in three intervals: one at the telomeric end, the others at 8p23.1, close to the beta-defensin gene cluster and olfactory receptor low-copy repeats. Hypothetical mechanisms for all cases are described. Our data extend the cohort of published patients with 8p aberrations and highlight the need to pay special attention to these sequences due to the risk of formation of new chromosomal aberrations with pathological effects.

Xq27.1 palindrome mediated interchromosomal insertion likely causes familial congenital bilateral laryngeal abductor paralysis (Plott syndrome)

Bilateral laryngeal abductor paralysis is a rare entity and the second most common cause of stridor in newborns. So far, no conclusive genetic or chromosomal aberration has been reported for X-linked isolated bilateral vocal cord paralysis, also referred to as Plott syndrome. Via whole genome sequencing (WGS), we identified a complex interchromosomal insertion in a large family with seven affected males. The 404 kb inserted fragment originates from chromosome 10q21.3, contains no genes and is inserted inversionally into the intergenic chromosomal region Xq27.1, 82 kb centromeric to the nearest gene SOX3. The patterns found at the breakpoint junctions resemble typical characteristics that arise in replication-based mechanisms with long-distance template switching. Non protein-coding insertions into the same genomic region have been described to result in different phenotypes, indicating that the phenotypic outcome likely depends on the introduction of regulatory elements. In conclusion, our data adds Plott syndrome as another entity, likely caused by the insertion of non-coding DNA into the intergenic chromosomal region Xq27.1. In this regard, we demonstrate the importance of WGS as a powerful diagnostic test in unsolved genetic diseases, as this genomic rearrangement has not been detected by current first-line diagnostic tests, i.e., exome sequencing and chromosomal microarray analysis.

De novo mutation rates at the single-mutation resolution in a human HBB gene-region associated with adaptation and genetic disease

Although it is known that the mutation rate varies across the genome, previous estimates were based on averaging across various numbers of positions. Here, we describe a method to measure the origination rates of target mutations at target base positions and apply it to a 6-bp region in the human hemoglobin subunit beta (HBB) gene and to the identical, paralogous hemoglobin subunit delta (HBD) region in sperm cells from both African and European donors. The HBB region of interest (ROI) includes the site of the hemoglobin S (HbS) mutation, which protects against malaria, is common in Africa, and has served as a classic example of adaptation by random mutation and natural selection. We found a significant correspondence between de novo mutation rates and past observations of alleles in carriers, showing that mutation rates vary substantially in a mutation-specific manner that contributes to the site frequency spectrum. We also found that the overall point mutation rate is significantly higher in Africans than in Europeans in the HBB region studied. Finally, the rate of the 20A→T mutation, called the “HbS mutation” when it appears in HBB, is significantly higher than expected from the genome-wide average for this mutation type. Nine instances were observed in the African HBB ROI, where it is of adaptive significance, representing at least three independent originations; no instances were observed elsewhere. Further studies will be needed to examine mutation rates at the single-mutation resolution across these and other loci and organisms and to uncover the molecular mechanisms responsible.

Mechanism for inverted-repeat recombination induced by a replication fork barrier

Replication stress and abundant repetitive sequences have emerged as primary conditions underlying genomic instability in eukaryotes. To gain insight into the mechanism of recombination between repeated sequences in the context of replication stress, we used a prokaryotic Tus/Ter barrier designed to induce transient replication fork stalling near inverted repeats in the budding yeast genome. Our study reveals that the replication fork block stimulates a unique recombination pathway dependent on Rad51 strand invasion and Rad52-Rad59 strand annealing activities, Mph1/Rad5 fork remodelers, Mre11/Exo1/Dna2 resection machineries, Rad1-Rad10 nuclease and DNA polymerase δ. Furthermore, we show recombination at stalled replication forks is limited by the Srs2 helicase and Mus81-Mms4/Yen1 nucleases. Physical analysis of the replication-associated recombinants revealed that half are associated with an inversion of sequence between the repeats. Based on our extensive genetic characterization, we propose a model for recombination of closely linked repeats that can robustly generate chromosome rearrangements.

Replication stress and abundant repetitive sequences have emerged as primary conditions underlying genomic instability in eukaryotes. Here the authors use a prokaryotic Tus/Ter barrier designed to induce transient replication fork stalling near inverted repeats in the budding yeast genome to support a model for recombination of closely linked repeats at stalled replication forks.

Establishment of a Dihydrofolate Reductase Gene Knock-In Zebrafish Strain to Aid Preliminary Analysis of Congenital Heart Disease Mechanisms

Background: The dihydrofolate reductase (DHFR) gene is imperative in development, therefore it is essential to explore its effects on heart development. Thus, here a dhfr zebrafish knock-in (KI) strain was constructed.

Methods: CRISPR/Cas9 technology was used to establish the dhfr KI zebrafish strain. This strain was hybridized with TgG fluorescent strain zebrafish to observe the phenotypes of heart shape, size, and circularization direction. Wild-type (WT) and KI zebrafish were then dissected and histologically stained to observe pathological changes. Western blot analysis was used to verify the increased expressions of zebrafish genes after KI. Hybridization experiments were used to confirm the presence of abnormal gonadal dysplasia.

Results: The zebrafish dhfr KI strain was successfully constructed through CRISPR/Cas9 technology. At 6 days post fertilization (dpf), microscopic examinations of KI (homozygous) specimens revealed pericardial effusions, heart compressions, and curled tails. Compared with WT, the Hematoxylin and Eosin (H&E) tissue sections of KI-homozygous zebrafish showed defects such as reduced atria and ventricles. Western blot analysis indicated that the expression of the DHFR protein increased in both heterozygotes and homozygotes of dhfr KI zebrafish. Hybridization experiments revealed that dhfr KI may affect gonadal function.

Conclusion: The DHFR gene plays an important regulatory role in the process of heart development, and copy number variations (CNVs) of this gene may constitute a new pathogenic mechanism of congenital heart disease (CHD).

Oxytocin and Oxytocin Receptor Gene Regulation in Williams Syndrome: A Systematic Review

Williams Syndrome (WS) is a rare genetic multisystem disorder that occurs because of a deletion of approximately 25 genes in the 7q11.23 chromosome region. This causes dysmorphic facial appearances, multiple congenital cardiovascular defects, delayed motor skills, and abnormalities in connective tissues and the endocrine system. The patients are mostly diagnosed with mild to moderate mental retardation, however, they have a hyper sociable, socially dis-inhibited, and outgoing personality, empathetic behavior, and are highly talkative. Oxytocin (OT), a neuropeptide synthesized at the hypothalamus, plays an important role in cognition and behavior, and is thought to be affecting WS patients' attitudes at its different amounts. Oxytocin receptor gene (OXTR), on chromosome 3p25.3, is considered regulating oxytocin receptors, via which OT exerts its effect. WS is a crucial disorder to understand gene, hormone, brain, and behavior associations in terms of sociality and neuropsychiatric conditions. Alterations to the WS gene region offer an opportunity to deepen our understandings of autism spectrum disorder, schizophrenia, anxiety, or depression. We aim to systematically present the data available of OT/OXTR regulation and expression, and the evidence for whether these mechanisms are dysregulated in WS. These results are important, as they predict strong epigenetic control over social behavior by methylation, single nucleotide polymorphisms, and other alterations. The comparison and collaboration of these studies may help to establish a better treatment or management approach for patients with WS if backed up with future research.

Clan genomics: From OMIM phenotypic traits to genes and biology

Clinical characterization of a patient phenotype has been the quintessential approach for elucidating a differential diagnosis and a hypothesis to explore a potential clinical diagnosis. This has resulted in a language of medicine and a semantic ontology, with both specialty- and subspecialty-specific lexicons, that can be challenging to translate and interpret. There is no 'Rosetta Stone' of clinical medicine such as the genetic code that can assist translation and interpretation of the language of genetics. Nevertheless, the information content embodied within a clinical diagnosis can guide management, therapeutic intervention, and potentially prognostic outlook of disease enabling anticipatory guidance for patients and families. Clinical genomics is now established firmly in medical practice. The granularity and informative content of a personal genome is immense. Yet, we are limited in our utility of much of that personal genome information by the lack of functional characterization of the overwhelming majority of computationally annotated genes in the haploid human reference genome sequence. Whereas DNA and the genetic code have provided a 'Rosetta Stone' to translate genetic variant information, clinical medicine, and clinical genomics provide the context to understand human biology and disease. A path forward will integrate deep phenotyping, such as available in a clinical synopsis in the Online Mendelian Inheritance in Man (OMIM) entries, with personal genome analyses.

Long-read technologies identify a hidden inverted duplication in a family with choroideremia

The lack of molecular diagnoses in rare genetic diseases can be explained by limitations of current standard genomic technologies. Upcoming long-read techniques have complementary strengths to overcome these limitations, with a particular strength in identifying structural variants. By using optical genome mapping and long-read sequencing, we aimed to identify the pathogenic variant in a large family with X-linked choroideremia. In this family, aberrant splicing of exon 12 of the choroideremia gene CHM was detected in 2003, but the underlying genomic defect remained elusive. Optical genome mapping and long-read sequencing approaches now revealed an intragenic 1,752 bp inverted duplication including exon 12 and surrounding regions, located downstream of the wild-type copy of exon 12. Both breakpoint junctions were confirmed with Sanger sequencing and segregate with the X-linked inheritance in the family. The breakpoint junctions displayed sequence microhomology suggestive for an erroneous replication mechanism as the origin of the structural variant. The inverted duplication is predicted to result in a hairpin formation of the pre-mRNA with the wild-type exon 12, leading to exon skipping in the mature mRNA. The identified inverted duplication is deemed the hidden pathogenic cause of disease in this family. Our study shows that optical genome mapping and long-read sequencing have significant potential for the identification of (hidden) structural variants in rare genetic diseases.

Genomic regions associated with microdeletion/microduplication syndromes exhibit extreme diversity of structural variation

Segmental duplications (SDs) are a class of long, repetitive DNA elements whose paralogs share a high level of sequence similarity with each other. SDs mediate chromosomal rearrangements that lead to structural variation in the general population as well as genomic disorders associated with multiple congenital anomalies, including the 7q11.23 (Williams-Beuren Syndrome, WBS), 15q13.3, and 16p12.2 microdeletion syndromes. Population-level characterization of SDs has generally been lacking because most techniques used for analyzing these complex regions are both labor and cost intensive. In this study, we have used a high-throughput technique to genotype complex structural variation with a single molecule, long-range optical mapping approach. We characterized SDs and identified novel structural variants (SVs) at 7q11.23, 15q13.3, and 16p12.2 using optical mapping data from 154 phenotypically normal individuals from 26 populations comprising five super-populations. We detected several novel SVs for each locus, some of which had significantly different prevalence between populations. Additionally, we localized the microdeletion breakpoints to specific paralogous duplicons located within complex SDs in two patients with WBS, one patient with 15q13.3, and one patient with 16p12.2 microdeletion syndromes. The population-level data presented here highlights the extreme diversity of large and complex SVs within SD-containing regions. The approach we outline will greatly facilitate the investigation of the role of inter-SD structural variation as a driver of chromosomal rearrangements and genomic disorders.

Clinical and genomic characterization of 8p cytogenomic disorders

PURPOSE

To provide a detailed clinical and cytogenomic summary of individuals with chromosome 8p rearrangements of invdupdel(8p), del(8p), and dup(8p).

METHODS

We enrolled 97 individuals with invdupdel(8p), del(8p), and dup(8p). Clinical and molecular data were collected to delineate and compare the clinical findings and rearrangement breakpoints. We included additional 5 individuals with dup(8p) from the literature for a total of 102 individuals.

RESULTS

Eighty-one individuals had recurrent rearrangements of invdupdel(8p) (n = 49), del(8p)_distal (n = 4), del(8p)_proximal (n = 9), del(8p)_proximal&distal (n = 12), and dup(8p)_proximal (n = 7). Twenty-one individuals had nonrecurrent rearrangements. While all individuals had neurodevelopmental features, the frequency and severity of clinical findings were higher in individuals with invdupdel(8p), and with larger duplications. All individuals with GATA4 deletion had structural congenital heart defects; however, the presence of structural heart defects in some individuals with normal GATA4 copy number suggests there are other potentially contributing gene(s) on 8p.

CONCLUSION

Our study may inform families and health-care providers about the associated clinical findings and severity in individuals with chromosome 8p rearrangements, and guide researchers in investigating the underlying molecular and biological mechanisms by providing detailed clinical and cytogenomic information about individuals with distinct 8p rearrangements.

Copy number variation and neuropsychiatric illness

Copy number variants (CNVs) at specific loci have been identified as important risk factors for several neuropsychiatric disorders, such as schizophrenia, autism spectrum disorder, intellectual disability (ID) and depression. These CNVs are individually rare (<0.5% frequency), have high effect sizes, and show pleiotropic effects for multiple neuropsychiatric disorders, which implies a shared aetiology. Neuropsychiatric CNVs are also associated with cognitive impairment and other medical morbidities, such as heart defects and obesity. As most neuropsychiatric CNVs are multigenic, it has been challenging to map their effects onto specific biological processes, although gene-set analyses have implicated genes related to the synapse and chromatin regulation. However, future whole-genome sequencing studies have potential for identifying novel single-gene CNV associations, which could provide insights into the pathophysiology underlying neuropsychiatric disorders.