Search

Figure S3. Illustration of modeling and correcting for dependence of each marker through the expression level of housekeeping genes.

Figure S5. Knockdown of MAPK7, CDC7, CDC2L5 and CSNKID induces Mesenchymal-Epithelial Transition (MET).

Supplementary Table S1: Patients and clinical information. Supplementary Table S2: Genes used to determine the AR signature.

Figure S2. Scatter-plots showing a linear relationship between expression of housekeeping genes.

Figure S1. Summary of kinome screen.

Figure S4. Heat map analysis of the kinome screen generated by the Wilcoxon method.

Supplementary Figure S1: RNA sequencing and pathway analysis of CTCs from Bone(+) versus Visceral(+) breast cancer patients. Supplementary Figure S2: Correlation between AR expression and AR signatures. Supplementary Figure S3: ER and PR expression in CTCs from Bone(+) versus Visceral(+) breast cancer patients. Supplementary Figure S4: AR antibody validation. Supplementary Figure S5: AR copy number assessment in breast cancer bone metastasis. Supplementary Figure S6: ER and PR expression in primary and metastatic breast cancer. Supplementary Figure S7: AR expression correlates with bone metastasis independently of ER status. Supplementary Figure S8: Enzalutamide treatment induces downregulation of a subset of AR-dependent genes in breast cancer cells. Supplementary Figure S9: AR shRNA validation. Supplementary Figure S10: AR inhibition does not alter primary tumor growth.

Figure S6. Knockdown of MAPK7 does not increase proliferation rate in primary tumor or lung metastasis.

Supplemental Methods, Tables S1-2, Table Legends and Figure Legends Table S1.List of qRT-PCR primers used for epithelial, mesenchymal, inducer and housekeeping genes. Table S2. RNA-ISH quantitation of expression scores used in Figure 5.

Purpose: Fulvestrant is known to promote the degradation of the estrogen receptor (ER) in the nucleus. However, fulvestrant also promotes the aggregation of the newly synthesized ER in the cytoplasm. Accumulation of protein aggregates leads to cell death but this effect is limited as a result of their elimination by the proteasome. We tested whether combining fulvestrant with the proteasome inhibitor, bortezomib, could enhance the accumulation of ER aggregates and cause apoptotic cell death.Experimental Design: The rate of aggregation of the ER was monitored in ER+ breast cancer cells lines, T47D, ZR-75.1, BT474, MDA-MB-361, MCF-7, fulvestrant resistance MCF-7, and tamoxifen-resistant T47D-cyclin D1 cells. Activation of the unfolded protein response, apoptosis, and metabolic rate were also monitored in these cell lines following treatment with fulvestrant, bortezomib, or bortezomib in combination with fulvestrant.Results: We found that bortezomib enhances the fulvestrant-mediated aggregation of the ER in the cytoplasm without blocking the degradation of the ER in the nucleus. Further, these aggregates activate a sustained unfolded protein response leading to apoptotic cell death. Further, we show that the combination induced tumor regression in a breast cancer mouse model of tamoxifen resistance.Conclusions: Adding bortezomib to fulvestrant enhances its efficacy by taking advantage of the unique ability of fulvestrant to promote cytoplasmic aggregates of the ER. As this effect of fulvestrant is independent of the transcriptional activity of the ER, these results suggest that this novel combination may be effective in breast cancers that are ER+ but estrogen independent. Clin Cancer Res; 17(8); 2292–300. ©2011 AACR.

Supplementary Figure S1.

Supplementary Figure 1 from Identification of Cyclin D1– and Estrogen-Regulated Genes Contributing to Breast Carcinogenesis and Progression

Supplementary Figure 8 from Genomic Alterations of Anaplastic Lymphoma Kinase May Sensitize Tumors to Anaplastic Lymphoma Kinase Inhibitors

Supplementary Figure 1 from Elevated CRAF as a Potential Mechanism of Acquired Resistance to BRAF Inhibition in Melanoma

Activating BRAF kinase mutations arise in ∼7% of all human tumors, and preclinical studies have validated the RAF–mitogen-activated protein/extracellular signal-regulated kinase (ERK) kinase–ERK signaling cascade as a potentially important therapeutic target in this setting. Selective RAF kinase inhibitors are currently undergoing clinical development, and based on the experience with other kinase-targeted therapeutics, it is expected that clinical responses to these agents, if observed, will lead to the eventual emergence of drug resistance in most cases. Thus, it is important to establish molecular mechanisms underlying such resistance to develop effective therapeutic strategies to overcome or prevent drug resistance. To anticipate potential mechanisms of acquired resistance to RAF inhibitors during the course of treatment, we established drug-resistant clones from a human melanoma-derived cell line harboring the recurrent V600E activating BRAF mutation, which exhibits exquisite sensitivity to AZ628, a selective RAF kinase inhibitor. We determined that elevated CRAF protein levels account for the acquisition of resistance to AZ628 in these cells, associated with a switch from BRAF to CRAF dependency in tumor cells. We also found that elevated CRAF protein levels may similarly contribute to primary insensitivity to RAF inhibition in a subset of BRAF mutant tumor cells. Interestingly, AZ628-resistant cells demonstrating either primary drug insensitivity or acquired drug resistance exhibit exquisite sensitivity to the HSP90 inhibitor geldanamycin. Geldanamycin effectively promotes the degradation of CRAF, thereby revealing a potential therapeutic strategy to overcome resistance to RAF inhibition in a subset of BRAF mutant tumors. [Cancer Res 2008;68(12):4853–61]

Supplementary Tables 1-2, Figure Legends 1-2 from Identification of Cyclin D1– and Estrogen-Regulated Genes Contributing to Breast Carcinogenesis and Progression

Supplementary Figure Legends 1-3, Table 1 from Elevated CRAF as a Potential Mechanism of Acquired Resistance to BRAF Inhibition in Melanoma

Supplementary Figure 2 from Elevated CRAF as a Potential Mechanism of Acquired Resistance to BRAF Inhibition in Melanoma

Supplementary Table 1 from Elevated CRAF as a Potential Mechanism of Acquired Resistance to BRAF Inhibition in Melanoma

Supplementary Figure 4 from Genomic Alterations of Anaplastic Lymphoma Kinase May Sensitize Tumors to Anaplastic Lymphoma Kinase Inhibitors

Tumors can become lethal when they progress from preinvasive lesions to invasive carcinomas. Here, we identify candidate tumor progression genes using gene array analysis of preinvasive and invasive tumors from mice, which were then evaluated in human cancers. Immediate early response protein IEX-1, small stress protein 1 (HSPB8), and tumor necrosis factor-associated factor–interacting protein mRNAs displayed higher expression levels in invasive lesions than in preinvasive lesions using samples obtained by laser capture microdissection (LCM) from transgenic erbB2, ras, and cyclin D1 mice. LCM-isolated tissues from patient-matched normal, ductal carcinoma in situ, and invasive ductal carcinoma revealed similar increased expression in invasive human cancers compared with preinvasive and normal samples. These genes induced anchorage independence, increased cell proliferation, and protected against apoptosis, singly or in collaboration with erbB2. Surprisingly, they were all up-regulated by 17β-estradiol and cyclin D1, and cyclin D1 overexpression increased p300/CBP binding to their promoters, supporting the model that cyclin D1-estrogen receptor (ER) coactivator interactions may be important to its role in ER-positive breast cancer. Additionally, an irreversible dual kinase inhibitor of ErbB signaling inhibited expression of the same genes. The up-regulation of genes contributing to increased invasiveness of ER-positive cancers offers a novel explanation for the contribution of cyclin D1 to a worse prognosis in ER-positive cancers. As targets of estrogen, cyclin D1, and erbB2 signaling, these candidates offer insights into the nature of the second events involved in breast cancer progression, regulatory events contributing to invasion, and potential targets of combined inhibition of hormone and growth factor signaling pathways. (Cancer Res 2006; 66(24): 11649-58)

Supplementary Legends for Table 1 and Figures 1-9 from Genomic Alterations of Anaplastic Lymphoma Kinase May Sensitize Tumors to Anaplastic Lymphoma Kinase Inhibitors

Supplementary Figure 3 from Genomic Alterations of Anaplastic Lymphoma Kinase May Sensitize Tumors to Anaplastic Lymphoma Kinase Inhibitors

Selective kinase inhibitors have had a substantial impact on the field of medical oncology. Whereas these agents can elicit dramatic clinical responses in some settings, their activity is generally limited to a subset of treated patients whose tumor cells harbor a specific genetic lesion. We have established an automated platform for examining the sensitivity to various molecularly targeted inhibitors across a large panel of human tumor-derived cell lines to identify additional genotype-correlated responses that may be clinically relevant. Among the inhibitors tested in a panel of 602 cell lines derived from a variety of human cancers, we found that a selective inhibitor of the anaplastic lymphoma kinase (ALK) potently suppressed growth of a small subset of tumor cells. This subset included lines derived from anaplastic large cell lymphomas, non–small-cell lung cancers, and neuroblastomas. ALK is a receptor tyrosine kinase that was first identified as part of a protein fusion derived from a chromosomal translocation detected in the majority of anaplastic large cell lymphoma patients, and has recently been implicated as an oncogene in a small fraction of non–small-cell lung cancers and neuroblastomas. Significantly, sensitivity in these cell lines was well correlated with specific ALK genomic rearrangements, including chromosomal translocations and gene amplification. Moreover, in such cell lines, ALK kinase inhibition can lead to potent suppression of downstream survival signaling and an apoptotic response. These findings suggest that a subset of lung cancers, lymphomas, and neuroblastomas that harbor genomic ALK alterations may be clinically responsive to pharmacologic ALK inhibition. [Cancer Res 2008;68(9):3389–95]

Supplementary Figure 7 from Genomic Alterations of Anaplastic Lymphoma Kinase May Sensitize Tumors to Anaplastic Lymphoma Kinase Inhibitors

Supplementary Figure 3 from Elevated CRAF as a Potential Mechanism of Acquired Resistance to BRAF Inhibition in Melanoma

Supplementary Figure 2 from Identification of Cyclin D1– and Estrogen-Regulated Genes Contributing to Breast Carcinogenesis and Progression

Supplementary Figure 1 from Genomic Alterations of Anaplastic Lymphoma Kinase May Sensitize Tumors to Anaplastic Lymphoma Kinase Inhibitors

Supplementary Figure 9 from Genomic Alterations of Anaplastic Lymphoma Kinase May Sensitize Tumors to Anaplastic Lymphoma Kinase Inhibitors

Supplementary Table 1 from Genomic Alterations of Anaplastic Lymphoma Kinase May Sensitize Tumors to Anaplastic Lymphoma Kinase Inhibitors

Supplementary Figure 2 from Genomic Alterations of Anaplastic Lymphoma Kinase May Sensitize Tumors to Anaplastic Lymphoma Kinase Inhibitors

Supplementary Figure 5 from Genomic Alterations of Anaplastic Lymphoma Kinase May Sensitize Tumors to Anaplastic Lymphoma Kinase Inhibitors

Supplementary Figure 6 from Genomic Alterations of Anaplastic Lymphoma Kinase May Sensitize Tumors to Anaplastic Lymphoma Kinase Inhibitors

In vitro expansion of human primordial germ cell-like cells (hPGCLCs), a pluripotent stem cell-derived PGC model, has proved challenging due to rapid loss of primordial germ cell (PGC)-like identity and limited cell survival/proliferation. Here, we describe long-term culture hPGCLCs (LTC-hPGCLCs), which actively proliferate in a serum-free, feeder-free condition without apparent limit as highly homogeneous diploid cell populations maintaining transcriptomic and epigenomic characteristics of hPGCLCs. Histone proteomics confirmed reduced H3K9me2 and increased H3K27me3 marks in LTC-hPGCLCs compared with induced pluripotent stem cells (iPSCs). LTC-hPGCLCs established from multiple human iPSC clones of both sexes were telomerase positive, senescence-free cells readily passaged with minimal cell death or deviation from the PGC-like identity. LTC-hPGCLCs are capable of differentiating to DAZL-positive M-spermatogonia-like cells in the xenogeneic reconstituted testis (xrTestis) organ culture milieu as well as efficiently producing fully pluripotent embryonic germ cell-like cells in the presence of stem cell factor and fibroblast growth factor 2. Thus, LTC-hPGCLCs provide convenient access to unlimited amounts of high-quality and homogeneous hPGCLCs.

An in vitro model of human ovarian follicles would greatly benefit the study of female reproduction. Ovarian development requires the combination of germ cells and several types of somatic cells. Among these, granulosa cells play a key role in follicle formation and support for oogenesis. Whereas efficient protocols exist for generating human primordial germ cell-like cells (hPGCLCs) from human induced pluripotent stem cells (hiPSCs), a method of generating granulosa cells has been elusive. Here, we report that simultaneous overexpression of two transcription factors (TFs) can direct the differentiation of hiPSCs to granulosa-like cells. We elucidate the regulatory effects of several granulosa-related TFs and establish that overexpression of NR5A1 and either RUNX1 or RUNX2 is sufficient to generate granulosa-like cells. Our granulosa-like cells have transcriptomes similar to human fetal ovarian cells and recapitulate key ovarian phenotypes including follicle formation and steroidogenesis. When aggregated with hPGCLCs, our cells form ovary-like organoids (ovaroids) and support hPGCLC development from the premigratory to the gonadal stage as measured by induction of DAZL expression. This model system will provide unique opportunities for studying human ovarian biology and may enable the development of therapies for female reproductive health.

Exposure of pregnant F0 mouse dams to the obesogen tributyltin (TBT) predisposes unexposed male descendants to obesity and diverts mesenchymal stem cells (MSCs) toward the adipocytic lineage. TBT promotes adipogenic commitment and differentiation of MSCs in vitro. To identify TBT-induced factors predisposing MSCs toward the adipocytic fate, we exposed mouse MSCs to TBT, the peroxisome proliferator activated receptor gamma (PPARγ)-selective agonist rosiglitazone, or the retinoid X receptor (RXR)-selective agonist LG-100268. Then we determined their transcriptomal profiles to determine candidate microRNAs (miR) regulating adipogenic commitment and differentiation. Of the top 10 candidate microRNAs predicted by Ingenuity Pathway Analysis, miR-21, miR-33, and miR-223 were expressed consistent with an ability to differentially regulate target genes during adipogenesis. We found that 24-hour exposure to 50nM TBT caused miR-223 levels in MSCs to increase; expression of its target genes ZEB1, NFIB, and FOXP1 was decreased. Rosiglitazone and TBT increased miR-223 levels. This induction was inhibited by the PPARγ antagonist T0070907 but not by the RXR antagonists HX531 or UVI3003, placing miR-223 downstream of PPARγ. Chromatin immunoprecipitation confirmed TBT-induced binding of PPARγ to regulatory elements in the miR-223 promoter. miR-223 levels were elevated in white adipose tissue of F2 and F3 male descendants of pregnant F0 mouse dams exposed to 50nM TBT throughout gestation. miR-223 levels were potentiated in males fed an increased fat diet. We infer that TBT induced miR-223 expression and increased adipogenesis in MSCs through the PPARγ pathway and that transgenerationally increased expression of miR-223 plays an important role in the development of obesity caused by TBT exposure.

Exposure of pregnant F0 mouse dams to the obesogen tributyltin (TBT) predisposes unexposed male descendants to obesity and diverts mesenchymal stem cells (MSCs) toward the adipocytic lineage. TBT also promotes adipogenic commitment and differentiation of MSCs, in vitro. We sought to identify TBT-induced factors predisposing MSCs toward the adipocytic fate. We exposed mouse MSCs to TBT, the PPARγ-selective agonist rosiglitazone or the RXR-selective agonist LG-100268 and determined their transcriptomal profiles to determine candidate microRNAs (miR) regulating adipogenic commitment and differentiation. Of the top 10 candidate microRNAs predicted by Ingenuity Pathway Analysis, miR-21, miR-33 and miR-223 were expressed in a manner consistent with an ability to differentially regulate target genes during adipogenesis. After 24 hours exposure to 50 nM TBT, miR-223 levels in MSCs were increased and expression of its target genes ZEB1, NFIB, and FOXP1 was decreased. Both ROSI and TBT increased miR-223 levels, and this induction was inhibited by the PPARγ antagonist T0070907 but not by the RXR antagonists HX531 or UVI3003, placing miR-223 downstream of PPARγ. Chromatin immunoprecipitation confirmed TBT-induced binding of PPARγ to regulatory elements in the miR-223 promoter. miR-223 levels were elevated in white adipose tissue of F2 and F3 male descendants of pregnant F0 mouse dams exposed to 50 nM TBT throughout gestation. miR-223 levels were further induced in males fed with an increased fat diet. We infer that TBT induced miR-223 expression and increased adipogenesis in MSCs through the PPARγ pathway and that transgenerationally increased expression of miR-223 plays an important role in the development of obesity caused by TBT exposure.

Burgeoning evidence demonstrates that responses to environmental exposures can be transmitted to subsequent generations through the germline without DNA mutations1,2. This is controversial because underlying mechanisms remain to be identified. Therefore, understanding how effects of environmental exposures are transmitted to unexposed generations without DNA mutations is a fundamental unanswered question in biology. Here, we used an established murine model of transgenerational obesity to show that direct or ancestral exposure to the obesogen tributyltin (TBT) elicited persistent changes in topologically associating domains (TADs) in primordial germ cells (PGCs) isolated from embryos of exposed and subsequent unexposed generations. New TAD boundaries were formed within theIdegene encoding insulin degrading enzyme in the exposed PGCs, then stably maintained in PGCs of the subsequent (unexposed) two generations. Concomitantly,IdemRNA expression was decreased in livers of male descendants from the exposed dams. These animals were hyperinsulinemic and hyperglycemic, phenocopyingIde-deficient mice that are predisposed to adult-onset obesity. Creation of new TAD boundaries in PGCs, suppression of hepaticIdemRNA, increased fat mass, hyperinsulinemia and hyperglycemia were male-specific. Our results provide a plausible molecular mechanism underlying transmission of the transgenerational predisposition to obesity caused by gestational exposure to an environmental obesogen. They also provide an entry point for future studies aimed at understanding how environmental exposures alter chromatin structure to influence physiology across multiple generations in mammals.

The hominoid-specific endogenous retrovirus LTR5_Hs is transcriptionally activated in human primordial germ cell-like cells (hPGCLCs), a pluripotent stem cell-derived cell culture model of PGCs. Here, taking the unique advantage of our novel cell culture method to obtain large amounts of pure hPGCLCs, we performed proteomics profiling of hPGCLCs and detected various viral proteins produced from the LTR5_Hs RNA via ribosomal frameshifting. We also present transmission electron microscopy images of 100-nm diameter virus-like particles (VLPs) assembled at the surface of hPGCLCs. Compared to hPGCLCs, expression of LTR5_Hs RNA is far weaker in human seminomas, the germ cell tumors resembling PGCs. Re-analysis of published single cell RNA-seq data of human embryos revealed strong activation of LTR5_Hs in migrating PGCs but suppressed in PGCs upon they reach the gonadal anlagen. In the microfluidics-supported polarized embryoids mimicking peri-implantation stages of human embryos, LTR5_Hs RNA was detected by RNA in situ hybridization in NANOG+/TFAP2C+/SOX17+ cells resembling freshly emerged PGCs. These results support that human germ cells produce LTR5_Hs proteins and VLPs during their earliest stages of normal development until their settlement in the gonadal anlagen.SUMMARY STATEMENTThe hominoid-specific endogenous retrovirus LTR5_Hs is activated in a cell culture model resembling early-stage human primordial germ cells, producing not only viral RNA but also retrovirus proteins and virus-like particles.

Germ cells are the vehicle of human reproduction, arising early in embryonic development and developing throughout adult life until menopause onset in women. Primordial germ cells are the common precursors of germline cells in both sexes, undergoing sexual specification into oogonia or gonocytes which further develop into oocytes or spermatocytes during development. Methods for recapitulation of primordial germ cell and oogonia formation have been developed extensively in recent decades, but fundamental technical limitations in their methodologies, throughput, and yield limit their utilization. Recently, transcription factor (TF)-based methods for human primordial germ cell-like cell (hPGCLC) formation, mouse meiotic entry, and mouse oocyte maturation have demonstrated the feasibility of gene overexpression screening in identifying potent regulators of germ cell development. Here we screened 47 folliculogenesis-regulating TFs for their role in hPGCLC and oogonia formation, identifying DLX5, HHEX, and FIGLA whose individual overexpression enhances hPGCLC formation from hiPSCs. Additionally, we identify a set of three TFs, ZNF281, LHX8, and SOHLH1, whose combinatorial overexpression drives direct oogonia-like formation from hiPSCs in a four-day, feeder-free monolayer culture condition with additional feeder-free culture capabilities post-isolation. We characterize these TF-based germ cells via gene and protein expression analyses, and demonstrate their broad similarity to in vivo germ cells. Together, these results identify novel regulators of human germ cell development and establish new TF-based tools for human in vitro oogenesis research.

An in vitro model of human ovarian follicles would greatly benefit the study of female reproduction. Ovarian development requires the combination of germ cells and their supporting somatic cells, known as granulosa cells. Whereas efficient protocols exist for generating human primordial germ cell-like cells (hPGCLCs) from human iPSCs, a method of generating granulosa cells has been elusive. Here we report that simultaneous overexpression of two transcription factors (TFs) can direct the differentiation of human iPSCs to granulosa-like cells. We elucidate the regulatory effects of several granulosa-related TFs, and establish that overexpression of NR5A1 and either RUNX1 or RUNX2 is necessary and sufficient to generate granulosa-like cells. Our granulosa-like cells form ovary-like organoids (ovaroids) when aggregated with hPGCLCs, and recapitulate key ovarian phenotypes including support of germ cell maturation, follicle formation, and steroidogenesis. This model system will provide unique opportunities for studying human ovarian biology, and may enable the development of therapies for female reproductive health.

Differential placental blood flow and nutrient transport can lead to both intrauterine growth restriction (IUGR) and macrosomia. Both conditions can lead to adult obesity and other conditions clustered as metabolic syndrome. We previously showed that pregnant hemi-ovariectomized mice have a crowded uterine horn, resulting in siblings whose birth weights differ by over 100% due to differential blood flow based on uterine position. We used this crowded uterus model to compare IUGR and macrosomic male mice and also identified IUGR males with rapid (IUGR-R) and low (IUGR-L) postweaning weight gain. At week 12 IUGR-R males were heavier than IUGR-L males and did not differ from macrosomic males. Rapid growth in IUGR-R males led to glucose intolerance compared to IUGR-L males and down-regulation of adipocyte signaling pathways for fat digestion and absorption and type II diabetes. Macrosomia led to increased fat mass and altered adipocyte size distribution compared to IUGR males, and down-regulation of signaling pathways for carbohydrate and fat digestion and absorption relative to IUGR-R. Clustering analysis of gonadal fat transcriptomes indicated more similarities than differences between IUGR-R and macrosomic males compared to IUGR-L males. Our findings suggest two pathways to adult metabolic disease: macrosomia and IUGR with rapid postweaning growth rate.

Polycyclic aromatic hydrocarbons like benzo[a]pyrene (BaP) are generated during incomplete combustion of organic materials. Prior research has demonstrated that BaP is a prenatal ovarian toxicant and carcinogen. However, the metabolic pathways active in the embryo and its developing gonads and the mechanisms by which prenatal exposure to BaP predisposes to ovarian tumors later in life remain to be fully elucidated. To address these data gaps, we orally dosed pregnant female mice with BaP from embryonic day (E) 6.5 to E11.5 (0, 0.2, or 2 mg/kg/day) for metabolite measurement or E9.5 to E11.5 (0 or 3.33 mg/kg/day) for embryonic gonad RNA sequencing. Embryos were harvested at E13.5 for both experiments. The sum of BaP metabolite concentrations increased significantly with dose in the embryos and placentas, and concentrations were significantly higher in female than male embryos and in embryos than placentas. RNA sequencing revealed that enzymes involved in metabolic activation of BaP are expressed at moderate to high levels in embryonic gonads and that greater transcriptomic changes occurred in the ovaries in response to BaP than in the testes. We identified 490 differentially expressed genes (DEGs) with false discovery rate P-values < 0.05 when comparing BaP-exposed to control ovaries but no statistically significant DEGs between BaP-exposed and control testes. Genes related to monocyte/macrophage recruitment and activity, prolactin family genes, and several keratin genes were among the most upregulated genes in the BaP-exposed ovaries. Results show that developing ovaries are more sensitive than testes to prenatal BaP exposure, which may be related to higher concentrations of BaP metabolites in female embryos.

Identification of dynamic changes in chromatin conformation is a fundamental task in genetics. In 2020, Galan et al.1 presented CHESS (Comparison of Hi-C Experiments using Structural Similarity), a novel computational algorithm designed for systematic identification of structural differences in chromatin-contact maps. Using CHESS, the same group recently reported that chromatin organization is largely maintained across tissues during dorsoventral patterning of fruit fly embryos despite tissue-specific chromatin states and gene expression2. However, here we show that the primary outputs of CHESS–namely, the structural similarity index (SSIM) profiles–are nearly identical regardless of the input matrices, even when query and reference reads were shuffled to destroy any significant differences. This issue stems from the dominance of the regional counting noise arising from stochastic sampling in chromatin-contact maps, reflecting a fundamentally incorrect assumption of the CHESS algorithm. Therefore, biological interpretation of SSIM profiles generated by CHESS requires considerable caution.

Primordial germ cells (PGCs) are common ancestors of all germline cells. In mammals, PGCs emerge in early-stage embryos around the timing of gastrulation at or near epiblast, and specification of PGCs from their precursor cells involves multiple growth factors secreted by adjacent cells. Recent advancements in germline stem cell biology have made it possible to generate PGC-like cell culture models (PGCLCs for PGC-like cells) from human and mouse pluripotent stem cells by mimicking the embryonic growth factor environment in vitro. Here we describe a method of producing human PGCLCs from primed-pluripotency induced pluripotent stem cells (iPSCs) via temporal conversion to naive pluripotency followed by formation of embryoid bodies (EBs) using the spin-EB method.

Primordial germ cells (PGCs) are common ancestors of all germline cells. In mammals, PGCs emerge in early-stage embryos around the timing of gastrulation at or near epiblast, and specification of PGCs from their precursor cells involves multiple growth factors secreted by adjacent cells. Recent advancements in germline stem cell biology have made it possible to generate PGC-like cell culture models (PGCLCs for PGC-like cells) from human and mouse pluripotent stem cells by mimicking the embryonic growth factor environment in vitro. Here we describe a method of producing human PGCLCs from primed-pluripotency induced pluripotent stem cells (iPSCs) via temporal conversion to naive pluripotency followed by formation of embryoid bodies (EBs) using the spin-EB method.

Whereas in ovo exposure of genetically male (ZZ) chicken embryos to exogenous estrogens temporarily feminizes gonads at the time of hatching, the morphologically ovarian ZZ-gonads (FemZZs for feminized ZZ gonads) are masculinized back to testes within 1 year. To identify the feminization-resistant “memory” of genetic male sex, FemZZs showing varying degrees of feminization were subjected to transcriptomic, DNA methylome, and immunofluorescence analyses. Protein-coding genes were classified based on their relative mRNA expression across normal ZZ-testes, genetically female (ZW) ovaries, and FemZZs. We identified a group of 25 genes that were strongly expressed in both ZZ-testes and FemZZs but dramatically suppressed in ZW-ovaries. Interestingly, 84% (21/25) of these feminization-resistant testicular marker genes, including the DMRT1 master masculinizing gene, were located in chromosome Z. Expression of representative marker genes of germline cells (eg, DAZL or DDX4/VASA) was stronger in FemZZs than normal ZZ-testes or ZW-ovaries. We also identified 231 repetitive sequences (RSs) that were strongly expressed in both ZZ-testes and FemZZs, but these RSs were not enriched in chromosome Z. Although 94% (165/176) of RSs exclusively expressed in ZW-ovaries were located in chromosome W, no feminization-inducible RS was detected in FemZZs. DNA methylome analysis distinguished FemZZs from normal ZZ- and ZW-gonads. Immunofluorescence analysis of FemZZ gonads revealed expression of DMRT1 protein in medullary SOX9+ somatic cells and apparent germline cell populations in both medulla and cortex. Taken together, our study provides evidence that both somatic and germline cell populations in morphologically feminized FemZZs maintain significant transcriptomic and epigenetic memories of genetic sex.

Human embryonic stem cells (hESCs) can be captured in a primed state in which they resemble the postimplantation epiblast, or in a naive state where they resemble the preimplantation epiblast. Naive-cell-specific culture conditions allow the study of preimplantation development ex vivo but reportedly lead to chromosomal abnormalities, which compromises their utility in research and potential therapeutic applications. Although MEK inhibition is essential for the naive state, here we show that reduced MEK inhibition facilitated the establishment and maintenance of naive hESCs that retained naive-cell-specific features, including global DNA hypomethylation, HERVK expression, and two active X chromosomes. We further show that hESCs cultured under these modified conditions proliferated more rapidly; accrued fewer chromosomal abnormalities; and displayed changes in the phosphorylation levels of MAPK components, regulators of DNA damage/repair, and cell cycle. We thus provide a simple modification to current methods that can enable robust growth and reduced genomic instability in naive hESCs.