Bridging sequence diversity and tissue-specific expression by DNA methylation in genes of the mouse prolactin superfamily

A diversity outbred F1 mouse model identifies host-intrinsic genetic regulators of response to immune checkpoint inhibitors

Immune checkpoint inhibitors (ICI) have improved outcomes for a variety of malignancies; however, many patients fail to benefit. While tumor-intrinsic mechanisms are likely involved in therapy resistance, it is unclear to what extent host genetic background influences response. To investigate this, we utilized the Diversity Outbred (DO) and Collaborative Cross (CC) mouse models. DO mice are an outbred stock generated by crossbreeding eight inbred founder strains, and CC mice are recombinant inbred mice generated from the same eight founders. We generated 207 DOB6F1 mice representing 48 DO dams and demonstrated that these mice reliably accept the C57BL/6-syngeneic B16F0 tumor and that host genetic background influences response to ICI. Genetic linkage analysis from 142 mice identified multiple regions including one within chromosome 13 that associated with therapeutic response. We utilized 6 CC strains bearing the positive (NZO) or negative (C57BL/6) driver genotype in this locus. We found that 2/3 of predicted responder CCB6F1 crosses show reproducible ICI response. The chromosome 13 locus contains the murine prolactin family, which is a known immunomodulating cytokine associated with various autoimmune disorders. To directly test whether prolactin influences ICI response rates, we implanted inbred C57BL/6 mice with subcutaneous slow-release prolactin pellets to induce mild hyperprolactinemia. Prolactin augmented ICI response against B16F0, with increased CD8 infiltration and 5/8 mice exhibiting slowed tumor growth relative to controls. This study highlights the role of host genetics in ICI response and supports the use of F1 crosses in the DO and CC mouse populations as powerful cancer immunotherapy models.

Blastocyst transfer in mice alters the placental transcriptome and growth

Assisted reproduction technologies (ART) are becoming increasingly common. Therefore, how these procedures influence gene regulation and feto-placental development are important to explore. Here, we assess the effects of blastocyst transfer on mouse placental growth and transcriptome. C57Bl/6 blastocysts were transferred into uteri of B6D2F1 pseudopregnant females and dissected at embryonic day 10.5 for analysis. Compared to non-transferred controls, placentas from transferred conceptuses weighed less even though the embryos were larger on average. This suggested a compensatory increase in placental efficiency. RNA-sequencing of whole male placentas revealed 543 differentially expressed genes (DEGs) after blastocyst transfer: 188 and 355 genes were down-regulated and up-regulated, respectively. DEGs were independently validated in male and female placentas. Bioinformatic analyses revealed that DEGs represented expression in all major placental cell types and included genes that are critical for placenta development and/or function. Furthermore, the direction of transcriptional change in response to blastocyst transfer implied an adaptive response to improve placental function to maintain fetal growth. Our analysis revealed that CpG methylation at regulatory regions of two DEGs was unchanged in female transferred placentas and that DEGs had fewer gene-associated CpG islands (within ~20 kb region) compared to the larger genome. These data suggested that altered methylation at proximal promoter regions might not lead to transcriptional disruption in transferred placentas. Genomic clustering of some DEGs warrants further investigation of long-range, cis-acting epigenetic mechanisms including histone modifications together with DNA methylation. We conclude that embryo transfer, a protocol required for ART, significantly impacts the placental transcriptome and growth.

Origin of a rapidly evolving homeostatic control system programming testis function

Mammals share common strategies for regulating reproduction, including a conserved hypothalamic-pituitary-gonadal axis; yet, individual species exhibit differences in reproductive performance. In this report, we describe the discovery of a species-restricted homeostatic control system programming testis growth and function. Prl3c1 is a member of the prolactin gene family and its protein product (PLP-J) was discovered as a uterine cytokine contributing to the establishment of pregnancy. We utilized mouse mutagenesis of Prl3c1 and revealed its involvement in the regulation of the male reproductive axis. The Prl3c1-null male reproductive phenotype was characterized by testiculomegaly and hyperandrogenism. The larger testes in the Prl3c1-null mice were associated with an expansion of the Leydig cell compartment. Prl3c1 locus is a template for two transcripts (Prl3c1-v1 and Prl3c1-v2) expressed in a tissue-specific pattern. Prl3c1-v1 is expressed in uterine decidua, while Prl3c1-v2 is expressed in Leydig cells of the testis. 5'RACE, chromatin immunoprecipitation and DNA methylation analyses were used to define cell-specific promoter usage and alternative transcript expression. We examined the Prl3c1 locus in five murid rodents and showed that the testicular transcript and encoded protein are the result of a recent retrotransposition event at the Mus musculus Prl3c1 locus. Prl3c1-v1 encodes PLP-J V1 and Prl3c1-v2 encodes PLP-J V2. Each protein exhibits distinct intracellular targeting and actions. PLP-J V2 possesses Leydig cell-static actions consistent with the Prl3c1-null testicular phenotype. Analysis of the biology of the Prl3c1 gene has provided insight into a previously unappreciated homeostatic setpoint control system programming testicular growth and function.

Novel O-GlcNAcylation on Ser(40) of canonical H2A isoforms specific to viviparity

We report here newly discovered O-linked-N-acetylglucosamine (O-GlcNAc) modification of histone H2A at Ser⁴⁰ (H2AS40Gc). The mouse genome contains 18 H2A isoforms, of which 13 have Ser⁴⁰ and the other five have Ala⁴⁰. The combination of production of monoclonal antibody and mass spectrometric analyses with reverse-phase (RP)-high performance liquid chromatography (HPLC) fractionation indicated that the O-GlcNAcylation is specific to the Ser⁴⁰ isoforms. The H2AS40Gc site is in the L1 loop structure where two H2A molecules interact in the nucleosome. Targets of H2AS40Gc are distributed genome-wide and are dramatically changed during the process of differentiation in mouse trophoblast stem cells. In addition to the mouse, H2AS40Gc was also detected in humans, macaques and cows, whereas non-mammalian species possessing only the Ala⁴⁰ isoforms, such as silkworms, zebrafish and Xenopus showed no signal. Genome database surveys revealed that Ser⁴⁰ isoforms of H2A emerged in Marsupialia and persisted thereafter in mammals. We propose that the emergence of H2A Ser⁴⁰ and its O-GlcNAcylation linked a genetic event to genome-wide epigenetic events that correlate with the evolution of placental animals.

Cell Type-Specific Sexual Dimorphism in Rat Pituitary Gene Expression During Maturation

The most obvious functional differences between mammalian males and females are related to the control of reproductive physiology and include patterns of GnRH and gonadotropin release, the timing of puberty, sexual and social behavior, and the regulation of food intake and body weight. Using the rat as the best-studied mammalian model for maturation, we examined the expression of major anterior pituitary genes in five secretory cell types of developing males and females. Corticotrophs show comparable Pomc profiles in both sexes, with the highest expression occurring during the infantile period. Somatotrophs and lactotrophs also exhibit no difference in Gh1 and Prl profiles during embryonic to juvenile age but show the amplification of Prl expression in females and Gh1 expression in males during peripubertal and postpubertal ages. Gonadotrophs exhibit highly synchronized Lhb, Fshb, Cga, and Gnrhr expression in both sexes, but the peak of expression occurs during the infantile period in females and at the end of the juvenile period in males. Thyrotrophs also show different developmental Tshb profiles, which are synchronized with the expression of gonadotroph genes in males but not in females. These results indicate the lack of influence of sex on Pomc expression and the presence of two patterns of sexual dimorphism in the expression of other pituitary genes: a time shift in the peak expression during postnatal development, most likely reflecting the perinatal sex-specific brain differentiation, and modulation of the amplitude of expression during late development, which is secondary to the establishment of the hypothalamic-pituitary-gonadal and -thyroid axes.

DNA methylation and histone modification profiles of mouse organic anion transporting polypeptides

Organic anion transporting polypeptides (rodents, Oatps; human, OATPs) are primarily involved in the transmembrane transportation of a wide range of endogenous and exogenous compounds. Multiple mouse Oatp1 isoforms are closely located on chromosome 6, where each isoform shows distinct tissue distribution; Oatp1b2, Oatp1a6, and Oatp1c1 are expressed exclusively in the liver, kidney, and cerebrum, respectively; Oatp1a1 in the liver and kidney; and Oatp1a4 in the liver and cerebrum. We have identified tissue-dependent differentially methylated region (T-DMR) around the transcriptional start site (TSS) of Oatp1b2, which correlates with its liver-specific expression. Bisulfite sequencing also demonstrated the presence of T-DMRs around the TSS in other Oatp1 genes: CpG dinucleotides at +149 relative to the TSS for Oatp1c1; -48, +101, and +356 for Oatp1a4; -572 and -550 for Oatp1a1; and -122 and +216 for Oatp1a6 were differentially methylated among the liver, kidney, and cerebrum. These methylation profiles were largely consistent with the tissue distribution of Oatp1 mRNAs. Chromatin immunoprecipitation assay revealed that the mRNA expression of Oatp1 genes was accompanied by acetylated histone H3. Human OATP1B1 and OATP1B3 are located on chromosome 12p12 in the OATP1 cluster; both show predominant expression in the liver. These genes also contained T-DMRs that were hypomethylated in the liver, compared with kidney cortex: -511, -411, and +92 relative to the TSS for OATP1B1 and -331, +70, and +73 for OATP1B3. These results suggest that the difference in epigenetic profiles comprising DNA methylation and histone acetylation determines the distinct tissue distribution of Oatp/OATP mRNAs.