Deciphering D4Z4 CpG methylation gradients in fascioscapulohumeral muscular dystrophy using nanopore sequencing.

Fascioscapulohumeral muscular dystrophy (FSHD) is caused by a unique genetic mechanism that relies on contraction and hypomethylation of the D4Z4 macrosatellite array on the Chromosome 4q telomere allowing ectopic expression of the DUX4 gene in skeletal muscle. Genetic analysis is difficult because of the large size and repetitive nature of the array, a nearly identical array on the 10q telomere, and the presence of divergent D4Z4 arrays scattered throughout the genome. Here, we combine nanopore long-read sequencing with Cas9-targeted enrichment of 4q and 10q D4Z4 arrays for comprehensive genetic analysis including determination of the length of the 4q and 10q D4Z4 arrays with base-pair resolution. In the same assay, we differentiate 4q from 10q telomeric sequences, determine A/B haplotype, identify paralogous D4Z4 sequences elsewhere in the genome, and estimate methylation for all CpGs in the array. Asymmetric, length-dependent methylation gradients were observed in the 4q and 10q D4Z4 arrays that reach a hypermethylation point at approximately 10 D4Z4 repeat units, consistent with the known threshold of pathogenic D4Z4 contractions. High resolution analysis of individual D4Z4 repeat methylation revealed areas of low methylation near the CTCF/insulator region and areas of high methylation immediately preceding the DUX4 transcriptional start site. Within the DUX4 exons, we observed a waxing/waning methylation pattern with a 180-nucleotide periodicity, consistent with phased nucleosomes. Targeted nanopore sequencing complements recently developed molecular combing and optical mapping approaches to genetic analysis for FSHD by adding precision of the length measurement, base-pair resolution sequencing, and quantitative methylation analysis.

NSD2 drives t(4;14) myeloma cell dependence on adenylate kinase 2 by diverting one-carbon metabolism to the epigenome.

ABSTRACT: Chromosomal translocation (4;14), an adverse prognostic factor in multiple myeloma (MM), drives overexpression of the histone methyltransferase nuclear receptor binding SET domain protein 2 (NSD2). A genome-wide CRISPR screen in MM cells identified adenylate kinase 2 (AK2), an enzyme critical for high-energy phosphate transfer from the mitochondria, as an NSD2-driven vulnerability. AK2 suppression in t(4;14) MM cells decreased nicotinamide adenine dinucleotide phosphate (NADP[H]) critical for conversion of ribonucleotides to deoxyribonucleosides, leading to replication stress, DNA damage, and apoptosis. Driving a large genome-wide increase in chromatin methylation, NSD2 overexpression depletes S-adenosylmethionine, compromising the synthesis of creatine from its precursor, guanidinoacetate. Creatine supplementation restored NADP(H) levels, reduced DNA damage, and rescued AK2-deficient t(4;14) MM cells. As the creatine phosphate shuttle constitutes an alternative means for mitochondrial high-energy phosphate transport, these results indicate that NSD2-driven creatine depletion underlies the hypersensitivity of t(4;14) MM cells to AK2 loss. Furthermore, AK2 depletion in t(4;14) cells impaired protein folding in the endoplasmic reticulum, consistent with impaired use of mitochondrial adenosine triphosphate (ATP). Accordingly, AK2 suppression increased the sensitivity of MM cells to proteasome inhibition. These findings delineate a novel mechanism in which aberrant transfer of carbon to the epigenome creates a metabolic vulnerability, with direct therapeutic implications for t(4;14) MM.

FoxO1 promotes ovarian cancer by increasing transcription and METTL14-mediated m6A modification of SMC4.

The transcription factor forkhead box protein O1 (FoxO1) is closely related to the occurrence and development of ovarian cancer (OC), however its role and molecular mechanisms remain unclear. Herein, we found that FoxO1 was highly expressed in clinical samples of OC patients and was significantly correlated with poor prognosis. FoxO1 knockdown inhibited the proliferation of OC cells in vitro and in vivo. ChIP-seq combined with GEPIA2 and Kaplan-Meier database analysis showed that structural maintenance of chromosome 4 (SMC4) is a downstream target of FoxO1, and FoxO1 promotes SMC4 transcription by binding to its -1400/-1390 bp promoter. The high expression of SMC4 significantly blocked the tumor inhibition effect of FoxO1 knockdown. Furtherly, FoxO1 increased SMC4 mRNA abundance by transcriptionally activating methyltransferase-like 14 (METTL14) and increasing SMC4 m6A methylation on its coding sequence region. The Cancer Genome Atlas dataset analysis confirmed a significant positive correlation between FoxO1, SMC4, and METTL14 expression in OC. In summary, this study revealed the molecular mechanisms of FoxO1 regulating SMC4 and established a clinical link between the expression of FoxO1/METTL14/SMC4 in the occurrence of OC, thus providing a potential diagnostic target and therapeutic strategy.

High-resolution breakpoint junction mapping of proximally extended D4Z4 deletions in FSHD1 reveals evidence for a founder effect.

Facioscapulohumeral muscular dystrophy (FSHD) is an inherited myopathy clinically characterized by weakness in the facial, shoulder girdle and upper a muscles. FSHD is caused by chromatin relaxation of the D4Z4 macrosatellite repeat, mostly by a repeat contraction, facilitating ectopic expression of DUX4 in skeletal muscle. Genetic diagnosis for FSHD is generally based on the sizing and haplotyping of the D4Z4 repeat on chromosome 4 by Southern blotting (SB), molecular combing or single-molecule optical mapping, which is usually straight forward but can be complicated by atypical rearrangements of the D4Z4 repeat. One of these rearrangements is a D4Z4 proximally extended deletion (DPED) allele, where not only the D4Z4 repeat is partially deleted, but also sequences immediately proximal to the repeat are lost, which can impede accurate diagnosis in all genetic methods. Previously, we identified several DPED alleles in FSHD and estimated the size of the proximal deletions by a complex pulsed-field gel electrophoresis and SB strategy. Here, using the next-generation sequencing, we have defined the breakpoint junctions of these DPED alleles at the base pair resolution in 12 FSHD families and 4 control individuals facilitating a PCR-based diagnosis of these DPED alleles. Our resultsshow that half of the DPED alleles are derivates of an ancient founder allele. For some DPED alleles, we found that genetic elements are deleted such as DUX4c, FRG2, DBE-T and myogenic enhancers necessitating re-evaluation of their role in FSHD pathogenesis.

Prenatal Diagnosis of Fetuses with 4q35 Deletion: Case Series and Review of the Literature.

INTRODUCTION: 4q35 deletion is a rare chromosomal syndrome with a wide range of phenotypes, which can be challenging to detect through prenatal ultrasound. This study aimed to summarize the fetal phenotypes of patients with 4q35 deletion. CASE PRESENTATION: The study included four fetuses with 4q35 deletion, with detailed records of prenatal ultrasound and genetic testing results. These cases included following phenotypes, fetal growth restriction (FGR) (2/4), cystic hygroma (2/4), single umbilical artery (1/4), and fused kidney (1/4). One case was terminated, while the other three were born and showed no obvious abnormalities at the 1-year follow-up. Previous reports have described the fetal phenotype of 4q35 deletion in 6 patients from five families, with prenatal phenotypes including FGR (2/6), cardiac structural abnormalities (1/6), brain ventriculomegaly (1/6), oligohydramnios (1/6), and multicystic dysplastic kidneys (1/6). CONCLUSION: Overall, the phenotypes of fetuses with 4q35 deletion are diverse, with FGR potentially being a significant phenotype in these cases.

Copy number variants at 4q31.3 affecting the regulatory region of FBXW7 associated with neurodevelopmental delay.

Emerging research has demonstrated that genomic alterations disrupting topologically associated domains (TADs) and chromatin interactions underlie the pathogenic mechanisms of specific copy number variants (CNVs) in neurodevelopmental disorders. We report two patients with a de novo deletion and a duplication in chromosome 4q31, potentially causing FBX-related neurodevelopmental syndrome by affecting the regulatory region of FBXW7. High-throughput chromosome conformation capture (Hi-C) analysis using available capture data in neural progenitor cells revealed the rewiring of the TAD boundary close to FBXW7. Both patients exhibited facial dysmorphisms, cardiac and limb abnormalities, and neurodevelopmental delays, showing significant clinical overlap with previously reported FBXW7-related features. We also included an additional 10 patients with CNVs in the 4q31 region from the literature and the DECIPHER database for Hi-C analysis, which confirmed that disruption of the regulatory region of FBXW7 likely contributes to the developmental defects observed in these patients.

Characterization of D4Z4 alleles and assessment of de novo cases in Facioscapulohumeral dystrophy (FSHD) in a cohort of Italian families.

Facioscapulohumeral dystrophy (FSHD) is an autosomal dominant disease, although 10%-30% of cases are sporadic. However, this percentage may include truly de novo patients (carrying a reduced D4Z4 allele that is not present in either of the parents) and patients with apparently sporadic disease resulting from mosaicism, non-penetrance, or complex genetic situations in either patients or parents. In this study, we characterized the D4Z4 Reduced Alleles (DRA) and evaluated the frequency of truly de novo cases in FSHD1 in a cohort of DNA samples received consecutively for FSHD-diagnostic from 100 Italian families. The D4Z4 testing revealed that 60 families reported a DRA compatible with FSHD1 (1-10 RU). The DRA co-segregated with the disease in most cases. Five families with truly de novo cases were identified, suggesting that this condition may be slightly lower (8%) than previously reported. In addition, D4Z4 characterization in the investigated families showed 4% of mosaic cases and 2% with translocations. This study further highlighted the importance of performing family studies for clarifying apparently sporadic FSHD cases, with significant implications for genetic counseling, diagnosis, clinical management, and procreative choices for patients and families.

Epigenetic dysregulation of eukaryotic initiation factor 3 subunit E (eIF3E) by lysine methyltransferase REIIBP confers a pro-inflammatory phenotype in t(4;14) myeloma.

REIIBP is a lysine methyltransferase aberrantly expressed through alternative promoter usage of NSD2 locus in t(4;14)-translocated multiple myeloma (MM). Clinically, t(4;14) translocation is an adverse prognostic factor found in approximately 15% of MM patients. The contribution of REIIBP relative to other NSD2 isoforms as a dependency gene in t(4;14)-translocated MM remains to be evaluated. Here, we demonstrated that despite homology with NSD2, REIIBP displayed distinct substrate specificity by preferentially catalyzing H3K4me3 and H3K27me3, with little activity on H3K36me2. Furthermore, REIIBP was regulated through microRNA by EZH2 in a Dicer-dependent manner, exemplifying a role of REIIBP in SET-mediated H3K27me3. Chromatin immunoprecipitation sequencing revealed chromatin remodeling characterized by changes in genome-wide and loci-specific occupancy of these opposing histone marks, allowing a bidirectional regulation of its target genes. Transcriptomics indicated that REIIBP induced a pro-inflammatory gene signature through upregulation of TLR7, which in turn led to B-cell receptor-independent activation of BTK and driving NFkB-mediated production of cytokines such as IL-6. Activation of this pathway is targetable using Ibrutinib and partially mitigated bortezomib resistance in a REIIBP xenograft model. Mechanistically, REIIBP upregulated TLR7 through eIF3E, and this relied on eIF3E RNA-binding function instead of its canonical protein synthesis activity, as demonstrated by direct binding to the 3'UTR of TLR7 mRNA. Altogether, we provided a rationale that co-existence of different NSD2 isoforms induced diversified oncogenic programs that should be considered in the strategies for t(4;14)-targeted therapy.

Methylation of the 4q35 D4Z4 repeat defines disease status in facioscapulohumeral muscular dystrophy.

Genetic diagnosis of facioscapulohumeral muscular dystrophy (FSHD) remains a challenge in clinical practice as it cannot be detected by standard sequencing methods despite being the third most common muscular dystrophy. The conventional diagnostic strategy addresses the known genetic parameters of FSHD: the required presence of a permissive haplotype, a size reduction of the D4Z4 repeat of chromosome 4q35 (defining FSHD1) or a pathogenic variant in an epigenetic suppressor gene (consistent with FSHD2). Incomplete penetrance and epistatic effects of the underlying genetic parameters as well as epigenetic parameters (D4Z4 methylation) pose challenges to diagnostic accuracy and hinder prediction of clinical severity. In order to circumvent the known limitations of conventional diagnostics and to complement genetic parameters with epigenetic ones, we developed and validated a multistage diagnostic workflow that consists of a haplotype analysis and a high-throughput methylation profile analysis (FSHD-MPA). FSHD-MPA determines the average global methylation level of the D4Z4 repeat array as well as the regional methylation of the most distal repeat unit by combining bisulphite conversion with next-generation sequencing and a bioinformatics pipeline and uses these as diagnostic parameters. We applied the diagnostic workflow to a cohort of 148 patients and compared the epigenetic parameters based on FSHD-MPA to genetic parameters of conventional genetic testing. In addition, we studied the correlation of repeat length and methylation level within the most distal repeat unit with age-corrected clinical severity and age at disease onset in FSHD patients. The results of our study show that FSHD-MPA is a powerful tool to accurately determine the epigenetic parameters of FSHD, allowing discrimination between FSHD patients and healthy individuals, while simultaneously distinguishing FSHD1 and FSHD2. The strong correlation between methylation level and clinical severity indicates that the methylation level determined by FSHD-MPA accounts for differences in disease severity among individuals with similar genetic parameters. Thus, our findings further confirm that epigenetic parameters rather than genetic parameters represent FSHD disease status and may serve as a valuable biomarker for disease status.

Molecular cytogenetic characterization of isolated recurrent 4q35.2 microduplication in Chinese population: a seven-year single-center retrospective study.

BACKGROUND: With the extensive use of chromosomal microarray analysis (CMA), an increasing number of variants of uncertain significance (VOUS) have been detected. The objective of the present study was to elucidate the pathogenicity and clinical variability associated with isolated recurrent 4q35.2 microduplications within the Chinese population. METHODS: The present study involved 14 cases of isolated recurrent 4q35.2 microduplication (including 12 fetuses and 2 cases of pediatric patients) out of 5,188 subjects who sought genetic consultation at our hospital and received CMA detection. WES technology was subsequently utilized to identify additional sequence variants in a patient with multiple clinical anomalies. RESULTS: All 14 cases exhibited isolated recurrent 4q35.2 microduplications spanning a 1.0-Mb region encompassing the ZFP42 gene. Among the 12 fetuses, 11 displayed normal clinical features, while one was born with renal duplication and hydronephrosis. Additionally, in the two pediatric patients, WES was performed for Case 1, who presented with congenital cataracts, severe intellectual disability, and seizures. This patient inherited the 4q35.2 microduplication from his phenotypically normal mother. WES identified a novel NM_000276:c.2042G > T (p.G681V) variant in the OCRL gene, which is associated with Lowe syndrome and may account for the observed phenotypic variability within this family. CONCLUSION: A series of 14 cases with isolated recurrent 4q35.2 microduplications were investigated, highlighting a potential association with increased susceptibility to renal abnormalities. Further, the present findings may expand the mutation spectrum of the OCRL gene associated with Lowe syndrome and provide valuable insights for the genetic etiological diagnosis of patients with unexplained copy number variants.

Genome-Wide Identification and Characterization of MYB Gene Family and Analysis of Its Sex-Biased Expression Pattern in Spinacia oleracea L.

The members of the myeloblastosis (MYB) family of transcription factors (TFs) participate in a variety of biological regulatory processes in plants, such as circadian rhythm, metabolism, and flower development. However, the characterization of MYB genes across the genomes of spinach Spinacia oleracea L. has not been reported. Here, we identified 140 MYB genes in spinach and described their characteristics using bioinformatics approaches. Among the MYB genes, 54 were 1R-MYB, 80 were 2R-MYB, 5 were 3R-MYB, and 1 was 4R-MYB. Almost all MYB genes were located in the 0-30 Mb region of autosomes; however, the 20 MYB genes were enriched at both ends of the sex chromosome (chromosome 4). Based on phylogeny, conserved motifs, and the structure of genes, 2R-MYB exhibited higher conservation relative to 1R-MYB genes. Tandem duplication and collinearity of spinach MYB genes drive their evolution, enabling the functional diversification of spinach genes. Subcellular localization prediction indicated that spinach MYB genes were mainly located in the nucleus. Cis-acting element analysis confirmed that MYB genes were involved in various processes of spinach growth and development, such as circadian rhythm, cell differentiation, and reproduction through hormone synthesis. Furthermore, through the transcriptome data analysis of male and female flower organs at five different periods, ten candidate genes showed biased expression in spinach males, suggesting that these genes might be related to the development of spinach anthers. Collectively, this study provides useful information for further investigating the function of MYB TFs and novel insights into the regulation of sex determination in spinach.

[Prenatal diagnosis and genetic analysis of two fetuses with Wolf-Hirschhorn syndrome].

OBJETIVE: To explore the prenatal ultrasound phenotype and genetic basis of two fetuses with Wolf-Hirschhorn syndrome (WHS). METHODS: A retrospective analysis was conducted on the ultrasound imaging data of two fetuses suspected for WHS at the Prenatal Diagnostic Center of Qingyuan People's Hospital in July 2017 and August 2019, respectively. Amniotic fluid samples of the two fetuses were subjected to chromosomal karyotyping and chromosomal microarray analysis (CMA). This study was approved by the Qingyuan People's Hospital (Ethics No. IRB-2022-064). RESULTS: Prenatal ultrasound examination of the two fetuses had consistently revealed WHS-associated traits including intrauterine growth restriction (IUGR), craniofacial abnormalities and cardiovascular anomalies. Karyotyping analysis suggested that both fetuses had harbored cryptic chromosomal translocations involving partial deletion of 4p. And parental verification revealed that it was de novo for fetus 1 and paternal for fetus 2. CMA has confirmed that fetus 1 had an approximately 8.7 Mb deletion at 4p16.3p16.1 and a 6.8 Mb duplication at 8p23.1p23.1, whilst fetus 2 had a 20.05 Mb deletion at 4p16.3p15.31 and a 7.66 Mb duplication at 9p24.3p24.1. The karyotype of fetus 1 was determined as 46,XN,der(4)t(4;8)(p16.1;p23.1)dn.arr[hg19]4p16.3p16.1(68345_8721580)×1, 8p23.3p23.1(158048_6933745)×3, and that of fetus 2 was determined as 46,XN,der(4)t(4;9)(p15.3;p24)pat.arr[hg19]4p16.3p15.31(68345_20116061)×1, 9p24.3p24.1(208454_7868292)×3. CONCLUSION: The 4p deletion is probably the main cause for the WHS phenotype in both fetuses. WHS should be suspected when IUGR, renal anomalies, craniofacial and cardiovascular abnormalities are detected upon prenatal ultrasound screening.

De novo variants and recombination at 4q35: Hints for preimplantation genetic testing in facioscapulohumeral muscular dystrophy.

Facioscapulohumeral muscular dystrophy (FSHD) has been associated with the deletion of an integral number of 3.3 kb units of the polymorphic D4Z4 repeat array at 4q35. The prenatal identification of this defect can be carried out on chorionic villi or amniocytes, whereas preimplantation genetic testing for monogenic disorders (PGT-M) requires molecular markers linked to the D4Z4 allele of reduced size. In this context the reliability of this association is crucial. To test the informativeness of the nearby polymorphic markers we investigated recombination at 4q35 using the polymorphic markers D4S1523, D4S163 and D4S139 positioned at 0.55, 0.5 and 0.21 Mb proximal to the D4Z4 array respectively. We determined the probability of recombination events to occur in the D4Z4-D4S1523 interval considering 86 subjects belonging to 12 FSHD families and found a recombination frequency of 14% between D4Z4 and D4S1523. Our study also revealed the occurrence of de novo variants and germline mosaicism. These findings highlight the recombinogenic nature of the 4q subtelomere and indicate that caution should be taken when interpreting PGT-M results. It is advisable that a woman who underwent a PGT-M cycle undertakes a prenatal DNA analysis to confirm the size of the D4Z4 alleles carried by the fetus.

Chromosome 10q-linked FSHD identifies DUX4 as principal disease gene.

BACKGROUND: Facioscapulohumeral dystrophy (FSHD) is an inherited muscular dystrophy clinically characterised by muscle weakness starting with the facial and upper extremity muscles. A disease model has been developed that postulates that failure in somatic repression of the transcription factor DUX4 embedded in the D4Z4 repeat on chromosome 4q causes FSHD. However, due to the position of the D4Z4 repeat close to the telomere and the complex genetic and epigenetic aetiology of FSHD, there is ongoing debate about the transcriptional deregulation of closely linked genes and their involvement in FSHD. METHOD: Detailed genetic characterisation and gene expression analysis of patients with clinically confirmed FSHD and control individuals. RESULTS: Identification of two FSHD families in which the disease is caused by repeat contraction and DUX4 expression from chromosome 10 due to a de novo D4Z4 repeat exchange between chromosomes 4 and 10. We show that the genetic lesion causal to FSHD in these families is physically separated from other candidate genes on chromosome 4. We demonstrate that muscle cell cultures from affected family members exhibit the characteristic molecular features of FSHD, including DUX4 and DUX4 target gene expression, without showing evidence for transcriptional deregulation of other chromosome 4-specific candidate genes. CONCLUSION: This study shows that in rare situations, FSHD can occur on chromosome 10 due to an interchromosomal rearrangement with the FSHD locus on chromosome 4q. These findings provide further evidence that DUX4 derepression is the dominant disease pathway for FSHD. Hence, therapeutic strategies should focus on DUX4 as the primary target.

[Delineation of a mosaicism fetal supernumerary marker chromosome with combined genetic techniques].

OBJECTIVE: To delineate the origin and content of a mosaicism small supernumerary marker chromosome (sSMC) in a fetus with combined chromosomal karyotyping, chromosomal microarray analysis (CMA) and fluorescence in situ hybridization (FISH). METHODS: The fetus of a 31-year-old pregnant woman who had presented at the Maternal and Child Health Care Hospital of Longhua District of Shenzhen City in 2022 was selected as the study subject. Non-invasive prenatal testing suggested that the fetus has harbored a 8.75 Mb duplication in 4q12q13.1. With informed consent, amniotic fluid and peripheral blood samples were taken from the couple for chromosomal karyotyping analysis. The origin and content of a sSMC was identified by CMA, and its proportion in amniotic fluid was determined with a FISH assay. RESULTS: The karyotypes of the pregnant woman, her husband and the fetus were respectively determined as 46,XX, 46,XY,inv(9)(p12q12), and 47,XY,inv(9)(p12q12)pat,+mar[75]/ 46,XY,inv(9)(p12q12)pat[25]. CMA test of the amniotic fluid sample was arr[hg19]4p11q13.1(48978053_63145931)×3, which revealed no mosaicism. However, FISH analysis showed that 59% of interphase cells from the cultured amniotic fluid sample had contained three signals for the centromere of chromosome 4, whilst 65% of interphase cells from the re-sampled amniotic fluid had three such signals, which confirmed the existence of trisomy 8 mosaicism. CONCLUSION: Chromosomal structural abnormality combined with mosaicism can be delineated with combined chromosomal karyotyping and molecular techniques such as FISH and CMA, which has enabled more accurate counseling for the family.

PAX7 target gene repression associates with FSHD progression and pathology over 1 year.

Facioscapulohumeral muscular dystrophy (FSHD) is a prevalent, inherited skeletal myopathy linked to hypomethylation of the D4Z4 macrosatellite at chromosome 4q35. This epigenetic de-repression permits expression of the transcription factor DUX4, which may drive pathology by direct activation of target genes or through inhibition of the homologous transcription factor PAX7. We demonstrated that PAX7 target gene repression is a superior biomarker of FSHD status compared with DUX4 target gene expression. However, despite importance for clinical trials, there remains no transcriptomic biomarker for FSHD progression. A recent study by Wong et al. [Longitudinal measures of RNA expression and disease activity in FSHD muscle biopsies. Hum. Mol. Genet., 29, 1030-1043] performed MRI, muscle biopsy transcriptomics and histopathology on a cohort of FSHD patients with 1-year follow-up. No significant changes in any biomarkers were reported over this time period. However, the authors did not consider PAX7 target gene repression as a marker of FSHD progression. Here we demonstrate that PAX7 target gene repression increases in these paired FSHD samples from year 1 to year 2 and is thus a marker of FSHD progression over 1 year. Moreover, we show that three validated DUX4 target gene expression biomarkers are not associated with FSHD progression over 1 year. We further confirm that PAX7 target gene repression associates with clinical correlates of FSHD disease activity, measured by MRI and histopathology. Thus, PAX7 target gene repression is a uniquely sensitive biomarker of FSHD progression and pathology, valid over a 1 year time frame, implicating its use in clinical trials.

The Genetics and Epigenetics of Facioscapulohumeral Muscular Dystrophy.

Facioscapulohumeral muscular dystrophy (FSHD), a progressive myopathy that afflicts individuals of all ages, provides a powerful model of the complex interplay between genetic and epigenetic mechanisms of chromatin regulation. FSHD is caused by dysregulation of a macrosatellite repeat, either by contraction of the repeat or by mutations in silencing proteins. Both cases lead to chromatin relaxation and, in the context of a permissive allele, aberrant expression of the DUX4 gene in skeletal muscle. DUX4 is a pioneer transcription factor that activates a program of gene expression during early human development, after which its expression is silenced in most somatic cells. When misexpressed in FSHD skeletal muscle, the DUX4 program leads to accumulated muscle pathology. Epigenetic regulators of the disease locus represent particularly attractive therapeutic targets for FSHD, as many are not global modifiers of the genome, and altering their expression or activity should allow correction of the underlying defect.

Efficient Variant Set Mixed Model Association Tests for Continuous and Binary Traits in Large-Scale Whole-Genome Sequencing Studies.

With advances in whole-genome sequencing (WGS) technology, more advanced statistical methods for testing genetic association with rare variants are being developed. Methods in which variants are grouped for analysis are also known as variant-set, gene-based, and aggregate unit tests. The burden test and sequence kernel association test (SKAT) are two widely used variant-set tests, which were originally developed for samples of unrelated individuals and later have been extended to family data with known pedigree structures. However, computationally efficient and powerful variant-set tests are needed to make analyses tractable in large-scale WGS studies with complex study samples. In this paper, we propose the variant-set mixed model association tests (SMMAT) for continuous and binary traits using the generalized linear mixed model framework. These tests can be applied to large-scale WGS studies involving samples with population structure and relatedness, such as in the National Heart, Lung, and Blood Institute's Trans-Omics for Precision Medicine (TOPMed) program. SMMATs share the same null model for different variant sets, and a virtue of this null model, which includes covariates only, is that it needs to be fit only once for all tests in each genome-wide analysis. Simulation studies show that all the proposed SMMATs correctly control type I error rates for both continuous and binary traits in the presence of population structure and relatedness. We also illustrate our tests in a real data example of analysis of plasma fibrinogen levels in the TOPMed program (n = 23,763), using the Analysis Commons, a cloud-based computing platform.

4q-D4Z4 chromatin architecture regulates the transcription of muscle atrophic genes in facioscapulohumeral muscular dystrophy.

Despite increasing insights in genome structure organization, the role of DNA repetitive elements, accounting for more than two thirds of the human genome, remains elusive. Facioscapulohumeral muscular dystrophy (FSHD) is associated with deletion of D4Z4 repeat array below 11 units at 4q35.2. It is known that the deletion alters chromatin structure in cis, leading to gene up-regulation. Here we show a genome-wide role of 4q-D4Z4 array in modulating gene expression via 3D nuclear contacts. We have developed an integrated strategy of 4q-D4Z4-specific 4C-seq and chromatin segmentation analyses, showing that 4q-D4Z4 3D interactome and chromatin states of interacting genes are impaired in FSHD1 condition; in particular, genes that have lost the 4q-D4Z4 interaction and with a more active chromatin state are enriched for muscle atrophy transcriptional signature. Expression level of these genes is restored by the interaction with an ectopic 4q-D4Z4 array, suggesting that the repeat directly modulates the transcription of contacted targets. Of note, the up-regulation of atrophic genes is a common feature of several FSHD1 and FSHD2 patients, indicating that we have identified a core set of deregulated genes involved in FSHD pathophysiology.

Myeloid malignancies with translocation t(4;12)(q11-13;p13): molecular landscape, clonal hierarchy and clinical outcomes.

Translocation t(4;12)(q11-13;p13) is a recurrent but very rare chromosomal aberration in acute myeloid leukaemia (AML) resulting in the non-constant expression of a CHIC2/ETV6 fusion transcript. We report clinico-biological features, molecular characteristics and outcomes of 21 cases of t(4;12) including 19 AML and two myelodysplastic syndromes (MDS). Median age at the time of t(4;12) was 78 years (range, 56-88). Multilineage dysplasia was described in 10 of 19 (53%) AML cases and CD7 and/or CD56 expression in 90%. FISH analyses identified ETV6 and CHIC2 region rearrangements in respectively 18 of 18 and 15 of 17 studied cases. The t(4;12) was the sole cytogenetic abnormality in 48% of cases. The most frequent associated mutated genes were ASXL1 (n = 8/16, 50%), IDH1/2 (n = 7/16, 44%), SRSF2 (n = 5/16, 31%) and RUNX1 (n = 4/16, 25%). Interestingly, concurrent FISH and molecular analyses showed that t(4;12) can be, but not always, a founding oncogenic event. Median OS was 7.8 months for the entire cohort. In the 16 of 21 patients (76%) who received antitumoral treatment, overall response and first complete remission rates were 37% and 31%, respectively. Median progression-free survival in responders was 13.7 months. Finally, t(4;12) cases harboured many characteristics of AML with myelodysplasia-related changes (multilineage dysplasia, MDS-related cytogenetic abnormalities, frequent ASXL1 mutations) and a poor prognosis.