MEPE evolution in mammals reveals regions and residues of prime functional importance

A synthetic review: natural history of amniote reproductive modes in light of comparative evolutionary genomics

There is a current lack of consensus on whether the ancestral parity mode was oviparity (egg-laying) or viviparity (live-birth) in amniotes and particularly in squamates (snakes, lizards, and amphisbaenids). How transitions between parity modes occur at the genomic level has primary importance for how science conceptualises the origin of amniotes, and highly variable parity modes in Squamata. Synthesising literature from medicine, poultry science, reproductive biology, and evolutionary biology, I review the genomics and physiology of five broad processes (here termed the 'Main Five') expected to change during transitions between parity modes: eggshell formation, embryonic retention, placentation, calcium transport, and maternal-fetal immune dynamics. Throughout, I offer alternative perspectives and testable hypotheses regarding proximate causes of parity mode evolution in amniotes and squamates. If viviparity did evolve early in the history of lepidosaurs, I offer the nucleation site hypothesis as a proximate explanation. The framework of this hypothesis can be extended to amniotes to infer their ancestral state. I also provide a mechanism and hypothesis on how squamates may transition from viviparity to oviparity and make predictions about the directionality of transitions in three species. After considering evidence for differing perspectives on amniote origins, I offer a framework that unifies (i) the extended embryonic retention model and (ii) the traditional model which describes the amniote egg as an adaptation to the terrestrial environment. Additionally, this review contextualises the origin of amniotes and parity mode evolution within Medawar's paradigm. Medawar posited that pregnancy could be supported by immunosuppression, inertness, evasion, or immunological barriers. I demonstrate that this does not support gestation or gravidity across most amniotes but may be an adequate paradigm to explain how the first amniote tolerated internal fertilization and delayed egg deposition. In this context, the eggshell can be thought of as an immunological barrier. If serving as a barrier underpins the origin of the amniote eggshell, there should be evidence that oviparous gravidity can be met with a lack of immunological responses in utero. Rare examples of two species that differentially express very few genes during gravidity, suggestive of an absent immunological reaction to oviparous gravidity, are two skinks Lampropholis guichenoti and Lerista bougainvillii. These species may serve as good models for the original amniote egg. Overall, this review grounds itself in the historical literature while offering a modern perspective on the origin of amniotes. I encourage the scientific community to utilise this review as a resource in evolutionary and comparative genomics studies, embrace the complexity of the system, and thoughtfully consider the frameworks proposed.

Evolution of the Avian Eggshell Biomineralization Protein Toolkit - New Insights From Multi-Omics

The avian eggshell is a remarkable biomineral, which is essential for avian reproduction; its properties permit embryonic development in the desiccating terrestrial environment, and moreover, are critically important to preserve unfertilized egg quality for human consumption. This calcium carbonate (CaCO₃) bioceramic is made of 95% calcite and 3.5% organic matrix; it protects the egg contents against microbial penetration and mechanical damage, allows gaseous exchange, and provides calcium for development of the embryonic skeleton. In vertebrates, eggshell occurs in the Sauropsida and in a lesser extent in Mammalia taxa; avian eggshell calcification is one of the fastest known CaCO₃ biomineralization processes, and results in a material with excellent mechanical properties. Thus, its study has triggered a strong interest from the researcher community. The investigation of eggshell biomineralization in birds over the past decades has led to detailed characterization of its protein and mineral constituents. Recently, our understanding of this process has been significantly improved using high-throughput technologies (i.e., proteomics, transcriptomics, genomics, and bioinformatics). Presently, more or less complete eggshell proteomes are available for nine birds, and therefore, key proteins that comprise the eggshell biomineralization toolkit are beginning to be identified. In this article, we review current knowledge on organic matrix components from calcified eggshell. We use these data to analyze the evolution of selected matrix proteins and underline their role in the biological toolkit required for eggshell calcification in avian species. Amongst the panel of eggshell-associated proteins, key functional domains are present such as calcium-binding, vesicle-binding and protein-binding. These technical advances, combined with progress in mineral ultrastructure analyses, have opened the way for new hypotheses of mineral nucleation and crystal growth in formation of the avian eggshell, including transfer of amorphous CaCO₃ in vesicles from uterine cells to the eggshell mineralization site. The enrichment of multi-omics datasets for bird species is critical to understand the evolutionary context for development of CaCO₃ biomineralization in metazoans, leading to the acquisition of the robust eggshell in birds (and formerly dinosaurs).

Biological stenciling of mineralization in the skeleton: Local enzymatic removal of inhibitors in the extracellular matrix

Biomineralization is remarkably diverse and provides myriad functions across many organismal systems. Biomineralization processes typically produce hardened, hierarchically organized structures usually having nanostructured mineral assemblies that are formed through inorganic-organic (usually protein) interactions. Calcium‑carbonate biomineral predominates in structures of small invertebrate organisms abundant in marine environments, particularly in shells (remarkably it is also found in the inner ear otoconia of vertebrates), whereas calcium-phosphate biomineral predominates in the skeletons and dentitions of both marine and terrestrial vertebrates, including humans. Reconciliation of the interplay between organic moieties and inorganic crystals in bones and teeth is a cornerstone of biomineralization research. Key molecular determinants of skeletal and dental mineralization have been identified in health and disease, and in pathologic ectopic calcification, ranging from small molecules such as pyrophosphate, to small membrane-bounded matrix vesicles shed from cells, and to noncollagenous extracellular matrix proteins such as osteopontin and their derived bioactive peptides. Beyond partly knowing the regulatory role of the direct actions of inhibitors on vertebrate mineralization, more recently the importance of their enzymatic removal from the extracellular matrix has become increasingly understood. Great progress has been made in deciphering the relationship between mineralization inhibitors and the enzymes that degrade them, and how adverse changes in this physiologic pathway (such as gene mutations causing disease) result in mineralization defects. Two examples of this are rare skeletal diseases having osteomalacia/odontomalacia (soft bones and teeth) - namely hypophosphatasia (HPP) and X-linked hypophosphatemia (XLH) - where inactivating mutations occur in the gene for the enzymes tissue-nonspecific alkaline phosphatase (TNAP, TNSALP, ALPL) and phosphate-regulating endopeptidase homolog X-linked (PHEX), respectively. Here, we review and provide a concept for how existing and new information now comes together to describe the dual nature of regulation of mineralization - through systemic mineral ion homeostasis involving circulating factors, coupled with molecular determinants operating at the local level in the extracellular matrix. For the local mineralization events in the extracellular matrix, we present a focused concept in skeletal mineralization biology called the Stenciling Principle - a principle (building upon seminal work by Neuman and Fleisch) describing how the action of enzymes to remove tissue-resident inhibitors defines with precision the location and progression of mineralization.

Evolutionary Analysis Predicts Sensitive Positions of MMP20 and Validates Newly- and Previously-Identified MMP20 Mutations Causing Amelogenesis Imperfecta

Amelogenesis imperfecta (AI) designates a group of genetic diseases characterized by a large range of enamel disorders causing important social and health problems. These defects can result from mutations in enamel matrix proteins or protease encoding genes. A range of mutations in the enamel cleavage enzyme matrix metalloproteinase-20 gene (MMP20) produce enamel defects of varying severity. To address how various alterations produce a range of AI phenotypes, we performed a targeted analysis to find MMP20 mutations in French patients diagnosed with non-syndromic AI. Genomic DNA was isolated from saliva and MMP20 exons and exon-intron boundaries sequenced. We identified several homozygous or heterozygous mutations, putatively involved in the AI phenotypes. To validate missense mutations and predict sensitive positions in the MMP20 sequence, we evolutionarily compared 75 sequences extracted from the public databases using the Datamonkey webserver. These sequences were representative of mammalian lineages, covering more than 150 million years of evolution. This analysis allowed us to find 324 sensitive positions (out of the 483 MMP20 residues), pinpoint functionally important domains, and build an evolutionary chart of important conserved MMP20 regions. This is an efficient tool to identify new- and previously-identified mutations. We thus identified six functional MMP20 mutations in unrelated families, finding two novel mutated sites. The genotypes and phenotypes of these six mutations are described and compared. To date, 13 MMP20 mutations causing AI have been reported, making these genotypes and associated hypomature enamel phenotypes the most frequent in AI.

Habitat variability does not generally promote metabolic network modularity in flies and mammals

The evolution of species habitat range is an important topic over a wide range of research fields. In higher organisms, habitat range evolution is generally associated with genetic events such as gene duplication. However, the specific factors that determine habitat variability remain unclear at higher levels of biological organization (e.g., biochemical networks). One widely accepted hypothesis developed from both theoretical and empirical analyses is that habitat variability promotes network modularity; however, this relationship has not yet been directly tested in higher organisms. Therefore, I investigated the relationship between habitat variability and metabolic network modularity using compound and enzymatic networks in flies and mammals. Contrary to expectation, there was no clear positive correlation between habitat variability and network modularity. As an exception, the network modularity increased with habitat variability in the enzymatic networks of flies. However, the observed association was likely an artifact, and the frequency of gene duplication appears to be the main factor contributing to network modularity. These findings raise the question of whether or not there is a general mechanism for habitat range expansion at a higher level (i.e., above the gene scale). This study suggests that the currently widely accepted hypothesis for habitat variability should be reconsidered.

Evolutionary analysis of selective constraints identifies ameloblastin (AMBN) as a potential candidate for amelogenesis imperfecta

Background

Ameloblastin (AMBN) is a phosphorylated, proline/glutamine-rich protein secreted during enamel formation. Previous studies have revealed that this enamel matrix protein was present early in vertebrate evolution and certainly plays important roles during enamel formation although its precise functions remain unclear. We performed evolutionary analyses of AMBN in order to (i) identify residues and motifs important for the protein function, (ii) predict mutations responsible for genetic diseases, and (iii) understand its molecular evolution in mammals.

Results

In silico searches retrieved 56 complete sequences in public databases that were aligned and analyzed computationally. We showed that AMBN is globally evolving under moderate purifying selection in mammals and contains a strong phylogenetic signal. In addition, our analyses revealed codons evolving under significant positive selection. Evidence for positive selection acting on AMBN was observed in catarrhine primates and the aye-aye. We also found that (i) an additional translation initiation site was recruited in the ancestral placental AMBN, (ii) a short exon was duplicated several times in various species including catarrhine primates, and (iii) several polyadenylation sites are present.

Conclusions

AMBN possesses many positions, which have been subjected to strong selective pressure for 200 million years. These positions correspond to several cleavage sites and hydroxylated, O-glycosylated, and phosphorylated residues. We predict that these conserved positions would be potentially responsible for enamel disorder if substituted. Some motifs that were previously identified as potentially important functionally were confirmed, and we found two, highly conserved, new motifs, the function of which should be tested in the near future. This study illustrates the power of evolutionary analyses for characterizing the functional constraints acting on proteins with yet uncharacterized structure.

Electronic supplementary material

The online version of this article (doi:10.1186/s12862-015-0431-0) contains supplementary material, which is available to authorized users.

Molecular evolution of the tissue-nonspecific alkaline phosphatase allows prediction and validation of missense mutations responsible for hypophosphatasia

ALPL encodes the tissue nonspecific alkaline phosphatase (TNSALP), which removes phosphate groups from various substrates. Its function is essential for bone and tooth mineralization. In humans, ALPL mutations lead to hypophosphatasia, a genetic disorder characterized by defective bone and/or tooth mineralization. To date, 275 ALPL mutations have been reported to cause hypophosphatasia, of which 204 were simple missense mutations. Molecular evolutionary analysis has proved to be an efficient method to highlight residues important for the protein function and to predict or validate sensitive positions for genetic disease. Here we analyzed 58 mammalian TNSALP to identify amino acids unchanged, or only substituted by residues sharing similar properties, through 220 millions years of mammalian evolution. We found 469 sensitive positions of the 524 residues of human TNSALP, which indicates a highly constrained protein. Any substitution occurring at one of these positions is predicted to lead to hypophosphatasia. We tested the 204 missense mutations resulting in hypophosphatasia against our predictive chart, and validated 99% of them. Most sensitive positions were located in functionally important regions of TNSALP (active site, homodimeric interface, crown domain, calcium site, …). However, some important positions are located in regions, the structure and/or biological function of which are still unknown. Our chart of sensitive positions in human TNSALP (i) enables to validate or invalidate at low cost any ALPL mutation, which would be suspected to be responsible for hypophosphatasia, by contrast with time consuming and expensive functional tests, and (ii) displays higher predictive power than in silico models of prediction.

Contribution of nonohnologous duplicated genes to high habitat variability in mammals

The mechanism by which genetic systems affect environmental adaptation is a focus of considerable attention in the fields of ecology, evolution, and conservation. However, the genomic characteristics that constrain adaptive evolution have remained unknown. A recent study showed that the proportion of duplicated genes in whole Drosophila genomes correlated with environmental variability within habitat, but it remains unclear whether the correlation is observed even in vertebrates whose genomes including a large number of duplicated genes generated by whole-genome duplication (WGD). Here, we focus on fully sequenced mammalian genomes that experienced WGD in early vertebrate lineages and show that the proportion of small-scale duplication (SSD) genes in the genome, but not that of WGD genes, is significantly correlated with habitat variability. Moreover, species with low habitat variability have a higher proportion of lost duplicated genes, particularly SSD genes, than those with high habitat variability. These results indicate that species that inhabit variable environments may maintain more SSD genes in their genomes and suggest that SSD genes are important for adapting to novel environments and surviving environmental changes. These insights may be applied to predicting invasive and endangered species.

Hen uterine gene expression profiling during eggshell formation reveals putative proteins involved in the supply of minerals or in the shell mineralization process

BACKGROUND

The chicken eggshell is a natural mechanical barrier to protect egg components from physical damage and microbial penetration. Its integrity and strength is critical for the development of the embryo or to ensure for consumers a table egg free of pathogens. This study compared global gene expression in laying hen uterus in the presence or absence of shell calcification in order to characterize gene products involved in the supply of minerals and / or the shell biomineralization process.

RESULTS

Microarrays were used to identify a repertoire of 302 over-expressed genes during shell calcification. GO terms enrichment was performed to provide a global interpretation of the functions of the over-expressed genes, and revealed that the most over-represented proteins are related to reproductive functions. Our analysis identified 16 gene products encoding proteins involved in mineral supply, and allowed updating of the general model describing uterine ion transporters during eggshell calcification. A list of 57 proteins potentially secreted into the uterine fluid to be active in the mineralization process was also established. They were classified according to their potential functions (biomineralization, proteoglycans, molecular chaperone, antimicrobials and proteases/antiproteases).

CONCLUSIONS

Our study provides detailed descriptions of genes and corresponding proteins over-expressed when the shell is mineralizing. Some of these proteins involved in the supply of minerals and influencing the shell fabric to protect the egg contents are potentially useful biological markers for the genetic improvement of eggshell quality.

Extracellular matrix mineralization in periodontal tissues: Noncollagenous matrix proteins, enzymes, and relationship to hypophosphatasia and X-linked hypophosphatemia

As broadly demonstrated for the formation of a functional skeleton, proper mineralization of periodontal alveolar bone and teeth - where calcium phosphate crystals are deposited and grow within an extracellular matrix - is essential for dental function. Mineralization defects in tooth dentin and cementum of the periodontium invariably lead to a weak (soft or brittle) dentition in which teeth become loose and prone to infection and are lost prematurely. Mineralization of the extremities of periodontal ligament fibers (Sharpey's fibers) where they insert into tooth cementum and alveolar bone is also essential for the function of the tooth-suspensory apparatus in occlusion and mastication. Molecular determinants of mineralization in these tissues include mineral ion concentrations (phosphate and calcium), pyrophosphate, small integrin-binding ligand N-linked glycoproteins and matrix vesicles. Amongst the enzymes important in regulating these mineralization determinants, two are discussed at length here, with clinical examples given, namely tissue-nonspecific alkaline phosphatase and phosphate-regulating gene with homologies to endopeptidases on the X chromosome. Inactivating mutations in these enzymes in humans and in mouse models lead to the soft bones and teeth characteristic of hypophosphatasia and X-linked hypophosphatemia, respectively, where the levels of local and systemic circulating mineralization determinants are perturbed. In X-linked hypophosphatemia, in addition to renal phosphate wasting causing low circulating phosphate levels, phosphorylated mineralization-regulating small integrin-binding ligand N-linked glycoproteins, such as matrix extracellular phosphoglycoprotein and osteopontin, and the phosphorylated peptides proteolytically released from them, such as the acidic serine- and aspartate-rich-motif peptide, may accumulate locally to impair mineralization in this disease.

The dentin matrix acidic phosphoprotein 1 (DMP1) in the light of mammalian evolution

Dentin matrix acidic phosphoprotein 1 (DMP1) is an acidic, highly phosphorylated, noncollagenous protein secreted during dentin and bone formation. Previous functional studies of DMP1 have revealed various motifs playing a role in either mineralization or cell differentiation. We performed an evolutionary analysis of DMP1 to identify residues and motifs that were conserved during 220 millions years (Ma) of mammalian evolution, and hence have an important function. In silico search provided us with 41 sequences that were aligned and analyzed using the Hyphy program. We showed that DMP1 contains 55 positions that were kept unchanged for 220 Ma. We also defined in a more precise manner some motifs that were already known (i.e., cleavage sites, RGD motif, ASARM peptide, glycosaminoglycan chain attachment site, nuclear localization signal sites, and dentin sialophosphoprotein-binding site), and we found five, highly conserved, new functional motifs. In the near future, functional studies could be performed to understand the role played by them.

Evolutionary conservation of the oocyte transcriptome among vertebrates and its implications for understanding human reproductive function

Cross-phylum and cross-species comparative transcriptomic analyses provide an evolutionary perspective on how specific tissues use genomic information. A significant mRNA subset present in the oocytes of most vertebrates is stabilized or stored for post-LH surge use. Since transcription is arrested in the oocyte before ovulation, this RNA is important for completing maturation and sustaining embryo development until zygotic genome activation. We compared the human oocyte transcriptome with an oocyte-enriched subset of mouse, bovine and frog (Xenopus laevis) genes in order to evaluate similarities between species. Graded temperature stringency hybridization on a multi-species oocyte cDNA array was used to measure the similarity of preferentially expressed sequences to the human oocyte library. Identity analysis of 679 human orthologs compared with each identified official gene symbol found in the subtractive (somatic-oocyte) libraries comprising our array revealed that bovine/human similarity was greater than mouse/human or frog/human similarity. However, based on protein sequence, mouse/human similarity was greater than bovine/human similarity. Among the genes over-expressed in oocytes relative to somatic tissue in Xenopus, Mus and Bos, a high level of conservation was found relative to humans, especially for genes involved in early embryonic development.

The chicken or the egg: PHEX, FGF23 and SIBLINGs unscrambled

The eggshell is an ancient innovation that helped the vertebrates' transition from the oceans and gain dominion over the land. Coincident with this conquest, several new eggshell and noncollagenous bone-matrix proteins (NCPs) emerged. The protein ovocleidin-116 is one of these proteins with an ancestry stretching back to the Triassic. Ovocleidin-116 is an avian homolog of Matrix Extracellular Phosphoglycoprotein (MEPE) and belongs to a group of proteins called Small Integrin-Binding Ligand Interacting Glycoproteins (SIBLINGs). The genes for these NCPs are all clustered on chromosome 5q in mice and chromosome 4q in humans. A unifying feature of the SIBLING proteins is an Acidic Serine Aspartate-Rich MEPE (ASARM)-associated motif. The ASARM motif and the released ASARM peptide play roles in mineralization, bone turnover, mechanotransduction, phosphate regulation and energy metabolism. ASARM peptides and motifs are physiological substrates for phosphate-regulating gene with homologies to endopeptidases on the X chromosome (PHEX), a Zn metalloendopeptidase. Defects in PHEX are responsible for X-linked hypophosphatemic rickets. PHEX interacts with another ASARM motif containing SIBLING protein, Dentin Matrix Protein-1 (DMP1). DMP1 mutations cause bone-renal defects that are identical with the defects caused by loss of PHEX function. This results in autosomal recessive hypophosphatemic rickets (ARHR). In both X-linked hypophosphatemic rickets and ARHR, increased fibroblast growth factor 23 (FGF23) expression occurs, and activating mutations in FGF23 cause autosomal dominant hypophosphatemic rickets (ADHR). ASARM peptide administration in vitro and in vivo also induces increased FGF23 expression. This review will discuss the evidence for a new integrative pathway involved in bone formation, bone-renal mineralization, renal phosphate homeostasis and energy metabolism in disease and health.

Regulation of bone-renal mineral and energy metabolism: the PHEX, FGF23, DMP1, MEPE ASARM pathway

More than 300 million years ago, vertebrates emerged from the vast oceans to conquer gravity and the dry land. With this transition, new adaptations occurred that included ingenious changes in reproduction, waste secretion, and bone physiology. One new innovation, the egg shell, contained an ancestral protein (ovocleidin-116) that likely first appeared with the dinosaurs and was preserved through the theropod lineage in modern birds and reptiles. Ovocleidin-116 is an avian homolog of matrix extracellular phosphoglycoprotein (MEPE) and belongs to a group of proteins called short integrin-binding ligand-interacting glycoproteins (SIBLINGs). These proteins are all localized to a defined region on chromosome 5q in mice and chromosome 4q in humans. A unifying feature of SIBLING proteins is an acidic serine aspartate-rich MEPE-associated motif (ASARM). Recent research has shown that the ASARM motif and the released ASARM peptide have regulatory roles in mineralization (bone and teeth), phosphate regulation, vascularization, soft-tissue calcification, osteoclastogenesis, mechanotransduction, and fat energy metabolism. The MEPE ASARM motif and peptide are physiological substrates for PHEX, a zinc metalloendopeptidase. Defects in PHEX are responsible for X-linked hypophosphatemic rickets (HYP). There is evidence that PHEX interacts with another ASARM motif containing SIBLING protein, dentin matrix protein-1 (DMP1). DMP1 mutations cause bone and renal defects that are identical with the defects caused by a loss of PHEX function. This results in autosomal recessive hypophosphatemic rickets (ARHR). In both HYP and ARHR, increased FGF23 expression plays a major role in the disease and in autosomal dominant hypophosphatemic rickets (ADHR), FGF23 half-life is increased by activating mutations. ASARM peptide administration in vitro and in vivo also induces increased FGF23 expression. FGF23 is a member of the fibroblast growth factor (FGF) family of cytokines, which surfaced 500 million years ago with the boney fish (i.e., teleosts) that do not contain SIBLING proteins. In terrestrial vertebrates, FGF23, like SIBLING proteins, is expressed in the osteocyte. The boney fish, however, are an-osteocytic, so a physiological bone-renal link with FGF23 and the SIBLINGs was cemented when life ventured from the oceans to the land during the Triassic period, approximately 300 million years ago. This link has been revealed by recent research that indicates a competitive displacement of a PHEX-DMP1 interaction by an ASARM peptide that leads to increased FGF23 expression. This review discusses the new discoveries that reveal a novel PHEX, DMP1, MEPE, ASARM peptide, and FGF23 bone-renal pathway. This pathway impacts not only bone formation, bone-renal mineralization, and renal phosphate homeostasis but also energy metabolism. The study of this new pathway is relevant for developing therapies for several diseases: bone-teeth mineral loss disorders, renal osteodystrophy, chronic kidney disease and bone mineralization disorders (CKD-MBD), end-stage renal diseases, ectopic arterial-calcification, cardiovascular disease renal calcification, diabetes, and obesity.

Adaptive evolution of the matrix extracellular phosphoglycoprotein in mammals

Background

Matrix extracellular phosphoglycoprotein (MEPE) belongs to a family of small integrin-binding ligand N-linked glycoproteins (SIBLINGs) that play a key role in skeleton development, particularly in mineralization, phosphate regulation and osteogenesis. MEPE associated disorders cause various physiological effects, such as loss of bone mass, tumors and disruption of renal function (hypophosphatemia). The study of this developmental gene from an evolutionary perspective could provide valuable insights on the adaptive diversification of morphological phenotypes in vertebrates.

Results

Here we studied the adaptive evolution of the MEPE gene in 26 Eutherian mammals and three birds. The comparative genomic analyses revealed a high degree of evolutionary conservation of some coding and non-coding regions of the MEPE gene across mammals indicating a possible regulatory or functional role likely related with mineralization and/or phosphate regulation. However, the majority of the coding region had a fast evolutionary rate, particularly within the largest exon (1467 bp). Rodentia and Scandentia had distinct substitution rates with an increased accumulation of both synonymous and non-synonymous mutations compared with other mammalian lineages. Characteristics of the gene (e.g. biochemical, evolutionary rate, and intronic conservation) differed greatly among lineages of the eight mammalian orders. We identified 20 sites with significant positive selection signatures (codon and protein level) outside the main regulatory motifs (dentonin and ASARM) suggestive of an adaptive role. Conversely, we find three sites under selection in the signal peptide and one in the ASARM motif that were supported by at least one selection model. The MEPE protein tends to accumulate amino acids promoting disorder and potential phosphorylation targets.

Conclusion

MEPE shows a high number of selection signatures, revealing the crucial role of positive selection in the evolution of this SIBLING member. The selection signatures were found mainly outside the functional motifs, reinforcing the idea that other regions outside the dentonin and the ASARM might be crucial for the function of the protein and future studies should be undertaken to understand its importance.

OC-116, the chicken ortholog of mammalian MEPE found in eggshell, is also expressed in bone cells

In chicken, ovocleidin 116 (OC-116) is found in the eggshell matrix and its encoding gene, OC-116, is expressed in uterine cells. In mammals, its orthologue MEPE encodes the matrix extracellular phosphoglycoprotein (MEPE), which has been shown to be involved in bone mineralization. Using RT-PCR and in situ hybridization on sections, we have checked whether OC-116 was also expressed in osteoblasts and osteocytes during bone development and mineralization in chicken embryos. We monitored OC-116 expression in the tibia and mandible of a growth series of chicken embryos from E3 to E19. Transcripts were identified in the osteoblasts as early as E5 in the tibia and E7 in the mandible, before matrix mineralization, then from these stages onwards in both the osteoblasts lining the mineralized bone matrix and the osteocytes. Therefore, early in chicken ontogeny and as soon as osteogenesis begins, OC-116 is involved. Its function, which remains still unknown, is maintained during further bone growth and mineralization, and later in adult, in which it is recruited for eggshell formation. We hypothesize that the ancestral OC-116/MEPE in a stem amniote was involved in these two functions and that the loss of eggshell in the mammalian lineage has probably favored the recruitment of some MEPE domains toward new functions in osteogenesis and mineralization, and in phosphatemia regulation.

Lipid-binding surfaces of membrane proteins: evidence from evolutionary and structural analysis

Membrane proteins function in the diverse environment of the lipid bilayer. Experimental evidence suggests that some lipid molecules bind tightly to specific sites on the membrane protein surface. These lipid molecules often act as co-factors and play important functional roles. In this study, we have assessed the evolutionary selection pressure experienced at lipid-binding sites in a set of α-helical and β-barrel membrane proteins using posterior probability analysis of the ratio of synonymous vs. nonsynonymous substitutions (ω-ratio). We have also carried out a geometric analysis of the membrane protein structures to identify residues in close contact with co-crystallized lipids. We found that residues forming cholesterol-binding sites in both β(2)-adrenergic receptor and Na(+)-K(+)-ATPase exhibit strong conservation, which can be characterized by an expanded cholesterol consensus motif for GPCRs. Our results suggest the functional importance of aromatic stacking interactions and interhelical hydrogen bonds in facilitating protein-cholesterol interactions, which is now reflected in the expanded motif. We also find that residues forming the cardiolipin-binding site in formate dehydrogenase-N γ-subunit and the phosphatidylglycerol binding site in KcsA are under strong purifying selection pressure. Although the lipopolysaccharide (LPS)-binding site in ferric hydroxamate uptake receptor (FhuA) is only weakly conserved, we show using a statistical mechanical model that LPS binds to the least stable FhuA β-strand and protects it from the bulk lipid. Our results suggest that specific lipid binding may be a general mechanism employed by β-barrel membrane proteins to stabilize weakly stable regions. Overall, we find that the residues forming specific lipid binding sites on the surfaces of membrane proteins often experience strong purifying selection pressure.

What makes an egg unique? Clues from evolutionary scenarios of egg-specific genes

The avian egg, which contains the egg yolk, the egg white, and the eggshell, represents the mostly advanced amniotic egg in oviparous vertebrates. In mammals, this reproductive strategy of laying egg has gradually evolved toward placentation. In order to better understand the unique status of the avian egg in the evolution of the vertebrate reproduction, we investigated the evolution of some Gallus gallus egg-specific protein-coding genes. Based on our finding and other recent research, we have summarized here that gene formation (such as ovalbumin genes, ovocalyxin-36 and apovitellenin-1 encoding genes in the G. gallus), gene divergence between G. gallus and mammals (such as the ovocalyxin-32 gene with its ortholog, the mammalian RARRES1, and the ovocleidin-116 with its ortholog, the mammalian MEPE), and gene loss (egg-expressed genes lost during the evolution of the mammals, such as vitellogenin and RBP encoding genes) play significant roles in the evolution of egg-specific genes.