Predicting human genes susceptible to genomic instability associated with Alu/Alu-mediated rearrangements

Compared to other NHEJ factors, DNA-PK protein and RNA levels are markedly increased in all higher primates, but not in prosimians or other mammals

The DNA dependent protein kinase (DNA-PK) initiates non-homologous recombination (NHEJ), the predominate DNA double-strand break (DSBR) pathway in higher vertebrates. It has been known for decades that the enzymatic activity of DNA-PK [that requires its three component polypeptides, Ku70, Ku80 (that comprise the DNA-end binding Ku heterodimer), and the catalytic subunit (DNA-PKcs)] is present in humans at 10-50 times the level observed in other mammals. Here, we show that the high level of DNA-PKcs protein expression appears evolutionarily in mammals between prosimians and higher primates. Moreover, the RNAs encoding the three component polypeptides of DNA-PK are present at similarly high levels in hominids, new-, and old-world monkeys, but expression of these RNAs in prosimians is ∼5-50 fold less, analogous to the levels observed in other non-primate species. This is reminiscent of the appearance of Alu repeats in primate genomes -- abundant in higher primates, but present at much lower density in prosimians. Alu repeats are well-known for their capacity to promote non-allelic homologous recombination (NAHR) a process known to be inhibited by DNA-PK. Nanopore sequence analyses of cultured cells proficient or deficient in DNA-PK revealed an increase of inter-chromosomal translocations caused by NAHR. Although the high levels of DNA-PK in primates may have many functions, we posit that high levels of DNA-PK may function to restrain deleterious NAHR events between Alu elements.

Inverted triplications formed by iterative template switches generate structural variant diversity at genomic disorder loci

The duplication-triplication/inverted-duplication (DUP-TRP/INV-DUP) structure is a complex genomic rearrangement (CGR). Although it has been identified as an important pathogenic DNA mutation signature in genomic disorders and cancer genomes, its architecture remains unresolved. Here, we studied the genomic architecture of DUP-TRP/INV-DUP by investigating the DNA of 24 patients identified by array comparative genomic hybridization (aCGH) on whom we found evidence for the existence of 4 out of 4 predicted structural variant (SV) haplotypes. Using a combination of short-read genome sequencing (GS), long-read GS, optical genome mapping, and single-cell DNA template strand sequencing (strand-seq), the haplotype structure was resolved in 18 samples. The point of template switching in 4 samples was shown to be a segment of ∼2.2-5.5 kb of 100% nucleotide similarity within inverted repeat pairs. These data provide experimental evidence that inverted low-copy repeats act as recombinant substrates. This type of CGR can result in multiple conformers generating diverse SV haplotypes in susceptible dosage-sensitive loci.

Novel and Recurrent Copy Number Variants in ABCA4-Associated Retinopathy

ABCA4 is the most frequently mutated gene leading to inherited retinal disease (IRD) with over 2200 pathogenic variants reported to date. Of these, ~1% are copy number variants (CNVs) involving the deletion or duplication of genomic regions, typically >50 nucleotides in length. An in-depth assessment of the current literature based on the public database LOVD, regarding the presence of known CNVs and structural variants in ABCA4, and additional sequencing analysis of ABCA4 using single-molecule Molecular Inversion Probes (smMIPs) for 148 probands highlighted recurrent and novel CNVs associated with ABCA4-associated retinopathies. An analysis of the coverage depth in the sequencing data led to the identification of eleven deletions (six novel and five recurrent), three duplications (one novel and two recurrent) and one complex CNV. Of particular interest was the identification of a complex defect, i.e., a 15.3 kb duplicated segment encompassing exon 31 through intron 41 that was inserted at the junction of a downstream 2.7 kb deletion encompassing intron 44 through intron 47. In addition, we identified a 7.0 kb tandem duplication of intron 1 in three cases. The identification of CNVs in ABCA4 can provide patients and their families with a genetic diagnosis whilst expanding our understanding of the complexity of diseases caused by ABCA4 variants.

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.

NODAL variants are associated with a continuum of laterality defects from simple D-transposition of the great arteries to heterotaxy

BACKGROUND

NODAL signaling plays a critical role in embryonic patterning and heart development in vertebrates. Genetic variants resulting in perturbations of the TGF-β/NODAL signaling pathway have reproducibly been shown to cause laterality defects in humans. To further explore this association and improve genetic diagnosis, the study aims to identify and characterize a broader range of NODAL variants in a large number of individuals with laterality defects.

METHODS

We re-analyzed a cohort of 321 proband-only exomes of individuals with clinically diagnosed laterality congenital heart disease (CHD) using family-based, rare variant genomic analyses. To this cohort we added 12 affected subjects with known NODAL variants and CHD from institutional research and clinical cohorts to investigate an allelic series. For those with candidate contributory variants, variant allele confirmation and segregation analysis were studied by Sanger sequencing in available family members. Array comparative genomic hybridization and droplet digital PCR were utilized for copy number variants (CNV) validation and characterization. We performed Human Phenotype Ontology (HPO)-based quantitative phenotypic analyses to dissect allele-specific phenotypic differences.

RESULTS

Missense, nonsense, splice site, indels, and/or structural variants of NODAL were identified as potential causes of heterotaxy and other laterality defects in 33 CHD cases. We describe a recurrent complex indel variant for which the nucleic acid secondary structure predictions implicate secondary structure mutagenesis as a possible mechanism for formation. We identified two CNV deletion alleles spanning NODAL in two unrelated CHD cases. Furthermore, 17 CHD individuals were found (16/17 with known Hispanic ancestry) to have the c.778G > A:p.G260R NODAL missense variant which we propose reclassification from variant of uncertain significance (VUS) to likely pathogenic. Quantitative HPO-based analyses of the observed clinical phenotype for all cases with p.G260R variation, including heterozygous, homozygous, and compound heterozygous cases, reveal clustering of individuals with biallelic variation. This finding provides evidence for a genotypic-phenotypic correlation and an allele-specific gene dosage model.

CONCLUSION

Our data further support a role for rare deleterious variants in NODAL as a cause for sporadic human laterality defects, expand the repertoire of observed anatomical complexity of potential cardiovascular anomalies, and implicate an allele specific gene dosage model.

HMZDupFinder: a robust computational approach for detecting intragenic homozygous duplications from exome sequencing data

Homozygous duplications contribute to genetic disease by altering gene dosage or disrupting gene regulation and can be more deleterious to organismal biology than heterozygous duplications. Intragenic exonic duplications can result in loss-of-function (LoF) or gain-of-function (GoF) alleles that when homozygosed, i.e. brought to homozygous state at a locus by identity by descent or state, could potentially result in autosomal recessive (AR) rare disease traits. However, the detection and functional interpretation of homozygous duplications from exome sequencing data remains a challenge. We developed a framework algorithm, HMZDupFinder, that is designed to detect exonic homozygous duplications from exome sequencing (ES) data. The HMZDupFinder algorithm can efficiently process large datasets and accurately identifies small intragenic duplications, including those associated with rare disease traits. HMZDupFinder called 965 homozygous duplications with three or less exons from 8,707 ES with a recall rate of 70.9% and a precision of 16.1%. We experimentally confirmed 8/10 rare homozygous duplications. Pathogenicity assessment of these copy number variant alleles allowed clinical genomics contextualization for three homozygous duplications alleles, including two affecting known OMIM disease genes EDAR (MIM# 224900), TNNT1(MIM# 605355), and one variant in a novel candidate disease gene: PAAF1.

LINE-1 retrotransposition and its deregulation in cancers: implications for therapeutic opportunities

Long interspersed element 1 (LINE-1) is the only protein-coding transposon that is active in humans. LINE-1 propagates in the genome using RNA intermediates via retrotransposition. This activity has resulted in LINE-1 sequences occupying approximately one-fifth of our genome. Although most copies of LINE-1 are immobile, ∼100 copies are retrotransposition-competent. Retrotransposition is normally limited via epigenetic silencing, DNA repair, and other host defense mechanisms. In contrast, LINE-1 overexpression and retrotransposition are hallmarks of cancers. Here, we review mechanisms of LINE-1 regulation and how LINE-1 may promote genetic heterogeneity in tumors. Finally, we discuss therapeutic strategies to exploit LINE-1 biology in cancers.

Copy Number Variations in Hereditary Spastic Paraplegia-Related Genes: Evaluation of an Iranian Hereditary Spastic Paraplegia Cohort and Literature Review

Introduction

In human genetic disorders, copy number variations (CNVs) are considered a considerable underlying cause. CNVs are generally detected by array-based methods but can also be discovered by read-depth analysis of whole-exome sequencing (WES) data. We performed WES-based CNV identification in a cohort of 35 Iranian families with hereditary spastic paraplegia (HSP) patients.

Methods

Thirty-five patients whose routine single-nucleotide variants (SNVs) and insertion/deletion analyses from exome data were unrevealing underwent a pipeline of CNV analysis using the read-depth detection method. Subsequently, a comprehensive search about the existence of CNVs in all 84 known HSP-causing genes was carried out in all reported HSP cases, so far.

Results and Discussion

CNV analysis of exome data indicated that 1 patient harbored a heterozygous deletion in exon 17 of the SPAST gene. Multiplex ligation-dependent probe amplification analysis confirmed this deletion in the proband and his affected father. Literature review demonstrated that, to date, pathogenic CNVs have been identified in 30 out of 84 HSP-causing genes (∼36%). However, CNVs in only 17 of these genes were specifically associated with the HSP phenotype. Among them, CNVs were more common in L1CAM, PLP1, SPAST, SPG7, SPG11, and REEP1 genes. The identification of the CNV in 1 of our patients suggests that WES allows the detection of both SNVs and CNVs from a single method without additional costs and execution time. However, because of intrinsic issues of WES in the detection of large rearrangements, it may not yet be exploited to replace the CNV detection methods in standard clinical practice.

Break-induced replication underlies formation of inverted triplications and generates unexpected diversity in haplotype structures

Background

The duplication-triplication/inverted-duplication (DUP-TRP/INV-DUP) structure is a type of complex genomic rearrangement (CGR) hypothesized to result from replicative repair of DNA due to replication fork collapse. It is often mediated by a pair of inverted low-copy repeats (LCR) followed by iterative template switches resulting in at least two breakpoint junctions in cis . Although it has been identified as an important mutation signature of pathogenicity for genomic disorders and cancer genomes, its architecture remains unresolved and is predicted to display at least four structural variation (SV) haplotypes.

Results

Here we studied the genomic architecture of DUP-TRP/INV-DUP by investigating the genomic DNA of 24 patients with neurodevelopmental disorders identified by array comparative genomic hybridization (aCGH) on whom we found evidence for the existence of 4 out of 4 predicted SV haplotypes. Using a combination of short-read genome sequencing (GS), long- read GS, optical genome mapping and StrandSeq the haplotype structure was resolved in 18 samples. This approach refined the point of template switching between inverted LCRs in 4 samples revealing a DNA segment of ∼2.2-5.5 kb of 100% nucleotide similarity. A prediction model was developed to infer the LCR used to mediate the non-allelic homology repair.

Conclusions

These data provide experimental evidence supporting the hypothesis that inverted LCRs act as a recombinant substrate in replication-based repair mechanisms. Such inverted repeats are particularly relevant for formation of copy-number associated inversions, including the DUP-TRP/INV-DUP structures. Moreover, this type of CGR can result in multiple conformers which contributes to generate diverse SV haplotypes in susceptible loci .

Somatic cell structural variant mutagenesis and neurologic disease

Detection of organismal mosaic states for variant alleles faces technical and analytical challenges, as does the association of such variant alleles with susceptibility to neurologic disease. In this issue of Cell Genomics, Maury et al.1 reanalyze genotyping arrays of a schizophrenia cohort providing evidence for the contribution of somatic structural variant mutagenesis and rare variant alleles.

Schizophrenia-associated somatic copy-number variants from 12,834 cases reveal recurrent NRXN1 and ABCB11 disruptions

While germline copy-number variants (CNVs) contribute to schizophrenia (SCZ) risk, the contribution of somatic CNVs (sCNVs)-present in some but not all cells-remains unknown. We identified sCNVs using blood-derived genotype arrays from 12,834 SCZ cases and 11,648 controls, filtering sCNVs at loci recurrently mutated in clonal blood disorders. Likely early-developmental sCNVs were more common in cases (0.91%) than controls (0.51%, p = 2.68e-4), with recurrent somatic deletions of exons 1-5 of the NRXN1 gene in five SCZ cases. Hi-C maps revealed ectopic, allele-specific loops forming between a potential cryptic promoter and non-coding cis-regulatory elements upon 5' deletions in NRXN1. We also observed recurrent intragenic deletions of ABCB11, encoding a transporter implicated in anti-psychotic response, in five treatment-resistant SCZ cases and showed that ABCB11 is specifically enriched in neurons forming mesocortical and mesolimbic dopaminergic projections. Our results indicate potential roles of sCNVs in SCZ risk.

Rare variant enrichment analysis supports GREB1L as a contributory driver gene in the etiology of Mayer-Rokitansky-Küster-Hauser syndrome

Mayer-Rokitansky-Küster-Hauser (MRKH) syndrome is characterized by aplasia of the female reproductive tract; the syndrome can include renal anomalies, absence or dysgenesis, and skeletal anomalies. While functional models have elucidated several candidate genes, only WNT4 (MIM: 603490) variants have been definitively associated with a subtype of MRKH with hyperandrogenism (MIM: 158330). DNA from 148 clinically diagnosed MRKH probands across 144 unrelated families and available family members from North America, Europe, and South America were exome sequenced (ES) and by family-based genomics analyzed for rare likely deleterious variants. A replication cohort consisting of 442 Han Chinese individuals with MRKH was used to further reproduce GREB1L findings in diverse genetic backgrounds. Proband and OMIM phenotypes annotated using the Human Phenotype Ontology were analyzed to quantitatively delineate the phenotypic spectrum associated with GREB1L variant alleles found in our MRKH cohort and those previously published. This study reports 18 novel GREB1L variant alleles, 16 within a multiethnic MRKH cohort and two within a congenital scoliosis cohort. Cohort-wide analyses for a burden of rare variants within a single gene identified likely damaging variants in GREB1L (MIM: 617782), a known disease gene for renal hypoplasia and uterine abnormalities (MIM: 617805), in 16 of 590 MRKH probands. GREB1L variant alleles, including a CNV null allele, were found in 8 MRKH type 1 probands and 8 MRKH type II probands. This study used quantitative phenotypic analyses in a worldwide multiethnic cohort to identify and strengthen the association of GREB1L to isolated uterine agenesis (MRKH type I) and syndromic MRKH type II.

Disruption of FBN1 by an Alu element insertion: A novel genetic cause of Marfan syndrome

Alu elements are retrotransposons with ubiquitous presence in the human genome that have contributed to human genomic diversity and health. These approximately 300-bp sequences can cause or mediate disease by disrupting coding/splicing regions in the germline, by insertional mutagenesis in somatic cells, and in promoting formation of copy-number variants. Alu elements may also disrupt epigenetic regulation by affecting non-coding regulatory regions. There are increasing reports of apparently sporadic and inherited genetic disorders caused by Alu-related gene disruption, but Marfan syndrome resulting from Alu element insertion has not been previously described. We report a family with classic features of Marfan syndrome whose previous FBN1 genetic testing was inconclusive. Using contemporary next-generation sequencing and bioinformatics analysis, a pathogenic/disruptive Alu insertion occurring in the coding region of the FBN1 gene was identified (c.6564_6565insAlu; p. Glu2189fs) and was confirmed and specified further with Sanger sequencing. This identified the molecular basis of disease in the family that was missed using previous genetic testing technologies and highlights a novel pathogenic mechanism for Marfan syndrome. This case adds to the growing literature of Mendelian diseases caused by Alu retrotransposition, and it also shows the growing capability of genomic technologies for detecting atypical mutation events.

Systematic assessment of the contribution of structural variants to inherited retinal diseases

Despite increasing success in determining genetic diagnosis for patients with inherited retinal diseases (IRDs), mutations in about 30% of the IRD cases remain unclear or unsettled after targeted gene panel or whole exome sequencing. In this study, we aimed to investigate the contributions of structural variants (SVs) to settling the molecular diagnosis of IRD with whole-genome sequencing (WGS). A cohort of 755 IRD patients whose pathogenic mutations remain undefined were subjected to WGS. Four SV calling algorithms including include MANTA, DELLY, LUMPY and CNVnator were used to detect SVs throughout the genome. All SVs identified by any one of these four algorithms were included for further analysis. AnnotSV was used to annotate these SVs. SVs that overlap with known IRD-associated genes were examined with sequencing coverage, junction reads and discordant read pairs. Polymerase Chain Reaction (PCR) followed by Sanger sequencing was used to further confirm the SVs and identify the breakpoints. Segregation of the candidate pathogenic alleles with the disease was performed when possible. A total of 16 candidate pathogenic SVs were identified in 16 families, including deletions and inversions, representing 2.1% of patients with previously unsolved IRDs. Autosomal dominant, autosomal recessive and X-linked inheritance of disease-causing SVs were observed in 12 different genes. Among these, SVs in CLN3, EYS and PRPF31 were found in multiple families. Our study suggests that the contribution of SVs detected by short-read WGS is about 0.25% of our IRD patient cohort and is significantly lower than that of single nucleotide changes and small insertions and deletions.

Adversary of DNA integrity: A long non-coding RNA stimulates driver oncogenic chromosomal rearrangement in human thyroid cells

The flurry of publications devoted to the functions of long non-coding RNAs (lncRNAs) published in the last decade leaves no doubt about the exceptional importance of lncRNAs in various areas including tumor biology. However, contribution of lncRNAs to the early stages of oncogenesis remains poorly understood. In this study we explored a new role for lncRNAs: stimulation of specific chromosomal rearrangements upon DNA damage. We demonstrated that lncRNA CASTL1 (ENSG00000269945) stimulates the formation of the CCDC6-RET inversion (RET/PTC1) in human thyroid cells subjected to radiation or chemical DNA damage. Facilitation of chromosomal rearrangement requires lncRNA to contain regions complementary to the introns of both CCDC6 and RET genes as deletion of these regions deprives CASTL1 of the ability to stimulate the gene fusion. We found that CASTL1 expression is elevated in tumors with CCDC6-RET fusion which is the most frequent rearrangement in papillary thyroid carcinoma. Our results open a new venue for the studies of early oncogenesis in various tumor types, especially those associated with physical or chemical DNA damage.

Shuffling the yeast genome using CRISPR/Cas9-generated DSBs that target the transposable Ty1 elements

Although homologous recombination between transposable elements can drive genomic evolution in yeast by facilitating chromosomal rearrangements, the details of the underlying mechanisms are not fully clarified. In the genome of the yeast Saccharomyces cerevisiae, the most common class of transposon is the retrotransposon Ty1. Here, we explored how Cas9-induced double-strand breaks (DSBs) directed to Ty1 elements produce genomic alterations in this yeast species. Following Cas9 induction, we observed a significant elevation of chromosome rearrangements such as deletions, duplications and translocations. In addition, we found elevated rates of mitotic recombination, resulting in loss of heterozygosity. Using Southern analysis coupled with short- and long-read DNA sequencing, we revealed important features of recombination induced in retrotransposons. Almost all of the chromosomal rearrangements reflect the repair of DSBs at Ty1 elements by non-allelic homologous recombination; clustered Ty elements were hotspots for chromosome rearrangements. In contrast, a large proportion (about three-fourths) of the allelic mitotic recombination events have breakpoints in unique sequences. Our analysis suggests that some of the latter events reflect extensive processing of the broken ends produced in the Ty element that extend into unique sequences resulting in break-induced replication. Finally, we found that haploid and diploid strain have different preferences for the pathways used to repair double-stranded DNA breaks. Our findings demonstrate the importance of DNA lesions in retrotransposons in driving genome evolution.

Systematic assessment of the contribution of structural variants to inherited retinal diseases

Despite increasing success in determining genetic diagnosis for patients with inherited retinal diseases (IRDs), mutations in about 30% of the IRD cases remain unclear or unsettled after targeted gene panel or whole exome sequencing. In this study, we aimed to investigate the contributions of structural variants (SVs) to settling the molecular diagnosis of IRD with whole-genome sequencing (WGS). A cohort of 755 IRD patients whose pathogenic mutations remain undefined was subjected to WGS. Four SV calling algorithms including include MANTA, DELLY, LUMPY, and CNVnator were used to detect SVs throughout the genome. All SVs identified by any one of these four algorithms were included for further analysis. AnnotSV was used to annotate these SVs. SVs that overlap with known IRD-associated genes were examined with sequencing coverage, junction reads, and discordant read pairs. PCR followed by Sanger sequencing was used to further confirm the SVs and identify the breakpoints. Segregation of the candidate pathogenic alleles with the disease was performed when possible. In total, sixteen candidate pathogenic SVs were identified in sixteen families, including deletions and inversions, representing 2.1% of patients with previously unsolved IRDs. Autosomal dominant, autosomal recessive, and X-linked inheritance of disease-causing SVs were observed in 12 different genes. Among these, SVs in CLN3, EYS, PRPF31 were found in multiple families. Our study suggests that the contribution of SVs detected by short-read WGS is about 0.25% of our IRD patient cohort and is significantly lower than that of single nucleotide changes and small insertions and deletions.

Comprehensive analysis of the PRPF31 gene in retinitis pigmentosa patients: Four novel Alu-mediated copy number variations at the PRPF31 locus

Retinitis pigmentosa (RP) is a monogenic disease characterized by irreversible degeneration of the retina. PRPF31, the second most common causative gene of autosomal dominant RP, frequently harbors copy number variations (CNVs), but the underlying mechanism is unclear. In this study, we summarized the phenotypic and genotypic characteristics of 18 RP families (F01-F18) with variants in PRPF31. The prevalence of PRPF31 variants in our cohort of Chinese RP families was 1.7% (18/1024). Seventeen different variants in PRPF31 were detected, including eight novel variants. Notably, four novel CNVs encompassing PRPF31, with a proportion of 22.2% (4/18), were validated to harbor gross deletions involving Alu/Alu-mediated rearrangements (AAMRs) in the same orientation. Among a total of 12 CNVs of PRPF31 with breakpoints mapped on nucleotide-resolution, 10 variants (83.3%) were presumably mediated by Alu elements. Furthermore, we described the correlation between the genotypes and phenotypes in PRPF31-related RP. Our findings expand the mutational spectrum of the PRPF31 gene and provide strong evidence that Alu elements of PRPF31 probably contribute to the susceptibility to genomic rearrangement in this locus.

Transposable element-mediated rearrangements are prevalent in human genomes

Transposable elements constitute about half of human genomes, and their role in generating human variation through retrotransposition is broadly studied and appreciated. Structural variants mediated by transposons, which we call transposable element-mediated rearrangements (TEMRs), are less well studied, and the mechanisms leading to their formation as well as their broader impact on human diversity are poorly understood. Here, we identify 493 unique TEMRs across the genomes of three individuals. While homology directed repair is the dominant driver of TEMRs, our sequence-resolved TEMR resource allows us to identify complex inversion breakpoints, triplications or other high copy number polymorphisms, and additional complexities. TEMRs are enriched in genic loci and can create potentially important risk alleles such as a deletion in TRIM65, a known cancer biomarker and therapeutic target. These findings expand our understanding of this important class of structural variation, the mechanisms responsible for their formation, and establish them as an important driver of human diversity.

Genetic genealogy uncovers a founder deletion mutation in the cerebral cavernous malformations 2 gene

Cerebral cavernous malformations (CCM) are vascular malformations consisting of collections of enlarged capillaries occurring in the brain or spinal cord. These vascular malformations can occur sporadically or susceptibility to develop these can be inherited as an autosomal dominant trait due to mutation in one of three genes. Over a decade ago, we described a 77.6 Kb germline deletion spanning exons 2-10 in the CCM2 gene found in multiple affected individuals from seemingly unrelated families. Segregation analysis using linked, microsatellite markers indicated that this deletion may have arisen at least twice independently. In the ensuing decades, many more CCM patients have been identified with this deletion. In this present study we examined 27 reportedly unrelated affected individuals with this deletion. To investigate the origin of the deletion at base pair level resolution, we sequenced approximately 10 Kb upstream and downstream from the recombination junction on the deleted allele. All patients showed the identical SNP haplotype across this combined 20 Kb interval. In parallel, genealogical records have traced 11 of these individuals to five separate pedigrees dating as far back as the 1600-1700s. These haplotype and genealogical data suggest that these families and the remaining "unrelated" samples converge on a common ancestor due to a founder mutation occurring centuries ago on the North American continent. We also note that another gene, NACAD, is included in this deletion. Although patient self-reporting does not indicate an apparent phenotypic consequence for heterozygous deletion of NACAD, further investigation is warranted for these patients.

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.

A reverse genetics and genomics approach to gene paralog function and disease: Myokymia and the juxtaparanode

The leucine-rich glioma-inactivated (LGI) family consists of four highly conserved paralogous genes, LGI1-4, that are highly expressed in mammalian central and/or peripheral nervous systems. LGI1 antibodies are detected in subjects with autoimmune limbic encephalitis and peripheral nerve hyperexcitability syndromes (PNHSs) such as Isaacs and Morvan syndromes. Pathogenic variations of LGI1 and LGI4 are associated with neurological disorders as disease traits including familial temporal lobe epilepsy and neurogenic arthrogryposis multiplex congenita 1 with myelin defects, respectively. No human disease has been reported associated with either LGI2 or LGI3. We implemented exome sequencing and family-based genomics to identify individuals with deleterious variants in LGI3 and utilized GeneMatcher to connect practitioners and researchers worldwide to investigate the clinical and electrophysiological phenotype in affected subjects. We also generated Lgi3-null mice and performed peripheral nerve dissection and immunohistochemistry to examine the juxtaparanode LGI3 microarchitecture. As a result, we identified 16 individuals from eight unrelated families with loss-of-function (LoF) bi-allelic variants in LGI3. Deep phenotypic characterization showed LGI3 LoF causes a potentially clinically recognizable PNHS trait characterized by global developmental delay, intellectual disability, distal deformities with diminished reflexes, visible facial myokymia, and distinctive electromyographic features suggestive of motor nerve instability. Lgi3-null mice showed reduced and mis-localized Kv1 channel complexes in myelinated peripheral axons. Our data demonstrate bi-allelic LoF variants in LGI3 cause a clinically distinguishable disease trait of PNHS, most likely caused by disturbed Kv1 channel distribution in the absence of LGI3.

Pathogenic variants detected by RNA sequencing in Cornelia de Lange syndrome

Recent studies suggest that transcript isoforms significantly overlap (approximately 60%) between brain tissue and Epstein-Barr virus-transformed lymphoblastoid cell lines (LCLs). Interestingly, 14 cohesion-related genes with variants that cause Cornelia de Lange Syndrome (CdLS) are highly expressed in the brain and LCLs. In this context, we first performed RNA sequencing of LCLs from 22 solved (with pathogenic variants) and 19 unsolved (with no confirmed variants) CdLS cases. Next, an RNA sequencing pipeline was developed using solved cases with two different methods: short variant analysis (for single-nucleotide and indel variants) and aberrant splicing detection analysis. Then, 19 unsolved cases were subsequently applied to our pipeline, and four pathogenic variants in NIPBL (one inframe deletion and three intronic variants) were newly identified. Two of three intronic variants were located at Alu elements in deep-intronic regions, creating cryptic exons. RNA sequencing with LCLs was useful for identifying hidden variants in exome-negative cases.

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.

Phenotypic and mutational spectrum of ROR2-related Robinow syndrome

Robinow syndrome is characterized by a triad of craniofacial dysmorphisms, disproportionate-limb short stature and genital hypoplasia. A significant degree of phenotypic variability seems to correlate with different genes/loci. Disturbances of the non-canonical WNT-pathway have been identified as the main cause of the syndrome. Biallelic variants in ROR2 cause an autosomal recessive form of the syndrome with distinctive skeletal findings. Twenty-two patients with a clinical diagnosis of autosomal recessive Robinow syndrome were screened for variants in ROR2 using multiple molecular approaches. We identified 25 putatively pathogenic ROR2 variants, 16 novel, including single nucleotide variants and exonic deletions. Detailed phenotypic analyses revealed that all subjects presented with a prominent forehead, hypertelorism, short nose, abnormality of the nasal tip, brachydactyly, mesomelic limb shortening, short stature and genital hypoplasia in male patients. A total of 19 clinical features were present in more than 75% of the subjects, thus pointing to an overall uniformity of the phenotype. Disease-causing variants in ROR2, contribute to a clinically recognizable AR trait phenotype with multiple skeletal defects. A comprehensive quantitative clinical evaluation this cohort delineated the phenotypic spectrum of ROR2-related Robinow syndrome. The identification of exonic deletion variant alleles further supports the contention of a loss-of-function mechanism in the etiology of the syndrome. This article is protected by copyright. All rights reserved.

Novel pathogenic variants and quantitative phenotypic analyses of Robinow syndrome: WNT signaling perturbation and phenotypic variability

Robinow syndrome (RS) is a genetically heterogeneous disorder with six genes that converge on the WNT/planar cell polarity (PCP) signaling pathway implicated (DVL1, DVL3, FZD2, NXN, ROR2, and WNT5A). RS is characterized by skeletal dysplasia and distinctive facial and physical characteristics. To further explore the genetic heterogeneity, paralog contribution, and phenotypic variability of RS, we investigated a cohort of 22 individuals clinically diagnosed with RS from 18 unrelated families. Pathogenic or likely pathogenic variants in genes associated with RS or RS phenocopies were identified in all 22 individuals, including the first variant to be reported in DVL2. We retrospectively collected medical records of 16 individuals from this cohort and extracted clinical descriptions from 52 previously published cases. We performed Human Phenotype Ontology (HPO) based quantitative phenotypic analyses to dissect allele-specific phenotypic differences. Individuals with FZD2 variants clustered into two groups with demonstrable phenotypic differences between those with missense and truncating alleles. Probands with biallelic NXN variants clustered together with the majority of probands carrying DVL1, DVL2, and DVL3 variants, demonstrating no phenotypic distinction between the NXN-autosomal recessive and dominant forms of RS. While phenotypically similar diseases on the RS differential matched through HPO analysis, clustering using phenotype similarity score placed RS-associated phenotypes in a unique cluster containing WNT5A, FZD2, and ROR2 apart from non-RS-associated paralogs. Through human phenotype analyses of this RS cohort and OMIM clinical synopses of Mendelian disease, this study begins to tease apart specific biologic roles for non-canonical WNT-pathway proteins.

CRISPR/Cas9-induced gene conversion between ATAD3 paralogs

Paralogs and pseudogenes are abundant within the human genome, and can mediate non-allelic homologous recombination (NAHR) or gene conversion events. The ATAD3 locus contains three paralogs situated in tandem, and is therefore prone to NAHR-mediated deletions and duplications associated with severe neurological phenotypes. To study this locus further, we aimed to generate biallelic loss-of-function variants in ATAD3A by CRISPR/Cas9 genome editing. Unexpectedly, two of the generated clones underwent gene conversion, as evidenced by replacement of the targeted sequence of ATAD3A by a donor sequence from its paralog ATAD3B. We highlight the complexity of CRISPR/Cas9 design, end-product formation, and recombination repair mechanisms for CRISPR/Cas9 delivery as a nucleic acid molecular therapy when targeting genes that have paralogs or pseudogenes, and advocate meticulous evaluation of resultant clones in model organisms. In addition, we suggest that endogenous gene conversion may be used to repair missense variants in genes with paralogs or pseudogenes.

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.

High prevalence of multilocus pathogenic variation in neurodevelopmental disorders in the Turkish population

Neurodevelopmental disorders (NDDs) are clinically and genetically heterogenous; many such disorders are secondary to perturbation in brain development and/or function. The prevalence of NDDs is > 3%, resulting in significant sociocultural and economic challenges to society. With recent advances in family-based genomics, rare-variant analyses, and further exploration of the Clan Genomics hypothesis, there has been a logarithmic explosion in neurogenetic "disease-associated genes" molecular etiology and biology of NDDs; however, the majority of NDDs remain molecularly undiagnosed. We applied genome-wide screening technologies, including exome sequencing (ES) and whole-genome sequencing (WGS), to identify the molecular etiology of 234 newly enrolled subjects and 20 previously unsolved Turkish NDD families. In 176 of the 234 studied families (75.2%), a plausible and genetically parsimonious molecular etiology was identified. Out of 176 solved families, deleterious variants were identified in 218 distinct genes, further documenting the enormous genetic heterogeneity and diverse perturbations in human biology underlying NDDs. We propose 86 candidate disease-trait-associated genes for an NDD phenotype. Importantly, on the basis of objective and internally established variant prioritization criteria, we identified 51 families (51/176 = 28.9%) with multilocus pathogenic variation (MPV), mostly driven by runs of homozygosity (ROHs) - reflecting genomic segments/haplotypes that are identical-by-descent. Furthermore, with the use of additional bioinformatic tools and expansion of ES to additional family members, we established a molecular diagnosis in 5 out of 20 families (25%) who remained undiagnosed in our previously studied NDD cohort emanating from Turkey.

Bi-allelic loss-of-function variants in BCAS3 cause a syndromic neurodevelopmental disorder

BCAS3 microtubule-associated cell migration factor (BCAS3) is a large, highly conserved cytoskeletal protein previously proposed to be critical in angiogenesis and implicated in human embryogenesis and tumorigenesis. Here, we established BCAS3 loss-of-function variants as causative for a neurodevelopmental disorder. We report 15 individuals from eight unrelated families with germline bi-allelic loss-of-function variants in BCAS3. All probands share a global developmental delay accompanied by pyramidal tract involvement, microcephaly, short stature, strabismus, dysmorphic facial features, and seizures. The human phenotype is less severe compared with the Bcas3 knockout mouse model and cannot be explained by angiogenic defects alone. Consistent with being loss-of-function alleles, we observed absence of BCAS3 in probands' primary fibroblasts. By comparing the transcriptomic and proteomic data based on probands' fibroblasts with those of the knockout mouse model, we identified similar dysregulated pathways resulting from over-representation analysis, while the dysregulation of some proposed key interactors could not be confirmed. Together with the results from a tissue-specific Drosophila loss-of-function model, we demonstrate a vital role for BCAS3 in neural tissue development.

A novel homozygous SLC13A5 whole-gene deletion generated by Alu/Alu-mediated rearrangement in an Iraqi family with epileptic encephalopathy

Biallelic loss-of-function (LoF) of SLC13A5 (solute carrier family 13, member 5) induced deficiency in sodium/citrate transporter (NaCT) causes autosomal recessive developmental epileptic encephalopathy 25 with hypoplastic amelogenesis imperfecta (DEE25; MIM #615905). Many pathogenic SLC13A5 single nucleotide variants (SNVs) and small indels have been described; however, no cases with copy number variants (CNVs) have been sufficiently investigated. We describe a consanguineous Iraqi family harboring an 88.5 kb homozygous deletion including SLC13A5 in Chr17p13.1. The three affected male siblings exhibit neonatal-onset epilepsy with fever-sensitivity, recurrent status epilepticus, global developmental delay/intellectual disability (GDD/ID), and other variable neurological findings as shared phenotypical features of DEE25. Two of the three affected subjects exhibit hypoplastic amelogenesis imperfecta (AI), while the proband shows no evidence of dental abnormalities or AI at 2 years of age with apparently unaffected primary dentition. Characterization of the genomic architecture at this locus revealed evidence for genomic instability generated by an Alu/Alu-mediated rearrangement; confirmed by break-point junction Sanger sequencing. This multiplex family from a distinct population elucidates the phenotypic consequence of complete LoF of SLC13A5 and illustrates the importance of read-depth-based CNV detection in comprehensive exome sequencing analysis to solve cases that otherwise remain molecularly unsolved.

Genetic-variant hotspots and hotspot clusters in the human genome facilitating adaptation while increasing instability

Background

Genetic variants, underlining phenotypic diversity, are known to distribute unevenly in the human genome. A comprehensive understanding of the distributions of different genetic variants is important for insights into genetic functions and disorders.

Methods

Herein, a sliding-window scan of regional densities of eight kinds of germline genetic variants, including single-nucleotide-polymorphisms (SNPs) and four size-classes of copy-number-variations (CNVs) in the human genome has been performed.

Results

The study has identified 44,379 hotspots with high genetic-variant densities, and 1135 hotspot clusters comprising more than one type of hotspots, accounting for 3.1% and 0.2% of the genome respectively. The hotspots and clusters are found to co-localize with different functional genomic features, as exemplified by the associations of hotspots of middle-size CNVs with histone-modification sites, work with balancing and positive selections to meet the need for diversity in immune proteins, and facilitate the development of sensory-perception and neuroactive ligand-receptor interaction pathways in the function-sparse late-replicating genomic sequences. Genetic variants of different lengths co-localize with retrotransposons of different ages on a “long-with-young” and “short-with-all” basis. Hotspots and clusters are highly associated with tumor suppressor genes and oncogenes (p < 10⁻¹⁰), and enriched with somatic tumor CNVs and the trait- and disease-associated SNPs identified by genome-wise association studies, exceeding tenfold enrichment in clusters comprising SNPs and extra-long CNVs.

Conclusions

In conclusion, the genetic-variant hotspots and clusters represent two-edged swords that spearhead both positive and negative genomic changes. Their strong associations with complex traits and diseases also open up a potential “Common Disease-Hotspot Variant” approach to the missing heritability problem.

Supplementary Information

The online version contains supplementary material available at 10.1186/s40246-021-00318-3.

Altered DNA repair creates novel Alu/Alu repeat-mediated deletions

Alu elements are the most abundant source of nonallelic homology that influences genetic instability in the human genome. When there is a DNA double-stranded break, the Alu element's high copy number, moderate length and distance and mismatch between elements uniquely influence recombination processes. We utilize a reporter-gene assay to show the complex influence of Alu mismatches on Alu-related repeat-mediated deletions (RMDs). The Alu/Alu heteroduplex intermediate can result in a nonallelic homologous recombination (HR). Alternatively, the heteroduplex can result in various DNA breaks around the Alu elements caused by competing nucleases. These breaks can undergo Alt-nonhomologous end joining to cause deletions focused around the Alu elements. Formation of these heteroduplex intermediates is largely RAD52 dependent. Cells with low ERCC1 levels utilize more of these alternatives resolutions, while cells with MSH2 defects tend to have more RMDs with a specific increase in the HR events. Therefore, Alu elements are expected to create different forms of deletions in various cancers depending on a number of these DNA repair defects.

BLM has Contrary Effects on Repeat-Mediated Deletions, based on the Distance of DNA DSBs to a Repeat and Repeat Divergence

Repeat-mediated deletions (RMDs) often involve repetitive elements (e.g., short interspersed elements) with sequence divergence that is separated by several kilobase pairs (kbps). We have examined RMDs induced by DNA double-strand breaks (DSBs) under varying conditions of repeat sequence divergence (identical versus 1% and 3% divergent) and DSB/repeat distance (16 bp–28.4 kbp). We find that the BLM helicase promotes RMDs with long DSB/repeat distances (e.g., 28.4 kbp), which is consistent with a role in extensive DSB end resection, because the resection nucleases EXO1 and DNA2 affect RMDs similarly to BLM. In contrast, BLM suppresses RMDs with sequence divergence and intermediate (e.g., 3.3 kbp) DSB/repeat distances, which supports a role in heteroduplex rejection. The role of BLM in heteroduplex rejection is not epistatic with MSH2 and is independent of the annealing factor RAD52. Accordingly, the role of BLM on RMDs is substantially affected by DSB/repeat distance and repeat sequence divergence.

Mendez-Dorantes et al. identify the BLM helicase as a key regulator of repeat-mediated deletions (RMDs). BLM, EXO1, and DNA2 mediate RMDs with remarkably long DNA break/repeat distances. BLM suppresses RMDs with sequence divergence that is optimal with a long non-homologous tail and is independent of MSH2 and RAD52.

Palindromes in DNA-A Risk for Genome Stability and Implications in Cancer

A palindrome in DNA consists of two closely spaced or adjacent inverted repeats. Certain palindromes have important biological functions as parts of various cis-acting elements and protein binding sites. However, many palindromes are known as fragile sites in the genome, sites prone to chromosome breakage which can lead to various genetic rearrangements or even cell death. The ability of certain palindromes to initiate genetic recombination lies in their ability to form secondary structures in DNA which can cause replication stalling and double-strand breaks. Given their recombinogenic nature, it is not surprising that palindromes in the human genome are involved in genetic rearrangements in cancer cells as well as other known recurrent translocations and deletions associated with certain syndromes in humans. Here, we bring an overview of current understanding and knowledge on molecular mechanisms of palindrome recombinogenicity and discuss possible implications of DNA palindromes in carcinogenesis. Furthermore, we overview the data on known palindromic sequences in the human genome and efforts to estimate their number and distribution, as well as underlying mechanisms of genetic rearrangements specific palindromic sequences cause.

Intronic Breakpoint Signatures Enhance Detection and Characterization of Clinically Relevant Germline Structural Variants

The relevance of large copy number variants (CNVs) to hereditary disorders has been long recognized, and population sequencing efforts have chronicled many common structural variants (SVs). However, limited data are available on the clinical contribution of rare germline SVs. Here, a detailed characterization of SVs identified using targeted next-generation sequencing was performed. Across 50 genes associated with hereditary cancer and cardiovascular disorders, a minimum of 828 unique SVs were reported, including 584 fully characterized SVs. Almost 40% of CNVs were <5 kb, with one in three deletions impacting a single exon. Additionally, 36 mid-range deletions/duplications (50 to 250 bp), 21 mobile element insertions, 6 inversions, and 27 complex rearrangements were detected. This data set was used to model SV detection in a bioinformatics pipeline solely relying on read depth, which revealed that genome sequencing (30×) allows detection of 71%, a 500× panel only targeting coding regions 53%, and exome sequencing (100×) <20% of characterized SVs. SVs accounted for 14.1% of all unique pathogenic variants, supporting the importance of SVs in hereditary disorders. Robust SV detection requires an ensemble of variant-calling algorithms that utilize sequencing of intronic regions. These algorithms should use distinct data features representative of each class of mutational mechanism, including recombination between two sequences sharing high similarity, covariants inserted between CNV breakpoints, and complex rearrangements containing inverted sequences.

Identification of the first Alu-mediated gross deletion involving the BCKDHA gene in a compound heterozygous patient with maple syrup urine disease

AIMS

To investigate a family with clinical symptoms of maple syrup urine disease and reveal a genetic cause underlying this disease.

METHODS

Targeted capture sequencing was used to screen for mutations in the patient. Real-Time PCR was carried out to perform exon 1, 5, 9 CNV analysis of samples from the patient's father, mother and sister. Whole genome sequencing was performed to map the approximate location of the break points of the gross deletion. Long-range PCR and Sanger sequencing were performed to identify the length of the deletion and to locate the break points.

RESULTS

The patient is a compound heterozygous mutation including a small deletion mutation (c.1227_1229del chr19: 41930402) and a gross novel deletion including exon1-9 in BCKDHA. The junction site of the gross deletion was localized within a microhomologous sequence in two Alu elements.

CONCLUSIONS

This study is the first time report on rearrangement sequences in BCKDHA mediated by Alu element, which resulted in MSUD. Our results may also offer new insights into the formation and pathogenicity of MSUD, and may be useful to genetic counseling and genetic testing.

Chromoanagenesis, the mechanisms of a genomic chaos

Designated under the name of chromoanagenesis, the phenomena of chromothripsis, chromanasynthesis and chromoplexy constitute new types of complex rearrangements, including many genomic alterations localized on a few chromosomal regions, and whose discovery over the last decade has changed our perception about the formation of chromosomal abnormalities and their etiology. Although exhibiting specific features, these new catastrophic mechanisms generally occur within a single cell cycle and their emergence is closely linked to genomic instability. Various non-exclusive exogenous or cellular mechanisms capable of generating chromoanagenesis have been evoked. However, recent experimental data shed light on 2 major processes, which following a defect in the mitotic segregation of chromosomes, can generate a cascade of cellular events leading to chromoanagenesis. These mechanisms are the formation of micronuclei integrating isolated chromosomal material, and the occurrence of chromatin bridges around chromosomal material resulting from telomeric fusions. In both cases, the cellular and molecular mechanisms of fragmentation, repair and transmission of damaged chromosomal material are consistent with the features of chromoanagenesis-related complex chromosomal rearrangements. In this review, we introduce each type of chromoanagenesis, and describe the experimental models that have allowed to validate the existence of chromoanagenesis events and to better understand their cellular mechanisms of formation and transmission, as well as their impact on the stability and the plasticity of the genome.

Hotspots of Human Mutation

Mutation of the human genome results in three classes of genomic variation: single nucleotide variants; short insertions or deletions; and large structural variants (SVs). Some mutations occur during normal processes, such as meiotic recombination or B cell development, and others result from DNA replication or aberrant repair of breaks in sequence-specific contexts. Regardless of mechanism, mutations are subject to selection, and some hotspots can manifest in disease. Here, we discuss genomic regions prone to mutation, mechanisms contributing to mutation susceptibility, and the processes leading to their accumulation in normal and somatic genomes. With further, more accurate human genome sequencing, additional mutation hotspots, mechanistic details of their formation, and the relevance of hotspots to evolution and disease are likely to be discovered.

Clinical genomics and contextualizing genome variation in the diagnostic laboratory

INTRODUCTION

The human genome contains the instructions for the development and biological homeostasis of the human organism and the genetic transmission of traits. Genome variation in human populations is the basis of evolution; individual or personal genomes vary tremendously, making each of us truly unique.

AREAS COVERED

Assaying this individual variation using genomic technologies has many applications in clinical medicine, from elucidating the biology of disease to designing strategies to ameliorate perturbations from homeostasis. Detecting pathogenic rare variation in a genome may provide a molecular diagnosis that can be informative for patient management and family healthcare.

EXPERT OPINION

Despite the increasing clinical use of unbiased genomic testing, including chromosome microarray analysis (CMA) with array comparative genomic hybridization (aCGH) or SNP arrays, clinical exome sequencing (cES), and whole-genome sequencing (WGS), to survey genome-wide for molecular aberrations, clinical acumen paired with an understanding of the limitations of each testing type will be needed to achieve molecular diagnoses. Potential opportunities for improving case solved rates, functionally annotating the majority of genes in the human genome, and further understanding genetic contributions to disease will empower clinical genomics and the precision medicine initiative.

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.

Evaluating the probability of CRISPR-based gene drive contaminating another species

The probability D that a given clustered regularly interspaced short palindromic repeats (CRISPR)-based gene drive element contaminates another, nontarget species can be estimated by the following Drive Risk Assessment Quantitative Estimate (DRAQUE) Equation: D = h y b + t r a n s f × e x p r e s s × c u t × f l a n k × i m m u n e × n o n e x t i n c t with hyb = probability of hybridization between the target species and a nontarget species; transf = probability of horizontal transfer of a piece of DNA containing the gene drive cassette from the target species to a nontarget species (with no hybridization); express = probability that the Cas9 and guide RNA genes are expressed; cut = probability that the CRISPR-guide RNA recognizes and cuts at a DNA site in the new host; flank = probability that the gene drive cassette inserts at the cut site; immune = probability that the immune system does not reject Cas9-expressing cells; nonextinct = probability of invasion of the drive within the population. We discuss and estimate each of the seven parameters of the equation, with particular emphasis on possible transfers within insects, and between rodents and humans. We conclude from current data that the probability of a gene drive cassette to contaminate another species is not insignificant. We propose strategies to reduce this risk and call for more work on estimating all the parameters of the formula.

Deficiency of Adenosine Deaminase 2 (DADA2): Hidden Variants, Reduced Penetrance, and Unusual Inheritance

PURPOSE

Deficiency of adenosine deaminase 2 (DADA2) is an autosomal recessive disorder that manifests with fever, early-onset vasculitis, strokes, and hematologic dysfunction. This study aimed to identify disease-causing variants by conventional Sanger and whole exome sequencing in two families suspected to have DADA2 and non-confirmatory genotypes. ADA2 enzymatic assay confirmed the clinical diagnosis of DADA2. Molecular diagnosis was important to accurately identify other family members at risk.

METHODS

We used a variety of sequencing technologies, ADA2 enzymatic testing, and molecular methods including qRT-PCR and MLPA.

RESULTS

Exome sequencing identified heterozygosity for the known pathogenic variant ADA2: c.1358A>G, p.Tyr453Cys in a 14-year-old female with a history of ischemic strokes, livedo, and vasculitis. No second pathogenic variant could be identified. ADA2 enzymatic testing in combination with quantitative RT-PCR suggested a loss-of-function allele. Subsequent genome sequencing identified a canonical splice site variant, c.-47+2T>C, within the 5'UTR of ADA2. Two of her unaffected siblings were found to carry the same two pathogenic variants. A homozygous 800-bp duplication comprising exon 7 of ADA2 was identified in a 5-year-old female with features consistent with Diamond-Blackfan anemia (DBA). The duplication was missed by Sanger sequencing of ADA2, chromosomal microarray, and exome sequencing but was detected by MLPA in combination with long-read PCR sequencing. The exon 7 duplication was also identified in her non-symptomatic father and younger sister.

CONCLUSIONS

ADA2 pathogenic variants may not be detected by conventional sequencing and genetic testing and may require the incorporation of additional diagnostic methods. A definitive molecular diagnosis is crucial for all family members to make informed treatment decisions.

A single mutation in the ACTR8 gene associated with lineage-specific expression in primates

Background

Alternative splicing (AS) generates various transcripts from a single gene and thus plays a significant role in transcriptomic diversity and proteomic complexity. Alu elements are primate-specific transposable elements (TEs) and can provide a donor or acceptor site for AS. In a study on TE-mediated AS, we recently identified a novel AluSz6-exonized ACTR8 transcript of the crab-eating monkey (Macaca fascicularis). In the present study, we sought to determine the molecular mechanism of AluSz6 exonization of the ACTR8 gene and investigate its evolutionary and functional consequences in the crab-eating monkey.

Results

We performed RT-PCR and genomic PCR to analyze AluSz6 exonization in the ACTR8 gene and the expression of the AluSz6-exonized transcript in nine primate samples, including prosimians, New world monkeys, Old world monkeys, and hominoids. AluSz6 integration was estimated to have occurred before the divergence of simians and prosimians. The Alu-exonized transcript obtained by AS was lineage-specific and expressed only in Old world monkeys and apes, and humans. This lineage-specific expression was caused by a single G duplication in AluSz6, which provides a new canonical 5′ splicing site. We further identified other alternative transcripts that were unaffected by the AluSz6 insertion. Finally, we observed that the alternative transcripts were transcribed into new isoforms with C-terminus deletion, and in silico analysis showed that these isoforms do not have a destructive function.

Conclusions

The single G duplication in the TE sequence is the source of TE exonization and AS, and this mutation may suffer a different fate of ACTR8 gene expression during primate evolution.

ALU non-B-DNA conformations, flipons, binary codes and evolution

ALUs contribute to genetic diversity by altering DNA's linear sequence through retrotransposition, recombination and repair. ALUs also have the potential to form alternative non-B-DNA conformations such as Z-DNA, triplexes and quadruplexes that alter the read-out of information from the genome. I suggest here these structures enable the rapid reprogramming of cellular pathways to offset DNA damage and regulate inflammation. The experimental data supporting this form of genetic encoding is presented. ALU sequence motifs that form non-B-DNA conformations under physiological conditions are called flipons. Flipons are binary switches. They are dissipative structures that trade energy for information. By efficiently targeting cellular machines to active genes, flipons expand the repertoire of RNAs compiled from a gene. Their action greatly increases the informational capacity of linearly encoded genomes. Flipons are programmable by epigenetic modification, synchronizing cellular events by altering both chromatin state and nucleosome phasing. Different classes of flipon exist. Z-flipons are based on Z-DNA and modify the transcripts compiled from a gene. T-flipons are based on triplexes and localize non-coding RNAs that direct the assembly of cellular machines. G-flipons are based on G-quadruplexes and sense DNA damage, then trigger the appropriate protective responses. Flipon conformation is dynamic, changing with context. When frozen in one state, flipons often cause disease. The propagation of flipons throughout the genome by ALU elements represents a novel evolutionary innovation that allows for rapid change. Each ALU insertion creates variability by extracting a different set of information from the neighbourhood in which it lands. By elaborating on already successful adaptations, the newly compiled transcripts work with the old to enhance survival. Systems that optimize flipon settings through learning can adapt faster than with other forms of evolution. They avoid the risk of relying on random and irreversible codon rewrites.