Search

The identification of subnetworks of interest-or active modules-by integrating biological networks with molecular profiles is a key resource to inform on the processes perturbed in different cellular conditions. We here propose MOGAMUN, a Multi-Objective Genetic Algorithm to identify active modules in MUltiplex biological Networks. MOGAMUN optimizes both the density of interactions and the scores of the nodes (e.g., their differential expression). We compare MOGAMUN with state-of-the-art methods, representative of different algorithms dedicated to the identification of active modules in single networks. MOGAMUN identifies dense and high-scoring modules that are also easier to interpret. In addition, to our knowledge, MOGAMUN is the first method able to use multiplex networks. Multiplex networks are composed of different layers of physical and functional relationships between genes and proteins. Each layer is associated to its own meaning, topology, and biases; the multiplex framework allows exploiting this diversity of biological networks. We applied MOGAMUN to identify cellular processes perturbed in Facio-Scapulo-Humeral muscular Dystrophy, by integrating RNA-seq expression data with a multiplex biological network. We identified different active modules of interest, thereby providing new angles for investigating the pathomechanisms of this disease. Availability: MOGAMUN is available at https://github.com/elvanov/MOGAMUN and as a Bioconductor package at https://bioconductor.org/packages/release/bioc/html/MOGAMUN.html. Contact: anais.baudot@univ-amu.fr., Competing Interests: The authors have declared that no competing interests exist.

Background: Facioscapulohumeral muscular dystrophy (FSHD) is an autosomal-dominant disorder and is one of the most common forms of muscular dystrophy. We have recently shown that some hallmarks of FSHD are already expressed in fetal FSHD biopsies, thus opening a new field of investigation for mechanisms leading to FSHD. As microRNAs (miRNAs) play an important role in myogenesis and muscle disorders, in this study we compared miRNAs expression levels during normal and FSHD muscle development., Methods: Muscle biopsies were obtained from quadriceps of both healthy control and FSHD1 fetuses with ages ranging from 14 to 33 weeks of development. miRNA expression profiles were analyzed using TaqMan Human MicroRNA Arrays., Results: During human skeletal muscle development, in control muscle biopsies we observed changes for 4 miRNAs potentially involved in secondary muscle fiber formation and 5 miRNAs potentially involved in fiber maturation. When we compared the miRNA profiles obtained from control and FSHD biopsies, we did not observe any differences in the muscle specific miRNAs. However, we identified 8 miRNAs exclusively expressed in FSHD1 samples (miR-330, miR-331-5p, miR-34a, miR-380-3p, miR-516b, miR-582-5p, miR-517* and miR-625) which could represent new biomarkers for this disease. Their putative targets are mainly involved in muscle development and morphogenesis. Interestingly, these FSHD1 specific miRNAs do not target the genes previously described to be involved in FSHD., Conclusions: This work provides new candidate mechanisms potentially involved in the onset of FSHD pathology. Whether these FSHD specific miRNAs cause deregulations during fetal development, or protect against the appearance of the FSHD phenotype until the second decade of life still needs to be investigated.

Generation of skeletal muscles with forms adapted to their function is essential for normal movement. Muscle shape is patterned by the coordinated polarity of collectively migrating myoblasts. Constitutive inactivation of the protocadherin gene Fat1 uncoupled individual myoblast polarity within chains, altering the shape of selective groups of muscles in the shoulder and face. These shape abnormalities were followed by early onset regionalised muscle defects in adult Fat1-deficient mice. Tissue-specific ablation of Fat1 driven by Pax3-cre reproduced muscle shape defects in limb but not face muscles, indicating a cell-autonomous contribution of Fat1 in migrating muscle precursors. Strikingly, the topography of muscle abnormalities caused by Fat1 loss-of-function resembles that of human patients with facioscapulohumeral dystrophy (FSHD). FAT1 lies near the critical locus involved in causing FSHD, and Fat1 mutant mice also show retinal vasculopathy, mimicking another symptom of FSHD, and showed abnormal inner ear patterning, predictive of deafness, reminiscent of another burden of FSHD. Muscle-specific reduction of FAT1 expression and promoter silencing was observed in foetal FSHD1 cases. CGH array-based studies identified deletion polymorphisms within a putative regulatory enhancer of FAT1, predictive of tissue-specific depletion of FAT1 expression, which preferentially segregate with FSHD. Our study identifies FAT1 as a critical determinant of muscle form, misregulation of which associates with FSHD., Competing Interests: The authors have declared that no competing interests exist.

Mammalian telomeres stabilize chromosome ends as a result of their assembly into a peculiar form of chromatin comprising a complex of non-histone proteins named shelterin. TRF2, one of the shelterin components, binds to the duplex part of telomeric DNA and is essential to fold the telomeric chromatin into a protective cap. Although most of the human telomeric DNA is organized into tightly spaced nucleosomes, their role in telomere protection and how they interplay with telomere-specific factors in telomere organization is still unclear. In this study we investigated whether TRF2 can regulate nucleosome assembly at telomeres.By means of chromatin immunoprecipitation (ChIP) and Micrococcal Nuclease (MNase) mapping assay, we found that the density of telomeric nucleosomes in human cells was inversely proportional to the dosage of TRF2 at telomeres. This effect was not observed in the G1 phase of the cell cycle but appeared coincident of late or post-replicative events. Moreover, we showed that TRF2 overexpression altered nucleosome spacing at telomeres increasing internucleosomal distance. By means of an in vitro nucleosome assembly system containing purified histones and remodeling factors, we reproduced the short nucleosome spacing found in telomeric chromatin. Importantly, when in vitro assembly was performed in the presence of purified TRF2, nucleosome spacing on a telomeric DNA template increased, in agreement with in vivo MNase mapping.Our results demonstrate that TRF2 negatively regulates the number of nucleosomes at human telomeres by a cell cycle-dependent mechanism that alters internucleosomal distance. These findings raise the intriguing possibility that telomere protection is mediated, at least in part, by the TRF2-dependent regulation of nucleosome organization.

The mechanisms governing telomere replication in humans are still poorly understood. To fill this gap, we investigated the timing of replication of single telomeres in human cells. Using in situ hybridization techniques, we have found that specific telomeres have preferential time windows for replication during the S-phase and that these intervals do not depend upon telomere length and are largely conserved between homologous chromosomes and between individuals, even in the presence of large subtelomeric segmental polymorphisms. Importantly, we show that one copy of the 3.3 kb macrosatellite repeat D4Z4, present in the subtelomeric region of the late replicating 4q35 telomere, is sufficient to confer both a more peripheral localization and a later-replicating property to a de novo formed telomere. Also, the presence of beta-satellite repeats next to a newly created telomere is sufficient to delay its replication timing. Remarkably, several native, non-D4Z4-associated, late-replicating telomeres show a preferential localization toward the nuclear periphery, while several early-replicating telomeres are associated with the inner nuclear volume. We propose that, in humans, chromosome arm-specific subtelomeric sequences may influence both the spatial distribution of telomeres in the nucleus and their replication timing., Competing Interests: The authors have declared that no competing interests exist.

Both genetic and epigenetic alterations contribute to Facio-Scapulo-Humeral Dystrophy (FSHD), which is linked to the shortening of the array of D4Z4 repeats at the 4q35 locus. The consequence of this rearrangement remains enigmatic, but deletion of this 3.3-kb macrosatellite element might affect the expression of the FSHD-associated gene(s) through position effect mechanisms. We investigated this hypothesis by creating a large collection of constructs carrying 1 to >11 D4Z4 repeats integrated into the human genome, either at random sites or proximal to a telomere, mimicking thereby the organization of the 4q35 locus. We show that D4Z4 acts as an insulator that interferes with enhancer-promoter communication and protects transgenes from position effect. This last property depends on both CTCF and A-type Lamins. We further demonstrate that both anti-silencing activity of D4Z4 and CTCF binding are lost upon multimerization of the repeat in cells from FSHD patients compared to control myoblasts from healthy individuals, suggesting that FSHD corresponds to a gain-of-function of CTCF at the residual D4Z4 repeats. We propose that contraction of the D4Z4 array contributes to FSHD physio-pathology by acting as a CTCF-dependent insulator in patients., Competing Interests: The authors have declared that no competing interests exist.