Induced Pluripotent Stem Cells for Modeling Physiological and Pathological Striated Muscle Complexity.

Neuromuscular disorders (NMDs) are a large group of diseases associated with either alterations of skeletal muscle fibers, motor neurons or neuromuscular junctions. Most of these diseases is characterized with muscle weakness or wasting and greatly alter the life of patients. Animal models do not always recapitulate the phenotype of patients. The development of innovative and representative human preclinical models is thus strongly needed for modeling the wide diversity of NMDs, characterization of disease-associated variants, investigation of novel genes function, or the development of therapies. Over the last decade, the use of patient's derived induced pluripotent stem cells (hiPSC) has resulted in tremendous progress in biomedical research, including for NMDs. Skeletal muscle is a complex tissue with multinucleated muscle fibers supported by a dense extracellular matrix and multiple cell types including motor neurons required for the contractile activity. Major challenges need now to be tackled by the scientific community to increase maturation of muscle fibers in vitro, in particular for modeling adult-onset diseases affecting this tissue (neuromuscular disorders, cachexia, sarcopenia) and the evaluation of therapeutic strategies. In the near future, rapidly evolving bioengineering approaches applied to hiPSC will undoubtedly become highly instrumental for investigating muscle pathophysiology and the development of therapeutic strategies.

miR-140-5p and miR-140-3p: Key Actors in Aging-Related Diseases?

microRNAs (miRNAs) are small single strand non-coding RNAs and powerful gene expression regulators. They mainly bind to the 3'UTR sequence of targeted mRNA, leading to their degradation or translation inhibition. miR-140 gene encodes the pre-miR-140 that generates the two mature miRNAs miR-140-5p and miR-140-3p. miR-140-5p/-3p have been associated with the development and progression of cancers, but also non-neoplastic diseases. In aging-related diseases, miR-140-5p and miR-140-3p expressions are modulated. The seric levels of these two miRNAs are used as circulating biomarkers and may represent predictive tools. They are also considered key actors in the pathophysiology of aging-related diseases. miR-140-5p/-3p repress targets regulating cell proliferation, apoptosis, senescence, and inflammation. This work focuses on the roles of miR-140-3p and miR-140-5p in aging-related diseases, details their regulation (i.e., by long non-coding RNA), and reviews the molecular targets of theses miRNAs involved in aging pathophysiology.

Mesenchymal stem cells derived from patients with premature aging syndromes display hallmarks of physiological aging.

Progeroid syndromes are rare genetic diseases with most of autosomal dominant transmission, the prevalence of which is less than 1/10,000,000. These syndromes caused by mutations in the <i>LMNA</i> gene encoding A-type lamins belong to a group of disorders called laminopathies. Lamins are implicated in the architecture and function of the nucleus and chromatin. Patients affected with progeroid laminopathies display accelerated aging of mesenchymal stem cells (MSCs)-derived tissues associated with nuclear morphological abnormalities. To identify pathways altered in progeroid patients' MSCs, we used induced pluripotent stem cells (hiPSCs) from patients affected with classical Hutchinson-Gilford progeria syndrome (HGPS, c.1824C&gt;T-p.G608G), HGPS-like syndrome (HGPS-L; c.1868C&gt;G-p.T623S) associated with farnesylated prelamin A accumulation, or atypical progeroid syndromes (APS; homozygous c.1583C&gt; T-p.T528M; heterozygous c.1762T&gt;C-p.C588R; compound heterozygous c.1583C&gt;T and c.1619T&gt;C-p.T528M and p.M540T) without progerin accumulation. By comparative analysis of the transcriptome and methylome of hiPSC-derived MSCs, we found that patient's MSCs display specific DNA methylation patterns and modulated transcription at early stages of differentiation. We further explored selected biological processes deregulated in the presence of <i>LMNA</i> variants and confirmed alterations of age-related pathways during MSC differentiation. In particular, we report the presence of an altered mitochondrial pattern; an increased response to double-strand DNA damage; and telomere erosion in HGPS, HGPS-L, and APS MSCs, suggesting converging pathways, independent of progerin accumulation, but a distinct DNA methylation profile in HGPS and HGPS-L compared with APS cells.

HRAS germline mutations impair LKB1/AMPK signaling and mitochondrial homeostasis in Costello syndrome models.

Germline mutations that activate genes in the canonical RAS/MAPK signaling pathway are responsible for rare human developmental disorders known as RASopathies. Here, we analyzed the molecular determinants of Costello syndrome (CS) using a mouse model expressing HRAS p.G12S, patient skin fibroblasts, hiPSC-derived human cardiomyocytes, a HRAS p.G12V zebrafish model, and human fibroblasts expressing lentiviral constructs carrying HRAS p.G12S or HRAS p.G12A mutations. The findings revealed alteration of mitochondrial proteostasis and defective oxidative phosphorylation in the heart and skeletal muscle of CS mice that were also found in the cell models of the disease. The underpinning mechanisms involved the inhibition of the AMPK signaling pathway by mutant forms of HRAS, leading to alteration of mitochondrial proteostasis and bioenergetics. Pharmacological activation of mitochondrial bioenergetics and quality control restored organelle function in HRAS p.G12A and p.G12S cell models, reduced left ventricle hypertrophy in CS mice, and diminished the occurrence of developmental defects in the CS zebrafish model. Collectively, these findings highlight the importance of mitochondrial proteostasis and bioenergetics in the pathophysiology of RASopathies and suggest that patients with CS may benefit from treatment with mitochondrial modulators.

Facioscapulohumeral dystrophy weakened sarcomeric contractility is mimicked in induced pluripotent stem cells-derived innervated muscle fibres.

BACKGROUND: Facioscapulohumeral dystrophy (FSHD) is a late-onset autosomal dominant form of muscular dystrophy involving specific groups of muscles with variable weakness that precedes inflammatory response, fat infiltration, and muscle atrophy. As there is currently no cure for this disease, understanding and modelling the typical muscle weakness in FSHD remains a major milestone towards deciphering the disease pathogenesis as it will pave the way to therapeutic strategies aimed at correcting the functional muscular defect in patients. METHODS: To gain further insights into the specificity of the muscle alteration in this disease, we derived induced pluripotent stem cells from patients affected with Types 1 and 2 FSHD but also from patients affected with Bosma arhinia and microphthalmia. We differentiated these cells into contractile innervated muscle fibres and analysed their transcriptome by RNA Seq in comparison with cells derived from healthy donors. To uncover biological pathways altered in the disease, we applied MOGAMUN, a multi-objective genetic algorithm that integrates multiplex complex networks of biological interactions (protein-protein interactions, co-expression, and biological pathways) and RNA Seq expression data to identify active modules. RESULTS: We identified 132 differentially expressed genes that are specific to FSHD cells (false discovery rate < 0.05). In FSHD, the vast majority of active modules retrieved with MOGAMUN converges towards a decreased expression of genes encoding proteins involved in sarcomere organization (P value 2.63e-12 ), actin cytoskeleton (P value 9.4e-5 ), myofibril (P value 2.19e-12 ), actin-myosin sliding, and calcium handling (with P values ranging from 7.9e-35 to 7.9e-21 ). Combined with in vivo validations and functional investigations, our data emphasize a reduction in fibre contraction (P value < 0.0001) indicating that the muscle weakness that is typical of FSHD clinical spectrum might be associated with dysfunction of calcium release (P value < 0.0001), actin-myosin interactions, motor activity, mechano-transduction, and dysfunctional sarcomere contractility. CONCLUSIONS: Identification of biomarkers of FSHD muscle remain critical for understanding the process leading to the pathology but also for the definition of readouts to be used for drug design, outcome measures, and monitoring of therapies. The different pathways identified through a system biology approach have been largely overlooked in the disease. Overall, our work opens new perspectives in the definition of biomarkers able to define the muscle alteration but also in the development of novel strategies to improve muscle function as it provides functional parameters for active molecule screening.

AKT Signaling Modifies the Balance between Cell Proliferation and Migration in Neural Crest Cells from Patients Affected with Bosma Arhinia and Microphthalmia Syndrome.

Over the recent years, the SMCHD1 (Structural Maintenance of Chromosome flexible Hinge Domain Containing 1) chromatin-associated factor has triggered increasing interest after the identification of variants in three rare and unrelated diseases, type 2 Facio Scapulo Humeral Dystrophy (FSHD2), Bosma Arhinia and Microphthalmia Syndrome (BAMS), and the more recently isolated hypogonadotrophic hypogonadism (IHH) combined pituitary hormone deficiency (CPHD) and septo-optic dysplasia (SOD). However, it remains unclear why certain mutations lead to a specific muscle defect in FSHD while other are associated with severe congenital anomalies. To gain further insights into the specificity of SMCHD1 variants and identify pathways associated with the BAMS phenotype and related neural crest defects, we derived induced pluripotent stem cells from patients carrying a mutation in this gene. We differentiated these cells in neural crest stem cells and analyzed their transcriptome by RNA-Seq. Besides classical differential expression analyses, we analyzed our data using MOGAMUN, an algorithm allowing the extraction of active modules by integrating differential expression data with biological networks. We found that in BAMS neural crest cells, all subnetworks that are associated with differentially expressed genes converge toward a predominant role for AKT signaling in the control of the cell proliferation-migration balance. Our findings provide further insights into the distinct mechanism by which defects in neural crest migration might contribute to the craniofacial anomalies in BAMS.

SORBS2 transcription is activated by telomere position effect-over long distance upon telomere shortening in muscle cells from patients with facioscapulohumeral dystrophy.

DNA is organized into complex three-dimensional chromatin structures, but how this spatial organization regulates gene expression remains a central question. These DNA/chromatin looping structures can range in size from 10-20 kb (enhancers/repressors) to many megabases during intra- and inter-chromosomal interactions. Recently, the influence of telomere length on chromatin organization prior to senescence has revealed the existence of long-distance chromatin loops that dictate the expression of genes located up to 10 Mb from the telomeres (Telomere Position Effect-Over Long Distances [TPE-OLD]). Here, we demonstrate the existence of a telomere loop at the 4q35 locus involving the sorbin and SH3 domain-containing protein 2 gene, SORBS2, a skeletal muscle protein using a modification of the chromosome conformation capture method. The loop reveals a cis-acting mechanism modifying SORBS2 transcription. The expression of this gene is altered by TPE-OLD in myoblasts from patients affected with the age-associated genetic disease, facioscapulohumeral muscular dystrophy (FSHD1A, MIM 158900). SORBS2 is expressed in FSHD myoblasts with short telomeres, while not detectable in FSHD myoblasts with long telomeres or in healthy myoblasts regardless of telomere length. This indicates that TPE-OLD may modify the regulation of the 4q35 locus in a pathogenic context. Upon differentiation, both FSHD and healthy myotubes express SORBS2, suggesting that SORBS2 is normally up-regulated by maturation/differentiation of skeletal muscle and is misregulated by TPE-OLD-dependent variegation in FSHD myoblasts. These findings provide additional insights for the complexity and age-related symptoms of FSHD.

miRNA expression in control and FSHD fetal human muscle biopsies.

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.

Low penetrance in facioscapulohumeral muscular dystrophy type 1 with large pathological D4Z4 alleles: a cross-sectional multicenter study.

BACKGROUND: Facioscapulohumeral muscular dystrophy type 1(FSHD1) is an autosomal dominant disorder associated with the contraction of D4Z4 less than 11 repeat units (RUs) on chromosome 4q35. Penetrance in the range of the largest alleles is poorly known. Our objective was to study the penetrance of FSHD1 in patients carrying alleles ranging between 6 to10 RUs and to evaluate the influence of sex, age, and several environmental factors on clinical expression of the disease. METHODS: A cross-sectional multicenter study was conducted in six French and one Swiss neuromuscular centers. 65 FSHD1 affected patients carrying a 4qA allele of 6-10 RUs were identified as index cases (IC) and their 119 at-risk relatives were included. The age of onset was recorded for IC only. Medical history, neurological examination and manual muscle testing were performed for each subject. Genetic testing determined the allele size (number of RUs) and the 4qA/4qB allelic variant. The clinical status of relatives was established blindly to their genetic testing results. The main outcome was the penetrance defined as the ratio between the number of clinically affected carriers and the total number of carriers. RESULTS: Among the relatives, 59 carried the D4Z4 contraction. At the clinical level, 34 relatives carriers were clinically affected and 25 unaffected. Therefore, the calculated penetrance was 57% in the range of 6-10 RUs. Penetrance was estimated at 62% in the range of 6-8 RUs, and at 47% in the range of 9-10 RUs. Moreover, penetrance was lower in women than men. There was no effect of drugs, anesthesia, surgery or traumatisms on the penetrance. CONCLUSIONS: Penetrance of FSHD1 is low for largest alleles in the range of 9-10 RUs, and lower in women than men. This is of crucial importance for genetic counseling and clinical management of patients and families.

Telomere position effect: regulation of gene expression with progressive telomere shortening over long distances.

While global chromatin conformation studies are emerging, very little is known about the chromatin conformation of human telomeres. Most studies have focused on the role of telomeres as a tumor suppressor mechanism. Here we describe how telomere length regulates gene expression long before telomeres become short enough to produce a DNA damage response (senescence). We directly mapped the interactions adjacent to specific telomere ends using a Hi-C (chromosome capture followed by high-throughput sequencing) technique modified to enrich for specific genomic regions. We demonstrate that chromosome looping brings the telomere close to genes up to 10 Mb away from the telomere when telomeres are long and that the same loci become separated when telomeres are short. Furthermore, expression array analysis reveals that many loci, including noncoding RNAs, may be regulated by telomere length. We report three genes (ISG15 [interferon-stimulated gene 15 kd], DSP [Desmoplakin], and C1S [complement component 1s subcomplement]) located at three different subtelomeric ends (1p, 6p, and 12p) whose expressions are altered with telomere length. Additionally, we confirmed by in situ analysis (3D-FISH [three-dimensional fluorescence in situ hybridization]) that chromosomal looping occurs between the loci of those genes and their respective telomere ends. We term this process TPE-OLD for "telomere position effect over long distances." Our results suggest a potential novel mechanism for how telomere shortening could contribute to aging and disease initiation/progression in human cells long before the induction of a critical DNA damage response.

Accelerated senescence of cord blood endothelial progenitor cells in premature neonates is driven by SIRT1 decreased expression.

Epidemiological and experimental studies indicate that early vascular dysfunction occurs in low-birth-weight subjects, especially preterm (PT) infants. We recently reported impaired angiogenic activity of endothelial colony-forming cells (ECFCs) in this condition. We hypothesized that ECFC dysfunction in PT might result from premature senescence and investigated the underlying mechanisms. Compared with ECFCs from term neonates (n = 18), ECFCs isolated from PT (n = 29) display an accelerated senescence sustained by growth arrest and increased senescence-associated ß-galactosidase activity. Increased p16(INK4a) expression, in the absence of telomere shortening, indicates that premature PT-ECFC aging results from stress-induced senescence. SIRT1 level, a nicotinamide adenine dinucleotide-dependent deacetylase with anti-aging activities, is dramatically decreased in PT-ECFCs and correlated with gestational age. SIRT1 deficiency is subsequent to epigenetic silencing of its promoter. Transient SIRT1 overexpression or chemical induction by resveratrol treatment reverses senescence phenotype, and rescues in vitro PT-ECFC angiogenic defect in a SIRT1-dependent manner. SIRT1 overexpression also restores PT-ECFC capacity for neovessel formation in vivo. We thus demonstrate that decreased expression of SIRT1 drives accelerated senescence of PT-ECFCs, and acts as a critical determinant of the PT-ECFC angiogenic defect. These findings lay new grounds for understanding the increased cardiovascular risk in individuals born prematurely and open perspectives for therapeutic strategy.

Acacetin and chrysin, two polyphenolic compounds, alleviate telomeric position effect in human cells.

We took advantage of the ability of human telomeres to silence neighboring genes (telomere position effect or TPE) to design a high-throughput screening assay for drugs altering telomeres. We identified, for the first time, that two dietary flavones, acacetin and chrysin, are able to specifically alleviate TPE in human cells. We further investigated their influence on telomere integrity and showed that both drugs drastically deprotect telomeres against DNA damage response. However, telomere deprotection triggered by shelterin dysfunction does not affect TPE, indicating that acacetin and chrysin target several functions of telomeres. These results show that TPE-based screening assays represent valuable methods to discover new compounds targeting telomeres.Molecular Therapy-Nucleic Acids (2013) 2, e116; doi:10.1038/mtna.2013.42; published online 20 August 2013.

Telomere protection and TRF2 expression are enhanced by the canonical Wnt signalling pathway.

The DNA-binding protein TRF2 is essential for telomere protection and chromosome stability in mammals. We show here that TRF2 expression is activated by the Wnt/ß-catenin signalling pathway in human cancer and normal cells as well as in mouse intestinal tissues. Furthermore, ß-catenin binds to TRF2 gene regulatory regions that are functional in a luciferase transactivating assay. Reduced ß-catenin expression in cancer cells triggers a marked increase in telomere dysfunction, which can be reversed by TRF2 overexpression. We conclude that the Wnt/ß-catenin signalling pathway maintains a level of TRF2 critical for telomere protection. This is expected to have an important role during development, adult stem cell function and oncogenesis.

TRF2 controls telomeric nucleosome organization in a cell cycle phase-dependent manner.

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.

TRF2 and apollo cooperate with topoisomerase 2alpha to protect human telomeres from replicative damage.

Human telomeres are protected from DNA damage by a nucleoprotein complex that includes the repeat-binding factor TRF2. Here, we report that TRF2 regulates the 5' exonuclease activity of its binding partner, Apollo, a member of the metallo-beta-lactamase family that is required for telomere integrity during S phase. TRF2 and Apollo also suppress damage to engineered interstitial telomere repeat tracts that were inserted far away from chromosome ends. Genetic data indicate that DNA topoisomerase 2alpha acts in the same pathway of telomere protection as TRF2 and Apollo. Moreover, TRF2, which binds preferentially to positively supercoiled DNA substrates, together with Apollo, negatively regulates the amount of TOP1, TOP2alpha, and TOP2beta at telomeres. Our data are consistent with a model in which TRF2 and Apollo relieve topological stress during telomere replication. Our work also suggests that cellular senescence may be caused by topological problems that occur during the replication of the inner portion of telomeres.

Platination of telomeric DNA by cisplatin disrupts recognition by TRF2 and TRF1.

Telomeres, the nucleoprotein complexes located at the ends of chromosomes, are involved in chromosome protection and genome stability. Telomeric repeat binding factor 1 (TRF1) and telomeric repeat binding factor 2 (TRF2) are the two telomeric proteins that bind to duplex telomeric DNA through interactions between their C-terminal domain and several guanines of the telomeric tract. Since the antitumour drug cisplatin binds preferentially to two adjacent guanines, we have investigated whether cisplatin adducts could affect the binding of TRF1 and TRF2 to telomeric DNA and the property of TRF2 to stimulate telomeric invasion, a process that is thought to participate in the formation of the t-loop. We show that the binding of TRF1 and TRF2 to telomeric sequences selectively modified by one GG chelate of cisplatin is markedly affected by cisplatin but that the effect is more drastic for TRF2 than for TRF1 (3-5-fold more sensitivity for TRF2 than for TRF1). We also report that platinum adducts cause a decrease in TRF2-dependent stimulation of telomeric invasion in vitro. Finally, in accordance with in vitro data, analysis of telomeric composition after cisplatin treatment reveals that 60% of TRF2 dissociate from telomeres.

Identification of a perinuclear positioning element in human subtelomeres that requires A-type lamins and CTCF.

The localization of genes within the nuclear space is of paramount importance for proper genome functions. However, very little is known on the cis-acting elements determining subnuclear positioning of chromosome segments. We show here that the D4Z4 human subtelomeric repeat localizes a telomere at the nuclear periphery. This perinuclear activity lies within an 80 bp sequence included within a region known to interact with CTCF and A-type Lamins. We further show that a reduced level of either CTCF or A-type Lamins suppresses the perinuclear activities of D4Z4 and that an array of multimerized D4Z4 sequence, which has lost its ability to bind CTCF and A-type Lamins, is not localized at the periphery. Overall, these findings reveal the existence of an 80 bp D4Z4 sequence that is sufficient to position an adjacent telomere to the nuclear periphery in a CTCF and A-type lamins-dependent manner. Strikingly, this sequence includes a 30 bp GA-rich motif, which binds CTCF and is present at several locations in the human genome.

Specific binding of the methyl binding domain protein 2 at the BRCA1-NBR2 locus.

The methyl-CpG binding domain (MBD) proteins are key molecules in the interpretation of DNA methylation signals leading to gene silencing. We investigated their binding specificity at the constitutively methylated region of a CpG island containing the bidirectional promoter of the Breast cancer predisposition gene 1, BRCA1, and the Near BRCA1 2 (NBR2) gene. In HeLa cells, quantitative chromatin immunoprecipitation assays indicated that MBD2 is associated with the methylated region, while MeCP2 and MBD1 were not detected at this locus. MBD2 depletion (approximately 90%), mediated by a transgene expressing a small interfering RNA (siRNA), did not induce MeCP2 or MBD1 binding at the methylated area. Furthermore, the lack of MBD2 at the BRCA1-NBR2 CpG island is associated with an elevated level of NBR2 transcripts and with a significant reduction of induced-DNA-hypomethylation response. In MBD2 knockdown cells, transient expression of a Mbd2 cDNA, refractory to siRNA-mediated decay, shifted down the NBR2 mRNA level to that observed in unmodified HeLa cells. Variations in MBD2 levels did not affect BRCA1 expression despite its stimulation by DNA hypomethylation. Collectively, our data indicate that MBD2 has specific targets and its presence at these targets is indispensable for gene repression.

Epigenetic marks at BRCA1 and p53 coding sequences in early human embryogenesis.

In the vertebrate genome, methylation of deoxycytosine residues of CpGs dinucleotide has been associated with transcriptional silencing of genes, parental imprinting, X-inactivation and chromatin remodelling. In human somatic tissues, the 5' end of the BRCA1 CpG island is methylated, whereas this region is unmethylated in mature germ cells and early embryos. In gametes, as in somatic tissues, the CpG sites in the coding region are methylated. We took advantage of this bimodal distribution as a model to analyse the epigenetic reprogramming of coding regions during early human embryogenesis using the bisulphite-based genomic sequencing method. During preimplantation divisions, exon 11 of BRCA1 was slowly demethylated and retained approximately 30% of its methylated residues at the blastocyst stage. Moreover, the change in the distribution of methylated residues was not restricted to the BRCA1 gene, since for another gene, p53, a relatively high level of methylation (50%) of exon 4 was observed in blastocysts. Taken together, these data suggest that a significant part of the methylated residues of coding sequences might be conserved during preimplantation development.