Variants in Candidate Genes for Phenotype Heterogeneity in Patients with the 22q11.2 Deletion Syndrome

22q11.2 deletion syndrome (22q11.2DS) is a microdeletion syndrome with a broad and heterogeneous phenotype, even though most of the deletions present similar sizes, involving ∼3 Mb of DNA. In a relatively large population of a Brazilian 22q11.2DS cohort (60 patients), we investigated genetic variants that could act as genetic modifiers and contribute to the phenotypic heterogeneity, using a targeted NGS (Next Generation Sequencing) with a specific Ion AmpliSeq panel to sequence nine candidate genes (CRKL, MAPK1, HIRA, TANGO2, PI4KA, HDAC1, ZDHHC8, ZFPM2, and JAM3), mapped in and outside the 22q11.2 hemizygous deleted region. In silico prediction was performed, and the whole-genome sequencing annotation analysis package (WGSA) was used to predict the possible pathogenic effect of single nucleotide variants (SNVs). For the in silico prediction of the indels, we used the genomic variants filtered by a deep learning model in NGS (GARFIELD-NGS). We identified six variants, 4 SNVs and 2 indels, in MAPK1, JAM3, and ZFPM2 genes with possibly synergistic deleterious effects in the context of the 22q11.2 deletion. Our results provide the opportunity for the discovery of the co-occurrence of genetic variants with 22q11.2 deletions, which may influence the patients´ phenotype.

Premature ovarian insufficiency is associated with global alterations in the regulatory landscape and gene expression in balanced X-autosome translocations

BACKGROUND

Patients with balanced X-autosome translocations and premature ovarian insufficiency (POI) constitute an interesting paradigm to study the effect of chromosome repositioning. Their breakpoints are clustered within cytobands Xq13-Xq21, 80% of them in Xq21, and usually, no gene disruption can be associated with POI phenotype. As deletions within Xq21 do not cause POI, and since different breakpoints and translocations with different autosomes lead to this same gonadal phenotype, a "position effect" is hypothesized as a possible mechanism underlying POI pathogenesis.

OBJECTIVE AND METHODS

To study the effect of the balanced X-autosome translocations that result in POI, we fine-mapped the breakpoints in six patients with POI and balanced X-autosome translocations and addressed gene expression and chromatin accessibility changes in four of them.

RESULTS

We observed differential expression in 85 coding genes, associated with protein regulation, multicellular regulation, integrin signaling, and immune response pathways, and 120 differential peaks for the three interrogated histone marks, most of which were mapped in high-activity chromatin state regions. The integrative analysis between transcriptome and chromatin data pointed to 12 peaks mapped less than 2 Mb from 11 differentially expressed genes in genomic regions not related to the patients' chromosomal rearrangement, suggesting that translocations have broad effects on the chromatin structure.

CONCLUSION

Since a wide impact on gene regulation was observed in patients, our results observed in this study support the hypothesis of position effect as a pathogenic mechanism for premature ovarian insufficiency associated with X-autosome translocations. This work emphasizes the relevance of chromatin changes in structural variation, since it advances our knowledge of the impact of perturbations in the regulatory landscape within interphase nuclei, resulting in the position effect pathogenicity.

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.

Phenotypic heterogeneity in 22q11.2 deletion syndrome: Copy Number Variants as genetic modifiers for congenital heart disease in a Brazilian cohort

The clinical heterogeneity in 22q11.2 deletion syndrome (22q11.2DS) underlies complex genetic mechanisms including variants in other regions of the genome, known as genetic modifiers. Congenital heart disease (CHD) is one of the most relevant phenotypes in the syndrome and copy number variants (CNVs) outside the 22q11.2 region could play a role in its variable expressivity. Since those described loci account for a small proportion of the variability, the CNV analysis in new cohorts from different ancestry-based populations constitutes a valuable resource to identify a wider range of modifiers. We performed SNP-array in 117 Brazilian patients with 22q11.2DS, with and without CHD, and leveraged genome-wide CNV analysis. After quality control, we selected 50 CNVs in 38 patients for downstream analysis. CNVs' genetic content and implicated biological pathways were compared between patients with and without CHD. CNV-affected genes in patients with CHD were enriched for several functional terms related to ubiquitination, transcription factor binding sites and miRNA targets, highlighting the complexity of the phenotype's expressivity. Cardiac-related genes were identified in both groups of patients suggesting that increasing risk and protective mechanisms could be involved. These genes and enriched pathways could indicate new modifiers to the cardiac phenotype in 22q11.2DS patients.

The impact of 22q11.2 copy-number variants on human traits in the general population

While extensively studied in clinical cohorts, the phenotypic consequences of 22q11.2 copy-number variants (CNVs) in the general population remain understudied. To address this gap, we performed a phenome-wide association scan in 405,324 unrelated UK Biobank (UKBB) participants by using CNV calls from genotyping array. We mapped 236 Human Phenotype Ontology terms linked to any of the 90 genes encompassed by the region to 170 UKBB traits and assessed the association between these traits and the copy-number state of 504 genotyping array probes in the region. We found significant associations for eight continuous and nine binary traits associated under different models (duplication-only, deletion-only, U-shape, and mirror models). The causal effect of the expression level of 22q11.2 genes on associated traits was assessed through transcriptome-wide Mendelian randomization (TWMR), revealing that increased expression of ARVCF increased BMI. Similarly, increased DGCR6 expression causally reduced mean platelet volume, in line with the corresponding CNV effect. Furthermore, cross-trait multivariable Mendelian randomization (MVMR) suggested a predominant role of genuine (horizontal) pleiotropy in the CNV region. Our findings show that within the general population, 22q11.2 CNVs are associated with traits previously linked to genes in the region, and duplications and deletions act upon traits in different fashions. We also showed that gain or loss of distinct segments within 22q11.2 may impact a trait under different association models. Our results have provided new insights to help further the understanding of the complex 22q11.2 region.

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.

Multicentric Carpotarsal Osteolysis Syndrome in a Mother and Daughter with a MAFB Missense Variant and Natural History of the Disease

Multicentric carpotarsal osteolysis syndrome (MCTO; MIM #166300) is a rare skeletal disorder characterized by osteolysis affecting particularly the carpal, metacarpal, and tarsal bones, although other bones might be involved. MCTO is an autosomal dominant disease caused by heterozygous variants in the MAFB gene, frequently misdiagnosed as juvenile rheumatoid arthritis due to similar clinical manifestations. This study reports the first Brazilian family diagnosed with MCTO with progressive osteolysis of the carpal and tarsal bones, presenting a c.161C>T (p.Ser54Leu) heterozygous variant in the MAFB gene, describing the clinical, radiological, and molecular findings, compared with literature data, and discussing the different clinical and molecular diagnosis, as well as the natural history of the disease. Since MCTO is a disorder with progressive symptoms, an early diagnosis is important to avoid unnecessary investigations and treatments and to provide the proper follow-up.

Novel MYT1 variants expose the complexity of oculo-auriculo-vertebral spectrum genetic mechanisms

Oculo-auriculo-vertebral spectrum (OAVS) is a developmental disorder characterized by anomalies mainly involving the structures derived from the first and second pharyngeal arches. The spectrum presents with heterogeneous clinical features and complex etiology with genetic factors not yet completely understood. To date, MYT1 is the most important gene unambiguously associated with the spectrum and with functional data confirmation. In this work, we aimed to identify new single nucleotide variants (SNVs) affecting MYT1 in a cohort of 73 Brazilian patients diagnosed with OAVS. In addition, we investigated copy number variations (CNVs) encompassing this gene or its cis-regulatory elements and compared the frequency of these events in patients versus a cohort of 455 Brazilian control individuals. A new SNV, predicted as likely deleterious, was identified in five unrelated patients with OAVS. All five patients presented hearing impairment and orbital asymmetry suggesting an association with the variant. CNVs near MYT1, located in its neighboring topologically associating domain (TAD), were found to be enriched in patients when compared to controls, indicating a possible involvement of this region with OAVS pathogenicity. Our findings highlight the genetic complexity of the spectrum that seems to involve more than one variant type and inheritance patterns.

Disruption of PCDH10 and TNRC18 Genes due to a Balanced Translocation

Balanced chromosomal rearrangements are usually associated with a normal phenotype, although in some individuals, phenotypic alterations are observed. In these patients, molecular characterization of the breakpoints can reveal the pathogenic mechanism, providing the annotation of disease-associated loci and a better genotype-phenotype correlation. In this study, we describe a patient with a balanced reciprocal translocation between 4q27 and 7p22 associated with neurodevelopmental delay. We performed cytogenetic evaluation, next-generation sequencing of microdissected derivative chromosomes, and Sanger sequencing of the junction points to define the translocation's breakpoints at base pair resolution. We found that the PCDH10 and TNRC18 genes were disrupted by the breakpoints at chromosomes 4 and 7, respectively, with the formation of chimeric genes at the junction points. Gene expression studies in the patient's peripheral blood showed reduced expression of TNRC18, a gene with unknown function and clinical significance. PCDH10 plays a role in the development of the nervous system and might be involved with the patient's neurodevelopmental delay. In this study, the full molecular characterization of the junction points was shown as an efficient tool for fine breakpoint mapping in balanced translocations in order to unmask gene disruptions and investigate the potential pathogenic role of the disrupted genes.

Rare single-nucleotide variants in oculo-auriculo-vertebral spectrum (OAVS)

Background

Oculo‐auriculo‐vertebral spectrum (OAVS) is a craniofacial developmental disorder that affects structures derived from the first and second pharyngeal arches. The clinically heterogeneous phenotype involves mandibular, oral, and ear development anomalies. Etiology is complex and poorly understood. Genetic factors have been associated, evidenced by chromosomal abnormalities affecting different genomic regions and genes. However, known pathogenic single‐nucleotide variants (SNVs) have only been identified in MYT1 in a restricted number of patients. Therefore, investigations of SNVs on candidate genes may reveal other pathogenic mechanisms.

Methods

In a cohort of 73 patients, coding and untranslated regions (UTR) of 10 candidate genes (CRKL, YPEL1, MAPK1, NKX3‐2, HMX1, MYT1, OTX2, GSC, PUF60, HOXA2) were sequenced. Rare SNVs were selected and in silico predictions were performed to ascertain pathogenicity. Likely pathogenic variants were validated by Sanger sequencing and heritability was assessed when possible.

Results

Four likely pathogenic variants in heterozygous state were identified in different patients. Two SNVs were located in the 5’UTR of YPEL1; one in the 3’UTR of CRKL and one in the 3’UTR of OTX2.

Conclusion

Our work described variants in candidate genes for OAVS and supported the genetic heterogeneity of the spectrum.

Inactivation of AMMECR1 is associated with growth, bone, and heart alterations

We report five individuals with loss-of-function of the X-linked AMMECR1: a girl with a balanced X-autosome translocation and inactivation of the normal X-chromosome; two boys with maternally inherited and de novo nonsense variants; and two half-brothers with maternally inherited microdeletion variants. They present with short stature, cardiac and skeletal abnormalities, and hearing loss. Variants of unknown significance in AMMECR1 in four male patients from two families with partially overlapping phenotypes were previously reported. AMMECR1 is coexpressed with genes implicated in cell cycle regulation, five of which were previously associated with growth and bone alterations. Our knockdown of the zebrafish orthologous gene resulted in phenotypes reminiscent of patients' features. The increased transcript and encoded protein levels of AMMECR1L, an AMMECR1 paralog, in the t(X;9) patient's cells indicate a possible partial compensatory mechanism. AMMECR1 and AMMECR1L proteins dimerize and localize to the nucleus as suggested by their nucleic acid-binding RAGNYA folds. Our results suggest that AMMECR1 is potentially involved in cell cycle control and linked to a new syndrome with growth, bone, heart, and kidney alterations with or without elliptocytosis.