Comparative genomics using <Emphasis Type="Italic">Fugu</Emphasis> reveals insights into regulatory subfunctionalization

Background  

A major mechanism for the preservation of gene duplicates in the genome is thought to be mediated via loss or modification of cis-regulatory subfunctions between paralogs following duplication (a process known as regulatory subfunctionalization). Despite a number of gene expression studies that support this mechanism, no comprehensive analysis of regulatory subfunctionalization has been undertaken at the level of the distal cis-regulatory modules involved. We have exploited fish-mammal genomic alignments to identify and compare more than 800 conserved non-coding elements (CNEs) that associate with genes that have undergone fish-specific duplication and retention.     

The transmembrane domain of N –acetylglucosaminyltransferase I is the key determinant for its Golgi subcompartmentation

Golgi‐resident type–II membrane proteins are asymmetrically distributed across the Golgi stack. The intrinsic features of the protein that determine its subcompartment‐specific concentration are still largely unknown. Here, we used a series of chimeric proteins to investigate the contribution of the cytoplasmic, transmembrane and stem region of Nicotiana benthamiana N–acetylglucosaminyltransferase I (GnTI) for its cis/medial‐Golgi localization and for protein–protein interaction in the Golgi. The individual GnTI protein domains were replaced with those from the well‐known trans‐Golgi enzyme α2,6–sialyltransferase (ST) and transiently expressed in Nicotiana benthamiana. Using co‐localization analysis and N–glycan profiling, we show that the transmembrane domain of GnTI is the major determinant for its cis/medial‐Golgi localization. By contrast, the stem region of GnTI contributes predominately to homomeric and heteromeric protein complex formation. Importantly, in transgenic Arabidopsis thaliana, a chimeric GnTI variant with altered sub‐Golgi localization was not able to complement the GnTI‐dependent glycosylation defect. Our results suggest that sequence‐specific features in the transmembrane domain of GnTI account for its steady‐state distribution in the cis/medial‐Golgi in plants, which is a prerequisite for efficient N–glycan processing in vivo.     

Evolution of multiple phosphodiesterase isoforms in stickleback involved in cAMP signal transduction pathway

Background  

Duplicate genes are considered to have evolved through the partitioning of ancestral functions among duplicates (subfunctionalization) and/or the acquisition of novel functions from a beneficial mutation (neofunctionalization). Additionally, an increase in gene dosage resulting from duplication may also confer an advantageous effect, as has been suggested for histone, tRNA, and rRNA genes. Currently, there is little understanding of the effect of increased gene dosage on subcellular networks like signal transduction pathways. Addressing this issue may provide further insights into the evolution by gene duplication.     

Differential gene-expression patterns in genital fibroblasts of normal males and 46,XY females with androgen insensitivity syndrome: evidence for early programming involving the androgen receptor

Background  

Androgen insensitivity syndrome (AIS) comprises a range of phenotypes from male infertility to complete feminization. Most individuals with AIS carry germline mutations of the androgen receptor (AR) that interfere with or ablate its function. As genital fibroblasts retain expression of the AR in vitro, we used genital skin fibroblasts from normal males and 46,XY females with complete AIS due to known AR mutations to gain insights into the role of the AR in human genital differentiation.     

Loss of Androgen Receptor-Dependent Growth Suppression by Prostate Cancer Cells Can Occur Independently from Acquiring Oncogenic Addiction to Androgen Receptor Signaling

The conversion of androgen receptor (AR) signaling as a mechanism of growth suppression of normal prostate epithelial cells to that of growth stimulation in prostate cancer cells is often associated with AR mutation, amplification and over-expression. Thus, down-regulation of AR signaling is commonly therapeutic for prostate cancer. The E006AA cell line was established from a hormone naïve, localized prostate cancer. E006AA cells are genetically aneuploid and grow equally well when xenografted into either intact or castrated male NOG but not nude mice. These cells exhibit: 1) X chromosome duplication and AR gene amplification, although paradoxically not coupled with increased AR expression, and 2) somatic, dominant-negative Serine-599-Glycine loss-of-function mutation within the dimerization surface of the DNA binding domain of the AR gene. No effect on the growth of E006AA cells is observed using targeted knockdown of endogenous mutant AR, ectopic expression of wild-type AR, or treatment with androgens or anti-androgens. E006AA cells represent a prototype for a newly identified subtype of prostate cancer cells that exhibit a dominant-negative AR loss-of-function in a hormonally naïve patient. Such loss-of-function eliminates AR-mediated growth suppression normally induced by normal physiological levels of androgens, thus producing a selective growth advantage for these malignant cells in hormonally naïve patients. These data highlight that loss of AR-mediated growth suppression is an independent process, and that, without additional changes, is insufficient for acquiring oncogene addiction to AR signaling. Thus, patients with prostate cancer cells harboring such AR loss-of-function mutations will not benefit from aggressive hormone or anti-AR therapies even though they express AR protein.     

Three novel and two known androgen receptor gene mutations associated with androgen insensitivity syndrome in sex-reversed XY female patients

Molecular characterization of 23 cytogenetically confirmed XY females was attempted by screening coding regions of SRY and androgen receptor (AR) genes. Five of the index cases showed sequence variations in various exons of the AR gene: a deletion (n.1911delG) and substitutions n.1761G >A and n.1317C >T in exon 1; n.3510C >T transition in exon 6 and deletion mutation (n.3672delT) in exon 7. Four mutations identified here lead to the formation of truncated receptor protein, involving a substantial loss of AR functional domains which explains the phenotype in the subjects. The n.1761G >A substitution has been previously reported in cases with mild androgen insensitivity. Although the ligand-binding domain was considered as the mutational hot spot in AR gene, we report here 3/5 variations in the N-terminal domain emphasizing the significance of considering the N-terminal domain of AR as well for mutation screening. Our present observation also strengthens the role of AR gene and its direct association with AIS.     

Alternative Splicing and Subfunctionalization Generates Functional Diversity in Fungal Proteomes

Alternative splicing is commonly used by the Metazoa to generate more than one protein from a gene. However, such diversification of the proteome by alternative splicing is much rarer in fungi. We describe here an ancient fungal alternative splicing event in which these two proteins are generated from a single alternatively spliced ancestral SKI7/HBS1 gene retained in many species in both the Ascomycota and Basidiomycota. While the ability to express two proteins from a single SKI7/HBS1 gene is conserved in many fungi, the exact mechanism by which they achieve this varies. The alternative splicing was lost in Saccharomyces cerevisiae following the whole-genome duplication event as these two genes subfunctionalized into the present functionally distinct HBS1 and SKI7 genes. When expressed in yeast, the single gene from Lachancea kluyveri generates two functionally distinct proteins. Expression of one of these proteins complements hbs1, but not ski7 mutations, while the other protein complements ski7, but not hbs1. This is the first known case of subfunctionalization by loss of alternative splicing in yeast. By coincidence, the ancestral alternatively spliced gene was also duplicated in Schizosaccharomyces pombe with subsequent subfunctionalization and loss of splicing. Similar subfunctionalization by loss of alternative splicing in fungi also explains the presence of two PTC7 genes in the budding yeast Tetrapisispora blattae, suggesting that this is a common mechanism to preserve duplicate alternatively spliced genes.     

Analysing gene function after duplication

After gene duplication, mutations cause the gene copies to diverge. The classical model predicts that these mutations will generally lead to the loss of function of one gene copy; rarely, new functions will be created and both duplicate genes are conserved. In contrast, under the subfunctionalization model both duplicates are preserved due to the partition of different functions between the duplicates. A recent study provides support for the subfunctionalization model, identifying several expressed gene duplicates common to humans and mice that contain regions conserved in one duplicate but variable in the other (and vice versa). We discuss both the methodology used in this study and also how gene phylogeny may lead to additional evidence for the importance of subfunctionalization in the evolution of new genes.     

Gene duplication and the evolution of phenotypic diversity in insect societies

Gene duplication is an important evolutionary process thought to facilitate the evolution of phenotypic diversity. We investigated if gene duplication was associated with the evolution of phenotypic differences in a highly social insect, the honeybee Apis mellifera. We hypothesized that the genetic redundancy provided by gene duplication could promote the evolution of social and sexual phenotypes associated with advanced societies. We found a positive correlation between sociality and rate of gene duplications across the Apoidea, indicating that gene duplication may be associated with sociality. We also discovered that genes showing biased expression between A. mellifera alternative phenotypes tended to be found more frequently than expected among duplicated genes than singletons. Moreover, duplicated genes had higher levels of caste‐, sex‐, behavior‐, and tissue‐biased expression compared to singletons, as expected if gene duplication facilitated phenotypic differentiation. We also found that duplicated genes were maintained in the A. mellifera genome through the processes of conservation, neofunctionalization, and specialization, but not subfunctionalization. Overall, we conclude that gene duplication may have facilitated the evolution of social and sexual phenotypes, as well as tissue differentiation. Thus this study further supports the idea that gene duplication allows species to evolve an increased range of phenotypic diversity.     

Molecular evolution of the COX7A gene family in primates.

COX VIIa is one of 10 nuclear-encoded subunits of the COX holoenzyme, and one of three that have isoforms with tissue-specific differences in expression. Analysis of nucleotide substitution rates revealed an accelerated rate of nonsynonymous substitutions relative to that of synonymous substitutions for the heart isoform gene (COX7AH) in six primate lineages. Rate accelerations have been noted for four other COX-related genes in this time period, suggesting that the COX holoenzyme has experienced an episode of adaptive evolution. A third member of the gene family, COX7AR, has recently been described. Although its function is currently unknown, low nonsynonymous substitution/synonymous substitution (N/S) ratios in mammalian evolution suggest that COX7AR is of functional importance. When the COX7A isoforms were divided into domains, examination of nucleotide substitution rates suggested that mitochondrial targeting residues experienced an accelerated nonsynonymous substitution rate in the period following gene duplication. In contrast, paralogous comparisons of the targeting residues of each isoform show they have been relatively conserved in mammalian evolution. This pattern is consistent with the evolution of tissue-specific function.     

Characterization of a novel amphiregulin-related molecule in 12-O-tetradecanoylphorbol-13-acetate-treated breast cancer cells

Amphiregulin (AR) can be induced at the mRNA level by 17-β-estradiol (E₂) or the phorbol ester tumor promoter 12-O-tetradecanoylphorbol-13-acetate (TPA). This study compares the effects of TPA and E₂ on the regulation of processing of AR isoforms and on subcellular localization in human MCF-7 breast cancer cells. AR was localized in the nucleus of MCF-7 cells after E₂ treatment, whereas it was predominantly secreted after TPA treatment. AR isoforms of 28, 18, and 10 kDa and an additional species of approximately 55–60 kDa were detected in the cellular conditioned media after TPA stimulation. Expression of this unusual AR isoform was inhibited by protein kinase C (PKC) inhibitors such as bryostatin or H-7. The biochemical properties of this isoform are consistent with it being an N-linked glycosylated form of the AR precursor that contains unprocessed mannose residues. The size of this large isoform is reduced to approximately 40 kDa after treating the TPA-induced MCF-7 cells with tunicamycin or treating the conditioned media of such cells with N-glycosidase F or with endoglycosidase H. Moreover, this isoform is able to bind several lectins with specificity for mannose residues. The 55–60 kDa glycosylated AR isoform, like lower Mr AR isoforms, is able to bind to heparin and to stimulate the growth of MCF-10A cells by interacting with the EGF receptor. These data suggest that TPA activation of PKC may be involved in post-translational modifications of AR, such as glycosylation, and in alteration of its subcellular routing to predominantly a secretory pathway. © 1996 Wiley-Liss, Inc.     

Adenosine A2A receptor and vascular response: role of soluble epoxide hydrolase,adenosine A1 receptor and angiotensin-II

Previously, we have reported that the coronary reactive hyperemic response was reduced in adenosine A_2A receptor-null (A_2AAR^?/?) mice, and it was reversed by the soluble epoxide hydrolase (sEH) inhibitor. However, it is unknown in aortic vascular response, therefore, we hypothesized that A_2AAR-gene deletion in mice (A_2AAR^?/?) affects adenosine-induced vascular response by increase in sEH and adenosine A₁ receptor (A₁AR) activities. A_2AAR^?/? mice showed an increase in sEH, A_I AR and CYP450-4A protein expression but decrease in CYP450-2C compared to C57Bl/6 mice. NECA (adenosine-analog) and CCPA (adenosine A₁ receptor-agonist)-induced dose-dependent vascular response was tested with t-AUCB (sEH-inhibitor) and angiotensin-II (Ang-II) in A_2AAR^?/? vs. C57Bl/6 mice. In A_2AAR^?/?, NECA and CCPA-induced increase in dose-dependent vasoconstriction compared to C57Bl/6 mice. However, NECA and CCPA-induced dose-dependent vascular contraction in A_2AAR^?/? was reduced by t-AUCB with NECA. Similarly, dose-dependent vascular contraction in A_2AAR^?/? was reduced by t-AUCB with CCPA. In addition, Ang-II enhanced NECA and CCPA-induced dose-dependent vascular contraction in A_2AAR^?/? with NECA. Similarly, the dose-dependent vascular contraction in A_2AAR^?/? was also enhanced by Ang-II with CCPA. Further, t-AUCB reduced Ang-II-enhanced NECA and CCPA-induced dose-dependent vascular contraction in A_2AAR^?/? mice. Our data suggest that the dose-dependent vascular contraction in A_2AAR^?/? mice depends on increase in sEH, A₁AR and CYP4A protein expression.

     

Two alloalleles of <Emphasis Type="Italic">Xenopus laevis hairy2</Emphasis> gene—evolution of duplicated gene function from a developmental perspective

Gene duplication is a fundamental source of a new gene in the process of evolution. A duplicated gene is able to accept many kinds of mutations that could lead to loss of function or novel phenotypic diversity. Alternatively, the duplicated genes complementarily lose part of their functions to play original roles as a set of genes, a process called subfunctionalization. Pseudotetraploid frog Xenopus laevis has four sets of genes, and it is generally thought that the alloalleles in X. laevis have mutually indistinguishable functions. In this paper, we report differences and similarities between Xhairy2a and Xhairy2b in the neural crest, floor plate, and prechordal plate. Knockdown studies showed that Xhairy2a seems not to function in the neural crest, although both of them are required in the floor plate and the prechordal plate. Temporal expression pattern analysis revealed that Xhairy2a is a maternal factor having lower zygotic expression than Xhairy2b, while Xhairy2b is not loaded in the egg but has high zygotic expression. Spatial expression pattern analysis demonstrated that future floor plate expression is shared by both alloalleles, but Xhairy2b expression in the neural crest is much higher than Xhairy2a expression, consistent with the results of individual knockdown experiments. Therefore, our data suggest that subfunctionalization occurs in Xhairy2. Electronic supplementary material The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

Complex arrangement of genes within a 220-kb region of double-duplicated DNA on human 2q37.1

Gene duplication events are followed by divergence of initially identical gene copies, due to the subsequent accumulation of mutations. These mutations tend to be degenerative and may lead to either nonfunctionalization or subfunctionalization of the gene copies. Here we report the molecular characterization of a 220-kb genomic DNA fragment from human 2q37.1, in which a double duplication and a partial triplication event has taken place. As a result, this region contains four copies of alkaline phosphatase (P), four copies of the ECEL1 gene (X), two copies of a newly identified gene (N), and two copies of a cholinergic receptor subunit (R), in the order N-P-X-P-X-P-X-N-P-X-R-R. While three of the four ECEL1 copies, one copy of the phosphatase gene and one copy of the newly identified gene have lost their function, three phosphatase gene copies and the two receptor subunits are still functionally active and thus may provide an example for subfunctionalization of duplicated genes.     

The limits of subfunctionalization

Background  

The duplication-degeneration-complementation (DDC) model has been proposed as an explanation for the unexpectedly high retention of duplicate genes. The hypothesis proposes that, following gene duplication, the two gene copies degenerate to perform complementary functions that jointly match that of the single ancestral gene, a process also known as subfunctionalization. We distinguish between subfunctionalization at the regulatory level and at the product level (e.g within temporal or spatial expression domains).     

Molecular and genomic studies of IMMP2L and mutation screening in autism and Tourette syndrome

We recently reported the disruption of the inner mitochondrial membrane peptidase 2-like (IMMP2L) gene by a chromosomal breakpoint in a patient with Gilles de la Tourette syndrome (GTS). In the present study we sought to identify genetic variation in IMMP2L, which, through alteration of protein function or level of expression might contribute to the manifestation of GTS. We screened 39 GTS patients, and, due to the localization of IMMP2L in the critical region for the autistic disorder (AD) locus on chromosome 7q (AUTS1), 95 multiplex AD families; however, no coding mutations were found in either GTS or AD patients. In addition, no parental-specific expression of IMMP2L was detected in somatic cell hybrids containing human chromosome 7 and human cell lines carrying a maternal uniparental disomy for chromosome 7 (mUPD7). Despite the fact that no deleterious mutations in IMMPL2 (other than the inverted duplication identified previously) were identified in either GTS or AD, this gene cannot be excluded as a possible rare cause of either disorder.