Genome-wide detection of gene extinction in early mammalian evolution

Enigmatic Nodal and Lefty gene repertoire discrepancy: Latent evolutionary history revealed by vertebrate-wide phylogeny

Homology in vertebrate body plans is traditionally ascribed to the high-level conservation of regulatory components within the genetic programs governing them, particularly during the "phylotypic stage." However, advancements in embryology and molecular phylogeny have unveiled the dynamic nature of gene repertoires responsible for early development. Notably, the Nodal and Lefty genes, members of the transforming growth factor-beta superfamily producing intercellular signaling molecules and crucial for left-right (L-R) symmetry breaking, exhibit distinctive features within their gene repertoires. These features encompass among-species gene repertoire variations resulting from gene gain and loss, as well as gene conversion. Despite their significance, these features have been largely unexplored in a phylogenetic context, but accumulating genome-wide sequence information is allowing the scrutiny of these features. It has exposed hidden paralogy between Nodal1 and Nodal2 genes resulting from differential gene loss in amniotes. In parallel, the tandem cluster of Lefty1 and Lefty2 genes, which was thought to be confined to mammals, is observed in sharks and rays, with an unexpected phylogenetic pattern. This article provides a comprehensive review of the current understanding of the origins of these vertebrate gene repertoires and proposes a revised nomenclature based on the elucidated history of vertebrate genome evolution.

Evolutionary Turnover in Wnt Gene Expression but Conservation of Wnt Signaling during Ovary Determination in a TSD Reptile

Temperature-dependent sex determination (TSD) is a well-known characteristic of many reptilian species. However, the molecular processes linking ambient temperature to determination of gonad fate remain hazy. Here, we test the hypothesis that Wnt expression and signaling differ between female- and male-producing temperatures in the snapping turtle Chelydra serpentina. Canonical Wnt signaling involves secretion of glycoproteins called WNTs, which bind to and activate membrane bound receptors that trigger β-catenin stabilization and translocation to the nucleus where β-catenin interacts with TCF/LEF transcription factors to regulate expression of Wnt targets. Non-canonical Wnt signaling occurs via 2 pathways that are independent of β-catenin: one involves intracellular calcium release (the Wnt/Ca2+ pathway), while the other involves activation of RAC1, JNK, and RHOA (the Wnt/planar cell polarity pathway). We screened 20 Wnt genes for differential expression between female- and male-producing temperatures during sex determination in the snapping turtle. Exposure of embryos to the female-producing temperature decreased expression of 7 Wnt genes but increased expression of 2 Wnt genes and Rspo1 relative to embryos at the male-producing temperature. Temperature also regulated expression of putative Wnt target genes in vivo and a canonical Wnt reporter (6x TCF/LEF sites drive H2B-GFP expression) in embryonic gonadal cells in vitro. Results indicate that Wnt signaling was higher at the female- than at the male-producing temperature. Evolutionary analyses of all 20 Wnt genes revealed that thermosensitive Wnts, as opposed to insensitive Wnts, were less likely to show evidence of positive selection and experienced stronger purifying selection within TSD species.

Shark and ray genomics for disentangling their morphological diversity and vertebrate evolution

Developmental studies of sharks and rays (elasmobranchs) have provided much insight into the process of morphological evolution of vertebrates. Although those studies are supposedly fueled by large-scale molecular sequencing information, whole-genome sequences of sharks and rays were made available only recently. One compelling difficulty of elasmobranch developmental biology is the low accessibility to embryonic study materials and their slow development. Another limiting factor is the relatively large size of their genomes. Moreover, their large body sizes restrict sustainable captive breeding, while their high body fluid osmolarity prevents reproducible cell culturing for in vitro experimentation, which has also limited our knowledge of their chromosomal organization for validation of genome sequencing products. This article focuses on egg-laying elasmobranch species used in developmental biology and provides an overview of the characteristics of the shark and ray genomes revealed to date. Developmental studies performed on a gene-by-gene basis are also reviewed from a whole-genome perspective. Among the popular regulatory genes studied in developmental biology, I scrutinize shark homologs of Wnt genes that highlight vanishing repertoires in many other vertebrate lineages, as well as Hox genes that underwent an unexpected modification unique to the elasmobranch lineage. These topics are discussed together with insights into the reconstruction of developmental programs in the common ancestor of vertebrates and its subsequent evolutionary trajectories that mark the features that are unique to, and those characterizing the diversity among, cartilaginous fishes.

Discovery of four Noggin genes in lampreys suggests two rounds of ancient genome duplication

The secreted protein Noggin1 was the first discovered natural embryonic inducer produced by cells of the Spemann organizer. Thereafter, it was shown that vertebrates have a whole family of Noggin genes with different expression patterns and functional properties. For example, Noggin1 and Noggin2 inhibit the activity of BMP, Nodal/Activin and Wnt-beta-catenin signalling, while Noggin4 cannot suppress BMP but specifically modulates Wnt signalling. In this work, we described and investigated phylogeny and expression patterns of four Noggin genes in lampreys, which represent the most basally divergent group of extant vertebrates, the cyclostomes, belonging to the superclass Agnatha. Assuming that lampreys have Noggin homologues in all representatives of another superclass of vertebrates, the Gnathostomata, we propose a model for Noggin family evolution in vertebrates. This model is in agreement with the hypotheses suggesting two rounds of genome duplication in the ancestor of vertebrates before the divergence of Agnatha and Gnathostomata.

Ermakova et al. report four Noggin genes in lampreys and using phylogenetics, gene synteny analysis, and in situ hybridization, suggest that the Noggin gene underwent two rounds of duplication and evolved specific functions before the divergence of vertebrate and lamprey lineages. These findings offer insight into early vertebrate genome and developmental evolution.

Robustness to extinction and plasticity derived from mutualistic bipartite ecological networks

Understanding the response of ecological networks to perturbations and disruptive events is needed to anticipate the biodiversity loss and extinction cascades. Here, we study how network plasticity reshapes the topology of mutualistic networks in response to species loss. We analyze more than one hundred empirical mutualistic networks and considered random and targeted removal as mechanisms of species extinction. Network plasticity is modeled as either random rewiring, as the most parsimonious approach, or resource affinity-driven rewiring, as a proxy for encoding the phylogenetic similarity and functional redundancy among species. This redundancy should be positively correlated with the robustness of an ecosystem, as functions can be taken by other species once one of them is extinct. We show that effective modularity, i.e. the ability of an ecosystem to adapt or restructure, increases with increasing numbers of extinctions, and with decreasing the replacement probability. Importantly, modularity is mostly affected by the extinction rather than by rewiring mechanisms. These changes in community structure are reflected in the robustness and stability due to their positive correlation with modularity. Resource affinity-driven rewiring offers an increase of modularity, robustness, and stability which could be an evolutionary favored mechanism to prevent a cascade of co-extinctions.

Nodal paralogues underlie distinct mechanisms for visceral left-right asymmetry in reptiles and mammals

Unidirectional fluid flow generated by motile cilia at the left-right organizer (LRO) breaks left-right (L-R) symmetry during early embryogenesis in mouse, frog and zebrafish. The chick embryo, however, does not require motile cilia for L-R symmetry breaking. The diversity of mechanisms for L-R symmetry breaking among vertebrates and the trigger for such symmetry breaking in non-mammalian amniotes have remained unknown. Here we examined how L-R asymmetry is established in two reptiles, Madagascar ground gecko and Chinese softshell turtle. Both of these reptiles appear to lack motile cilia at the LRO. The expression of the Nodal gene at the LRO in the reptilian embryos was found to be asymmetric, in contrast to that in vertebrates such as mouse that are dependent on cilia for L-R patterning. Two paralogues of the Nodal gene derived from an ancient gene duplication are retained and expressed differentially in cilia-dependent and cilia-independent vertebrates. The expression of these two Nodal paralogues is similarly controlled in the lateral plate mesoderm but regulated differently at the LRO. Our in-depth analysis of reptilian embryos thus suggests that mammals and non-mammalian amniotes deploy distinct strategies dependent on different Nodal paralogues for rendering Nodal activity asymmetric at the LRO.

Discovery of high-confidence human protein-coding genes and exons by whole-genome PhyloCSF helps elucidate 118 GWAS loci

The most widely appreciated role of DNA is to encode protein, yet the exact portion of the human genome that is translated remains to be ascertained. We previously developed PhyloCSF, a widely used tool to identify evolutionary signatures of protein-coding regions using multispecies genome alignments. Here, we present the first whole-genome PhyloCSF prediction tracks for human, mouse, chicken, fly, worm, and mosquito. We develop a workflow that uses machine learning to predict novel conserved protein-coding regions and efficiently guide their manual curation. We analyze more than 1000 high-scoring human PhyloCSF regions and confidently add 144 conserved protein-coding genes to the GENCODE gene set, as well as additional coding regions within 236 previously annotated protein-coding genes, and 169 pseudogenes, most of them disabled after primates diverged. The majority of these represent new discoveries, including 70 previously undetected protein-coding genes. The novel coding genes are additionally supported by single-nucleotide variant evidence indicative of continued purifying selection in the human lineage, coding-exon splicing evidence from new GENCODE transcripts using next-generation transcriptomic data sets, and mass spectrometry evidence of translation for several new genes. Our discoveries required simultaneous comparative annotation of other vertebrate genomes, which we show is essential to remove spurious ORFs and to distinguish coding from pseudogene regions. Our new coding regions help elucidate disease-associated regions by revealing that 118 GWAS variants previously thought to be noncoding are in fact protein altering. Altogether, our PhyloCSF data sets and algorithms will help researchers seeking to interpret these genomes, while our new annotations present exciting loci for further experimental characterization.

Evolutionary Origin and Nomenclature of Vertebrate Wnt11-Family Genes

To adequately connect zebrafish medical models to human biology, it is essential that gene nomenclature reflects gene orthology. Analysis of gene phylogenies and conserved syntenies shows that the zebrafish gene currently called wnt11 (ENSDARG00000004256, ZFIN ID: ZDB-GENE-990603-12) is not the ortholog of the human gene called WNT11 (ENSG00000085741); instead, the gene currently called wnt11r (ENSDARG00000014796, ZFIN ID: ZDB-GENE-980526-249) is the zebrafish ortholog of human WNT11. Genomic analysis of Wnt11-family genes suggests a model for the birth of Wnt11-family gene ohnologs in genome duplication events, provides a mechanism for the death of a Wnt11-family ohnolog in mammals after they diverged from birds, and suggests revised nomenclature to better connect teleost disease models to human biology.

Evolution of nodal and nodal-related genes and the putative composition of the heterodimers that trigger the nodal pathway in vertebrates

Nodal is a signaling molecule that belongs to the transforming growth factor-β superfamily that plays key roles during the early stages of development of animals. In vertebrates Nodal forms an heterodimer with a GDF1/3 protein to activate the Nodal pathway. Vertebrates have a paralog of nodal in their genomes labeled Nodal-related, but the evolutionary history of these genes is a matter of debate, mainly because of the presence of a variable numbers of genes in the vertebrate genomes sequenced so far. Thus, the goal of this study was to investigate the evolutionary history of the Nodal and Nodal-related genes with an emphasis in tracking changes in the number of genes among vertebrates. Our results show the presence of two gene lineages (Nodal and Nodal-related) that can be traced back to the ancestor of jawed vertebrates. These lineages have undergone processes of differential retention and lineage-specific expansions. Our results imply that Nodal and Nodal-related duplicated at the latest in the ancestor of gnathostomes, and they still retain a significant level of functional redundancy. By comparing the evolution of the Nodal/Nodal-related with GDF1/3 gene family, it is possible to infer that there are several types of heterodimers that can trigger the Nodal pathway among vertebrates.

TGF-β Family Signaling in Early Vertebrate Development

TGF-β family ligands function in inducing and patterning many tissues of the early vertebrate embryonic body plan. Nodal signaling is essential for the specification of mesendodermal tissues and the concurrent cellular movements of gastrulation. Bone morphogenetic protein (BMP) signaling patterns tissues along the dorsal-ventral axis and simultaneously directs the cell movements of convergence and extension. After gastrulation, a second wave of Nodal signaling breaks the symmetry between the left and right sides of the embryo. During these processes, elaborate regulatory feedback between TGF-β ligands and their antagonists direct the proper specification and patterning of embryonic tissues. In this review, we summarize the current knowledge of the function and regulation of TGF-β family signaling in these processes. Although we cover principles that are involved in the development of all vertebrate embryos, we focus specifically on three popular model organisms: the mouse Mus musculus, the African clawed frog of the genus Xenopus, and the zebrafish Danio rerio, highlighting the similarities and differences between these species.

Madagascar ground gecko genome analysis characterizes asymmetric fates of duplicated genes

Background

Conventionally, comparison among amniotes – birds, mammals, and reptiles – has often been approached through analyses of mammals and, for comparison, birds. However, birds are morphologically and physiologically derived and, moreover, some parts of their genomes are recognized as difficult to sequence and/or assemble and are thus missing in genome assemblies. Therefore, sequencing the genomes of reptiles would aid comparative studies on amniotes by providing more comprehensive coverage to help understand the molecular mechanisms underpinning evolutionary changes.

Results

Herein, we present the whole genome sequences of the Madagascar ground gecko (Paroedura picta), a promising study system especially in developmental biology, and used it to identify changes in gene repertoire across amniotes. The genome-wide analysis of the Madagascar ground gecko allowed us to reconstruct a comprehensive set of gene phylogenies comprising 13,043 ortholog groups from diverse amniotes. Our study revealed 469 genes retained by some reptiles but absent from available genome-wide sequence data of both mammals and birds. Importantly, these genes, herein collectively designated as ‘elusive’ genes, exhibited high nucleotide substitution rates and uneven intra-genomic distribution. Furthermore, the genomic regions flanking these elusive genes exhibited distinct characteristics that tended to be associated with increased gene density, repeat element density, and GC content.

Conclusion

This highly continuous and nearly complete genome assembly of the Madagascar ground gecko will facilitate the use of this species as an experimental animal in diverse fields of biology. Gene repertoire comparisons across amniotes further demonstrated that the fate of a duplicated gene can be affected by the intrinsic properties of its genomic location, which can persist for hundreds of millions of years.

Electronic supplementary material

The online version of this article (10.1186/s12915-018-0509-4) contains supplementary material, which is available to authorized users.

Identification of a fish short-chain dehydrogenase/reductase associated with bone metabolism

Although human and mouse genetics have largely contributed to the better understanding of the mechanisms underlying skeletogenesis, much more remains to be uncovered. In this regard alternative and complementary systems have been sought and cell systems capable of in vitro calcification have been developed to study the mechanisms underlying bone formation. In gilthead seabream (Sparus aurata), a gene coding for an unknown protein that is strongly up-regulated during extracellular matrix (ECM) mineralization of a pre-osteoblast cell line was recently identified as a potentially important player in bone formation. In silico analysis of the deduced protein revealed the presence of domains typical of short-chain dehydrogenase/reductases (SDR). Closely related to carbonyl reductase 1, seabream protein belongs to a novel subfamily of SDR proteins with no orthologs in mammals. Analysis of gene expression by qPCR confirmed the strong up-regulation of sdr-like expression during in vitro mineralization but also revealed high expression levels in calcified tissues. A possible role for Sdr-like in osteoblast and bone metabolism was further evidenced through (i) the localization by in situ hybridization of sdr-like transcript in pre-osteoblasts of the operculum and (ii) the regulation of sdr-like gene transcription by Runx2 and retinoic acid receptor, two regulators of osteoblast differentiation and mineralization. Expression data also indicated a role for Sdr-like in gastrointestinal tract homeostasis and during gilthead seabream development at gastrulation and metamorphosis. This study reports a new subfamily of short-chain dehydrogenases/reductases in vertebrates and, for the first time, provides evidence of a role for SDRs in bone metabolism, osteoblast differentiation and/or tissue mineralization.

Stem Cells, Patterning and Regeneration in Planarians: Self-Organization at the Organismal Scale

The establishment of size and shape remains a fundamental challenge in biological research that planarian flatworms uniquely epitomize. Planarians can regenerate complete and perfectly proportioned animals from tiny and arbitrarily shaped tissue pieces; they continuously renew all organismal cell types from abundant pluripotent stem cells, yet maintain shape and anatomy in the face of constant turnover; they grow when feeding and literally degrow when starving, while scaling form and function over as much as a 40-fold range in body length or an 800-fold change in total cell numbers. This review provides a broad overview of the current understanding of the planarian stem cell system, the mechanisms that pattern the planarian body plan and how the interplay between patterning signals and cell fate choices orchestrates regeneration. What emerges is a conceptual framework for the maintenance and regeneration of the planarian body plan on basis of the interplay between pluripotent stem cells and self-organizing patterns and further, the general utility of planarians as model system for the mechanistic basis of size and shape.

Generation of animal form by the Chordin/Tolloid/BMP gradient: 100 years after D'Arcy Thompson

The classic book "On Growth and Form" by naturalist D'Arcy Thompson was published 100 years ago. To celebrate this landmark, we present experiments in the Xenopus embryo that provide a framework for understanding how simple, quantitative transformations of a morphogen gradient might have affected evolution and morphological diversity of organisms. D'Arcy Thompson proposed that different morphologies might be generated by modifying physical parameters in an underlying system of Cartesian coordinates that pre-existed in Nature and arose during evolutionary history. Chordin is a BMP antagonist secreted by the Spemann organizer located on the dorsal side of the gastrula. Chordin generates a morphogen gradient as first proposed by mathematician Alan Turing. The rate-limiting step of this dorsal-ventral (D-V) morphogen is the degradation of Chordin by the Tolloid metalloproteinase in the ventral side. Chordin is expressed at gastrula on the dorsal side where BMP signaling is low, while at the opposite side peak levels of BMP signaling are reached. In fishes, amphibians, reptiles and birds, high BMP signaling in the ventral region induces transcription of a secreted inhibitor of Tolloid called Sizzled. By depleting Sizzled exclusively in the ventral half of the embryo we were able to expand the ventro-posterior region in an otherwise normal embryo. Conversely, ventral depletion of Tolloid, which stabilizes Chordin, decreased ventral and tail structures, phenocopying the tolloid zebrafish mutation. We explain how historical constraints recorded in the language of DNA become subject to the universal laws of physics when an ancestral reaction-diffusion morphogen gradient dictates form.

Incorporating tree-thinking and evolutionary time scale into developmental biology

Phylogenetic approaches are indispensable in any comparative molecular study involving multiple species. These approaches are in increasing demand as the amount and availability of DNA sequence information continues to increase exponentially, even for organisms that were previously not extensively studied. Without the sound application of phylogenetic concepts and knowledge, one can be misled when attempting to infer ancestral character states as well as the timing and order of evolutionary events, both of which are frequently exerted in evolutionary developmental biology. The ignorance of phylogenetic approaches can also impact non-evolutionary studies and cause misidentification of the target gene or protein to be examined in functional characterization. This review aims to promote tree-thinking in evolutionary conjecture and stress the importance of a sense of time scale in cross-species comparisons, in order to enhance the understanding of phylogenetics in all biological fields including developmental biology. To this end, molecular phylogenies of several developmental regulatory genes, including those denoted as "cryptic pan-vertebrate genes", are introduced as examples.

Unexpected multiplicity of QRFP receptors in early vertebrate evolution

The neuropeptide QRFP, also called 26RFa, and its G protein-coupled receptor GPR103 have been identified in all vertebrates investigated. In mammals, this peptide-receptor pair has been found to have several effects including stimulation of appetite. Recently, we reported that a QRFP peptide is present in amphioxus, Branchiostoma floridae, and we also identified a QRFP receptor (QRFPR) that mediates a functional response to sub-nanomolar concentrations of the amphioxus peptide as well as short and long human QRFP (Xu et al., submitted). Because the ancestral vertebrate underwent two tetraploidizations, it might be expected that duplicates of the QRFP gene and its receptor gene may exist. Indeed, we report here the identification of multiple vertebrate QRFPR genes. Three QRFPR genes are present in the coelacanth Latimeria chalumnae, representing an early diverging sarcopterygian lineage. Three QRFPR genes are present in the basal actinopterygian fish, the spotted gar. Phylogenetic and chromosomal analyses show that only two of these receptor genes are orthologous between the two species, thus demonstrating a total of four distinct vertebrate genes. Three of the QRFPR genes resulted from the early vertebrate tetraploidizations and were copied along with syntenic neuropeptide Y receptor genes. The fourth QRFPR gene may be an even older and distinct lineage. Because mammals and birds have only a single QRFPR gene, this means that three genes have been lost in these lineages, and at least one of these was lost independently in mammals and birds because it is still present in a turtle. In conclusion, these results show that the QRFP system gained considerable complexity in the early stages of vertebrate evolution and still maintains much of this in some lineages, and that it has been secondarily reduced in mammals.

Development: scaling to size by protease inhibition

The dorsal half of bisected Xenopus laevis embryos can regenerate a well-proportioned organism on a smaller scale. A new study indicates that the removal of ventral tissue generates a steeper Chordin gradient by reducing Sizzled, a secreted inhibitor of Tolloid chordinases.

Evolution of mammalian Opn5 as a specialized UV-absorbing pigment by a single amino acid mutation

Opn5 is one of the recently identified opsin groups that is responsible for nonvisual photoreception in animals. We previously showed that a chicken homolog of mammalian Opn5 (Opn5m) is a Gi-coupled UV sensor having molecular properties typical of bistable pigments. Here we demonstrated that mammalian Opn5m evolved to be a more specialized photosensor by losing one of the characteristics of bistable pigments, direct binding of all-trans-retinal. We first confirmed that Opn5m proteins in zebrafish, Xenopus tropicalis, mouse, and human are also UV-sensitive pigments. Then we found that only mammalian Opn5m proteins lack the ability to directly bind all-trans-retinal. Mutational analysis showed that these characteristics were acquired by a single amino acid replacement at position 168. By comparing the expression patterns of Opn5m between mammals and chicken, we found that, like chicken Opn5m, mammalian Opn5m was localized in the ganglion cell layer and inner nuclear layer of the retina. However, the mouse and primate (common marmoset) opsins were distributed not in the posterior hypothalamus (including the region along the third ventricle) where chicken Opn5m is localized, but in the preoptic hypothalamus. Interestingly, RPE65, an essential enzyme for forming 11-cis-retinal in the visual cycle is expressed near the preoptic hypothalamus of the mouse and common marmoset brain but not near the region of the chicken brain where chicken Opn5m is expressed. Therefore, mammalian Opn5m may work exclusively as a short wavelength sensor in the brain as well as in the retina with the assistance of an 11-cis-retinal-supplying system.

Diversity and history of the long-chain acyl-CoA synthetase (Acsl) gene family in vertebrates

BACKGROUND

Fatty acids, a considerable fraction of lipid molecules, participate in fundamental physiological processes. They undergo activation into their corresponding CoA esters for oxidation or esterification into complex lipids (e.g. triglycerides, phospholipids and cholesterol esters), a process that is carried out by acyl-CoA synthases (ACS). Here we analyze the evolution of the gene family encoding for the long-chain acyl-CoA synthetases (Acsl) in vertebrates.

RESULTS

By means of phylogenetics and comparative genomics we show that genome duplications (2R) generated the diversity of Acsl genes in extant vertebrate lineages. In the vertebrate ancestor two separate genes originated the current Acsl1/5/6 and the Acsl3/4 gene families, and the extra gene duplicates in teleosts are a consequence of the teleost specific third round of genome duplication (3R). Moreover, the diversity of Acsl family members is broader than anticipated. Our strategy uncovered a novel uncharacterized Acsl-like gene found in teleosts, spotted gar, coelacanth and possibly lamprey, which we designate Acsl2. The detailed analysis of the Acsl2 teleost gene locus strongly supports the conclusion that it corresponds to a retained 2R paralogue, lost in tetrapods.

CONCLUSIONS

We provide here the first evolutionary analysis of the Acsl gene family in vertebrates, showing the specific contribution of 2R/3R to the diversity of this gene family. We find also that the division of ACSL enzymes into two groups predates at least the emergence of deuterostomes. Our study indicates that genome duplications significantly contributed to the elaboration of fatty acid activation metabolism in vertebrates.

Recurrent gene loss correlates with the evolution of stomach phenotypes in gnathostome history

The stomach, a hallmark of gnathostome evolution, represents a unique anatomical innovation characterized by the presence of acid- and pepsin-secreting glands. However, the occurrence of these glands in gnathostome species is not universal; in the nineteenth century the French zoologist Cuvier first noted that some teleosts lacked a stomach. Strikingly, Holocephali (chimaeras), dipnoids (lungfish) and monotremes (egg-laying mammals) also lack acid secretion and a gastric cellular phenotype. Here, we test the hypothesis that loss of the gastric phenotype is correlated with the loss of key gastric genes. We investigated species from all the main gnathostome lineages and show the specific contribution of gene loss to the widespread distribution of the agastric condition. We establish that the stomach loss correlates with the persistent and complete absence of the gastric function gene kit--H(+)/K(+)-ATPase (Atp4A and Atp4B) and pepsinogens (Pga, Pgc, Cym)--in the analysed species. We also find that in gastric species the pepsinogen gene complement varies significantly (e.g. two to four in teleosts and tens in some mammals) with multiple events of pseudogenization identified in various lineages. We propose that relaxation of purifying selection in pepsinogen genes and possibly proton pump genes in response to dietary changes led to the numerous independent events of stomach loss in gnathostome history. Significantly, the absence of the gastric genes predicts that reinvention of the stomach in agastric lineages would be highly improbable, in line with Dollo's principle.

Evolutionary history of chordate PAX genes: dynamics of change in a complex gene family

Paired box (PAX) genes are transcription factors that play important roles in embryonic development. Although the PAX gene family occurs in animals only, it is widely distributed. Among the vertebrates, its 9 genes appear to be the product of complete duplication of an original set of 4 genes, followed by an additional partial duplication. Although some studies of PAX genes have been conducted, no comprehensive survey of these genes across the entire taxonomic unit has yet been attempted. In this study, we conducted a detailed comparison of PAX sequences from 188 chordates, which revealed restricted variation. The absence of PAX4 and PAX8 among some species of reptiles and birds was notable; however, all 9 genes were present in all 74 mammalian genomes investigated. A search for signatures of selection indicated that all genes are subject to purifying selection, with a possible constraint relaxation in PAX4, PAX7, and PAX8. This result indicates asymmetric evolution of PAX family genes, which can be associated with the emergence of adaptive novelties in the chordate evolutionary trajectory.

A high resolution genome-wide scan for significant selective sweeps: an application to pooled sequence data in laying chickens

In most studies aimed at localizing footprints of past selection, outliers at tails of the empirical distribution of a given test statistic are assumed to reflect locus-specific selective forces. Significance cutoffs are subjectively determined, rather than being related to a clear set of hypotheses. Here, we define an empirical p-value for the summary statistic by means of a permutation method that uses the observed SNP structure in the real data. To illustrate the methodology, we applied our approach to a panel of 2.9 million autosomal SNPs identified from re-sequencing a pool of 15 individuals from a brown egg layer line. We scanned the genome for local reductions in heterozygosity, suggestive of selective sweeps. We also employed a modified sliding window approach that accounts for gaps in the sequence and increases scanning resolution by moving the overlapping windows by steps of one SNP only, and suggest to call this a “creeping window” strategy. The approach confirmed selective sweeps in the region of previously described candidate genes, i.e. TSHR, PRL, PRLHR, INSR, LEPR, IGF1, and NRAMP1 when used as positive controls. The genome scan revealed 82 distinct regions with strong evidence of selection (genome-wide p-value<0.001), including genes known to be associated with eggshell structure and immune system such as CALB1 and GAL cluster, respectively. A substantial proportion of signals was found in poor gene content regions including the most extreme signal on chromosome 1. The observation of multiple signals in a highly selected layer line of chicken is consistent with the hypothesis that egg production is a complex trait controlled by many genes.

GLADX: an automated approach to analyze the lineage-specific loss and pseudogenization of genes

A well-established ancestral gene can usually be found, in one or multiple copies, in different descendant species. Sometimes during the course of evolution, all the representatives of a well-established ancestral gene disappear in specific lineages; such gene losses may occur in the genome by deletion of a DNA fragment or by pseudogenization. The loss of an entire gene family in a given lineage may reflect an important phenomenon, and could be due either to adaptation, or to a relaxation of selection that leads to neutral evolution. Therefore, the lineage-specific gene loss analyses are important to improve the understanding of the evolutionary history of genes and genomes. In order to perform this kind of study from the increasing number of complete genome sequences available, we developed a unique new software module called GLADX in the DAGOBAH framework, based on a comparative genomic approach. The software is able to automatically detect, for all the species of a phylum, the presence/absence of a representative of a well-established ancestral gene, and by systematic steps of re-annotation, confirm losses, detect and analyze pseudogenes and find novel genes. The approach is based on the use of highly reliable gene phylogenies, of protein predictions and on the analysis of genomic mutations. All the evidence associated to evolutionary approach provides accurate information for building an overall view of the evolution of a given gene in a selected phylum. The reliability of GLADX has been successfully tested on a benchmark analysis of 14 reported cases. It is the first tool that is able to fully automatically study the lineage-specific losses and pseudogenizations. GLADX is available at http://ioda.univ-provence.fr/IodaSite/gladx/.

The evolution of pepsinogen C genes in vertebrates: duplication, loss and functional diversification

Background

Aspartic proteases comprise a large group of enzymes involved in peptide proteolysis. This collection includes prominent enzymes globally categorized as pepsins, which are derived from pepsinogen precursors. Pepsins are involved in gastric digestion, a hallmark of vertebrate physiology. An important member among the pepsinogens is pepsinogen C (Pgc). A particular aspect of Pgc is its apparent single copy status, which contrasts with the numerous gene copies found for example in pepsinogen A (Pga). Although gene sequences with similarity to Pgc have been described in some vertebrate groups, no exhaustive evolutionary framework has been considered so far.

Methodology/Principal Findings

By combining phylogenetics and genomic analysis, we find an unexpected Pgc diversity in the vertebrate sub-phylum. We were able to reconstruct gene duplication timings relative to the divergence of major vertebrate clades. Before tetrapod divergence, a single Pgc gene tandemly expanded to produce two gene lineages (Pgbc and Pgc2). These have been differentially retained in various classes. Accordingly, we find Pgc2 in sauropsids, amphibians and marsupials, but not in eutherian mammals. Pgbc was retained in amphibians, but duplicated in the ancestor of amniotes giving rise to Pgb and Pgc1. The latter was retained in mammals and probably in reptiles and marsupials but not in birds. Pgb was kept in all of the amniote clade with independent episodes of loss in some mammalian species. Lineage specific expansions of Pgc2 and Pgbc have also occurred in marsupials and amphibians respectively. We find that teleost and tetrapod Pgc genes reside in distinct genomic regions hinting at a possible translocation.

Conclusions

We conclude that the repertoire of Pgc genes is larger than previously reported, and that tandem duplications have modelled the history of Pgc genes. We hypothesize that gene expansion lead to functional divergence in tetrapods, coincident with the invasion of terrestrial habitats.