Protein-coding genes are epigenetically regulated in Arabidopsis polyploids

The fate of redundant genes resulting from genome duplication is poorly understood. Previous studies indicated that ribosomal RNA genes from one parental origin are epigenetically silenced during interspecific hybridization or polyploidization. Regulatory mechanisms for protein-coding genes in polyploid genomes are unknown, partly because of difficulty in studying expression patterns of homologous genes. Here we apply amplified fragment length polymorphism (AFLP)–cDNA display to perform a genome-wide screen for orthologous genes silenced in Arabidopsis suecica, an allotetraploid derived from Arabidopsis thaliana and Cardaminopsis arenosa. We identified ten genes that are silenced from either A. thaliana or C. arenosa origin in A. suecica and located in four of the five A. thaliana chromosomes. These genes represent a variety of RNA and predicted proteins including four transcription factors such as TCP3. The silenced genes in the vicinity of TCP3 are hypermethylated and reactivated by blocking DNA methylation, suggesting epigenetic regulation is involved in the expression of orthologous genes in polyploid genomes. Compared with classic genetic mutations, epigenetic regulation may be advantageous for selection and adaptation of polyploid species during evolution and development.

Horizontal gene transfer among genomes: The complexity hypothesis

Increasingly, studies of genes and genomes are indicating that considerable horizontal transfer has occurred between prokaryotes. Extensive horizontal transfer has occurred for operational genes (those involved in housekeeping), whereas informational genes (those involved in transcription, translation, and related processes) are seldomly horizontally transferred. Through phylogenetic analysis of six complete prokaryotic genomes and the identification of 312 sets of orthologous genes present in all six genomes, we tested two theories describing the temporal flow of horizontal transfer. We show that operational genes have been horizontally transferred continuously since the divergence of the prokaryotes, rather than having been exchanged in one, or a few, massive events that occurred early in the evolution of prokaryotes. In agreement with earlier studies, we found that differences in rates of evolution between operational and informational genes are minimal, suggesting that factors other than rate of evolution are responsible for the observed differences in horizontal transfer. We propose that a major factor in the more frequent horizontal transfer of operational genes is that informational genes are typically members of large, complex systems, whereas operational genes are not, thereby making horizontal transfer of informational gene products less probable (the complexity hypothesis).

Identification of a human nuclear receptor defines a new signaling pathway for CYP3A induction

Nuclear receptors regulate metabolic pathways in response to changes in the environment by appropriate alterations in gene expression of key metabolic enzymes. Here, a computational search approach based on iteratively built hidden Markov models of nuclear receptors was used to identify a human nuclear receptor, termed hPAR, that is expressed in liver and intestines. hPAR was found to be efficiently activated by pregnanes and by clinically used drugs including rifampicin, an antibiotic known to selectively induce human but not murine CYP3A expression. The CYP3A drug-metabolizing enzymes are expressed in gut and liver in response to environmental chemicals and clinically used drugs. Interestingly, hPAR is not activated by pregnenolone 16α-carbonitrile, which is a potent inducer of murine CYP3A genes and an activator of the mouse receptor PXR.1. Furthermore, hPAR was found to bind to and trans-activate through a conserved regulatory sequence present in human but not murine CYP3A genes. These results provide evidence that hPAR and PXR.1 may represent orthologous genes from different species that have evolved to regulate overlapping target genes in response to pharmacologically distinct CYP3A activators, and have potential implications for the in vitro identification of drug interactions important to humans.

Measuring genome evolution

The determination of complete genome sequences provides us with an opportunity to describe and analyze evolution at the comprehensive level of genomes. Here we compare nine genomes with respect to their protein coding genes at two levels: (i) we compare genomes as “bags of genes” and measure the fraction of orthologs shared between genomes and (ii) we quantify correlations between genes with respect to their relative positions in genomes. Distances between the genomes are related to their divergence times, measured as the number of amino acid substitutions per site in a set of 34 orthologous genes that are shared among all the genomes compared. We establish a hierarchy of rates at which genomes have changed during evolution. Protein sequence identity is the most conserved, followed by the complement of genes within the genome. Next is the degree of conservation of the order of genes, whereas gene regulation appears to evolve at the highest rate. Finally, we show that some genomes are more highly organized than others: they show a higher degree of the clustering of genes that have orthologs in other genomes.

A novel application of gene arrays: Escherichia coli array provides insight into the biology of the obligate endosymbiont of tsetse flies

Symbiotic associations with microorganisms are pivotal in many insects. Yet, the functional roles of obligate symbionts have been difficult to study because it has not been possible to cultivate these organisms in vitro. The medically important tsetse fly (Diptera: Glossinidae) relies on its obligate endosymbiont, Wigglesworthia glossinidia, a member of the Enterobacteriaceae, closely related to Escherichia coli, for fertility and possibly nutrition. We show here that the intracellular Wigglesworthia has a reduced genome size smaller than 770 kb. In an attempt to understand the composition of its genome, we used the gene arrays developed for E. coli. We were able to identify 650 orthologous genes in Wigglesworthia corresponding to ≈85% of its genome. The arrays were also applied for expression analysis using Wigglesworthia cDNA and 61 gene products were detected, presumably coding for some of its most abundant products. Overall, genes involved in cell processes, DNA replication, transcription, and translation were found largely retained in the small genome of Wigglesworthia. In addition, genes coding for transport proteins, chaperones, biosynthesis of cofactors, and some amino acids were found to comprise a significant portion, suggesting an important role for these proteins in its symbiotic life. Based on its expression profile, we predict that Wigglesworthia may be a facultative anaerobic organism that utilizes ammonia as its major source of nitrogen. We present an application of E. coli gene arrays to obtain broad genome information for a closely related organism in the absence of complete genome sequence data.

Comparative mapping of the human 22q11 chromosomal region and the orthologous region in mice reveals complex changes in gene organization

The region of human chromosome 22q11 is prone to rearrangements. The resulting chromosomal abnormalities are involved in Velo-cardio-facial and DiGeorge syndromes (VCFS and DGS) (deletions), “cat eye” syndrome (duplications), and certain types of tumors (translocations). As a prelude to the development of mouse models for VCFS/DGS by generating targeted deletions in the mouse genome, we examined the organization of genes from human chromosome 22q11 in the mouse. Using genetic linkage analysis and detailed physical mapping, we show that genes from a relatively small region of human 22q11 are distributed on three mouse chromosomes (MMU6, MMU10, and MMU16). Furthermore, although the region corresponding to about 2.5 megabases of the VCFS/DGS critical region is located on mouse chromosome 16, the relative organization of the region is quite different from that in humans. Our results show that the instability of the 22q11 region is not restricted to humans but may have been present throughout evolution. The results also underscore the importance of detailed comparative mapping of genes in mice and humans as a prerequisite for the development of mouse models of human diseases involving chromosomal rearrangements.

Comparative analysis of the gene-dense ACHE/TFR2 region on human chromosome 7q22 with the orthologous region on mouse chromosome 5

Chromosome 7q22 has been the focus of many cytogenetic and molecular studies aimed at delineating regions commonly deleted in myeloid leukemias and myelodysplastic syndromes. We have compared a gene-dense, GC-rich sub-region of 7q22 with the orthologous region on mouse chromosome 5. A physical map of 640 kb of genomic DNA from mouse chromosome 5 was derived from a series of overlapping bacterial artificial chromosomes. A 296 kb segment from the physical map, spanning Ache to Tfr2, was compared with 267 kb of human sequence. We identified a conserved linkage of 12 genes including an open reading frame flanked by Ache and Asr2, a novel cation-chloride cotransporter interacting protein Cip1, Ephb4, Zan and Perq1. While some of these genes have been previously described, in each case we present new data derived from our comparative sequence analysis. Adjacent unfinished sequence data from the mouse contains an orthologous block of 10 additional genes including three novel cDNA sequences that we subsequently mapped to human 7q22. Methods for displaying comparative genomic information, including unfinished sequence data, are becoming increasingly important. We supplement our printed comparative analysis with a new, Web-based program called Laj (local alignments with java). Laj provides interactive access to archived pairwise sequence alignments via the WWW. It displays synchronized views of a dot-plot, a percent identity plot, a nucleotide-level local alignment and a variety of relevant annotations. Our mouse–human comparison can be viewed at http://web.uvic.ca/~bioweb/laj.html. Laj is available at http://bio.cse.psu.edu/, along with online documentation and additional examples of annotated genomic regions.

Genomic analysis of orthologous mouse and human olfactory receptor loci

Olfactory receptor (OR) genes represent ≈1% of genomic coding sequence in mammals, and these genes are clustered on multiple chromosomes in both the mouse and human genomes. We have taken a comparative genomics approach to identify features that may be involved in the dynamic evolution of this gene family and in the transcriptional control that results in a single OR gene expressed per olfactory neuron. We sequenced ≈350 kb of the murine P2 OR cluster and used synteny, gene linkage, and phylogenetic analysis to identify and sequence ≈111 kb of an orthologous cluster in the human genome. In total, 18 mouse and 8 human OR genes were identified, including 7 orthologs that appear to be functional in both species. Noncoding homology is evident between orthologs and generally is confined within the transcriptional unit. We find no evidence for common regulatory features shared among paralogs, and promoter regions generally do not contain strong promoter motifs. We discuss these observations, as well as OR clustering, in the context of evolutionary expansion and transcriptional regulation of OR repertoires.

Evolution and divergence of sodium channel genes in vertebrates

Invertebrate species possess one or two Na+ channel genes, yet there are 10 in mammals. When did this explosive growth come about during vertebrate evolution? All mammalian Na+ channel genes reside on four chromosomes. It has been suggested that this came about by multiple duplications of an ancestral chromosome with a single Na+ channel gene followed by tandem duplications of Na+ channel genes on some of these chromosomes. Because a large-scale expansion of the vertebrate genome likely occurred before the divergence of teleosts and tetrapods, we tested this hypothesis by cloning Na+ channel genes in a teleost fish. Using an approach designed to clone all of the Na+ channel genes in a genome, we found six Na+ channel genes. Phylogenetic comparisons show that each teleost gene is orthologous to a Na+ channel gene or gene cluster on a different mammalian chromosome, supporting the hypothesis that four Na+ channel genes were present in the ancestors of teleosts and tetrapods. Further duplications occurred independently in the teleost and tetrapod lineages, with a greater number of duplications in tetrapods. This pattern has implications for the evolution of function and specialization of Na+ channel genes in vertebrates. Sodium channel genes also are linked to homeobox (Hox) gene clusters in mammals. Using our phylogeny of Na+ channel genes to independently test between two models of Hox gene evolution, we support the hypothesis that Hox gene clusters evolved as (AB) (CD) rather than {D[A(BC)]}.

Cloning the human and mouse MMS19 genes and functional complementation of a yeast mms19 deletion mutant

The MMS19 gene of the yeast Saccharomyces cerevisiae encodes a polypeptide of unknown function which is required for both nucleotide excision repair (NER) and RNA polymerase II (RNAP II) transcription. Here we report the molecular cloning of human and mouse orthologs of the yeast MMS19 gene. Both human and Drosophila MMS19 cDNAs correct thermosensitive growth and sensitivity to killing by UV radiation in a yeast mutant deleted for the MMS19 gene, indicating functional conservation between the yeast and mammalian gene products. Alignment of the translated sequences of MMS19 from multiple eukaryotes, including mouse and human, revealed the presence of several conserved regions, including a HEAT repeat domain near the C-terminus. The presence of HEAT repeats, coupled with functional complementation of yeast mutant phenotypes by the orthologous protein from higher eukaryotes, suggests a role of Mms19 protein in the assembly of a multiprotein complex(es) required for NER and RNAP II transcription. Both the mouse and human genes are ubiquitously expressed as multiple transcripts, some of which appear to derive from alternative splicing. The ratio of different transcripts varies in several different tissue types.

Contrasting histories of avian and mammalian Mhc genes revealed by class II B sequences from songbirds.

To explore the evolutionary dynamics of genes in the major histocompatibility complex (Mhc) in nonmammalian vertebrates, we have amplified complete sequences of the polymorphic second (beta1) and third (beta2) exons of class II beta chain genes of songbirds. The pattern of nucleotide substitution in the antigen-binding site of sequences cloned from three behaviorally and phylogenetically divergent songbirds [scrub jays Aphelocoma coerulescens), red-winged blackbirds (Agelaius phoeniceus), and house finches (Carpodacus mexicanus) reveals that class II B genes of songbirds are subject to the same types of diversifying forces as those observed at mammalian class II loci. By contrast, the tree of avian class II B genes reveals that orthologous relationships have not been retained as in placental mammals and that, unlike class II genes in mammals, genes in songbirds and chickens have had very recent common ancestors within their respective groups. Thus, whereas the selective forces diversifying class II B genes of birds are likely similar to those in mammals, their long-term evolutionary dynamics appear to be characterized by much higher rates of concerted evolution.