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

How Many Sirtuin Genes Are Out There? Evolution of Sirtuin Genes in Vertebrates With a Description of a New Family Member

Studying the evolutionary history of gene families is a challenging and exciting task with a wide range of implications. In addition to exploring fundamental questions about the origin and evolution of genes, disentangling their evolution is also critical to those who do functional/structural studies to allow a deeper and more precise interpretation of their results in an evolutionary context. The sirtuin gene family is a group of genes that are involved in a variety of biological functions mostly related to aging. Their duplicative history is an open question, as well as the definition of the repertoire of sirtuin genes among vertebrates. Our results show a well-resolved phylogeny that represents an improvement in our understanding of the duplicative history of the sirtuin gene family. We identified a new sirtuin gene family member (SIRT3.2) that was apparently lost in the last common ancestor of amniotes but retained in all other groups of jawed vertebrates. According to our experimental analyses, elephant shark SIRT3.2 protein is located in mitochondria, the overexpression of which leads to an increase in cellular levels of ATP. Moreover, in vitro analysis demonstrated that it has deacetylase activity being modulated in a similar way to mammalian SIRT3. Our results indicate that there are at least eight sirtuin paralogs among vertebrates and that all of them can be traced back to the last common ancestor of the group that existed between 676 and 615 millions of years ago.

The gastric proton pump in gobiid and mudskipper fishes. Evidence of stomach loss?

Stomach loss has occurred independently multiple times during gnathostome evolution with notable frequency within the Teleostei. Significantly, this loss of acid-peptic digestion has been found to correlate with the secondary genomic loss of the gastric proton pump subunits (atp4a, atp4b) and pepsinogens/pepsins (pga, pgc). Gastric glands produce gastric juice containing the acid and pepsin and thus their presence is a hallmark feature of a digestive system capable of acid-peptic digestion. However, in gobiid fishes although oesogaster and gastric glands have been identified histologically, their functional significance has been questioned. In the present study we address whether the gastric proton pump is present and expressed in gastric glands of the goby Neogobius species (Gobiidae) and in members of the family Oxudercidae, a group of amphibious gobiid fishes commonly known as mudskippers (genera: Periophthalmus, Boleophthalmus, Periophthalmodon and Scartelaos). We confirmed the presence of gastric glands and have immunohistochemically localized gastric proton pump expression to these glands in Neogobius fluviatilis and Periophthalmus novemradiatus, Periophthalmus barbarus and Boleophthalmus boddarti. Genome analysis in Neogobius melanostomus, Periophthalmus magnuspinnatus, Scartelaos histophorus, Boleophthalmus pectinirostris, and Periophthalmodon schlosseri revealed the presence of both atp4a and atp4b subunit orthologues in all species in a conserved genomic loci organization. Moreover, it was possible to deduce that the complete open reading frame and the key functional amino acid residues are present. The conserved expression of the gastric proton pump provides clear evidence of the potential for gastric acid secretion indicating that acid digestion is retained in these gobiid fishes and not lost.

Inhibition of gastric acid secretion with omeprazole affects fish specific dynamic action and growth rate: Implications for the development of phenotypic stomach loss

An acid-secreting stomach provides many selective advantages to fish and other vertebrates; however, phenotypic stomach loss has occurred independently multiple times and is linked to loss of expression of both the gastric proton pump and the protease pepsin. Reasons underpinning stomach loss remain uncertain. Understanding the importance of gastric acid-secretion to the metabolic costs of digestion and growth will provide information about the metabolic expense of acid-production and performance. In this study, omeprazole, a well characterized gastric proton pump inhibitor, was used to simulate the agastric phenotype by significantly inhibiting gastric acidification in Nile tilapia. The effects on post-prandial metabolic rate and growth were assessed using intermittent flow respirometry and growth trials, respectively. Omeprazole reduced the duration (34.4%) and magnitude (34.5%) of the specific dynamic action and specific growth rate (21.3%) suggesting a decrease in digestion and assimilation of the meal. Gastric pH was measured in control and omeprazole treated fish to confirm that gastric acid secretion was inhibited for up to 12 h post-treatment (p < 0.05). Gastric evacuation measurements confirm a more rapid emptying of the stomach in omeprazole treated fish. These findings reinforce the importance of stomach acidification in digestion and growth and present a novel way of determining costs of gastric digestion.

Evolution of digestive enzymes and dietary diversification in birds

As the most species-rich class of tetrapod vertebrates, Aves possesses diverse feeding habits, with multiple origins of insectivory, carnivory, frugivory, nectarivory, granivory and omnivory. Since digestive enzymes mediate and limit energy and nutrient uptake, we hypothesized that genes encoding digestive enzymes have undergone adaptive evolution in birds. To test this general hypothesis, we identified 16 digestive enzyme genes (including seven carbohydrase genes (hepatic amy, pancreatic amy, salivary amy, agl, g6pc, gaa and gck), three lipase genes (cyp7a1, lipf and pnlip), two protease genes (ctrc and pgc), two lysozyme genes (lyz and lyg) and two chitinase genes (chia and chit1)) from the available genomes of 48 bird species. Among these 16 genes, three (salivary amy, lipf and chit1) were not found in all 48 avian genomes, which was further supported by our synteny analysis. Of the remaining 13 genes, eight were single-copy and five (chia, gaa, lyz, lyg and pgc) were multi-copy. Moreover, the multi-copy genes gaa, lyg and pgc were predicted to exhibit functional divergence among copies. Positively selected sites were detected in all of the analyzed digestive enzyme genes, except agl, g6pc, gaa and gck, suggesting that different diets may have favored differences in catalytic capacities of these enzymes. Furthermore, the analysis also revealed that the pancreatic amylase gene and one of the lipase genes (cyp7a1) have higher ω (the ratio of nonsynonymous to the synonymous substitution rates) values in species consuming a larger amount of seeds and meat, respectively, indicating an intense selection. In addition, the gck carbohydrase gene in species consuming a smaller amount of seeds, fruits or nectar, and a lipase gene (pnlip) in species consuming less meat were found to be under relaxed selection. Thus, gene loss, gene duplication, functional divergence, positive selection and relaxed selection have collectively shaped the evolution of digestive enzymes in birds, and the evolutionary flexibility of these enzymes may have facilitated their dietary diversification.

Molecular ontogeny of the stomach in the catshark Scyliorhinus canicula

The origin of extracellular digestion in metazoans was accompanied by structural and physiological alterations of the gut. These adaptations culminated in the differentiation of a novel digestive structure in jawed vertebrates, the stomach. Specific endoderm/mesenchyme signalling is required for stomach differentiation, involving the growth and transcription factors: 1) Shh and Bmp4, required for stomach outgrowth; 2) Barx1, Sfrps and Sox2, required for gastric epithelium development and 3) Cdx1 and Cdx2, involved in intestinal versus gastric identity. Thus, modulation of endoderm/mesenchyme signalling emerges as a plausible mechanism linked to the origin of the stomach. In order to gain insight into the ancient mechanisms capable of generating this structure in jawed vertebrates, we characterised the development of the gut in the catshark Scyliorhinus canicula. As chondrichthyans, these animals retained plesiomorphic features of jawed vertebrates, including a well-differentiated stomach. We identified a clear molecular regionalization of their embryonic gut, characterised by the expression of barx1 and sox2 in the prospective stomach region and expression of cdx1 and cdx2 in the prospective intestine. Furthermore, we show that gastric gland development occurs close to hatching, accompanied by the onset of gastric proton pump activity. Our findings favour a scenario in which the developmental mechanisms involved in the origin of the stomach were present in the common ancestor of chondrichthyans and osteichthyans.

A time-resolved multi-omic atlas of the developing mouse stomach

The mammalian stomach is structurally highly diverse and its organ functionality critically depends on a normal embryonic development. Although there have been several studies on the morphological changes during stomach development, a system-wide analysis of the underlying molecular changes is lacking. Here, we present a comprehensive, temporal proteome and transcriptome atlas of the mouse stomach at multiple developmental stages. Quantitative analysis of 12,108 gene products allows identifying three distinct phases based on changes in proteins and RNAs and the gain of stomach functions on a longitudinal time scale. The transcriptome indicates functionally important isoforms relevant to development and identifies several functionally unannotated novel splicing junction transcripts that we validate at the peptide level. Importantly, many proteins differentially expressed in stomach development are also significantly overexpressed in diffuse-type gastric cancer. Overall, our study provides a resource to understand stomach development and its connection to gastric cancer tumorigenesis.

Evaluation of Pepsinogen I as a Biomarker of Drug-induced Gastric Mucosal Injury in Cynomolgus Monkeys

Gastric mucosal injury is frequently observed in nonclinical studies of nonhuman primates. Because microscopic evaluation of stomach is generally a terminal procedure, our objective was to determine whether serum pepsinogen I (PG I) could serve as a noninvasive biomarker for detection of gastric mucosal injury in monkey. Serum PG I was measured using a commercial human immunoassay in cynomolgus monkeys ( n = 166) prior to dosing and/or terminally in 11 studies of up to 1 month duration. Mean ( SD) PG I values (ug/L) for monkeys with ( n = 59) and without ( n = 100) gastric mucosal degeneration were 101 (215) and 28 (12.6), respectively. For monkeys with baseline and terminal PG I data, mean ( SD) fold change (ratio of terminal to baseline PG I) for monkeys with ( n = 57) and without ( n = 76) glandular degeneration were 4.1 (11.3) and 1 (0.28). Receiver operating characteristic area under the curve (AUC) data demonstrated moderate diagnostic accuracy for serum PG I for glandular degeneration, AUC ( SE) 0.789 (0.04), with improved diagnostic accuracy as a fold change of baseline, AUC ( SE) 0.816 (0.04), consistent with the large interindividual but low intraindividual variability of serum PG I values in control monkeys. These data demonstrate that serum PG I is a useful biomarker of drug-induced gastric mucosal injury in the cynomolgus monkey.

The Gastric Phenotype in the Cypriniform Loaches: A Case of Reinvention?

The stomach, which is characterized by acid peptic digestion in vertebrates, has been lost secondarily multiple times in the evolution of the teleost fishes. The Cypriniformes are largely seen as an agastric order; however, within the superfamily Cobitoidea, the closely related sister groups Nemacheilidae and Balitoridae have been identified as gastric families. The presence of these most recently diverged gastric families in an otherwise agastric clade indicates that either multiple (>2–3) loss events occurred with the Cyprinidae, Catostomidae and Cobitidae, or that gastric reinvention arose in a recent ancestor of the Nemacheilidae/Balitoridae sister clade. In the present study, the foregut regions of Cobitidae, Nemacheilidae/Balitoridae and the ancestral Botiidae family members were examined for the presence of gastric glands and gastric proton pump (Atp4a) α subunit expression by histology and immunohistochemistry respectively. Atp4a gene expression was assessed by reverse transcriptase-polymerase chain reaction (RT-PCR). Gastric glands expressing apical H⁺/K⁺-ATPase α subunit and isolated partial sequences of atp4a, identified using degenerate primers showing clear orthology to other vertebrate atp4a sequences, were detected in representative species from Nemacheilidae/ Balitoridae and Botiidae, but not Cobitidae (Misgurnus anguillicaudatus). In summary, we provide evidence for an uninterrupted gastric evolutionary lineage in the Cobitoidea, making it highly improbable that the stomach was reinvented in the Nemacheilidae/Balitoridae clade consistent with Dollo’s principle. These results also indicate that the gastric trait may be present elsewhere in the Cobitoidea.

Elevated Gene Copy Number Does Not Always Explain Elevated Amylase Activities in Fishes

Amylase activity variation in the guts of several model organisms appears to be explained by amylase gene copy number variation. We tested the hypothesis that amylase gene copy number is always elevated in animals with high amylolytic activity. We therefore sequenced the amylase genes and examined amylase gene copy number in prickleback fishes (family Stichaeidae) with different diets including two species of convergently evolved herbivores with the elevated amylase activity phenotype. We found elevated amylase gene copy number (six haploid copies) with sequence variation among copies in one herbivore (Cebidichthys violaceus) and modest gene copy number (two to three haploid copies) with little sequence variation in the remaining taxa, which included herbivores, omnivores, and a carnivore. Few functional differences in amylase biochemistry were observed, and previous investigations showed similar digestibility among the convergently evolved herbivores with differing amylase genetics. Hence, the phenotype of elevated amylase activity can be achieved by different mechanisms (i.e., elevated expression of fewer genes, increased gene copy number, or expression of more efficient amylase proteins) with similar results. Phylogenetic and comparative genomic analyses of available fish amylase genes show mostly lineage-specific duplication events leading to gene copy number variation, although a whole-genome duplication event or chromosomal translocation may have produced multiple amylase copies in the Ostariophysi, again showing multiple routes to the same result.

Purification and molecular cloning of aspartic proteinases from the stomach of adult Japanese fire belly newts, Cynops pyrrhogaster

Six aspartic proteinase precursors, a pro-cathepsin E (ProCatE) and five pepsinogens (Pgs), were purified from the stomach of adult newts (Cynops pyrrhogaster). On sodium dodecylsulfate-polyacrylamide gel electrophoresis, the molecular weights of the Pgs and active enzymes were 37-38 kDa and 31-34 kDa, respectively. The purified ProCatE was a dimer whose subunits were connected by a disulphide bond. cDNA cloning by polymerase chain reaction and subsequent phylogenetic analysis revealed that three of the purified Pgs were classified as PgA and the remaining two were classified as PgBC belonging to C-type Pg. Our results suggest that PgBC is one of the major constituents of acid protease in the urodele stomach. We hypothesize that PgBC is an amphibian-specific Pg that diverged during its evolutional lineage. PgBC was purified and characterized for the first time. The purified urodele pepsin A was completely inhibited by equal molar units of pepstatin A. Conversely, the urodele pepsin BC had low sensitivity to pepstatin A. In acidic condition, the activation rates of newt pepsin A and BC were similar to those of mammalian pepsin A and C1, respectively. Our results suggest that the enzymological characters that distinguish A- and C-type pepsins appear to be conserved in mammals and amphibians.

Comparative analysis of zebrafish bone morphogenetic proteins 2, 4 and 16: molecular and evolutionary perspectives

BMP2, BMP4 and BMP16 form a subfamily of bone morphogenetic proteins acting as pleiotropic growth factors during development and as bone inducers during osteogenesis. BMP16 is the most recent member of this subfamily and basic data regarding protein structure and function, and spatio-temporal gene expression is still scarce. In this work, insights on BMP16 were provided through the comparative analysis of structural and functional data for zebrafish BMP2a, BMP2b, BMP4 and BMP16 genes and proteins, determined from three-dimensional models, patterns of gene expression during development and in adult tissues, regulation by retinoic acid and capacity to activate BMP-signaling pathway. Structures of Bmp2a, Bmp2b, Bmp4 and Bmp16 were found to be remarkably similar; with residues involved in receptor binding being highly conserved. All proteins could activate the BMP-signaling pathway, suggesting that they share a common function. On the contrary, stage- and tissue-specific expression of bmp2, bmp4 and bmp16 suggested the genes might be differentially regulated (e.g. different transcription factors, enhancers and/or regulatory modules) but also that they are involved in distinct physiological processes, although with the same function. Retinoic acid, a morphogen known to interact with BMP-signaling during bone formation, was shown to down-regulate the expression of bmp2, bmp4 and bmp16, although to different extents. Taxonomic and phylogenetic analyses indicated that bmp16 diverged before bmp2 and bmp4, is not restricted to teleost fish lineage as previously reported, and that it probably arose from a whole genomic duplication event that occurred early in vertebrate evolution and disappeared in various tetrapod lineages through independent events.

Extensive expansion of A1 family aspartic proteinases in fungi revealed by evolutionary analyses of 107 complete eukaryotic proteomes

The A1 family of eukaryotic aspartic proteinases (APs) forms one of the 16 AP families. Although one of the best characterized families, the recent increase in genome sequence data has revealed many fungal AP homologs with novel sequence characteristics. This study was performed to explore the fungal AP sequence space and to obtain an in-depth understanding of fungal AP evolution. Using a comprehensive phylogeny of approximately 700 AP sequences from the complete proteomes of 87 fungi and 20 nonfungal eukaryotes, 11 major clades of APs were defined of which clade I largely corresponds to the A1A subfamily of pepsin-archetype APs. Clade II largely corresponds to the A1B subfamily of nepenthesin-archetype APs. Remarkably, the nine other clades contain only fungal APs, thus indicating that fungal APs have undergone a large sequence diversification. The topology of the tree indicates that fungal APs have been subject to both “birth and death” evolution and “functional redundancy and diversification.” This is substantiated by coclustering of certain functional sequence characteristics. A meta-analysis toward the identification of Cluster Determining Positions (CDPs) was performed in order to investigate the structural and biochemical basis for diversification. Seven CDPs contribute to the secondary structure of the enzyme. Three other CDPs are found in the vicinity of the substrate binding cleft. Tree topology, the large sequence variation among fungal APs, and the apparent functional diversification suggest that an amendment to update the current A1 AP classification based on a comprehensive phylogenetic clustering might contribute to refinement of the classification in the MEROPS peptidase database.

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.

Evolution of the human gastrokine locus and confounding factors regarding the pseudogenicity of GKN3

In a screen for genes expressed specifically in gastric mucous neck cells, we identified GKN3, the recently discovered third member of the gastrokine family. We present confirmatory mouse data and novel porcine data showing that mouse GKN3 expression is confined to mucous cells of the corpus neck and antrum base and is prominently expressed in metaplastic lesions. GKN3 was proposed originally to be expressed in some human populations and a pseudogene in others. To investigate that hypothesis, we studied human GKN3 evolution in the context of its paralogous genomic neighbors, GKN1 and GKN2. Haplotype analysis revealed that GKN3 mimics GKN2 in patterns of exonic SNP allocation, whereas GKN1 appeared to be more stringently selected. GKN3 showed signatures of both directional selection and population based selective sweeps in humans. One such selective sweep includes SNP rs10187256, originally identified as an ancestral tryptophan to premature STOP codon mutation. The derived (nonancestral) allele went to fixation in Asia. We show that another SNP, rs75578132, identified 5 bp downstream of rs10187256, exhibits a second selective sweep in almost all Europeans, some Latinos, and some Africans, possibly resulting from a reintroduction of European genes during African colonization. Finally, we identify a mutation that would destroy the splice donor site in the putative exon3-intron3 boundary, which occurs in all human genomes examined to date. Our results highlight a stomach-specific human genetic locus, which has undergone various selective sweeps across European, Asian, and African populations and thus reflects geographic and ethnic patterns in genome evolution.

Molecular cloning of pepsinogens A and C from adult newt (Cynops pyrrhogaster) stomach

The full-length cDNAs of three pepsinogens (Pgs) were cloned from the stomach of newt, Cynops pyrrhogaster, and nucleotide sequences of the full-length cDNAs were determined. Molecular phylogenetic analysis showed that two Pgs, named PgC1 and PgC2, belong to the pepsinogen C group, and one Pg, named PgA, belongs to the pepsinogen A group. The sequences contain an open reading frame (ORF) encoding 385 amino acid residues for PgC1, 383 amino acid residues for PgC2 and 377 amino acid residues for PgA. In addition, all of the three amino acid sequences conserve some unique characteristics such as six cysteine residues and putative active site two aspartic acid residues. All of the pepsinogen mRNAs were detected in the stomach by RT-PCR but not in other organs. Although a slight difference at the time of the start of expression was seen among the three pepsinogen genes, all of them were expressed in the larval stage after hatching. This is the first report on cloning of pepsinogens from urodele stomach.