Positive selection within sperm-egg adhesion domains of fertilin: an ADAM gene with a potential role in fertilization

Genes with a role in fertilization show a common pattern of rapid evolution. The role played by positive selection versus lack of selective constraints has been more difficult to establish. One problem arises from attempts to detect selection in an overall gene sequence analysis. I have analyzed the pattern of molecular evolution of fertilin, a gene coding for a heterodimeric sperm protein belonging to the ADAM (A disintegrin and A metalloprotease) gene family. A nonsynonymous to synonymous rate ratio (d(N)/d(S)) analysis for different protein domains of fertilin alpha and fertilin beta showed d(N)/d(S) < 1, suggesting that purifying selection has shaped fertilin's evolution. However, an analysis of the distribution of single positively selected codon sites using phylogentic analysis by maximum likelihood (PAML) showed sites within adhesion domains (disintegrin and cysteine-rich) of fertilin beta evolving under positive selection. The region 3' to the EGF-like domain of fertilin alpha, where the transmembrane and cytoplasmic tail regions are supposed to be localized, showed higher d(N) and d(S) than any other fertilin alpha region. However, it was not possible to identify positively selected codon sites due to ambiguous alignments of the carboxy-end region (ClustalX vs. DiAlign2). When this region was excluded from the PAML analysis, most single positively selected codon sites were concentrated within adhesion domains (cysteine-rich and EGF-like). The use of an ancestral sequence prior to a recent duplication event of fertilin alpha among non-Hominidae primates (Macaca, Papio, and Saguinus) revealed that the duplication is partially responsible for masking the detection of positively selected sites within the disintegrin domain. Finally, most ADAM genes with a potential role in sperm maturation and/or fertilization showed significantly higher d(N) estimates than other ADAM genes.     

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.     

Detecting the undetectable: uncovering duplicated segments in Arabidopsis by comparison with rice

Genome analysis shows that large-scale gene duplications have occurred in fungi, animals and plants, creating genomic regions that show similarity in gene content and order. However, the high frequency of gene loss reduces colinearity resulting in duplicated regions that, in the extreme, no longer share homologous genes. Here, we show that by comparison with an appropriate second genome, such paralogous regions can still be identified.     

Detecting the undetectable: uncovering duplicated segments in Arabidopsis by comparison with rice

Genome analysis shows that large-scale gene duplications have occurred in fungi, animals and plants, creating genomic regions that show similarity in gene content and order. However, the high frequency of gene loss reduces colinearity resulting in duplicated regions that, in the extreme, no longer share homologous genes. Here, we show that by comparison with an appropriate second genome, such paralogous regions can still be identified.     

Analysis of 41 cancer cell lines reveals excessive allelic loss and novel mutations in the SIRT1 gene

SIRT1 is an evolutionarily conserved protein deacetylase that modulates stress response, cellular metabolism and aging in model organisms. While SIRT1 exerts beneficial effects in protecting against age-related diseases, the role of SIRT1 in cancer has been controversial. SIRT1 promotes cell survival by deacetylating, and thereby negatively regulating the activity of important tumor suppressors such as p53. In this regard, SIRT1 has been considered to be a potential oncogene, and SIRT1 inhibitors have been studied for possible anticancer therapeutic effects. In contrast, it has been shown that SIRT1 deficiency leads to increased genomic instability and tumorigenesis, and that overexpression of SIRT1 attenuates cancer formation in mice, suggesting it may also act as a tumor suppressor. Based on this evidence, SIRT1-activating molecules could act as candidate chemotherapeutic drugs. In order to gain insight into the role of SIRT1 in cancer, we performed a comprehensive resequencing analysis of the SIRT1 gene in 41 tumor cell lines and found an unusually excessive homozygosity, which was confirmed to be allelic loss by microsatellite analysis. Furthermore, we found two novel SIRT1 mutations (D739Y and R65_A72del) in addition to the known, rare non-synonymous variation resulting in I731V. In vitro assays using purified SIRT1 protein showed that these mutations do not alter SIRT1 deacetylase activity or telomerase activity, which was shown to be regulated by SIRT1. We conclude that allelic loss or mutations in the SIRT1 gene occur prevalently during tumorigenesis, supporting the assertion that SIRT1 may serve as a tumor suppressor.     

A hotspot of gene order rearrangement by tandem duplication and random loss in the vertebrate mitochondrial genome

Most reported examples of change in vertebrate mitochondrial (mt) gene order could be explained by a tandem duplication followed by random loss of redundant genes (tandem duplication-random loss [TDRL] model). Under this model of evolution, independent loss of genes arising from a single duplication in an ancestral species are predicted, and remnant pseudogenes expected, intermediate states that may remain in rearranged genomes. However, evidence for this is rare and largely scattered across vertebrate lineages. Here, we report new derived mt gene orders in the vertebrate "WANCY" region of four closely related caecilian amphibians. The novel arrangements found in this genomic region (one of them is convergent with the derived arrangement of marsupials), presence of pseudogenes, and positions of intergenic spacers fully satisfy predictions from the TDRL model. Our results, together with comparative data for the available vertebrate complete mt genomes, provide further evidence that the WANCY genomic region is a hotspot for gene order rearrangements and support the view that TDRL is the dominant mechanism of gene order rearrangement in vertebrate mt genomes. Convergent gene rearrangements are not unlikely in hotspots of gene order rearrangement by TDRL.     

Rapid identification of the Leptosphaeria maculans avirulence gene AvrLm2 using an intraspecific comparative genomics approach

Five avirulence genes from Leptosphaeria maculans, the causal agent of blackleg of canola (Brassica napus), have been identified previously through map‐based cloning. In this study, a comparative genomic approach was used to clone the previously mapped AvrLm2. Given the lack of a presence–absence gene polymorphism coincident with the AvrLm2 phenotype, 36 L. maculans isolates were resequenced and analysed for single‐nucleotide polymorphisms (SNPs) in predicted small secreted protein‐encoding genes present within the map interval. Three SNPs coincident with the AvrLm2 phenotype were identified within LmCys1, previously identified as a putative effector‐coding gene. Complementation of a virulent isolate with LmCys1, as the candidate AvrLm2 allele, restored the avirulent phenotype on Rlm2‐containing B. napus lines. AvrLm2 encodes a small cysteine‐rich protein with low similarity to other proteins in the public databases. Unlike other avirulence genes, AvrLm2 resides in a small GC island within an AT‐rich isochore of the genome, and was never found to be deleted completely in virulent isolates.     

Evolution of trace amine associated receptor (TAAR) gene family in vertebrates: lineage-specific expansions and degradations of a second class of vertebrate chemosensory receptors expressed in the olfactory epithelium

The trace amine-associated receptors (TAARs) form a specific family of G protein-coupled receptors in vertebrates. TAARs were initially considered neurotransmitter receptors, but recent study showed that mouse TAARs function as chemosensory receptors in the olfactory epithelium. To clarify the evolutionary dynamics of the TAAR gene family in vertebrates, near-complete repertoires of TAAR genes and pseudogenes were identified from the genomic assemblies of 4 teleost fishes (zebrafish, fugu, stickleback, and medaka), western clawed frogs, chickens, 3 mammals (humans, mice, and opossum), and sea lampreys. Database searches revealed that fishes had many putatively functional TAAR genes (13-109 genes), whereas relatively small numbers of TAAR genes (3-22 genes) were identified in tetrapods. Phylogenetic analysis of these genes indicated that the TAAR gene family was subdivided into 5 subfamilies that diverged before the divergence of ray-finned fishes and tetrapods. In tetrapods, virtually all TAAR genes were located in 1 specific region of their genomes as a gene cluster; however, in fishes, TAAR genes were scattered throughout more than 2 genomic locations. This possibly reflects a whole-genome duplication that occurred in the common ancestor of ray-finned fishes. Expression analysis of zebrafish and stickleback TAAR genes revealed that many TAARs in these fishes were expressed in the olfactory organ, suggesting the relatively high importance of TAARs as chemosensory receptors in fishes. A possible evolutionary history of the vertebrate TAAR gene family was inferred from the phylogenetic and comparative genomic analyses.     

Identification of five human novel genes associated with cell proliferation by cell-based screening from an expressed cDNA ORF library

The development of functional profiling technologies provides opportunity for high-throughput functional genomics studies. We describe a cell-based screening system to identify novel human genes associated with cell proliferation. The method integrates luciferase reporter gene activity, fluorescence stain, automated microscopy and cellular phenotype assays. We successfully used the system to screen 409 novel human genes cloned by our lab and found that 27 genes significantly up-regulated promoter-Renilla luciferase reporter plasmid (pRL) activity. Among them, five genes, TRAF3IP3, ZNF306, ZNF250, SGOL1, and ZNF434, were determined through morphological observation, calcein AM fluorescence stain, MTT assay and cell cycle analysis to be associated with cell proliferation. Furthermore, we showed that the gene TRAF3IP3, which initially was identified to specifically interact with TRAF3, stimulated cell growth by modulating the c-Jun N-terminal kinase (JNK) pathway, and RNAi of TRAF3IP3 confirmed that the effect was physiological and necessary. In summary, we integrated a rapid and efficient system for screening novel growth regulatory genes. Using the new screening system we identified five genes associated with cell proliferation for the first time.     

Evolution of the vertebrate twist family and synfunctionalization: a mechanism for differential gene loss through merging of expression domains

Twist genes are essential for embryonic development and are conserved from jellyfish to human. To study the vertebrate twist family and its evolution, the entire complement of twist genes was obtained for 9 representative species. Phylogenetic analysis showed that a single protochordate twist gene was duplicated at least twice before the teleost-tetrapod split to give rise to 3 ancestral genes, which were further duplicated or deleted, resulting in fluctuating number of twist paralogs in different vertebrate lineages. To find whether changes in gene copy number were associated with changes in gene function, embryonic expression patterns of twist orthologs were evaluated against the number of twist paralogs in different species. The results showed evidence for both neo- and subfunctionalization, and, in addition, for loss of an ancestral regulatory gene. For example, in Xenopus, twist2 was lost, but the twist1 paralog acquired, and therefore preserved, twist2 function. A general model is proposed to explain the data. In this process, termed synfunctionalization, one paralog acquires the expression domain(s) of another. The merging may lead to function shuffle. Alternatively, it may leave one paralog redundant and thus subject to deletion--while its function is retained by the surviving paralog(s). Synfunctionalization is a mechanism that, together with neo- and subfunctionalization, may work to establish equilibrium in the number of genes that regulate developmental processes; it may regulate the complexity of regulatory regions as well as gene copy number and therefore may play a role in evolution of gene function and the structure of genome.     

Genome sequence comparison reveals independent inactivation of the caspase-15 gene in different evolutionary lineages of mammals

We have recently demonstrated that placental mammalian species such as pig and dog express a novel proapoptotic protease, caspase-15, whereas mouse and humans lack this enzyme. Here we investigated the evolutionary fate of the caspase-15 gene in different mammalian lineages by analyzing whole-genome shotgun sequences of 30 mammalian species for the presence of caspase-15 orthologs. Caspase-15 gene sequences were found in representatives of all major mammalian clades except for the superorders Afrotheria (tenrec, rock hyrax, and elephant) and Euarchontoglires (rodents, rabbit, tree shrew, and primates), which either lacked any caspase-15-like sequences or contained mutated remnants of the caspase-15 gene. Polymerase chain reaction screenings confirmed the results of the database searches and showed that the caspase-15 gene is expressed not only in various placental mammals but also in the marsupial, Monodelphis domestica. The observed species distribution implies that caspase-15 has originated in an early ancestor of modern mammals and has been conserved, over more than 180 Myr, in marsupials and many placental mammals, whereas it was independently lost in 2 phylogenetically distant clades of placental mammals, that is, Afrotheria and Euarchontoglires. Our data suggest that the inactivation of the caspase-15 gene was not counteracted by, and may even have been driven by, evolutionary constraints in these clades, and therefore, caution against the uncritical use of gene absence for the inference of phylogenetic relationships.     

On the incidence of intron loss and gain in paralogous gene families

Understanding gene duplication and gene structure evolution are fundamental goals of molecular evolutionary biology. A previous study by Babenko et al. (2004. Prevalence of intron gain over intron loss in the evolution of paralogous gene families. Nucleic Acids Res. 32:3724-3733) employed Dollo parsimony to infer spliceosomal intron losses and gains in paralogous gene families and concluded that there was a general excess of gains over losses. This result contrasts with patterns in orthologous genes, in which most lineages show an excess of intron losses over gains, suggesting the possibility of fundamentally different modes of intron evolution between orthologous and paralogous genes. We further studied the data and found a low level of intron position conservation with outgroups, and this led to problems with using Dollo parsimony to analyze the data. Statistical reanalysis of the data suggests, instead, that intron losses have outnumbered intron gains in paralogous gene families.     

Direct inhibition of cytochrome c-induced caspase activation in vitro by Toxoplasma gondii reveals novel mechanisms of interference with host cell apoptosis

The intracellular parasite Toxoplasma gondii is known to inhibit apoptosis of its host cell. The molecular mechanisms of this interference are, however, not yet completely understood. We show here that viable parasites prominently inhibited the activation of caspase 3/7 induced by cytochrome c, dATP and dithiothreitol in cytosolic extracts of human-derived Jurkat leukemic T cells. In contrast, granzyme B-induced caspase activity was only slightly diminished. De novo protein biosynthesis by T. gondii was dispensable for the inhibition of cytochrome c-induced caspase activation. Furthermore, a complete parasite lysate or, more importantly, molecules released by extracellular parasites mediated the interaction with the caspase cascade. The cell-free system applied here is thus a valuable tool to study the interaction of T. gondii and possibly other intracellular pathogens with host cell apoptosis.     

Identification of a second egfr gene in Xiphophorus uncovers an expansion of the epidermal growth factor receptor family in fish

The epidermal growth factor receptor (EGFR) gives name to a family of receptors formed by four members in mammals (EGFR, ErbB2, ErbB3, and ErbB4). Members of this family can be activated to become potent oncogenes, and many human and animal tumors express qualitatively or quantitatively altered receptors from this group. We have isolated and characterized a second egfr gene in the melanoma model fish Xiphophorus. Both Xiphophorus egfra and egfrb duplicates are co-orthologs of the mammalian egfr gene. Database analysis showed that not only egfr but also erbB3 and erbB4 are present as duplicates in some fish species. They originated from ancient duplication events that might be consistent with the hypothesis of a fish-specific genome duplication. In Xiphophorus, the egfrb gene underwent a second duplication that generated the melanoma-inducing oncogene Xmrk. The study and comparison of some of the functional characteristics of both Xiphophorus EGF receptors, including expression profile, ligand-binding abilities, and intracellular signal transduction revealed that Xiphophorus Egfra not only shares common features with Egfrb and the human EGFR but also shows significant differences in its functional characteristics. The mechanism of maintenance of these duplicates remains to be clarified.     

Chromosomal localization of the chicken and mammalian orthologues of the orphan phosphatase PHOSPHO1 gene

PHOSPHO1 is a recently identified phosphatase expressed at high levels in the chicken growth plate and which may be involved in generating inorganic phosphate for skeletal matrix mineralization. Using a degenerate RT-PCR approach a fragment of human PHOSPHO1 was cloned. This enabled the identification of the human orthologue on HSA17q21, and the mouse orthologue on a region of MMU11 that exhibits conservation of synteny with HSA17q21. Chicken PHOSPHO1 was mapped by SSCP analysis to position 44 cM on GGA27, adjacent to the HOXB@ (44 cM) and COL1A1 (36 cM) loci. Comparison of genes on GGA27 with their orthologues on the preliminary draft of the human genome identifies regions of conserved synteny equivalent to 25 Mb on HSA17q21.2-23.3 and approximately 20 Mb on GGA27 in which the gene order appears to be conserved. Mapping of the PHOSPHO1 genes to regions of HSA17q21.3, MMU11 and GGA27 that exhibit conservation of synteny provides strong evidence that they are orthologous.     

Adaptive loss of an old duplicated gene during incipient speciation

To probe the role of natural selection in species origin, we performed a DNA polymorphism survey of the Drosophila melanogaster desaturase2 (ds2) locus. ds2 is responsible for a cuticular hydrocarbon difference between two behaviorally isolated races--Zimbabwe (Z) and Cosmopolitan (M). The ds2 allele prevalent in the Z populations is functional, while the allele from the M populations harbors a 16-bp deletion upstream of the gene which knocks out its expression. We find a signature of positive selection in the ds2 promoter, but not in the control gene, sas. This signature appears to be confined to the derived M population. We also find that the selection has been recent because the gene retains a signature of a selective sweep evidenced by the departure of Fay and Wu's H test from neutral expectation. We also find that ds2, as well as its duplicate pair ds1, has been maintained in the Drosophila genus for at least 40 Myr without any sign of adaptive change. Taken together with previous molecular genetic evidence, our results suggest that ds2 is one of the genes responsible for adaptive divergence of the Z and M races of D. melanogaster.