In vitro reconstitution of epigenetic reprogramming in the human germ line.

Epigenetic reprogramming resets parental epigenetic memories and differentiates primordial germ cells (PGCs) into mitotic pro-spermatogonia or oogonia. This process ensures sexually dimorphic germ cell development for totipotency1. In vitro reconstitution of epigenetic reprogramming in humans remains a fundamental challenge. Here we establish a strategy for inducing epigenetic reprogramming and differentiation of pluripotent stem-cell-derived human PGC-like cells (hPGCLCs) into mitotic pro-spermatogonia or oogonia, coupled with their extensive amplification (about >1010-fold). Bone morphogenetic protein (BMP) signalling is a key driver of these processes. BMP-driven hPGCLC differentiation involves attenuation of the MAPK (ERK) pathway and both de novo and maintenance DNA methyltransferase activities, which probably promote replication-coupled, passive DNA demethylation. hPGCLCs deficient in TET1, an active DNA demethylase abundant in human germ cells2,3, differentiate into extraembryonic cells, including amnion, with de-repression of key genes that bear bivalent promoters. These cells fail to fully activate genes vital for spermatogenesis and oogenesis, and their promoters remain methylated. Our study provides a framework for epigenetic reprogramming in humans and an important advance in human biology. Through the generation of abundant mitotic pro-spermatogonia and oogonia-like cells, our results also represent a milestone for human in vitro gametogenesis research and its potential translation into reproductive medicine.

Reconstituting human somitogenesis in vitro.

The segmented body plan of vertebrates is established during somitogenesis, a well-studied process in model organisms; however, the details of this process in humans remain largely unknown owing to ethical and technical limitations. Despite recent advances with pluripotent stem cell-based approaches1-5, models that robustly recapitulate human somitogenesis in both space and time remain scarce. Here we introduce a pluripotent stem cell-derived mesoderm-based 3D model of human segmentation and somitogenesis-which we termed 'axioloid'-that captures accurately the oscillatory dynamics of the segmentation clock and the morphological and molecular characteristics of sequential somite formation in vitro. Axioloids show proper rostrocaudal patterning of forming segments and robust anterior-posterior FGF-WNT signalling gradients and retinoic acid signalling components. We identify an unexpected critical role of retinoic acid signalling in the stabilization of forming segments, indicating distinct, but also synergistic effects of retinoic acid and extracellular matrix on the formation and epithelialization of somites. Comparative analysis demonstrates marked similarities of axioloids to the human embryo, further validated by the presence of a Hox code in axioloids. Finally, we demonstrate the utility of axioloids for studying the pathogenesis of human congenital spine diseases using induced pluripotent stem cells with mutations in HES7 and MESP2. Our results indicate that axioloids represent a promising platform for the study of axial development and disease in humans.

Induction of fetal meiotic oocytes from embryonic stem cells in cynomolgus monkeys.

Human in vitro oogenesis provides a framework for clarifying the mechanism of human oogenesis. To create its benchmark, it is vital to promote in vitro oogenesis using a model physiologically close to humans. Here, we establish a foundation for in vitro oogenesis in cynomolgus (cy) monkeys (Macaca fascicularis): cy female embryonic stem cells harboring one active and one inactive X chromosome (Xa and Xi, respectively) differentiate robustly into primordial germ cell-like cells, which in xenogeneic reconstituted ovaries develop efficiently into oogonia and, remarkably, further into meiotic oocytes at the zygotene stage. This differentiation entails comprehensive epigenetic reprogramming, including Xi reprogramming, yet Xa and Xi remain epigenetically asymmetric with, as partly observed in vivo, incomplete Xi reactivation. In humans and monkeys, the Xi epigenome in pluripotent stem cells functions as an Xi-reprogramming determinant. We further show that developmental pathway over-activations with suboptimal up-regulation of relevant meiotic genes impede in vitro meiotic progression. Cy in vitro oogenesis exhibits critical homology with the human system, including with respect to bottlenecks, providing a salient model for advancing human in vitro oogenesis.

Mechanistic insights into the evolution of the differential expression of tandemly arrayed cone opsin genes in zebrafish.

The genome of many organisms contains several loci consisting of duplicated genes that are arrayed in tandem. The daughter genes produced by duplication typically exhibit differential expression patterns with each other or otherwise experience pseudogenization. Remarkably, opsin genes in fish are preserved after many duplications in different lineages. This fact indicates that fish opsin genes are characterized by a regulatory mechanism that could intrinsically facilitate the differentiation of the expression patterns. However, little is known about the mechanisms that underlie the differential expression patterns or how they were established during evolution. The loci of green (RH2)- and red (LWS)-sensitive cone opsin genes in zebrafish have been used as model systems to study the differential regulation of tandemly arrayed opsin genes. Over a decade of studies have uncovered several mechanistic features that might have assisted the differentiation and preservation of duplicated genes. Furthermore, recent progress in the understanding of the transcriptional process in general has added essential insights. In this article, the current understanding of the transcriptional regulation of differentially expressed tandemly arrayed cone opsin genes in zebrafish is summarized and a possible evolutionary scenario that could achieve this differentiation is discussed.

A discrete transition zone organizes the topological and regulatory autonomy of the adjacent tfap2c and bmp7 genes.

Despite the well-documented role of remote enhancers in controlling developmental gene expression, the mechanisms that allocate enhancers to genes are poorly characterized. Here, we investigate the cis-regulatory organization of the locus containing the Tfap2c and Bmp7 genes in vivo, using a series of engineered chromosomal rearrangements. While these genes lie adjacent to one another, we demonstrate that they are independently regulated by distinct sets of enhancers, which in turn define non-overlapping regulatory domains. Chromosome conformation capture experiments reveal a corresponding partition of the locus in two distinct structural entities, demarcated by a discrete transition zone. The impact of engineered chromosomal rearrangements on the topology of the locus and the resultant gene expression changes indicate that this transition zone functionally organizes the structural partition of the locus, thereby defining enhancer-target gene allocation. This partition is, however, not absolute: we show that it allows competing interactions across it that may be non-productive for the competing gene, but modulate expression of the competed one. Altogether, these data highlight the prime role of the topological organization of the genome in long-distance regulation of gene expression.

Functional and topological characteristics of mammalian regulatory domains.

Long-range regulatory interactions play an important role in shaping gene-expression programs. However, the genomic features that organize these activities are still poorly characterized. We conducted a large operational analysis to chart the distribution of gene regulatory activities along the mouse genome, using hundreds of insertions of a regulatory sensor. We found that enhancers distribute their activities along broad regions and not in a gene-centric manner, defining large regulatory domains. Remarkably, these domains correlate strongly with the recently described TADs, which partition the genome into distinct self-interacting blocks. Different features, including specific repeats and CTCF-binding sites, correlate with the transition zones separating regulatory domains, and may help to further organize promiscuously distributed regulatory influences within large domains. These findings support a model of genomic organization where TADs confine regulatory activities to specific but large regulatory domains, contributing to the establishment of specific gene expression profiles.

Spatially differentiated expression of quadruplicated green-sensitive RH2 opsin genes in zebrafish is determined by proximal regulatory regions and gene order to the locus control region.

BACKGROUND: Fish are remarkably diverse in repertoires of visual opsins by gene duplications. Differentiation of their spatiotemporal expression patterns and absorption spectra enables fine-tuning of feature detection in spectrally distinct regions of the visual field during ontogeny. Zebrafish have quadruplicated green-sensitive (RH2) opsin genes in tandem (RH2-1, -2, -3, -4), which are expressed in the short member of the double cones (SDC). The shortest wavelength RH2 subtype (RH2-1) is expressed in the central to dorsal area of the adult retina. The second shortest wave subtype (RH2-2) is expressed overlapping with RH2-1 but extending outside of it. The second longest wave subtype (RH2-3) is expressed surrounding the RH2-2 area, and the longest wave subtype (RH2-4) is expressed outside of the RH2-3 area broadly occupying the ventral area. Expression of the four RH2 genes in SDC requires a single enhancer (RH2-LCR), but the mechanism of their spatial differentiation remains elusive. RESULTS: Functional comparison of the RH2-LCR with its counterpart in medaka revealed that the regulatory role of the RH2-LCR in SDC-specific expression is evolutionarily conserved. By combining the RH2-LCR and the proximal upstream region of each RH2 gene with fluorescent protein reporters, we show that the RH2-LCR and the RH2-3 proximal regulatory region confer no spatial selectivity of expression in the retina. But those of RH2-1, -2 and -4 are capable of inducing spatial differentiation of expression. Furthermore, by analyzing transgenic fish with a series of arrays consisting of the RH2-LCR and multiple upstream regions of the RH2 genes in different orders, we show that a gene expression pattern related to an upstream region is greatly influenced by another flanking upstream region in a relative position-dependent manner. CONCLUSIONS: The zebrafish RH2 genes except RH2-3 acquired differential cis-elements in the proximal upstream regions to specify the differential expression patterns. The input from these proximal elements collectively dictates the actual gene expression pattern of the locus, context-dependently. Importantly, competition for the RH2-LCR activity among the replicates is critical in this collective regulation, facilitating differentiation of expression among them. This combination of specificity and generality enables seemingly complicated spatial differentiation of duplicated opsin genes characteristic in fish.

TRACER: a resource to study the regulatory architecture of the mouse genome.

BACKGROUND: Mammalian genes are regulated through the action of multiple regulatory elements, often distributed across large regions. The mechanisms that control the integration of these diverse inputs into specific gene expression patterns are still poorly understood. New approaches enabling the dissection of these mechanisms in vivo are needed. RESULTS: Here, we describe TRACER (http://tracerdatabase.embl.de), a resource that centralizes information from a large on-going functional exploration of the mouse genome with different transposon-associated regulatory sensors. Hundreds of insertions have been mapped to specific genomic positions, and their corresponding regulatory potential has been documented by analysis of the expression of the reporter sensor gene in mouse embryos. The data can be easily accessed and provides information on the regulatory activities present in a large number of genomic regions, notably in gene-poor intervals that have been associated with human diseases. CONCLUSIONS: TRACER data enables comparisons with the expression pattern of neighbouring genes, activity of surrounding regulatory elements or with other genomic features, revealing the underlying regulatory architecture of these loci. TRACER mouse lines can also be requested for in vivo transposition and chromosomal engineering, to analyse further regions of interest.

A single enhancer regulating the differential expression of duplicated red-sensitive opsin genes in zebrafish.

A fundamental step in the evolution of the visual system is the gene duplication of visual opsins and differentiation between the duplicates in absorption spectra and expression pattern in the retina. However, our understanding of the mechanism of expression differentiation is far behind that of spectral tuning of opsins. Zebrafish (Danio rerio) have two red-sensitive cone opsin genes, LWS-1 and LWS-2. These genes are arrayed in a tail-to-head manner, in this order, and are both expressed in the long member of double cones (LDCs) in the retina. Expression of the longer-wave sensitive LWS-1 occurs later in development and is thus confined to the peripheral, especially ventral-nasal region of the adult retina, whereas expression of LWS-2 occurs earlier and is confined to the central region of the adult retina, shifted slightly to the dorsal-temporal region. In this study, we employed a transgenic reporter assay using fluorescent proteins and P1-artificial chromosome (PAC) clones encompassing the two genes and identified a 0.6-kb "LWS-activating region" (LAR) upstream of LWS-1, which regulates expression of both genes. Under the 2.6-kb flanking upstream region containing the LAR, the expression pattern of LWS-1 was recapitulated by the fluorescent reporter. On the other hand, when LAR was directly conjugated to the LWS-2 upstream region, the reporter was expressed in the LDCs but also across the entire outer nuclear layer. Deletion of LAR from the PAC clones drastically lowered the reporter expression of the two genes. These results suggest that LAR regulates both LWS-1 and LWS-2 by enhancing their expression and that interaction of LAR with the promoters is competitive between the two genes in a developmentally restricted manner. Sharing a regulatory region between duplicated genes could be a general way to facilitate the expression differentiation in duplicated visual opsins.

Controlling gene activation by enhancers through a drug-inducible topological insulator.

While regulation of gene-enhancer interaction is intensively studied, its application remains limited. Here, we reconstituted arrays of CTCF-binding sites and devised a synthetic topological insulator with tetO for chromatin-engineering (STITCH). By coupling STITCH with tetR linked to the KRAB domain to induce heterochromatin and disable the insulation, we developed a drug-inducible system to control gene activation by enhancers. In human induced pluripotent stem cells, STITCH inserted between MYC and the enhancer down-regulated MYC. Progressive mutagenesis of STITCH led to a preferential escalation of the gene-enhancer interaction, corroborating the strong insulation ability of STITCH. STITCH also altered epigenetic states around MYC. Time-course analysis by drug induction uncovered deposition and removal of H3K27me3 repressive marks follows and reflects, but does not precede and determine, the expression change. Finally, STITCH inserted near NEUROG2 impaired the gene activation in differentiating neural progenitor cells. Thus, STITCH should be broadly useful for functional genetic studies.

Rescue from oculocutaneous albinism type 4 using medaka slc45a2 cDNA driven by its own promoter.

Patients and vertebrate mutants with oculocutaneous albinism type 4 (OCA4) have mutations in the solute carrier family 45 member 2 (slc45a2) gene. However, there is no empirical evidence for this gene-phenotype relationship. There is a unique OCA4 mutant in medaka (b) that exhibits albinism only in the skin, but the mechanism underlying this phenotype is also unknown. In this study, we rescued medaka OCA4 phenotypes, in both the eyes and the skin, by micro-injection of an slc45a2-containing genomic fragment or slc45a2 cDNA driven by its own 0.9-kb promoter. We also identified a spontaneous nucleotide change of 339 bp in the promoter as the b mutation. There are multiple transcription start sites in medaka slc45a2, as in its human ortholog, and only the shortest and eye-specific mRNA is transcribed with the b mutation. Interestingly, we further revealed a conserved pyrimidine (Py)-rich sequence of approximately 10 bp in the promoter by medaka-pufferfish comparative genomics and verified that it plays an indispensable role for expression of slc45a2 in the skin. Further studies of the 0.9-kb promoter identified in this study should provide insights into the cis/trans-regulatory mechanisms underlying the ocular and cutaneous expression of slc45a2.

Control of directionality of chromatin folding for the inter- and intra-domain contacts at the Tfap2c-Bmp7 locus.

BACKGROUND: Contact domains of chromatin serve as a fundamental unit to regulate action of enhancers for target genes. Looping between a pair of CCCTC-binding factor (CTCF)-binding sites in convergent orientations underlies the formation of contact domains, while those in divergent orientations establish domain boundaries. However, every CTCF site is not necessarily engaged in loop or boundary structures, leaving functions of CTCF in varied genomic contexts still elusive. The locus containing Tfap2c and Bmp7 encompasses two contact domains separated by a region between the two genes, termed transition zone (TZ), characterized by two arrays of CTCF sites in divergent configuration. In this study, we created deletion and inversion alleles of these and other regions across the locus and investigated how they impinge on the conformation. RESULTS: Deletion of the whole two CTCF arrays with the CRISPR/Cas9 system resulted in impairment of blocking of chromatin contacts by the TZ, as assessed by the circular chromatin conformation capture assay (4C-seq). Deletion and inversion of either of the two arrays similarly, but less pronouncedly, led to reduction in the blocking activity. Thus, the divergent configuration provides the TZ with the strong boundary activity. Uniquely, we show the TZ harbors a 50-kb region within one of the two arrays that contacts broadly with the both flanking intervals, regardless of the presence or orientation of the other CTCF array. Further, we show the boundary CTCF array has little impact on intra-domain folding; instead, locally associating CTCF sites greatly affect it. CONCLUSIONS: Our results show that the TZ not only separates the two domains, but also bears a wide interval that shows isotropic behavior of chromatin folding, indicating a potentially complex nature of actual boundaries in the genome. We also show that CTCF-binding sites inside a domain greatly contribute to the intra-domain folding of chromatin. Thus, the study reveals diverse and context-dependent roles of CTCF in organizing chromatin conformation at different levels.

Long-term effects of low calcium dialysates on the serum calcium levels during maintenance hemodialysis treatments: A systematic review and meta-analysis.

Hypercalcemia and hyperparathyroidism in patients receiving maintenance hemodialysis (MHD) can cause the progression of cardiovascular diseases (CVD) and mineral bone disorders (MBD). The KDIGO recommends the dialysates with a calcium (Ca) concentration of 1.25-1.5 mmol/L for MHD treatments, but the optimal concentration remains controversial. Here, we conducted a systematic review and a meta-analysis of seven randomized controlled trials examining a total of 622 patients to investigate the optimal concentration for MHD for 6 months or longer. The dialysates with a low Ca concentration (1.125 or 1.25 mmol/L) significantly lowered the serum Ca and raised the intact parathyroid hormone levels by 0.52 mg/dL (95% confidence interval, 0.20-0.85) and 39.59 pg/mL (14.80-64.38), respectively, compared with a high Ca concentration (1.50 or 1.75 mmol/L). Three studies showed that a low concentration was preferred for lowering arterial calcifications or atherosclerosis in different arteries, but one study showed that coronary arterial calcifications increased with a low concentration. Two studies showed contradictory outcomes in terms of MBD. Our meta-analysis showed that a dialysate with a low Ca concentration lowered the serum Ca levels in patients receiving long-term MHD, but further studies are needed to determine the optimal Ca concentration in terms of CVD and MBD.

Roles and regulation of bone morphogenetic protein-7 in kidney development and diseases.

The gene encoding bone morphogenetic protein-7 (Bmp7) is expressed in the developing kidney in embryos and also in the mature organ in adults. During kidney development, expression of Bmp7 is essential to determine the final number of nephrons in and proper size of the organ. The secreted BMP7 acts on the nephron progenitor cells to exert its dual functions: To maintain and expand the progenitor population and to provide them with competence to respond to differentiation cues, each relying on distinct signaling pathways. Intriguingly, in the adult organ, BMP7 has been implicated in protection against and regeneration from injury. Exogenous administration of recombinant BMP7 to animal models of kidney diseases has shown promising effects in counteracting inflammation, apoptosis and fibrosis evoked upon injury. Although the expression pattern of Bmp7 has been well described, the mechanisms by which it is regulated have remained elusive and the processes by which the secretion sites of BMP7 impinge upon its functions in kidney development and diseases have not yet been assessed. Understanding the regulatory mechanisms will pave the way towards gaining better insight into the roles of BMP7, and to achieving desired control of the gene expression as a therapeutic strategy for kidney diseases.

Evolutionarily conserved and divergent regulatory sequences in the fish rod opsin promoter.

Fish have multiple types and subtypes of opsin genes that are expressed in a highly regulated manner in retinal photoreceptor cells. In the rod opsin proximal promoter region (RPPR) of zebrafish (Danio rerio), the BAT 1 regulatory region contains highly conserved OTX (GATTA) and OTX-like (TATTA) sequences that can be recognized by the mammalian cone-rod homeobox (CRX) protein. However, binding of zebrafish crx to the OTX sequence has remained elusive. In contrast to the BAT 1 region, the Ret 1 region, located approximately 20 bp upstream of the BAT 1 region in mammals, is not conserved in zebrafish. In the Ret 1 region, even the core OTX-like sequence (AATTA sequence in mammals) is destructed. We show in this study that a region between Ret 1 and BAT 1 (denoted IRB, Inter-Ret 1-BAT 1) is highly conserved among fish species. Using electrophoretic mobility shift assay (EMSA), we show that zebrafish crx binds to the conserved OTX sequence and that the fish-specific IRB region specifically binds elements present in both retinal and brain nuclear extracts of zebrafish. These results imply that the regulatory mechanisms of opsin gene expression consist not only of evolutionarily conserved but also of divergent machinery among different animal taxa.

Adult stem-like cells in kidney.

Human pluripotent cells are promising for treatment for kidney diseases, but the protocols for derivation of kidney cell types are still controversial. Kidney tissue regeneration is well confirmed in several lower vertebrates such as fish, and the repair of nephrons after tubular damages is commonly observed after renal injury. Even in adult mammal kidney, renal progenitor cell or system is reportedly presents suggesting that adult stem-like cells in kidney can be practical clinical targets for kidney diseases. However, it is still unclear if kidney stem cells or stem-like cells exist or not. In general, stemness is defined by several factors such as self-renewal capacity, multi-lineage potency and characteristic gene expression profiles. The definite use of stemness may be obstacle to understand kidney regeneration, and here we describe the recent broad findings of kidney regeneration and the cells that contribute regeneration.

Characterization of transgenic zebrafish lines that express GFP in the retina, pineal gland, olfactory bulb, hatching gland, and optic tectum.

Transgenic animals are powerful tools to study gene function invivo. Here we characterize several transgenic zebrafish lines that express green fluorescent protein (GFP) under the control of the LCR(RH2)-RH2-1 or LCR(RH2)-RH2-2 green opsin regulatory elements. Using confocal immunomicroscopy, stereo-fluorescence microscopy, and Western blotting, we show that the Tg(LCR(RH2)-RH2-1:GFP)(pt112) and Tg(LCR(RH2)-RH2-2:GFP)(pt115) transgenic zebrafish lines express GFP in the pineal gland and certain types of photoreceptors. In addition, some of these lines also express GFP in the hatching gland, optic tectum, or olfactory bulb. Some of the expression patterns differ significantly from previously published similar transgenic fish lines, making them useful tools for studying the development of the corresponding tissues and organs. In addition, the variations of GFP expression among different lines corroborate the notion that transgenic expression is often subjected to position effect, thus emphasizing the need for careful verification of expression patterns when transgenic animal models are utilized for research.

Bipolar cell-photoreceptor connectivity in the zebrafish (Danio rerio) retina.

Bipolar cells convey luminance, spatial, and color information from photoreceptors to amacrine and ganglion cells. We studied the photoreceptor connectivity of 321 bipolar cells in the adult zebrafish retina. 1,1'-Dioctadecyl-3,3,3',3'-tetramethylindocarbocyanine perchlorate (DiI) was inserted into whole-mounted transgenic zebrafish retinas to label bipolar cells. The photoreceptors that connect to these DiI-labeled cells were identified by transgenic fluorescence or their positions relative to the fluorescent cones, as cones are arranged in a highly ordered mosaic: rows of alternating blue- (B) and ultraviolet-sensitive (UV) single cones alternate with rows of red-(R) and green-sensitive (G) double cones. Rod terminals intersperse among cone terminals. As many as 18 connectivity subtypes were observed, 9 of which-G, GBUV, RG, RGB, RGBUV, RGRod, RGBRod, RGBUVRod, and RRod bipolar cells-accounted for 96% of the population. Based on their axon terminal stratification, these bipolar cells could be further subdivided into ON, OFF, and ON-OFF cells. The dendritic spread size, soma depth and size, and photoreceptor connections of the 308 bipolar cells within the nine common connectivity subtypes were determined, and their dendritic tree morphologies and axonal stratification patterns compared. We found that bipolar cells with the same axonal stratification patterns could have heterogeneous photoreceptor connectivity whereas bipolar cells with the same dendritic tree morphology usually had the same photoreceptor connectivity, although their axons might stratify on different levels.

Identification of a locus control region for quadruplicated green-sensitive opsin genes in zebrafish.

Duplication of opsin genes has a crucial role in the evolution of visual system. Zebrafish have four green-sensitive (RH2) opsin genes (RH2-1, RH2-2, RH2-3, and RH2-4) arrayed in tandem. They are expressed in the short member of the double cones (SDC) but differ in expression areas in the retina and absorption spectra of their encoding photopigments. The shortest and the second shortest wavelength subtypes, RH2-1 and RH2-2, are expressed in the central-to-dorsal retina. The longer wavelength subtype, RH2-3, is expressed circumscribing the RH2-1/RH2-2 area, and the longest subtype, RH2-4, is expressed further circumscribing the RH2-3 area and mainly occupying the ventral retina. The present report shows that a 0.5-kb region located 15 kb upstream of the RH2 gene array is an essential regulator for their expression. When the 0.5-kb region was deleted from a P1-artificial chromosome (PAC) clone encompassing the four RH2 genes and when one of these genes was replaced with a reporter GFP gene, the GFP expression in SDCs was abolished in the zebrafish to which a series of the modified PAC clones were introduced. Transgenic studies also showed that the 0.5-kb region conferred the SDC-specific expression for promoters of a non-SDC (UV opsin) and a nonretinal (keratin 8) gene. Changing the location of the 0.5-kb region in the PAC clone conferred the highest expression for its proximal gene. The 0.5-kb region was thus designated as RH2-LCR analogous to the locus control region of the L-M opsin genes of primates.