Expression diversity and evolutionary dynamics of rice duplicate genes

Duplicate genes are believed to be a major source of new gene functions over evolutionary time. In order to evaluate the evolutionary dynamics of rice duplicate genes, formed principally by paleoployploidization prior to the speciation of the Poaceae family, we have employed a public microarray dataset including 155 gene expression omnibus sample plates and bioinformatics tools. At least 57.4% of old ~70 million years ago (MYA) duplicate gene pairs exhibit divergences in expression over the given experimental set, whereas at least 50.9% of young ~7.7-MYA duplicate gene pairs were shown to be divergent. When grouping the rice duplicate genes according to functional categories, we noted a striking and significant enrichment of divergent duplicate metabolism-associated genes, as compared to that observed in non-divergent duplicate genes. While both non-synonymous substitution (Ka) and synonymous substitution (Ks) values between non- and divergent duplicate gene pairs evidenced significant differences, the Ka/Ks values between them exhibited no significant differences. Interestingly, the average numbers of conserved motifs of the duplicate gene pairs revealed a pattern of decline along with an increase in expression diversity, partially supporting the subfunctionalization model with degenerative complementation in regulatory motifs. Duplicate gene pairs with high local similarity (HLS) segments, which might be formed via conversion between rice paleologs, evidenced higher expression correlations than were observed in the gene pairs without the HLS segments; this probably resulted in an increased likelihood of gene conversion in promoters of the gene pairs harboring HLS segments. More than 60% of the rice gene families exhibited similar high expression diversity between members as compared to that of randomly selected gene pairs. These findings are likely reflective of the evolutionary dynamics of rice duplicate genes for gene retention. An erratum to this article can be found at     

The yan gene is highly conserved in Drosophila and its expression suggests a complex role throughout development

 Competence for cell fate determination and cellular differentiation is under tight control of regulatory genes. Yan, a nuclear target of receptor tyrosine kinase (RTK) signaling, is an E twenty six (ETS) DNA-binding protein that functions as a negative regulator of cell differentiation and proliferation in Drosophila. Most members of RTK signaling pathways are highly conserved through evolution, yet no yan orthologues have been identified to date in vertebrates. To investigate the degree of yan conservation during evolution, we have characterized a yan homologue from a sibling species of D. melanogaster, D. virilis. Our results show that the organization, primary structure and expression pattern of yan are highly conserved. Both genes span over 20 kb and contain four exons with introns at identical positions. The areas with highest amino acid similarity include the Pointed and ETS domain but there are other discrete regions with a high degree of similarity. Phylogenetic analysis reveals that yan’s closest relative is the human tel gene, a negative regulator of differentiation in hematopoetic precursors. In both species, Yan is dynamically expressed beginning as early as stage 4/5 and persisting throughout embryogenesis. In third instar larvae, Yan is expressed in and behind the morphogenetic furrow of the eye imaginal disc as well as in the laminar precursor cells of the brain. Ovarian follicle cells also contain Yan protein. Conservation of the structure and expression patterns of yan genes strongly suggests that regulatory mechanisms for their expression are also conserved in these two species. Received: 3 November 1998 / Accepted: 9 December 1998     

An IQGAP-related gene is activated during tentacle formation in the simple metazoan Hydra

Differentiation of body column epithelial cells into tentacle epithelial cells in Hydra is accompanied by changes in both cell shape and cell-cell contact. The molecular mechanism by which epithelial cells acquire tentacle cell characteristics is unknown. Here we report that expression of a Hydra homologue of the mammalian IQGAP1 protein is strongly upregulated during tentacle formation. Like mammalian IQGAP, Hydra IQGAP1 contains an N-terminal calponin-homology domain, IQ repeats and a conserved C terminus. In adult polyps a high level of Hydra IQGAP1 mRNA is detected at the basis of tentacles. Consistent with a role in tentacle formation, IQGAP1 expression is activated during head regeneration and budding at a time when tentacles are emerging. The observations support the previous hypothesis that IQGAP proteins are involved in cytoskeletal as well as cell-cell contact rearrangements. Received: 25 January 2000 / Accepted: 2 May 2000     

Molecular and biological characterization of a zonula occludens-1 homologue in Hydra vulgaris, named HZO-1

Zonula occludens-1 (ZO-1) is one of the earliest identified molecular components of tight junctions. Sequence analysis has placed ZO-1 into the broader membrane-associated guanylate kinase (MAGUK) protein family that contains such diverse members as postsynaptic density 95 (PSD-95), Drosophila discs large tumor suppressor gene product (dlg-A), p55, and TamA. Studies in both vertebrates and invertebrates have established that the MAGUK family is involved in a wide variety of cellular functions. These functions involve the regulation of such cellular processes as: (1) tight junction formation, (2) cell proliferation, (3) cell differentiation, and (4) neuronal synapse transmission. Extending these studies, we report the presence of a ZO-1 homologue in Hydra vulgaris, a member of the Cnidaria, the second oldest phylum of the animal kingdom. Hydra ZO-1 (HZO-1) is encoded by a single messenger RNA (mRNA) of approximately 6.0 kb that contains an open reading frame of 5,085 bp. The 191 kDa predicted protein consists of a characteristic MAGUK domain structure, including three PSD-95/SAP90, discs-large, ZO-1 (PDZ) domains, a src homology (SH₃) domain, and a guanylate kinase (GUK) domain. Western blot analysis using an antibody generated from a synthetic peptide designed from the HZO-1 sequence confirmed the presence of a Hydra protein of the appropriate mass. While whole mount in situ hybridization determined that HZO-1 mRNA was expressed along the entire longitudinal axis of Hydra, cross-sectional analysis established that HZO-1 mRNA expression was restricted to the ectoderm or outer cell layer of the organism’s epithelial bilayer. Consistent with this mRNA expression pattern, immunofluorescence studies localized HZO-1 protein to the apical plasma membrane of ectodermal cells. It is unclear what role HZ0-1 has in the cellular physiology of Hydra; however, immunolocalization studies indicate a conserved plasma membrane-associated function(s), as reported for its counterparts in other invertebrate and vertebrate species. These studies establish that the MAGUK family of proteins with a membrane-associated function arose early during metazoan evolution, even before the divergence of protostomes and deuterostomes. Received: 11 May 2000 / Accepted: 26 July 2000     

Morphogenesis of the atrichous isorhiza, a type of nematocyst, in Hydra observed with a monoclonal antibody

 We produced a monoclonal antibody, AE03, which recognized mucous granules in the basal disk gland cells in Hydra and the secreted mucus with which they stick onto substrate. AE03 also recognized atrichous isorhizas, one of the four types of nematocyst present in tentacles, and their nematoblasts present in the body column. With this monoclonal antibody, we could observe the detailed morphogenesis of the atrichous isorhiza from the beginning of its formation. The elongation and invagination processes of external tubes and correspondence between the external tubes and the thread of discharged nematocysts were confirmed. Received: 20 May 1997 / Accepted: 21 September 1997     

Response to insulin and the expression pattern of a gene encoding an insulin receptor homologue suggest a role for an insulin-like molecule in regulating growth and patterning in Hydra

 A gene encoding a receptor protein-tyrosine kinase closely related to the vertebrate insulin receptor has been identified in the Cnidarian Hydra vulgaris. The gene is expressed in both epithelial layers of the adult polyp. A particularly high level of expression is seen in the ectoderm of the proximal portions of the tentacles and in a ring of ectodermal cells at the border between the foot basal disk and body column. The expression pattern of the gene in asexual buds is dynamic; expression is high throughout the newly emerging bud but the area of high expression becomes restricted to the apex as the bud lengthens. When the bud begins hypostome and tentacle formation, a high level of expression appears at the bases of the emerging tentacles. Finally, a ring of high expression appears just above the foot of the bud, completing the pattern seen in the adult polyp. The presence of this receptor and its pattern of expression suggested that an endogenous molecule related to insulin plays a role in regulating cell division in the body column and in differentiation of the tentacle and foot cells in Hydra, with the switch between the two being determined by the level of the receptor. Treatment of Hydra polyps with mammalian insulin caused an increase in the number of ectodermal and endodermal cells undergoing DNA synthesis. Received: 19 April 1996 / Accepted: 5 July 1996     

The level of expression of a protein kinase C gene may be an important component of the patterning process in Hydra

 Several studies have provided strong, but indirect evidence that signalling through pathways involving protein kinase C (PKC) plays an important role in morphogenesis and patterning in Hydra. We have cloned a gene (HvPKC2) from Hydra vulgaris which encodes a member of the nPKC subfamily. In adult polyps, HvPKC2 is expressed at high levels in two locations, the endoderm of the foot and the endoderm of the hypostomal tip. Increased expression of HvPKC2 is an early event during head and foot regeneration, with the rise in expression being restricted to the endodermal cells underlying the regenerating ends. No upregulation is observed if regenerates are cut too close to the head to form a foot. Elevated expression of HvPKC2 is also observed in the endoderm underlying lithium-induced ectopic feet. A dynamic and complex pattern of expression is seen in developing buds. Regeneration of either head or foot is accompanied by an increase in the amount of PKC in both soluble and particulate fractions. An increase in the fraction of PKC activity which is membrane-bound is specifically associated with head regeneration. Taken together these data suggest that patterning of the head and foot in Hydra is controlled in part by the level of HvPKC2 expression, whilst head formation is accompanied by an in vivo activation of both calcium-dependent and independent PKC isoforms. Received: 10 July 1997 / Accepted: 8 November 1997     

In the multiheaded strain (mh-1) of Hydra magnipapillata the ectodermal epithelial cells are responsible for the formation of additional heads and the endodermal epithelial cells for the reduced ability to regenerate a foot

Hydra consist of three self-renewing cell lineages: the ectodermal epithelial, endodermal epithelial and interstitial cell lineages. The role of these cell lineages in head formation and foot regeneration in Hydra magnipapillata was studied by comparing the multiheaded strain mh-1 with the wild-type. Adult polyps of this strain show a reduced ability to regenerate a foot in the apical body half several days before additional heads are formed there. Cell lineage chimeras were produced, and it was found that in mh-1, the ectodermal epithelial cell lineage is responsible for the formation of additional heads, whereas the endodermal epithelial cell lineage and, to a lesser extent, the derivatives of the interstitial cell lineage, are responsible for the reduced ability of foot regeneration.     

The mapping of phytochrome genes and photomorphogenic mutants of tomato

The map positions of five previously described phytochrome genes have been determined in tomato (Lycopersicon esculentum Mill.) The position of the yg-2 gene on chromosome 12 has been confirmed and the classical map revised. The position of the phytochrome A (phy A)-deficient fri mutants has been refined by revising the classical map of chromosome 10. The position of the PhyA gene is indistinguishable from that of the fri locus. The putative phyB1-deficient tri mutants were mapped by classical and RFLP analysis to chromosome 1. The PhyB1 gene, as predicted, was located at the same position. Several mutants with the high pigment (hp) phenotype, which exaggerates phytochrome responses, have been reported. Allelism tests confirmed that the hp-2 mutant is not allelic to other previously described hp (proposed here to be called hp-1) mutants and a second stronger hp-2 allele (hp-2 ^j ) was identified. The hp-2 gene was mapped to the classical, as well as the RFLP, map of chromosome 1. Received: 24 May 1996 / Accepted: 14 June 1996     

Comparative analysis of essential genes and nonessential genes in Escherichia coli K12

Genes can be classified as essential or nonessential based on their indispensability for a living organism. Previous researches have suggested that essential genes evolve more slowly than nonessential genes and the impact of gene dispensability on a gene’s evolutionary rate is not as strong as expected. However, findings have not been consistent and evidence is controversial regarding the relationship between the gene indispensability and the rate of gene evolution. Understanding how different classes of genes evolve is essential for a full understanding of evolutionary biology, and may have medical relevance in the design of new antibacterial agents. We therefore performed an investigation into the properties of essential and nonessential genes. Analysis of evolutionary conservation, protein length distribution and amino acid usage between essential and nonessential genes in Escherichia coli K12 demonstrated that essential genes are relatively preserved throughout the bacterial kingdom when compared to nonessential genes. Furthermore, results show that essential genes, compared to nonessential genes, have a significantly higher proportion of large (>534 amino acids) and small proteins (<139 amino acids) relative to medium-sized proteins. The pattern of amino acids usage shows a similar trend for essential and nonessential genes, although some notable exceptions are observed. These findings help to clarify our understanding of the evolutionary mechanisms of essential and nonessential genes, relevant to the study of mutagenesis and possibly allowing prediction of gene properties in other poorly understood organisms.     

Mitochondrial DNA of Hydra attenuata (Cnidaria): A Sequence That Includes an End of One Linear Molecule and the Genes for l-rRNA, tRNAf-Met, tRNATrp, COII, and ATPase8

The 3231-nucleotide-pair (ntp) sequence of one end of one of the two linear mitochondrial (mt) DNA molecules of Hydra attenuata (phylum Cnidaria, class Hydrozoa, order Anthomedusae) has been determined. This segment contains complete genes for tRNA^f-Met, l-rRNA, tRNA^Trp, subunit 2 of cytochrome c oxidase (COII), subunit 8 of ATP synthetase (ATPase8), and the 5′ 136 ntp of ATPase6. These genes are arranged in the order given and are transcribed from the same strand of the molecule. As in two other cnidarians, the hexacorallian anthozoan Metridium senile and the octocorallian anthozoan Sarcophyton glaucum, the mt-genetic code of H. attenuata is near standard. The only modification appears to be that TGA specifies tryptophan rather than termination. Also as in M. senile and S. glaucum, the encoded H. attenuata mt-tRNA^f-Met has primary and secondary structural features resembling those of Escherichia coli initiator tRNA^t-Met. As the encoded mt-tRNA^Trp cannot be folded into a totally orthodox secondary structure, two alternative forms are suggested. The encoded H. attenuata mt-l-rRNA is 1738 nt, which is 451 nt shorter than the M. senile mt-l-rRNA. Comparisons of secondary structure models of these two mt-l-rRNAs indicate that most of the size difference results from loss of nucleotides in the H. attenuata molecule at a minimum of 46 locations, which includes elimination of six distinct helical elements. Received: 9 March 2000 / Accepted: 24 July 2000     

Discovery of genes expressed in Hydra embryogenesis

Hydra's remarkable capacity to regenerate, to proliferate asexually by budding, and to form a pattern de novo from aggregates allows studying complex cellular and molecular processes typical for embryonic development. The underlying assumption is that patterning in adult hydra tissue relies on factors and genes which are active also during early embryogenesis. Previously, we reported that in Hydra the timing of expression of conserved regulatory genes, known to be involved in adult patterning, differs greatly in adults and embryos (Fr?bius, A.C., Genikhovich, G., Kürn, U., Anton-Erxleben, F. and Bosch, T.C.G., 2003. Expression of developmental genes during early embryogenesis of Hydra. Dev. Genes Evol. 213, 445-455). Here, we describe an unbiased screening strategy to identify genes that are relevant to Hydra vulgaris embryogenesis. The approach yielded two sets of differentially expressed genes: one set was expressed exclusively or nearly exclusively in the embryos, while the second set was upregulated in embryos in comparison to adult polyps. Many of the genes identified in hydra embryos had no matches in the database. Among the conserved genes upregulated in embryos is the Hydra orthologue of Embryonic Ectoderm Development (HyEED). The expression pattern of HyEED in developing embryos suggests that interstitial stem cells in Hydra originate in the endoderm. Importantly, the observations uncover previously unknown differences in genes expressed by embryos and polyps and indicate that not only the timing of expression of developmental genes but also the genetic context is different in Hydra embryos compared to adults.     

MHC class I genes of the channel catfish: sequence analysis and expression

 Four cDNAs encoding the major histocompatibility complex (MHC) class I α chain were isolated from a channel catfish clonal B-cell cDNA library. Sequence analysis suggests these cDNAs represent three different MHC class I loci. All cDNAs encoded conserved residues characteristic of the MHC class I α chain: namely, those involved in peptide binding, salt bridges, disulfide bond formation, and glycosylation. Southern blot analyses of individual outbred and second-generation gynogenetic fish indicated the existence of both polygenic and polymorphic loci. Northern blot studies demonstrated that catfish B, T, and macrophage cell lines transcribed markedly higher levels of class I α and β₂-microglobulin (β₂m) mRNA than fibroblast cell lines. In addition, immunoprecipitation data showed that a 41 000 M _r glycoprotein (presumably class I α) was associated with β₂m on the surface of catfish B cells. This latter finding is the first direct evidence for the cell surface association of β₂m with the MHC class I α chain on teleost cells and supports the notion that functional MHC class I proteins exist in teleosts. Received: 25 March 1998 / Revised: 28 July 1998