MUC1-C Is a Common Driver of Acquired Osimertinib Resistance in NSCLC

INTRODUCTION

Osimertinib is an irreversible EGFR tyrosine kinase inhibitor approved for the first-line treatment of patients with metastatic NSCLC harboring EGFR exon 19 deletions or L858R mutations. Patients treated with osimertinib invariably develop acquired resistance by mechanisms involving additional EGFR mutations, MET amplification, and other pathways. There is no known involvement of the oncogenic MUC1-C protein in acquired osimertinib resistance.

METHODS

H1975/EGFR (L858R/T790M) and patient-derived NSCLC cells with acquired osimertinib resistance were investigated for MUC1-C dependence in studies of EGFR pathway activation, clonogenicity, and self-renewal capacity.

RESULTS

We reveal that MUC1-C is up-regulated in H1975 osimertinib drug-tolerant persister cells and is necessary for activation of the EGFR pathway. H1975 cells selected for stable osimertinib resistance (H1975-OR) and MGH700-2D cells isolated from a patient with acquired osimertinib resistance are found to be dependent on MUC1-C for induction of (1) phospho (p)-EGFR, p-ERK, and p-AKT, (2) EMT, and (3) the resistant phenotype. We report that MUC1-C is also required for p-EGFR, p-ERK, and p-AKT activation and self-renewal capacity in acquired osimertinib-resistant (1) MET-amplified MGH170-1D #2 cells and (2) MGH121 Res#2/EGFR (T790M/C797S) cells. Importantly, targeting MUC1-C in these diverse models reverses osimertinib resistance. In support of these results, high MUC1 mRNA and MUC1-C protein expression is associated with a poor prognosis for patients with EGFR-mutant NSCLCs.

CONCLUSIONS

Our findings reveal that MUC1-C is a common effector of osimertinib resistance and is a potential target for the treatment of osimertinib-resistant NSCLCs.

Modeling embryo-endometrial interface recapitulating human embryo implantation

The initiation of human pregnancy is marked by the implantation of an embryo into the uterine environment; however, the underlying mechanisms remain largely elusive. To address this knowledge gap, we developed hormone-responsive endometrial organoids (EMO), termed apical-out (AO)-EMO, which emulate the in vivo architecture of endometrial tissue. The AO-EMO comprise an exposed apical epithelium surface, dense stromal cells, and a self-formed endothelial network. When cocultured with human embryonic stem cell-derived blastoids, the three-dimensional feto-maternal assembloid system recapitulates critical implantation stages, including apposition, adhesion, and invasion. Endometrial epithelial cells were subsequently disrupted by syncytial cells, which invade and fuse with endometrial stromal cells. We validated this fusion of syncytiotrophoblasts and stromal cells using human blastocysts. Our model provides a foundation for investigating embryo implantation and feto-maternal interactions, offering valuable insights for advancing reproductive medicine.

The impact of selective HDAC inhibitors on the transcriptome of early mouse embryos

BACKGROUND

Histone acetylation, which is regulated by histone acetyltransferases (HATs) and histone deacetylases (HDACs), plays a crucial role in the control of gene expression. HDAC inhibitors (HDACi) have shown potential in cancer therapy; however, the specific roles of HDACs in early embryos remain unclear. Moreover, although some pan-HDACi have been used to maintain cellular undifferentiated states in early embryos, the specific mechanisms underlying their effects remain unknown. Thus, there remains a significant knowledge gap regarding the application of selective HDACi in early embryos.

RESULTS

To address this gap, we treated early embryos with two selective HDACi (MGCD0103 and T247). Subsequently, we collected and analyzed their transcriptome data at different developmental stages. Our findings unveiled a significant effect of HDACi treatment during the crucial 2-cell stage of zygotes, leading to a delay in embryonic development after T247 and an arrest at 2-cell stage after MGCD0103 administration. Furthermore, we elucidated the regulatory targets underlying this arrested embryonic development, which pinpointed the G2/M phase as the potential period of embryonic development arrest caused by MGCD0103. Moreover, our investigation provided a comprehensive profile of the biological processes that are affected by HDACi, with their main effects being predominantly localized in four aspects of zygotic gene activation (ZGA): RNA splicing, cell cycle regulation, autophagy, and transcription factor regulation. By exploring the transcriptional regulation and epigenetic features of the genes affected by HDACi, we made inferences regarding the potential main pathways via which HDACs affect gene expression in early embryos. Notably, Hdac7 exhibited a distinct response, highlighting its potential as a key player in early embryonic development.

CONCLUSIONS

Our study conducted a comprehensive analysis of the effects of HDACi on early embryonic development at the transcriptional level. The results demonstrated that HDACi significantly affected ZGA in embryos, elucidated the distinct actions of various selective HDACi, and identified specific biological pathways and mechanisms via which these inhibitors modulated early embryonic development.

Exome-wide benchmark of difficult-to-sequence regions using short-read next-generation DNA sequencing

Next-generation DNA sequencing (NGS) in short-read mode has recently been used for genetic testing in various clinical settings. NGS data accuracy is crucial in clinical settings, and several reports regarding quality control of NGS data, primarily focusing on establishing NGS sequence read accuracy, have been published thus far. Variant calling is another critical source of NGS errors that remains unexplored at the single-nucleotide level despite its established significance. In this study, we used a machine-learning-based method to establish an exome-wide benchmark of difficult-to-sequence regions at the nucleotide-residue resolution using 10 genome sequence features based on real-world NGS data accumulated in The Genome Aggregation Database (gnomAD) of the human reference genome sequence (GRCh38/hg38). The newly acquired metric, designated the 'UNMET score,' along with additional lines of structural information from the human genome, allowed us to assess the sequencing challenges within the exonic region of interest using conventional short-read NGS. Thus, the UNMET score could provide a basis for addressing potential sequential errors in protein-coding exons of the human reference genome sequence GRCh38/hg38 in clinical sequencing.

CRISPR screening in human trophoblast stem cells reveals both shared and distinct aspects of human and mouse placental development

The placenta serves as the interface between the mother and fetus, facilitating the exchange of gases and nutrients between their separate blood circulation systems. Trophoblasts in the placenta play a central role in this process. Our current understanding of mammalian trophoblast development relies largely on mouse models. However, given the diversification of mammalian placentas, findings from the mouse placenta cannot be readily extrapolated to other mammalian species, including humans. To fill this knowledge gap, we performed CRISPR knockout screening in human trophoblast stem cells (hTSCs). We targeted genes essential for mouse placental development and identified more than 100 genes as critical regulators in both human hTSCs and mouse placentas. Among them, we further characterized in detail two transcription factors, DLX3 and GCM1, and revealed their essential roles in hTSC differentiation. Moreover, a gene function-based comparison between human and mouse trophoblast subtypes suggests that their relationship may differ significantly from previous assumptions based on tissue localization or cellular function. Notably, our data reveal that hTSCs may not be analogous to mouse TSCs or the extraembryonic ectoderm (ExE) in which in vivo TSCs reside. Instead, hTSCs may be analogous to progenitor cells in the mouse ectoplacental cone and chorion. This finding is consistent with the absence of ExE-like structures during human placental development. Our data not only deepen our understanding of human trophoblast development but also facilitate cross-species comparison of mammalian placentas.

Pan-cancer analysis of whole-genome doubling and its association with patient prognosis

BACKGROUND

Whole-genome doubling (WGD) is a common mutation in cancer. Various studies have suggested that WGD is associated with a poor prognosis in cancer. However, the detailed association between WGD occurrence and prognosis remains unclear. In this study, we aimed to elucidate the mechanism by which WGD affects prognosis using sequencing data from the Pan-Cancer Analysis of Whole Genomes (PCAWG) and The Cancer Genome Atlas.

METHODS

Whole-genome sequencing data of 23 cancer types were downloaded from PCAWG project. We defined the WGD event in each sample using the WGD status annotated using PCAWG. We used MutationTimeR to predict the relative timings of mutations and loss of heterozygosity (LOH) in WGD, thus evaluating their association with WGD. We also analyzed the association between WGD-associated factors and patient prognosis.

RESULTS

WGD was associated with several factors, e.g., length of LOH regions. Survival analysis using WGD-associated factors revealed that longer LOH regions and LOH in chr17 were associated with poor prognosis in samples with WGD (WGD samples) and samples without WGD (nWGD samples). In addition to these two factors, nWGD samples showed that the number of mutations in tumor suppressor genes was associated with prognosis. Moreover, we explored the genes associated with prognosis in both samples separately.

CONCLUSION

The prognosis-related factors in WGD samples differed significantly compared with those in nWGD samples. This study emphasizes the need for different treatment strategies for WGD and nWGD samples.

Possible involvement of silent mutations in cancer pathogenesis and evolution

Recent studies have shown that some silent mutations can be harmful to various processes. In this study, we performed a comprehensive in silico analysis to elucidate the effects of silent mutations on cancer pathogenesis using exome sequencing data derived from the Cancer Genome Atlas. We focused on the codon optimality scores of silent mutations, which were defined as the difference between the optimality of synonymous codons, calculated using the codon usage table. The relationship between cancer evolution and silent mutations showed that the codon optimality score of the mutations that occurred later in carcinogenesis was significantly higher than of those that occurred earlier. In addition, mutations with higher scores were enriched in genes involved in the cell cycle and cell division, while those with lower scores were enriched in genes involved in apoptosis and cellular senescence. Our results demonstrate that some silent mutations can be involved in cancer pathogenesis.

A case of breast milk-acquired cytomegalovirus infection in an extremely low birth weight infant

INTRODUCTION

Although breast milk is considered the optimal nutrition for infants, it is also the primary cause of postnatal cytomegalovirus (CMV) infection. Preterm infants with postnatal CMV infections are susceptible to a variety of life-threatening conditions.

CASE SUMMARY

Twin male infants were delivered via emergency caesarian section at 27 weeks' gestation secondary to maternal complete uterine rupture. The Apgar scores at 1 and 5 min were 1 and 1 for the older twin (Twin A) and 0 and 3 for the younger twin (Twin B). Their birth weights were 1203 g (+ 0.65SD) and 495 g (- 3.79SD) respectively. On day 41, laboratory blood test results for Twin B showed a moderate elevation in C-reactive protein (CRP), thrombocytopenia. CMV quantitative polymerase chain reaction (qPCR) tests in Twin B's urine and blood as well as in the mother's breast milk were positive, but stored, dried umbilical cord CMV qPCR tests were negative. Twin B was diagnosed with a postnatal CMV infection secondary to infected breast milk and ganciclovir was commenced on day 52. Treatment was switched to valganciclovir at 74 days of age, but a negative CMV-DNA level in the blood was not achieved. Postnatal CMV infection in this infant led to an exacerbation of pre-existing bronchopulmonary dysplasia (BPD) and he demised at 182 days of age.

CONCLUSION

Postnatal cytomegalovirus infections may lead to exacerbations of BPD. Early use of raw breast milk in preterm infants should be done with careful consideration of this potential complication.

Genome-wide identification of copy neutral loss of heterozygosity reveals its possible association with spatial positioning of chromosomes

Loss of heterozygosity (LOH) is a genetic alteration that results from the loss of one allele at a heterozygous locus. In particular, copy neutral LOH (CN-LOH) events are generated, for example, by mitotic homologous recombination after monoallelic defection or gene conversion, resulting in novel homozygous locus having two copies of the normal counterpart allele. This phenomenon can serve as a source of genome diversity and is associated with various diseases. To clarify the nature of the CN-LOH such as the frequency, genomic distribution and inheritance pattern, we made use of whole-genome sequencing data of the three-generation CEPH/Utah family cohort, with the pedigree consisting of grandparents, parents and offspring. We identified an average of 40.7 CN-LOH events per individual taking advantage of 285 healthy individuals from 33 families in the cohort. On average 65% of them were classified as gonosomal-mosaicism-associated CN-LOH, which exists in both germline and somatic cells. We also confirmed that the incidence of the CN-LOH has little to do with the parents' age and sex. Furthermore, through the analysis of the genomic region including the CN-LOH, we found that the chance of the occurrence of the CN-LOH tends to increase at the GC-rich locus and/or on the chromosome having a relatively close inter-homolog distance. We expect that these results provide significant insights into the association between genetic alteration and spatial position of chromosomes as well as the intrinsic genetic property of the CN-LOH.

Increased neutrophils in inflammatory bowel disease accelerate the accumulation of amyloid plaques in the mouse model of Alzheimer's disease

BACKGROUND

Alzheimer's disease (AD) is one of the neurodegenerative diseases and characterized by the appearance and accumulation of amyloid-β (Aβ) aggregates and phosphorylated tau with aging. The aggregation of Aβ, which is the main component of senile plaques, is closely associated with disease progression. App^NL-G-F mice, a mouse model of AD, have three familial AD mutations in the amyloid-β precursor gene and exhibit age-dependent AD-like symptoms and pathology. Gut-brain interactions have attracted considerable attention and inflammatory bowel disease (IBD) has been associated with a higher risk of dementia, especially AD, in humans. However, the underlying mechanisms and the effects of intestinal inflammation on the brain in AD remain largely unknown. Therefore, we aimed to investigate the effects of intestinal inflammation on AD pathogenesis.

METHODS

Wild-type and App^NL-G-F mice at three months of age were fed with water containing 2% dextran sulfate sodium (DSS) to induce colitis. Immune cells in the brain were analyzed using single-cell RNA sequencing (scRNA-seq) analysis, and the aggregation of Aβ protein in the brain was analyzed via immunohistochemistry.

RESULTS

An increase in aggregated Aβ was observed in the brains of App^NL-G-F mice with acute intestinal inflammation. Detailed scRNA-seq analysis of immune cells in the brain showed that neutrophils in the brain increased after acute enteritis. Eliminating neutrophils by antibodies suppressed the accumulation of Aβ, which increased because of intestinal inflammation.

CONCLUSION

These results suggest that neutrophils infiltrate the AD brain parenchyma when acute colitis occurs, and this infiltration is significantly related to disease progression. Therefore, we propose that neutrophil-targeted therapies could reduce Aβ accumulation observed in early AD and prevent the increased risk of AD due to colitis.

Tmsb10 triggers fetal Leydig differentiation by suppressing the RAS/ERK pathway

Leydig cells in fetal testes play crucial roles in masculinizing fetuses through androgen production. Gene knockout studies have revealed that growth factors are implicated in fetal Leydig cell (FLC) differentiation, but little is known about the mechanisms regulating this process. We investigate this issue by characterizing FLC progenitor cells using single-cell RNA sequencing. The sequence datasets suggest that thymosin β10 (Tmsb10) is transiently upregulated in the progenitors. While studying the function of Tmsb10, we reveal that platelet-derived growth factor (PDGF) regulates ciliogenesis through the RAS/ERK and PI3K/AKT pathways, and thereby promotes desert hedgehog (DHH)-dependent FLC differentiation. Tmsb10 expressed in the progenitor cells induces their differentiation into FLCs by suppressing the RAS/ERK pathway. Through characterizing the transiently expressed Tmsb10 in the FLC progenitors, this study unveils the molecular process of FLC differentiation and shows that it is cooperatively induced by DHH and PDGF.

Investigation of fetal Leydig progenitors shows that thymosin β10 (Tmsb10) suppresses the RAS/ERK pathway, inducing progenitor differentiation into fetal Leydig cells.

Mapping of promoter usage QTL using RNA-seq data reveals their contributions to complex traits

Genomic variations are associated with gene expression levels, which are called expression quantitative trait loci (eQTL). Most eQTL may affect the total gene expression levels by regulating transcriptional activities of a specific promoter. However, the direct exploration of genomic loci associated with promoter activities using RNA-seq data has been challenging because eQTL analyses treat the total expression levels estimated by summing those of all isoforms transcribed from distinct promoters. Here we propose a new method for identifying genomic loci associated with promoter activities, called promoter usage quantitative trait loci (puQTL), using conventional RNA-seq data. By leveraging public RNA-seq datasets from the lymphoblastoid cell lines of 438 individuals from the GEUVADIS project, we obtained promoter activity estimates and mapped 2,592 puQTL at the 10% FDR level. The results of puQTL mapping enabled us to interpret the manner in which genomic variations regulate gene expression. We found that 310 puQTL genes (16.1%) were not detected by eQTL analysis, suggesting that our pipeline can identify novel variant–gene associations. Furthermore, we identified genomic loci associated with the activity of “hidden” promoters, which the standard eQTL studies have ignored. We found that most puQTL signals were concordant with at least one genome-wide association study (GWAS) signal, enabling novel interpretations of the molecular mechanisms of complex traits. Our results emphasize the importance of the re-analysis of public RNA-seq datasets to obtain novel insights into gene regulation by genomic variations and their contributions to complex traits.

Many variations exist in the human genome, creating phenotypic diversity among individuals. It is well known that they are associated with the risk of disease development by affecting the expression levels of genes. Genes are transcribed from regulatory elements called promoters. Although some genes are transcribed from multiple promoters and translated into proteins with different functions, the relationship between genomic variations and promoter activities has not been investigated in depth compared to the relationship between genomic variations and gene expression levels. In this study, we proposed a new method to detect the association between genomic variations and promoter activities. Our method identified the associations between many variations and promoters using genomic and promoter activity data from blood cells of 438 individuals. This study allowed us to identify new functional associations between genomic variations and genes. Furthermore, we identified previously undiscovered variation-gene-disease associations. Our results will help to elucidate the molecular mechanisms of diseases in which genetic factors are involved.

In Vitro Generation of Brain Regulatory T Cells by Co-culturing With Astrocytes

Regulatory T cells (Tregs) are normally born in the thymus and activated in secondary lymphoid tissues to suppress immune responses in the lymph node and at sites of inflammation. Tregs are also resident in various tissues or accumulate in damaged tissues, which are now called tissue Tregs, and contribute to homeostasis and tissue repair by interacting with non-immune cells. We have shown that Tregs accumulate in the brain during the chronic phase in a mouse cerebral infarction model, and these Tregs acquire the characteristic properties of brain Tregs and contribute to the recovery of neurological damage by interacting with astrocytes. However, the mechanism of tissue Treg development is not fully understood. We developed a culture method that confers brain Treg characteristics in vitro. Naive Tregs from the spleen were activated and efficiently amplified by T-cell receptor (TCR) stimulation in the presence of primary astrocytes. Furthermore, adding IL-33 and serotonin could confer part of the properties of brain Tregs, such as ST2, peroxisome proliferator-activated receptor γ (PPARγ), and serotonin receptor 7 (Htr7) expression. Transcriptome analysis revealed that in vitro generated brain Treg-like Tregs (induced brain Tregs; iB-Tregs) showed similar gene expression patterns as those in in vivo brain Tregs, although they were not identical. Furthermore, in Parkinson’s disease models, in which T cells have been shown to be involved in disease progression, iB-Tregs infiltrated into the brain more readily and ameliorated pathological symptoms more effectively than splenic Tregs. These data indicate that iB-Tregs contribute to our understanding of brain Treg development and could also be therapeutic for inflammatory brain diseases.

Identification of a drug-response gene in multiple myeloma through longitudinal single-cell transcriptome sequencing

Despite recent therapeutic advances for multiple myeloma (MM), relapse is very common. Here, we conducted longitudinal single-cell transcriptome sequencing (scRNA-seq) of MM cells from a patient with relapsed MM, treated with multiple anti-myeloma drugs. We observed five subclusters of MM cells, which appeared and/or disappeared in response to the therapeutic pressure, and identified cluster 3 which emerged during lenalidomide treatment and disappeared after proteasome inhibitor (PI) treatment. Among the differentially expressed genes in cluster 3, we found a candidate drug-response gene; pellino E3 ubiquitin-protein ligase family member 2 (PELI2), which is responsible for PI-induced cell death in in vitro assay. Kaplan-Meier survival analysis of database revealed that higher expression of PELI2 is associated with a better prognosis. Our integrated strategy combining longitudinal scRNA-seq analysis, in vitro functional assay, and database analysis would facilitate the understanding of clonal dynamics of MM in response to anti-myeloma drugs and identification of drug-response genes.

•

Longitudinal scRNA-seq reveals clonal dynamics of MM under therapeutic pressure

•

PELI2 is identified as a candidate proteasome inhibitors (PI)-sensitive gene from the PI-sensitive cluster

•

Overexpression of PELI2 sensitizes PI to an MM cell line in the cytotoxic assay

•

In database analysis, high expression of PELI2 is associated with a better prognosis

The microRNA cluster C19MC confers differentiation potential into trophoblast lineages upon human pluripotent stem cells

The first cell fate commitment during mammalian development is the specification of the inner cell mass and trophectoderm. This irreversible cell fate commitment should be epigenetically regulated, but the precise mechanism is largely unknown in humans. Here, we show that naïve human embryonic stem (hES) cells can transdifferentiate into trophoblast stem (hTS) cells, but primed hES cells cannot. Our transcriptome and methylome analyses reveal that a primate-specific miRNA cluster on chromosome 19 (C19MC) is active in naïve hES cells but epigenetically silenced in primed ones. Moreover, genome and epigenome editing using CRISPR/Cas systems demonstrate that C19MC is essential for hTS cell maintenance and C19MC-reactivated primed hES cells can give rise to hTS cells. Thus, we reveal that C19MC activation confers differentiation potential into trophoblast lineages on hES cells. Our findings are fundamental to understanding the epigenetic regulation of human early development and pluripotency.

Little is known about the epigenetic mechanisms of the first cell fate commitment in humans. Here, the authors show that activation of the miRNA cluster C19MC confers differentiation potential into trophoblast lineages on human embryonic stem cells.

Pan-cancer analysis of mutations in open chromatin regions and their possible association with cancer pathogenesis

BACKGROUND

Open chromatin is associated with gene transcription. Previous studies have shown that the density of mutations in open chromatin regions is lower than that in flanking regions because of the higher accessibility of DNA repair machinery. However, in several cancer types, open chromatin regions show an increased local density of mutations in activated regulatory regions. Although the mutation distribution within open chromatin regions in cancer cells has been investigated, only few studies have focused on their functional implications in cancer. To reveal the impact of highly mutated open chromatin regions on cancer, we investigated the association between mutations in open chromatin regions and their possible functions.

METHODS

Whole-genome sequencing data of 18 cancer types were downloaded from the PanCancer Analysis of Whole Genomes and Catalog of Somatic Mutations in Cancer. We quantified the mutations located in open chromatin regions defined by The Cancer Genome Atlas and classified open chromatin regions into three categories based on the number of mutations. Then, we investigated the chromatin state, amplification, and possible target genes of the open chromatin regions with a high number of mutations. We also analyzed the association between the number of mutations in open chromatin regions and patient prognosis.

RESULTS

In some cancer types, the proportion of promoter or enhancer chromatin state in open chromatin regions with a high number of mutations was significantly higher than that in the regions with a low number of mutations. The possible target genes of open chromatin regions with a high number of mutations were more strongly associated with cancer than those of other open chromatin regions. Moreover, a high number of mutations in open chromatin regions was significantly associated with a poor prognosis in some cancer types.

CONCLUSIONS

These results suggest that highly mutated open chromatin regions play an important role in cancer pathogenesis and can be effectively used to predict patient prognosis.

Kastor and Polluks polypeptides encoded by a single gene locus cooperatively regulate VDAC and spermatogenesis

Although several long noncoding RNAs (lncRNAs) have recently been shown to encode small polypeptides, those in testis remain largely uncharacterized. Here we identify two sperm-specific polypeptides, Kastor and Polluks, encoded by a single mouse locus (Gm9999) previously annotated as encoding a lncRNA. Both Kastor and Polluks are inserted in the outer mitochondrial membrane and directly interact with voltage-dependent anion channel (VDAC), despite their different amino acid sequences. Male VDAC3-deficient mice are infertile as a result of reduced sperm motility due to an abnormal mitochondrial sheath in spermatozoa, and deficiency of both Kastor and Polluks also severely impaired male fertility in association with formation of a similarly abnormal mitochondrial sheath. Spermatozoa lacking either Kastor or Polluks partially recapitulate the phenotype of those lacking both. Cooperative function of Kastor and Polluks in regulation of VDAC3 may thus be essential for mitochondrial sheath formation in spermatozoa and for male fertility.

A number of testes-specific lncRNAs have been annotated but their roles remain largely unexplored. Here the authors identify two small peptides encoded by the lncRNA Gm9999, Kastor and Polluks, both of which are required for male fertility in mice.

Genome-wide identification of exon extension/shrinkage events induced by splice-site-creating mutations

Mutations that affect phenotypes have been identified primarily as those that directly alter amino acid sequences or disrupt splice sites. However, some mutations not located in functionally important sites can also affect phenotypes, such as splice-site-creating mutations (SCMs). To investigate how frequent exon extension/shrinkage events induced by SCMs occur in normal individuals, we used personal genome sequencing data and transcriptome data of the corresponding individuals and identified 371 exon extension/shrinkage events in normal individuals. This number was about three times higher than the number of pseudo-exon activation events identified in the previous study. The average numbers of exon extension and exon shrinkage events in each sample were 3.3 and 11.2, respectively. We also evaluated the impact of exon extension/shrinkage events on the resulting transcripts and their protein products and found that 40.2% of the identified events may have possible functional impacts by either generating premature termination codons in transcripts or affecting protein domains. Our results indicated that a certain fraction of SCMs identified in this study can be pathogenic mutations by creating novel splice sites.

Pan-cancer analysis of non-coding recurrent mutations and their possible involvement in cancer pathogenesis

Cancer-related mutations have been mainly identified in protein-coding regions. Recent studies have demonstrated that mutations in non-coding regions of the genome could also be a risk factor for cancer. However, the non-coding regions comprise 98% of the total length of the human genome and contain a huge number of mutations, making it difficult to interpret their impacts on pathogenesis of cancer. To comprehensively identify cancer-related non-coding mutations, we focused on recurrent mutations in non-coding regions using somatic mutation data from COSMIC and whole-genome sequencing data from The Cancer Genome Atlas (TCGA). We identified 21 574 recurrent mutations in non-coding regions that were shared by at least two different samples from both COSMIC and TCGA databases. Among them, 580 candidate cancer-related non-coding recurrent mutations were identified based on epigenomic and chromatin structure datasets. One of such mutation was located in RREB1 binding site that is thought to interact with TEAD1 promoter. Our results suggest that mutations may disrupt the binding of RREB1 to the candidate enhancer region and increase TEAD1 expression levels. Our findings demonstrate that non-coding recurrent mutations and coding mutations may contribute to the pathogenesis of cancer.

In silico identification of pseudo-exon activation events in personal genome and transcriptome data

Causative mutations for human genetic disorders have mainly been identified in exonic regions that code for amino acid sequences. Recently, however, it has been reported that mutations in deep intronic regions can also cause certain human genetic disorders by creating novel splice sites, leading to pseudo-exon activation. To investigate how frequently pseudo-exon activation events occur in normal individuals, we conducted in silico identification of such events using personal genome data and corresponding high-quality transcriptome data. With rather stringent conditions, on average, 2.6 pseudo-exon activation events per individual were identified. More pseudo-exon activation events were found in 5' donor splice sites than in 3' acceptor splice sites. Although pseudo-exon activation events have sporadically been reported as causative mutations in genetic disorders, it is revealed in this study that such events can be observed in normal individuals at a certain frequency. We estimate that human genomes typically contain on average at least 10 pseudo-exon activation events. The actual number should be higher than this, because we used stringent criteria to identify pseudo-exon activation events. This suggests that it is worth considering the possibility of pseudo-exon activation when searching for causative mutations of genetic disorders if candidate mutations are not identified in coding regions or RNA splice sites.

Direct reprogramming of human umbilical vein- and peripheral blood-derived endothelial cells into hepatic progenitor cells

Recent advances have enabled the direct induction of human tissue-specific stem and progenitor cells from differentiated somatic cells. However, it is not known whether human hepatic progenitor cells (hHepPCs) can be generated from other cell types by direct lineage reprogramming with defined transcription factors. Here, we show that a set of three transcription factors, FOXA3, HNF1A, and HNF6, can induce human umbilical vein endothelial cells to directly acquire the properties of hHepPCs. These induced hHepPCs (hiHepPCs) propagate in long-term monolayer culture and differentiate into functional hepatocytes and cholangiocytes by forming cell aggregates and cystic epithelial spheroids, respectively, under three-dimensional culture conditions. After transplantation, hiHepPC-derived hepatocytes and cholangiocytes reconstitute damaged liver tissues and support hepatic function. The defined transcription factors also induce hiHepPCs from endothelial cells circulating in adult human peripheral blood. These expandable and bipotential hiHepPCs may be useful in the study and treatment of human liver diseases.

The conditions to induce human hepatic progenitor cells from other cell types are unclear. Here, the authors reprogram human endothelial cells to hepatic progenitor cells by expressing FOXA3, HNF1A and HNF6, capable of giving rise to hepatocytes and cholangiocytes that reconstitute damaged liver tissues on transplantation.

Zfp281 Shapes the Transcriptome of Trophoblast Stem Cells and Is Essential for Placental Development

Placental development is a key event in mammalian reproduction and embryogenesis. However, the molecular basis underlying placental development is not fully understood. Here, we conduct a forward genetic screen to identify regulators for extraembryonic development and identify Zfp281 as a key factor. Zfp281 overexpression in mouse embryonic stem cells facilitates the induction of trophoblast stem-like cells. Zfp281 is preferentially expressed in the undifferentiated trophoblast stem cell population in an FGF-dependent manner, and disruption of Zfp281 in mice causes severe defects in early placental development. Consistently, Zfp281-depleted trophoblast stem cells exhibit defects in maintaining the transcriptome and differentiation capacity. Mechanistically, Zfp281 interacts with MLL or COMPASS subunits and occupies the promoters of its target genes. Importantly, ZNF281, the human ortholog of this factor, is required to stabilize the undifferentiated status of human trophoblast stem cells. Thus, we identify Zfp281 as a conserved factor for the maintenance of trophoblast stem cell plasticity.

Augmented oxidative stress increases 8-oxoguanine preferentially in the transcriptionally active genomic regions

8-Oxoguanine (8-oxoG) is the most common DNA base modification in the mammalian genome, associated with oxidative stress. Here we analysed the alterations in the distribution of 8-oxoG across the entire murine genome, before and after an elevation of oxidative stress by the use of ferric nitrilotriacetate (Fe-NTA) as an oxidative stress inducer in the renal proximal tubules. We isolated DNA fragments containing 8-oxoGs with immunoprecipitation from the murine genome, and amplified them by PCR for a distribution analysis with microarray-based comparative genomic hybridisation. The distribution profiles revealed that frequencies of 8-oxoG fluctuated with a cycle of 1-10 Mb along the chromosomes and the amplitude of the fluctuation was reduced after Fe-NTA administration. The distributions of 8-oxoG along the entire genome in the control and oxidatively stressed conditions were negatively correlated with that of gene density but positively correlated with that of Lamin B1 interaction, which corresponds to lamina-associated domains. These results on the murine genome were consistent with those on the rat genome we previously reported. We further discovered a negative correlation between the distributions of 8-oxoG and transcriptional activity along the genome. Finally, a comparison of the distributions before and after Fe-NTA administration suggested that 8-oxoGs are generated in response to the augmented oxidative stress preferentially in the transcriptionally active genomic regions, where 8-oxoGs have been less accumulated in the control condition.

Comparative genomic analysis of inbred rat strains reveals the existence of ancestral polymorphisms

In an alignment of closely related genomic sequences, the existence of discordant mutation sites, which do not reflect the phylogenetic relationship of the genomes, is often observed. Although these discordant mutation sites are thought to have emerged by ancestral polymorphism or gene flow, their frequency and distribution in the genome have not yet been analyzed in detail. Using the genome sequences of all protein coding genes of 25 inbred rat strains, we analyzed the frequency and genome-wide distribution of the discordant mutation sites. From the comparison of different substrains, it was found that these loci are not substrain specific, but are common among different groups of substrains, suggesting that the discordant sites might have mainly emerged through ancestral polymorphism. It was also revealed that the discordant sites are not uniformly distributed along chromosomes, but are concentrated at certain genomic loci, such as RT1, major histocompatibility complex of rats, and olfactory receptors, indicating that genes known to be highly polymorphic tend to have more discordant sites. Our results also showed that loci with a high density of discordant sites are also rich in heterozygous variants, even though these are inbred strains.

An integrated analysis of public genomic data unveils a possible functional mechanism of psoriasis risk via a long-range ERRFI1 enhancer

BACKGROUND

Psoriasis is a chronic inflammatory skin disease, for which genome-wide association studies (GWAS) have identified many genetic variants as risk markers. However, the details of underlying molecular mechanisms, especially which variants are functional, are poorly understood.

METHODS

We utilized a computational approach to survey psoriasis-associated functional variants that might affect protein functions or gene expression levels. We developed a pipeline by integrating publicly available datasets provided by GWAS Catalog, FANTOM5, GTEx, SNP2TFBS, and DeepBlue. To identify functional variants on exons or splice sites, we used a web-based annotation tool in the Ensembl database. To search for noncoding functional variants within promoters or enhancers, we used eQTL data calculated by GTEx. The data of variants lying on transcription factor binding sites provided by SNP2TFBS were used to predict detailed functions of the variants.

RESULTS

We discovered 22 functional variant candidates, of which 8 were in noncoding regions. We focused on the enhancer variant rs72635708 (T > C) in the 1p36.23 region; this variant is within the enhancer region of the ERRFI1 gene, which regulates lipid metabolism in the liver and skin morphogenesis via EGF signaling. Further analysis showed that the ERRFI1 promoter spatially contacts with the enhancer, despite the 170 kb distance between them. We found that this variant lies on the AP-1 complex binding motif and may modulate binding levels.

CONCLUSIONS

The minor allele rs72635708 (rs72635708-C) might affect the ERRFI1 promoter activity, which results in unstable expression of ERRFI1, enhancing the risk of psoriasis via disruption of lipid metabolism and skin cell proliferation. Our study represents a successful example of predicting molecular pathogenesis by integration and reanalysis of public data.

DBTSS/DBKERO for integrated analysis of transcriptional regulation

DBTSS (Database of Transcriptional Start Sites)/DBKERO (Database of Kashiwa Encyclopedia for human genome mutations in Regulatory regions and their Omics contexts) is the database originally initiated with the information of transcriptional start sites and their upstream transcriptional regulatory regions. In recent years, we updated the database to assist users to elucidate biological relevance of the human genome variations or somatic mutations in cancers which may affect the transcriptional regulation. In this update, we facilitate interpretations of disease associated genomic variation, using the Japanese population as a model case. We enriched the genomic variation dataset consisting of the 13,368 individuals collected for various genome-wide association studies and the reference epigenome information in the surrounding regions using a total of 455 epigenome datasets (four tissue types from 67 healthy individuals) collected for the International Human Epigenome Consortium (IHEC). The data directly obtained from the clinical samples was associated with that obtained from various model systems, such as the drug perturbation datasets using cultured cancer cells. Furthermore, we incorporated the results obtained using the newly developed analytical methods, Nanopore/10x Genomics long-read sequencing of the human genome and single cell analyses. The database is made publicly accessible at the URL (http://dbtss.hgc.jp/).

CLEC3A, MMP7, and LCN2 as novel markers for predicting recurrence in resected G1 and G2 pancreatic neuroendocrine tumors

Although the postoperative recurrence rate for pancreatic neuroendocrine tumors (PNETs) is reported to be 13.5%-30%, the paucity of valuable biomarkers to predict recurrence poses a problem for the early detection of relapse. Hence, this study aimed to identify new biomarkers to predict the recurrence of PNETs. We performed RNA sequencing (RNA-Seq) on RNA isolated from frozen primary tumors sampled from all localized G1/G2 PNETs resected curatively from 1998 to 2015 in our institution. We calculated differentially expressed genes (DEGs) in tumor with and without recurrence (≥3 years) for the propensity-matched cohort. Gene ontology analysis for the identified DEGs was also performed. Furthermore, we evaluated the expression levels of candidate genes as recurrence predictors via immunostaining. Comparison of transcriptional levels in tumors with and without recurrence identified 166 DEGs. Up- and downregulated genes with high significance in these tumors were mainly related to extracellular organization and cell adhesion, respectively. We observed the top three upregulated genes, C-type lectin domain family 3 member A (CLEC3A), matrix metalloproteinase-7 (MMP7), and lipocalin2 (LCN2) immunohistochemically and compared their levels in recurrent and nonrecurrent tumors. Significantly higher recurrence rate was shown in patients with positive expression of CLEC3A (P = 0.028), MMP7 (P = 0.003), and LCN2 (P = 0.040) than that with negative expression. We identified CLEC3A, MMP7, and LCN2 known to be associated with the phosphatidylinositol-3-kinase/Akt pathway, as potential novel markers to predict the postoperative recurrence of PNETs.

Mouse polycomb group gene Cbx2 promotes osteoblastic but suppresses adipogenic differentiation in postnatal long bones

A set of key developmental genes is essential for skeletal growth from multipotent progenitor cells at weaning. Polycomb group proteins, which regulate such genes contributes to the cell lineage commitment and subsequent differentiation via epigenetic chromatin modification and remodeling. However, it is unclear which cell lineage and gene sets are targeted by polycomb proteins during skeletal growth. We now report that mice deficient in a polycomb group gene Cbx2^cterm/cterm exhibited skeletal hypoplasia in the tibia, femur, and cranium. Long bone cavities in these mice contained fewer multipotent mesenchymal stromal cells. RNA-sequencing of bone marrow cells showed downregulation and upregulation of osteoblastic and adipogenic genes, respectively. Furthermore, the expression levels of genes specifically expressed in B-cell precursors were decreased. Forced expression of Cbx2 in Cbx2^cterm/cterm bone marrow stromal cell recovered fibroblastic colony formation and suppressed adipogenic differentiation. Collectively, our results suggest that Cbx2 controls the maintenance and adipogenic differentiation of mesenchymal stromal cells in the bone marrow.

Pan-cancer analysis of intratumor heterogeneity associated with patient prognosis using multidimensional measures

Human cancers accumulate various mutations during development and consist of highly heterogeneous cell populations. This phenomenon is called intratumor heterogeneity (ITH). ITH is known to be involved in tumor growth, progression, invasion, and metastasis, presenting obstacles to accurate diagnoses and effective treatments. Numerous studies have explored the dynamics of ITH, including constructions of phylogenetic trees in cancer samples using multiregional ultradeep sequencing and simulations of evolution using statistical models. Although ITH is associated with prognosis, it is still challenging to use the characteristics of ITH as prognostic factors because of difficulties in quantifying ITH precisely. In this study, we analyzed the relationship between patient prognosis and the distribution of variant allele frequencies (VAFs) in cancer samples (n = 6,064) across 16 cancer types registered in The Cancer Genome Atlas. To measure VAF distributions multidimensionally, we adopted parameters that define the shape of VAF distributions and evaluated the relationships between these parameters and prognosis. In seven cancer types, we found significant relationships between prognosis and VAF distributions. Moreover, we observed that samples with a larger amount of mutations were not necessarily linked to worse prognosis. By evaluating the ITH from multidimensional viewpoints, it will be possible to provide a more accurate prediction of cancer prognosis.

Fetal Leydig cells dedifferentiate and serve as adult Leydig stem cells

Previous studies have established that fetal Leydig cells (FLCs) and adult Leydig cells (ALCs) show distinct functional characteristics. However, the lineage relationship between FLCs and ALCs has not been clarified yet. Here, we reveal that a subset of FLCs dedifferentiate at fetal stages to give rise to ALCs at the pubertal stage. Moreover, the dedifferentiated cells contribute to the peritubular myoid cell and vascular pericyte populations in the neonatal testis, and these non-steroidogenic cells serve as potential ALC stem cells. We generated FLC lineage-specific Nr5a1 (Ad4BP/SF-1) gene-disrupted mice and mice lacking the fetal Leydig enhancer (FLE) of the Nr5a1 gene. Phenotypes of these mice support the conclusion that most of the ALCs arise from dedifferentiated FLCs, and that the FLE of the Nr5a1 gene is essential for both initial FLC differentiation and pubertal ALC redifferentiation.

Intratumor heterogeneity of HMCN1 mutant alleles associated with poor prognosis in patients with breast cancer

Human breast cancers comprise a complex and highly heterogeneous population of tumor cells. Intratumor heterogeneity is an underlying cause of resistance to effective therapies and disease recurrence. To explore prognostic factors based on intratumor heterogeneity, we analyzed genomic mutations in breast cancer patients registered in The Cancer Genome Atlas. We calculated the variant allele frequency (VAF) at each mutation site and evaluated the associations of VAFs with the prognosis of breast cancer. VAFs of HMCN1 correlated with the prognosis and lymph node status. Although the detailed function of HMCN1 remains unknown, it is located in extracellular matrix and the mutation in the gene might be associated with cancer cell invasion and metastasis. This finding suggests that HMCN1 is a potential metastatic factor and can be a candidate gene for targeted breast cancer therapy.

Ad4BP/SF-1 regulates cholesterol synthesis to boost the production of steroids

Housekeeping metabolic pathways such as glycolysis are active in all cell types. In addition, many types of cells are equipped with cell-specific metabolic pathways. To properly perform their functions, housekeeping and cell-specific metabolic pathways must function cooperatively. However, the regulatory mechanisms that couple metabolic pathways remain largely unknown. Recently, we showed that the steroidogenic cell-specific nuclear receptor Ad4BP/SF-1, which regulates steroidogenic genes, also regulates housekeeping glycolytic genes. Here, we identify cholesterogenic genes as the targets of Ad4BP/SF-1. Further, we reveal that Ad4BP/SF-1 regulates Hummr, a candidate mediator of cholesterol transport from endoplasmic reticula to mitochondria. Given that cholesterol is the starting material for steroidogenesis and is synthesized from acetyl-CoA, which partly originates from glucose, our results suggest that multiple biological processes involved in synthesizing steroid hormones are governed by Ad4BP/SF-1. To our knowledge, this study provides the first example where housekeeping and cell-specific metabolism are coordinated at the transcriptional level.

Identification of a candidate enhancer for DMRT3 involved in spastic cerebral palsy pathogenesis

Cerebral palsy (CP) is a major neuronal disease and the most common movement disorder in children. Although environmental factors leading to CP have been greatly investigated, the genetic mechanism underlying CP is not well understood. Here we focused on two clinical reports that characterized a deletion involving the KANK1 gene locus in the 9p24.3 region. One report shows spastic CP and the other shows no spastic CP phenotype. Based on the epigenetic status and evolutionary conservation, we first found a functional genomic element at the noncoding region that was deleted only in patients with spastic CP. This element contains the retinoic acid receptor/retinoid X receptor (RAR/RXR) complex-binding motif that is widely conserved among placental mammals. RAR/RXR ChIP-seq data from mouse F9 embryonal carcinoma cells that were treated with trans-retinoic acids showed that the element has a binding ability. In addition, data regarding chromosome conformation capture from mouse neural progenitor and ES cells suggested that the element spatially interacts with the Doublesex and mab-3 related transcription factor 3 (Dmrt3) gene promoter that is located approximately 120 kb downstream of the RAR/RXR-binding site. Dmrt3 is detected in the developing mouse forebrain and in some interneurons in the spinal cord, and it works as a locomotion coordinator in horses and mice. Thus, the deletion of the cis-regulatory element for DMRT3 in humans may cause impaired development of the forebrain and gait abnormalities, resulting in spastic CP. In conclusion, this study provides new mechanistic insights into the genetic basis of CP.

Derivation of Human Trophoblast Stem Cells

Trophoblast cells play an essential role in the interactions between the fetus and mother. Mouse trophoblast stem (TS) cells have been derived and used as the best in vitro model for molecular and functional analysis of mouse trophoblast lineages, but attempts to derive human TS cells have so far been unsuccessful. Here we show that activation of Wingless/Integrated (Wnt) and EGF and inhibition of TGF-β, histone deacetylase (HDAC), and Rho-associated protein kinase (ROCK) enable long-term culture of human villous cytotrophoblast (CT) cells. The resulting cell lines have the capacity to give rise to the three major trophoblast lineages, which show transcriptomes similar to those of the corresponding primary trophoblast cells. Importantly, equivalent cell lines can be derived from human blastocysts. Our data strongly suggest that the CT- and blastocyst-derived cell lines are human TS cells, which will provide a powerful tool to study human trophoblast development and function.

Identification of Secretory Leukoprotease Inhibitor As an Endogenous Negative Regulator in Allergic Effector Cells

Mast cells, basophils, and eosinophils are central effectors in allergic inflammatory disorders. These cells secrete abundant serine proteases as well as chemical mediators and cytokines; however, the expression profiles and functions of their endogenous inhibitors remain elusive. We found that murine secretory leukoprotease inhibitor (SLPI) is expressed in basophils and eosinophils but in not in mast cells. SLPI-deficient (Slpi^−/−) basophils produce more cytokines than wild-type mice after IgE stimulation. Although the deletion of SLPI in basophils did not affect the release of chemical mediators upon IgE stimulation, the enzymatic activity of the serine protease tryptase was increased in Slpi^−/− basophils. Mice transferred with Slpi^−/− basophils were highly sensitive to IgE-mediated chronic allergic inflammation. Eosinophils lacking SLPI showed greater interleukin-6 secretion and invasive activity upon lipopolysaccharide stimulation, and the expression of matrix metalloproteinase-9 by these eosinophils was increased without stimulation. The absence of SLPI increases JNK1 phosphorylation at the steady state, and augments the serine phosphorylation of JNK1-downstream ETS transcriptional factor Elk-1 in eosinophils upon stimulation. Of note, SLPI interacts with a scaffold protein, JNK-interacting protein 3 (JIP3), that constitutively binds to the cytoplasmic domain of toll-like receptor (TLR) 4, suggesting that SLPI controls Elk-1 activation via binding to JIP3 in eosinophils. Mice transferred with Slpi^−/− eosinophils showed the exacerbation of chitin-induced allergic inflammation. These findings showed that SLPI is a negative regulator in allergic effector cells and suggested a novel inhibitory role of SLPI in the TLR4 signaling pathways.

Histone methyltransferase G9a is a key regulator of the starvation-induced behaviors in Drosophila melanogaster

Organisms have developed behavioral strategies to defend themselves from starvation stress. Despite of their importance in nature, the underlying mechanisms have been poorly understood. Here, we show that Drosophila G9a (dG9a), one of the histone H3 Lys 9-specific histone methyltransferases, functions as a key regulator for the starvation-induced behaviors. RNA-sequencing analyses utilizing dG9a null mutant flies revealed that the expression of some genes relating to gustatory perception are regulated by dG9a under starvation conditions. Reverse transcription quantitative-PCR analyses showed that the expression of gustatory receptor genes for sensing sugar are up-regulated in starved dG9a null mutant. Consistent with this, proboscis extension reflex tests indicated that dG9a depletion increased the sensitivity to sucrose under starvation conditions. Furthermore, the locomotion activity was promoted in starved dG9a null mutant. We also found that dG9a depletion down-regulates the expression of insulin-like peptide genes that are required for the suppression of starvation-induced hyperactivity. Furthermore, refeeding of wild type flies after starvation conditions restores the hyperactivity and increased sensitivity to sucrose as well as dG9a expression level. These data suggest that dG9a functions as a key regulator for the decision of behavioral strategies under starvation conditions.

Role of Ad4-binding protein/steroidogenic factor 1 in regulating NADPH production in adrenocortical Y-1 cells

Ad4-binding protein/steroidogenic factor 1 (Ad4BP/SF-1), a member of the nuclear receptor superfamily, is expressed in steroidogenic cells and regulates all steroidogenic gene expression. We recently employed mRNA and chromatin immunoprecipitation sequence (ChIP-seq) to demonstrate that Ad4BP/SF-1 directly regulates the expression of nearly all glycolytic genes. The pentose phosphate pathway (PPP) contributes to the production of nicotinamide adenine dinucleotide phosphate (NADPH). Although the expression of PPP genes and intracellular NADPH were decreased by Ad4BP/SF-1 knockdown, these genes were not the direct targets of Ad4BP/SF-1. This study therefore investigates whether Ad4BP/SF-1 directly regulates genes implicated in NADPH production. Examination of previously published data sets of mRNA sequence (mRNA-seq) and ChIP-seq strongly suggested a possibility that other NADPH-producing genes, such as malic enzyme 1 (Me1) and methylenetetrahydrofolate dehydrogenase 2 (Mthfd2), are the direct targets of Ad4BP/SF-1. Reporter gene assays and determination of intracellular NADPH concentration supported the notion that Ad4BP/SF-1 regulates NADPH production by regulating these genes. NADPH is required for macromolecule synthesis of compounds such as steroids, and for detoxification of reactive oxygen species. When synthesizing steroid hormones, steroidogenic cells consume NADPH through enzymatic reactions mediated by steroidogenic P450s. NADPH is also consumed through elimination of reactive oxygen species produced as the byproducts of the P450 reactions. Overall, Ad4BP/SF-1 potentially maintains the intracellular NADPH level through cooperative regulation of genes involved in the biological processes for consumption and supply.

A deletion in the intergenic region upstream of Ednrb causes head spot in the rat strain KFRS4/Kyo

Background

Head spot is one of the phenotypes identified in the KFRS4/Kyo rat strain. Although previous linkage analysis suggested that Ednrb, which is frequently involved in coat color variations in various animals, could be the gene responsible for this phenotype, no mutations have been identified in its coding region.

Results

To identify mutations causative of this phenotype in KFRS4/Kyo, we analyzed target capture sequencing data that we recently generated. Our target capture method has a unique feature, i.e., it covers not only exonic regions but also conserved non-coding sequences (CNSs) among vertebrates; therefore, it has the potential to detect regulatory mutations. We identified a deletion of approximately 50 kb in length approximately 50 kb upstream of Ednrb. A comparative analysis with the epigenomic data in the corresponding region in humans and mice showed that one of the CNSs might be an enhancer. Further comparison with Hi-C data, which provide information about chromosome conformation, indicated that the putative enhancer is spatially close to the promoter of Ednrb, suggesting that it acts as an enhancer of Ednrb.

Conclusions

These in silico data analyses strongly suggest that the identified deletion in the intergenic region upstream of Ednrb, which might contain a melanocyte-specific enhancer, is the mutation causative of the head spot phenotype in the KFRS4/Kyo rat strain.

Electronic supplementary material

The online version of this article (doi:10.1186/s12863-017-0497-3) contains supplementary material, which is available to authorized users.

Control of tissue size and development by a regulatory element in the yorkie 3'UTR

Regulation of the Hippo pathway via phosphorylation of Yorkie (Yki), the Drosophila homolog of human Yes-associated protein 1, is conserved from Drosophila to humans. Overexpression of a non-phosphorylatable form of Yki induces severe overgrowth in adult fly eyes. Here, we show that yki mRNA associates with microsomal fractions and forms foci that partially colocalize to processing bodies in the vicinity of endoplasmic reticulum. This localization is dependent on a stem-loop (SL) structure in the 3' untranslated region of yki. Surprisingly, expression of SL deleted yki in eye imaginal discs also results in severe overgrowth phenotypes. When the structure of the SL is disrupted, Yki protein levels increase without a significant effect on RNA levels. When the SL is completely removed, protein levels drastically increase, but in this case, due to increased RNA stability. In the latter case, we show that the increased RNA accumulation is due to removal of a putative miR-8 seed sequence in the SL. These data demonstrate the function of two novel regulatory mechanisms, both controlled by the yki SL element, that are essential for proper Hippo pathway mediated growth regulation.

Software updates in the Illumina HiSeq platform affect whole-genome bisulfite sequencing

BACKGROUND

Methylation of cytosine in genomic DNA is a well-characterized epigenetic modification involved in many cellular processes and diseases. Whole-genome bisulfite sequencing (WGBS), such as MethylC-seq and post-bisulfite adaptor tagging sequencing (PBAT-seq), uses the power of high-throughput DNA sequencers and provides genome-wide DNA methylation profiles at single-base resolution. However, the accuracy and consistency of WGBS outputs in relation to the operating conditions of high-throughput sequencers have not been explored.

RESULTS

We have used the Illumina HiSeq platform for our PBAT-based WGBS, and found that different versions of HiSeq Control Software (HCS) and Real-Time Analysis (RTA) installed on the system provided different global CpG methylation levels (approximately 5% overall difference) for the same libraries. This problem was reproduced multiple times with different WGBS libraries and likely to be associated with the low sequence diversity of bisulfite-converted DNA. We found that HCS was the major determinant in the observed differences. To determine which version of HCS is most suitable for WGBS, we used substrates with predetermined CpG methylation levels, and found that HCS v2.0.5 is the best among the examined versions. HCS v2.0.12 showed the poorest performance and provided artificially lower CpG methylation levels when 5-methylcytosine is read as guanine (first read of PBAT-seq and second read of MethylC-seq). In addition, paired-end sequencing of low diversity libraries using HCS v2.2.38 or the latest HCS v2.2.58 was greatly affected by cluster densities.

CONCLUSIONS

Software updates in the Illumina HiSeq platform can affect the outputs from low-diversity sequencing libraries such as WGBS libraries. More recent versions are not necessarily the better, and HCS v2.0.5 is currently the best for WGBS among the examined HCS versions. Thus, together with other experimental conditions, special care has to be taken on this point when CpG methylation levels are to be compared between different samples by WGBS.

Alterations in Fetal Leydig Cell Gene Expression during Fetal and Adult Development

Fetal Leydig cells (FLCs) and adult Leydig cells (ALCs) develop in the mammalian prenatal and postnatal testes, respectively. In mice, FLCs emerge in the interstitial space of the testis as early as embryonic day 12.5 and thereafter increase in number during the fetal stage. We previously established a transgenic mouse line in which FLCs are labeled with EGFP and demonstrated that the EGFP-labeled FLCs were present even in adult testes. However, the characteristics of FLCs during postnatal stages remained unclear. In the present study, a comparison of the transcriptomes of FLCs from prenatal and postnatal testes and of ALCs from adult testes revealed that FLCs gradually alter their characteristics across developmental stages and come to roughly resemble ALCs. Many cholesterogenic genes simultaneously expressed a unique alternation pattern, while many oxidative phosphorylation and β-oxidation (both mitochondrial functions) genes showed a different unique pattern. These metabolic gene expression alterations might be triggered by milieu changes, such as nutrient and oxygen supply, from the prenatal to the postnatal period.

Application of target capture sequencing of exons and conserved non-coding sequences to 20 inbred rat strains

We report sequence data obtained by our recently devised target capture method TargetEC applied to 20 inbred rat strains. This method encompasses not only all annotated exons but also highly conserved non-coding sequences shared among vertebrates. The total length of the target regions covers 146.8 Mb. On an average, we obtained 31.7 × depth of target coverage and identified 154,330 SNVs and 24,368 INDELs for each strain. This corresponds to 470,037 unique SNVs and 68,652 unique INDELs among the 20 strains. The sequence data can be accessed at DDBJ/EMBL/GenBank under accession number PRJDB4648, and the identified variants have been deposited at http://bioinfo.sls.kyushu-u.ac.jp/rat_target_capture/20_strains.vcf.gz.

Linkage disequilibrium analysis of allelic heterogeneity in DNA methylation

Heterogeneity of DNA methylation status among alleles is observed in various cell types and is involved in epigenetic gene regulation and cancer biology. However, the individual methylation profile within each allele has not yet been examined at the whole-genome level. In the present study, we applied linkage disequilibrium analysis to the DNA methylation data obtained from whole-genome bisulfite sequencing studies in mouse germline and other types of cells. We found that the methylation status of 2 consecutive CpG sites showed deviation from equilibrium frequency toward concordant linkage (both methylated or both unmethylated) in germline cells. In the imprinting loci where methylation of constituent alleles is known, our analysis detected the deviation toward the concordant linkage as expected. In addition, we applied this analysis to the transitional zone between methylated and unmethylated regions and to the cells undergoing epigenetic reprogramming. In both cases, deviation to the concordant-linked alleles was conspicuous, indicating that the methylation pattern is not random but rather concordant within each allele. These results will provide the key to understanding the mechanism underlying allelic heterogeneity.

Allele-Specific Methylome and Transcriptome Analysis Reveals Widespread Imprinting in the Human Placenta

DNA methylation is globally reprogrammed after fertilization, and as a result, the parental genomes have similar DNA-methylation profiles after implantation except at the germline differentially methylated regions (gDMRs). We and others have previously shown that human blastocysts might contain thousands of transient maternally methylated gDMRs (transient mDMRs), whose maternal methylation is lost in embryonic tissues after implantation. In this study, we performed genome-wide allelic DNA methylation analyses of purified trophoblast cells from human placentas and, surprisingly, found that more than one-quarter of the transient-in-embryo mDMRs maintained their maternally biased DNA methylation. RNA-sequencing-based allelic expression analyses revealed that some of the placenta-specific mDMRs were associated with expression of imprinted genes (e.g., TIGAR, SLC4A7, PROSER2-AS1, and KLHDC10), and three imprinted gene clusters were identified. This approach also identified some X-linked gDMRs. Comparisons of the data with those from other mammals revealed that genomic imprinting in the placenta is highly variable. These findings highlight the incomplete erasure of germline DNA methylation in the human placenta; understanding this erasure is important for understanding normal placental development and the pathogenesis of developmental disorders with imprinting effects.

Design and application of a target capture sequencing of exons and conserved non-coding sequences for the rat

Background

Target capture sequencing is an efficient approach to directly identify the causative mutations of genetic disorders. To apply this strategy to laboratory rats exhibiting various phenotypes, we developed a novel target capture probe set, TargetEC (target capture for exons and conserved non-coding sequences), which can identify mutations not only in exonic regions but also in conserved non-coding sequences and thus can detect regulatory mutations.

Results

TargetEC covers 1,078,129 regions spanning 146.8 Mb of the genome. We applied TargetEC to four inbred rat strains (WTC/Kyo, WTC-swh/Kyo, PVG/Seac, and KFRS4/Kyo) maintained by the National BioResource Project for the Rat in Japan, and successfully identified mutations associated with these phenotypes, including one mutation detected in a conserved non-coding sequence.

Conclusions

The method developed in this study can be used to efficiently identify regulatory mutations, which cannot be detected using conventional exome sequencing, and will help to deepen our understanding of the relationships between regulatory mutations and associated phenotypes.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-016-2975-9) contains supplementary material, which is available to authorized users.