Cross-species Alignment Research Articles

Cross-species alignment along the chronological axis reveals evolutionary effect on structural development of the human brain.

Disentangling the evolution mysteries of the human brain has always been an imperative endeavor in neuroscience. Although many previous comparative studies revealed genetic, brain structural and connectivity distinctness between human and other nonhuman primates, the brain evolutional mechanism is still largely unclear. Here, we proposed to embed the brain anatomy of human and macaque in the developmental chronological axis to construct cross-species predictive model to quantitatively characterize brain evolution using two large public human and macaque datasets. We observed that applying the trained models within-species could well predict the chronological age. Interestingly, we found the model trained in macaque showed a higher accuracy in predicting the chronological age of human than the model trained in human in predicting the chronological age of macaque. The cross-application of the trained model introduced an individual brain cross-species age gap index to quantify the cross-species discrepancy along the temporal axis of brain development and was found to be associated with the behavioral performance in visual acuity test and picture vocabulary test in human. Taken together, our study situated the cross-species brain development along the chronological axis, which highlighted the disproportionately anatomical development in human brain to extend our understanding of the potential evolutionary effects.

Approximate IsoRank for Scalable and Functionally Meaningful Cross-Species Alignments of Protein Interaction Networks.

The IsoRank algorithm of Singh, Xu, and Berger was a pioneering algorithmic advance that applied spectral methods to the problem of cross-species global alignment of biological networks. We develop a new IsoRank approximation that exploits the mathematical properties of IsoRank's linear system to solve the problem in quadratic time with respect to the maximum size of the two protein-protein interaction (PPI) networks. We further propose a refinement to this initial approximation so that the updated result is even closer to the original IsoRank formulation while remaining computationally inexpensive. In experiments on synthetic and real PPI networks with various proposed metrics to measure alignment quality, we find the results of our approximate IsoRank are nearly as accurate as the original IsoRank. In fact, for functional enrichment-based measures of global network alignment quality, our approximation performs better than the exact IsoRank, which is doubtless because it is more robust to the noise of missing or incorrect edges. It also performs competitively against two more recent global network alignment algorithms. We also present an analogous approximation to IsoRankN, which extends the network alignment to more than two species.

The NCBI Comparative Genome Viewer (CGV) is an interactive visualization tool for the analysis of whole-genome eukaryotic alignments.

We report a new visualization tool for analysis of whole-genome assembly-assembly alignments, the Comparative Genome Viewer (CGV) (https://ncbi.nlm.nih.gov/genome/cgv/). CGV visualizes pairwise same-species and cross-species alignments provided by National Center for Biotechnology Information (NCBI) using assembly alignment algorithms developed by us and others. Researchers can examine large structural differences spanning chromosomes, such as inversions or translocations. Users can also navigate to regions of interest, where they can detect and analyze smaller-scale deletions and rearrangements within specific chromosome or gene regions. RefSeq or user-provided gene annotation is displayed where available. CGV currently provides approximately 800 alignments from over 350 animal, plant, and fungal species. CGV and related NCBI viewers are undergoing active development to further meet needs of the research community in comparative genome visualization.

Saccharomyces cerevisiae cis-acting DNA sequences curation pipeline (Sc-cADSs-CP): Master transcription factors prediction in yeasts

Designing promoter architectures hinges on genomic and functional annotation. Saccharomyces cerevisiae is the first model yeast whose databases host genomic and functional annotation information. To predict transcription factors (TFs) regulating central pathways in yeasts, we first introduce S. cerevisiae cis-acting DNA sequences/sites (cADSs) curation pipeline (Sc-cADSs-CP). The promoters of the genes involved in the central pathways of S. cerevisiae were retrieved from the genome sequences. We processed the binding frequency matrices of TFs with the following two criteria. First, we extracted cADSs based on the TF motifs in the TRANSFAC database; then, if there were more than one frequency matrix for a TF, the longest one with the maximum sensitivity was used. Next, we developed the direct scanning algorithm ScanAlgo (uses the Biopython-library), scanning DNA motifs for pairing cross-species alignments. We used the tools Sc-cADSs-CP and ScanAlgo to predict master TFs for Pichia pastoris, which lacks extensive functional annotation studies. The phylogenetic footprinting results were obtained by aligning the scanned S. cerevisiae promoters against orthologous P. pastoris promoters. The predicted cADSs were summed into position weight matrices unique to P. pastoris. We annotated 116 TFs based on the phylogenetic footprinting predictions of cADSs regulating the central pathways in P. pastoris. The presented methodology with the tools Sc-cADSs-CP and ScanAlgo enables the prediction of master TFs and cADSs in yeasts.

Heterozygous FMN2 missense variant found in a family case of premature ovarian insufficiency

BackgroundPremature ovarian insufficiency (POI) plagues 1% of women under 40, while quite a few remain an unknown cause. The development of sequencing has helped find pathogenic genes and reveal the relationship between DNA repair and ovarian reserve. Through the exome sequencing, our study targets screening out the possible POI pathogenic gene and variants in a Chinese family and 20 sporadic POI patients, preliminarily exploring the functional impact and finding out potential linkages between the gene and POI.ResultsThe whole exome sequencing suggested a novel FMN2 heterozygous variant c.1949C > T (p.Ser650Leu) carried by all three patients in a Chinese family and another c.1967G > A(p.Arg656His) variant in a sporadic case. Since no FMN2 missense mutation is reported for causing human POI, we preliminarily assessed p.Ser650Leu variant via cross-species alignment and 3D modeling and found it possibly deleterious. A series of functional evidence was consistent with our hypothesis. We proved the expression of FMN2 in different stages of oocytes and observed a statistical difference of chromosomal breakages between the POI patient carrying p.Arg656His variant and the health control (p = 0.0013). Western Blot also suggested a decrease in FMN2 and P21 in the mutant type and an associated increase in H2AX. The p.Arg656His variant with an extremely low frequency also indicated that the gene FMN2 might play an essential role in the genetic etiology of POI. To the best of our knowledge, this is the first POI report on missense variants of FMN2.ConclusionThis finding indicates a novel gene possibly related to POI and sheds lights on the study of FMN2.

Exploiting genomic synteny in Felidae: cross-species genome alignments and SNV discovery can aid conservation management

BackgroundWhile recent advances in genomics has enabled vast improvements in the quantification of genome-wide diversity and the identification of adaptive and deleterious alleles in model species, wildlife and non-model species have largely not reaped the same benefits. This has been attributed to the resources and infrastructure required to develop essential genomic datasets such as reference genomes. In the absence of a high-quality reference genome, cross-species alignments can provide reliable, cost-effective methods for single nucleotide variant (SNV) discovery. Here, we demonstrated the utility of cross-species genome alignment methods in gaining insights into population structure and functional genomic features in cheetah (Acinonyx jubatas), snow leopard (Panthera uncia) and Sumatran tiger (Panthera tigris sumatrae), relative to the domestic cat (Felis catus).ResultsAlignment of big cats to the domestic cat reference assembly yielded nearly complete sequence coverage of the reference genome. From this, 38,839,061 variants in cheetah, 15,504,143 in snow leopard and 13,414,953 in Sumatran tiger were discovered and annotated. This method was able to delineate population structure but limited in its ability to adequately detect rare variants. Enrichment analysis of fixed and species-specific SNVs revealed insights into adaptive traits, evolutionary history and the pathogenesis of heritable diseases.ConclusionsThe high degree of synteny among felid genomes enabled the successful application of the domestic cat reference in high-quality SNV detection. The datasets presented here provide a useful resource for future studies into population dynamics, evolutionary history and genetic and disease management of big cats. This cross-species method of variant discovery provides genomic context for identifying annotated gene regions essential to understanding adaptive and deleterious variants that can improve conservation outcomes.

Cross-species functional alignment reveals evolutionary hierarchy within the connectome

Evolution provides an important window into how cortical organization shapes function and vice versa. The complex mosaic of changes in brain morphology and functional organization that have shaped the mammalian cortex during evolution, complicates attempts to chart cortical differences across species. It limits our ability to fully appreciate how evolution has shaped our brain, especially in systems associated with unique human cognitive capabilities that lack anatomical homologues in other species. Here, we develop a function-based method for cross-species alignment that enables the quantification of homologous regions between humans and rhesus macaques, even when their location is decoupled from anatomical landmarks. Critically, we find cross-species similarity in functional organization reflects a gradient of evolutionary change that decreases from unimodal systems and culminates with the most pronounced changes in posterior regions of the default mode network (angular gyrus, posterior cingulate and middle temporal cortices). Our findings suggest that the establishment of the default mode network, as the apex of a cognitive hierarchy, has changed in a complex manner during human evolution – even within subnetworks.

Transcriptomic evidence that von Economo neurons are regionally specialized extratelencephalic-projecting excitatory neurons

von Economo neurons (VENs) are bipolar, spindle-shaped neurons restricted to layer 5 of human frontoinsula and anterior cingulate cortex that appear to be selectively vulnerable to neuropsychiatric and neurodegenerative diseases, although little is known about other VEN cellular phenotypes. Single nucleus RNA-sequencing of frontoinsula layer 5 identifies a transcriptomically-defined cell cluster that contained VENs, but also fork cells and a subset of pyramidal neurons. Cross-species alignment of this cell cluster with a well-annotated mouse classification shows strong homology to extratelencephalic (ET) excitatory neurons that project to subcerebral targets. This cluster also shows strong homology to a putative ET cluster in human temporal cortex, but with a strikingly specific regional signature. Together these results suggest that VENs are a regionally distinctive type of ET neuron. Additionally, we describe the first patch clamp recordings of VENs from neurosurgically-resected tissue that show distinctive intrinsic membrane properties relative to neighboring pyramidal neurons.

Few SINEs of life: Alu elements have little evidence for biological relevance despite elevated translation.

Transposable elements colonize genomes and with time may end up being incorporated into functional regions. SINE Alu elements, which appeared in the primate lineage, are ubiquitous in the human genome and more than a thousand overlap annotated coding exons. Although almost all Alu-derived coding exons appear to be in alternative transcripts, they have been incorporated into the main coding transcript in at least 11 genes. The extent to which Alu regions are incorporated into functional proteins is unclear, but we detected reliable peptide evidence to support the translation to protein of 33 Alu-derived exons. All but one of the Alu elements for which we detected peptides were frame-preserving and there was proportionally seven times more peptide evidence for Alu elements as for other primate exons. Despite this strong evidence for translation to protein we found no evidence of selection, either from cross species alignments or human population variation data, among these Alu-derived exons. Overall, our results confirm that SINE Alu elements have contributed to the expansion of the human proteome, and this contribution appears to be stronger than might be expected over such a relatively short evolutionary timeframe. Despite this, the biological relevance of these modifications remains open to question.

Biological process activity transformation of single cell gene expression for cross-species alignment

The maintenance and transition of cellular states are controlled by biological processes. Here we present a gene set-based transformation of single cell RNA-Seq data into biological process activities that provides a robust description of cellular states. Moreover, as these activities represent species-independent descriptors, they facilitate the alignment of single cell states across different organisms.

P013 Development of the IPD-MHC database

The IPD-MHC Database was released in 2003 to provide a curated database of MHC sequences from a number of non-human species. The system was modeled on the IPD-IMGT/HLA database but expanded to cover multiple species and nomenclature systems. The initial release contained data on non-human primates, felines and canines. Since this time further data on horse, sheep, swine, cattle, fish and rats has been added. Since its release, the database has grown rapidly and now contains nearly 7000 alleles covering over 70 non-human species. The site averages ∼ 1500 unique visitors a month viewing ∼ 5,000 pages. This growth has led to a number of challenges in maintaining and developing the site. The bioinformatics challenges of running a locus specific database across multiple nomenclature systems and locations have impacted on the performance and sustainability of the current IPD-MHC model. In 2015 a BBSRC Bioinformatics and Biological Resources (BBR) grant (UK) was awarded to fund both essential upgrades and future expansion of IPD-MHC. Work is currently focused on the underlying database, the public website and the submission procedures. The key aims of this first phase are development of a universal cross species data submission and display tool and streamlining and standardising the work of the nomenclature committees in curating the data. This will allow for simpler and more frequent species updates. Future developments aim to incorporate cross-species alignment, primer design tools and the incorporation of NGS data. Ultimately this will create an improved system capable of coping with large number of sequence variants, genes and species within the dataset. The project is overseen by a steering committee representing key stakeholders and non-human MHC nomenclature committees. This group is also helping to develop new MHC nomenclature standards and guidelines. IPD-MHC: https://www.ebi.ac.uk/ipd/mhc/ .

A Novel Missense Mutation of GATA4 in a Chinese Family with Congenital Heart Disease.

BackgroundCongenital heart disease (CHD) is the most prevalent type of birth defect in human, with high morbidity in infant. Several genes essential for heart development have been identified. GATA4 is a pivotal transcription factor that can regulate the cardiac development. Many GATA4 mutations have been identified in patients with different types of CHD.AimsIn this study, the NKX2-5, HAND1 and GATA4 coding regions were sequenced in a family spanning three generations in which seven patients had CHD. Disease-causing potential variation in this family was evaluated by bioinformatics programs and the transcriptional activity of mutant protein was analyzed by the dual luciferase reporter assay.ResultsA novel GATA4 mutation, c.C931T (p.R311W), was identified and co-segregated with the affected patients in this family. The bioinformatics programs predicted this heterozygous mutation to be deleterious and the cross-species alignment of GATA4 sequences showed that the mutation occurred within a highly conserved amino acid. Even though it resided in the nuclear localization signal domain, the mutant protein didn’t alter its intracellular distribution. Nevertheless, further luciferase reporter assay demonstrated that the p.R311W mutation reduced the ability of GATA4 to activate its downstream target gene.ConclusionsOur study identified a novel mutation in GATA4 that likely contributed to the CHD in this family. This finding expanded the spectrum of GATA4 mutations and underscored the pathogenic correlation between GATA4 mutations and CHD.

ConBind: motif-aware cross-species alignment for the identification of functional transcription factor binding sites.

Eukaryotic gene expression is regulated by transcription factors (TFs) binding to promoter as well as distal enhancers. TFs recognize short, but specific binding sites (TFBSs) that are located within the promoter and enhancer regions. Functionally relevant TFBSs are often highly conserved during evolution leaving a strong phylogenetic signal. While multiple sequence alignment (MSA) is a potent tool to detect the phylogenetic signal, the current MSA implementations are optimized to align the maximum number of identical nucleotides. This approach might result in the omission of conserved motifs that contain interchangeable nucleotides such as the ETS motif (IUPAC code: GGAW). Here, we introduce ConBind, a novel method to enhance alignment of short motifs, even if their mutual sequence similarity is only partial. ConBind improves the identification of conserved TFBSs by improving the alignment accuracy of TFBS families within orthologous DNA sequences. Functional validation of the Gfi1b + 13 enhancer reveals that ConBind identifies additional functionally important ETS binding sites that were missed by all other tested alignment tools. In addition to the analysis of known regulatory regions, our web tool is useful for the analysis of TFBSs on so far unknown DNA regions identified through ChIP-sequencing.

Sperm Lysozyme-Like Protein 1 (SLLP1), an intra-acrosomal oolemmal-binding sperm protein, reveals filamentous organization in protein crystal form.

Sperm lysozyme-like protein 1 (SLLP1) is one of the lysozyme-like proteins predominantly expressed in mammalian testes that lacks bacteriolytic activity, localizes in the sperm acrosome, and exhibits high affinity for an oolemmal receptor, SAS1B. The crystal structure of mouse SLLP1 (mSLLP1) was determined at 2.15Å resolution. mSLLP1 monomer adopts a structural fold similar to that of chicken/mouse lysozymes retaining all four canonical disulfide bonds. mSLLP1 is distinct from c-lysozyme by substituting two essential catalytic residues (E35T/D52N), exhibiting different surface charge distribution, and by forming helical filaments approximately 75Å in diameter with a 25Å central pore comprised of six monomers per helix turn repeating every 33Å. Cross-species alignment of all reported SLLP1 sequences revealed a set of invariant surface regions comprising a characteristic fingerprint uniquely identifying SLLP1 from other c-lysozyme family members. The fingerprint surface regions reside around the lips of the putative glycan-binding groove including three polar residues (Y33/E46/H113). A flexible salt bridge (E46-R61) was observed covering the glycan-binding groove. The conservation of these regions may be linked to their involvement in oolemmal protein binding. Interaction between SLLP1 monomer and its oolemmal receptor SAS1B was modeled using protein-protein docking algorithms, utilizing the SLLP1 fingerprint regions along with the SAS1B conserved surface regions. This computational model revealed complementarity between the conserved SLLP1/SAS1B interacting surfaces supporting the experimentally observed SLLP1/SAS1B interaction involved in fertilization.

Annotation of the Protein Coding Regions of the Equine Genome.

Current gene annotation of the horse genome is largely derived from in silico predictions and cross-species alignments. Only a small number of genes are annotated based on equine EST and mRNA sequences. To expand the number of equine genes annotated from equine experimental evidence, we sequenced mRNA from a pool of forty-three different tissues. From these, we derived the structures of 68,594 transcripts. In addition, we identified 301,829 positions with SNPs or small indels within these transcripts relative to EquCab2. Interestingly, 780 variants extend the open reading frame of the transcript and appear to be small errors in the equine reference genome, since they are also identified as homozygous variants by genomic DNA resequencing of the reference horse. Taken together, we provide a resource of equine mRNA structures and protein coding variants that will enhance equine and cross-species transcriptional and genomic comparisons.

Harnessing cross-species alignment to discover SNPs and generate a draft genome sequence of a bighorn sheep (Ovis canadensis).

BackgroundWhole genome sequences (WGS) have proliferated as sequencing technology continues to improve and costs decline. While many WGS of model or domestic organisms have been produced, a growing number of non-model species are also being sequenced. In the absence of a reference, construction of a genome sequence necessitates de novo assembly which may be beyond the ability of many labs due to the large volumes of raw sequence data and extensive bioinformatics required. In contrast, the presence of a reference WGS allows for alignment which is more tractable than assembly. Recent work has highlighted that the reference need not come from the same species, potentially enabling a wide array of species WGS to be constructed using cross-species alignment. Here we report on the creation a draft WGS from a single bighorn sheep (Ovis canadensis) using alignment to the closely related domestic sheep (Ovis aries).ResultsTwo sequencing libraries on SOLiD platforms yielded over 865 million reads, and combined alignment to the domestic sheep reference resulted in a nearly complete sequence (95% coverage of the reference) at an average of 12x read depth (104 SD). From this we discovered over 15 million variants and annotated them relative to the domestic sheep reference. We then conducted an enrichment analysis of those SNPs showing fixed differences between the reference and sequenced individual and found significant differences in a number of gene ontology (GO) terms, including those associated with reproduction, muscle properties, and bone deposition.ConclusionOur results demonstrate that cross-species alignment enables the creation of novel WGS for non-model organisms. The bighorn sheep WGS will provide a resource for future resequencing studies or comparative genomics.Electronic supplementary materialThe online version of this article (doi:10.1186/s12864-015-1618-x) contains supplementary material, which is available to authorized users.

APPRIS WebServer and WebServices.

This paper introduces the APPRIS WebServer (http://appris.bioinfo.cnio.es) and WebServices (http://apprisws.bioinfo.cnio.es). Both the web servers and the web services are based around the APPRIS Database, a database that presently houses annotations of splice isoforms for five different vertebrate genomes. The APPRIS WebServer and WebServices provide access to the computational methods implemented in the APPRIS Database, while the APPRIS WebServices also allows retrieval of the annotations. The APPRIS WebServer and WebServices annotate splice isoforms with protein structural and functional features, and with data from cross-species alignments. In addition they can use the annotations of structure, function and conservation to select a single reference isoform for each protein-coding gene (the principal protein isoform). APPRIS principal isoforms have been shown to agree overwhelmingly with the main protein isoform detected in proteomics experiments. The APPRIS WebServer allows for the annotation of splice isoforms for individual genes, and provides a range of visual representations and tools to allow researchers to identify the likely effect of splicing events. The APPRIS WebServices permit users to generate annotations automatically in high throughput mode and to interrogate the annotations in the APPRIS Database. The APPRIS WebServices have been implemented using REST architecture to be flexible, modular and automatic.

Tracing the evolution of lineage-specific transcription factor binding sites in a birth-death framework.

Changes in cis-regulatory element composition that result in novel patterns of gene expression are thought to be a major contributor to the evolution of lineage-specific traits. Although transcription factor binding events show substantial variation across species, most computational approaches to study regulatory elements focus primarily upon highly conserved sites, and rely heavily upon multiple sequence alignments. However, sequence conservation based approaches have limited ability to detect lineage-specific elements that could contribute to species-specific traits. In this paper, we describe a novel framework that utilizes a birth-death model to trace the evolution of lineage-specific binding sites without relying on detailed base-by-base cross-species alignments. Our model was applied to analyze the evolution of binding sites based on the ChIP-seq data for six transcription factors (GATA1, SOX2, CTCF, MYC, MAX, ETS1) along the lineage toward human after human-mouse common ancestor. We estimate that a substantial fraction of binding sites (∼58–79% for each factor) in humans have origins since the divergence with mouse. Over 15% of all binding sites are unique to hominids. Such elements are often enriched near genes associated with specific pathways, and harbor more common SNPs than older binding sites in the human genome. These results support the ability of our method to identify lineage-specific regulatory elements and help understand their roles in shaping variation in gene regulation across species.

Predicting Long Non-coding RNAs Based on Genomic Sequence Information

The binary classification of coding and non-coding genes is simplified near to 50 years. Genome-wide transcriptome studies have revealed that there exist tens of thousands of long non-coding RNAs (lncRNAs), while the functions are being uncovered slowly. Accurate identification of lncRNAs is the initial step to the systematic characterization of lncRNAs. The diversity of transcription patterns for lncRNAs challenges the available non-coding RNA prediction algorithms. Until now, prediction of lncRNAs mostly relies on genomic sequence and cross-species alignment information. Here, we introduce the main strategies that can discriminate lncRNA from protein-coding transcripts. Especially, recently available machine learning algorithms are shown efficient to the rapid and accurate identification of lncRNAs from a large number of putative lncRNAs based on transcriptome assembled transcripts, which would provide the basis of understanding of lncRNA biology.

Detection and clinical value of novel GATA6 mutations associated with congenital atrial septal defects

Objective To identify novel mutations in the GATA6 gene associated with congenital atrial septal defects (ASD). Methods This was a case-control study.A cohort of 220 unrelated Han-race patients with congenital ASD and 200 unrelated ethnically matched healthy individuals used as controls, who were admitted to Tongji University Affiliated Tongji Hospital from January, 2007 to October, 2011, were recruited.The peripheral venous blood samples from the participants were prepared.All the coding exons and their flanking sequences of the GATA6 gene were amplified by polymerase chain reaction and sequenced using the di-deoxynucleotide chain termination technique.The acquired sequences were aligned with the sequences derived from GenBank by BLAST to identify the sequence variations.The software ClustalW was used to analyze the conservation of the altered amino acids. Results Three novel heterozygous missense GATA6 mutations, c.250G>A (p.A84T), c.649G>C (p.G217R) and c.1270A>C (p.S424R), were identified in 3 of 220 ASD patients, respectively.None of the three mutations was detected in 200 healthy control individuals.A cross-species alignment of GATA6 encoded protein sequences showed that the mutated amino acids were relatively conserved evolutionarily. Conclusion The identification of novel GATA6 mutations associated with ASD contributes to the reveal of the mechanism involved in the pathogenesis of ASD.（Chin J Lab Med,2012,35:1108-1111） Key words: Heart septal defects; atrial; GATA6 transcription factor; Mutation