Search

Enhancers are often studied as noncoding regulatory elements that modulate the precise spatiotemporal expression of genes in a highly tissue-specific manner. This paradigm has been challenged by recent evidence of individual enhancers acting in multiple tissues or developmental contexts. However, the frequency of these enhancers with high degrees of “pleiotropy” out of all putative enhancers is not well understood. Consequently, it is unclear how the variation of enhancer pleiotropy corresponds to the variation in expression breadth of target genes. Here, we use multi-tissue chromatin maps from diverse human tissues to investigate the enhancer–gene interaction architecture while accounting for 1) the distribution of enhancer pleiotropy, 2) the variations of regulatory links from enhancers to target genes, and 3) the expression breadth of target genes. We show that most enhancers are tissue-specific and that highly pleiotropy enhancers account for

Data for the article 'Transmembrane helices are an over-presented and evolutionarily conserved source of major histocompatibility complex class I and II epitopes' by Rich��l J.C. Bilderbeek, Maksim V. Baranov, Geert van den Bogaart and Frans Bianchi. Link to article: [https://github.com/richelbilderbeek/bbbq_article](https://github.com/richelbilderbeek/bbbq_article), Frans Bianchi is funded by a Veni grant from the Netherlands Organization for Scientific Research (016.Veni.192.026) and an Off-Road Grant from the Dutch Medical Science Foundation (ZonMW 04510011910005). Geert van den Bogaart is funded by a Young Investigator Grant from the Human Frontier Science Program (HFSP; RGY0080/2018), and a Vidi grant from the Netherlands Organization for Scientific Research (NWO-ALW VIDI 864.14.001). Geert van den Bogaart has received funding from the European Research Council (ERC) under the European Union's Horizon 2020 research and innovation programme (grant agreement No. 862137.

In the past few years, there has been an increasing neuroscientific interest in understanding the function of mammalian chromodomains helicase DNA-binding (CHD) proteins due to their association with severe developmental syndromes. Mammalian CHDs include nine members (CHD1 to CHD9), grouped into subfamilies according to the presence of specific functional domains, generally highly conserved in evolutionary terms. Mutations affecting these domains hold great potential to disrupt protein function, leading to meaningful pathogenic scenarios, such as embryonic defects incompatible with life. Here, we analysed the evolution of CHD proteins by performing a comparative study of the functional domains of CHD proteins between orthologous and paralogous protein sequences. Our findings show that the highest degree of inter-species conservation was observed at Group II (CHD3, CHD4, and CHD5) and that most of the pathological variations documented in humans involve amino acid residues that are conserved not only between species but also between paralogs. The parallel analysis of both orthologous and paralogous proteins, in cases where gene duplications have occurred, provided extra information showing patterns of flexibility as well as interchangeability between amino acid positions. This added complexity needs to be considered when the impact of novel mutations is assessed in terms of evolutionary conservation.

It is now more than 20 years since the FRA16D common chromosomal fragile site was characterised and the WWOX gene spanning this site was identified. In this time, much information has been discovered about its contribution to disease; however, the normal biological role of WWOX is not yet clear. Experiments leading to the identification of the WWOX gene are recounted, revealing enigmatic relationships between the fragile site, its gene and the encoded protein. We also highlight research mainly using the genetically tractable model organism Drosophila melanogaster that has shed light on the integral role of WWOX in metabolism. In addition to this role, there are some particularly outstanding questions that remain regarding WWOX, its gene and its chromosomal location. This review, therefore, also aims to highlight two unanswered questions. Firstly, what is the biological relationship between the WWOX gene and the FRA16D common chromosomal fragile site that is located within one of its very large introns? Secondly, what is the actual substrate and product of the WWOX enzyme activity? It is likely that understanding the normal role of WWOX and its relationship to chromosomal fragility are necessary in order to understand how the perturbation of these normal roles results in disease.

Topologically associating domains (TADs) are fundamental units of three-dimensional (3D) nuclear organization. The regions bordering TADs-TAD boundaries-contribute to the regulation of gene expression by restricting interactions of cis-regulatory sequences to their target genes. TAD and TAD-boundary disruption have been implicated in rare-disease pathogenesis; however, we have a limited framework for integrating TADs and their variation across cell types into the interpretation of common-trait-associated variants. Here, we investigate an attribute of 3D genome architecture-the stability of TAD boundaries across cell types-and demonstrate its relevance to understanding how genetic variation in TADs contributes to complex disease. By synthesizing TAD maps across 37 diverse cell types with 41 genome-wide association studies (GWASs), we investigate the differences in disease association and evolutionary pressure on variation in TADs versus TAD boundaries. We demonstrate that genetic variation in TAD boundaries contributes more to complex-trait heritability, especially for immunologic, hematologic, and metabolic traits. We also show that TAD boundaries are more evolutionarily constrained than TADs. Next, stratifying boundaries by their stability across cell types, we find substantial variation. Compared to boundaries unique to a specific cell type, boundaries stable across cell types are further enriched for complex-trait heritability, evolutionary constraint, CTCF binding, and housekeeping genes. Thus, considering TAD boundary stability across cell types provides valuable context for understanding the genome's functional landscape and enabling variant interpretation that takes 3D structure into account.

Understanding how protein function has evolved and diversified is of great importance for human genetics and medicine. Here, we tackle the problem of describing the whole transcript variability observed in several species by generalising the definition of splicing graph. We provide a practical solution to building parsimonious evolutionary splicing graphs where each node is a minimal transcript building block defined across species. We show a clear link between the functional relevance, tissue-regulation and conservation of AS events on a set of 50 genes. By scaling up to the whole human protein-coding genome, we identify a few thousands of genes where alternative splicing modulates the number and composition of pseudo-repeats. We have implemented our approach in ThorAxe, an efficient, versatile, and robust computational tool freely available at https://github.com/PhyloSofS-Team/thoraxe. The results are accessible and can be browsed interactively at http://www.lcqb.upmc.fr/ThorAxe.

The opioid receptor (OPR) family comprises the mu-, delta-, and kappa-opioid, and nociceptin receptors that belong to the superfamily of 7-transmembrane spanning G protein-coupled receptors (GPCRs). The mu-opioid receptor is the main target for clinically used opioid analgesics, and its biology has been extensively studied. The N-terminally truncated 6TM receptors isoform produced through alternative splicing of the OPRM1 gene displays unique signaling and analgesic properties, but it is unclear if other OPRs have the same ability. In this study, we have built a comprehensive map of alternative splicing events that produce 6TM receptor variants in all the OPRs and demonstrated their evolutionary conservation. We then obtained evidence for their translation through ribosomal footprint analysis. We discovered that N-terminally truncated 6TM GPCRs are rare in the human genome and OPRs are overrepresented in this group. Finally, we also observed a significant enrichment of 6TM GPCR genes among genes associated with pain, psychiatric disorders, and addiction. Understanding the biology of 6TM receptors and leveraging this knowledge for drug development should pave the way for novel therapies. Electronic supplementary material The online version of this article (10.1007/s10571-020-00971-7) contains supplementary material, which is available to authorized users.

Background Structural homology modeling supported by bioinformatics analysis plays a key role in uncovering new molecular interactions within gene regulatory networks. Here, we have applied this powerful approach to analyze the molecular interactions orchestrating death receptor signaling networks. In particular, we focused on the molecular mechanisms of CD95-mediated NF-κB activation and the role of c-FLIP/NEMO interaction in the induction of this pathway. Results To this end, we have created the homology model of the c-FLIP/NEMO complex using the reported structure of the v-FLIP/NEMO complex, and rationally designed peptides targeting this complex. The designed peptides were based on the NEMO structure. Strikingly, the experimental in vitro validation demonstrated that the best inhibitory effects on CD95-mediated NF-κB activation are exhibited by the NEMO-derived peptides with the substitution D242Y of NEMO. Furthermore, we have assumed that the c-FLIP/NEMO complex is recruited to the DED filaments formed upon CD95 activation and validated this assumption in silico. Further insight into the function of c-FLIP/NEMO complex was provided by the analysis of evolutionary conservation of interacting regions which demonstrated that this interaction is common in distinct mammalian species. Conclusions Taken together, using a combination of bioinformatics and experimental approaches we obtained new insights into CD95-mediated NF-κB activation, providing manifold possibilities for targeting the death receptor network.

Significance Although gamma delta (γδ) T cells compose an evolutionarily conserved third lineage of diversified lymphocytes, alongside αβ T cells and B cells, they can seem overtly different across species and tissues. Thus, human blood γδ cells show butyrophilin (BTN)3A1-dependent responses to metabolites (“phosphoantigens”) not seen by rodent γδ cells, whereas some rodent, γδ-rich compartments, notably in the skin, lack obvious human counterparts. Recently, however, mouse and human intraepithelial gut γδ cells were found to be regulated by pairings of BTN-like genes. This study now shows that BTN3A1 also functions as a pairing, with its subcellular trafficking and optimal activity both regulated by BTN3A2. Hence, seemingly diverse γδ cell biologies across species and tissues are underpinned by conserved mechanisms., The long-held view that gamma delta (γδ) T cells in mice and humans are fundamentally dissimilar, as are γδ cells in blood and peripheral tissues, has been challenged by emerging evidence of the cells’ regulation by butyrophilin (BTN) and butyrophilin-like (BTNL) molecules. Thus, murine Btnl1 and the related gene, Skint1, mediate T cell receptor (TCR)-dependent selection of murine intraepithelial γδ T cell repertoires in gut and skin, respectively; BTNL3 and BTNL8 are TCR-dependent regulators of human gut γδ cells; and BTN3A1 is essential for TCR-dependent activation of human peripheral blood Vγ9Vδ2+ T cells. However, some observations concerning BTN/Btnl molecules continue to question the extent of mechanistic conservation. In particular, murine and human gut γδ cell regulation depends on pairings of Btnl1 and Btnl6 and BTNL3 and BTNL8, respectively, whereas blood γδ cells are reported to be regulated by BTN3A1 independent of other BTNs. Addressing this paradox, we show that BTN3A2 regulates the subcellular localization of BTN3A1, including functionally important associations with the endoplasmic reticulum (ER), and is specifically required for optimal BTN3A1-mediated activation of Vγ9Vδ2+ T cells. Evidence that BTNL3/BTNL8 and Btnl1/Btnl6 likewise associate with the ER reinforces the prospect of broadly conserved mechanisms underpinning the selection and activation of γδ cells in mice and humans, and in blood and extralymphoid sites.

Studies of regulatory activity and gene expression have revealed an intriguing dichotomy: There is substantial turnover in the regulatory activity of orthologous sequences between species; however, the expression level of orthologous genes is largely conserved. Understanding how distal regulatory elements, for example, enhancers, evolve and function is critical, as alterations in gene expression levels can drive the development of both complex disease and functional divergence between species. In this study, we investigated determinants of the conservation of regulatory enhancer activity for orthologous sequences across mammalian evolution. Using liver enhancers identified from genome-wide histone modification profiles in ten diverse mammalian species, we compared orthologous sequences that exhibited regulatory activity in all species (conserved-activity enhancers) to shared sequences active only in a single species (species-specific-activity enhancers). Conserved-activity enhancers have greater regulatory potential than species-specific-activity enhancers, as quantified by both the density and diversity of transcription factor binding motifs. Consistent with their greater regulatory potential, conserved-activity enhancers have greater regulatory activity in humans than species-specific-activity enhancers: They are active across more cellular contexts, and they regulate more genes than species-specific-activity enhancers. Furthermore, the genes regulated by conserved-activity enhancers are expressed in more tissues and are less tolerant of loss-of-function mutations than those targeted by species-specific-activity enhancers. These consistent results across various stages of gene regulation demonstrate that conserved-activity enhancers are more pleiotropic than their species-specific-activity counterparts. This suggests that pleiotropy is associated with the conservation of regulatory across mammalian evolution.

Side effects from targeted drugs remain a serious concern. One reason is the nonselective binding of a drug to unintended proteins such as its paralogs, which are highly homologous in sequences and have similar structures and drug-binding pockets. To identify targetable differences between paralogs, we analyzed two types (type-I and type-II) of functional divergence between two paralogs in the known target protein receptor family G-protein coupled receptors (GPCRs) at the amino acid level. Paralogous protein receptors in glucagon-like subfamily, glucagon receptor (GCGR) and glucagon-like peptide-1 receptor (GLP-1R), exhibit divergence in ligands and are clinically validated drug targets for type 2 diabetes. Our data showed that type-II amino acids were significantly enriched in the binding sites of antagonist MK-0893 to GCGR, which had a radical shift in physicochemical properties between GCGR and GLP-1R. We also examined the role of type-I amino acids between GCGR and GLP-1R. The divergent features between GCGR and GLP-1R paralogs may be helpful in their discrimination, thus enabling the identification of binding sites to reduce undesirable side effects and increase the target specificity of drugs.

Contains fulltext : 182589.pdf (Publisher’s version ) (Open Access) Whole exomes of patients with a genetic disorder are nowadays routinely sequenced but interpretation of the identified genetic variants remains a major challenge. The increased availability of population-based human genetic variation has given rise to measures of genetic tolerance that have been used, for example, to predict disease-causing genes in neurodevelopmental disorders. Here, we investigated whether combining variant information from homologous protein domains can improve variant interpretation. For this purpose, we developed a framework that maps population variation and known pathogenic mutations onto 2,750 "meta-domains." These meta-domains consist of 30,853 homologous Pfam protein domain instances that cover 36% of all human protein coding sequences. We find that genetic tolerance is consistent across protein domain homologues, and that patterns of genetic tolerance faithfully mimic patterns of evolutionary conservation. Furthermore, for a significant fraction (68%) of the meta-domains high-frequency population variation re-occurs at the same positions across domain homologues more often than expected. In addition, we observe that the presence of pathogenic missense variants at an aligned homologous domain position is often paired with the absence of population variation and vice versa. The use of these meta-domains can improve the interpretation of genetic variation.

Despite the tremendous growth of the DNA sequencing data in the last decade, our understanding of the human genome is still in its infancy. To understand the implications of genetic variants in the light of population genetics and molecular evolution, we developed a database, PGG.SNV (https://www.pggsnv.org), which gives much higher weight to previously under-investigated indigenous populations in Asia. PGG.SNV archives 265 million SNVs across 220,147 present-day genomes and 1018 ancient genomes, including 1009 newly sequenced genomes, representing 977 global populations. Moreover, estimation of population genetic diversity and evolutionary parameters is available in PGG.SNV, a unique feature compared with other databases.

The growing body of experimental and computational data describing how proteins interact with each other has emphasized the multiplicity of protein interactions and the complexity underlying protein surface usage and deformability. In this work, we propose new concepts and methods toward deciphering such complexity. We introduce the notion of interacting region to account for the multiple usage of a protein's surface residues by several partners and for the variability of protein interfaces coming from molecular flexibility. We predict interacting patches by crossing evolutionary, physicochemical and geometrical properties of the protein surface with information coming from complete cross‐docking (CC‐D) simulations. We show that our predictions match well interacting regions and that the different sources of information are complementary. We further propose an indicator of whether a protein has a few or many partners. Our prediction strategies are implemented in the dynJET2 algorithm and assessed on a new dataset of 262 protein on which we performed CC‐D. The code and the data are available at: http://www.lcqb.upmc.fr/dynJET2/.

Inflammasomes are cytosolic multiprotein complexes that regulate inflammatory responses to danger stimuli and infection, and their dysregulation is associated with an increasing number of autoinflammatory diseases. In recent years, zebrafish models of human pathologies to study inflammasome function in vivo have started to emerge. Here, we discuss inflammasome research in zebrafish in light of current knowledge about mammalian inflammasomes. We summarize the evolutionary conservation of inflammasome components between zebrafish and mammals, highlighting the similarities and possible divergence in functions of these components.We present new insights into the evolution of the caspase-1 family in the teleost lineage, and how its evolutionaryorigin may help contextualize its functions. We also review existing infectious and non-infectious models in zebrafish in which inflammasomes have been directly implicated. Finally, we discuss the advantages of zebrafish larvae for intravital imaging of inflammasome activation and summarize available tools that will help to advance inflammasome research.

Summary FOXO transcription factors (FOXOs) are central regulators of lifespan across species, yet they also have cell‐specific functions, including adult stem cell homeostasis and immune function. Direct targets of FOXOs have been identified genome‐wide in several species and cell types. However, whether FOXO targets are specific to cell types and species or conserved across cell types and throughout evolution remains uncharacterized. Here, we perform a meta‐analysis of direct FOXO targets across tissues and organisms, using data from mammals as well as Caenorhabditis elegans and Drosophila. We show that FOXOs bind cell type‐specific targets, which have functions related to that particular cell. Interestingly, FOXOs also share targets across different tissues in mammals, and the function and even the identity of these shared mammalian targets are conserved in invertebrates. Evolutionarily conserved targets show enrichment for growth factor signaling, metabolism, stress resistance, and proteostasis, suggesting an ancestral, conserved role in the regulation of these processes. We also identify candidate cofactors at conserved FOXO targets that change in expression with age, including CREB and ETS family factors. This meta‐analysis provides insight into the evolution of the FOXO network and highlights downstream genes and cofactors that may be particularly important for FOXO's conserved function in adult homeostasis and longevity.

Although evidence indicates that drug target genes share some common evolutionary features, there have been few studies analyzing evolutionary features of drug targets from an overall level. Therefore, we conducted an analysis which aimed to investigate the evolutionary characteristics of drug target genes. We compared the evolutionary conservation between human drug target genes and non-target genes by combining both the evolutionary features and network topological properties in human protein-protein interaction network. The evolution rate, conservation score and the percentage of orthologous genes of 21 species were included in our study. Meanwhile, four topological features including the average shortest path length, betweenness centrality, clustering coefficient and degree were considered for comparison analysis. Then we got four results as following: compared with non-drug target genes, 1) drug target genes had lower evolutionary rates; 2) drug target genes had higher conservation scores; 3) drug target genes had higher percentages of orthologous genes and 4) drug target genes had a tighter network structure including higher degrees, betweenness centrality, clustering coefficients and lower average shortest path lengths. These results demonstrate that drug target genes are more evolutionarily conserved than non-drug target genes. We hope that our study will provide valuable information for other researchers who are interested in evolutionary conservation of drug targets.

Classification of rare BRCA1 missense variants presents a major challenge for the counseling and treatment of patients. Variant classification can be complicated by conflicting lines of evidence. BRCA1 c.5309G>T p.(Gly1770Val) has been shown to abrogate BRCA1 protein homologous DNA repair; however, multiple sequence alignment demonstrates a lack of sequence conservation at this position, suggesting that glycine at position 1770 may not be essential for cellular maintenance in humans. We analyzed clinical information to resolve the classification of BRCA1 c.5309G>T p.(Gly1770Val). We performed multifactorial likelihood analysis combining segregation data for 14 informative families, and breast tumor histopathological data for 17 variant carriers, ascertained through the ENIGMA consortium. Bayes segregation analysis gave a likelihood ratio of 101:1 in favor of pathogenicity. The vast majority of breast tumors showed features indicative of pathogenic variant carrier status, resulting in a likelihood ratio of 15800794:1 towards pathogenicity. Despite a low prior probability of pathogenicity (0.03) based on bioinformatic prediction, multifactorial likelihood analysis including segregation and histopathology analysis gave a posterior probability of > 0.99 and final classification of Pathogenic. We provide evidence that BRCA1 c.5309G>T p.(Gly1770Val), previously described as a Moroccan founder variant, should be treated as a disease-causing variant despite a lack of evolutionary conservation at this amino acid position. Additionally, we stress that bioinformatic information should be used in combination with other data, either direct clinical evidence or some form of clinical calibration, to arrive at a final clinical classification.

International audience; Human ribonuclease 4 (RNase 4) is the most evolutionarily conserved member of the 8 canonical human pancreatic-like RNases, showing more than 90 % identity with bovine and porcine homologues. The enzyme displays ribonucleolytic activity with a strong preference for uracil-containing RNA substrates, a feature only shared with human eosinophil derived-neurotoxin (EDN, or RNase 2) and eosinophil cationic protein (ECP, or RNase 3). It is also the shortest member of the human family, with a significantly truncated C-terminal tail. Its unique active-site pocket and high degree of conservation among vertebrates suggest that the enzyme plays a crucial biological function. Here, we report on the (1)H, (13)C and (15)N backbone resonance assignments of RNase 4, providing means to characterize its molecular function at the atomic level by NMR.

Assigning pathogenicity to mt-tRNA variants requires multiple strands of evidence. Evolutionary conservation is often considered mandatory, but lack of a standard panel of organisms to assess conservation complicates comparison between reports and undermines the value of conservation-based evidence. We demonstrate that intra-species MTT sequence variation is sufficiently low for sequence data from a single organism to adequately represent a species. On this basis, we propose a standardised panel of organisms for conservation assessment and describe integration of this conservation panel into a pathogenicity scoring system designed to assess mt-tRNA variation associated with mitochondrial disease., Highlights ► Assigning pathogenicity to mt-tRNA variants requires multiple veins of evidence, including evolutionary conservation; however, there is no “standard sequence panel”. ► Species used for conservation assessment previously, noting alarming variation. ► Whether a panel with a single sequence from each organism can be used to assess conservation, given intra species variation. ► We set out robust criteria for sequence selection and suggest a “standard sequence panel”. ► Aid decisions on when to carry out functional studies upon identification of novel mutations and improve consistency of mt-tRNA mutation reporting.

Recent publications have revealed that the evolution of phosphosites is influenced by the local protein structures and whether the phosphosites have characterized functions or not. With knowledge of the wide functional range of phosphorylation, we attempted to clarify whether the evolutionary conservation of phosphosites is different among distinct functional modules. We grouped the phosphosites in the human genome into the modules according to the functional categories of KEGG (Kyoto Encyclopedia of Genes and Genomes) and investigated their evolutionary conservation in vertebrate genomes from mouse to zebrafish. We have found that the phosphosites in the vertebrate-specific functional modules (VFMs), such as cellular signaling processes and responses to stimuli, are evolutionarily more conserved than those in the basic functional modules (BFMs), such as metabolic and genetic processes. The phosphosites in the VFMs are also significantly more conserved than their flanking regions, whereas those in the BFMs are not. These results hold for both serine/threonine and tyrosine residues, although the fraction of phosphorylated tyrosine residues is increased in the VFMs. Moreover, the difference in the evolutionary conservation of the phosphosites between the VFMs and BFMs could not be explained by the difference in the local protein structures. There is also a higher fraction of phosphosites with known functions in the VFMs than BFMs. Based on these findings, we have concluded that protein phosphorylation may play more dominant roles for the VFMs than BFMs during the vertebrate evolution. As phosphorylation is a quite rapid biological reaction, the VFMs that quickly respond to outer stimuli and inner signals might heavily depend on this regulatory mechanism. Our results imply that phosphorylation may have an essential role in the evolution of vertebrates.

Hemophilia B has been considered as a classical monogenic disorder with complete penetrance. Here, we observed that its allelic spectrum encompasses mutations that show least phenotypic variation as well as those showing phenotypic plasticity, thereby establishing hemophilia B as a quasi-quantitative condition with variable expressivity. Thus, we relocate it from simple monogenic diseases with complete penetrance into a space that marks the conceptual continuum of Mendelian to complex diseases. By computational analysis, we show that mutations showing phenotypic variation were characterized by relatively less conserved mutant sites and more tolerated conservative substitutions. We also demonstrate that consideration of continuous phenotypes renders quantitative rigor for interrogation and hence higher predictive value while analyzing for differential properties than classifications based on clinical end points. Certain mutations have been consistently reported to cause mild phenotype, calling for a check in indiscriminate termination of fetuses following prenatal diagnosis.

It has been postulated that noncoding RNAs (ncRNAs) are involved in the posttranscriptional control of gene expression, and may have contributed to the emergence of the complex attributes observed in mammalians. We show here that the complement of ncRNAs expressed from intronic regions of the human and mouse genomes comprises at least 78,147 and 39,660 transcriptional units, respectively. To identify conserved intronic sequences expressed in both humans and mice, we used custom-designed human cDNA microarrays to separately interrogate RNA from mouse and human liver, kidney, and prostate tissues. An overlapping tissue expression signature was detected for both species, comprising 198 transcripts; among these, 22 RNAs map to intronic regions with evidence of evolutionary conservation in humans and mice. Transcription of selected human–mouse intronic ncRNAs was confirmed using strand-specific RT-PCR. Altogether, these results support an evolutionarily conserved role of intronic ncRNAs in human and mouse, which are likely to be involved in the fine tuning of gene expression regulation in different mammalian tissues.

Identifying the genetic input for fetal growth will help to understand common, serious complications of pregnancy such as fetal growth restriction. Genomic imprinting is an epigenetic process that silences one parental allele, resulting in monoallelic expression. Imprinted genes are important in mammalian fetal growth and development. Evidence has emerged showing that genes that are paternally expressed promote fetal growth, whereas maternally expressed genes suppress growth. We have assessed whether the expression levels of key imprinted genes correlate with fetal growth parameters during pregnancy, either early in gestation, using chorionic villus samples (CVS), or in term placenta. We have found that the expression of paternally expressing insulin-like growth factor 2 (IGF2), its receptor IGF2R, and the IGF2/IGF1R ratio in CVS tissues significantly correlate with crown-rump length and birthweight, whereas term placenta expression shows no correlation. For the maternally expressing pleckstrin homology-like domain family A, member 2 (PHLDA2), there is no correlation early in pregnancy in CVS but a highly significant negative relationship in term placenta. Analysis of the control of imprinted expression of PHLDA2 gave rise to a maternally and compounded grand-maternally controlled genetic effect with a birthweight increase of 93/155 g, respectively, when one copy of the PHLDA2 promoter variant is inherited. Expression of the growth factor receptor-bound protein 10 (GRB10) in term placenta is significantly negatively correlated with head circumference. Analysis of the paternally expressing delta-like 1 homologue (DLK1) shows that the paternal transmission of type 1 diabetes protective G allele of rs941576 single nucleotide polymorphism (SNP) results in significantly reduced birth weight (2132 g). In conclusion, we have found that the expression of key imprinted genes show a strong correlation with fetal growth and that for both genetic and genomics data analyses, it is important not to overlook parent-of-origin effects.

Members of the platelet-derived growth factor/vascular endothelial growth factor (PDGF/VEGF) family have been implicated in a variety of functions in vertebrates, especially angiogenesis. Here we identify and characterize a PDGF/VEGF-like factor (named PVF-1) from the nematode C. elegans. We show that PVF-1 has biochemical properties similar to vertebrate PDGF/VEGF growth factors. More important, PVF-1 binds to the human receptors VEGFR-1 (Flt-1) and VEGFR-2 (KDR) and is able to induce angiogenesis in two model systems derived from vertebrates. Our results highlight the widespread evolutionary conservation of this important class of growth factors and raise the possibility that C. elegans can provide a simple experimental system in which to investigate how these factors function.—Tarsitano, M., De Falco, S., Colonna, V., McGhee, J. D., Persico, M. G. The C. elegans pvf-1 gene encodes a PDGF/VEGF-like factor able to bind mammalian VEGF receptors and to induce angiogenesis.

Background Nucleic acids containing guanine tracts can form quadruplex structures via non-Watson-Crick base pairing. Formation of G-quadruplexes is associated with the regulation of important biological functions such as transcription, genetic instability, DNA repair, DNA replication, epigenetic mechanisms, regulation of translation, and alternative splicing. G-quadruplexes play important roles in human diseases and are being considered as targets for a variety of therapies. Identification of functional G-quadruplexes and the study of their overall distribution in genomes and transcriptomes is an important pursuit. Traditional computational methods map sequence motifs capable of forming G-quadruplexes but have difficulty in distinguishing motifs that occur by chance from ones which fold into G-quadruplexes. Results We present Quadruplex forming ‘G’-rich sequences (QGRS)-Conserve, a computational method for calculating motif conservation across exomes and supports filtering to provide researchers with more precise methods of studying G-quadruplex distribution patterns. Our method quantitatively evaluates conservation between quadruplexes found in homologous nucleotide sequences based on several motif structural characteristics. QGRS-Conserve also efficiently manages overlapping G-quadruplex sequences such that the resulting datasets can be analyzed effectively. Conclusions We have applied QGRS-Conserve to identify a large number of G-quadruplex motifs in the human exome conserved across several mammalian and non-mammalian species. We have successfully identified multiple homologs of many previously published G-quadruplexes that play post-transcriptional regulatory roles in human genes. Preliminary large-scale analysis identified many homologous G-quadruplexes in the 5′- and 3′-untranslated regions of mammalian species. An expectedly smaller set of G-quadruplex motifs was found to be conserved across larger phylogenetic distances. QGRS-Conserve provides means to build datasets that can be filtered and categorized in a variety of biological dimensions for more targeted studies in order to better understand the roles that G-quadruplexes play.

Background The existence of introns in eukaryotic genes is believed to provide an evolutionary advantage by increasing protein diversity through exon shuffling and alternative splicing. However, this eukaryotic feature is associated with the necessity of exclusion of intronic sequences, which requires considerable energy expenditure and can lead to splicing errors. The relationship between intronic burden and evolution is poorly understood. The goal of this study was to analyze the relationship between the intronic burden and the level of evolutionary conservation of the gene. Results We found a positive correlation between the level of evolutionary conservation of a gene and its intronic burden. The level of evolutionary conservation was estimated using the conservation index (CI). The CI value was determined on the basis of the most distant ortholog of the human protein sequence and ranged from 0 (the gene was unique to the human genome) to 9 (an ortholog of the human gene was detected in plants). In multivariable model, both the number of introns and total intron size remained significant predictors of CI. We also found that the number of alternative splice variants was positively correlated with CI. The expression level of a gene was negatively correlated with the number of introns and total size of intronic region. Genes with a greater intronic burden had lower density of missense and nonsense mutations in the coding regions of the gene, which suggests that they are under a stronger pressure from purifying selection. Conclusions We identified a positive association between intronic burden and CI. One of the possible explanations of this is the idea of a cost-benefits balance. Evolutionarily conserved (functionally important) genes can “afford” the negative consequences of maintaining multiple introns because these consequences are outweighed by the benefit of maintaining the gene. Evolutionarily conserved and functionally important genes may use introns to create novel splice variants to tune the gene function to developmental stage and tissue type.

Experiments on seven vertebrates suggest that identifying the locations of islands of non-methylated DNA provides more insights into evolutionarily-conserved epigenetic regulatory elements than studies of CpG islands.

Two-thirds of gene promoters in mammals are associated with regions of non-methylated DNA, called CpG islands (CGIs), which counteract the repressive effects of DNA methylation on chromatin. In cold-blooded vertebrates, computational CGI predictions often reside away from gene promoters, suggesting a major divergence in gene promoter architecture across vertebrates. By experimentally identifying non-methylated DNA in the genomes of seven diverse vertebrates, we instead reveal that non-methylated islands (NMIs) of DNA are a central feature of vertebrate gene promoters. Furthermore, NMIs are present at orthologous genes across vast evolutionary distances, revealing a surprising level of conservation in this epigenetic feature. By profiling NMIs in different tissues and developmental stages we uncover a unifying set of features that are central to the function of NMIs in vertebrates. Together these findings demonstrate an ancient logic for NMI usage at gene promoters and reveal an unprecedented level of epigenetic conservation across vertebrate evolution. DOI: http://dx.doi.org/10.7554/eLife.00348.001, eLife digest DNA methylation—the addition of a methyl group to cytosine, one of the four bases found in DNA—is a central process in genetics. By preventing genes from being expressed as proteins, DNA methylation is one of a number of epigenetic mechanisms that can determine which proteins are made in different cell types without changing the underlying DNA sequence. In warm-blooded vertebrates such as mammals most of the genome is methylated, however short regions of non-methylated DNA are known to be associated with gene promoters (regions of DNA that act as binding sites for the enzymes and transcription factors that transcribe the DNA in the gene into RNA). Much of our current understanding of the role of these islands of non-methylated DNA is based on computational predictions rather than experimental data. In cold-blooded vertebrates, for example, computer models often predict that non-methylated islands are not associated with gene promoters, which potentially suggests an evolutionary divergence in the role of methylation amongst vertebrates. However, this idea has not been confirmed by experimental data. Long et al. have performed experiments to compare the location of non-methylated islands in seven different vertebrate species. In general they find that computational models are not a reliable method for identifying non-methylated islands. Moreover they find that non-methylated islands are a central epigenetic feature of gene promoters in all vertebrates analysed–including three mammals, a bird, a lizard, a frog and a fish—and not just in warm-blooded vertebrates as suggested by computational models. This shows that the epigenetic function of these non-methylated islands has been conserved over more than 450 million years of evolution. In addition to the non-methylated islands associated with gene promoters, Long et al. identify two other types: intergenic non-methylated islands that are found away from gene promoters and are said to be ‘plastic’ because the DNA in these islands can acquire methyl groups, and ‘broad’ non-methylated islands that span many of the genes that are involved in embryonic development. By showing that the epigenetic role of non-methylated islands has been conserved over time, and identifying three specific types of island, the work of Long et al. marks an important change in our understanding of epigenetics in vertebrates. DOI: http://dx.doi.org/10.7554/eLife.00348.002

Background Polyadenylation, an essential step in eukaryotic gene expression, requires both cis-elements and a plethora of trans- acting polyadenylation factors. The polyadenylation factors are largely conserved across mammals and fungi. The conservation seems also extended to plants based on the analyses of Arabidopsis polyadenylation factors. To extend this observation, we systemically identified the orthologs of yeast and human polyadenylation factors from 10 plant species chosen based on both the availability of their genome sequences and their positions in the evolutionary tree, which render them representatives of different plant lineages. Results The evolutionary trajectories revealed several interesting features of plant polyadenylation factors. First, the number of genes encoding plant polyadenylation factors was clearly increased from “lower” to “higher” plants. Second, the gene expansion in higher plants was biased to some polyadenylation factors, particularly those involved in RNA binding. Finally, while there are clear commonalities, the differences in the polyadenylation apparatus were obvious across different species, suggesting an ongoing process of evolutionary change. These features lead to a model in which the plant polyadenylation complex consists of a conserved core, which is rather rigid in terms of evolutionary conservation, and a panoply of peripheral subunits, which are less conserved and associated with the core in various combinations, forming a collection of somewhat distinct complex assemblies. Conclusions The multiple forms of plant polyadenylation complex, together with the diversified polyA signals may explain the intensive alternative polyadenylation (APA) and its regulatory role in biological functions of higher plants.

The impact of gene silencing on cellular phenotypes is difficult to establish due to the complexity of interactions in the associated biological processes and pathways. A recent genome-wide RNA knock-down study both identified and phenotypically characterized a set of important genes for the cell cycle in HeLa cells. Here, we combine a molecular interaction network analysis, based on physical and functional protein interactions, in conjunction with evolutionary information, to elucidate the common biological and topological properties of these key genes. Our results show that these genes tend to be conserved with their corresponding protein interactions across several species and are key constituents of the evolutionary conserved molecular interaction network. Moreover, a group of bistable network motifs is found to be conserved within this network, which are likely to influence the network stability and therefore the robustness of cellular functioning. They form a cluster, which displays functional homogeneity and this cluster is significantly enriched in genes phenotypically relevant for mitosis. Additional results reveal a relationship between specific cellular processes and the phenotypic outcomes induced by gene silencing. This study introduces new ideas regarding the relationship between genotype and phenotype in the context of the cell cycle. We show that the analysis of molecular interaction networks can result in the identification of genes relevant to cellular processes, which is a promising avenue for future research.

Genomic imprinting is an important epigenetic phenomenon, which on the phenotypic level can be detected by the difference between the two heterozygote classes of a gene. Imprinted genes are important in both the development of the placenta and the embryo, and we hypothesized that imprinted genes might be involved in female fertility traits. We therefore performed an association study for imprinted genes related to female fertility traits in two commercial pig populations. For this purpose, 309 SNPs in fifteen evolutionary conserved imprinted regions were genotyped on 689 and 1050 pigs from the two pig populations. A single SNP association study was used to detect additive, dominant and imprinting effects related to four reproduction traits; total number of piglets born, the number of piglets born alive, the total weight of the piglets born and the total weight of the piglets born alive. Several SNPs showed significant (q-value < 0.10) additive and dominant effects and one SNP showed a significant imprinting effect. The SNP with a significant imprinting effect is closely linked to DIO3, a gene involved in thyroid metabolism. The imprinting effect of this SNP explained approximately 1.6% of the phenotypic variance, which corresponded to approximately 15.5% of the additive genetic variance. In the other population, the imprinting effect of this QTL was not significant (q-value > 0.10), but had a similar effect as in the first population. The results of this study indicate a possible association between the imprinted gene DIO3 and female fertility traits in pigs.

The kidney and inner ear CLC-K chloride channels, which are involved in salt absorption and endolymph production, are regulated by extracellular Ca2+ in the millimolar concentration range. Recently, Gradogna et al. (2010. J. Gen. Physiol. http://dx.doi.org/10.1085/jgp.201010455) identified a pair of acidic residues (E261 and D278) located in the loop between helices I and J as forming a putative intersubunit Ca2+-binding site in hClC-Ka. In this study, we sought to explore the properties of the binding site in more detail. First, we verified that the site is conserved in hClC-Kb and rClC-K1. In addition, we could confer Ca2+ sensitivity to the Torpedo marmorata ClC-0 channel by exchanging its I–J loop with that from ClC-Ka, demonstrating a direct role of the loop in Ca2+ binding. Based on a structure of a bacterial CLC and a new sequence alignment, we built homology models of ClC-Ka. The models suggested additional amino acids involved in Ca2+ binding. Testing mutants of these residues, we could restrict the range of plausible models and positively identify two more residues (E259 and E281) involved in Ca2+ coordination. To investigate cation specificity, we applied extracellular Zn2+, Mg2+, Ba2+, Sr2+, and Mn2+. Zn2+ blocks ClC-Ka as well as its Ca2+-insensitive mutant, suggesting that Zn2+ binds to a different site. Mg2+ does not activate CLC-Ks, but the channels are activated by Ba2+, Sr2+, and Mn2+ with a rank order of potency of Ca2+ > Ba2+ > Sr2+ = Mn2+ for the human CLC-Ks. Dose–response analysis indicates that the less potent Ba2+ has a lower affinity rather than a lower efficacy. Interestingly, rClC-K1 shows an altered rank order (Ca2+ > Sr2+ >> Ba2+), but homology models suggest that residues outside the I–J loop are responsible for this difference. Our detailed characterization of the regulatory Ca2+-binding site provides a solid basis for the understanding of the physiological modulation of CLC-K channel function in the kidney and inner ear.

Monoacylglycerol lipases (MGLs) catalyse the hydrolysis of monoacylglycerol into free fatty acid and glycerol. MGLs have been identified throughout all genera of life and have adopted different substrate specificities depending on their physiological role. In humans, MGL plays an integral part in lipid metabolism affecting energy homeostasis, signalling processes and cancer cell progression. In bacteria, MGLs degrade short-chain monoacylglycerols which are otherwise toxic to the organism. We report the crystal structures of MGL from the bacterium Bacillus sp. H257 (bMGL) in its free form at 1.2 Å and in complex with phenylmethylsulfonyl fluoride at 1.8 Å resolution. In both structures, bMGL adopts an α/β hydrolase fold with a cap in an open conformation. Access to the active site residues, which were unambiguously identified from the protein structure, is facilitated by two different channels. The larger channel constitutes the highly hydrophobic substrate binding pocket with enough room to accommodate monoacylglycerol. The other channel is rather small and resembles the proposed glycerol exit hole in human MGL. Molecular dynamics simulation of bMGL yielded open and closed states of the entrance channel and the glycerol exit hole. Despite differences in the number of residues, secondary structure elements, and low sequence identity in the cap region, this first structure of a bacterial MGL reveals striking structural conservation of the overall cap architecture in comparison with human MGL. Thus it provides insight into the structural conservation of the cap amongst MGLs throughout evolution and provides a framework for rationalising substrate specificities in each organism., Highlights ► First crystal structure of monoacylglycerol lipase from bacterial species. ► Small angle X-ray scattering shows that the protein is a monomer in solution. ► A large hydrophobic channel enables access of the substrate to the active site. ► Molecular dynamic simulations reveal open and closed states of the cap region. ► The cap architecture is conserved on a structural but not on a sequence level.

CD1d-restricted invariant natural killer T cells (NKT cells) have been well characterized in humans and mice, but it is unknown whether they are present in other species. Here we describe the invariant TCR alpha chain and the full length CD1d transcript of pig and horse. Molecular modeling predicts that porcine (po) invariant TCR alpha chain/poCD1d/alpha-GalCer and equine (eq) invariant TCR alpha chain/eqCD1d/alpha-GalCer form complexes that are highly homologous to the human complex. Since a prerequisite for the presence of NKT cells is the expression of CD1d protein, we performed searches for CD1D genes and CD1d transcripts in multiple species. Previously, cattle and guinea pig have been suggested to lack CD1D genes. The CD1D genes of European taurine cattle (Bos taurus) are known to be pseudogenes because of disrupting mutations in the start codon and in the donor splice site of the first intron. Here we show that the same mutations are found in six other ruminants: African buffalo, sheep, bushbuck, bongo, N'Dama cattle, and roe deer. In contrast, intact CD1d transcripts were found in guinea pig, African elephant, horse, rabbit, and pig. Despite the discovery of a highly homologous NKT/CD1d system in pig and horse, our data suggest that functional CD1D and CD1d-restricted NKT cells are not universally present in mammals. (C) 2009 Elsevier Ltd. All rights reserved.

Divergent morphology of species has largely been ascribed to genetic differences in the tissue-specific expression of proteins, which could be achieved by divergence in cis-regulatory elements or by altering the binding specificity of transcription factors (TFs). The relative importance of the latter has been difficult to assess, as previous systematic analyses of TF binding specificity have been performed using different methods in different species. To address this, we determined the binding specificities of 242 Drosophila TFs, and compared them to human and mouse data. This analysis revealed that TF binding specificities are highly conserved between Drosophila and mammals, and that for orthologous TFs, the similarity extends even to the level of very subtle dinucleotide binding preferences. The few human TFs with divergent specificities function in cell types not found in fruit flies, suggesting that evolution of TF specificities contributes to emergence of novel types of differentiated cells. DOI: http://dx.doi.org/10.7554/eLife.04837.001, eLife digest Flies look very different from humans, but both are descended from a common ancestor that existed over 600 million years ago. Some differences between animal species are due to them having different genes: stretches of DNA that contain the instructions to make proteins and other molecules. However, often differences are caused by the same or similar genes being switched on and off at different times and in different tissues in each species. The instructions that control when and where a gene is expressed are written in the sequence of DNA bases located in the regulatory region of the gene. These instructions are written in a language that is often called the ‘gene regulatory code’. This code is read and interpreted by proteins called transcription factors that bind to specific sequences of DNA (or ‘DNA words’) and increase or decrease gene expression. Changes in gene expression between species could therefore be due to changes in the transcription factors and/or changes in the instructions within the regulatory regions of specific genes. Gene regulatory regions are not well conserved between species. However, it is unclear if the instructions in these regions are written using the same gene regulatory code, and whether transcription factors found in different species recognize different DNA words. Nitta et al. have now used high-throughput methods to identify the DNA words recognized by 242 transcription factors from a fruit fly called Drosophila melanogaster. Nitta et al. then used new computational tools to find motifs, or collections of DNA words, that are recognized by each of the transcription factors. By comparing the motifs, they observed that, in spite of more than 600 million years of evolution, almost all known motifs found in humans and mice were recognized by fruit fly transcription factors. Nitta et al. noted that both fruit flies and humans have transcription factors that recognize a few unique motifs, and confer properties that are specific to each species. For example, some of the transcription factors that control the development of the fruit fly wing are not present in humans. Moreover, fruit flies lack both mucus-producing goblet cells and the ability to recognize a motif read by the transcription factor that controls the development of these cells in humans. The findings of Nitta et al. also indicate that transcription factors do not evolve to recognize subtly different DNA motifs, but instead appear constrained to recognize the same motifs. Thus, much like the genetic code that instructs how to build proteins, the gene regulatory code that determines how DNA sequences direct gene expression is also highly conserved in animals. The language used to guide the development of animals has, as such, remained very similar for millions of years. What makes animals different is differences in the content and length of the instructions that are written using this language into the regulatory regions of their genes. DOI: http://dx.doi.org/10.7554/eLife.04837.002

Submitted by Paulo Gavaia (pgavaia@ualg.pt) on 2014-05-30T10:20:11Z No. of bitstreams: 1 Identification of a Promoter Element within the Zebrafish colXα1 Gene Responsive to Runx2 Isoforms Osf2 Cbfa1 and til-1 but not to pebp2αA2.pdf: 688008 bytes, checksum: 267bf20255e71f87b6bfae9859067b5d (MD5) Approved for entry into archive by Pedro Sousa (a19661@ualg.pt) on 2014-07-14T10:44:50Z (GMT) No. of bitstreams: 2 Identification of a Promoter Element within the Zebrafish colXα1 Gene Responsive to Runx2 Isoforms Osf2 Cbfa1 and til-1 but not to pebp2αA2.pdf: 688008 bytes, checksum: 267bf20255e71f87b6bfae9859067b5d (MD5) 8753041680192686.zip: 577095 bytes, checksum: 020c5747d54e09c1fe2342995eec56b5 (MD5) Made available in DSpace on 2014-07-14T10:44:51Z (GMT). No. of bitstreams: 2 Identification of a Promoter Element within the Zebrafish colXα1 Gene Responsive to Runx2 Isoforms Osf2 Cbfa1 and til-1 but not to pebp2αA2.pdf: 688008 bytes, checksum: 267bf20255e71f87b6bfae9859067b5d (MD5) 8753041680192686.zip: 577095 bytes, checksum: 020c5747d54e09c1fe2342995eec56b5 (MD5) Previous issue date: 2006

The human gene RPMS12 encodes a protein similar to bacterial ribosomal protein S12 and is proposed to represent the human mitochondrial orthologue, RPMS12 reporter gene expression in cultured human cells supports the idea that the gene product is mitochondrial and is localized to the inner membrane. Human cells contain at least four structurally distinct RPMS12 mRNAs that differ in their 5'-untranslated region (5'-UTR) as a result of alternate splicing and of 5' end heterogeneity. All of them encode the same polypeptide, The full 5'-UTR contains two types of sequence element implicated elsewhere in translational regulation as follows: a short upstream open reading frame and an oligopyrimidine tract similar to that found at the 5' end of mRNAs encoding other growth-regulated proteins, including those of cytosolic ribosomes, The fully spliced (short) mRNA is the predominant form in all cell types studied and is translationally down-regulated in cultured cells in response to serum starvation, even though it lacks both of the putative translational regulatory elements. By contrast, other splice variants containing one or both of these elements are not translationally regulated by growth status but are translated poorly in both growing and non-growing cells. Reporter analysis identified a 26-nucleotide tract of the 5'-UTR of the short mRNA that is essential for translational down-regulation in growth-inhibited cells. Such experiments also confirmed that the 5'-UTR of the longer mRNA variants contains negative regulatory elements for translation. Tissue representation of RPMS12 mRNA is highly variable, following a typical mitochondrial pattern, but the relative levels of the different splice variants are similar in different tissues. These findings indicate a complex, multilevel regulation of RPMS12 gene expression in response to signals mediating growth, tissue specialization, and probably metabolic needs. ispartof: Journal of biological chemistry vol:274 issue:45 pages:31853-31862 ispartof: location:United States status: published

Drosophila type 2 Btk29A reveals the highest homology to Btk among mammalian Tec kinases. In Btk29AficP mutant males, the apodeme holding the penis split into two pieces. Human Btk rescued this phenotype in 39% of Btk29AficP males, while the Drosophila transgenes did so in 90–100% of mutants. The Btk29AficP mutation reduced adult longevity to 11% that of wild-type. This effect was counteracted by Drosophila type 2, yielding 76% of the wild-type lifespan. Human Btk extended the lifespan of Btk29AficP mutants only to 20% that of wild-type. Thus human Btk can partially replace Drosophila Btk29A+ in male genital development and survival.