The animal in the genome: comparative genomics and evolution

Deciphering a shared transcriptomic regulation and the relative contribution of each regulator type through endometrial gene expression signatures

BACKGORUND

While various endometrial biomarkers have been characterized at the transcriptomic and functional level, there is generally a poor overlap among studies, making it unclear to what extent their upstream regulators (e.g., ovarian hormones, transcription factors (TFs) and microRNAs (miRNAs)) realistically contribute to menstrual cycle progression and function. Unmasking the intricacies of the molecular interactions in the endometrium from a novel systemic point of view will help gain a more accurate perspective of endometrial regulation and a better explanation the molecular etiology of endometrial-factor infertility.

METHODS

An in-silico analysis was carried out to identify which regulators consistently target the gene biomarkers proposed in studies related to endometrial progression and implantation failure (19 gene lists/signatures were included). The roles of these regulators, and of genes related to progesterone and estrogens, were then analysed in transcriptomic datasets compiled from samples collected throughout the menstrual cycle (n = 129), and the expression of selected TFs were prospectively validated in an independent cohort of healthy participants (n = 19).

RESULTS

A total of 3,608 distinct genes from the 19 gene lists were associated with endometrial progression and implantation failure. The lists' regulation was significantly favoured by TFs (89% (17/19) of gene lists) and progesterone (47% (8 /19) of gene lists), rather than miRNAs (5% (1/19) of gene lists) or estrogen (0% (0/19) of gene lists), respectively (FDR < 0.05). Exceptionally, two gene lists that were previously associated with implantation failure and unexplained infertility were less hormone-dependent, but primarily regulated by estrogen. Although endometrial progression genes were mainly targeted by hormones rather than non-hormonal contributors (odds ratio = 91.94, FDR < 0.05), we identified 311 TFs and 595 miRNAs not previously associated with ovarian hormones. We highlight CTCF, GATA6, hsa-miR-15a-5p, hsa-miR-218-5p, hsa-miR-107, hsa-miR-103a-3p, and hsa-miR-128-3p, as overlapping novel master regulators of endometrial function. The gene expression changes of selected regulators throughout the menstrual cycle (FDR < 0.05), dually validated in-silico and through endometrial biopsies, corroborated their potential regulatory roles in the endometrium.

CONCLUSIONS

This study revealed novel hormonal and non-hormonal regulators and their relative contributions to endometrial progression and pathology, providing new leads for the potential causes of endometrial-factor infertility.

Diversification and Functional Evolution of HOX Proteins

Gene duplication and divergence is a major contributor to the generation of morphological diversity and the emergence of novel features in vertebrates during evolution. The availability of sequenced genomes has facilitated our understanding of the evolution of genes and regulatory elements. However, progress in understanding conservation and divergence in the function of proteins has been slow and mainly assessed by comparing protein sequences in combination with in vitro analyses. These approaches help to classify proteins into different families and sub-families, such as distinct types of transcription factors, but how protein function varies within a gene family is less well understood. Some studies have explored the functional evolution of closely related proteins and important insights have begun to emerge. In this review, we will provide a general overview of gene duplication and functional divergence and then focus on the functional evolution of HOX proteins to illustrate evolutionary changes underlying diversification and their role in animal evolution.

Transcriptional changes through menstrual cycle reveal a global transcriptional derepression underlying the molecular mechanism involved in the window of implantation

The human endometrium is a dynamic tissue that only is receptive to host the embryo during a brief time in the middle secretory phase, called the window of implantation (WOI). Despite its importance, regulation of the menstrual cycle remains incompletely understood. The aim of this study was to characterize the gene cooperation and regulation of menstrual cycle progression, to dissect the molecular complexity underlying acquisition of endometrial receptivity for a successful pregnancy, and to provide the scientific community with detailed gene co-expression information throughout the menstrual cycle on a user-friendly web-tool database. A retrospective gene co-expression analysis was performed based on the endometrial receptivity array (ERarray) gene signature from 523 human endometrial samples collected across the menstrual cycle, including during the WOI. Gene co-expression analysis revealed the WOI as having the significantly smallest proportion of negative correlations for transcriptional profiles associated with successful pregnancies compared to other cycle stages, pointing to a global transcriptional derepression being involved in acquisition of endometrial receptivity. Regulation was greatest during the transition between proliferative and secretory endometrial phases. Further, we prioritized nuclear hormone receptors as major regulators of this derepression and proved that some genes and transcription factors involved in this process were dysregulated in patients with recurrent implantation failure. We also compiled the wealth of gene co-expression data to stimulate hypothesis-driven single-molecule endometrial studies in a user-friendly database: Menstrual Cycle Gene Co-expression Network (www.menstrualcyclegcn.com). This study revealed a global transcriptional repression across the menstrual cycle, which relaxes when the WOI opens for transcriptional profiles associated with successful pregnancies. These findings suggest that a global transcriptional derepression is needed for embryo implantation and early development.

Multivariate analysis of genomic variables, effective population size, and mutation rate

Objective

The relationship between genomic variables (genome size, gene number, intron size, and intron number) and evolutionary forces has two implications. First, they help to unravel the mechanism underlying genome evolution. Second, they provide a solution to the debate over discrepancy between genome size variation and organismal complexity. Previously, a clear correlation between genomic variables and effective population size and mutation rate (Neu) led to an important hypothesis to consider random genetic drift as a major evolutionary force during evolution of genome size and complexity. But recent reports also support natural selection as the leading evolutionary force. As such, the debate remains unresolved.

Results

Here, we used a multivariate method to explore the relationship between genomic variables and Neu in order to understand the evolution of genome. Previously reported patterns between genomic variables and Neu were not observed in our multivariate study. We found only one association between intron number and Neu, but no relationships were observed between genome size, intron size, gene number, and Neu, suggesting that Neu of the organisms solely does not influence genome evolution. We, therefore, concluded that Neu influences intron evolution, while it may not be the only force that provides mechanistic insights into genome evolution and complexity.

Electronic supplementary material

The online version of this article (10.1186/s13104-019-4097-3) contains supplementary material, which is available to authorized users.

Perspectives on the history of evo-devo and the contemporary research landscape in the genomics era

A fundamental question in biology is how the extraordinary range of living organisms arose. In this theme issue, we celebrate how evolutionary studies on the origins of morphological diversity have changed over the past 350 years since the first publication of the Philosophical Transactions of The Royal Society Current understanding of this topic is enriched by many disciplines, including anatomy, palaeontology, developmental biology, genetics and genomics. Development is central because it is the means by which genetic information of an organism is translated into morphology. The discovery of the genetic basis of development has revealed how changes in form can be inherited, leading to the emergence of the field known as evolutionary developmental biology (evo-devo). Recent approaches include imaging, quantitative morphometrics and, in particular, genomics, which brings a new dimension. Articles in this issue illustrate the contemporary evo-devo field by considering general principles emerging from genomics and how this and other approaches are applied to specific questions about the evolution of major transitions and innovations in morphology, diversification and modification of structures, intraspecific morphological variation and developmental plasticity. Current approaches enable a much broader range of organisms to be studied, thus building a better appreciation of the origins of morphological diversity.This article is part of the themed issue 'Evo-devo in the genomics era, and the origins of morphological diversity'.

The Unicellular Ancestry of Groucho-Mediated Repression and the Origins of Metazoan Transcription Factors

Groucho is a co-repressor that interacts with many transcription factors playing a crucial role in animal development. The evolutionary origins of Groucho are not clear. It is generally regarded as being a distinct animal-specific protein, although with similarities to the yeast Tup-like proteins. Here, it is shown that Groucho has true orthologs in unicellular relatives of animals. Based on their phylogenetic distribution, and an analysis of ligand-binding residues, these genes are unlikely to be orthologs of the fungal Tup-like genes. By identifying conserved candidate Groucho interaction motifs (GIMs) in nonmetazoan transcription factors, it is demonstrated that the details of molecular interactions between Groucho and transcription factors are likely to have been established prior to the origin of animals, but that the association of GIMs with many transcription factor types can be regarded as a metazoan innovation.

Metabolic Field (Schrodinger); an explanatory platform for biology: Based on lecture at Trinity College, Dublin, Ireland, July 18, 2012

Metabolism represents the nexus of fundamental physical forces, which while present in all structure and function require new explanatory emergent principles, which, so far, cannot be predicted or derived solely from description of chemistry and physics. Metabolism is essentially concerned with the transduction of energy flows with respect to time, space, and matter. Language models and metaphors contribute to construction of scientific explanation within biology. The concept of a metabolic field yields a deeper, broader, more quantitative integrated theoretical framework leading to novel predictive models of systems biology.

The relationship between gene isoform multiplicity, number of exons and protein divergence

At present we know that phenotypic differences between organisms arise from a variety of sources, like protein sequence divergence, regulatory sequence divergence, alternative splicing, etc. However, we do not have yet a complete view of how these sources are related. Here we address this problem, studying the relationship between protein divergence and the ability of genes to express multiple isoforms. We used three genome-wide datasets of human-mouse orthologs to study the relationship between isoform multiplicity co-occurrence between orthologs (the fact that two orthologs have more than one isoform) and protein divergence. In all cases our results showed that there was a monotonic dependence between these two properties. We could explain this relationship in terms of a more fundamental one, between exon number of the largest isoform and protein divergence. We found that this last relationship was present, although with variations, in other species (chimpanzee, cow, rat, chicken, zebrafish and fruit fly). In summary, we have identified a relationship between protein divergence and isoform multiplicity co-occurrence and explained its origin in terms of a simple gene-level property. Finally, we discuss the biological implications of these findings for our understanding of inter-species phenotypic differences.

Ancestral and conserved cis-regulatory architectures in developmental control genes

Among developmental control genes, transcription factor-target gene "linkages"--the direct connections between target genes and the factors that control their patterns of expression--can show remarkable evolutionary stability. However, the specific binding sites that mediate and define these regulatory connections are themselves often subject to rapid turnover. Here we describe several instances in which particular transcription factor binding motif combinations have evidently been conserved upstream of orthologous target genes for extraordinarily long evolutionary periods. This occurs against a backdrop in which other binding sites for the same factors are coming and going rapidly. Our examples include a particular Dpp Silencer Element upstream of insect brinker genes, in combination with a novel motif we refer to as the Downstream Element; combinations of a Suppressor of Hairless Paired Site (SPS) and a specific proneural protein binding site associated with arthropod Notch pathway target genes; and a three-motif combination, also including an SPS, upstream of deuterostome Hes repressor genes, which are also Notch targets. We propose that these stable motif architectures have been conserved intact from a deep ancestor, in part because they mediate a special mode of regulation that cannot be supplied by the other, unstable motif instances.

Correlating traits of gene retention, sequence divergence, duplicability and essentiality in vertebrates, arthropods, and fungi

Delineating ancestral gene relations among a large set of sequenced eukaryotic genomes allowed us to rigorously examine links between evolutionary and functional traits. We classified 86% of over 1.36 million protein-coding genes from 40 vertebrates, 23 arthropods, and 32 fungi into orthologous groups and linked over 90% of them to Gene Ontology or InterPro annotations. Quantifying properties of ortholog phyletic retention, copy-number variation, and sequence conservation, we examined correlations with gene essentiality and functional traits. More than half of vertebrate, arthropod, and fungal orthologs are universally present across each lineage. These universal orthologs are preferentially distributed in groups with almost all single-copy or all multicopy genes, and sequence evolution of the predominantly single-copy orthologous groups is markedly more constrained. Essential genes from representative model organisms, Mus musculus, Drosophila melanogaster, and Saccharomyces cerevisiae, are significantly enriched in universal orthologs within each lineage, and essential-gene-containing groups consistently exhibit greater sequence conservation than those without. This study of eukaryotic gene repertoire evolution identifies shared fundamental principles and highlights lineage-specific features, it also confirms that essential genes are highly retained and conclusively supports the "knockout-rate prediction" of stronger constraints on essential gene sequence evolution. However, the distinction between sequence conservation of single- versus multicopy orthologs is quantitatively more prominent than between orthologous groups with and without essential genes. The previously underappreciated difference in the tolerance of gene duplications and contrasting evolutionary modes of "single-copy control" versus "multicopy license" may reflect a major evolutionary mechanism that allows extended exploration of gene sequence space.

Functional genomics and networks: new approaches in the extraction of complex gene modules

The engine that makes the cell work is made of an intricate network of molecular interactions. Nowadays, the elements and relationships of this complex network can be studied with several types of high-throughput techniques. The dream of having a global picture of the cell from different perspectives that can jointly explain cell behavior is, at least technically, feasible. However, this task can only be accomplished by filling the gap between data and information. The availability of methods capable of accurately managing, integrating and analyzing the results from these experiments is crucial for this purpose. Here, we review the new challenges raised by the availability of different genomic data, as well as the new proposals presented to cope with the increasing data complexity. Special emphasis is given to approaches that explore the transcriptome trying to describe the modules of genes that account for the traits studied.

Alzheimer's disease: insights from Drosophila melanogaster models

The power of fruit fly genetics is being deployed against some of the most intractable and economically significant problems in modern medicine, the neurodegenerative diseases. Fly models of Alzheimer's disease can be exposed to the rich diversity of biological techniques that are available to the community and are providing new insights into disease mechanisms, and assisting in the identification of novel targets for therapy. Similar approaches might also help us to interpret the results of genome-wide association studies of human neurodegenerative diseases by allowing us to triage gene “hits” according to whether a candidate risk factor gene has a modifying effect on the disease phenotypes in fly model systems.

Evidence of uneven selective pressure on different subsets of the conserved human genome; implications for the significance of intronic and intergenic DNA

Background

Human genetic variation produces the wide range of phenotypic differences that make us individual. However, little is known about the distribution of variation in the most conserved functional regions of the human genome. We examined whether different subsets of the conserved human genome have been subjected to similar levels of selective constraint within the human population. We used set theory and high performance computing to carry out an analysis of the density of Single Nucleotide Polymorphisms (SNPs) within the evolutionary conserved human genome, at three different selective stringencies, intersected with exonic, intronic and intergenic coordinates.

Results

We demonstrate that SNP density across the genome is significantly reduced in conserved human sequences. Unexpectedly, we further demonstrate that, despite being conserved to the same degree, SNP density differs significantly between conserved subsets. Thus, both the conserved exonic and intronic genomes contain a significantly reduced density of SNPs compared to the conserved intergenic component. Furthermore the intronic and exonic subsets contain almost identical densities of SNPs indicating that they have been constrained to the same degree.

Conclusion

Our findings suggest the presence of a selective linkage between the exonic and intronic subsets and ascribes increased significance to the role of introns in human health. In addition, the identification of increased plasticity within the conserved intergenic subset suggests an important role for this subset in the adaptation and diversification of the human population.

Early origin of the bilaterian developmental toolkit

Whole-genome sequences from the choanoflagellate Monosiga brevicollis, the placozoan Trichoplax adhaerens and the cnidarian Nematostella vectensis have confirmed results from comparative evolutionary developmental studies that much of the developmental toolkit once thought to be characteristic of bilaterians appeared much earlier in the evolution of animals. The diversity of transcription factors and signalling pathway genes in animals with a limited number of cell types and a restricted developmental repertoire is puzzling, particularly in light of claims that such highly conserved elements among bilaterians provide evidence of a morphologically complex protostome-deuterostome ancestor. Here, I explore the early origination of elements of what became the bilaterian toolkit, and suggest that placozoans and cnidarians represent a depauperate residue of a once more diverse assemblage of early animals, some of which may be represented in the Ediacaran fauna (c. 585-542 Myr ago).

Apprehending multicellularity: regulatory networks, genomics, and evolution

The genomic revolution has provided the first glimpses of the architecture of regulatory networks. Combined with evolutionary information, the "network view" of life processes leads to remarkable insights into how biological systems have been shaped by various forces. This understanding is critical because biological systems, including regulatory networks, are not products of engineering but of historical contingencies. In this light, we attempt a synthetic overview of the natural history of regulatory networks operating in the development and differentiation of multicellular organisms. We first introduce regulatory networks and their organizational principles as can be deduced using ideas from the graph theory. We then discuss findings from comparative genomics to illustrate the effects of lineage-specific expansions, gene-loss, and nonprotein-coding DNA on the architecture of networks. We consider the interaction between expansions of transcription factors, and cis regulatory and more general chromatin state stabilizing elements in the emergence of morphological complexity. Finally, we consider a case study of the Notch subnetwork, which is present throughout Metazoa, to examine how such a regulatory system has been pieced together in evolution from new innovations and pre-existing components that were originally functionally distinct.

Functional and comparative genomics analyses of pmp22 in medaka fish

Background

Pmp22, a member of the junction protein family Claudin/EMP/PMP22, plays an important role in myelin formation. Increase of pmp22 transcription causes peripheral neuropathy, Charcot-Marie-Tooth disease type1A (CMT1A). The pathophysiological phenotype of CMT1A is aberrant axonal myelination which induces a reduction in nerve conduction velocity (NCV). Several CMT1A model rodents have been established by overexpressing pmp22. Thus, it is thought that pmp22 expression must be tightly regulated for correct myelin formation in mammals. Interestingly, the myelin sheath is also present in other jawed vertebrates. The purpose of this study is to analyze the evolutionary conservation of the association between pmp22 transcription level and vertebrate myelin formation, and to find the conserved non-coding sequences for pmp22 regulation by comparative genomics analyses between jawed fishes and mammals.

Results

A transgenic pmp22 over-expression medaka fish line was established. The transgenic fish had approximately one fifth the peripheral NCV values of controls, and aberrant myelination of transgenic fish in the peripheral nerve system (PNS) was observed. We successfully confirmed that medaka fish pmp22 has the same exon-intron structure as mammals, and identified some known conserved regulatory motifs. Furthermore, we found novel conserved sequences in the first intron and 3'UTR.

Conclusion

Medaka fish undergo abnormalities in the PNS when pmp22 transcription increases. This result indicates that an adequate pmp22 transcription level is necessary for correct myelination of jawed vertebrates. Comparison of pmp22 orthologs between distantly related species identifies evolutionary conserved sequences that contribute to precise regulation of pmp22 expression.

Comparative functional analysis of the Caenorhabditis elegans and Drosophila melanogaster proteomes

The nematode Caenorhabditis elegans is a popular model system in genetics, not least because a majority of human disease genes are conserved in C. elegans. To generate a comprehensive inventory of its expressed proteome, we performed extensive shotgun proteomics and identified more than half of all predicted C. elegans proteins. This allowed us to confirm and extend genome annotations, characterize the role of operons in C. elegans, and semiquantitatively infer abundance levels for thousands of proteins. Furthermore, for the first time to our knowledge, we were able to compare two animal proteomes (C. elegans and Drosophila melanogaster). We found that the abundances of orthologous proteins in metazoans correlate remarkably well, better than protein abundance versus transcript abundance within each organism or transcript abundances across organisms; this suggests that changes in transcript abundance may have been partially offset during evolution by opposing changes in protein abundance.

Proteins are the active players that execute the genetic program of a cell, and their levels and interactions are precisely controlled. Routinely monitoring thousands of proteins is difficult, as they can be present at vastly different abundances, come with various sizes, shapes, and charge, and have a more complex alphabet of twenty “letters,” in contrast to the four letters of the genome itself. Here, we used mass spectrometry to extensively characterize the proteins of a popular model organism, the nematode Caenorhabditis elegans. Together with previous data from the fruit fly Drosophila melanogaster, this allows us to compare the protein levels of two animals on a global scale. Surprisingly, we find that individual protein abundance is highly conserved between the two species. So, although worms and flies look very different, they need similar amounts of each conserved, orthologous protein. Because many C. elegans and D. melanogaster proteins also have counterparts in humans, our results suggest that similar rules may apply to our own proteins.

A quantitative comparison of two animal proteomes shows a striking correlation of protein abundance levels, a better correlation than transcript levels. Are the latter more variable during evolution?

Estimating the size of the human interactome

After the completion of the human and other genome projects it emerged that the number of genes in organisms as diverse as fruit flies, nematodes, and humans does not reflect our perception of their relative complexity. Here, we provide reliable evidence that the size of protein interaction networks in different organisms appears to correlate much better with their apparent biological complexity. We develop a stable and powerful, yet simple, statistical procedure to estimate the size of the whole network from subnet data. This approach is then applied to a range of eukaryotic organisms for which extensive protein interaction data have been collected and we estimate the number of interactions in humans to be approximately 650,000. We find that the human interaction network is one order of magnitude bigger than the Drosophila melanogaster interactome and approximately 3 times bigger than in Caenorhabditis elegans.

The evolution of the animals: introduction to a Linnean tercentenary celebration

Celebrating 300 years since the birth of Carl Linnaeus (1707-1778), a meeting was held in June 2007 to review recent progress made in understanding the origins and evolutionary radiation of the animals. The year 2008 celebrates the 250th anniversary of the publication of the 10th edition of Linnaeus' Systema Naturae, generally considered to be the starting point of zoological nomenclature. With subsequent advances in comparative taxonomic and systematic studies, Darwin's discovery of evolution by natural selection, the birth of phylogenetic systematics, and the wider interest in biodiversity, it is salutary to consider that many of the major advances in our understanding of animal evolution have been made in recent years. Phylogenetic systematics, drawing from evidence provided by genotype, phenotype and an understanding of the link between them through comparative embryological and evolutionary developmental studies, has provided a wide consensus of the major branching patterns of the tree of life. More importantly, the integrated approaches discussed in the 16 contributions to this volume highlight the identity and nature of problematic taxa, the missing data, errors in existing analytical procedures and the promise of a wealth of additional characters from genomes that need to be accumulated and assessed in providing a definitive Systema Naturae.