Genome architecture, rearrangements and genomic disorders

Genome-wide maps of highly-similar intrachromosomal repeats that can mediate ectopic recombination in three human genome assemblies

Repeated sequences spread throughout the genome play important roles in shaping the structure of chromosomes and facilitating the generation of new genomic variation through structural rearrangements. Several mechanisms of structural variation formation use shared nucleotide similarity between repeated sequences as substrate for ectopic recombination. We performed genome-wide analyses of direct and inverted intrachromosomal repeated sequence pairs with 200 bp or more and 80% or greater sequence identity in three human genome assemblies, GRCh37, GRCh38, and T2T-CHM13. Overall, the composition and distribution of direct and inverted repeated sequences identified was similar among the three assemblies involving 13%-15% of the haploid genome, with an increased, albeit not significant, number of repeated sequences in T2T-CHM13. Interestingly, the majority of repeated sequences are below 1 kb in length with a median of 84.2% identity, highlighting the potential relevance of smaller, less identical repeats, such as Alu-Alu pairs, for ectopic recombination. We cross-referenced the identified repeated sequences with protein-coding genes to identify those at risk for being involved in genomic rearrangements. Olfactory receptors and immune response genes were enriched among those impacted.

Leveraging the T2T assembly to resolve rare and pathogenic inversions in reference genome gaps

Chromosomal inversions (INVs) are particularly challenging to detect due to their copy-number neutral state and association with repetitive regions. Inversions represent about 1/20 of all balanced structural chromosome aberrations and can lead to disease by gene disruption or altering regulatory regions of dosage-sensitive genes in cis Short-read genome sequencing (srGS) can only resolve ∼70% of cytogenetically visible inversions referred to clinical diagnostic laboratories, likely due to breakpoints in repetitive regions. Here, we study 12 inversions by long-read genome sequencing (lrGS) (n = 9) or srGS (n = 3) and resolve nine of them. In four cases, the inversion breakpoint region was missing from at least one of the human reference genomes (GRCh37, GRCh38, T2T-CHM13) and a reference agnostic analysis was needed. One of these cases, an INV9 mappable only in de novo assembled lrGS data using T2T-CHM13 disrupts EHMT1 consistent with a Mendelian diagnosis (Kleefstra syndrome 1; MIM#610253). Next, by pairwise comparison between T2T-CHM13, GRCh37, and GRCh38, as well as the chimpanzee and bonobo, we show that hundreds of megabases of sequence are missing from at least one human reference, highlighting that primate genomes contribute to genomic diversity. Aligning population genomic data to these regions indicated that these regions are variable between individuals. Our analysis emphasizes that T2T-CHM13 is necessary to maximize the value of lrGS for optimal inversion detection in clinical diagnostics. These results highlight the importance of leveraging diverse and comprehensive reference genomes to resolve unsolved molecular cases in rare diseases.

A fungal plant pathogen overcomes mlo-mediated broad-spectrum disease resistance by rapid gene loss

Hosts and pathogens typically engage in a coevolutionary arms race. This also applies to phytopathogenic powdery mildew fungi, which can rapidly overcome plant resistance and perform host jumps. Using experimental evolution, we show that the powdery mildew pathogen Blumeria hordei is capable of breaking the agriculturally important broad-spectrum resistance conditioned by barley loss-of-function mlo mutants. Partial mlo virulence of evolved B. hordei isolates is correlated with a distinctive pattern of adaptive mutations, including small-sized (c. 8-40 kb) deletions, of which one is linked to the de novo insertion of a transposable element. Occurrence of the mutations is associated with a transcriptional induction of effector protein-encoding genes that is absent in mlo-avirulent isolates on mlo mutant plants. The detected mutational spectrum comprises the same loci in at least two independently isolated mlo-virulent isolates, indicating convergent multigenic evolution. The mutational events emerged in part early (within the first five asexual generations) during experimental evolution, likely generating a founder population in which incipient mlo virulence was later stabilized by additional events. This work highlights the rapid dynamic genome evolution of an obligate biotrophic plant pathogen with a transposon-enriched genome.

AGAP duplicons associate with structural diversity at Chromosome 10q11.22

The 10q11.22 chromosomal region is a duplication-rich interval of the human genome and one of the last to be fully assembled. It carries copy number-variable genes associated with intellectual disability, bipolar disorder, and obesity. In this study, we characterized the structural diversity at this locus by analyzing 64 haploid assemblies produced by the Human Pangenome Reference Consortium. We identified 11 alternative haplotypes that differ in the copy number and/or orientation of large genomic segments, ranging from hundreds of kilobase pairs (kbp) to over one megabase pair (Mbp). We uncovered a 2.4 Mbp size difference between the shortest and longest haplotypes. Breakpoint analysis revealed that genomic instability results from nonallelic homologous recombination between segmental duplication (SD) pairs with varying similarity (94.4%-99.6%). Nonetheless, these pairs generally recombine at positions where their identity is higher (>99.6%). Recurrent inversions occur with different breakpoints within the same inverted SD pair. Inversion polymorphisms shuffle the entire SD arrangement, creating new predispositions to copy-number variations. The SD architecture is associated with a catarrhine-specific subgroup of the AGAP gene family, which likely triggered the accumulation of SDs at this locus over the past 25 million years of human evolution. Our results reveal extensive structural diversity and genomic instability at the 10q11.22 locus, and expand the general understanding of the mutational mechanisms behind SD-mediated rearrangements.

The origin and maintenance of supergenes contributing to ecological adaptation in Atlantic herring

Chromosomal inversions are associated with local adaptation in many species. However, questions regarding how they are formed, maintained and impact various other evolutionary processes remain elusive. Here, using a large genomic dataset of long-read and short-read sequencing, we ask these questions in one of the most abundant vertebrates on Earth, the Atlantic herring. This species has four megabase-sized inversions associated with ecological adaptation that correlate with water temperature. The S and N inversion alleles at these four loci dominate in the southern and northern parts, respectively, of the species distribution in the North Atlantic Ocean. By determining breakpoint coordinates of the four inversions and the structural variations surrounding them, we hypothesize that these inversions are formed by ectopic recombination between duplicated sequences immediately outside of the inversions. We show that these are old inversions (>1 MY), albeit formed after the split between the Atlantic herring and its sister species, the Pacific herring. There is evidence for extensive gene flux between inversion alleles at all four loci. The large Ne of herring combined with the common occurrence of opposite homozygotes across the species distribution has allowed effective purifying selection to prevent the accumulation of genetic load and repeats within the inversions.

Topological stress triggers persistent DNA lesions in ribosomal DNA with ensuing formation of PML-nucleolar compartment

PML, a multifunctional protein, is crucial for forming PML-nuclear bodies involved in stress responses. Under specific conditions, PML associates with nucleolar caps formed after RNA polymerase I (RNAPI) inhibition, leading to PML-nucleolar associations (PNAs). This study investigates PNAs-inducing stimuli by exposing cells to various genotoxic stresses. We found that the most potent inducers of PNAs introduced topological stress and inhibited RNAPI. Doxorubicin, the most effective compound, induced double-strand breaks (DSBs) in the rDNA locus. PNAs co-localized with damaged rDNA, segregating it from active nucleoli. Cleaving the rDNA locus with I-PpoI confirmed rDNA damage as a genuine stimulus for PNAs. Inhibition of ATM, ATR kinases, and RAD51 reduced I-PpoI-induced PNAs, highlighting the importance of ATM/ATR-dependent nucleolar cap formation and homologous recombination (HR) in their triggering. I-PpoI-induced PNAs co-localized with rDNA DSBs positive for RPA32-pS33 but deficient in RAD51, indicating resected DNA unable to complete HR repair. Our findings suggest that PNAs form in response to persistent rDNA damage within the nucleolar cap, highlighting the interplay between PML/PNAs and rDNA alterations due to topological stress, RNAPI inhibition, and rDNA DSBs destined for HR. Cells with persistent PNAs undergo senescence, suggesting PNAs help avoid rDNA instability, with implications for tumorigenesis and aging.

Recurrent evolution and selection shape structural diversity at the amylase locus

The adoption of agriculture triggered a rapid shift towards starch-rich diets in human populations1. Amylase genes facilitate starch digestion, and increased amylase copy number has been observed in some modern human populations with high-starch intake2, although evidence of recent selection is lacking3,4. Here, using 94 long-read haplotype-resolved assemblies and short-read data from approximately 5,600 contemporary and ancient humans, we resolve the diversity and evolutionary history of structural variation at the amylase locus. We find that amylase genes have higher copy numbers in agricultural populations than in fishing, hunting and pastoral populations. We identify 28 distinct amylase structural architectures and demonstrate that nearly identical structures have arisen recurrently on different haplotype backgrounds throughout recent human history. AMY1 and AMY2A genes each underwent multiple duplication/deletion events with mutation rates up to more than 10,000-fold the single-nucleotide polymorphism mutation rate, whereas AMY2B gene duplications share a single origin. Using a pangenome-based approach, we infer structural haplotypes across thousands of humans identifying extensively duplicated haplotypes at higher frequency in modern agricultural populations. Leveraging 533 ancient human genomes, we find that duplication-containing haplotypes (with more gene copies than the ancestral haplotype) have rapidly increased in frequency over the past 12,000 years in West Eurasians, suggestive of positive selection. Together, our study highlights the potential effects of the agricultural revolution on human genomes and the importance of structural variation in human adaptation.

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.

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.

Non-homologous end joining shapes the genomic rearrangement landscape of chromothripsis from mitotic errors

Mitotic errors generate micronuclei entrapping mis-segregated chromosomes, which are susceptible to catastrophic fragmentation through chromothripsis. The reassembly of fragmented chromosomes by error-prone DNA double-strand break (DSB) repair generates diverse genomic rearrangements associated with human diseases. How specific repair pathways recognize and process these lesions remains poorly understood. Here we use CRISPR/Cas9 to systematically inactivate distinct DSB repair pathways and interrogate the rearrangement landscape of fragmented chromosomes. Deletion of canonical non-homologous end joining (NHEJ) components substantially reduces complex rearrangements and shifts the rearrangement landscape toward simple alterations without the characteristic patterns of chromothripsis. Following reincorporation into the nucleus, fragmented chromosomes localize within sub-nuclear micronuclei bodies (MN bodies) and undergo ligation by NHEJ within a single cell cycle. In the absence of NHEJ, chromosome fragments are rarely engaged by alternative end-joining or recombination-based mechanisms, resulting in delayed repair kinetics, persistent 53BP1-labeled MN bodies, and cell cycle arrest. Thus, we provide evidence supporting NHEJ as the exclusive DSB repair pathway generating complex rearrangements from mitotic errors.

Diversity of ribosomes at the level of rRNA variation associated with human health and disease

Ribosomal DNA and RNA (rDNA and rRNA) sequences are usually discarded from sequencing analyses. But with hundreds of copies of rDNA genes it is unknown whether they possess sequence variations that form different types of ribosomes that affect human physiology and disease. Here, we developed an algorithm for variant-calling between paralog genes (termed RGA) and compared rDNA variations found in short- and long-read sequencing data from the 1,000 Genomes Project (1KGP) and Genome In A Bottle (GIAB). We additionally developed a novel protocol for long-read sequencing full-length rRNA (RIBO-RT) from actively translating ribosomes. Our analyses identified hundreds of rDNA variants, most of which, surprisingly, are short insertion-deletions (indels) and dozens of highly abundant rRNA variants that are incorporated into translationally active ribosomes. To visualize variant ribosomes at the single cell level, we developed an in-situ rRNA sequencing method (SWITCH-seq) which revealed that variants are co-expressed within individual cells. Strikingly, by analyzing rDNA, we found that variants assemble into distinct ribosome subtypes. We discovered that these subtypes acquire different rRNA structures by successfully employing dimethyl sulfate (DMS) probing of full length rRNA. With this atlas we investigated rRNA variation changes across human tissues and cancer types. This revealed tissue-specific rRNA subtype expression in endoderm/ectoderm-derived tissues. In cancer, low abundant rRNA variants can become highly expressed, which suggests the presence of cancer-specific ribosomes. Together, this study identifies and comprehensively characterizes the diversity of ribosomes at the level of rRNA variants which is dominated by indel variants, their chromosomal location and unique structure as well as the association of ribosome variation with tissue-specific biology and cancer.

X-linked congenital adrenal hypoplasia: Report of long clinical follow-up and description of a new complex variant in the NR0B1 gene

Adrenal hypoplasia congenita, attributed to NR0B1 pathogenic variants, accounts for more than 50% of the incidence of primary adrenal insufficiency in children. Although more than 250 different deleterious variations have been described, no genotype-phenotype correlation has been defined to date. We report a case of an adopted boy who reported the onset of an adrenal crisis at 2 weeks of age, requiring replacement therapy with mineralocorticoids and glucocorticoids for 4 months. For 3 years, he did well without treatment. At almost 4 years of age, the disorder was restarted. A long follow-up showed the evolution of hypogonadotropic hypogonadism. Molecular studies on NR0B1 revealed a novel and deleterious deletion-insertion-inversion-deletion complex rearrangement sorted in the 5'-3' direction, which is described as follows: (1) deletion of the intergenic region (between TASL and NR0B1 genes) and 5' region, (2) insertion of a sequence containing 37 bp at the junction of the intergenic region of the TASL gene and a part of exon 1 of the NR0B1 gene, (3) inversion of a part of exon 1, (4) deletion of the final portion of exon 1 and exon 2 and beginning of the 3'UTR region, (5) maintenance of part of the intergenic sequence (between genes MAGEB1 and NR0B1, telomeric sense), (6) large posterior deletion, in the same sense. The path to molecular diagnosis was challenging and involved several molecular biology techniques. Evaluating the breakpoints in our patient, we assumed that it was a nonrecurrent rearrangement that had not yet been described. It may involve a repair mechanism known as nonhomologous end-joining (NHEJ), which joins two ends of DNA in an imprecise manner, generating an "information scar," represented herein by the 37 bp insertion. In addition, the local Xp21 chromosome architecture with sequences capable of modifying the DNA structure could impact the formation of complex rearrangements.

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

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

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.

Family case of Potocki-Lupski syndrome

BACKGROUND

Potocki-Lupski syndrome (PTLS, OMIM # 610883) is a rare genetic developmental disorder resulting from a partial heterozygous microduplication at chromosome 17p11.2. The condition is characterized by a wide variability of clinical expression, which can make its clinical and molecular diagnosis challenging.

CASE PRESENTATION

We report here a family (mother and her two children) diagnosed with PTLS. When examining children, neurological and psychological (neuropsychiatric) manifestations (speech delay, mild mental retardation), motor disorders, craniofacial dysmorphism (microcephaly, dolichocephaly, triangular face, wide bulging forehead, long chin, antimongoloid slant, "elfin" ears) were revealed. The suspected clinical diagnosis was confirmed by MLPA and CMA molecular genetic testing which revealed the presence of a segmental aneusomy; microduplication in the 17p11.2 region.

CONCLUSIONS

Children with PTLS can have a clinically recognizable and specific phenotype: craniofacial dysmorphism, motor and neurological manifestations, which may implicate a possible genetic disease to the attending physician. Moreover, each child with this syndrome is unique and may have a different clinical picture. The management of such patients requires a multidisciplinary team approach, including medical genetic counseling.

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.

Case report: A new de novo 6q21q22.1 interstitial deletion case in a girl with cerebellar vermis hypoplasia and developmental delay and literature review

Interstitial deletions involving 6q chromosomal region are rare. Less than 30 patients have been described to date, and fewer have been characterized by high-resolution techniques, such as chromosomal microarray. Deletions involving 6q21q22.1 region are associated with an extremely wide and heterogeneous clinical spectrum, thus genotype-phenotype correlation based on the size of the rearranged region and on the involved genes is complex, even among individuals with overlapping deletions. Here we describe the phenotypic and molecular characterization of a new 6q interstitial deletion in a girl with developmental delay, intellectual disability, cerebellar vermis hypoplasia, facial peculiar characteristics, ataxia and ocular abnormalities. Microarray analysis of the proposita revealed a 7.9 Mb interstitial de novo deletion at 6q21q22.1 chromosomal region, which spanned from nucleotides 108,337,770 to 116,279,453 (GRCh38/hg38). The present case, alongside with a systematic review of the literature, provides further evidence that could aid to the definition of the Smallest Region of Overlap and of the genomic traits that are associated with particular phenotypes, focusing on neurological findings and especially on cerebellar anomalies.

Genome-wide maps of highly-similar intrachromosomal repeats that mediate ectopic recombination in three human genome assemblies

Repeated sequences spread throughout the genome play important roles in shaping the structure of chromosomes and facilitating the generation of new genomic variation. Through a variety of mechanisms, repeats are involved in generating structural rearrangements such as deletions, duplications, inversions, and translocations, which can have the potential to impact human health. Despite their significance, repetitive regions including tandem repeats, transposable elements, segmental duplications, and low-copy repeats remain a challenge to characterize due to technological limitations inherent to many sequencing methodologies. We performed genome-wide analyses and comparisons of direct and inverted repeated sequences in the latest available human genome reference assemblies including GRCh37 and GRCh38 and the most recent telomere-to-telomere alternate assembly (T2T-CHM13). Overall, the composition and distribution of direct and inverted repeats identified remains similar among the three assemblies but we observed an increase in the number of repeated sequences detected in the T2T-CHM13 assembly versus the reference assemblies. As expected, there is an enrichment of repetitive regions in the short arms of acrocentric chromosomes, which had been previously unresolved in the human genome reference assemblies. We cross-referenced the identified repeats with protein-coding genes across the genome to identify those at risk for being involved in genomic disorders. We observed that certain gene categories, such as olfactory receptors and immune response genes, are enriched among those impacted by repeated sequences likely contributing to human diversity and adaptation. Through this analysis, we have produced a catalogue of direct and inversely oriented repeated sequences across the currently three most widely used human genome assemblies. Bioinformatic analyses of these repeats and their contribution to genome architecture can reveal regions that are most susceptible to genomic instability. Understanding how the architectural genomic features of repeat pairs such as their homology, size and distance can lead to complex genomic rearrangement formation can provide further insights into the molecular mechanisms leading to genomic disorders and genome evolution.

Author summary

This study focused on the characterization of intrachromosomal repeated sequences in the human genome that can play important roles in shaping chromosome structure and generating new genomic variation in three human genome assemblies. We observed an increase in the number of repeated sequence pairs detected in the most recent telomere-to-telomere alternate assembly (T2T-CHM13) compared to the reference assemblies (GRCh37 and GRCh38). We observed an enrichment of repeats in the T2T-CHM13 acrocentric chromosomes, which had been previously unresolved. Importantly, our study provides a catalogue of direct and inverted repeated sequences across three commonly used human genome assemblies, which can aid in the understanding of genomic architecture instability, evolution, and disorders. Our analyses provide insights into repetitive regions in the human genome that may contribute to complex genomic rearrangements.

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.

Expanding Genotype-Phenotype Correlation of CLCNKA and CLCNKB Variants Linked to Hearing Loss

The ClC-K channels CLCNKA and CLCNKB are crucial for the transepithelial transport processes required for sufficient urinary concentrations and sensory mechanoelectrical transduction in the cochlea. Loss-of-function alleles in these channels are associated with various clinical phenotypes, ranging from hypokalemic alkalosis to sensorineural hearing loss (SNHL) accompanied by severe renal conditions, i.e., Bartter's syndrome. Using a stepwise genetic approach encompassing whole-genome sequencing (WGS), we identified one family with compound heterozygous variants in the ClC-K channels, specifically a truncating variant in CLCNKA in trans with a contiguous deletion of CLCNKA and CLCNKB. Breakpoint PCR and Sanger sequencing elucidated the breakpoint junctions derived from WGS, and allele-specific droplet digital PCR confirmed one copy loss of the CLCNKA_CLCNKB contiguous deletion. The proband that harbors the CLCNKA_CLCNKB variants is characterized by SNHL without hypokalemic alkalosis and renal anomalies, suggesting a distinct phenotype in the ClC-K channels in whom SNHL predominantly occurs. These results expanded genotypes and phenotypes associated with ClC-K channels, including the disease entities associated with non-syndromic hearing loss. Repeated identification of deletions across various extents of CLCNKA_CLCNKB suggests a mutational hotspot allele, highlighting the need for an in-depth analysis of the CLCNKA_CLCNKB intergenic region, especially in undiagnosed SNHL patients with a single hit in CLCNKA.

Comparative genomic diversity analysis of copy number variations (CNV) in indicine and taurine cattle thriving in Europe and Indian subcontinent

Copy number variations (CNVs) include deletions, duplications, and insertions that are larger than 50 bp in size causing structural variation responsible for diversity, adaptation, and breed development. Indian cattle breeds are highly diverse from the taurine breeds. The pattern of CNVRs in 191 animals belonging to 39 cattle breeds (four Indicine and 35 Taurine) was studied based on Illumina 777K BovineHD chip data. The Indicine breeds revealed 2590 CNVs and 335 copy number variation regions (CNVRs) in autosomes. Out of the identified CNVs, 50 were found to be novel. Structure analysis revealed admixed nature of Siri. Neighbor joining tree from CNVR data showed that hot (Kankrej and Hallikar) and cold (Ladakhi and Siri) adapted cattle breeds clustered separately. CNVR of Indian and European breeds revealed that Balkan and Italian breeds of Podolian group are admixed with Indian cattle breeds corroborating indicine introgression (6.1-13.5%). CNVRs spanning the regions of olfactory receptors and immune system genes were identified. AMOVA revealed 9% variation among populations which is 2% greater than SNP based studies showing higher inclusion of variation by CNVR. Detailed analysis of CNVs/CNVRs in Indian cattle adapted to hot and cold climate, and their diversity among worldwide cattle is presented in this study.

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.

Two sides of the same coin: a complex presentation of autosomal dominant tubulointerstitial kidney diseases: a literature review and case reports

Introduction

Genetic kidney diseases are underdiagnosed; namely, from 7% to 40% of patients suffering from chronic kidney disease (CKD) can carry a pathogenic variant, depending on population characteristics. Hereditary tubulointerstitial kidney diseases, including autosomal dominant tubulointerstitial kidney diseases (ADTKD), are even more challenging to diagnose. ADTKD is a rare form of genetic kidney disease resulting from pathogenic variants in the MUC1, UMOD, HNF1B, REN, SEC61A1, and DNAJB11 genes. There is no typical clinical or histopathological sign of ADTKD, it is characterized by progressive CKD, an autosomal dominant inheritance pattern, and tubular atrophy with interstitial fibrosis on kidney biopsy. There is no significant proteinuria, and the urinary sediment is bland. The patients usually do not have severe arterial hypertension. There can be a history of early gout, especially when compared to the UMOD gene variants. Children can have enuresis due to a loss of renal concentration. On ultrasound, the kidneys can appear normal or small in size. Renal cysts are not pathognomonic for any of the named diseases. End-stage renal disease (ESRD) develops at the average age of 45, but this can be very variable. Family history that suggests autosomal dominant inheritance and CKD fulfilling the aforementioned characteristics of tubulointerstitial kidney disease should raise suspicion of ADTKD. In the setting of a negative family history for CKD, clinical suspicion should be raised based on clinical characteristics, including early onset of hyperuricemia or gout and compatible histology on the kidney biopsy. Contrary to the aforementioned characteristics of ADTKD, in the case of HNF1B-related disease, there is a more complex clinical presentation with extrarenal manifestations of the disease (diabetes mellitus, hypomagnesemia, neurologic and psychiatric disturbances, etc.). The diagnosis of ADTKD is based on a positive family history and a detection of the pathogenic variant in one of the genes in an affected individual.

Aim

The aim of our study is to present two case reports of ADTKD with different characteristics (slowly progressive CKD vs. complex clinical presentation with an extrarenal manifestation of the disease) with a literature review.

Methods

A 34-year-old patient with CKD and a positive family history of CKD in whom kidney biopsy showed nonspecific chronic changes, with only genetic analysis confirming the diagnosis of MUC1-related ADTKD. Our second case is of a 17-year-old patient with an unremarkable family history who was initially referred to genetic counseling due to cognitive and motor impairment with long-lasting epilepsy. Extensive workup revealed increased serum creatinine levels with no proteinuria and bland urinary sediment, along with hypomagnesemia. His genetic analysis revealed 17q12 deletion syndrome, causing the loss of one copy of the HNF1B gene, the AATF, and the LHX1 gene.

Conclusion

Autosomal dominant tubulointerstitial kidney diseases are challenging to diagnose due to a lack of typical clinical or histopathological signs as well as an uncharacteristic and versatile clinical presentation. Increased clinical awareness is crucial for the detection of these diseases.

Structural Variation Evolution at the 15q11-q13 Disease-Associated Locus

The impact of segmental duplications on human evolution and disease is only just starting to unfold, thanks to advancements in sequencing technologies that allow for their discovery and precise genotyping. The 15q11-q13 locus is a hotspot of recurrent copy number variation associated with Prader-Willi/Angelman syndromes, developmental delay, autism, and epilepsy and is mediated by complex segmental duplications, many of which arose recently during evolution. To gain insight into the instability of this region, we characterized its architecture in human and nonhuman primates, reconstructing the evolutionary history of five different inversions that rearranged the region in different species primarily by accumulation of segmental duplications. Comparative analysis of human and nonhuman primate duplication structures suggests a human-specific gain of directly oriented duplications in the regions flanking the GOLGA cores and HERC segmental duplications, representing potential genomic drivers for the human-specific expansions. The increasing complexity of segmental duplication organization over the course of evolution underlies its association with human susceptibility to recurrent disease-associated rearrangements.

Genomic Complexity and Complex Chromosomal Rearrangements in Genetic Diagnosis: Two Illustrative Cases on Chromosome 7

Complex chromosomal rearrangements are rare events compatible with survival, consisting of an imbalance and/or position effect of one or more genes, that contribute to a range of clinical presentations. The investigation and diagnosis of these cases are often difficult. The interpretation of the pattern of pairing and segregation of these chromosomes during meiosis is important for the assessment of the risk and the type of imbalance in the offspring. Here, we investigated two unrelated pediatric carriers of complex rearrangements of chromosome 7. The first case was a 2-year-old girl with a severe phenotype. Conventional cytogenetics evidenced a duplication of part of the short arm of chromosome 7. By array-CGH analysis, we found a complex rearrangement with three discontinuous trisomy regions (7p22.1p21.3, 7p21.3, and 7p21.3p15.3). The second case was a newborn investigated for hypodevelopment and dimorphisms. The karyotype analysis promptly revealed a structurally altered chromosome 7. The array-CGH analysis identified an even more complex rearrangement consisting of a trisomic region at 7q11.23q22 and a tetrasomic region of 4.5 Mb spanning 7q21.3 to q22.1. The mother's karyotype examination revealed a complex rearrangement of chromosome 7: the 7q11.23q22 region was inserted in the short arm at 7p15.3. Finally, array-CGH analysis showed a trisomic region that corresponds to the tetrasomic region of the son. Our work proved that the integration of several technical solutions is often required to appropriately analyze complex chromosomal rearrangements in order to understand their implications and offer appropriate genetic counseling.

Cross-species analysis identifies mitochondrial dysregulation as a functional consequence of the schizophrenia-associated 3q29 deletion

The 1.6-megabase deletion at chromosome 3q29 (3q29Del) is the strongest identified genetic risk factor for schizophrenia, but the effects of this variant on neurodevelopment are not well understood. We interrogated the developing neural transcriptome in two experimental model systems with complementary advantages: isogenic human cortical organoids and isocortex from the 3q29Del mouse model. We profiled transcriptomes from isogenic cortical organoids that were aged for 2 and 12 months, as well as perinatal mouse isocortex, all at single-cell resolution. Systematic pathway analysis implicated dysregulation of mitochondrial function and energy metabolism. These molecular signatures were supported by analysis of oxidative phosphorylation protein complex expression in mouse brain and assays of mitochondrial function in engineered cell lines, which revealed a lack of metabolic flexibility and a contribution of the 3q29 gene PAK2. Together, these data indicate that metabolic disruption is associated with 3q29Del and is conserved across species.

Non-homologous end joining shapes the genomic rearrangement landscape of chromothripsis from mitotic errors

Errors in mitosis can generate micronuclei that entrap mis-segregated chromosomes, which are susceptible to catastrophic fragmentation through a process termed chromothripsis. The reassembly of fragmented chromosomes by error-prone DNA double-strand break (DSB) repair generates a spectrum of simple and complex genomic rearrangements that are associated with human cancers and disorders. How specific DSB repair pathways recognize and process these lesions remains poorly understood. Here we used CRISPR/Cas9 to systematically inactivate distinct DSB processing or repair pathways and interrogated the rearrangement landscape of fragmented chromosomes from micronuclei. Deletion of canonical non-homologous end joining (NHEJ) components, including DNA-PKcs, LIG4, and XLF, substantially reduced the formation of complex rearrangements and shifted the rearrangement landscape toward simple alterations without the characteristic patterns of cancer-associated chromothripsis. Following reincorporation into the nucleus, fragmented chromosomes localize within micronuclei bodies (MN bodies) and undergo successful ligation by NHEJ within a single cell cycle. In the absence of NHEJ, chromosome fragments were rarely engaged by polymerase theta-mediated alternative end-joining or recombination-based mechanisms, resulting in delayed repair kinetics and persistent 53BP1-labeled MN bodies in the interphase nucleus. Prolonged DNA damage signaling from unrepaired fragments ultimately triggered cell cycle arrest. Thus, we provide evidence supporting NHEJ as the exclusive DSB repair pathway generating complex rearrangements following chromothripsis from mitotic errors.

Genomic structural variation: A complex but important driver of human evolution

Structural variants (SVs)-including duplications, deletions, and inversions of DNA-can have significant genomic and functional impacts but are technically difficult to identify and assay compared with single-nucleotide variants. With the aid of new genomic technologies, it has become clear that SVs account for significant differences across and within species. This phenomenon is particularly well-documented for humans and other primates due to the wealth of sequence data available. In great apes, SVs affect a larger number of nucleotides than single-nucleotide variants, with many identified SVs exhibiting population and species specificity. In this review, we highlight the importance of SVs in human evolution by (1) how they have shaped great ape genomes resulting in sensitized regions associated with traits and diseases, (2) their impact on gene functions and regulation, which subsequently has played a role in natural selection, and (3) the role of gene duplications in human brain evolution. We further discuss how to incorporate SVs in research, including the strengths and limitations of various genomic approaches. Finally, we propose future considerations in integrating existing data and biospecimens with the ever-expanding SV compendium propelled by biotechnology advancements.

Applications of advanced technologies for detecting genomic structural variation

Chromosomal structural variation (SV) encompasses a heterogenous class of genetic variants that exerts strong influences on human health and disease. Despite their importance, many structural variants (SVs) have remained poorly characterized at even a basic level, a discrepancy predicated upon the technical limitations of prior genomic assays. However, recent advances in genomic technology can identify and localize SVs accurately, opening new questions regarding SV risk factors and their impacts in humans. Here, we first define and classify human SVs and their generative mechanisms, highlighting characteristics leveraged by various SV assays. We next examine the first-ever gapless assembly of the human genome and the technical process of assembling it, which required third-generation sequencing technologies to resolve structurally complex loci. The new portions of that "telomere-to-telomere" and subsequent pangenome assemblies highlight aspects of SV biology likely to develop in the near-term. We consider the strengths and limitations of the most promising new SV technologies and when they or longstanding approaches are best suited to meeting salient goals in the study of human SV in population-scale genomics research, clinical, and public health contexts. It is a watershed time in our understanding of human SV when new approaches are expected to fundamentally change genomic applications.

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.

Cross-species transcriptomic analysis identifies mitochondrial dysregulation as a functional consequence of the schizophrenia-associated 3q29 deletion

Recent advances in the genetics of schizophrenia (SCZ) have identified rare variants that confer high disease risk, including a 1.6 Mb deletion at chromosome 3q29 with a staggeringly large effect size (O.R. > 40). Understanding the impact of the 3q29 deletion (3q29Del) on the developing CNS may therefore lead to insights about the pathobiology of schizophrenia. To gain clues about the molecular and cellular perturbations caused by the 3q29 deletion, we interrogated transcriptomic effects in two experimental model systems with complementary advantages: isogenic human forebrain cortical organoids and isocortex from the 3q29Del mouse model. We first created isogenic lines by engineering the full 3q29Del into an induced pluripotent stem cell line from a neurotypical individual. We profiled transcriptomes from isogenic cortical organoids that were aged for 2 months and 12 months, as well as day p7 perinatal mouse isocortex, all at single cell resolution. Differential expression analysis by genotype in each cell-type cluster revealed that more than half of the differentially expressed genes identified in mouse cortex were also differentially expressed in human cortical organoids, and strong correlations were observed in mouse-human differential gene expression across most major cell-types. We systematically filtered differentially expressed genes to identify changes occurring in both model systems. Pathway analysis on this filtered gene set implicated dysregulation of mitochondrial function and energy metabolism, although the direction of the effect was dependent on developmental timepoint. Transcriptomic changes were validated at the protein level by analysis of oxidative phosphorylation protein complexes in mouse brain tissue. Assays of mitochondrial function in human heterologous cells further confirmed robust mitochondrial dysregulation in 3q29Del cells, and these effects are partially recapitulated by ablation of the 3q29Del gene PAK2 . Taken together these data indicate that metabolic disruption is associated with 3q29Del and is conserved across species. These results converge with data from other rare SCZ-associated variants as well as idiopathic schizophrenia, suggesting that mitochondrial dysfunction may be a significant but overlooked contributing factor to the development of psychotic disorders. This cross-species scRNA-seq analysis of the SCZ-associated 3q29 deletion reveals that this copy number variant may produce early and persistent changes in cellular metabolism that are relevant to human neurodevelopment.

Diagnostic Value of Chromosomal Microarray Analysis for Fetal Congenital Heart Defects with Different Cardiac Phenotypes and Extracardiac Abnormalities

(1) Background: The objective of this study was to investigate the diagnostic value of chromosomal microarray analysis (CMA) for congenital heart defects (CHDs) with different cardiac phenotypes and extracardiac abnormalities (ECAs) and to explore the pathogenic genetic factors of CHDs. (2) Methods: We collected fetuses diagnosed with CHDs by echocardiography at our hospital from January 2012 to December 2021. We analyzed the CMA results of 427 fetuses with CHDs. We then categorized the CHD into different groups according to two dimensions: different cardiac phenotypes and whether it was combined with ECAs. The correlation between the numerical chromosomal abnormalities (NCAs) and copy number variations (CNVs) with CHDs was analyzed. Statistical analyses, including Chi-square tests and t-tests, were performed on the data using IBM SPSS and GraphPad Prism. (3) Results: In general, CHDs with ECAs increased the detection rate for CA, especially the conotruncal defects. CHD combined with the thoracic and abdominal walls and skeletal, thymic and multiple ECAs, were more likely to exhibit CA. Among the CHD phenotypes, VSD and AVSD were associated with NCA, while DORV may be associated with NCA. The cardiac phenotypes associated with pCNVs were IAA (type A and B), RAA, TAPVC, CoA and TOF. In addition, IAA, B, RAA, PS, CoA and TOF were also associated with 22q11.2DS. The length distribution of the CNV was not significantly different between each CHD phenotype. We detected twelve CNV syndromes, of which six syndromes may be related to CHDs. The pregnancy outcome in this study suggests that termination of pregnancy with fetal VSD and vascular abnormality is more dependent on genetic diagnosis, whereas the outcome in other phenotypes of CHDs may be associated with other additional factors. (4) Conclusions: CMA examination for CHDs is still necessary. We should identify the existence of fetal ECAs and specific cardiac phenotypes, which are helpful for genetic counseling and prenatal diagnosis.

Quantitative trait locus for calving traits on Bos taurus autosome 18 in Holstein cattle is embedded in a complex genomic region

Although the quantitative trait locus (QTL) on chromosome 18 (BTA18) associated with paternal calving ease and stillbirth in Holstein Friesian cattle and its cross has been known for over 20 years, to our knowledge, the exact causal genetic sequence has yet escaped identification. The aim of this study was to re-examine the region of the published QTL on BTA18 and to investigate the possible reasons behind this elusiveness. For this purpose, we carried out a combined linkage disequilibrium and linkage analysis using genotyping data of 2,697 German Holstein Friesian (HF) animals and subsequent whole-genome sequencing (WGS) data analyses and genome assembly of HF samples. We confirmed the known QTL in the 95% confidence interval of 1.089 Mbp between 58.34 and 59.43 Mbp on BTA18. Additionally, these 4 SNPs in the near-perfect linkage disequilibrium with the QTL haplotype were identified: rs381577268 (on 57,816,137 bp, C/T), rs381878735 (on 59,574,329 bp, A/T), rs464221818 (on 59,329,176 bp, C/T), and rs472502785 (on 59,345,689 bp, T/C). Search for the causal mutation using short and long-read sequences, and methylation data of the BTA18 QTL region did not reveal any candidates though. The assembly showed problems in the region, as well as an abundance of segmental duplications within and around the region. Taking the QTL of BTA18 in Holstein cattle as an example, the data presented in this study comprehensively characterize the genomic features that could also be relevant for other such elusive QTL in various other cattle breeds and livestock species as well.

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.

Accurate and Effective Detection of Recurrent Copy Number Variants in Large SNP Genotype Datasets

Structural variations, including recurrent Copy Number Variants (CNVs) at specific genomic loci, have been found to be associated with increased risk of several diseases and syndromes. CNV carrier status can be determined in large collections of samples using SNP arrays and, more recently, sequencing data. Although there is some consensus among researchers about the essential steps required in such analysis (i.e., CNV calling, filtering of putative carriers, and visual validation using intensity data plots of the genomic region), standard methodologies and processes to control the quality and consistency of the results are lacking. Here, we present a comprehensive and user-friendly protocol that we have refined from our extensive research experience in the field. We cover every aspect of the analysis, from input data curation to final results. For each step, we highlight which parameters affect the analysis the most and how different settings may lead to different results. We provide a pipeline to run the complete analysis with effective (but customizable) pre-sets. We present software that we developed to better handle and filter putative CNV carriers and perform visual inspection to validate selected candidates. Finally, we describe methods to evaluate the critical sections and actions to counterbalance potential problems. The current implementation is focused on Illumina SNP array data. All the presented software is freely available and provided in a ready-to-use docker container. © 2022 The Authors. Current Protocols published by Wiley Periodicals LLC. Basic Protocol 1: From raw intensity data files to CNV calls Basic Protocol 2: From CNV calls to validated CNV carriers. Basic Protocol 3: Quality control and quality assessment Basic Protocol 4: Install the necessary software.

Identifying Genetic Etiology in Patients with Intellectual Disability: An Experience in Public Health Services in Northeastern Brazil

Intellectual disability (ID) is considered a common neuropsychiatric disorder that affects up to 3% of the population. The etiologic origin of ID may be genetic, environmental, and multifactorial. Chromosomopathies are relatively common among the genetic causes of ID, especially in the most severe cases and those associated with dysmorphic features. Currently, the application of new molecular cytogenetics technologies has increasingly allowed the identification of microdeletions, microduplications, and unbalanced translocations as causes of ID. The objective of this study was to investigate the etiology of ID in patients admitted to a public hospital in Northeastern Brazil. In total, 119 patients with ID who had normal karyotypes and fragile X exams participated in this study. The patients were initially physically examined for microdeletion syndromes and then tested using fluorescence in situ hybridization (FISH), multiplex ligation-dependent probe amplification (MLPA), methylation-sensitive polymerase chain reaction (MS-PCR), and chromosome microarray analysis (CMA), according to clinical suspicion. Patients with no diagnoses after FISH, MLPA, and/or MS-PCR evaluations were subsequently tested by CMA. The rate of etiologic diagnoses of ID in the current study was 28%. FISH diagnosed 25 out of 79 tested (31%), MLPA diagnosed 26 out of 79 tested (32%), MS-PCR diagnosed 7 out of 20 tested (35%), and the single nucleotide polymorphism array diagnosed 6 out of 27 tested (22%). Although the CMA is the most complete and recommended tool for the diagnosis of microdeletions, microduplications, and unbalance translocations in patients with ID, FISH, MLPA, and MS-PCR testing can be used as the first tests for specific syndromes, as long as the patients are first physically screened clinically, especially in the public health networks system in Brazil, where resources are scarce.

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.

A DNA Replication Mechanism Can Explain Structural Variation at the Pigeon Recessive Red Locus

For species to adapt to their environment, evolution must act upon genetic variation that is present in the population. Elucidating the molecular mechanisms that give rise to this variation is thus of crucial importance for understanding how organisms evolve. In addition to variation caused by point mutations, structural variation (deletions, duplications, inversions, translocations) is also an important source of variety. Mechanisms involving recombination, transposition and retrotransposition, and replication have been proposed for generating structural variation, and each are capable of explaining certain rearrangements. In this study, we conduct a detailed analysis of two partially overlapping rearrangements (e1 and e2 allele) in domestic rock pigeon (Columba livia) which are both associated with the recessive red phenotype. We find that a replicative mechanism is best able to explain the complex architecture of the e1 allele, and is also compatible with the simpler architecture of the e2 allele as well.

Mutation of Proteolipid Protein 1 Gene: From Severe Hypomyelinating Leukodystrophy to Inherited Spastic Paraplegia

Pelizaeus-Merzbacher Disease (PMD) is an inherited leukodystrophy affecting the central nervous system (CNS)-a rare disorder that especially concerns males. Its estimated prevalence is 1.45-1.9 per 100,000 individuals in the general population. Patients affected by PMD exhibit a drastic reduction or absence of myelin sheaths in the white matter areas of the CNS. The Proteolipid Protein 1 (PLP1) gene encodes a transmembrane proteolipid protein. PLP1 is the major protein of myelin, and it plays a key role in the compaction, stabilization, and maintenance of myelin sheaths. Its function is predominant in oligodendrocyte development and axonal survival. Mutations in the PLP1 gene cause the development of a wide continuum spectrum of leukopathies from the most severe form of PMD for whom patients exhibit severe CNS hypomyelination to the relatively mild late-onset type 2 spastic paraplegia, leading to the concept of PLP1-related disorders. The genetic diversity and the biochemical complexity, along with other aspects of PMD, are discussed to reveal the obstacles that hinder the development of treatments. This review aims to provide a clinical and mechanistic overview of this spectrum of rare diseases.

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.

Comprehensive evaluation of structural variant genotyping methods based on long-read sequencing data

BACKGROUND

Structural variants (SVs) play a crucial role in gene regulation, trait association, and disease in humans. SV genotyping has been extensively applied in genomics research and clinical diagnosis. Although a growing number of SV genotyping methods for long reads have been developed, a comprehensive performance assessment of these methods has yet to be done.

RESULTS

Based on one simulated and three real SV datasets, we performed an in-depth evaluation of five SV genotyping methods, including cuteSV, LRcaller, Sniffles, SVJedi, and VaPoR. The results show that for insertions and deletions, cuteSV and LRcaller have similar F1 scores (cuteSV, insertions: 0.69-0.90, deletions: 0.77-0.90 and LRcaller, insertions: 0.67-0.87, deletions: 0.74-0.91) and are superior to other methods. For duplications, inversions, and translocations, LRcaller yields the most accurate genotyping results (0.84, 0.68, and 0.47, respectively). When genotyping SVs located in tandem repeat region or with imprecise breakpoints, cuteSV (insertions and deletions) and LRcaller (duplications, inversions, and translocations) are better than other methods. In addition, we observed a decrease in F1 scores when the SV size increased. Finally, our analyses suggest that the F1 scores of these methods reach the point of diminishing returns at 20× depth of coverage.

CONCLUSIONS

We present an in-depth benchmark study of long-read SV genotyping methods. Our results highlight the advantages and disadvantages of each genotyping method, which provide practical guidance for optimal application selection and prospective directions for tool improvement.

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.

Segmental duplication as potential biomarkers for non-invasive prenatal testing of aneuploidies

Background

Segmental duplication (SD) regions are distinct targets for aneuploidy detection owing to the virtual elimination of amplification bias. The difficulty of searching SD sequences for assay design has hampered their applications.

Methods

We developed a computational program, ChAPDes, which integrates SD searching, refinement, and design of specific PCR primer/probe sets in a pipeline to remove most of the manual work. The generated primer/probe sets were first tested in a multiplex multicolour melting curve analysis for the detection of five common aneuploidies. The primer/probe sets were then tested in a digital PCR assay for the detection of trisomy 21. Finally, a digital PCR protocol was established to quantify maternal plasma DNA sequences for the non-invasive prenatal detection of fetal trisomy 21.

Findings

ChAPDes could output 21,772 candidate primer/probe sets for trisomy 13, 18, 21 and sex chromosome aneuploidies within 2 working days. Clinical evaluation of the multiplex multicolour melting curve analysis involving 463 fetal genomic DNA samples revealed a sensitivity of 100% and specificity of 99.64% in comparison with the reference methods. Using the established digital PCR protocol, we correctly identified two trisomy 21 fetuses and thirteen euploid foetuses from the maternal plasma samples.

Interpretation

The combination of ChAPDes with digital PCR detection could facilitate the use of SD as potential biomarkers for the non-invasive prenatal testing of fetal chromosomal aneuploidies.

Clinicoradiological aspects of pontine tegmental cap dysplasia: Case report of a rare hindbrain malformation

Malformations involving the brainstem are very rare and present with a varied spectrum of clinical symptoms due to multiple cranial nerve palsies and pyramidal tract involvement. Of these, pontine tegmental cap dysplasia is a very unusual malformation, characterized by ventral pons hypoplasia and an ectopic dorsal band of tissue, projecting into the fourth ventricle, from dorsal pontine tegmentum. A 4-year-old male child, presenting with left facial nerve palsy, revealed hypoplastic ventral pons and an ectopic structure on magnetic resonance imaging (MRI). The ectopic structure was isointense to pons, arose from the left side of dorsal pontine tegmentum, at pontomedullary junction and protruded into the fourth ventricle, impinging upon the left seventh and eighth cranial nerves. Diffusion tensor imaging (DTI) depicted abnormal white matter tracts in ectopic tissue with absent transverse pontine fibres and abnormal middle and superior cerebellar peduncles. The typical MRI appearance, coupled with DTI, helped us reach an accurate diagnosis of pontine tegmental cap dysplasia, in a setting of neurological dysfunction.

Modelling segmental duplications in the human genome

Background

Segmental duplications (SDs) are long DNA sequences that are repeated in a genome and have high sequence identity. In contrast to repetitive elements they are often unique and only sometimes have multiple copies in a genome. There are several well-studied mechanisms responsible for segmental duplications: non-allelic homologous recombination, non-homologous end joining and replication slippage. Such duplications play an important role in evolution, however, we do not have a full understanding of the dynamic properties of the duplication process.

Results

We study segmental duplications through a graph representation where nodes represent genomic regions and edges represent duplications between them. The resulting network (the SD network) is quite complex and has distinct features which allow us to make inference on the evolution of segmantal duplications. We come up with the network growth model that explains features of the SD network thus giving us insights on dynamics of segmental duplications in the human genome. Based on our analysis of genomes of other species the network growth model seems to be applicable for multiple mammalian genomes.

Conclusions

Our analysis suggests that duplication rates of genomic loci grow linearly with the number of copies of a duplicated region. Several scenarios explaining such a preferential duplication rates were suggested.

Supplementary Information

The online version contains supplementary material available at (10.1186/s12864-021-07789-7).

Study of complex structural variations of X-linked deafness-2 based on single-molecule sequencing

X-linked deafness-2 (DFNX2) is cochlear incomplete partition type III (IP-III), one of inner ear malformations characterized by an abnormally wide opening in the bone separating the basal turn of the cochlea from the internal auditory canal, fixation of the stapes and cerebrospinal fluid (CSF) gusher upon stapedectomy or cochleostomy. The causative gene of DFNX2 was POU3F4. To investigate the genetic causes of DFNX2 and compare the efficiency of different sequencing methods, 12 unrelated patients were enrolled in the present study. Targeted next-generation sequencing (NGS) and long-read sequencing were used to analyze the genetic etiology of DFNX2. Six variants of POU3F4 were identified in this cohort by NGS. Three patients with a negative diagnosis based on NGS were enrolled in further long-read sequencing. Two of them were all found to carry structural variations (SVs) on chromosome X, consisting of an 870-kb deletion (DEL) at upstream of POU3F4 and an 8-Mb inversion (INV). The 870-kb DEL may have been due to non-homologous end joining (NHEJ), while non-allelic homologous recombination (NAHR) within a single chromatid may have accounted for the 8-Mb INV. Common POU3F4 mutations in DFNX2 included point mutations, small insertions and deletions (INDELs), and exon mutations, which can be detected by Sanger sequencing and NGS. Single-molecule long-read sequencing constitutes an additional and valuable method for accurate detection of pathogenic SVs in IP-III patients with negative NGS results.

The Role of Structural Variation in Adaptation and Evolution of Yeast and Other Fungi

Mutations in DNA can be limited to one or a few nucleotides, or encompass larger deletions, insertions, duplications, inversions and translocations that span long stretches of DNA or even full chromosomes. These so-called structural variations (SVs) can alter the gene copy number, modify open reading frames, change regulatory sequences or chromatin structure and thus result in major phenotypic changes. As some of the best-known examples of SV are linked to severe genetic disorders, this type of mutation has traditionally been regarded as negative and of little importance for adaptive evolution. However, the advent of genomic technologies uncovered the ubiquity of SVs even in healthy organisms. Moreover, experimental evolution studies suggest that SV is an important driver of evolution and adaptation to new environments. Here, we provide an overview of the causes and consequences of SV and their role in adaptation, with specific emphasis on fungi since these have proven to be excellent models to study SV.

Electro-clinical features in epileptic children with chromosome 15q duplication syndrome

OBJECTIVE

We aimed to describe epilepsy and EEG patterns related to vigilance states and age, in chromosome15-long-arm-duplication-syndrome (dup15q) children with epilepsy, in both duplication types: interstitial (intdup15) and isodicentric (idic15).

METHODS

Clinical data and 70 EEGs of 12 patients (5 intdup15, 7 idic15), followed from 4.5 m.o to 17y4m (median follow-up 8y3m), were retrospectively reviewed. EEGs were analyzed visually and using power spectrum analysis.

RESULTS

Seventy video-EEGs were analyzed (1-16 per patient, median 6), follow-up lasting up to 8y10m (median 4y2m): 25 EEGs in intdup15 (8 m.o to 12y.o, median 4y6m) and 45 EEGs in idic15 (7 m.o to 12 y.o, median 15 m). Epilepsy: 6 West syndrome (WS) (2intdup15, 4idic15); 4 Lennox-Gastaut syndromes (LGS) (1 intdup15, 3 idic15), 2 evolving from WS; focal epilepsy (3 intdup15). In idic15, WS displayed additional myoclonic seizures (3), atypical (4) or no hypsarrhythmia (2) and posterior predominant spike and polyspike bursts (4). Beta-band rapid-rhythms (RR): present in 11 patients, power decreased during non-REM-sleep, localization shifted from diffuse to anterior, peak frequency increased with age.

CONCLUSION

WS with peculiar electro-clinical features and LGS, along with beta-band RR decreasing in non-REM-sleep and shifting from diffuse to anterior localization with age are recognizable features pointing towards dup15q diagnosis in children with autism spectrum disorder and developmental delay.

SIGNIFICANCE

This study describes electroclinical features in both interstitial and isodicentric duplications of chromosome 15q, in epileptic children, including some recent extensions regarding sleep features; and illustrates how the temporo-spatial organization of beta oscillations can be of significant help in directing towards dup15q diagnosis hypothesis.