Search

Background There has been a rapid increase in the brain size relative to body size during mammalian evolutionary history. In particular, the enlarged and globular brain is the most distinctive anatomical feature of modern humans that set us apart from other extinct and extant primate species. Genetic basis of large brain size in modern humans has largely remained enigmatic. Genes associated with the pathological reduction of brain size (primary microcephaly-MCPH) have the characteristics and functions to be considered ideal candidates to unravel the genetic basis of evolutionary enlargement of human brain size. For instance, the brain size of microcephaly patients is similar to the brain size of Pan troglodyte and the very early hominids like the Sahelanthropus tchadensis and Australopithecus afarensis. Results The present study investigates the molecular evolutionary history of subset of autosomal recessive primary microcephaly (MCPH) genes; CEP135, ZNF335, PHC1, SASS6, CDK6, MFSD2A, CIT, and KIF14 across 48 mammalian species. Codon based substitutions site analysis indicated that ZNF335, SASS6, CIT, and KIF14 have experienced positive selection in eutherian evolutionary history. Estimation of divergent selection pressure revealed that almost all of the MCPH genes analyzed in the present study have maintained their functions throughout the history of placental mammals. Contrary to our expectations, human-specific adoptive evolution was not detected for any of the MCPH genes analyzed in the present study. Conclusion Based on these data it can be inferred that protein-coding sequence of MCPH genes might not be the sole determinant of increase in relative brain size during primate evolutionary history.

The zinc finger-containing transcription factor Gli3 is a key mediator of Hedgehog (Hh) signaling pathway. In vertebrates, Gli3 has widespread expression pattern during early embryonic development. Along the anteroposterior axes of the central nervous system (CNS), dorsoventral neural pattern elaboration is achieved through Hh mediated spatio-temporal deployment of Gli3 transcripts. Previously, we and others uncovered a set of enhancers that mediate many of the known aspects of Gli3 expression during neurogenesis. However, the potential role of Gli3 associated enhancers in trait evolution has not yet received any significant attention. Here, we investigate the evolutionary patterns of Gli3 associated CNS-specific enhancers that have been reported so far. A subset of these enhancers has undergone an accelerated rate of molecular evolution in the human lineage in comparison to other primates/mammals. These fast-evolving enhancers have acquired human-specific changes in transcription factor binding sites (TFBSs). These human-unique changes within subset of Gli3 associated CNS-specific enhancers were further validated as single nucleotide polymorphisms through 1000 Genome Project Phase 3 data. This work not only infers the molecular evolutionary patterns of Gli3 associated enhancers but also provides clues for putative genetic basis of the population-specificity of gene expression regulation.

Background The zinc finger-containing transcription factor Gli2, is a key mediator of Hedgehog (Hh) signaling and participates in embryonic patterning of various organs including the central nervous system (CNS) and limbs. Abnormal expression of Gli2 can impede the transcription of Hh target genes through disruption of proper balance between Gli2 and Gli3 functions. Therefore, delineation of enhancers that are required for complementary roles of Glis would allow the interrogation of those pathogenic variants that cause gene dysregulation, and a corresponding abnormal phenotype. Previously, we reported tissue-specific enhancers for Gli family including Gli2 through direct tetrapod-teleost comparisons. Results Here, we employed the sequence alignments of slowly evolving spotted gar and elephant shark and have identified six novel conserved noncoding elements in human GLI2 containing locus. Zebrafish-based transgenic assays revealed that combined action of these autonomous CNEs reflects many aspects of Gli2 specific endogenous transcriptional activity, including CNS and pectoral fins. Conclusion Taken together with our previous findings, this study suggests that Hh-signaling controlled deployment of Gli2 activity in embryonic patterning arose in the common ancestor of gnathostomes. These GLI2 specific cis-regulatory modules will help to identify DNA variants that probably reside outside of coding intervals and are associated with congenital anomalies.

Background Parkinson’s disease (PD) is the second most common neurodegenerative disorder. PD associated human UCHL1 (Ubiquitin C-terminal hydrolase L1) gene belongs to the family of deubiquitinases and is known to be highly expressed in neurons (1–2% in soluble form). Several functions of UCHL1 have been proposed including ubiquitin hydrolyze activity, ubiquitin ligase activity and stabilization of the mono-ubiquitin. Mutations in human UCHL1 gene have been associated with PD and other neurodegenerative disorders. The present study aims to decipher the sequence evolutionary pattern and structural dynamics of UCHL1. Furthermore, structural and interactional analysis of UCHL1 was performed to help elucidate the pathogenesis of PD. Results The phylogenetic tree topology suggests that the UCHL1 gene had originated in early gnathostome evolutionary history. Evolutionary rate analysis of orthologous sequences reveals strong purifying selection on UCHL1. Comparative structural analysis of UCHL1 pinpoints an important protein segment spanning amino acid residues 32 to 39 within secretion site with crucial implications in evolution and PD pathogenesis through a well known phenomenon called intragenic epistasis. Identified critical protein segment appears to play an indispensable role in protein stability, proper protein conformation as well as harboring critical interaction sites. Conclusions Conclusively, the critical protein segment of UCHL1 identified in the present study not only demonstrates the relevant role of intraprotein conformational epistasis in the pathophysiology of PD but also offers a novel therapeutic target for the disease.

The outbreak of coronavirus disease 2019 (COVID-19) caused by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is spreading fast worldwide. There is a pressing need to understand how the virus counteracts host innate immune responses. Deleterious clinical manifestations of coronaviruses have been associated with virus-induced direct dysregulation of innate immune responses occurring via viral macrodomains located within nonstructural protein-3 (Nsp3). However, no substantial information is available concerning the relationship of macrodomains to the unusually high pathogenicity of SARS-CoV-2. Here, we show that structural evolution of macrodomains may impart a critical role to the unique pathogenicity of SARS-CoV-2. Using sequence, structural, and phylogenetic analysis, we identify a specific set of historical substitutions that recapitulate the evolution of the macrodomains that counteract host immune response. These evolutionary substitutions may alter and reposition the secondary structural elements to create new intra-protein contacts and, thereby, may enhance the ability of SARS-CoV-2 to inhibit host immunity. Further, we find that the unusual virulence of this virus is potentially the consequence of Darwinian selection‐driven epistasis in protein evolution. Our findings warrant further characterization of macrodomain-specific evolutionary substitutions in in vitro and in vivo models to determine their inhibitory effects on the host immune system.

The BIG Data Center at Beijing Institute of Genomics (BIG) of the Chinese Academy of Sciences provides a suite of database resources in support of worldwide research activities in both academia and industry. With the vast amounts of multi-omics data generated at unprecedented scales and rates, the BIG Data Center is continually expanding, updating and enriching its core database resources through big data integration and value-added curation. Resources with significant updates in the past year include BioProject (a biological project library), BioSample (a biological sample library), Genome Sequence Archive (GSA, a data repository for archiving raw sequence reads), Genome Warehouse (GWH, a centralized resource housing genome-scale data), Genome Variation Map (GVM, a public repository of genome variations), Science Wikis (a catalog of biological knowledge wikis for community annotations) and IC4R (Information Commons for Rice). Newly released resources include EWAS Atlas (a knowledgebase of epigenome-wide association studies), iDog (an integrated omics data resource for dog) and RNA editing resources (for editome-disease associations and plant RNA editosome, respectively). To promote biodiversity and health big data sharing around the world, the Open Biodiversity and Health Big Data (BHBD) initiative is introduced. All of these resources are publicly accessible at http://bigd.big.ac.cn.

Empirical assessments of human accelerated noncoding DNA frgaments have delineated presence of many cis-regulatory elements. Enhancers make up an important category of such accelerated cis-regulatory elements that efficiently control the spatiotemporal expression of many developmental genes. Establishing plausible reasons for accelerated enhancer sequence divergence in Homo sapiens has been termed significant in various previously published studies. This acceleration by including closely related primates and archaic human data has the potential to open up evolutionary avenues for deducing present-day brain structure. This study relied on empirically confirmed brain exclusive enhancers to avoid any misjudgments about their regulatory status and categorized among them a subset of enhancers with an exceptionally accelerated rate of lineage specific divergence in humans. In this assorted set, 13 distinct transcription factor binding sites were located that possessed unique existence in humans. Three of 13 such sites belonging to transcription factors SOX2, RUNX1/3, and FOS/JUND possessed single nucleotide variants that made them unique to H. sapiens upon comparisons with Neandertal and Denisovan orthologous sequences. These variants modifying the binding sites in modern human lineage were further substantiated as single nucleotide polymorphisms via exploiting 1000 Genomes Project Phase3 data. Long range haplotype based tests laid out evidence of positive selection to be governing in African population on two of the modern human motif modifying alleles with strongest results for SOX2 binding site. In sum, our study acknowledges acceleration in noncoding regulatory landscape of the genome and highlights functional parts within it to have undergone accelerated divergence in present-day human population.

The vertebrate genome contains several closely spaced sets of paralogous genes from distinct gene families on typically two, three or four different chromosomes (paralogons). These four fold paralogy regions have been considered as historical remnants of whole genome duplication events (WGDs/2R hypothesis). To examine the 2R hypothesis, a robust phylogenetic analysis of 40 multigene families with triplicated or quadruplicated distribution on human MHC bearing chromosomes (1/6/9/19) was conducted. Topology comparison approach categorized the members of 40 families into six distinct co-duplicated groups. Genes belonging to a particular co-duplicated group are duplicated concurrently, whereas genes of two different co-duplicated groups do not share their evolutionary history and have not duplicated in harmony. Our results based on this large scale phylogenetic data set contradict the polyploidization model and are indicative of small-scale duplications and rearrangement events that cover the entire span of animal history.

The National Genomics Data Center (NGDC) provides a suite of database resources to support worldwide research activities in both academia and industry. With the rapid advancements in higher-throughput and lower-cost sequencing technologies and accordingly the huge volume of multi-omics data generated at exponential scales and rates, NGDC is continually expanding, updating and enriching its core database resources through big data integration and value-added curation. In the past year, efforts for update have been mainly devoted to BioProject, BioSample, GSA, GWH, GVM, NONCODE, LncBook, EWAS Atlas and IC4R. Newly released resources include three human genome databases (PGG.SNV, PGG.Han and CGVD), eLMSG, EWAS Data Hub, GWAS Atlas, iSheep and PADS Arsenal. In addition, four web services, namely, eGPS Cloud, BIG Search, BIG Submission and BIG SSO, have been significantly improved and enhanced. All of these resources along with their services are publicly accessible at https://bigd.big.ac.cn.

In recent years, the number of human long noncoding RNAs (lncRNAs) that have been identified has increased exponentially. However, these lncRNAs are poorly annotated compared to protein-coding genes, posing great challenges for a better understanding of their functional significance and elucidating their complex functioning molecular mechanisms. Here we employ both community and expert curation to yield a comprehensive collection of human lncRNAs and their annotations. Specifically, LncRNAWiki (http://lncrna.big.ac.cn/index.php/Main_Page) uses a wiki-based community curation model, thus showing great promise in dealing with the flood of biological knowledge, while LncBook (http://bigd.big.ac.cn/lncbook) is an expert curation-based database that provides a complement to LncRNAWiki. LncBook features a comprehensive collection of human lncRNAs and a systematic curation of lncRNAs by multi-omics data integration, functional annotation, and disease association. These protocols provide step-by-step instructions on how to browse and search a specific lncRNA and how to obtain a range of related information including expression, methylation, variation, function, and disease association. © 2019 by John Wiley & Sons, Inc.

Background The hypothesis that vertebrates have experienced two ancient, whole genome duplications (WGDs) is of central interest to evolutionary biology and has been implicated in evolution of developmental complexity. Three-way and Four-way paralogy regions in human and other vertebrate genomes are considered as vital evidence to support this hypothesis. Alternatively, it has been proposed that such paralogy regions are created by small-scale duplications that occurred at different intervals over the evolution of life. Results To address this debate, the present study investigates the evolutionary history of multigene families with at least three-fold representation on human chromosomes 1, 2, 8 and 20. Phylogenetic analysis and the tree topology comparisons classified the members of 36 multigene families into four distinct co-duplicated groups. Gene families falling within the same co-duplicated group might have duplicated together, whereas genes belong to different co-duplicated groups might have distinct evolutionary origins. Conclusion Taken together with previous investigations, the current study yielded no proof in favor of WGDs hypothesis. Rather, it appears that the vertebrate genome evolved as a result of small-scale duplication events, that cover the entire span of the animals’ history. Electronic supplementary material The online version of this article (10.1186/s12862-019-1441-0) contains supplementary material, which is available to authorized users.

Background Frizzled family members belong to G-protein coupled receptors and encode proteins accountable for cell signal transduction, cell proliferation and cell death. Members of Frizzled receptor family are considered to have critical roles in causing various forms of cancer, cardiac hypertrophy, familial exudative vitreoretinopathy (FEVR) and schizophrenia. Results This study investigates the evolutionary and structural aspects of Frizzled receptors, with particular focus on FEVR associated FZD4 gene. The phylogenetic tree topology suggests the diversification of Frizzled receptors at the root of metazoans history. Moreover, comparative structural data reveals that FEVR associated missense mutations in FZD4 effect the common protein region (amino acids 495–537) through a well-known phenomenon called epistasis. This critical protein region is present at the carboxyl-terminal domain and encompasses the K-T/S-XXX-W, a PDZ binding motif and S/T-X-V PDZ recognition motif. Conclusion Taken together these results demonstrate that during the course of evolution, FZD4 has acquired new functions or epistasis via complex patter of gene duplications, sequence divergence and conformational remodeling. In particular, amino acids 495–537 at the C-terminus region of FZD4 protein might be crucial in its normal function and/or pathophysiology. This critical region of FZD4 protein may offer opportunities for the development of novel therapeutics approaches for human retinal vascular disease. Electronic supplementary material The online version of this article (10.1186/s12862-019-1400-9) contains supplementary material, which is available to authorized users.

The human nose has a diversified morphology among different geographical distributions and serves as a natural conditioning system for the inhaled air. Several studies have put forth statistically sound observations of the presence of large-bulbous noses in hot and humid Africa and narrower-aquiline noses in cold and dry Europe. These diverse nasal shapes have been termed as local adaptations to a region's climate. Over the past couple of decades, several GWA studies have linked potential loci with several diseases and morphological traits including the human nasal forms. In this study, highly significant SNPs associated with eight nasal morphological traits from previously reported GWA studies were collected. SNP data was obtained from the 1000 Genomes Project Phase III. By employing tests such as extended haplotype homozygosity (EHH), haplotype bifurcation diagrams and Fst measure, a subset of 6 nasal-associated SNPs were found to be significantly differentiated among the human sub-populations indicating that positive selection has acted on human nasal-associated loci.

To elucidate important cellular and molecular interactions that regulate patterning and skeletal development, vertebrate limbs served as a model organ. A growing body of evidence from detailed studies on a subset of limb regulators like the HOXD cluster or SHH, reveals the importance of enhancers in limb related developmental and disease processes. Exploiting the recent genome-wide availability of functionally confirmed enhancer dataset, this study establishes regulatory interactions for dozens of human limb developmental genes. From these data, it appears that the long-range regulatory interactions are fairly common during limb development. This observation highlights the significance of chromosomal breaks/translocations in human limb deformities. Transcriptional factor (TF) analysis predicts that the differentiation of early nascent limb-bud into future territories entail distinct TF interaction networks. Conclusively, an important motivation for annotating the human limb specific regulatory networks is to pave way for the systematic exploration of their role in disease and evolution.

Human brain evolution is characterized by dramatic expansion in cerebral cortex size. WDR62 (WD repeat domain 62) is one of the important gene in controlling human cortical development. Mutations in WDR62 lead to primary microcephaly, a neurodevelopmental disease characterized by three to four fold reduction in cerebral cortex size of affected individuals. This study analyzes comparative protein evolutionary rate to provide a useful insight into the molecular evolution of WDR62 and hence pinpointed human specific amino acid replacements. Comparative analysis of human WDR62 with two archaic humans (Neanderthals and Denisovans) and modern human populations revealed that five hominin specific amino acid residues (human specific amino acids shared with two archaic humans) might have been accumulated in the common ancestor of extinct archaic humans and modern humans about 550,000–765,000 years ago. Collectively, the data demonstrates an acceleration of WDR62 sequence evolution in hominin lineage and suggests that the ability of WDR62 protein to mediate the neurogenesis has been altered in the course of hominin evolution., Graphical abstract Image 1, Highlights • We trace the evolutionary history of WDR62 and its putative paralogs. • We identify accelerated sequence evolution in human WDR62. • We pinpoint eight human specific amino acid sites that reside on the C-terminal. • Out of eight, six sites are shared with archaic humans.

Motivation The significance of long non-coding RNAs (lncRNAs) in many biological processes and diseases has gained intense interests over the past several years. However, computational identification of lncRNAs in a wide range of species remains challenging; it requires prior knowledge of well-established sequences and annotations or species-specific training data, but the reality is that only a limited number of species have high-quality sequences and annotations. Results Here we first characterize lncRNAs in contrast to protein-coding RNAs based on feature relationship and find that the feature relationship between open reading frame length and guanine-cytosine (GC) content presents universally substantial divergence in lncRNAs and protein-coding RNAs, as observed in a broad variety of species. Based on the feature relationship, accordingly, we further present LGC, a novel algorithm for identifying lncRNAs that is able to accurately distinguish lncRNAs from protein-coding RNAs in a cross-species manner without any prior knowledge. As validated on large-scale empirical datasets, comparative results show that LGC outperforms existing algorithms by achieving higher accuracy, well-balanced sensitivity and specificity, and is robustly effective (>90% accuracy) in discriminating lncRNAs from protein-coding RNAs across diverse species that range from plants to mammals. To our knowledge, this study, for the first time, differentially characterizes lncRNAs and protein-coding RNAs based on feature relationship, which is further applied in computational identification of lncRNAs. Taken together, our study represents a significant advance in characterization and identification of lncRNAs and LGC thus bears broad potential utility for computational analysis of lncRNAs in a wide range of species. Availability and implementation LGC web server is publicly available at http://bigd.big.ac.cn/lgc/calculator. The scripts and data can be downloaded at http://bigd.big.ac.cn/biocode/tools/BT000004. Supplementary information Supplementary data are available at Bioinformatics online.

HOX genes encode transcriptional factors that play a pivotal role in specifying regional identity in nearly every bilateral animal. The birth of HOX gene cluster and its subsequent evolution, either in regulation or function, underlie the evolution of many bilaterian features and hence to the evolutionary radiation of this group. Despite of this importance, evolution of HOX cluster in vertebrates remains largely obscure because the phylogenetic history of these genes is poorly resolved. This has led to the controversy about whether four HOX clusters in human originated through two rounds (2R) of whole-genome duplications or instead evolved by small-scale events early in vertebrate evolution. Recently, the large-scale phylogenetic analysis of triplicate and quadruplicate paralogous regions residing on human HOX-bearing chromosomes provided an unprecedented insight into events that shaped vertebrate genome early in their history. Based on these data and comparative genomic analysis of fruit fly, red floor beetle, and human, this study infers the genic content of minimal HOX locus in the Urbilaterian and reconstructs its duplication history. It appears that four HOX clusters of humans are not remnants of polyploidy events in vertebrate ancestry. Rather, current evidence suggests that one-to-four transition in HOX cluster number occurred by three-step sequential process involving regional duplication events. Therefore, it is concluded that the evolutionary origin of vertebrate novelties, including the complexity of their body, is the consequence of small-scale genetic changes at widely different times over their history.

The zinc-finger transcription factor GLI3 acts as a primary transducer of Sonic hedgehog (Shh) signaling in a context-dependent combinatorial fashion. GLI3 participates in the patterning and growth of many organs, including the central nervous system (CNS) and limbs. Previously, we reported a subset of human intronic cis-regulators controlling many known aspects of endogenous Gli3 expression in mouse and zebrafish. Here we demonstrate in a transgenic zebrafish assay the potential of two novel tetrapod-teleost conserved non-coding elements (CNEs) docking within GLI3 intronic intervals (intron 3 and 4) to induce reporter gene expression at known sites of endogenous Gli3 transcription in embryonic domains such as the central nervous system (CNS) and limbs. Interestingly, the cell culture based assays reveal harmony with the context dependent dual nature of intra-GLI3 conserved elements. Furthermore, a transgenic zebrafish assay of previously reported limb-specific GLI3 transcriptional enhancers (previously tested in mice and chicken limb buds) induced reporter gene expression in zebrafish blood precursor cells and notochord instead of fin. These results demonstrate that the appendage-specific activity of a subset of GLI3-associated enhancers might be a tetrapod innovation. Taken together with our recent data, these results suggest that during the course of vertebrate evolution Gli3 expression control acquired a complex cis-regulatory landscape for spatiotemporal patterning of CNS and limbs. Comparative data from fish and mice suggest that the functional aspects of a subset of these cis-regulators have diverged significantly between these two lineages.

Enhancers lie at the heart of transcriptional and developmental gene regulation. Therefore, changes in enhancer sequences usually disrupt the target gene expression and result in disease phenotypes. Despite the well-established role of enhancers in development and disease, evolutionary sequence studies are lacking. The current study attempts to unravel the puzzle of bony vertebrates’ conserved noncoding elements (CNE) enhancer evolution. Bayesian phylogenetics of enhancer sequences spotlights promising interordinal relationships among placental mammals, proposing a closer relationship between humans and laurasiatherians while placing rodents at the basal position. Clock-based estimates of enhancer evolution provided a dynamic picture of interspecific rate changes across the bony vertebrate lineage. Moreover, coelacanth in the study augmented our appreciation of the vertebrate cis-regulatory evolution during water–land transition. Intriguingly, we observed a pronounced upsurge in enhancer evolution in land-dwelling vertebrates. These novel findings triggered us to further investigate the evolutionary trend of coding as well as CNE nonenhancer repertoires, to highlight the relative evolutionary dynamics of diverse genomic landscapes. Surprisingly, the evolutionary rates of enhancer sequences were clearly at odds with those of the coding and the CNE nonenhancer sequences during vertebrate adaptation to land, with land vertebrates exhibiting significantly reduced rates of coding sequence evolution in comparison to their fast evolving regulatory landscape. The observed variation in tetrapod cis-regulatory elements caused the fine-tuning of associated gene regulatory networks. Therefore, the increased evolutionary rate of tetrapods’ enhancer sequences might be responsible for the variation in developmental regulatory circuits during the process of vertebrate adaptation to land.

Background: GLI2, a zinc finger transcription factor, mediates Sonic hedgehog signaling, a critical pathway in vertebrate embryogenesis. GLI2 has been implicated in diverse set of embryonic developmental processes, including patterning of central nervous system and limbs. In humans, mutations in GLI2 are associated with several developmental defects, including holoprosencephaly and polydactyly. Results: Here, we demonstrate in transient transgenic zebrafish assays, the potential of a subset of tetrapod-teleost conserved non-coding elements (CNEs) residing within human GLI2 intronic intervals to induce reporter gene expression at known regions of endogenous GLI2 transcription. The regulatory activities of these elements are observed in several embryonic domains, including neural tube and pectoral fin. Moreover, our data reveal an overlapping expression profile of duplicated copies of an enhancer during zebrafish evolution. Conclusions: Our data suggest that during vertebrate history GLI2 acquired a high level of complexity in the genetic mechanisms regulating its expression during spatiotemporal patterning of the central nervous system (CNS) and limbs. Developmental Dynamics 244:681–692, 2015. © 2015 Wiley Periodicals, Inc.

Background The vertebrate genome often contains closely spaced set of paralogous genes from distinct gene families on typically two, three or four different chromosomes (paralogons). This type of genome architecture is widely considered to be remnants of whole genome duplication events (WGD/2R). Results Taking advantage of the well-annotated and high-quality human genomic sequence map as well as the ever-increasing accessibility of large-scale genomic sequence data from a diverse range of animal species, we investigated the evolutionary history of potential quadruplicated regions residing on human HOX-cluster bearing chromosomes (chromosomes 2/7/12/17). For this purpose a detailed phylogenetic analysis was performed for those multigene families, including members of at least three of the four HOX-bearing chromosomes. Topology comparison approach categorized the members of 63 families into distinct co-duplicated groups. Distinct gene families belonging to a particular co-duplicated group, exhibit similar evolutionary history and hence have duplicated concurrently, whereas genes of two different co-duplicated groups do not share their history and have not duplicated in concert with each other. Conclusions These results based on large-scale phylogenetic dataset yielded no evidence in favor of polyploidization events; instead it appears that triplicated and quadruplicated genomic segments on the human HOX-bearing chromosomes arose by small-scale duplication events that occurred at widely different time points in animal evolution.

The zinc finger transcription factor Gli3 is an important mediator of Sonic hedgehog (Shh) signaling. During early embryonic development Gli3 participates in patterning and growth of the central nervous system, face, skeleton, limb, tooth and gut. Precise regulation of the temporal and spatial expression of Gli3 is crucial for the proper specification of these structures in mammals and other vertebrates. Previously we reported a set of human intronic cis-regulators controlling almost the entire known repertoire of endogenous Gli3 expression in mouse neural tube and limbs. However, the genetic underpinning of GLI3 expression in other embryonic domains such as craniofacial structures and internal organs remain elusive. Here we demonstrate in a transgenic mice assay the potential of a subset of human/fish conserved non-coding sequences (CNEs) residing within GLI3 intronic intervals to induce reporter gene expression at known regions of endogenous Gli3 transcription in embryonic domains other than central nervous system (CNS) and limbs. Highly specific reporter expression was observed in craniofacial structures, eye, gut, and genitourinary system. Moreover, the comparison of expression patterns directed by these intronic cis-acting regulatory elements in mouse and zebrafish embryos suggests that in accordance with sequence conservation, the target site specificity of a subset of these elements remains preserved among these two lineages. Taken together with our recent investigations, it is proposed here that during vertebrate evolution the Gli3 expression control acquired multiple, independently acting, intronic enhancers for spatiotemporal patterning of CNS, limbs, craniofacial structures and internal organs.

The zinc-finger transcription factor GLI3 acts during vertebrate development in a combinatorial, context-dependent fashion as a primary transducer of sonic hedgehog (SHH) signaling. In humans, mutations affecting this key regulator of development are associated with GLI3-morphopathies, a group of congenital malformations in which forebrain and limb development are preferentially affected. We show that a non-coding element from intron two of GLI3, ultraconserved in mammals and highly conserved in the pufferfish Fugu, is a transcriptional enhancer. In transient transfection assays, it activates reporter gene transcription in human cell cultures expressing endogenous GLI3 but not in GLI3 negative cells. The identified enhancer element is predicted to contain conserved binding sites for transcription factors crucial for developmental steps in which GLI3 is involved. The regulatory potential of this element is conserved and was used to direct tissue-specific expression of a green fluorescent protein reporter gene in zebrafish embryos and of a beta-galactosidase reporter in transgenic mouse embryos. Time, location, and quantity of reporter gene expression are congruent with part of the pattern previously reported for endogenous GLI3 transcription.

After Alzheimer, Parkinson’s disease (PD) is the second most common neurodegenerative disorder. Alpha synuclein (SNCA) is deemed as a major component of Lewy bodies, a neuropathological feature of PD. Five point mutations in SNCA have been reported so far, responsible for autosomal dominant PD. This study aims to decipher evolutionary and structural insights of SNCA by revealing its sequence and structural evolutionary patterns among sarcopterygians and its paralogous counterparts (SNCB and SNCG). Rate analysis detected strong purifying selection on entire synuclein family. Structural dynamics divulges that during the course of sarcopterygian evolutionary history, the region encompassed 32 to 58 of N-terminal domain of SNCA has acquired its critical functional significance through the epistatic influence of the lineage specific substitutions. In sum, these findings provide an evidence that the region from 32 to 58 of N-terminal lipid binding alpha helix domain of SNCA is the most critical region, not only from the evolutionary perspective but also for the stability and the proper conformation of the protein as well as crucial for the disease pathogenesis, harboring critical interaction sites.

Evolution of organismal complexity and origin of novelties during vertebrate history has been widely explored in context of both regulation of gene expression and gene duplication events. Ohno (1970) for the first time put forward the idea of two rounds whole genome duplication events as the most plausible explanation for evolutionarizing the vertebrate lineage (2R hypothesis). To test the validity of 2R hypothesis, a robust phylogenomic analysis of multigene families with triplicated or quadruplicated representation on human FGFR bearing chromosomes (4/5/8/10) was performed. Topology comparison approach categorized members of 80 families into five distinct co-duplicated groups. Genes belonging to one co-duplicated group are duplicated concurrently, whereas genes of two different co-duplicated groups do not share their duplication history and have not duplicated in congruency. Our findings contradict the 2R model and are indicative of small-scale duplications and rearrangements that cover the entire span of animal's history.

Syndactyly type II or synpolydactyly (SPD) is the second most frequent syndactyly type and is inherited in an autosomal dominant fashion. The cardinal features of this malformation are the cutaneous or bony fusion of third and fourth fingers, and fourth and fifth toes associated with additional digital elements within the web. It shows incomplete penetrance and high inter- and intrafamilial phenotypic variability. Two loci are known for SPD (MIM 186000, MIM 608180) associated with mutations in HOXD13 and FBLN1, respectively. Here, we report further genetic heterogeneity for SDP. Employing a whole genomic screen, we demonstrate, in a large Pakistani kindred, that the classical phenotype of SPD maps on a new locus at chromosome 14q11.2-q12. The highest LOD score (Z(max) = 4.06) was obtained with microsatellite marker D14S264, and the multipoint LOD score reached a maximum of 5.01. Haplotype analysis revealed that the disease interval is flanked by microsatellite markers D14S283 and D14S1060, encompassing a physical distance of 10.72 Mb. We propose to allocate to this locus the symbol SPD3 (synpolydactyly 3), and to name the loci associated with HOXD13 or FBLN1 mutations SPD1 and SPD2, respectively.

This article describes the identification of a novel locus (DFNB39) responsible for an autosomal recessive form of hearing loss segregating in a Pakistani consanguineous family. The hearing impaired members of this family present with profound prelingual sensorineural hearing impairment and use sign language for communications. Linkage was established to microsatellite markers located on chromosome 7q with a maximum multipoint lod score of 3.8. The region of homozygosity spans a 19 cM region that is bounded by markers D7S3046 and D7S644.

GLI2, a zinc finger transcription factor, mediates Sonic hedgehog signaling, a critical pathway in vertebrate embryogenesis. GLI2 has been implicated in diverse set of embryonic developmental processes, including patterning of central nervous system and limbs. In humans, mutations in GLI2 are associated with several developmental defects, including holoprosencephaly and polydactyly.Here, we demonstrate in transient transgenic zebrafish assays, the potential of a subset of tetrapod-teleost conserved non-coding elements (CNEs) residing within human GLI2 intronic intervals to induce reporter gene expression at known regions of endogenous GLI2 transcription. The regulatory activities of these elements are observed in several embryonic domains, including neural tube and pectoral fin. Moreover, our data reveal an overlapping expression profile of duplicated copies of an enhancer during zebrafish evolution.Our data suggest that during vertebrate history GLI2 acquired a high level of complexity in the genetic mechanisms regulating its expression during spatiotemporal patterning of the central nervous system (CNS) and limbs.

Background Understanding the genetic mechanisms underlying the organismal complexity and origin of novelties during vertebrate history is one of the central goals of evolutionary biology. Ohno (1970) was the first to postulate that whole genome duplications (WGD) have played a vital role in the evolution of new gene functions: permitting an increase in morphological, physiological and anatomical complexity during early vertebrate history. Results Here, we analyze the evolutionary history of human FGFR-bearing paralogon (human autosome 4/5/8/10) by the phylogenetic analysis of multigene families with triplicate and quadruplicate distribution on these chromosomes. Our results categorized the histories of 21 families into discrete co-duplicated groups. Genes of a particular co-duplicated group exhibit identical evolutionary history and have duplicated in concert with each other, whereas genes belonging to different groups have dissimilar histories and have not duplicated concurrently. Conclusion Taken together with our previously published data, we submit that there is sufficient empirical evidence to disprove the 1R/2R hypothesis and to support the general prediction that vertebrate genome evolved by relatively small-scale, regional duplication events that spread across the history of life.

Background Human genome is enriched with thousands of conserved non-coding elements (CNEs). Recently, a medium throughput strategy was employed to analyze the ability of human CNEs to drive tissue specific expression during mouse embryogenesis. These data led to the establishment of publicly available genome wide catalog of functionally defined human enhancers. Scattering of enhancers over larger regions in vertebrate genomes seriously impede attempts to pinpoint their precise target genes. Such associations are prerequisite to explore the significance of this in vivo characterized catalog of human enhancers in development, disease and evolution. Results This study is an attempt to systematically identify the target gene-bodies for functionally defined human CNE-enhancers. For the purpose we adopted the orthology/paralogy mapping approach and compared the CNE induced reporter expression with reported endogenous expression pattern of neighboring genes. This procedure pinpointed specific target gene-bodies for the total of 192 human CNE-enhancers. This enables us to gauge the maximum genomic search space for enhancer hunting: 4 Mb of genomic sequence around the gene of interest (2 Mb on either side). Furthermore, we used human-rodent comparison for a set of 159 orthologous enhancer pairs to infer that the central nervous system (CNS) specific gene expression is closely associated with the cooperative interaction among at least eight distinct transcription factors: SOX5, HFH, SOX17, HNF3β, c-FOS, Tal1beta-E47S, MEF and FREAC. Conclusions In conclusion, the systematic wiring of cis-acting sites and their target gene bodies is an important step to unravel the role of in vivo characterized catalog of human enhancers in development, physiology and medicine.

Background Susumu Ohno’s idea that modern vertebrates are degenerate polyploids (concept referred as 2R hypothesis) has been the subject of intense debate for past four decades. It was proposed that intra-genomic synteny regions (paralogons) in human genome are remains of ancient polyploidization events that occurred early in the vertebrate history. The quadruplicated paralogon centered on human HOX clusters is taken as evidence that human HOX-bearing chromosomes were structured by two rounds of whole genome duplication (WGD) events. Results Evolutionary history of human HOX-bearing chromosomes (chromosomes 2/7/12/17) was evaluated by the phylogenetic analysis of multigene families with triplicated or quadruplicated distribution on these chromosomes. Topology comparison approach categorized the members of 44 families into four distinct co-duplicated groups. Distinct gene families belonging to a particular co-duplicated group, exhibit similar evolutionary history and hence have duplicated simultaneously, whereas genes of two distinct co-duplicated groups do not share their evolutionary history and have not duplicated in concert with each other. Conclusion The recovery of co-duplicated groups suggests that “ancient segmental duplications and rearrangements” is the most rational model of evolutionary events that have generated the triplicated and quadruplicated paralogy regions seen on the human HOX-bearing chromosomes.

Hair is a unique mammalian trait that is absent in all other animal forms. Hairlessness is rare in mammals and humans are exceptional among primates in lacking dense layer of hair covering. HR was the first gene identified to be implicated in hair-cycle regulation. Point mutations in HR lead to congenital human hair loss, which results in the complete loss of body and scalp hairs. HR functions are indispensable for initiation of postnatal hair follicular cycling. This study investigates the phylogenetic history and analyzes the protein evolutionary rate to provide useful insight into the molecular evolution of HR. The data demonstrates an acceleration of HR sequence evolution in human branch and suggests that the ability of HR protein to mediate postnatal hair-cycling has been altered in the course of human evolution. In particular those residues were pinpointed which should be regarded as target of positive Darwinian selection during human evolution.

The new body of evidence from fossils and comparative-developmental analysis of subset of appendicular patterning genes has revealed that limb elements seen in tetrapods are assembled in fish fin over evolutionary time. However, despite of deep homology in basic structure and underlying developmental system, there remains a large morphological gap between distal elements of tetrapod limb and distal fin skeleton of tetrapodomorph fish. Understanding the genetic basis of major transformations in distal-limb morphology is the next challenge for evolutionary developmental biologists. Here by integrating data from fossils, comparative-developmental and genetic studies, models are proposed describing the evolution of cis-regulatory elements as a basis for diversification of appendicular architecture. Instead of emphasizing the subset of developmental genes, for instance Hoxd genes, the focus here is on the significance of elucidating cis-regulatory elements for multiple other key molecular players of limb/fin development and genetic/molecular interactions among them, for a better understanding of the developmental and genetic basis of limb evolution.

Background The zinc-finger transcription factor GLI3 is an important mediator of Sonic hedgehog signaling and crucial for patterning of many aspects of the vertebrate body plan. In vertebrates, the mechanism of SHH signal transduction and its action on target genes by means of activating or repressing forms of GLI3 have been studied most extensively during limb development and the specification of the central nervous system. From these studies it has emerged, that Gli3 expression must be subject to a tight spatiotemporal regulation. However, the genetic mechanisms and the cis-acting elements controlling the expression of Gli3 remained largely unknown. Results Here, we demonstrate in chicken and mouse transgenic embryos that human GLI3-intronic conserved non-coding sequence elements (CNEs) autonomously control individual aspects of Gli3 expression. Their combined action shows many aspects of a Gli3-specific pattern of transcriptional activity. In the mouse limb bud, different CNEs enhance Gli3-specific expression in evolutionary ancient stylopod and zeugopod versus modern skeletal structures of the autopod. Limb bud specificity is also found in chicken but had not been detected in zebrafish embryos. Three of these elements govern central nervous system specific gene expression during mouse embryogenesis, each targeting a subset of endogenous Gli3 transcription sites. Even though fish, birds, and mammals share an ancient repertoire of gene regulatory elements within Gli3, the functions of individual enhancers from this catalog have diverged significantly. During evolution, ancient broad-range regulatory elements within Gli3 attained higher specificity, critical for patterning of more specialized structures, by abolishing the potential for redundant expression control. Conclusion These results not only demonstrate the high level of complexity in the genetic mechanisms controlling Gli3 expression, but also reveal the evolutionary significance of cis-acting regulatory networks of early developmental regulators in vertebrates.

Background Vertebrate genomes contain extensive intra-genomic conserved synteny, which is the presence of similar set of genes on two or more chromosomes (paralogons). The existence of these paralogons has led to the proposal that vertebrate genome was structured by one or more rounds of ancient whole genome duplications (2R hypothesis). Results The 2R hypothesis was tested by phylogenetic analysis of gene families residing on human HOX-bearing chromosomes (HOX-cluster paralogons). These results revealed that, based on their duplication history, 23 gene families with representation on three or four of the human HOX-bearing chromosomes can be partitioned into four discrete co-duplicated groups. The distinct genes within each co-duplicated group share the same evolutionary history and are duplicated in concert with each other, while the constituent genes of two different co-duplicated groups do not share their evolutionary history and are not duplicated simultaneously. These co-duplicated groups are large constituting members from 3 to 8 gene families and suggest that human HOX-cluster paralogons were shaped by ancient segmental duplications (SDs) and rearrangement events that occurred at least as early as before the divergence of bony fishes and tetrapods. Conclusions Based on the recovery of ancient SD events in this analysis and given the widespread evidence in favor of the fact that recent SD events played a pivotal role in changing genome architecture of primates and other recently diverged animals, it is concluded that a more realistic model of ancient vertebrate genome evolutionary history can be deduced by tracing the evolutionary trajectory of the genomes of recently diverged vertebrate species.

Background The human chromosomes 2q, 7, 12q and 17q show extensive intra-genomic homology, containing duplicate, triplicate and quadruplicate paralogous regions centered on the HOX gene clusters. The fact that two or more representatives of different gene families are linked with HOX clusters is taken as evidence that these paralogous gene sets might have arisen from a single chromosomal segment through block or whole chromosome duplication events. This would imply that the constituent genes including the HOX clusters reflect the architecture of a single ancestral block (before vertebrate origin) where all of these genes were linked in a single copy. Results In the present study we have employed the currently available set of protein data for a wide variety of vertebrate and invertebrate genomes to analyze the phylogenetic history of 11 multigene families with three or more of their representatives linked to human HOX clusters. A topology comparison approach revealed four discrete co-duplicated groups: group 1 involves the genes from GLI, HH, INHB, IGFBP (cluster-1), and SLC4A families; group 2 involves ERBB, ZNFN1A, and IGFBP (cluster-2) gene families; group 3 involves the HOX clusters and the SP gene family; group 4 involves the integrin beta chain and myosine light chain families. The distinct genes within each co-duplicated group share the same evolutionary history and are duplicated in concert with each other, while the constituent genes of two different co-duplicated groups may not share their evolutionary history and may not have duplicated simultaneously. Conclusion We conclude that co-duplicated groups may themselves be remnants of ancient small-scale duplications (involving chromosomal segments or gene-clusters) which occurred at different time points during chordate evolution. Whereas the recent combination of genes from distinct co-duplicated groups on different chromosomal regions (human chromosomes 2q, 7, 12q, and 17q) is probably the outcome of subsequent rearrangement of genomic segments, including syntenic groups of genes.

The zinc-finger transcription factor GLI3 is a key regulator of development, acting as a primary transducer of Sonic hedgehog (SHH) signaling in a combinatorial context dependent fashion controlling multiple patterning steps in different tissues/organs. A tight temporal and spatial control of gene expression is indispensable, however, cis-acting sequence elements regulating GLI3 expression have not yet been reported. We show that 11 ancient genomic DNA signatures, conserved from the pufferfish Takifugu (Fugu) rubripes to man, are distributed throughout the introns of human GLI3. They map within larger conserved non-coding elements (CNEs) that are found in the tetrapod lineage. Full length CNEs transiently transfected into human cell cultures acted as cell type specific enhancers of gene transcription. The regulatory potential of these elements is conserved and was exploited to direct tissue specific expression of a reporter gene in zebrafish embryos. Assays of deletion constructs revealed that the human-Fugu conserved sequences within the GLI3 intronic CNEs were essential but not sufficient for full-scale transcriptional activation. The enhancer activity of the CNEs is determined by a combinatorial effect of a core sequence conserved between human and teleosts (Fugu) and flanking tetrapod-specific sequences, suggesting that successive clustering of sequences with regulatory potential around an ancient, highly conserved nucleus might be a possible mechanism for the evolution of cis-acting regulatory elements.