The evolution of vitamin C biosynthesis and transport in animals

Background

Vitamin C (VC) is an indispensable antioxidant and co-factor for optimal function and development of eukaryotic cells. In animals, VC can be synthesized by the organism, acquired through the diet, or both. In the single VC synthesis pathway described in animals, the penultimate step is catalysed by Regucalcin, and the last step by l-gulonolactone oxidase (GULO). The GULO gene has been implicated in VC synthesis only, while Regucalcin has been shown to have multiple functions in mammals.

Results

Both GULO and Regucalcin can be found in non-bilaterian, protostome and deuterostome species. Regucalcin, as here shown, is involved in multiple functions such as VC synthesis, calcium homeostasis, and the oxidative stress response in both Deuterostomes and Protostomes, and in insects in receptor-mediated uptake of hexamerin storage proteins from haemolymph. In Insecta and Nematoda, however, there is no GULO gene, and in the latter no Regucalcin gene, but species from these lineages are still able to synthesize VC, implying at least one novel synthesis pathway. In vertebrates, SVCT1, a gene that belongs to a family with up to five members, as here shown, is the only gene involved in the uptake of VC in the gut. This specificity is likely the result of a subfunctionalization event that happened at the base of the Craniata subphylum. SVCT-like genes present in non-Vertebrate animals are likely involved in both VC and nucleobase transport. It is also shown that in lineages where GULO has been lost, SVCT1 is now an essential gene, while in lineages where SVCT1 gene has been lost, GULO is now an essential gene.

Conclusions

The simultaneous study, for the first time, of GULO, Regucalcin and SVCTs evolution provides a clear picture of VC synthesis/acquisition and reveals very different selective pressures in different animal taxonomic groups.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12862-022-02040-7.

Advances in Novel Animal Vitamin C Biosynthesis Pathways and the Role of Prokaryote-Based Inferences to Understand Their Origin

Vitamin C (VC) is an essential nutrient required for the optimal function and development of many organisms. VC has been studied for many decades, and still today, the characterization of its functions is a dynamic scientific field, mainly because of its commercial and therapeutic applications. In this review, we discuss, in a comparative way, the increasing evidence for alternative VC synthesis pathways in insects and nematodes, and the potential of myo-inositol as a possible substrate for this metabolic process in metazoans. Methodological approaches that may be useful for the future characterization of the VC synthesis pathways of Caenorhabditis elegans and Drosophila melanogaster are here discussed. We also summarize the current distribution of the eukaryote aldonolactone oxidoreductases gene lineages, while highlighting the added value of studies on prokaryote species that are likely able to synthesize VC for both the characterization of novel VC synthesis pathways and inferences on the complex evolutionary history of such pathways. Such work may help improve the industrial production of VC.

Predictive Models of within- and between-Species SARS-CoV-2 Transmissibility

Viruses from the Coronaviridae family have been reported to infect a large range of hosts, including humans. The latest human-infecting coronavirus, SARS-CoV-2, turned into a pandemic and subtypes with different transmissibility have appeared since then. The SARS-CoV-2 Spike (S) protein interacts with the angiotensin-converting enzyme 2 (ACE2) host receptor, and thus, in silico models, based on the structural features of the SARS-CoV-2 S protein–ACE2 receptor complex, as well as ACE2 amino acid patterns, may be used to predict the within- and between-species transmissibility of SARS-CoV-2 subtypes. Here, it is shown that at the beginning of the pandemic, the SARS-CoV-2 S protein was, as expected for a virus that just jumped the species barrier, ill-adapted to the human ACE2 receptor, and that the replacement of one SARS-CoV-2 variant by another is partially due to a better fitting of the S protein–human ACE2 complex. Moreover, it is shown that mutations that are predicted to lead to a better fit have increased in the population due to positive selection. It is also shown that the number of ACE2-interfacing residues is positively correlated with the transmissibility rate of SARS-CoV-2 variants. Finally, it is shown that the number of species that are susceptible to infection by SARS-CoV-2, and that could be a reservoir for this virus, is likely higher than previously thought.

Target of Rapamycin Complex 1 (TORC1), Protein Kinase A (PKA) and Cytosolic pH Regulate a Transcriptional Circuit for Lipid Droplet Formation

Lipid droplets (LDs) are ubiquitous organelles that fulfill essential roles in response to metabolic cues. The identification of several neutral lipid synthesizing and regulatory protein complexes have propelled significant advance on the mechanisms of LD biogenesis in the endoplasmic reticulum (ER). However, our understanding of signaling networks, especially transcriptional mechanisms, regulating membrane biogenesis is very limited. Here, we show that the nutrient-sensing Target of Rapamycin Complex 1 (TORC1) regulates LD formation at a transcriptional level, by targeting DGA1 expression, in a Sit4-, Mks1-, and Sfp1-dependent manner. We show that cytosolic pH (pHc), co-regulated by the plasma membrane H+-ATPase Pma1 and the vacuolar ATPase (V-ATPase), acts as a second messenger, upstream of protein kinase A (PKA), to adjust the localization and activity of the major transcription factor repressor Opi1, which in turn controls the metabolic switch between phospholipid metabolism and lipid storage. Together, this work delineates hitherto unknown molecular mechanisms that couple nutrient availability and pHc to LD formation through a transcriptional circuit regulated by major signaling transduction pathways.

BDBM 1.0: A Desktop Application for Efficient Retrieval and Processing of High-Quality Sequence Data and Application to the Identification of the Putative Coffea S-Locus

Nowadays, bioinformatics is one of the most important areas in modern biology and the creation of high-quality scientific software supporting this recent research area is one of the core activities of many researchers. In this context, high-quality sequence datasets are needed to perform inferences on the evolution of species, genes, and gene families, or to get evidence for adaptive amino acid evolution, among others. Nevertheless, sequence data are very often spread over several databases, many useful genomes and transcriptomes are non-annotated, the available annotation is not for the desired coding sequence isoform, and/or is unlikely to be accurate. Moreover, although the FASTA text-based format is quite simple and usable by most software applications, there are a number of issues that may be critical depending on the software used to analyse such files. Therefore, researchers without training in informatics often use a fraction of all available data. The above issues can be addressed using already available software applications, but there is no easy-to-use single piece of software that allows performing all these tasks within the same graphical interface, such as the one here presented, named BDBM (Blast DataBase Manager). BDBM can be used to efficiently get gene sequences from annotated and non-annotated genomes and transcriptomes. Moreover, it can be used to look for alternatives to existing annotations and to easily create reliable custom databases. Such databases are essential to prepare high-quality datasets. The analyses that we have performed on the Coffea canephora genome using BDBM aimed at the identification of the S-locus region (that harbours the genes involved in gametophytic self-incompatibility) led to the conclusion that there are two likely regions, one on chromosome 2 (around region 6600000-6650000), and another on chromosome 5 (around 15830000-15930000). Such findings are discussed in the context of the Rubiaceae gametophytic self-incompatibility evolution.

Bioinformatics Protocols for Quickly Obtaining Large-Scale Data Sets for Phylogenetic Inferences

Useful insight into the evolution of genes and gene families can be provided by the analysis of all available genome datasets rather than just a few, which are usually those of model species. Handling and transforming such datasets into the desired format for downstream analyses is, however, often a difficult and time-consuming task for researchers without a background in informatics. Therefore, we present two simple and fast protocols for data preparation, using an easy-to-install, open-source, cross-platform software application with user-friendly, rich graphical user interface (SEDA; http://www.sing-group.org/seda/index.html ). The first protocol is a substantial improvement over one recently published (López-Fernández et al. Practical applications of computational biology and bioinformatics, 12th International conference. Springer, Cham, pp 88-96 (2019)[1]), which was used to study the evolution of GULO, a gene that encodes the enzyme responsible for the last step of vitamin C synthesis. In this paper, we show how the sequence data file used for the phylogenetic analyses can now be obtained much faster by changing the way coding sequence isoforms are removed, using the newly implemented SEDA operation "Remove isoforms". This protocol can be used to easily show that putative functional GULO genes are present in several Prostotomian groups such as Molluscs, Priapulida and Arachnida. Such findings could have been easily missed if only a few Protostomian model species had been used. The second protocol allowed us to identify positively selected amino acid sites in a set of 19 primate HLA immunity genes. Interestingly, the proteins encoded by MHC class II genes can show just as many positively selected amino acid sites as those encoded by classical MHC class I genes. Although a significant percentage of codons, which can be as high as 14.8%, are evolving under positive selection, the main mode of evolution of HLA immunity genes is purifying selection. Using a large number of primate species, the probability of missing the identification of positively selected amino acid sites is lower. Both projects were performed in less than one week, and most of the time was spent running the analyses rather than preparing the files. Such protocols can be easily adapted to answer many other questions using a phylogenetic approach.

GC4S: A bioinformatics-oriented Java software library of reusable graphical user interface components

Modern bioinformatics and computational biology are fields of study driven by the availability of effective software required for conducting appropriate research tasks. Apart from providing reliable and fast implementations of different data analysis algorithms, these software applications should also be clear and easy to use through proper user interfaces, providing appropriate data management and visualization capabilities. In this regard, the user experience obtained by interacting with these applications via their Graphical User Interfaces (GUI) is a key factor for their final success and real utility for researchers. Despite the existence of different packages and applications focused on advanced data visualization, there is a lack of specific libraries providing pertinent GUI components able to help scientific bioinformatics software developers. To that end, this paper introduces GC4S, a bioinformatics-oriented collection of high-level, extensible, and reusable Java GUI elements specifically designed to speed up bioinformatics software development. Within GC4S, developers of new applications can focus on the specific GUI requirements of their projects, relying on GC4S for generalities and abstractions. GC4S is free software distributed under the terms of GNU Lesser General Public License and both source code and documentation are publicly available at http://www.sing-group.org/gc4s.

Inferences on specificity recognition at the Malus×domestica gametophytic self-incompatibility system

In Malus × domestica (Rosaceae) the product of each SFBB gene (the pollen component of the gametophytic self-incompatibility (GSI) system) of a S-haplotype (the combination of pistil and pollen genes that are linked) interacts with a sub-set of non-self S-RNases (the pistil component), but not with the self S-RNase. To understand how the Malus GSI system works, we identified 24 SFBB genes expressed in anthers, and determined their gene sequence in nine M. domestica cultivars. Expression of these SFBBs was not detected in the petal, sepal, filament, receptacle, style, stigma, ovary or young leaf. For all SFBBs (except SFBB15), identical sequences were obtained only in cultivars having the same S-RNase. Linkage with a particular S-RNase was further established using the progeny of three crosses. Such data is needed to understand how other genes not involved in GSI are affected by the S-locus region. To classify SFBBs specificity, the amino acids under positive selection obtained when performing intra-haplotypic analyses were used. Using this information and the previously identified S-RNase positively selected amino acid sites, inferences are made on the S-RNase amino acid properties (hydrophobicity, aromatic, aliphatic, polarity, and size), at these positions, that are critical features for GSI specificity determination.

The Drosophila melanogaster Muc68E Mucin Gene Influences Adult Size, Starvation Tolerance, and Cold Recovery

Mucins have been implicated in many different biological processes, such as protection from mechanical damage, microorganisms, and toxic molecules, as well as providing a luminal scaffold during development. Nevertheless, it is conceivable that mucins have the potential to modulate food absorption as well, and thus contribute to the definition of several important phenotypic traits. Here we show that the Drosophila melanogaster Muc68E gene is 40- to 60-million-yr old, and is present in Drosophila species of the subgenus Sophophora only. The central repeat region of this gene is fast evolving, and shows evidence for repeated expansions/contractions. This and/or frequent gene conversion events lead to the homogenization of its repeats. The amino acid pattern P[ED][ED][ST][ST][ST] is found in the repeat region of Muc68E proteins from all Drosophila species studied, and can occur multiple times within a single conserved repeat block, and thus may have functional significance. Muc68E is a nonessential gene under laboratory conditions, but Muc68E mutant flies are smaller and lighter than controls at birth. However, at 4 d of age, Muc68E mutants are heavier, recover faster from chill-coma, and are more resistant to starvation than control flies, although they have the same percentage of lipids as controls. Mutant flies have enlarged abdominal size 1 d after chill-coma recovery, which is associated with higher lipid content. These results suggest that Muc68E has a role in metabolism modulation, food absorption, and/or feeding patterns in larvae and adults, and under normal and stress conditions. Such biological function is novel for mucin genes.

No evidence for Fabaceae Gametophytic self-incompatibility being determined by Rosaceae, Solanaceae, and Plantaginaceae S-RNase lineage genes

BACKGROUND

Fabaceae species are important in agronomy and livestock nourishment. They have a long breeding history, and most cultivars have lost self-incompatibility (SI), a genetic barrier to self-fertilization. Nevertheless, to improve legume crop breeding, crosses with wild SI relatives of the cultivated varieties are often performed. Therefore, it is fundamental to characterize Fabaceae SI system(s). We address the hypothesis of Fabaceae gametophytic (G)SI being RNase based, by recruiting the same S-RNase lineage gene of Rosaceae, Solanaceae or Plantaginaceae SI species.

RESULTS

We first identify SSK1 like genes (described only in species having RNase based GSI), in the Trifolium pratense, Medicago truncatula, Cicer arietinum, Glycine max, and Lupinus angustifolius genomes. Then, we characterize the S-lineage T2-RNase genes in these genomes. In T. pratense, M. truncatula, and C. arietinum we identify S-RNase lineage genes that in phylogenetic analyses cluster with Pyrinae S-RNases. In M. truncatula and C. arietinum genomes, where large scaffolds are available, these sequences are surrounded by F-box genes that in phylogenetic analyses also cluster with S-pollen genes. In T. pratense the S-RNase lineage genes show, however, expression in tissues not involved in GSI. Moreover, levels of diversity are lower than those observed for other S-RNase genes. The M. truncatula and C. arietinum S-RNase and S-pollen like genes phylogenetically related to Pyrinae S-genes, are also expressed in tissues other than those involved in GSI. To address if other T2-RNases could be determining Fabaceae GSI, here we obtained a style with stigma transcriptome of Cytisus striatus, a species that shows significant difference on the percentage of pollen growth in self and cross-pollinations. Expression and polymorphism analyses of the C. striatus S-RNase like genes revealed that none of these genes, is the S-pistil gene.

CONCLUSION

We find no evidence for Fabaceae GSI being determined by Rosaceae, Solanaceae, and Plantaginaceae S-RNase lineage genes. There is no evidence that T2-RNase lineage genes could be determining GSI in C. striatus. Therefore, to characterize the Fabaceae S-pistil gene(s), expression analyses, levels of diversity, and segregation analyses in controlled crosses are needed for those genes showing high expression levels in the tissues where GSI occurs.

Convergent evolution at the gametophytic self-incompatibility system in Malus and Prunus

S-RNase-based gametophytic self-incompatibility (GSI) has evolved once before the split of the Asteridae and Rosidae. This conclusion is based on the phylogenetic history of the S-RNase that determines pistil specificity. In Rosaceae, molecular characterizations of Prunus species, and species from the tribe Pyreae (i.e., Malus, Pyrus, Sorbus) revealed different numbers of genes determining S-pollen specificity. In Prunus only one pistil and pollen gene determine GSI, while in Pyreae there is one pistil but multiple pollen genes, implying different specificity recognition mechanisms. It is thus conceivable that within Rosaceae the genes involved in GSI in the two lineages are not orthologous but possibly paralogous. To address this hypothesis we characterised the S-RNase lineage and S-pollen lineage genes present in the genomes of five Rosaceae species from three genera: M. × domestica (apple, self-incompatible (SI); tribe Pyreae), P. persica (peach, self-compatible (SC); Amygdaleae), P. mume (mei, SI; Amygdaleae), Fragaria vesca (strawberry, SC; Potentilleae), and F. nipponica (mori-ichigo, SI; Potentilleae). Phylogenetic analyses revealed that the Malus and Prunus S-RNase and S-pollen genes belong to distinct gene lineages, and that only Prunus S-RNase and SFB-lineage genes are present in Fragaria. Thus, S-RNase based GSI system of Malus evolved independently from the ancestral system of Rosaceae. Using expression patterns based on RNA-seq data, the ancestral S-RNase lineage gene is inferred to be expressed in pistils only, while the ancestral S-pollen lineage gene is inferred to be expressed in tissues other than pollen.

The unique regulation of iron-sulfur cluster biogenesis in a Gram-positive bacterium

Iron-sulfur clusters function as cofactors of a wide range of proteins, with diverse molecular roles in both prokaryotic and eukaryotic cells. Dedicated machineries assemble the clusters and deliver them to the final acceptor molecules in a tightly regulated process. In the prototypical Gram-negative bacterium Escherichia coli, the two existing iron-sulfur cluster assembly systems, iron-sulfur cluster (ISC) and sulfur assimilation (SUF) pathways, are closely interconnected. The ISC pathway regulator, IscR, is a transcription factor of the helix-turn-helix type that can coordinate a [2Fe-2S] cluster. Redox conditions and iron or sulfur availability modulate the ligation status of the labile IscR cluster, which in turn determines a switch in DNA sequence specificity of the regulator: cluster-containing IscR can bind to a family of gene promoters (type-1) whereas the clusterless form recognizes only a second group of sequences (type-2). However, iron-sulfur cluster biogenesis in Gram-positive bacteria is not so well characterized, and most organisms of this group display only one of the iron-sulfur cluster assembly systems. A notable exception is the unique Gram-positive dissimilatory metal reducing bacterium Thermincola potens, where genes from both systems could be identified, albeit with a diverging organization from that of Gram-negative bacteria. We demonstrated that one of these genes encodes a functional IscR homolog and is likely involved in the regulation of iron-sulfur cluster biogenesis in T. potens. Structural and biochemical characterization of T. potens and E. coli IscR revealed a strikingly similar architecture and unveiled an unforeseen conservation of the unique mechanism of sequence discrimination characteristic of this distinctive group of transcription regulators.

The founding of Mauritian endemic coffee trees by a synchronous long-distance dispersal event

The stochastic process of long-distance dispersal is the exclusive means by which plants colonize oceanic islands. Baker's rule posits that self-incompatible plant lineages are unlikely to successfully colonize oceanic islands because they must achieve a coordinated long-distance dispersal of sufficiently numerous individuals to establish an outcrossing founder population. Here, we show for the first time that Mauritian Coffea species are self-incompatible and thus represent an exception to Baker's rule. The genus Coffea (Rubiaceae) is composed of approximately 124 species with a paleotropical distribution. Phylogenetic evidence strongly supports a single colonization of the oceanic island of Mauritius from either Madagascar or Africa. We employ Bayesian divergence time analyses to show that the colonization of Mauritius was not a recent event. We genotype S-RNase alleles from Mauritian endemic Coffea, and using S-allele gene genealogies, we show that the Mauritian allelic diversity is confined to just seven deeply divergent Coffea S-RNase allelic lineages. Based on these data, we developed an individual-based model and performed a simulation study to estimate the most likely number of founding individuals involved in the colonization of Mauritius. Our simulations show that to explain the observed S-RNase allelic diversity, the founding population was likely composed of fewer than 31 seeds that were likely synchronously dispersed from an ancestral mainland species.

Evolution stings: the origin and diversification of scorpion toxin peptide scaffolds

The episodic nature of natural selection and the accumulation of extreme sequence divergence in venom-encoding genes over long periods of evolutionary time can obscure the signature of positive Darwinian selection. Recognition of the true biocomplexity is further hampered by the limited taxon selection, with easy to obtain or medically important species typically being the subject of intense venom research, relative to the actual taxonomical diversity in nature. This holds true for scorpions, which are one of the most ancient terrestrial venomous animal lineages. The family Buthidae that includes all the medically significant species has been intensely investigated around the globe, while almost completely ignoring the remaining non-buthid families. Australian scorpion lineages, for instance, have been completely neglected, with only a single scorpion species (Urodacus yaschenkoi) having its venom transcriptome sequenced. Hence, the lack of venom composition and toxin sequence information from an entire continent’s worth of scorpions has impeded our understanding of the molecular evolution of scorpion venom. The molecular origin, phylogenetic relationships and evolutionary histories of most scorpion toxin scaffolds remain enigmatic. In this study, we have sequenced venom gland transcriptomes of a wide taxonomical diversity of scorpions from Australia, including buthid and non-buthid representatives. Using state-of-art molecular evolutionary analyses, we show that a majority of CSα/β toxin scaffolds have experienced episodic influence of positive selection, while most non-CSα/β linear toxins evolve under the extreme influence of negative selection. For the first time, we have unraveled the molecular origin of the major scorpion toxin scaffolds, such as scorpion venom single von Willebrand factor C-domain peptides (SV-SVC), inhibitor cystine knot (ICK), disulphide-directed beta-hairpin (DDH), bradykinin potentiating peptides (BPP), linear non-disulphide bridged peptides and antimicrobial peptides (AMP). We have thus demonstrated that even neglected lineages of scorpions are a rich pool of novel biochemical components, which have evolved over millions of years to target specific ion channels in prey animals, and as a result, possess tremendous implications in therapeutics.

The Drosophila melanogaster methuselah gene: a novel gene with ancient functions

The Drosophila melanogaster G protein-coupled receptor gene, methuselah (mth), has been described as a novel gene that is less than 10 million years old. Nevertheless, it shows a highly specific expression pattern in embryos, larvae, and adults, and has been implicated in larval development, stress resistance, and in the setting of adult lifespan, among others. Although mth belongs to a gene subfamily with 16 members in D. melanogaster, there is no evidence for functional redundancy in this subfamily. Therefore, it is surprising that a novel gene influences so many traits. Here, we explore the alternative hypothesis that mth is an old gene. Under this hypothesis, in species distantly related to D. melanogaster, there should be a gene with features similar to those of mth. By performing detailed phylogenetic, synteny, protein structure, and gene expression analyses we show that the D. virilis GJ12490 gene is the orthologous of mth in species distantly related to D. melanogaster. We also show that, in D. americana (a species of the virilis group of Drosophila), a common amino acid polymorphism at the GJ12490 orthologous gene is significantly associated with developmental time, size, and lifespan differences. Our results imply that GJ12490 orthologous genes are candidates for developmental time and lifespan differences in Drosophila in general.

Patterns of evolution at the gametophytic self-incompatibility Sorbus aucuparia (Pyrinae) S pollen genes support the non-self recognition by multiple factors model

S-RNase-based gametophytic self-incompatibility evolved once before the split of the Asteridae and Rosidae. In Prunus (tribe Amygdaloideae of Rosaceae), the self-incompatibility S-pollen is a single F-box gene that presents the expected evolutionary signatures. In Malus and Pyrus (subtribe Pyrinae of Rosaceae), however, clusters of F-box genes (called SFBBs) have been described that are expressed in pollen only and are linked to the S-RNase gene. Although polymorphic, SFBB genes present levels of diversity lower than those of the S-RNase gene. They have been suggested as putative S-pollen genes, in a system of non-self recognition by multiple factors. Subsets of allelic products of the different SFBB genes interact with non-self S-RNases, marking them for degradation, and allowing compatible pollinations. This study performed a detailed characterization of SFBB genes in Sorbus aucuparia (Pyrinae) to address three predictions of the non-self recognition by multiple factors model. As predicted, the number of SFBB genes was large to account for the many S-RNase specificities. Secondly, like the S-RNase gene, the SFBB genes were old. Thirdly, amino acids under positive selection-those that could be involved in specificity determination-were identified when intra-haplotype SFBB genes were analysed using codon models. Overall, the findings reported here support the non-self recognition by multiple factors model.