Asynchronous replication and autosome-pair non-equivalence in human embryonic stem cells

Allele-specific control of replication timing and genome organization during development

DNA replication occurs in a defined temporal order known as the replication-timing (RT) program. RT is regulated during development in discrete chromosomal units, coordinated with transcriptional activity and 3D genome organization. Here, we derived distinct cell types from F1 hybrid musculus × castaneus mouse crosses and exploited the high single-nucleotide polymorphism (SNP) density to characterize allelic differences in RT (Repli-seq), genome organization (Hi-C and promoter-capture Hi-C), gene expression (total nuclear RNA-seq), and chromatin accessibility (ATAC-seq). We also present HARP, a new computational tool for sorting SNPs in phased genomes to efficiently measure allele-specific genome-wide data. Analysis of six different hybrid mESC clones with different genomes (C57BL/6, 129/sv, and CAST/Ei), parental configurations, and gender revealed significant RT asynchrony between alleles across ∼12% of the autosomal genome linked to subspecies genomes but not to parental origin, growth conditions, or gender. RT asynchrony in mESCs strongly correlated with changes in Hi-C compartments between alleles but not as strongly with SNP density, gene expression, imprinting, or chromatin accessibility. We then tracked mESC RT asynchronous regions during development by analyzing differentiated cell types, including extraembryonic endoderm stem (XEN) cells, four male and female primary mouse embryonic fibroblasts (MEFs), and neural precursor cells (NPCs) differentiated in vitro from mESCs with opposite parental configurations. We found that RT asynchrony and allelic discordance in Hi-C compartments seen in mESCs were largely lost in all differentiated cell types, accompanied by novel sites of allelic asynchrony at a considerably smaller proportion of the genome, suggesting that genome organization of homologs converges to similar folding patterns during cell fate commitment.

Programming asynchronous replication in stem cells

Many regions of the genome replicate asynchronously and are expressed monoallelically. It is thought that asynchronous replication may be involved in choosing one allele over the other, but little is known about how these patterns are established during development. We show that, unlike somatic cells, which replicate in a clonal manner, embryonic and adult stem cells are programmed to undergo switching, such that daughter cells with an early-replicating paternal allele are derived from mother cells that have a late-replicating paternal allele. Furthermore, using ground-state embryonic stem (ES) cells, we demonstrate that in the initial transition to asynchronous replication, it is always the paternal allele that is chosen to replicate early, suggesting that primary allelic choice is directed by preset gametic DNA markers. Taken together, these studies help define a basic general strategy for establishing allelic discrimination and generating allelic diversity throughout the organism.

Random monoallelic expression: regulating gene expression one allele at a time

Monoallelic gene expression is a remarkable process in which transcription occurs from only one of two homologous alleles in a diploid cell. Interestingly, between 0.5% and 15% of autosomal genes exhibit random monoallelic gene expression, in which different cells express only one allele independently of the underlying genomic sequence, in a cell type-specific manner. Recently, genome-wide studies have increased our understanding of the cell type-specific incidence of random monoallelic gene expression, and how the imbalance in allelic expression is distinguished within the cell and potentially maintained across cell generations. Monoallelic gene expression is likely generated through stochastic independent regulation of the two alleles upon differentiation, and has varied implications for the cell and organism, in particular with respect to disease.

Random monoallelic gene expression increases upon embryonic stem cell differentiation

Random autosomal monoallelic gene expression refers to the transcription of a gene from one of two homologous alleles. We assessed the dynamics of monoallelic expression during development through an allele-specific RNA-sequencing screen in clonal populations of hybrid mouse embryonic stem cells (ESCs) and neural progenitor cells (NPCs). We identified 67 and 376 inheritable autosomal random monoallelically expressed genes in ESCs and NPCs, respectively, a 5.6-fold increase upon differentiation. Although DNA methylation and nuclear positioning did not distinguish the active and inactive alleles, specific histone modifications were differentially enriched between the two alleles. Interestingly, expression levels of 8% of the monoallelically expressed genes remained similar between monoallelic and biallelic clones. These results support a model in which random monoallelic expression occurs stochastically during differentiation and, for some genes, is compensated for by the cell to maintain the required transcriptional output of these genes.

Genome-wide association study implicates testis-sperm specific FKBP6 as a susceptibility locus for impaired acrosome reaction in stallions

Impaired acrosomal reaction (IAR) of sperm causes male subfertility in humans and animals. Despite compelling evidence about the genetic control over acrosome biogenesis and function, the genomics of IAR is as yet poorly understood, providing no molecular tools for diagnostics. Here we conducted Equine SNP50 Beadchip genotyping and GWAS using 7 IAR–affected and 37 control Thoroughbred stallions. A significant (P<6.75E-08) genotype–phenotype association was found in horse chromosome 13 in FK506 binding protein 6 (FKBP6). The gene belongs to the immunophilins FKBP family known to be involved in meiosis, calcium homeostasis, clathrin-coated vesicles, and membrane fusions. Direct sequencing of FKBP6 exons in cases and controls identified SNPs g.11040315G>A and g.11040379C>A (p.166H>N) in exon 4 that were significantly associated with the IAR phenotype both in the GWAS cohort (n = 44) and in a large multi-breed cohort of 265 horses. All IAR stallions were homozygous for the A-alleles, while this genotype was found only in 2% of controls. The equine FKBP6 was exclusively expressed in testis and sperm and had 5 different transcripts, of which 4 were novel. The expression of this gene in AC/AG heterozygous controls was monoallelic, and we observed a tendency for FKBP6 up-regulation in IAR stallions compared to controls. Because exon 4 SNPs had no effect on the protein structure, it is likely that FKBP6 relates to the IAR phenotype via regulatory or modifying functions. In conclusion, FKBP6 was considered a susceptibility gene of incomplete penetrance for IAR in stallions and a candidate gene for male subfertility in mammals. FKBP6 genotyping is recommended for the detection of IAR–susceptible individuals among potential breeding stallions. Successful use of sperm as a source of DNA and RNA propagates non-invasive sample procurement for fertility genomics in animals and humans.

Impaired acrosomal reaction (IAR) of sperm causes male subfertility in humans and animals, and currently the molecular causes of the condition are not known. Here we report the mapping, identification, and functional analysis of a susceptibility locus for IAR in stallions. The candidate region was mapped to horse chromosome 13 by SNP genotyping and GWAS of 7 IAR affected and 44 control Thoroughbred stallions. Re-sequencing and case-control analysis of functionally relevant candidate genes in the region identified FKBP6 gene as a significantly associated locus. The association was confirmed by genotyping 265 male horses of multiple breeds. FKBP6 belongs to the immunophilins FKBP family known to be involved in meiosis, calcium homeostasis, clathrin-coated vesicles, and membrane fusions. We showed that the equine FKBP6 is exclusively and monoallelically expressed in testis and sperm and has 5 different transcripts, of which 4 were novel. Overall, FKBP6 was considered a susceptibility gene of incomplete penetrance for IAR in stallions and a candidate gene for male subfertility in other mammals. Successful use of sperm as a source of DNA and RNA propagates non-invasive sample procurement for fertility genomics in animals and humans.

Epigenetics of haematopoietic cell development

Cells of the immune system are generated through a developmental cascade that begins in haematopoietic stem cells. During this process, gene expression patterns are programmed in a series of stages that bring about the restriction of cell potential, ultimately leading to the formation of specialized innate immune cells and mature lymphocytes that express antigen receptors. These events involve the regulation of both gene expression and DNA recombination, mainly through the control of chromatin accessibility. In this Review, we describe the epigenetic changes that mediate this complex differentiation process and try to understand the logic of the programming mechanism.

Epigenetic control of recombination in the immune system

Immune receptor gene expression is regulated by a series of developmental events that modify their accessibility in a locus, cell type, stage and allele-specific manner. This is carried out by a programmed combination of many different molecular mechanisms, including region-wide replication timing, changes in nuclear localization, chromatin contraction, histone modification, nucleosome positioning and DNA methylation. These modalities ultimately work by controlling steric interactions between receptor loci and the recombination machinery.

Allelic inactivation of rDNA loci

Human cells contain several hundred ribosomal genes (rDNA) that are clustered into nucleolar organizer regions (NORs) on the short arms of five different acrocentric chromosomes. Only approximately 50% of the gene copies are actually expressed in somatic cells. Here, we used a new cytological technique to demonstrate that rDNA is regulated allelically in a regional manner, with one parental copy of each NOR being repressed in any individual cell. This process is similar to that of X-chromosome inactivation in females. Early in development, one copy of each NOR becomes late-replicating, thus probably marking it for inactivation and subsequent targeted de novo methylation at rDNA promoter regions. Once established, this multichromosomal allelic pattern is then maintained clonally in somatic cells. This pathway may serve as an epigenetic mechanism for controlling the number of available rDNA copies during development.

Predictable dynamic program of timing of DNA replication in human cells

The organization of mammalian DNA replication is poorly understood. We have produced high-resolution dynamic maps of the timing of replication in human erythroid, mesenchymal, and embryonic stem (ES) cells using TimEX, a method that relies on gaussian convolution of massive, highly redundant determinations of DNA copy-number variations during S phase to produce replication timing profiles. We first obtained timing maps of 3% of the genome using high-density oligonucleotide tiling arrays and then extended the TimEX method genome-wide using massively parallel sequencing. We show that in untransformed human cells, timing of replication is highly regulated and highly synchronous, and that many genomic segments are replicated in temporal transition regions devoid of initiation, where replication forks progress unidirectionally from origins that can be hundreds of kilobases away. Absence of initiation in one transition region is shown at the molecular level by single molecule analysis of replicated DNA (SMARD). Comparison of ES and erythroid cells replication patterns revealed that these cells replicate about 20% of their genome in different quarters of S phase. Importantly, we detected a strong inverse relationship between timing of replication and distance to the closest expressed gene. This relationship can be used to predict tissue-specific timing of replication profiles from expression data and genomic annotations. We also provide evidence that early origins of replication are preferentially located near highly expressed genes, that mid-firing origins are located near moderately expressed genes, and that late-firing origins are located far from genes.