Mónica Lopes-Marques

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

BackgroundFatty 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.ResultsBy 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.ConclusionsWe 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.

Evolutionary functional elaboration of the Elovl2/5 gene family in chordates

The biosynthesis of long-chain polyunsaturated fatty acids (LC-PUFA) provides an intriguing example on how multi-enzymatic cascades evolve. Essential LC-PUFA, such as arachidonic, eicosapentaenoic, and docosahexaenoic acids (DHA), can be acquired from the diet but are also endogenously retailored from C18 precursors through consecutive elongations and desaturations catalyzed, respectively, by fatty acyl elongase and desaturase enzymes. The molecular wiring of this enzymatic pathway defines the ability of a species to biosynthesize LC-PUFA. Exactly when and how in animal evolution a functional LC-PUFA pathway emerged is still elusive. Here we examine key components of the LC-PUFA cascade, the Elovl2/Elovl5 elongases, from amphioxus, an invertebrate chordate, the sea lamprey, a representative of agnathans, and the elephant shark, a basal jawed vertebrate. We show that Elovl2 and Elovl5 emerged from genome duplications in vertebrate ancestry. The single Elovl2/5 from amphioxus efficiently elongates C18 and C20 and, to a marked lesser extent, C22 LC-PUFA. Lamprey is incapable of elongating C22 substrates. The elephant shark Elovl2 showed that the ability to efficiently elongate C22 PUFA and thus to synthesize DHA through the Sprecher pathway, emerged in the jawed vertebrate ancestor. Our findings illustrate how non-integrated “metabolic islands” evolve into fully wired pathways upon duplication and neofunctionalization.

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.

Statins: An undesirable class of aquatic contaminants?

Emerging pollutants, such as pharmaceuticals, may pose a considerable environment risk. Hypocholesterolaemic drugs such as statins are among the most prescribed human pharmaceuticals in western European countries. In vertebrates, this therapeutic class disrupts the cholesterol synthesis by inhibiting the enzyme 3-hydroxy-3-methyl-glutaryl-CoA reductase (HMGR), responsible for the limiting step in the mevalonate pathway. Recently, functional studies have shown that statins competitively inhibit HMGR in vertebrates and arthropods, two taxa that have diverged over 450 million years ago. Importantly, chronic simvastatin exposure disrupts crustacean reproduction and development at environmentally relevant concentrations. Hence, a fundamental question emerges: what is the taxonomic scope of statins-induced HMGR inhibition across metazoans? Here, we address this central question in a large sampling of metazoans using comparative genomics, homology modelling and molecular docking. Sequence alignment of metazoan HMGRs allowed the annotation of highly conserved catalytic, co-factor and substrate binding sites, including residues highjacked for statin binding. Furthermore, molecular docking shows that the catalytic domains of metazoan HMGRs are highly conserved regarding interactions, not only with HMG-CoA, but also with both simvastatin and atorvastatin, the top prescribed statins in Europe and USA. Hence, the data indicates that both statins are expected to competitively inhibit metazoans HMGRs, and therefore all metazoan taxa might be at risk. The environmental relevance of these findings are discussed and research priorities established. We believe that the conceptual framework used in this study can be applied to other emerging pollutants and assist in the design of toxicity testing and risk assessment.

Molecular and functional characterization of a fads2 orthologue in the Amazonian teleost, Arapaima gigas

The Brazilian teleost Arapaima gigas is an iconic species of the Amazon. In recent years a significant effort has been put into the farming of arapaima to mitigate overfishing threats. However, little is known regarding the nutritional requirements of A. gigas in particular those for essential fatty acids including the long-chain polyunsaturated fatty acids (LC-PUFA) eicosapentaenoic acid (EPA) and docosahexaenoic acid (DHA). The ability to biosynthesize LC-PUFA is dependent upon the gene repertoire of fatty acyl desaturases (Fads) and elongases (Elovl), as well as their fatty acid specificities. In the present study we characterized both molecularly and functionally an orthologue of the desaturase fatty acid desaturase 2 (fads2) from A. gigas. The isolated sequence displayed the typical desaturase features, a cytochrome b5-domain with the heme-binding motif, two transmembrane domains and three histidine-rich regions. Functional characterization of A. gigas fads2 showed that, similar to other teleosts, the A. gigas fads2 exhibited a predominant Δ6 activity complemented with some capacity for Δ8 desaturation. Given that A. gigas belongs to one of the oldest teleostei lineages, the Osteoglossomorpha, these findings offer a significant insight into the evolution LC-PUFA biosynthesis in teleosts.

A novel Acetyl-CoA synthetase short-chain subfamily member 1 (Acss1) gene indicates a dynamic history of paralogue retention and loss in vertebrates.

Acetyl-CoA short chain synthetases (ACSSs) are key enzymes in the activation of fatty acids through the formation of thioesters with CoA. Three subfamily members are currently recognized in the human genome, ACSS1, ACSS2 and ACSS3, all single copy genes. The mitochondrial isoform, Acss1, plays a key role in the metabolism of acetate for energy production. While the single copy condition has been accurately established in humans, the evolutionary history of the Acss1 subfamily in vertebrates has yet to be elucidated, in particular, the isoform diversity, origin and function. Through genome database mining we analyzed the diversity of Acss1 isoforms in vertebrate classes. We detected the presence of a novel Acss1 isoform, which we name Acss1B. This new gene, Acss1B, has a curious phylogenetic distribution being found in teleosts (except zebrafish), sauropsids (birds and reptiles) and probably chondrichthyes. In contrast Acss1A is found in all the investigated species, except the teleost medaka. By means of comparative genomics and phylogenetics we show that Acss1A and Acss1B were generated in the quadruplication of the vertebrate genome. In effect, we find that amphioxus, a pre-genome duplication chordate, has a single Acss1 gene in a genomic region equally related to a quadrupled vertebrate genomic set. Consequently, Acss1B has been lost in some teleosts, amphibians and mammals, while Acss1A is probably absent in medaka. The reported findings illustrate an especially dynamic pattern of paralogue retention and independent loss in vertebrate species involving the Acss1 subfamily, whose functional consequences in energy metabolism are as yet unknown.

A mollusk VDR/PXR/CAR-like (NR1J) nuclear receptor provides insight into ancient detoxification mechanisms.

The origin and diversification of the metazoan endocrine systems represents a fundamental research issue in biology. Nuclear receptors are critical components of these systems. A particular group named VDR/PXR/CAR (NR1I/J) is central in the mediation of detoxification responses. While orthologues have been thoroughly characterized in vertebrates, a sparse representation is currently available for invertebrates. Here, we provide the first isolation and characterization of a lophotrochozoan protostome VDR/PXR/CAR nuclear receptor (NR1J), in the estuarine bivalve the peppery furrow shell (Scrobicularia plana). Using a reporter gene assay, we evaluated the xenobiotic receptor plasticity comparing the human PXR with the S. plana NR1Jβ. Our results show that the molluscan receptor responds to a natural toxin (okadaic acid) in a similar fashion to that reported for other invertebrates. In contrast, the pesticide esfenvalerate displayed a unique response, since it down regulated transactivation at higher concentrations, while for triclosan no response was observed. Additionally, we uncovered lineage specific gene duplications and gene loss in the gene group encoding NRs in protostomes with likely impacts on the complexity of detoxification mechanisms across different phyla. Our findings pave the way for the development of multi-specific sensor tools to screen xenobiotic compounds acting via the NR1I/J group.

Characterization of the human ornithine transcarbamylase 3' untranslated regulatory region.

Mutations in the untranslated regulatory regions of genes may result in abnormal gene expression or transcriptional regulation. In this study, we characterize the ornithine transcarbamylase (OTC) mRNA isoforms of the X-linked OTC gene involved in the urea formation in the liver. Our data revealed that two major transcripts (OTC-t1 and OTC-t2) are more highly expressed than any of the other isoforms in all the tissues analyzed, though a longer transcript (OTC-t3) was also isolated and characterized from the brain sample. The OTC-t2 sequence fully matches the OTC mRNA reference sequence (NM_000531.5). All three isoforms use a canonical AAUAAA hexamer that is predicted to fold into a hairpin secondary structure which might be exposed to the cleavage and polyadenylation specificity factor. In addition, we observed that the OTC-t1 and OTC-t2 transcripts display heterogeneity at the cleavage sites in a tissue-dependent manner. Taken together, our data demonstrate that several mRNA isoforms are transcribed from the OTC gene, thereby indicating a wide degree of variability in post-transcriptional regulation.

LXRα and LXRβ nuclear receptors evolved in the common ancestor of gnathostomes

Nuclear receptors (NRs) regulate numerous aspects of the endocrine system. They mediate endogenous and exogenous cues, ensuring a homeostatic control of development and metabolism. Gene duplication, loss and mutation have shaped the repertoire and function of NRs in metazoans. Here, we examine the evolution of a pivotal orchestrator of cholesterol metabolism in vertebrates, the liver X receptors (LXRs). Previous studies suggested that LXRα and LXRβ genes emerged in the mammalian ancestor. However, we show through genome analysis and functional assay that bona fide LXRα and LXRβ orthologues are present in reptiles, coelacanth and chondrichthyans but not in cyclostomes. These findings show that LXR duplicated before gnathostome radiation, followed by asymmetric paralogue loss in some lineages. We suggest that a tighter control of cholesterol levels in vertebrates was achieved through the exploitation of a wider range of oxysterols, an ability contingent on ligand-binding pocket remodeling.

The Origin and Diversity of Cpt1 Genes in Vertebrate Species

The Carnitine palmitoyltransferase I (Cpt1) gene family plays a crucial role in energy homeostasis since it is required for the occurrence of fatty acid β-oxidation in the mitochondria. The exact gene repertoire in different vertebrate lineages is variable. Presently, four genes are documented: Cpt1a, also known as Cpt1a1, Cpt1a2; Cpt1b and Cpt1c. The later is considered a mammalian innovation resulting from a gene duplication event in the ancestor of mammals, after the divergence of sauropsids. In contrast, Cpt1a2 has been found exclusively in teleosts. Here, we reassess the overall evolutionary relationships of Cpt1 genes using a combination of approaches, including the survey of the gene repertoire in basal gnathostome lineages. Through molecular phylogenetics and synteny studies, we find that Cpt1c is most likely a rapidly evolving orthologue of Cpt1a2. Thus, Cpt1c is present in other lineages such as cartilaginous fish, reptiles, amphibians and the coelacanth. We show that genome duplications (2R) and variable rates of sequence evolution contribute to the history of Cpt1 genes in vertebrates. Finally, we propose that loss of Cpt1b is the likely cause for the unusual energy metabolism of elasmobranch.