Mammalian Genome Sequences Research Articles

OrthoMaM v12: a database of curated single-copy ortholog alignments and trees to study mammalian evolutionary genomics.

To date, the databases built to gather information on gene orthology do not provide end-users with descriptors of the molecular evolution information and phylogenetic pattern of these orthologues. In this context, we developed OrthoMaM, a database of ORTHOlogous MAmmalian Markers describing the evolutionary dynamics of coding sequences in mammalian genomes. OrthoMaM version 12 includes 15,868 alignments of orthologous coding sequences (CDS) from the 190 complete mammalian genomes currently available. All annotations and 1-to-1 orthology assignments are based on NCBI. Orthologous CDS can be mined for potential informative markers at the different taxonomic levels of the mammalian tree. To this end, several evolutionary descriptors of DNA sequences are provided for querying purposes (e.g. base composition and relative substitution rate). The graphical web interface allows the user to easily browse and sort the results of combined queries. The corresponding multiple sequence alignments and ML trees, inferred using state-of-the art approaches, are available for download both at the nucleotide and amino acid levels. OrthoMaM v12 can be used by researchers interested either in reconstructing the phylogenetic relationships of mammalian taxa or in understanding the evolutionary dynamics of coding sequences in their genomes. OrthoMaM is available for browsing, querying and complete or filtered download at https://orthomam.mbb.cnrs.fr/.

Transposable elements in mammalian chromatin organization.

Transposable elements (TEs) are mobile DNA elements that comprise almost 50% of mammalian genomic sequence. TEs are capable of making additional copies of themselves that integrate into new positions in host genomes. This unique property has had an important impact on mammalian genome evolution and on the regulation of gene expression because TE-derived sequences can function as cis-regulatory elements such as enhancers, promoters and silencers. Now, advances in our ability to identify and characterize TEs have revealed that TE-derived sequences also regulate gene expression by both maintaining and shaping 3D genome architecture. Studies are revealing how TEs contribute raw sequence that can give rise to the structures that shape chromatin organization, and thus gene expression, allowing for species-specific genome innovation and evolutionary novelty.

EvalDNA: a machine learning-based tool for the comprehensive evaluation of mammalian genome assembly quality

BackgroundTo select the most complete, continuous, and accurate assembly for an organism of interest, comprehensive quality assessment of assemblies is necessary. We present a novel tool, called Evaluation of De Novo Assemblies (EvalDNA), which uses supervised machine learning for the quality scoring of genome assemblies and does not require an existing reference genome for accuracy assessment.ResultsEvalDNA calculates a list of quality metrics from an assembled sequence and applies a model created from supervised machine learning methods to integrate various metrics into a comprehensive quality score. A well-tested, accurate model for scoring mammalian genome sequences is provided as part of EvalDNA. This random forest regression model evaluates an assembled sequence based on continuity, completeness, and accuracy, and was able to explain 86% of the variation in reference-based quality scores within the testing data. EvalDNA was applied to human chromosome 14 assemblies from the GAGE study to rank genome assemblers and to compare EvalDNA to two other quality evaluation tools. In addition, EvalDNA was used to evaluate several genome assemblies of the Chinese hamster genome to help establish a better reference genome for the biopharmaceutical manufacturing community. EvalDNA was also used to assess more recent human assemblies from the QUAST-LG study completed in 2018, and its ability to score bacterial genomes was examined through application on bacterial assemblies from the GAGE-B study.ConclusionsEvalDNA enables scientists to easily identify the best available genome assembly for their organism of interest without requiring a reference assembly. EvalDNA sets itself apart from other quality assessment tools by producing a quality score that enables direct comparison among assemblies from different species.

Evolution of the mammalian insulin (Ins) gene; Changes in proteolytic processing

Disruption of insulin signaling in humans leads to diabetes yet changes in insulin function is tolerated in some species. Taking advantage of the large number of publicly available mammalian genome sequences I identified insulin gene (Ins) in the genomes of 151 of 156 mammalian species with well-annotated genomes, of which 141 had complete Ins coding sequences. Complete Ins coding sequences were identified from 8 additional species that lack complete genomes. Duplicated Ins genes were found in 12 rodents (9 with complete genomes) resulting in the identification of a total of 161 complete mammalian Ins coding sequences. While all 161 proinsulin protein sequences were predicted to have functional signal peptides, which should allow secretion of the hormone, unexpectedly, substitutions were found at prohormone convertase processing sites in sequences from 6 species, 2 from Chiroptera (Myotis brandtii and M. lucifugus) and 4 from Afrotheria (Chrysochloris asiatica, Echinops telfairi, Elephantulus edwardii, and Orycteropus afer). Both basic residues at the C-peptide-A-chain junction in the bats M. brandtii and M. lucifugas are replaced, which should prevent processing. Replacements of a single basic residue are found at the B-chain-C-peptide junction, in the two bats, and at the C-peptide-A-chain junction, in 4 species of Afrotheria, processing sites that suggest impaired processing. In addition, a large number of substitutions at sites that interact with the insulin receptor were found in the insulin sequences from M. brandtii and M. lucifugas suggesting a change in biological function.

Roles of Noncoding RNAs in Islet Biology.

The discovery that most mammalian genome sequences are transcribed to ribonucleic acids (RNA) has revolutionized our understanding of the mechanisms governing key cellular processes and of the causes of human diseases, including diabetes mellitus. Pancreatic islet cells were found to contain thousands of noncoding RNAs (ncRNAs), including micro-RNAs (miRNAs), PIWI-associated RNAs, small nucleolar RNAs, tRNA-derived fragments, long non-coding RNAs, and circular RNAs. While the involvement of miRNAs in islet function and in the etiology of diabetes is now well documented, there is emerging evidence indicating that other classes of ncRNAs are also participating in different aspects of islet physiology. The aim of this article will be to provide a comprehensive and updated view of the studies carried out in human samples and rodent models over the past 15 years on the role of ncRNAs in the control of α- and β-cell development and function and to highlight the recent discoveries in the field. We not only describe the role of ncRNAs in the control of insulin and glucagon secretion but also address the contribution of these regulatory molecules in the proliferation and survival of islet cells under physiological and pathological conditions. It is now well established that most cells release part of their ncRNAs inside small extracellular vesicles, allowing the delivery of genetic material to neighboring or distantly located target cells. The role of these secreted RNAs in cell-to-cell communication between β-cells and other metabolic tissues as well as their potential use as diabetes biomarkers will be discussed. © 2020 American Physiological Society. Compr Physiol 10:893-932, 2020.

Transposable elements as a potent source of diverse cis-regulatory sequences in mammalian genomes.

Eukaryotic gene regulation is mediated by cis-regulatory elements, which are embedded within the vast non-coding genomic space and recognized by the transcription factors in a sequence- and context-dependent manner. A large proportion of eukaryotic genomes, including at least half of the human genome, are composed of transposable elements (TEs), which in their ancestral form carried their own cis-regulatory sequences able to exploit the host trans environment to promote TE transcription and facilitate transposition. Although not all present-day TE copies have retained this regulatory function, the preexisting regulatory potential of TEs can provide a rich source of cis-regulatory innovation for the host. Here, we review recent evidence documenting diverse contributions of TE sequences to gene regulation by functioning as enhancers, promoters, silencers and boundary elements. We discuss how TE-derived enhancer sequences can rapidly facilitate changes in existing gene regulatory networks and mediate species- and cell-type-specific regulatory innovations, and we postulate a unique contribution of TEs to species-specific gene expression divergence in pluripotency and early embryogenesis. With advances in genome-wide technologies and analyses, systematic investigation of TEs' cis-regulatory potential is now possible and our understanding of the biological impact of genomic TEs is increasing.This article is part of a discussion meeting issue ‘Crossroads between transposons and gene regulation’.

Evolution of the RNA N 6-Methyladenosine Methylome Mediated by Genomic Duplication.

RNA N 6-methyladenosine (m6A) modification is the most abundant form of RNA epigenetic modification in eukaryotes. Given that m6A evolution is associated with the selective constraints of nucleotide sequences in mammalian genomes, we hypothesize that m6A evolution can be linked, at least in part, to genomic duplication events in complex polyploid plant genomes. To test this hypothesis, we presented the maize (Zea mays) m6A modification landscape in a transcriptome-wide manner and identified 11,968 m6A peaks carried by 5,893 and 3,811 genes from two subgenomes (maize1 and maize2, respectively). Each of these subgenomes covered over 2,200 duplicate genes. Within these duplicate genes, those carrying m6A peaks exhibited significant differences in retention rate. This biased subgenome fractionation of m6A-methylated genes is associated with multiple sequence features and is influenced by asymmetric evolutionary rates. We also characterized the coevolutionary patterns of m6A-methylated genes and transposable elements, which can be mediated by whole genome duplication and tandem duplication. We revealed the evolutionary conservation and divergence of duplicated m6A functional factors and the potential role of m6A modification in maize responses to drought stress. This study highlights complex interplays between m6A modification and gene duplication, providing a reference for understanding the mechanisms underlying m6A evolution mediated by genome duplication events.

Identification and Characterization of the OCT4 Upstream Regulatory Region in Sus scrofa.

OCT4 plays pivotal roles in maintaining pluripotency during early mammalian embryonic development and in embryonic stem cells. It is essential to establish a reporter system based on the OCT4 promoter region to study pluripotency. However, there is still a lack of information about the porcine OCT4 upstream reporter system. To improve our understanding of the porcine OCT4 regulatory region, we identified conserved regions in the porcine OCT4 promoter upstream region by sequence-based comparative analysis using various mammalian genome sequences. The similarity of nucleotide sequences in the 5′ upstream region was low among mammalian species. However, the OCT4 promoter and four regulatory regions, including distal and proximal enhancer elements, had high similarity. Next, a functional analysis of the porcine OCT4 promoter region was conducted. Luciferase reporter assay results indicated that the porcine OCT4 distal enhancer and proximal enhancer were highly activated in mouse embryonic stem cells and embryonic carcinoma cells, respectively. A comparison analysis of naïve and primed state marker gene expression in a dual-reporter assay showed that the expression levels of naïve and primed markers differed in fluorescence signal between high-expressing cells and low-expressing cells. Similar to OCT4 upstream-based reporter systems derived from other species, the porcine OCT4 upstream region-based reporter constructs showed exclusive expression patterns depending on the state of pluripotency. This work provides basic information about the porcine OCT4 upstream region and various porcine OCT4 fluorescence reporter constructs, which can be applied to study species-specific pluripotency in early embryo development and the establishment of embryonic stem cells in pigs.

124 Functional analysis of porcine OCT4 transcriptional regulatory region-based reporter system

Gene OCT4 plays pivotal roles in maintaining pluripotency of early mammalian embryonic development and embryonic stem cells. It is essential to establish a reporter system based on the OCT4 promoter region for the study of pluripotency. However, there is still a lack of sufficient information about the porcine OCT4 upstream reporter system. To improve our understanding of the porcine OCT4 regulatory region, first, we conducted an investigation to find conserved regions in the porcine OCT4 promoter upstream region by sequence-based comparative analysis using various mammalian genome sequences. A similarity of nucleotide sequences of the 5′ upstream region was low among mammalian species. However, the OCT4 promoter and 4 regulatory regions including distal and proximal enhancer elements have a high similarity. Next, a functional analysis of the porcine OCT4 promoter region was conducted. Luciferase reporter assay indicated that the porcine OCT4 distal enhancer and proximal enhancer are highly activated in mouse embryonic stem cells and embryonic carcinoma cells, respectively (n=3). Comparison analysis of naïve (Tbx3, Nr0b1, Rex1, Esrrb, Nanog, Klf2) or primed (Gata6, Mixl1, Fgf5, Otx2) state marker gene expression in a dual-reporter assay using pOCT4-DE-eGFP and pOCT4-PE-DsRed2 showed that expression of naïve and primed markers were up-regulated in cells with high green fluorescent protein and red fluorescent protein expression, respectively (n=3). Porcine OCT4-upstream region-based reporter constructs showed exclusive expression patterns depending on the state of pluripotency. This work could provide basic information for the porcine OCT4 upstream region and the various porcine OCT4-fluorescence reporter constructs, which can be applied to study species-specific pluripotency in early embryo development and for the establishment of embryonic stem cells in pigs. This work was supported by the Korea Institute of Planning and Evaluation for Technology in Food, Agriculture, Forestry and Fisheries (IPET) through the Development of High Value-Added Food Technology Program, funded by the Ministry of Agriculture, Food and Rural Affairs (MAFRA, 118042-03-1-HD020).

Genome classification improvements based on k-mer intervals in sequences

Given the vast amount of genomic data, alignment-free sequence comparison methods are required due to their low computational complexity. k-mer based methods can improve comparison accuracy by extracting an effective feature of the genome sequences. The aim of this paper is to extract k-mer intervals of a sequence as a feature of a genome for high comparison accuracy. In the proposed method, we calculated the distance between genome sequences by comparing the distribution of k-mer intervals. Then, we identified the classification results using phylogenetic trees. We used viral, mitochondrial (MT), microbial and mammalian genome sequences to perform classification for various genome sets. We confirmed that the proposed method provides a better classification result than other k-mer based methods. Furthermore, the proposed method could efficiently be applied to long sequences such as human and mouse genomes.

Peroxisome Proliferator Activated Receptor Agonists Modulate Transposable Element Expression in Brain and Liver.

Peroxisome proliferator activated receptors (PPARs) are nuclear hormone receptors that act as transcription factors in response to endogenous lipid messengers. The fibrates and thiazolidinediones are synthetic PPAR agonists used clinically to treat dyslipidemia and Type 2 Diabetes Mellitus, respectively, but also improve symptoms of several other diseases. Transposable elements (TEs), repetitive sequences in mammalian genomes, are implicated in many of the same conditions for which PPAR agonists are therapeutic, including neurodegeneration, schizophrenia, and drug addiction. We tested the hypothesis that there is a link between actions of PPAR agonists and TE expression. We developed an innovative application of microarray data by mapping Illumina mouse WG-6 microarray probes to areas of the mouse genome that contain TEs. Using this information, we assessed the effects of systemic administration of three PPAR agonists with different PPAR subtype selectivity: fenofibrate, tesaglitazar, and bezafibrate, on TE probe expression in mouse brain [prefrontal cortex (PFC) and amygdala] and liver. We found that fenofibrate, and bezafibrate to a lesser extent, up-regulated probes mapped to retrotransposons: Short-Interspersed Elements (SINEs) and Long-Interspersed Elements (LINEs), in the PFC. Conversely, all PPAR agonists down-regulated LINEs and tesaglitazar and bezafibrate also down-regulated SINEs in liver. We built gene coexpression networks that partitioned the diverse transcriptional response to PPAR agonists into groups of probes with highly correlated expression patterns (modules). Most of the differentially expressed retrotransposons were within the same module, suggesting coordinated regulation of their expression, possibly by PPAR signaling. One TE module was conserved across tissues and was enriched with genes whose products participate in epigenetic regulation, suggesting that PPAR agonists affect TE expression via epigenetic mechanisms. Other enriched functional categories included phenotypes related to embryonic development and learning and memory, suggesting functional links between these biological processes and TE expression. In summary, these findings suggest mechanistic relationships between retrotransposons and PPAR agonists and provide a basis for future exploration of their functional roles in brain and liver.

Lineage-associated underrepresented permutations (LAUPs) of mammalian genomic sequences based on a Jellyfish-based LAUPs analysis application (JBLA).

This study addresses several important questions related to naturally underrepresented sequences: (i) are there permutations of real genomic DNA sequences in a defined length (k-mer) and a given lineage that do not actually exist or underrepresented? (ii) If there are such sequences, what are their characteristics in terms of k-mer length and base composition? (iii) Are they related to CpG or TpA underrepresentation known for human sequences? We propose that the answers to these questions are of great significance for the study of sequence-associated regulatory mechanisms, such cytosine methylation and chromosomal structures in physiological or pathological conditions such as cancer. We empirically defined sequences that were not included in any well-known public databases as lineage-associated underrepresented permutations (LAUPs). Then, we developed a Jellyfish-based LAUPs analysis application (JBLA) to investigate LAUPs for 24 representative species. The present discoveries include: (i) lengths for the shortest LAUPs, ranging from 10 to 14, which collectively constitute a low proportion of the genome. (ii) Common LAUPs showing higher CG content over the analysed mammalian genome and possessing distinct CG*CG motifs. (iii) Neither CpG-containing LAUPs nor CpG island sequences are randomly structured and distributed over the genomes; some LAUPs and most CpG-containing sequences exhibit an opposite trend within the same k and n variants. In addition, we demonstrate that the JBLA algorithm is more efficient than the original Jellyfish for computing LAUPs. We developed a Jellyfish-based LAUP analysis (JBLA) application by integrating Jellyfish (Marçais and Kingsford, 2011), MEME (Bailey, et al., 2009) and the NCBI genome database (Pruitt, et al., 2007) applications, which are listed as Supplementary Material. Supplementary data are available at Bioinformatics online.

SinEx DB: a database for single exon coding sequences in mammalian genomes.

Eukaryotic genes are typically interrupted by intragenic, noncoding sequences termed introns. However, some genes lack introns in their coding sequence (CDS) and are generally known as ‘single exon genes’ (SEGs). In this work, a SEG is defined as a nuclear, protein-coding gene that lacks introns in its CDS. Whereas, many public databases of Eukaryotic multi-exon genes are available, there are only two specialized databases for SEGs. The present work addresses the need for a more extensive and diverse database by creating SinEx DB, a publicly available, searchable database of predicted SEGs from 10 completely sequenced mammalian genomes including human. SinEx DB houses the DNA and protein sequence information of these SEGs and includes their functional predictions (KOG) and the relative distribution of these functions within species. The information is stored in a relational database built with My SQL Server 5.1.33 and the complete dataset of SEG sequences and their functional predictions are available for downloading. SinEx DB can be interrogated by: (i) a browsable phylogenetic schema, (ii) carrying out BLAST searches to the in-house SinEx DB of SEGs and (iii) via an advanced search mode in which the database can be searched by key words and any combination of searches by species and predicted functions. SinEx DB provides a rich source of information for advancing our understanding of the evolution and function of SEGs.Database URL: www.sinex.cl

Promoter-distal RNA polymerase II binding discriminates active from inactive CCAAT/ enhancer-binding protein beta binding sites.

Transcription factors (TFs) bind to thousands of DNA sequences in mammalian genomes, but most of these binding events appear to have no direct effect on gene expression. It is unclear why only a subset of TF bound sites are actively involved in transcriptional regulation. Moreover, the key genomic features that accurately discriminate between active and inactive TF binding events remain ambiguous. Recent studies have identified promoter-distal RNA polymerase II (RNAP2) binding at enhancer elements, suggesting that these interactions may serve as a marker for active regulatory sequences. Despite these correlative analyses, a thorough functional validation of these genomic co-occupancies is still lacking. To characterize the gene regulatory activity of DNA sequences underlying promoter-distal TF binding events that co-occur with RNAP2 and TF sites devoid of RNAP2 occupancy using a functional reporter assay, we performed cis-regulatory element sequencing (CRE-seq). We tested more than 1000 promoter-distal CCAAT/enhancer-binding protein beta (CEBPB)-bound sites in HepG2 and K562 cells, and found that CEBPB-bound sites co-occurring with RNAP2 were more likely to exhibit enhancer activity. CEBPB-bound sites further maintained substantial cell-type specificity, indicating that local DNA sequence can accurately convey cell-type–specific regulatory information. By comparing our CRE-seq results to a comprehensive set of genome annotations, we identified a variety of genomic features that are strong predictors of regulatory element activity and cell-type–specific activity. Collectively, our functional assay results indicate that RNAP2 occupancy can be used as a key genomic marker that can distinguish active from inactive TF bound sites.

GREAM: A Web Server to Short-List Potentially Important Genomic Repeat Elements Based on Over-/Under-Representation in Specific Chromosomal Locations, Such as the Gene Neighborhoods, within or across 17 Mammalian Species.

BackgroundGenome-wide repeat sequences, such as LINEs, SINEs and LTRs share a considerable part of the mammalian nuclear genomes. These repeat elements seem to be important for multiple functions including the regulation of transcription initiation, alternative splicing and DNA methylation. But it is not possible to study all repeats and, hence, it would help to short-list before exploring their potential functional significance via experimental studies and/or detailed in silico analyses.ResultWe developed the ‘Genomic Repeat Element Analyzer for Mammals’ (GREAM) for analysis, screening and selection of potentially important mammalian genomic repeats. This web-server offers many novel utilities. For example, this is the only tool that can reveal a categorized list of specific types of transposons, retro-transposons and other genome-wide repetitive elements that are statistically over-/under-represented in regions around a set of genes, such as those expressed differentially in a disease condition. The output displays the position and frequency of identified elements within the specified regions. In addition, GREAM offers two other types of analyses of genomic repeat sequences: a) enrichment within chromosomal region(s) of interest, and b) comparative distribution across the neighborhood of orthologous genes. GREAM successfully short-listed a repeat element (MER20) known to contain functional motifs. In other case studies, we could use GREAM to short-list repetitive elements in the azoospermia factor a (AZFa) region of the human Y chromosome and those around the genes associated with rat liver injury. GREAM could also identify five over-represented repeats around some of the human and mouse transcription factor coding genes that had conserved expression patterns across the two species.ConclusionGREAM has been developed to provide an impetus to research on the role of repetitive sequences in mammalian genomes by offering easy selection of more interesting repeats in various contexts/regions. GREAM is freely available at http://resource.ibab.ac.in/GREAM/.

Identification of a Recently Active Mammalian SINE Derived from Ribosomal RNA

Complex eukaryotic genomes are riddled with repeated sequences whose derivation does not coincide with phylogenetic history and thus is often unknown. Among such sequences, the capacity for transcriptional activity coupled with the adaptive use of reverse transcription can lead to a diverse group of genomic elements across taxa, otherwise known as selfish elements or mobile elements. Short interspersed nuclear elements (SINEs) are nonautonomous mobile elements found in eukaryotic genomes, typically derived from cellular RNAs such as tRNAs, 7SL or 5S rRNA. Here, we identify and characterize a previously unknown SINE derived from the 3′-end of the large ribosomal subunit (LSU or 28S rDNA) and transcribed via RNA polymerase III. This new element, SINE28, is represented in low-copy numbers in the human reference genome assembly, wherein we have identified 27 discrete loci. Phylogenetic analysis indicates these elements have been transpositionally active within primate lineages as recently as 6 MYA while modern humans still carry transcriptionally active copies. Moreover, we have identified SINE28s in all currently available assembled mammalian genome sequences. Phylogenetic comparisons indicate that these elements are frequently rederived from the highly conserved LSU rRNA sequences in a lineage-specific manner. We propose that this element has not been previously recognized as a SINE given its high identity to the canonical LSU, and that SINE28 likely represents one of possibly many unidentified, active transposable elements within mammalian genomes.

Functions of noncoding sequences in mammalian genomes.

Most of the mammalian genome consists of nucleotide sequences not coding for proteins. Exons of genes make up only 3% of the human genome, while the significance of most other sequences remains unknown. Recent genome studies with high-throughput methods demonstrate that the so-called noncoding part of the genome may perform important functions. This hypothesis is supported by three groups of experimental data: 1) approximately 10% of the sequences, most of which are located in noncoding parts of the genome, is evolutionarily conserved and thus can be of functional importance; 2) up to 99% of the mammalian genome is being transcribed forming short and long noncoding RNAs in addition to common mRNA; and 3) mutations in noncoding parts of the genome can be accompanied by progression of pathological states of the organism. In the light of these data, in the review we consider the functional role of numerous known sequences of noncoding parts of the genome including introns, DNA methylation regions, enhancers and locus control regions, insulators, S/MAR sequences, pseudogenes, and genes of noncoding RNAs, as well as transposons and simple repeats of centromeric and telomeric regions of chromosomes. The assumption is made that the intergenic noncoding sequences without definite/clear functions can be involved in spatial organization of genetic loci in interphase nuclei.

Function-based identification of mammalian enhancers using site-specific integration.

The accurate and comprehensive identification of functional regulatory sequences in mammalian genomes remains a major challenge. Here we describe Site-specific Integration FACS-sequencing (SIF-seq), an unbiased, medium-throughput functional assay for the discovery of distant-acting enhancers. Pluripotent cell reporter assays, targeted single-copy genomic integration, and flow cytometry are coupled with high-throughput DNA sequencing to enable parallel screening of large numbers of DNA sequences. We demonstrate the utility of this method by functionally interrogating >500 kb of mouse and human sequence for enhancer activity and identifying embryonic stem (ES) cell enhancers at pluripotency loci including NANOG. We also demonstrate the effectiveness of the approach in differentiated cell populations through the identification of cardiac enhancers from cardiomyocytes and neuronal enhancers from neural progenitors. SIF-seq is a powerful and flexible method for the de novo functional identification of mammalian enhancers in a potentially wide variety of cell types.

Segmented K-mer and its application on similarity analysis of mitochondrial genome sequences

K-mer-based approach has been widely used in similarity analyses so as to discover similarity/dissimilarity among different biological sequences. In this study, we have improved the traditional K-mer method, and introduce a segmented K-mer approach (s-K-mer). After each primary sequence is divided into several segments, we simultaneously transform all these segments into corresponding K-mer-based vectors. In this approach, it is vital how to determine the optimal combination of distance metric with the number of K and the number of segments, i.e., (K⁎, s⁎, and d⁎). Based on the cascaded feature vectors transformed from s⁎ segmented sequences, we analyze 34 mammalian genome sequences using the proposed s-K-mer approach. Meanwhile, we compare the results of s-K-mer with those of traditional K-mer. The contrastive analysis results demonstrate that s-K-mer approach outperforms the traditionally K-mer method on similarity analysis among different species.

CpGislandEVO: A Database and Genome Browser for Comparative Evolutionary Genomics of CpG Islands

Hypomethylated, CpG-rich DNA segments (CpG islands, CGIs) are epigenome markers involved in key biological processes. Aberrant methylation is implicated in the appearance of several disorders as cancer, immunodeficiency, or centromere instability. Furthermore, methylation differences at promoter regions between human and chimpanzee strongly associate with genes involved in neurological/psychological disorders and cancers. Therefore, the evolutionary comparative analyses of CGIs can provide insights on the functional role of these epigenome markers in both health and disease. Given the lack of specific tools, we developed CpGislandEVO. Briefly, we first compile a database of statistically significant CGIs for the best assembled mammalian genome sequences available to date. Second, by means of a coupled browser front-end, we focus on the CGIs overlapping orthologous genes extracted from OrthoDB, thus ensuring the comparison between CGIs located on truly homologous genome segments. This allows comparing the main compositional features between homologous CGIs. Finally, to facilitate nucleotide comparisons, we lifted genome coordinates between assemblies from different species, which enables the analysis of sequence divergence by direct count of nucleotide substitutions and indels occurring between homologous CGIs. The resulting CpGislandEVO database, linking together CGIs and single-cytosine DNA methylation data from several mammalian species, is freely available at our website.