Protocol for generating high-quality genome-scale DNA methylation sequencing data from human cancer biospecimens

Aberrant DNA methylation is a universal feature of cancer. Here, we present a protocol for generating high-quality genome-scale DNA methylation sequencing data from a variety of human cancer biospecimens including immortalized cell lines, fresh-frozen surgical resections, and formalin-fixed paraffin-embedded tissues. We describe steps for DNA extraction considerations, reduced representation bisulfite sequencing, data processing and quality control, and downstream data analysis and integration. This protocol is also applicable for other human diseases and methylome profiling in other organisms. For complete details on the use and execution of this protocol, please refer to Rodger et al. (2023).1.

An epigenetic signature of advanced colorectal cancer metastasis

Colorectal cancer (CRC) is a leading cause of morbidity and mortality worldwide. The majority of CRC deaths are caused by tumor metastasis, even following treatment. There is strong evidence for epigenetic changes, such as DNA methylation, accompanying CRC metastasis and poorer patient survival. Earlier detection and a better understanding of molecular drivers for CRC metastasis are of critical clinical importance. Here, we identify a signature of advanced CRC metastasis by performing whole genome-scale DNA methylation and full transcriptome analyses of paired primary cancers and liver metastases from CRC patients. We observed striking methylation differences between primary and metastatic pairs. A subset of loci showed coordinated methylation-expression changes, suggesting these are potentially epigenetic drivers that control the expression of critical genes in the metastatic cascade. The identification of CRC epigenomic markers of metastasis has the potential to enable better outcome prediction and lead to the discovery of new therapeutic targets.

The transposable element-derived transcript of LIN28B has a placental origin and is not specific to tumours

Transposable elements (TEs) are genetic elements that have evolved as crucial regulators of human development and cancer, functioning as both genes and regulatory elements. When TEs become dysregulated in cancer cells, they can serve as alternate promoters to activate oncogenes, a process known as onco-exaptation. This study aimed to explore the expression and epigenetic regulation of onco-exaptation events in early human developmental tissues. We discovered co-expression of some TEs and oncogenes in human embryonic stem cells and first trimester and term placental tissues. Previous studies identified onco-exaptation events in various cancer types, including an AluJb SINE element-LIN28B interaction in lung cancer cells, and showed that the TE-derived LIN28B transcript is associated with poor patient prognosis in hepatocellular carcinoma. This study further characterized the AluJb-LIN28B transcript and confirmed that its expression is restricted to the placenta. Targeted DNA methylation analysis revealed differential methylation of the two LIN28B promoters between placenta and healthy somatic tissues, indicating that some TE-oncogene interactions are not cancer-specific but arise from the epigenetic reactivation of developmental TE-derived regulatory events. In conclusion, our findings provide evidence that some TE-oncogene interactions are not limited to cancer and may originate from the epigenetic reactivation of TE-derived regulatory events that are involved in early development. These insights broaden our understanding of the role of TEs in gene regulation and suggest the potential importance of targeting TEs in cancer therapy beyond their conventional use as cancer-specific markers.

Strategy for RNA-Seq Experimental Design and Data Analysis

Ribonucleic acids (RNAs) are fundamental molecules that control regulation and expression of the genome and therefore the function of a cell. Robust analysis and quantification of RNA transcripts hold critical importance in understanding cell function, altered phenotypes in different biological context, for understanding and targeting diseases. The development of RNA-sequencing (RNA-Seq) now provides opportunities to analyze the expression and function of RNA molecules at an unprecedented scale. However, the strategy for RNA-Seq experimental design and data analysis can substantially differ depending on the biological application. The design choice could also have significant impact for downstream results and interpretation of data. Here we describe key critical considerations required for RNA-Seq experimental design and also describe a step-by-step bioinformatics workflow detailing the different steps required for RNA-Seq data analysis. We believe this article will be a valuable guide for designing and analyzing RNA-Seq data to address a wide range of different biological questions.

Dynamic Epigenetic Changes during a Relapse and Recovery Cycle in Myalgic Encephalomyelitis/Chronic Fatigue Syndrome

Myalgic Encephalomyelitis/Chronic Fatigue Syndrome (ME/CFS) is a complex disease with variable severity. Patients experience frequent relapses where symptoms increase in severity, leaving them with a marked reduction in quality of life. Previous work has investigated molecular differences between ME/CFS patients and healthy controls, but not the dynamic changes specific to each individual patient. We applied precision medicine here to map genomic changes in two selected ME/CFS patients through a period that contained a relapse recovery cycle. DNA was isolated from two patients and a healthy age/gender matched control at regular intervals and captured the patient relapse in each case. Reduced representation DNA methylation sequencing profiles were obtained spanning the relapse recovery cycle. Both patients showed a significantly larger methylome variability (10–20-fold) through the period of sampling compared with the control. During the relapse, changes in the methylome profiles of the two patients were detected in regulatory-active regions of the genome that were associated, respectively, with 157 and 127 downstream genes, indicating disturbed metabolic, immune and inflammatory functions. Severe health relapses in the ME/CFS patients resulted in functionally important changes in their DNA methylomes that, while differing between the two patients, led to very similar compromised physiology. DNA methylation as a signature of disease variability in ongoing ME/CFS may have practical applications for strategies to decrease relapse frequency.

Innate immune checkpoint inhibitor resistance is associated with melanoma sub-types exhibiting invasive and de-differentiated gene expression signatures

Melanoma is a highly aggressive skin cancer, which, although highly immunogenic, frequently escapes the body’s immune defences. Immune checkpoint inhibitors (ICI), such as anti-PD1, anti-PDL1, and anti-CTLA4 antibodies lead to reactivation of immune pathways, promoting rejection of melanoma. However, the benefits of ICI therapy remain limited to a relatively small proportion of patients who do not exhibit ICI resistance. Moreover, the precise mechanisms underlying innate and acquired ICI resistance remain unclear. Here, we have investigated differences in melanoma tissues in responder and non-responder patients to anti-PD1 therapy in terms of tumour and immune cell gene-associated signatures. We performed multi-omics investigations on melanoma tumour tissues, which were collected from patients before starting treatment with anti-PD1 immune checkpoint inhibitors. Patients were subsequently categorized into responders and non-responders to anti-PD1 therapy based on RECIST criteria. Multi-omics analyses included RNA-Seq and NanoString analysis. From RNA-Seq data we carried out HLA phenotyping as well as gene enrichment analysis, pathway enrichment analysis and immune cell deconvolution studies. Consistent with previous studies, our data showed that responders to anti-PD1 therapy had higher immune scores (median immune score for responders = 0.1335, median immune score for non-responders = 0.05426, p-value = 0.01, Mann-Whitney U two-tailed exact test) compared to the non-responders. Responder melanomas were more highly enriched with a combination of CD8+ T cells, dendritic cells (p-value = 0.03) and an M1 subtype of macrophages (p-value = 0.001). In addition, melanomas from responder patients exhibited a more differentiated gene expression pattern, with high proliferative- and low invasive-associated gene expression signatures, whereas tumours from non-responders exhibited high invasive- and frequently neural crest-like cell type gene expression signatures. Our findings suggest that non-responder melanomas to anti-PD1 therapy exhibit a de-differentiated gene expression signature, associated with poorer immune cell infiltration, which establishes a gene expression pattern characteristic of innate resistance to anti-PD1 therapy. Improved understanding of tumour-intrinsic gene expression patterns associated with response to anti-PD1 therapy will help to identify predictive biomarkers of ICI response and may help to identify new targets for anticancer treatment, especially with a capacity to function as adjuvants to improve ICI outcomes.

Generating Sequencing-Based DNA Methylation Maps from Low DNA Input Samples

Reduced representation bisulfite sequencing (RRBS) is a technique used for assessing genome-wide DNA methylation patterns in eukaryotes. RRBS was introduced to focus on CpG-rich regions that are likely to be of most interest for epigenetic regulation, such as gene promoters and enhancer sequence elements (Meissner et al., Nature 454:766-770, 2008). This "reduced representation" lowers the cost of sequencing and also gives increased depth of coverage, facilitating the resolution of more subtle changes in methylation levels. Here, we describe a modified RRBS sequencing (RRBS-seq) library preparation. Our protocol is optimized for generating single base-resolution libraries when low input DNA is a concern (10-100 ng). Our protocol includes steps to optimize library preparation, such as using deparaffinization solution (when formalin-fixed material is used), and a replacement of gel size-selection with sample purification beads. The described protocol can be accomplished in 3 days and has been successfully applied to tissues or cells from different organisms, including formalin-fixed tissues, to yield robust and reproducible results.

The Mouse Papillomavirus Epigenetic Signature Is Characterised by DNA Hypermethylation after Lesion Regression

Papillomaviruses (PVs) are double-stranded DNA tumour viruses that can infect cutaneous and mucosal epidermis. Human papillomavirus (HPV) types have been linked to the causality of cutaneous squamous cell carcinoma (cSCC); however, HPV DNA is not always detected in the resultant tumour. DNA methylation is an epigenetic change that can contribute to carcinogenesis. We hypothesise that the DNA methylation pattern in cells is altered following PV infection. We tested if DNA methylation was altered by PV infection in the mouse papillomavirus (MmuPV1) model. Immunosuppressed mice were infected with MmuPV1 on cutaneous tail skin. Immunosuppression was withdrawn for some mice, causing lesions to spontaneously regress. Reduced representation bisulphite sequencing was carried out on DNA from the actively infected lesions, visibly regressed lesions, and mock-infected control mice. DNA methylation libraries were generated and analysed for differentially methylated regions throughout the genome. The presence of MmuPV1 sequences was also assessed. We identified 834 predominantly differentially hypermethylated fragments in regressed lesions, and no methylation differences in actively infected lesions. The promoter regions of genes associated with tumorigenicity, including the tumour suppressor protein DAPK1 and mismatch repair proteins MSH6 and PAPD7, were hypermethylated. Viral DNA was detected in active lesions and in some lesions that had regressed. This is the first description of the genome-wide DNA methylation landscape for active and regressed MmuPV1 lesions. We propose that the DNA hypermethylation in the regressed lesions that we report here may increase the susceptibility of cells to ultraviolet-induced cSCC.

RepExpress: A Novel Pipeline for the Quantification and Characterization of Transposable Element Expression from RNA-seq Data

Transposable elements (TEs) are key regulators of both development and disease; however, their repetitive nature presents substantial computational challenges to their analysis. Due to a lack of computational tools and suitable analysis frameworks, TE expression is often not quantified at the locus level. Therefore, we have developed RepExpress, a novel pipeline that enables locus-level TE quantification and characterization. RepExpress enables the characterization of TE expression in a genomic context, and is the first tool focusing on the identification of tissue-specific TE-derived and TE-regulated genes. RepExpress identifies expressed TEs overlapping with annotated genomic features and enables tissue-specific profiles of TE-derived genes. TEs that are expressed with no overlap with any known genomic features are characterized by the closest downstream genomic feature enabling identification of novel TE-gene regulatory relationships. RepExpress takes standard RNA-seq data as input and performs genomic alignment optimized for TEs. Our novel pipeline quantifies expression of both TEs and genes using featureCounts and Stringtie, respectively. RepExpress then filters expressed repeats and characterizes their genomic context, enabling the identification of TEs that overlap with genes, or that may be influencing gene expression. Here, we describe RepExpress, and provide a step-by-step protocol detailing its workflow. We also discuss other TE analysis tools and their applicability to addressing different biological questions. © 2021 Wiley Periodicals LLC. Basic Protocol: RepExpress workflow.

Comparison of Global DNA Methylation Patterns in Human Melanoma Tissues and Their Derivative Cell Lines

Simple Summary

Cancer cell lines are a defined population of cells, originally sourced from tumour tissue, that can be maintained in culture for an extended period of time. They are a critical laboratory-based model, and are frequently used to unravel mechanisms of cancer cell biology. In all cells, gene activity is in part regulated by DNA methylation, the epigenetic process by which methyl groups are added to DNA. In this study, we demonstrate that at a global level, DNA methylation profiles are globally well conserved, but we identify specific sites that are consistently more methylated in tumour-derived cell lines compared to the original tumour tissue. The genes associated with these common differentially methylated regions are involved in important cellular processes and are strongly enriched for epigenetic mechanisms associated with suppression of gene activity. This study provides a valuable resource for identifying false positives due to cell culture and for better interpretation of future cancer epigenetics studies.

Abstract

DNA methylation is a heritable epigenetic mark that is fundamental to mammalian development. Aberrant DNA methylation is an epigenetic hallmark of cancer cells. Cell lines are a critical in vitro model and very widely used to unravel mechanisms of cancer cell biology. However, limited data are available to assess whether DNA methylation patterns in tissues are retained when cell lines are established. Here, we provide the first genome-scale sequencing-based methylation map of metastatic melanoma tumour tissues and their derivative cell lines. We show that DNA methylation profiles are globally conserved in vitro compared to the tumour tissue of origin. However, we identify sites that are consistently hypermethylated in cell lines compared to their tumour tissue of origin. The genes associated with these common differentially methylated regions are involved in cell metabolism, cell cycle and apoptosis and are also strongly enriched for the H3K27me3 histone mark and PRC2 complex-related genes. Our data indicate that although global methylation patterns are similar between tissues and cell lines, there are site-specific epigenomic differences that could potentially impact gene expression. Our work provides a valuable resource for identifying false positives due to cell culture and for better interpretation of cancer epigenetics studies in the future.

Genome-wide DNA methylation and RNA expression differences correlate with invasiveness in melanoma cell lines

Aims & objectives: The aim of this study was to investigate the role of DNA methylation in invasiveness in melanoma cells. Materials & methods: The authors carried out genome-wide transcriptome (RNA sequencing) and reduced representation bisulfite sequencing methylome profiling between noninvasive (n = 4) and invasive melanoma cell lines (n = 5). Results: The integration of differentially expressed genes and differentially methylated fragments (DMFs) identified 12 DMFs (two in AVPI1, one in HMG20B, two in BCL3, one in NTSR1, one in SYNJ2, one in ROBO2 and four in HORMAD2) that overlapped with either differentially expressed genes (eight DMFs and six genes) or cis-targets of lncRNAs (five DMFs associated with cis-targets and four differentially expressed lncRNAs). Conclusions: DNA methylation changes are associated with a number of transcriptional differences observed in noninvasive and invasive phenotypes in melanoma.

Reawakening the Developmental Origins of Cancer Through Transposable Elements

Transposable elements (TEs) have an established role as important regulators of early human development, functioning as tissue-specific genes and regulatory elements. Functional TEs are highly active during early development, and interact with important developmental genes, some of which also function as oncogenes. Dedifferentiation is a hallmark of cancer, and is characterized by genetic and epigenetic changes that enable proliferation, self-renewal and a metabolism reminiscent of embryonic stem cells. There is also compelling evidence suggesting that the path to dedifferentiation in cancer can contribute to invasion and metastasis. TEs are frequently expressed in cancer, and recent work has identified a newly proposed mechanism involving extensive recruitment of TE-derived promoters to drive expression of oncogenes and subsequently promote oncogenesis—a process termed onco-exaptation. However, the mechanism by which this phenomenon occurs, and the extent to which it contributes to oncogenesis remains unknown. Initial hypotheses have proposed that onco-exaptation events are cancer-specific and arise randomly due to the dysregulated and hypomethylated state of cancer cells and abundance of TEs across the genome. However, we suspect that exaptation-like events may not just arise due to chance activation of novel regulatory relationships as proposed previously, but as a result of the reestablishment of early developmental regulatory relationships. Dedifferentiation in cancer is well-documented, along with expression of TEs. The known interactions between TEs and pluripotency factors such as NANOG and OCTt4 during early development, along with the expression of some placental-specific TE-derived transcripts in cancer support a possible link between TEs and dedifferentiation of tumor cells. Thus, we hypothesize that onco-exaptation events can be associated with the epigenetic reawakening of early developmental TEs to regulate expression of oncogenes and promote oncogenesis. We also suspect that activation of these early developmental regulatory TEs may promote dedifferentiation, although at this stage it is hard to predict whether TE activation is one of the initial drivers of dedifferentiation. We expect that developmental TE activation occurs as a result of the establishment of an epigenetic landscape in cancer that resembles that of early development and that developmental TE activation may also enable cancers to exploit early developmental pathways, repurposing them to promote malignancy.

Extensive Inter-Cyst DNA Methylation Variation in Autosomal Dominant Polycystic Kidney Disease Revealed by Genome Scale Sequencing

Autosomal dominant polycystic kidney disease (ADPKD) is a heritable disease characterized by bilateral renal enlargement due to the growth of cysts throughout the kidneys. Inheritance of a disease-causing mutation is required to develop ADPKD, which results in end-stage kidney disease and is associated with a high morbidity. The pathology underlying cyst formation is not well understood. To address this, we have previously shown the global methylome is altered in ADPKD tissue, suggesting a role of DNA methylation in disease-state renal tissue. As cysts are believed to arise independently, we hypothesize that DNA methylation changes vary accordingly. Here we further investigate the role of DNA methylation within independent cysts to characterize key intra-individual changes. We demonstrate that fragments within CpG islands and gene bodies harbor the greatest amount of variation across the ADPKD kidney, while intergenic fragments are comparatively stable. A proportion of variably methylated genes were also differentially methylated in ADPKD tissue. Our data provide evidence that individual molecular mechanisms are operating in the development of each cyst.

Characterisation of DNA methylation changes in EBF3 and TBC1D16 associated with tumour progression and metastasis in multiple cancer types

Background

Characteristic DNA methylation differences have been identified between primary and metastatic melanomas at EBF3 and/or TBC1D16 gene loci. To further evaluate whether these epigenetic changes may act more generally as drivers of tumour onset and metastasis, we have investigated DNA methylation changes involving EBF3 and TBC1D16 in additional publicly available data of multiple different tumour types.

Results

Promoter hypermethylation and gene body hypomethylation of EBF3 were observed in a number of metastatic tumour types, when compared to normal or primary tumour tissues, as well as in tumour vs normal tissues and in a colorectal primary/metastasis pair, although not all tumour samples or primary/metastasis cancer pairs exhibited altered patterns of EBF3 methylation. In addition, hypomethylation of TBC1D16 was observed in multiple tumours, including a breast cancer primary/metastasis pair, and to a lesser degree in melanoma, although again not all tumours or cancer primary/metastasis pairs exhibited altered patterns of methylation.

Conclusions

These findings suggest characteristic DNA methylation changes in EBF3 and TBC1D16 are relatively common tumour-associated epigenetic events in multiple tumour types, which is consistent with a potential role as more general drivers of tumour progression.

Abstract 826: Global DNA methylation levels regulate PD-L1 expression in melanoma

The programmed death-ligand 1 (PD-L1) receptor is an important immune checkpoint and is often upregulated in cancer cells to allow immune evasion. In melanoma, the patients with PD-L1 expression and absense of tumour infiltrating lymphocytes (TILs) (i.e. “constitutive PD-L1 or PD-L1CON”) show worse response rates and prognosis than patients with PD-L1 expression and the presence of TILs (i.e. “inducible PD-L1 or PD-L1IND”). However, how PD-L1 expression is regulated in melanoma cells remains elusive. Understanding the mechanisms of how PD-L1 is regulated is important for predicting responses for anti-PD-L1 treatment and for developing new combinatorial therapies. We hypothesised that epigenetic state regulates constitutive and inducible PD-L1 expression in melanoma.

To address the hypothesis, we have generated whole-genome scale DNA methylomes (using reduced representation bisulfite sequencing) and transcriptomes (RNA-Seq) for patient derived melanoma cell lines (in PD-L1IND and PD-L1CON groups). We discovered extensive global hypomethylation in the constitutive lines, particularly pronounced in intergenic repeat regions and gene bodies. A high proportion of hypomethylated regions exhibited dichotomous methylation patterns indicating a common regulatory mechanism between the PD-L1IND and PD-L1CON lines. RNA-Sequening data indicated that the hypomethylated state of the PD-L1CON cells was correlated with higher upregulation of the differentially expressed genes at a global-scale and the upregulated genes were associated with cancer hallmark properties. The upregulated genes exhibited signatures of viral mimicry and cytosolic sensing of dsRNA similar to what has been observed after DNA methyltransferase inhibitor (DNMTi) treatment in cancers. DNMTi treatment increased PD-L1 transcription in the PD-L1IND cell lines. PD-L1CON samples showed greater resistance to DNMTi compared to PD-L1 negative samples. This was observed at protein-coding genes, transposable elements and long non-coding RNAs. Moreover, genes involved in the innate immune pathway was strongly increased upon DNMTi mediated hypomethylation in the PD-L1IND cell lines (FDR-adjusted P-value = 0.02), however there was no expression change in the PD-L1 positive cells (Adjusted P-value = 1.0). Furthermore, we found that the IRF1 transcription factor, which plays a critical role in the innate immune pathway and also binds to the PD-L1 promoter, is significantly upregulated upon DNMTi in the PD-L1 negative cell lines (Adjusted P-value = 0.01) but not in the PD-L1 positive cells (Adjusted P-value = 0.9). These results show that global hypomethylation levels regulate PD-L1 expression in melanoma. We believe these results are the first to show that DNA methylation levels play a role in regulating PD-L1 on melanoma and suggest they may have important implications for combined treatments targeting methylation (DNMTi) and PD1/PD-L1 (anti-PD1 antibodies).

Citation Format: Aniruddha Chatterjee, Antonio Ahn, Euan J. Rodger, Peter A. Stockwell, Matthew Parry, Jyoti Motwani, Stuart J. Gallagher, Elena Shklovskaya, Jessamy Tiffen, Peter Hersey, Michael R. Eccles. Global DNA methylation levels regulate PD-L1 expression in melanoma [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2019; 2019 Mar 29-Apr 3; Atlanta, GA. Philadelphia (PA): AACR; Cancer Res 2019;79(13 Suppl):Abstract nr 826.

DNA methylation and gene expression alterations in zebrafish early-life stages exposed to the antibacterial agent triclosan

There is increasing evidence that toxicant exposure can alter DNA methylation profile, one of the main epigenetic mechanisms, particularly during embryogenesis when DNA methylation patterns are being established. In order to investigate the effects of the antibacterial agent Triclosan on DNA methylation and its correlation with gene expression, zebrafish embryos were exposed during 7 days post-fertilization (starting at maximum 8-cells stage) to 50 and 100 μg/l, two conditions for which increased sensitivity and acclimation have been respectively reported. Although global DNA methylation was not significantly affected, a total of 171 differentially methylated fragments were identified by Reduced Representation Bisulfite Sequencing. The majority of these fragments were found between the two exposed groups, reflecting dose-dependant specific responses. Gene ontology analysis revealed that pathways involved in TGF-β signaling were enriched in larvae exposed to 50 μg/l, while de novo pyrimidine biosynthesis functions were overrepresented in fish exposed to 100 μg/l. In addition, gene expression analysis revealed a positive correlation between mRNA levels and DNA methylation patterns in introns, together with significant alterations of the transcription of genes involved in nervous system development, transcriptional factors and histone methyltransferases. Overall this work provides evidence that Triclosan alters DNA methylation in zebrafish exposed during embryogenesis as well as related genes expression and proposes concentration specific modes of action. Further studies will investigate the possible long-term consequences of these alterations, i.e. latent defects associated with developmental exposure and transgenerational effects, and the possible implications in terms of fitness and adaptation to environmental pollutants.

Marked Global DNA Hypomethylation Is Associated with Constitutive PD-L1 Expression in Melanoma

Constitutive expression of the immune checkpoint, PD-L1, inhibits anti-tumor immune responses in cancer, although the factors involved in PD-L1 regulation are poorly understood. Here we show that loss of global DNA methylation, particularly in intergenic regions and repeat elements, is associated with constitutive (PD-L1_CON), versus inducible (PD-L1_IND), PD-L1 expression in melanoma cell lines. We further show this is accompanied by transcriptomic up-regulation. De novo epigenetic regulators (e.g., DNMT3A) are strongly correlated with PD-L1 expression and methylome status. Accordingly, decitabine-mediated inhibition of global methylation in melanoma cells leads to increased PD-L1 expression. Moreover, viral mimicry and immune response genes are highly expressed in lymphocyte-negative plus PD-L1-positive melanomas, versus PD-L1-negative melanomas in The Cancer Genome Atlas (TCGA). In summary, using integrated genomic analysis we identified that global DNA methylation influences PD-L1 expression in melanoma, and hence melanoma's ability to evade anti-tumor immune responses. These results have implications for combining epigenetic therapy with immunotherapy.

A streamlined method for analysing genome-wide DNA methylation patterns from low amounts of FFPE DNA

Background

Formalin fixed paraffin embedded (FFPE) tumor samples are a major source of DNA from patients in cancer research. However, FFPE is a challenging material to work with due to macromolecular fragmentation and nucleic acid crosslinking. FFPE tissue particularly possesses challenges for methylation analysis and for preparing sequencing-based libraries relying on bisulfite conversion. Successful bisulfite conversion is a key requirement for sequencing-based methylation analysis.

Methods

Here we describe a complete and streamlined workflow for preparing next generation sequencing libraries for methylation analysis from FFPE tissues. This includes, counting cells from FFPE blocks and extracting DNA from FFPE slides, testing bisulfite conversion efficiency with a polymerase chain reaction (PCR) based test, preparing reduced representation bisulfite sequencing libraries and massively parallel sequencing.

Results

The main features and advantages of this protocol are:

An optimized method for extracting good quality DNA from FFPE tissues.

An efficient bisulfite conversion and next generation sequencing library preparation protocol that uses 50 ng DNA from FFPE tissue.

Incorporation of a PCR-based test to assess bisulfite conversion efficiency prior to sequencing.

Conclusions

We provide a complete workflow and an integrated protocol for performing DNA methylation analysis at the genome-scale and we believe this will facilitate clinical epigenetic research that involves the use of FFPE tissue.

Electronic supplementary material

The online version of this article (10.1186/s12920-017-0290-1) contains supplementary material, which is available to authorized users.

Comparative assessment of DNA methylation patterns between reduced representation bisulfite sequencing and Sequenom EpiTyper methylation analysis

AIM

Validation of sequencing-based DNA methylation data is an important step for meaningful translation of findings. However, there has been limited assessment of different platforms to validate methylation data from next generation sequencing.

METHODS

We performed a comparative methylation analysis between the genome-wide platform of reduced representation bisulfite sequencing with a targeted, Sequenom EpiTyper platform (four genes were analyzed in 15 cell lines covering 52 CpG sites).

RESULTS

We show that the accuracy of validation substantially improves if results from multiple and adjacent CpG sites are combined rather than at single CpG sites. We demonstrate increased read number improves accuracy of reduced representation bisulfite sequencing results. Further, by using series of replicates, we document variation in samples analyzed by Sequenom EpiTyper, indicating the importance of including replicates to increase precision.

CONCLUSION

The results reveal potential sources of bias and provide a guideline for refining study design for DNA methylation analysis.

Genome-wide methylation sequencing of paired primary and metastatic cell lines identifies common DNA methylation changes and a role for EBF3 as a candidate epigenetic driver of melanoma metastasis

Epigenetic alterations are increasingly implicated in metastasis, whereas very few genetic mutations have been identified as authentic drivers of cancer metastasis. Yet, to date, few studies have identified metastasis-related epigenetic drivers, in part because a framework for identifying driver epigenetic changes in metastasis has not been established. Using reduced representation bisulfite sequencing (RRBS), we mapped genome-wide DNA methylation patterns in three cutaneous primary and metastatic melanoma cell line pairs to identify metastasis-related epigenetic drivers. Globally, metastatic melanoma cell lines were hypomethylated compared to the matched primary melanoma cell lines. Using whole genome RRBS we identified 75 shared (10 hyper- and 65 hypomethylated) differentially methylated fragments (DMFs), which were associated with 68 genes showing significant methylation differences. One gene, Early B Cell Factor 3 (EBF3), exhibited promoter hypermethylation in metastatic cell lines, and was validated with bisulfite sequencing and in two publicly available independent melanoma cohorts (n = 40 and 458 melanomas, respectively). We found that hypermethylation of the EBF3 promoter was associated with increased EBF3 mRNA levels in metastatic melanomas and subsequent inhibition of DNA methylation reduced EBF3 expression. RNAi-mediated knockdown of EBF3 mRNA levels decreased proliferation, migration and invasion in primary and metastatic melanoma cell lines. Overall, we have identified numerous epigenetic changes characterising metastatic melanoma cell lines, including EBF3-induced aggressive phenotypic behaviour with elevated EBF3 expression in metastatic melanoma, suggesting that EBF3 promoter hypermethylation may be a candidate epigenetic driver of metastasis.

Tools and Strategies for Analysis of Genome-Wide and Gene-Specific DNA Methylation Patterns

DNA methylation is a stable epigenetic mechanism that has important roles in the normal function of a cell and therefore also in disease etiology. Accurate measurements of normal and altered DNA methylation patterns are important to understand its role in regulating gene expression and cell phenotype. Remarkable progress has been made over the last decade in developing methodologies to investigate DNA methylation. The availability of next-generation sequencing has enabled the profiling of methylation marks at an unprecedented scale. Several methods that were previously used to profile locus-specific methylation have now been upgraded to a genome-wide scale using high-throughput sequencing or array platforms. However, because there are so many techniques available, researchers are faced with the challenge of assessing the potential merits or limitations of each technique and selecting the appropriate method for their analysis. In this review we discuss the strengths and weaknesses of genome-wide and gene-specific analysis tools for interrogating DNA methylation. We particularly focus on the design and analysis strategies involved. This review will provide a guideline for selecting the appropriate methods and tools for large-scale and locus-specific DNA methylation analysis.

Generating Multiple Base-Resolution DNA Methylomes Using Reduced Representation Bisulfite Sequencing

Reduced representation bisulfite sequencing (RRBS) is an effective technique for profiling genome-wide DNA methylation patterns in eukaryotes. RRBS couples size selection, bisulfite conversion, and second-generation sequencing to enrich for CpG-dense regions of the genome. The progressive improvement of second-generation sequencing technologies and reduction in cost provided an opportunity to examine the DNA methylation patterns of multiple genomes. Here, we describe a protocol for sequencing multiple RRBS libraries in a single sequencing reaction to generate base-resolution methylomes. Furthermore, we provide a brief guideline for base-calling and data analysis of multiplexed RRBS libraries. These strategies will be useful to perform large-scale, genome-wide DNA methylation analysis.

scan_tcga tools for integrated epigenomic and transcriptomic analysis of tumor subgroups

Aim: The Cancer Genome Atlas contains multiple levels of genomic data (mutation, gene expression, DNA methylation, copy number variation) for 33 cancer types for almost 11,000 patients. However, a dearth of appropriate software tools makes it difficult for bench scientists to use these data effectively. Materials & methods: Here, we present a suite of flexible, fast and command line-based scripts that will allow retrieval and analysis of DNA methylation (tool: scan_tcga_methylation.awk), mRNA (tool: scan_tcga_mRNA.awk) and miRNA expression (tool: scan_tcga_miRNAs.awk) from cancer genome atlas network level 3 data. Results: We demonstrate the utility of these tools by analyzing DNA methylation and mRNA expression signatures of 60 frequently deregulated cancer genes and also of 30 miRNAs in primary (n = 102) and metastatic melanoma patients (n = 367). Conclusion: Our analysis illustrates the validity of the scan_tcga tools and reveals the epigenomic signatures and importance of identifying smaller patient subgroups with distinct molecular profiles.

Sex differences in DNA methylation and expression in zebrafish brain: a test of an extended 'male sex drive' hypothesis

The sex drive hypothesis predicts that stronger selection on male traits has resulted in masculinization of the genome. Here we test whether such masculinizing effects can be detected at the level of the transcriptome and methylome in the adult zebrafish brain. Although methylation is globally similar, we identified 914 specific differentially methylated CpGs (DMCs) between males and females (435 were hypermethylated and 479 were hypomethylated in males compared to females). These DMCs were prevalent in gene body, intergenic regions and CpG island shores. We also discovered 15 distinct CpG clusters with striking sex-specific DNA methylation differences. In contrast, at transcriptome level, more female-biased genes than male-biased genes were expressed, giving little support for the male sex drive hypothesis. Our study provides genome-wide methylome and transcriptome assessment and sheds light on sex-specific epigenetic patterns and in zebrafish for the first time.

Genome-scale DNA methylome and transcriptome profiling of human neutrophils

Methylation of DNA molecules is a key mechanism associated with human disease, altered gene expression and phenotype. Using reduced representation bisulphite sequencing (RRBS) technology we have analysed DNA methylation patterns in healthy individuals and identified genes showing significant inter-individual variation. Further, using whole genome transcriptome analysis (RNA-Seq) on the same individuals we showed a local and specific relationship of exon inclusion and variable DNA methylation pattern. For RRBS, 363 million, 100-bp reads were generated from 13 samples using Illumina GAII and HiSeq2000 platforms. Here we also present additional RRBS data for a female pair of monozygotic twins that was not described in our original publication. Further, We performed RNA-Seq on four of these individuals, generating 174 million, 51-bp high quality reads on an Illumina HiSeq2000 platform. The current data set could be exploited as a comprehensive resource for understanding the nature and mechanism of variable phenotypic traits and altered disease susceptibility due to variable DNA methylation and gene expression patterns in healthy individuals.

Base-resolution DNA methylation landscape of zebrafish brain and liver

Zebrafish (Danio rerio) is a vertebrate model organism that is widely used for studying a plethora of biological questions, including developmental processes, effects of external cues on phenotype, and human disease modeling. DNA methylation is an important epigenetic mechanism that contributes to gene regulation, and is prevalent in all vertebrates. Reduced representation bisulfite sequencing (RRBS) is a cost-effective technique to generate genome-wide DNA methylation maps and has been used in mammalian genomes (e.g., human, mouse and rat) but not in zebrafish. High-resolution DNA methylation data in zebrafish are limited: increased availability of such data will enable us to model and better understand the roles, causes and consequences of changes in DNA methylation.

Here we present five high-resolution DNA methylation maps for wild-type zebrafish brain (two pooled male and two pooled female methylomes) and liver. These data were generated using the RRBS technique (includes 1.43 million CpG sites of zebrafish genome) on the Illumina HiSeq platform. Alignment to the reference genome was performed using the Zv9 genome assembly.

To our knowledge, these datasets are the only RRBS datasets and base-resolution DNA methylation data available at this time for zebrafish brain and liver. These datasets could serve as a resource for future studies to document the functional role of DNA methylation in zebrafish. In addition, these datasets could be used as controls while performing analysis on treated samples.

DMAP: differential methylation analysis package for RRBS and WGBS data

MOTIVATION

The rapid development of high-throughput sequencing technologies has enabled epigeneticists to quantify DNA methylation on a massive scale. Progressive increase in sequencing capacity present challenges in terms of processing analysis and the interpretation of the large amount of data; investigating differential methylation between genome-scale data from multiple samples highlights this challenge.

RESULTS

We have developed a differential methylation analysis package (DMAP) to generate coverage-filtered reference methylomes and to identify differentially methylated regions across multiple samples from reduced representation bisulphite sequencing and whole genome bisulphite sequencing experiments. We introduce a novel fragment-based approach for investigating DNA methylation patterns for reduced representation bisulphite sequencing data. Further, DMAP provides the identity of gene and CpG features and distances to the differentially methylated regions in a format that is easily analyzed with limited bioinformatics knowledge.

AVAILABILITY AND IMPLEMENTATION

The software has been implemented in C and has been written to ensure portability between different platforms. The source code and documentation is freely available (DMAP: as compressed TAR archive folder) from http://biochem.otago.ac.nz/research/databases-software/. Two test datasets are also available for download from the Web site. Test dataset 1 contains reads from chromosome 1 of a patient and a control, which is used for comparative analysis in the current article. Test dataset 2 contains reads from a part of chromosome 21 of three disease and three control samples for testing the operation of DMAP, especially for the analysis of variance. Example commands for the analyses are included.

Mapping the zebrafish brain methylome using reduced representation bisulfite sequencing

Reduced representation bisulfite sequencing (RRBS) has been used to profile DNA methylation patterns in mammalian genomes such as human, mouse and rat. The methylome of the zebrafish, an important animal model, has not yet been characterized at base-pair resolution using RRBS. Therefore, we evaluated the technique of RRBS in this model organism by generating four single-nucleotide resolution DNA methylomes of adult zebrafish brain. We performed several simulations to show the distribution of fragments and enrichment of CpGs in different in silico reduced representation genomes of zebrafish. Four RRBS brain libraries generated 98 million sequenced reads and had higher frequencies of multiple mapping than equivalent human RRBS libraries. The zebrafish methylome indicates there is higher global DNA methylation in the zebrafish genome compared with its equivalent human methylome. This observation was confirmed by RRBS of zebrafish liver. High coverage CpG dinucleotides are enriched in CpG island shores more than in the CpG island core. We found that 45% of the mapped CpGs reside in gene bodies, and 7% in gene promoters. This analysis provides a roadmap for generating reproducible base-pair level methylomes for zebrafish using RRBS and our results provide the first evidence that RRBS is a suitable technique for global methylation analysis in zebrafish.

Pseudomonas syringae pv. actinidiae from recent outbreaks of kiwifruit bacterial canker belong to different clones that originated in China

A recently emerged plant disease, bacterial canker of kiwifruit (Actinidia deliciosa and A. chinensis), is caused by Pseudomonas syringae pv. actinidiae (PSA). The disease was first reported in China and Japan in the 1980s. A severe outbreak of PSA began in Italy in 2008 and has spread to other European countries. PSA was found in both New Zealand and Chile in 2010. To study the evolution of the pathogen and analyse the transmission of PSA between countries, genomes of strains from China and Japan (where the genus Actinidia is endemic), Italy, New Zealand and Chile were sequenced. The genomes of PSA strains are very similar. However, all strains from New Zealand share several single nucleotide polymorphisms (SNPs) that distinguish them from all other PSA strains. Similarly, all the PSA strains from the 2008 Italian outbreak form a distinct clonal group and those from Chile form a third group. In addition to the rare SNPs present in the core genomes, there is abundant genetic diversity in a genomic island that is part of the accessory genome. The island from several Chinese strains is almost identical to the island present in the New Zealand strains. The island from a different Chinese strain is identical to the island present in the strains from the recent Italian outbreak. The Chilean strains of PSA carry a third variant of this island. These genomic islands are integrative conjugative elements (ICEs). Sequencing of these ICEs provides evidence of three recent horizontal transmissions of ICE from other strains of Pseudomonas syringae to PSA. The analyses of the core genome SNPs and the ICEs, combined with disease history, all support the hypothesis of an independent Chinese origin for both the Italian and the New Zealand outbreaks and suggest the Chilean strains also originate from China.

Comparison of alignment software for genome-wide bisulphite sequence data

Recent advances in next generation sequencing (NGS) technology now provide the opportunity to rapidly interrogate the methylation status of the genome. However, there are challenges in handling and interpretation of the methylation sequence data because of its large volume and the consequences of bisulphite modification. We sequenced reduced representation human genomes on the Illumina platform and efficiently mapped and visualized the data with different pipelines and software packages. We examined three pipelines for aligning bisulphite converted sequencing reads and compared their performance. We also comment on pre-processing and quality control of Illumina data. This comparison highlights differences in methods for NGS data processing and provides guidance to advance sequence-based methylation data analysis for molecular biologists.

Noncanonical translation initiation of the Arabidopsis flowering time and alternative polyadenylation regulator FCA

The RNA binding protein FCA regulates the floral transition and is required for silencing RNAs corresponding to specific noncoding sequences in the Arabidopsis thaliana genome. Through interaction with the canonical RNA 3' processing machinery, FCA affects alternative polyadenylation of many transcripts, including antisense RNAs at the locus encoding the floral repressor FLC. This potential for widespread alteration of gene regulation clearly needs to be tightly regulated, and we have previously shown that FCA expression is autoregulated through poly(A) site choice. Here, we show distinct layers of FCA regulation that involve sequences within the 5' region that regulate noncanonical translation initiation and alter the expression profile. FCA translation in vivo occurs exclusively at a noncanonical CUG codon upstream of the first in-frame AUG. We fully define the upstream flanking sequences essential for its selection, revealing features that distinguish this from other non-AUG start site mechanisms. Bioinformatic analysis identified 10 additional Arabidopsis genes that likely initiate translation at a CUG codon. Our findings reveal further unexpected complexity in the regulation of FCA expression with implications for its roles in regulating flowering time and gene expression and more generally show plant mRNA exceptions to AUG translation initiation.

Transterm: a database to aid the analysis of regulatory sequences in mRNAs

Messenger RNAs, in addition to coding for proteins, may contain regulatory elements that affect how the protein is translated. These include protein and microRNA-binding sites. Transterm (http://mRNA.otago.ac.nz/Transterm.html) is a database of regions and elements that affect translation with two major unique components. The first is integrated results of analysis of general features that affect translation (initiation, elongation, termination) for species or strains in Genbank, processed through a standard pipeline. The second is curated descriptions of experimentally determined regulatory elements that function as translational control elements in mRNAs. Transterm focuses on protein binding sites, particularly those in 3′-untranslated regions (3′-UTR). For this release the interface has been extensively updated based on user feedback. The data is now accessible by strain rather than species, for example there are 10 Escherichia coli strains (genomes) analysed separately. In addition to providing a repository of data, the database also provides tools for users to query their own mRNA sequences. Users can search sequences for Transterm or user defined regulatory elements, including protein or miRNA targets. Transterm also provides a central core of links to related resources for complementary analyses.

Transterm--extended search facilities and improved integration with other databases

Transterm has now been publicly available for >10 years. Major changes have been made since its last description in this database issue in 2002. The current database provides data for key regions of mRNA sequences, a curated database of mRNA motifs and tools to allow users to investigate their own motifs or mRNA sequences. The key mRNA regions database is derived computationally from Genbank. It contains 3' and 5' flanking regions, the initiation and termination signal context and coding sequence for annotated CDS features from Genbank and RefSeq. The database is non-redundant, enabling summary files and statistics to be prepared for each species. Advances include providing extended search facilities, the database may now be searched by BLAST in addition to regular expressions (patterns) allowing users to search for motifs such as known miRNA sequences, and the inclusion of RefSeq data. The database contains >40 motifs or structural patterns important for translational control. In this release, patterns from UTRsite and Rfam are also incorporated with cross-referencing. Users may search their sequence data with Transterm or user-defined patterns. The system is accessible at http://uther.otago.ac.nz/Transterm.html.

Transterm: a database of mRNAs and translational control elements

Transterm is a database that facilitates studies of translation and the translational control of protein synthesis. It contains a curated collection of elements in mRNAs that control translation, and biologically relevant mRNA regions extracted from GenBank. It is organised largely on a taxonomic basis with files and summaries for each species. Global patterns that may affect translation in particular species, for example bias in the context of initiation codons (Kozak's consensus or Shine-Dalgarno sequences) or termination codons, can be detected in the consensus and information content bias summaries. Several types of access are provided via a web browser interface. Transterm defined elements may be matched in a user's sequence or in the database. Alternatively, elements can be entered by the user to search specific sections of the database (for example, coding regions or 3' flanking regions or the 3'-UTRs) or the user's sequence. Each Transterm defined element has an associated biological description with references. The database is accessible at http://uther.otago.ac.nz/Transterm.html.

ToothPrint, a proteomic database for dental tissues

Increasing demand exists to disseminate and integrate proteomic data as proteome analysis assumes a commanding role in the postgenome era. Databases on the World Wide Web are an effective means to share information obtained from two-dimensional gels and allied proteomic approaches. Here we report the establishment of ToothPrint, a proteomic database for dental tissues accessed at http://toothprint.otago.ac.nz. Using developing rat enamel as a prototype, ToothPrint provides a variety of functionally relevant data (ligand binding, subcellular localisation, developmental regulation) in addition to protein identification maps. Features designed to enhance usability of the website and simplify its computing requirements are also outlined. Customized for mineralizing tissues, ToothPrint should prove to be an effective bioinformatic resource for investigations of dental biology.

Self-programmable, self-assembling two-dimensional genetic matter

Putative two-dimensional coding systems can be constructed from aqueous solutions of purine and pyrimidine nucleic acid bases evaporated at moderate temperatures on the surfaces of inorganic solids. The resultant structures are monolayers which are formed spontaneously by molecular self-assembly and they have been observed with molecular resolution by scanning tunnelling microscopy (STM). When formed from solutions of a single base, the monolayers of adenine and uracil have crystalline characteristics and the STM images can be interpreted in terms of the geometrical placement of planar arranged molecules that interact laterally by intermolecular hydrogen bonding. When formed from solutions containing a mixture of adenine and uracil, the monolayers have aperiodic structures. Small crystalline domains within these monolayers can be interpreted in terms of the single phase configurations of the molecules and the remaining aperiodic structures can presumably be interpreted, geometrically, in terms of the 21 theoretically possible adenine-adenine, uracil-uracil and adenine-uracil hydrogen bonding interactions. We propose that combinatorial arrangements of planar arranged purine and pyrimidine bases could provide the necessary complexity to act as a primitive genetic mechanism and may have relevance to the origin of life.

Transterm: a database of messenger RNA components and signals

Transterm facilitates studies of messenger RNAs and translational control signals. Each messenger RNA (mRNA) from GenBank is extracted and broken into its functional components, its coding sequence, initiation context, termination context, flanking sequence representing its 5' UTR (untranslated region), 3' UTR and translational signals. In addition, numerical parameters characterising each coding region in Transterm, including codon and GC bias, are available. For each species in Transterm, the initiation and termination regions are aligned by their start or stop codons and presented as base frequency matrices and tables of the information content of the bases in the alignments. Users can obtain summaries of characteristics of the mRNAs for species of their choice and search for translational signals both in the Transterm database and in their own sequence. The current release contains data from over 10 000 species, including the complete genomes of 20 prokaryotes and three eukaryotes. Both flat-file and relational database forms of Transterm are accessible via the WWW at http://biochem.otago.ac.nz/Transterm/

TransTerm, the translational signal database, extended to include full coding sequences and untranslated regions

TransTerm is a database of mRNA sequences and parameters useful for detecting translational control signals in general. TransTerm-98 has been expanded beyond previous years to include full coding sequences and UTRs, while retaining the original small contexts about the coding sequence start- and stop-codons. The database contains more than 130 000 non-redundant coding sequences with associated untranslated regions (UTRs) from over 450 species. This includes the complete genomes of 12 prokaryotic and one eukaryotic organism. Several coding sequence parameters are available: coding sequence length, Nc, GC3 and, when it is computable, Codon Adaptation Index (CAI). Codon usage tables and summaries of start- and stop-codon contexts are also included. TransTerm-98 has both a relational database form with a WWW interface and a flatfile format, also available by Internet browser. TransTerm is available at: http://biochem.otago.ac.nz:800/Transterm/homepage.h tml

Is gene deletion in eukaryotes sequence-dependent? A study of nine deletion junctions and nineteen other deletion breakpoints in intron 7 of the human dystrophin gene

Although large deletions comprise 65% of the mutations that underlie most cases of Duchenne and Becker muscular dystrophies, the DNA sequence characteristics of the deletions and the molecular processes leading to their formation are largely unknown. Intron 7 of the human dystrophin gene is unusually large (110 kb) and a substantial number of deletions have been identified with endpoints within this intron. The distribution of 28 deletion endpoints was mapped to local sequence elements by PCR. The break points were distributed among unique sequence, LINE-1, Alu, MIR, MER and microsatellite sequences with frequencies expected from the frequency of those sequences in the intron. Thus, deletions in this intron are not associated primarily with any one of those sequences in the intron. Nine deletion junctions were amplified and sequenced. Eight were deletions between DNA sequences with minimal homology (0-4 bp) and are therefore unlikely to be products of homologous recombination. In the ninth case, a complex rearrangement was found to be consistent with unequal recombinational exchange between two Alu sequences coupled with a duplication. We have hypothesized that a paucity of matrix attachment regions in this very large intron expanded by the insertion of many mobile elements might provoke a chromatin structure that stimulates deletions (McNaughton et al., 1997, Genomics 40, 294-304). The data presented here are consistent with that idea and demonstrate that the deletion sequences are not usually produced by homologous DNA misalignments.

The translational signal database, TransTerm, is now a relational database

TransTerm-97 contains more than 97 500 non-redundant coding-sequence initiation and termination contexts compiled from GenBank, release 101 (15-June-1997). In addition, several coding sequence parameters are available: coding sequence length, Nc, GC3, and, when it is computable, codon adaptation index (CAI). Codon usage tables and summaries of start and stop codon contexts are also included. The information covers more than 325 species and organelles, including seven complete bacterial genomes and one complete eukaryotic genome. To promote research in translational control of protein synthesis, TransTerm has been converted into a relational database to ease the process of making queries. The relational database manager, Postgresql, gives access to the database using SQL (Structured Query Language). A World Wide Web interface using forms is being completed to allow the casual user access to the database. Extensions are planned to include the full 5'-UTR, full coding sequence and 3'-UTR. TransTerm-97 is available on the World Wide Web at:http://biochem. otago.ac.nz:800/Transterm/homepage.html

The evolution of an intron: analysis of a long, deletion-prone intron in the human dystrophin gene

The sequence of a 112-kb region of the human dystrophin (DMD/BMD) gene encompassing the deletion prone intron 7 (110 kb) and the much shorter intron 8 (1.1 kb) has been determined. Recognizable insertion sequences account for approximately 40% of intron 7. LINE-1 and THE-1/LTR sequences occur in intron 7 with significantly higher frequency than would be expected statistically while Alu sequences are underrepresented. Intron 7 also contains numerous mammalian-wide interspersed repeats, a diverse range of medium reiteration repeats of unknown origin, and a sequence derived from a mariner transposon. By contrast, the shorter intron 8 contains no detectable insertion sequences. Dating of the LI and Alu sequences suggests that intron 7 has approximately doubled in size within the past 130 million years, and comparison with the corresponding intron from the pufferfish (Fugu rubripes) suggests that the intron has expanded some 44-fold over a period of 400 million years. The possible contribution of the insertion elements to the instability of intron 7 is discussed.

The translational signal database, TransTerm: more organisms, complete genomes

TransTerm is a database of initiation and termination sequence contexts from more than 250 organisms listed in GenBank, including the four complete genomes:Haemophilus influenzae, Methanococcus jannaschii, Mycoplasma genitalium,and Saccharomyces cerevisiae. For the current release, more than 60 000 coding sequences were analysed. The tabulated data include initiation and termination contexts organised by species along with quantitative parameters about individual coding sequences (length, %GC, GC3, Nc and CAI). There are also tables of initiation- and termination-region nucleotide-frequencies, codon usage tables and summaries of stop signal usage. TransTerm is available on the World Wide Web at: http://biochem.otago.ac.nz:800/Transterm/homepage.h tml

TransTerm: a database of translational signals

The TransTerm database of sequence contexts of stop and start codons has been expanded to include approximately 50% more species than last year's release. It now contains 148 organisms and >39 500 coding sequences; it is now available on the World Wide Web. The database includes: (i) initiation and termination sequence contexts organized by species; (ii) summary parameters about the individual sequences (sequence length, GC%, GC3, Nc, CAI) in addition to tables of base frequencies for each species' stop and start codon sequence context; (iii) species codon usage tables; and (iv) summary tables of stop signal frequency.

Phylogenetic relationships among transposon-like elements in human and primate DNA

A THE-1 sequence in intron 7 of the human dystrophin gene has been found to represent a new subfamily of THE-1 elements. The sequence is closely related to the MstII family of repetitive sequences and is more like single-copy sequences found in the galago genome than any other THE-1 sequence previously reported. This new THE-1 sequence has been compared with two other complete THE-1 sequences and three related long-terminal repeat elements that we have previously found in intron 7 of the dystrophin gene, and with members of the same family from elsewhere in the primate genome. Parsimony and deletion analysis show that the cluster of THE-1 sequences in intron 7 of the dystrophin gene has arisen from at least three individual insertion events, rather than from the insertion and duplication of a single progenitor sequence.

Sequence of a transposon identified as Tn1000 (gamma delta)

We report the complete sequence of a transposon found in a cosmid clone of a human DNA sequence. The transposon is identified as the Escherichia coli transposon Tn1000 (also known as gamma delta) on the basis of the identity of the restriction map of the new sequence with that previously recorded for Tn1000 and homology between parts of the new sequence and that of published fragments of Tn1000 sequence. The transposon, which comprises 5,981 nucleotides including two 35 bp inverted terminal repeat sequences (ITRs), contains three open reading frames. The sequence of the resolvase coding region (tnpR) is identical to that published by others. A second reading frame can be identified as the tnpA gene, coding for the transposase, on the grounds of its strong homology with the corresponding gene from transposon Tn3. The third reading frame has the potential to code for a protein of unknown function containing 698 amino acids.

The translational termination signal database (TransTerm) now also includes initiation contexts

The TransTerm database of termination codon contexts has been extended to include sense codon usage, and initiation codon contexts. The database was constructed from 23,721 coding sequences from 93 organisms. The database contains: a) the sequence around the termination codon (-10, +10); b) the sequence around the initiation codon (-20, +10); c) the length, 'G+C%' of the third position of codons (GC3), the 'codon adaptation index' (CAI) and the 'effective number of codons' statistic (Nc); d) summary tables for each organism including total codon usage, stop codon and tetranucleotide stop-signal usage, and matrices tallying base frequencies at each position around the initiation and termination codons. The data are arranged to facilitate investigation of the relationships between the three phases of protein synthesis. The database is available electronically from EMBL.

A cluster of transposon-like repetitive sequences in intron 7 of the human dystrophin gene

A 32 kilobase-pair fragment of intron 7 of the human dystrophin gene has been sequenced and analysed for the presence of repetitive elements and open reading frames. Two transposon-like human elements (THE-1 sequences), and three intervening, and related, long terminal repeat elements, together with three Alu sequences and a LINE sequence have been identified. These represent an unexpected clustering of highly-repetitive sequences within this single segment of intron DNA. Amplification of a region of chimpanzee genomic DNA by the polymerase chain reaction has provided evidence that at least one of the THE-1 sequences is present in the same position in the chimpanzee genome and the high homology between the human and chimpanzee sequences indicates that this element was fixed within the ancestral genome before the divergence of the two species. The possible role of repetitive, transposon-like sequences in natural mutagenesis of the dystrophin gene is discussed.

The translational termination signal database

The Translational Termination Database (TransTerm) consists of the immediate context sequences around the natural termination codons from 45 organisms, and summary tables. The influence of termination codon context on their effectivness as stop signals has been widely documented. The SPECIES--TRI.DAT table shows trinucleotide stop codon usage in each organism and for comparison the occurrence of these sequences in the noncoding region. The SPECIES--TETRA.DAT table contains is a similar table of tetranucleotide stop signal usage. The database is available from EMBL.

The putative single-stranded DNA-binding protein of the filamentous bacteriophage, Ifl. Amino acid sequence of the protein and structure of the gene

The protein product corresponding to the gene located in the region of the coliphage Ifl genome shown to contain the code for the single-stranded DNA (ssDNA)-binding proteins of all filamentous phages so far studied has been isolated from infected bacterial cells and its amino acid sequence determined. The mature protein contains 95 amino acids (calculated molecular mass 10553 Da). Its sequence corresponds to that predicted from the DNA sequence but lacks the initiating methionine residue. Although there is little direct sequence homology between the phage Ifl protein and the ssDNA-binding proteins of the other filamentous phages that have been studied, computer-based comparisons of various physical and structural parameters showed that the phage Ifl protein contains a domain that is closely related to domains in the coliphage T4 gene 32 protein and the Pseudomonas phage Pfl ssDNA-binding protein and suggest that the Ifl protein does have a ssDNA-binding function although we were unable to show this directly.

Sequence analysis suggests that tetra-nucleotides signal the termination of protein synthesis in eukaryotes

An increasing number of cases where tri-nucleotide stop codons do not signal the termination of protein synthesis are being reported. In order to identify what constitutes an efficient stop signal, we analysed the region around natural stop codons in genes from a wide variety of eukaryotic species and gene families. Certain stop codons and nucleotides following stop codons are over-represented, and this pattern is accentuated in highly expressed genes. For example, the preferred signal for Saccharomyces cerevisiae and Drosophila melanogaster highly expressed genes is UAAG, and generally the signals UAA(A/G) and UGA(A/G) are preferred in eukaryotes. The GC% of the organism or DNA region can affect whether there is A or G in the second or fourth positions. We suggest therefore, that the stop codon and the nucleotide following it comprise a tetra-nucleotide stop signal. A model is proposed in which the polypeptide chain release factor, a protein, recognises this sequence, but will tolerate some substitution, particularly A to G in the second or third positions.