Whole-genome sequence of the bovine blood fluke Schistosoma bovis supports interspecific hybridization with S. haematobium.

Mesenteric infection by the parasitic blood fluke Schistosoma bovis is a common veterinary problem in Africa and the Middle East and occasionally in the Mediterranean Region. The species also has the ability to form interspecific hybrids with the human parasite S. haematobium with natural hybridisation observed in West Africa, presenting possible zoonotic transmission. Additionally, this exchange of alleles between species may dramatically influence disease dynamics and parasite evolution. We have generated a 374 Mb assembly of the S. bovis genome using Illumina and PacBio-based technologies. Despite infecting different hosts and organs, the genome sequences of S. bovis and S. haematobium appeared strikingly similar with 97% sequence identity. The two species share 98% of protein-coding genes, with an average sequence identity of 97.3% at the amino acid level. Genome comparison identified large continuous parts of the genome (up to several 100 kb) showing almost 100% sequence identity between S. bovis and S. haematobium. It is unlikely that this is a result of genome conservation and provides further evidence of natural interspecific hybridization between S. bovis and S. haematobium. Our results suggest that foreign DNA obtained by interspecific hybridization was maintained in the population through multiple meiosis cycles and that hybrids were sexually reproductive, producing viable offspring. The S. bovis genome assembly forms a highly valuable resource for studying schistosome evolution and exploring genetic regions that are associated with species-specific phenotypic traits.

Whole-genome sequencing of human malignant mesothelioma tumours and cell lines.

Pleural mesothelioma is a cancer of serosal surfaces caused by environmental exposure to asbestos. Clinical outcome remains poor and while trials of new treatments are ongoing it remains an understudied cancer. Mesothelioma cell lines can readily be grown from primary tumour and from tumour cells shed into pleural effusion with the latter representing a particularly valuable source of DNA in clinical settings, procurable without the need for additional invasive procedures. However, it is not well understood how accurately patient-derived cultured tumour cells represent the molecular characteristics of their primary tumour. We used whole-genome sequencing of primary tumour and matched cultured cells to comprehensively characterize mutations and structural alterations. Most cases had complex rearranged genomes with evidence of chromoanagenesis and rearrangements reminiscent of chromoplexy. Many of the identified driver mutations were structural, indicating that mesothelioma is often caused by structural alterations and catastrophic genomic events, rather than point mutations. Because the majority of genomic changes detected in tumours were also displayed by the genomes of cultured tumour cells, we conclude that low-passage cultured tumour cells are generally suitable for molecular characterization of mesothelioma and may be particularly useful where tissue samples with high tumour cell content are not available. However, the subclonal compositions of the cell lines did not fully recapitulate the subclonal diversity of the primary tumours. Furthermore, longitudinal acquisition of major alterations in subclonal cell populations was observed after long-term passaging. These two factors define limitations of tumour-derived cell lines as genomic substrate for clinical purposes.

Trim33 Binds and Silences a Class of Young Endogenous Retroviruses in the Mouse Testis; a Novel Component of the Arms Race between Retrotransposons and the Host Genome.

Transposable elements (TEs) have been active in the mammalian genome for millions of years and the silencing of these elements in the germline is important for the survival of the host. Mice carrying reporter transgenes can be used to model transcriptional silencing. A mutagenesis screen for modifiers of epigenetic gene silencing produced a line with a mutation in Trim33; the mutants displayed increased expression of the reporter transgene. ChIP-seq of Trim33 in testis revealed 9,109 peaks, mostly at promoters. This is the first report of ChIP-seq for Trim33 in any tissue. Comparison with ENCODE datasets showed that regions of high read density for Trim33 had high read density for histone marks associated with transcriptional activity and mapping to TE consensus sequences revealed Trim33 enrichment at RLTR10B, the LTR of one of the youngest retrotransposons in the mouse genome, MMERVK10C. We identified consensus sequences from the 266 regions at which Trim33 ChIP-seq peaks overlapped RLTR10B elements and found a match to the A-Myb DNA-binding site. We found that TRIM33 has E3 ubiquitin ligase activity for A-MYB and regulates its abundance. RNA-seq revealed that mice haploinsufficient for Trim33 had altered expression of a small group of genes in the testis and the gene with the most significant increase was found to be transcribed from an upstream RLTR10B. These studies provide the first evidence that A-Myb has a role in the actions of Trim33 and suggest a role for both A-Myb and Trim33 in the arms race between the transposon and the host. This the first report of any factor specifically regulating RLTR10B and adds to the current literature on the silencing of MMERVK10C retrotransposons. This is also the first report that A-Myb has a role in the transcription of any retrotransposon.

Identification of the CIMP-like subtype and aberrant methylation of members of the chromosomal segregation and spindle assembly pathways in esophageal adenocarcinoma.

The incidence of esophageal adenocarcinoma (EAC) has risen significantly over recent decades. Although survival has improved, cure rates remain poor, with <20% of patients surviving 5 years. This is the first study to explore methylome, transcriptome and ENCODE data to characterize the role of methylation in EAC. We investigate the genome-wide methylation profile of 250 samples including 125 EAC, 19 Barrett's esophagus (BE), 85 squamous esophagus and 21 normal stomach. Transcriptome data of 70 samples (48 EAC, 4 BE and 18 squamous esophagus) were used to identify changes in methylation associated with gene expression. BE and EAC showed similar methylation profiles, which differed from squamous tissue. Hypermethylated sites in EAC and BE were mainly located in CpG-rich promoters. A total of 18575 CpG sites associated with 5538 genes were differentially methylated, 63% of these genes showed significant correlation between methylation and mRNA expression levels. Pathways involved in tumorigenesis including cell adhesion, TGF and WNT signaling showed enrichment for genes aberrantly methylated. Genes involved in chromosomal segregation and spindle formation were aberrantly methylated. Given the recent evidence that chromothripsis may be a driver mechanism in EAC, the role of epigenetic perturbation of these pathways should be further investigated. The methylation profiles revealed two EAC subtypes, one associated with widespread CpG island hypermethylation overlapping H3K27me3 marks and binding sites of the Polycomb proteins. These subtypes were supported by an independent set of 89 esophageal cancer samples. The most hypermethylated tumors showed worse patient survival.

The recently identified modifier of murine metastable epialleles, Rearranged L-Myc Fusion, is involved in maintaining epigenetic marks at CpG island shores and enhancers.

BACKGROUND: We recently identified a novel protein, Rearranged L-myc fusion (Rlf), that is required for DNA hypomethylation and transcriptional activity at two specific regions of the genome known to be sensitive to epigenetic gene silencing. To identify other loci affected by the absence of Rlf, we have now analysed 12 whole genome bisulphite sequencing datasets across three different embryonic tissues/stages from mice wild-type or null for Rlf. RESULTS: Here we show that the absence of Rlf results in an increase in DNA methylation at thousands of elements involved in transcriptional regulation and many of the changes occur at enhancers and CpG island shores. ChIP-seq for H3K4me1, a mark generally found at regulatory elements, revealed associated changes at many of the regions that are differentially methylated in the Rlf mutants. RNA-seq showed that the numerous effects of the absence of Rlf on the epigenome are associated with relatively subtle effects on the mRNA population. In vitro studies suggest that Rlf's zinc fingers have the capacity to bind DNA and that the protein interacts with other known epigenetic modifiers. CONCLUSION: This study provides the first evidence that the epigenetic modifier Rlf is involved in the maintenance of DNA methylation at enhancers and CGI shores across the genome.

Identification of novel hypomorphic and null mutations in Klf1 derived from a genetic screen for modifiers of α-globin transgene variegation.

Position-effect variegation of transgene expression is sensitive to the chromatin state. We previously reported a forward genetic screen in mice carrying a variegated α-globin GFP transgene to find novel genes encoding epigenetic regulators. We named the phenovariant strains "Mommes" for modifiers of murine metastable epialleles. Here we report positional cloning of mutations in two Momme strains which result in suppression of variegation. Both strains harbour point mutations in the erythroid transcription factor, Klf1. One (D11) generates a stop codon in the zinc finger domain and a homozygous null phenotype. The other (D45) generates an amino acid transversion (H350R) within a conserved linker between zinc fingers two and three. Homozygous MommeD45 mice have chronic microcytic anaemia which models the phenotype in a recently described family. This is the first genetic evidence that the linkers between the zinc fingers of transcription factors have a function beyond that of a simple spacer.

SMCHD1 has separable roles in chromatin architecture and gene silencing that could be targeted in disease.

The interplay between 3D chromatin architecture and gene silencing is incompletely understood. Here, we report a novel point mutation in the non-canonical SMC protein SMCHD1 that enhances its silencing capacity at endogenous developmental targets. Moreover, it also results in enhanced silencing at the facioscapulohumeral muscular dystrophy associated macrosatellite-array, D4Z4, resulting in enhanced repression of DUX4 encoded by this repeat. Heightened SMCHD1 silencing perturbs developmental Hox gene activation, causing a homeotic transformation in mice. Paradoxically, the mutant SMCHD1 appears to enhance insulation against other epigenetic regulators, including PRC2 and CTCF, while depleting long range chromatin interactions akin to what is observed in the absence of SMCHD1. These data suggest that SMCHD1's role in long range chromatin interactions is not directly linked to gene silencing or insulating the chromatin, refining the model for how the different levels of SMCHD1-mediated chromatin regulation interact to bring about gene silencing in normal development and disease.

The characterisation of piRNA-related 19mers in the mouse.

BACKGROUND: Piwi interacting RNA, or piRNA, is a class of small RNA almost exclusively expressed in the germline where they serve essential roles in retrotransposon silencing. There are two types, primary and secondary piRNA, and the latter is a product of enzymatic cleavage of retrotransposons' transcripts directed by the former. Recently, a new class of 19nt long RNA was discovered that is specific to testis and appears to be linked to secondary piRNA biogenesis. RESULTS: We locate clusters of the testis-specific 19mers, which we call piRNA-related 19mers (pr19RNA), and characterise the transcripts from which they are derived. Most pr19RNA clusters were associated with retrotransposons and unannotated antisense transcripts overlapping piRNA clusters. At these loci the abundance of 19mers was found to be greater than that of secondary piRNAs. CONCLUSION: We find that pr19RNAs are distinguished from other RNA populations by their length and flanking sequence, allowing their identification without requiring overlapping piRNAs. Using such sequence features allows identification of the source transcripts, and we suggest that these likely represent the substrates of primary piRNA-guided RNA cleavage events. While pr19RNAs appear not to bind directly to Miwi or Mili, their abundance relative to secondary piRNAs, in combination with their precise length, suggests they may be more than by-products of secondary piRNA biogenesis.

The role of MORC3 in silencing transposable elements in mouse embryonic stem cells.

BACKGROUND: Microrchidia proteins (MORCs) are involved in epigenetic gene silencing in a variety of eukaryotic organisms. Deletion of MORCs result in several developmental abnormalities and their dysregulation has been implicated in developmental disease and multiple cancers. Specifically, mammalian MORC3 mutations are associated with immune system defects and human cancers such as bladder, uterine, stomach, lung, and diffuse large B cell lymphomas. While previous studies have shown that MORC3 binds to H3K4me3 in vitro and overlaps with H3K4me3 ChIP-seq peaks in mouse embryonic stem cells, the mechanism by which MORC3 regulates gene expression is unknown. RESULTS: In this study, we identified that mutation in Morc3 results in a suppressor of variegation phenotype in a Modifiers of murine metastable epialleles Dominant (MommeD) screen. We also find that MORC3 functions as an epigenetic silencer of transposable elements (TEs) in mouse embryonic stem cells (mESCs). Loss of Morc3 results in upregulation of TEs, specifically those belonging to the LTR class of retrotransposons also referred to as endogenous retroviruses (ERVs). Using ChIP-seq we found that MORC3, in addition to its known localization at H3K4me3 sites, also binds to ERVs, suggesting a direct role in regulating their expression. Previous studies have shown that these ERVs are marked by the repressive histone mark H3K9me3 which plays a key role in their silencing. However, we found that levels of H3K9me3 showed only minor losses in Morc3 mutant mES cells. Instead, we found that loss of Morc3 resulted in increased chromatin accessibility at ERVs as measured by ATAC-seq. CONCLUSIONS: Our results reveal MORC3 as a novel regulator of ERV silencing in mouse embryonic stem cells. The relatively minor changes of H3K9me3 in the Morc3 mutant suggests that MORC3 acts mainly downstream of, or in a parallel pathway with, the TRIM28/SETDB1 complex that deposits H3K9me3 at these loci. The increased chromatin accessibility of ERVs in the Morc3 mutant suggests that MORC3 may act at the level of chromatin compaction to effect TE silencing.

Radiomics Biomarkers Correlate with CD8 Expression and Predict Immune Signatures in Melanoma Patients.

Treatment for metastatic melanoma includes targeted and/or immunotherapy. Although many patients respond, only a subset has complete response. As late-stage patients often have multiple tumors in difficult access sites, non-invasive techniques are necessary for the development of predictive/prognostic biomarkers. PET/CT scans from 52 patients with stage III/IV melanoma were assessed and CT image parameters were evaluated as prognostic biomarkers. Analysis indicated patients with high standard deviation or high mean of positive pixels (MPP) had worse progression-free survival (P = 0.00047 and P = 0.0014, respectively) and worse overall survival (P = 0.0223 and P = 0.0465, respectively). Whole-exome sequencing showed high MPP was associated with BRAF mutation status (P = 0.0389). RNA-sequencing indicated patients with immune "cold" signatures had worse survival, which was associated with CT biomarker, MPP4 (P = 0.0284). Multiplex immunofluorescence confirmed a correlation between CD8 expression and image biomarkers (P = 0.0028). IMPLICATIONS: CT parameters have the potential to be cost-effective biomarkers of survival in melanoma, and reflect the tumor immune-microenvironment. VISUAL OVERVIEW: http://mcr.aacrjournals.org/content/molcanres/19/6/950/F1.large.jpg.

Targeting an adenosine-mediated "don't eat me signal" augments anti-lymphoma immunity by anti-CD20 monoclonal antibody.

A growing body of evidence suggests that macrophage immune checkpoint molecules are potential targets in the era of cancer immunotherapy. Here we showed that extracellular adenosine, an abundant metabolite in the tumor microenvironment, critically impedes the therapeutic efficacy of anti-CD20 monoclonal antibodies (mAbs) against B-cell lymphoma. Using a syngeneic B-cell lymphoma model, we showed that host deficiency of adenosine 2A receptor (A2AR), but not A2BR, remarkably improved lymphoma control by anti-CD20 mAb therapy. Conditional deletion of A2AR in myeloid cells, and to a lesser extent in NK cells, augmented therapeutic efficacy of anti-CD20 mAb. Indeed, adenosine signaling impaired antibody-mediated cellular phagocytosis (ADCP) by macrophages and limited the generation of anti-lymphoma CD8+ T cells. Pharmacological inhibition of A2AR overcame the adenosine-mediated negative regulation of ADCP by rituximab in a xeno-transplanted lymphoma model. Moreover, aberrant overexpression of CD39, an apical ecto-enzyme for adenosine generation, showed a negative impact on prognosis in patients with diffuse large B-cell lymphoma, as well as on preclinical efficacy of rituximab. Together, adenosine acts as a "don't eat me signal", and may be a potential target to harness anti-lymphoma immunity.

Pathogenic germline variants are associated with poor survival in stage III/IV melanoma patients.

Patients with late stage resected cutaneous melanoma have poor overall survival (OS) and experience irreversible adverse events from systemic therapy. There is a clinical need to identify biomarkers to predict outcome. Performing germline/tumour whole-exome sequencing of 44 stage III/IV melanoma patients we identified pathogenic germline mutations in CDKN2A, CDK4, ATM, POLH, MRE11A, RECQL4 and XPC, affecting 7/44 patients. These mutations were associated with poor OS (p = 0.0082). We confirmed our findings in The Cancer Genome Atlas (TCGA) human skin cutaneous melanoma cohort where we identified pathogenic variants in 40/455 patients (p = 0.0203). Combining these cohorts (n = 499) further strengthened these findings showing germline carriers had worse OS (p = 0.0009). Additionally, we determined whether tumour mutation burden (TMB) or BRAF status were prognostic markers of survival. Low TMB rate (< 20 Mut/Mb; p = 0.0034) and BRAF p.V600 mutation (p = 0.0355) were associated with worse progression-free survival. Combining these biomarkers indicated that V600 mutant patients had significantly lower TMB (p = 0.0155). This was confirmed in the TCGA (n = 443, p = 0.0007). Integrative analysis showed germline mutation status conferred the highest risk (HR 5.2, 95% CI 1.72-15.7). Stage IV (HR 2.5, 0.74-8.6) and low TMB (HR 2.3, 0.57-9.4) were similar, whereas BRAF V600 status was the weakest prognostic biomarker (HR 1.5, 95% CI 0.44-5.2).

Chromosome arm aneuploidies shape tumour evolution and drug response.

Chromosome arm aneuploidies (CAAs) are pervasive in cancers. However, how they affect cancer development, prognosis and treatment remains largely unknown. Here, we analyse CAA profiles of 23,427 tumours, identifying aspects of tumour evolution including probable orders in which CAAs occur and CAAs predicting tissue-specific metastasis. Both haematological and solid cancers initially gain chromosome arms, while only solid cancers subsequently preferentially lose multiple arms. 72 CAAs and 88 synergistically co-occurring CAA pairs multivariately predict good or poor survival for 58% of 6977 patients, with negligible impact of whole-genome doubling. Additionally, machine learning identifies 31 CAAs that robustly alter response to 56 chemotherapeutic drugs across cell lines representing 17 cancer types. We also uncover 1024 potential synthetic lethal pharmacogenomic interactions. Notably, in predicting drug response, CAAs substantially outperform  mutations and focal deletions/amplifications combined. Thus, CAAs predict cancer prognosis, shape tumour evolution, metastasis and drug response, and may advance precision oncology.

Whole-genome sequence of the oriental lung fluke Paragonimus westermani.

Background: Foodborne infections caused by lung flukes of the genus Paragonimus are a significant and widespread public health problem in tropical areas. Approximately 50 Paragonimus species have been reported to infect animals and humans, but Paragonimus westermani is responsible for the bulk of human disease. Despite their medical and economic importance, no genome sequence for any Paragonimus species is available. Results: We sequenced and assembled the genome of P. westermani, which is among the largest of the known pathogen genomes with an estimated size of 1.1 Gb. A 922.8 Mb genome assembly was generated from Illumina and Pacific Biosciences (PacBio) sequence data, covering 84% of the estimated genome size. The genome has a high proportion (45%) of repeat-derived DNA, particularly of the long interspersed element and long terminal repeat subtypes, and the expansion of these elements may explain some of the large size. We predicted 12,852 protein coding genes, showing a high level of conservation with related trematode species. The majority of proteins (80%) had homologs in the human liver fluke Opisthorchis viverrini, with an average sequence identity of 64.1%. Assembly of the P. westermani mitochondrial genome from long PacBio reads resulted in a single high-quality circularized 20.6 kb contig. The contig harbored a 6.9 kb region of non-coding repetitive DNA comprised of three distinct repeat units. Our results suggest that the region is highly polymorphic in P. westermani, possibly even within single worm isolates. Conclusions: The generated assembly represents the first Paragonimus genome sequence and will facilitate future molecular studies of this important, but neglected, parasite group.

Wiz binds active promoters and CTCF-binding sites and is required for normal behaviour in the mouse.

We previously identified Wiz in a mouse screen for epigenetic modifiers. Due to its known association with G9a/GLP, Wiz is generally considered a transcriptional repressor. Here, we provide evidence that it may also function as a transcriptional activator. Wiz levels are high in the brain, but its function and direct targets are unknown. ChIP-seq was performed in adult cerebellum and Wiz peaks were found at promoters and transcription factor CTCF binding sites. RNA-seq in Wiz mutant mice identified genes differentially regulated in adult cerebellum and embryonic brain. In embryonic brain most decreased in expression and included clustered protocadherin genes. These also decreased in adult cerebellum and showed strong Wiz ChIP-seq enrichment. Because a precise pattern of protocadherin gene expression is required for neuronal development, behavioural tests were carried out on mutant mice, revealing an anxiety-like phenotype. This is the first evidence of a role for Wiz in neural function.

Hypomethylation of ERVs in the sperm of mice haploinsufficient for the histone methyltransferase Setdb1 correlates with a paternal effect on phenotype.

The number of reports of paternal epigenetic influences on the phenotype of offspring in rodents is increasing but the molecular events involved remain unclear. Here, we show that haploinsufficiency for the histone 3 lysine 9 methyltransferase Setdb1 in the sire can influence the coat colour phenotype of wild type offspring. This effect occurs when the allele that directly drives coat colour is inherited from the dam, inferring that the effect involves an "in trans" step. The implication of this finding is that epigenetic state of the sperm can alter the expression of genes inherited on the maternally derived chromosomes. Whole genome bisulphite sequencing revealed that Setdb1 mutant mice show DNA hypomethylation at specific classes of transposable elements in the sperm. Our results identify Setdb1 as a paternal effect gene in the mouse and suggest that epigenetic inheritance may be more likely in individuals with altered levels of epigenetic modifiers.

Genetic and epigenetic variation among inbred mouse littermates: identification of inter-individual differentially methylated regions.

BACKGROUND: Phenotypic variability among inbred littermates reared in controlled environments remains poorly understood. Metastable epialleles refer to loci that intrinsically behave in this way and a few examples have been described. They display differential methylation in association with differential expression. For example, inbred mice carrying the agouti viable yellow (A (vy) ) allele show a range of coat colours associated with different DNA methylation states at the locus. The availability of next-generation sequencing, in particular whole genome sequencing of bisulphite converted DNA, allows us, for the first time, to search for metastable epialleles at base pair resolution. RESULTS: Using whole genome bisulphite sequencing of DNA from the livers of five mice from the A (vy) colony, we searched for sites at which DNA methylation differed among the mice. A small number of loci, 356, were detected and we call these inter-individual Differentially Methylated Regions, iiDMRs, 55 of which overlap with endogenous retroviral elements (ERVs). Whole genome resequencing of two mice from the colony identified very few differences and these did not occur at or near the iiDMRs. Further work suggested that the majority of ERV iiDMRs are metastable epialleles; the level of methylation was maintained in tissue from other germ layers and the level of mRNA from the neighbouring gene inversely correlated with methylation state. Most iiDMRs that were not overlapping ERV insertions occurred at tissue-specific DMRs and it cannot be ruled out that these are driven by changes in the ratio of cell types in the tissues analysed. CONCLUSIONS: Using the most thorough genome-wide profiling technologies for differentially methylated regions, we find very few intrinsically epigenetically variable regions that we term iiDMRs. The most robust of these are at retroviral elements and appear to be metastable epialleles. The non-ERV iiDMRs cannot be described as metastable epialleles at this stage but provide a novel class of variably methylated elements for further study.

An ENU mutagenesis screen identifies novel and known genes involved in epigenetic processes in the mouse.

BACKGROUND: We have used a sensitized ENU mutagenesis screen to produce mouse lines that carry mutations in genes required for epigenetic regulation. We call these lines Modifiers of murine metastable epialleles (Mommes). RESULTS: We report a basic molecular and phenotypic characterization for twenty of the Momme mouse lines, and in each case we also identify the causative mutation. Three of the lines carry a mutation in a novel epigenetic modifier, Rearranged L-myc fusion (Rlf), and one gene, Rap-interacting factor 1 (Rif1), has not previously been reported to be involved in transcriptional regulation in mammals. Many of the other lines are novel alleles of known epigenetic regulators. For two genes, Rlf and Widely-interspaced zinc finger (Wiz), we describe the first mouse mutants. All of the Momme mutants show some degree of homozygous embryonic lethality, emphasizing the importance of epigenetic processes. The penetrance of lethality is incomplete in a number of cases. Similarly ,abnormalities in phenotype seen in the heterozygous individuals of some lines occur with incomplete penetrance. CONCLUSIONS: Recent advances in sequencing enhance the power of sensitized mutagenesis screens to identify the function of previously uncharacterized factors and to discover additional functions for previously characterized proteins. The observation of incomplete penetrance of phenotypes in these inbred mutant mice, at various stages of development, is of interest. Overall, the Momme collection of mouse mutants provides a valuable resource for researchers across many disciplines.

A forward genetic screen identifies eukaryotic translation initiation factor 3, subunit H (eIF3h), as an enhancer of variegation in the mouse.

We have used a forward genetic screen to identify genes required for transgene silencing in the mouse. Previously these genes were found using candidate-based sequencing, a slow and labor-intensive process. Recently, whole-exome deep sequencing has accelerated our ability to find the causative point mutations, resulting in the discovery of novel and sometimes unexpected genes. Here we report the identification of translation initiation factor 3, subunit H (eIF3h) in two modifier of murine metastable epialleles (Mommes) lines. Mice carrying mutations in this gene have not been reported previously, and a possible involvement of eIF3h in transcription or epigenetic regulation has not been considered.