Clustering of housekeeping genes provides a unified model of gene order in the human genome

It is often supposed that, except for tandem duplicates, genes are randomly distributed throughout the human genome. However, recent analyses suggest that when all the genes expressed in a given tissue (notably placenta and skeletal muscle) are examined, these genes do not map to random locations but instead resolve to clusters. We have asked three questions: (i) is this clustering true for most tissues, or are these the exceptions; (ii) is any clustering simply the result of the expression of tandem duplicates and (iii) how, if at all, does this relate to the observed clustering of genes with high expression rates? We provide a unified model of gene clustering that explains the previous observations. We examined Serial Analysis of Gene Expression (SAGE) data for 14 tissues and found significant clustering, in each tissue, that persists even after the removal of tandem duplicates. We confirmed clustering by analysis of independent expressed-sequence tag (EST) data. We then tested the possibility that the human genome is organized into subregions, each specializing in genes needed in a given tissue. By comparing genes expressed in different tissues, we show that this is not the case: those genes that seem to be tissue-specific in their expression do not, as a rule, cluster. We report that genes that are expressed in most tissues (housekeeping genes) show strong clustering. In addition, we show that the apparent clustering of genes with high expression rates is a consequence of the clustering of housekeeping genes.     

Genome rearrangement by replication-directed translocation

Gene order in bacteria is poorly conserved during evolution. For example, although many homologous genes are shared by the proteobacteria Escherichia coli, Haemophilus influenzae and Helicobacter pylori, their relative positions are very different in each genome, except local functional clusters such as operons. The complete sequences of the more closely related bacterial genomes, such as pairs of Chlamydia, H. pylori and Mycobacterium species, now allow identification of the processes and mechanisms involved in genome evolution. Here we provide evidence that a substantial proportion of rearrangements in gene order results from recombination sites that are determined by the positions of the replication forks. Our observations suggest that replication has a major role in directing genome evolution.     

Evidence for substantial fine-scale variation in recombination rates across the human genome

Characterizing fine-scale variation in human recombination rates is important, both to deepen understanding of the recombination process and to aid the design of disease association studies. Current genetic maps show that rates vary on a megabase scale, but studying finer-scale variation using pedigrees is difficult. Sperm-typing experiments have characterized regions where crossovers cluster into 1-2-kb hot spots, but technical difficulties limit the number of studies. An alternative is to use population variation to infer fine-scale characteristics of the recombination process. Several surveys reported 'block-like' patterns of diversity, which may reflect fine-scale recombination rate variation, but limitations of available methods made this impossible to assess. Here, we applied a new statistical method, which overcomes these limitations, to infer patterns of fine-scale recombination rate variation in 74 genes. We found extensive rate variation both within and among genes. In particular, recombination hot spots are a common feature of the human genome: 47% (35 of 74) of genes showed substantive evidence for a hot spot, and many more showed evidence for some rate variation. No primary sequence characteristics are consistently associated with precise hot-spot location, although G+C content and nucleotide diversity are correlated with local recombination rate.     

Disruption of an imprinted gene cluster by a targeted chromosomal translocation in mice

Genomic imprinting is an epigenetic process in which the activity of a gene is determined by its parent of origin. Mechanisms governing genomic imprinting are just beginning to be understood. However, the tendency of imprinted genes to exist in chromosomal clusters suggests a sharing of regulatory elements. To better understand imprinted gene clustering, we disrupted a cluster of imprinted genes on mouse distal chromosome 7 using the Cre/loxP recombination system. In mice carrying a site-specific translocation separating Cdkn1c and Kcnq1, imprinting of the genes retained on chromosome 7, including Kcnq1, Kcnq1ot1, Ascl2, H19 and Igf2, is unaffected, demonstrating that these genes are not regulated by elements near or telomeric to Cdkn1c. In contrast, expression and imprinting of the translocated Cdkn1c, Slc22a1l and Tssc3 on chromosome 11 are affected, consistent with the hypothesis that elements regulating both expression and imprinting of these genes lie within or proximal to Kcnq1. These data support the proposal that chromosomal abnormalities, including translocations, within KCNQ1 that are associated with the human disease Beckwith-Wiedemann syndrome (BWS) may disrupt CDKN1C expression. These results underscore the importance of gene clustering for the proper regulation of imprinted genes.     

Deletion of Peg10, an imprinted gene acquired from a retrotransposon, causes early embryonic lethality

By comparing mammalian genomes, we and others have identified actively transcribed Ty3/gypsy retrotransposon-derived genes with highly conserved DNA sequences and insertion sites. To elucidate the functions of evolutionarily conserved retrotransposon-derived genes in mammalian development, we produced mice that lack one of these genes, Peg10 (paternally expressed 10), which is a paternally expressed imprinted gene on mouse proximal chromosome 6. The Peg10 knockout mice showed early embryonic lethality owing to defects in the placenta. This indicates that Peg10 is critical for mouse parthenogenetic development and provides the first direct evidence of an essential role of an evolutionarily conserved retrotransposon-derived gene in mammalian development.     

Linkage disequilibrium mapping in isolated founder populations: diastrophic dysplasia in Finland.

Linkage disequilibrium mapping in isolated populations provides a powerful tool for fine structure localization of disease genes. Here, Luria and Delbrück's classical methods for analysing bacterial cultures are adapted to the study of human isolated founder populations in order to estimate (i) the recombination fraction between a disease locus and a marker; (ii) the expected degree of allelic homogeneity in a population; and (iii) the mutation rate of marker loci. Using these methods, we report striking linkage disequilibrium for diastrophic dysplasia (DTD) in Finland indicating that the DTD gene should lie within 0.06 centimorgans (or about 60 kilobases) of the CSF1R gene. Predictions about allelic homogeneity in Finland and mutation rates in simple sequence repeats are confirmed by independent observations.     

Intragenic tandem repeats generate functional variability

Tandemly repeated DNA sequences are highly dynamic components of genomes. Most repeats are in intergenic regions, but some are in coding sequences or pseudogenes. In humans, expansion of intragenic triplet repeats is associated with various diseases, including Huntington chorea and fragile X syndrome. The persistence of intragenic repeats in genomes suggests that there is a compensating benefit. Here we show that in the genome of Saccharomyces cerevisiae, most genes containing intragenic repeats encode cell-wall proteins. The repeats trigger frequent recombination events in the gene or between the gene and a pseudogene, causing expansion and contraction in the gene size. This size variation creates quantitative alterations in phenotypes (e.g., adhesion, flocculation or biofilm formation). We propose that variation in intragenic repeat number provides the functional diversity of cell surface antigens that, in fungi and other pathogens, allows rapid adaptation to the environment and elusion of the host immune system.     

Adaptation shapes patterns of genome evolution on sexual and asexual chromosomes in Drosophila

What advantage might sexual recombination confer? Population genetics theory predicts that asexual genomes are less efficient at eliminating deleterious mutations and incorporating beneficial alleles. Here, I compare patterns of genome evolution in a 40-kb gene-rich region on homologous neo-sex chromosomes of Drosophila miranda. Genes on the non-recombining neo-Y show various signs of degeneration, including transposable-element insertions, frameshift mutations and a higher rate of amino-acid substitution. In contrast, loci on the recombining neo-X show intact open reading frames and generally low rates of amino-acid substitution. One exceptional gene on the neo-X shows evidence for adaptive protein evolution, affecting patterns of variability at neighboring regions along the chromosome. These findings illustrate the limits to natural selection in an asexual genome. Deleterious mutations, including repetitive DNA, accumulate on a non-recombining chromosome, whereas rapid protein evolution due to positive selection is confined to the recombining homolog.     

A conserved microRNA module exerts homeotic control over Petunia hybrida and Antirrhinum majus floral organ identity

It is commonly thought that deep phylogenetic conservation of plant microRNAs (miRNAs) and their targets indicates conserved regulatory functions. We show that the blind (bl) mutant of Petunia hybrida and the fistulata (fis) mutant of Antirrhinum majus, which have similar homeotic phenotypes, are recessive alleles of two homologous miRNA-encoding genes. The BL and FIS genes control the spatial restriction of homeotic class C genes to the inner floral whorls, but their ubiquitous early floral expression patterns are in contradiction with a potential role in patterning C gene expression. We provide genetic evidence for the unexpected function of the MIRFIS and MIRBL genes in the center of the flower and propose a dynamic mechanism underlying their regulatory role. Notably, Arabidopsis thaliana, a more distantly related species, also contains this miRNA module but does not seem to use it to confine early C gene expression to the center of the flower.     

Fine-scale recombination patterns differ between chimpanzees and humans

Recombination rates seem to vary extensively along the human genome. Pedigree analysis suggests that rates vary by an order of magnitude when measured at the megabase scale, and at a finer scale, sperm typing studies point to the existence of recombination hotspots. These are short regions (1-2 kb) in which recombination rates are 10-1,000 times higher than the background rate. Less is known about how recombination rates change over time. Here we determined to what degree recombination rates are conserved among closely related species by estimating recombination rates from 14 Mb of linkage disequilibrium data in central chimpanzee and human populations. The results suggest that recombination hotspots are not conserved between the two species and that recombination rates in larger (50 kb) genomic regions are only weakly conserved. Therefore, the recombination landscape has changed markedly between the two species.     

A high-resolution survey of deletion polymorphism in the human genome

Recent work has shown that copy number polymorphism is an important class of genetic variation in human genomes. Here we report a new method that uses SNP genotype data from parent-offspring trios to identify polymorphic deletions. We applied this method to data from the International HapMap Project to produce the first high-resolution population surveys of deletion polymorphism. Approximately 100 of these deletions have been experimentally validated using comparative genome hybridization on tiling-resolution oligonucleotide microarrays. Our analysis identifies a total of 586 distinct regions that harbor deletion polymorphisms in one or more of the families. Notably, we estimate that typical individuals are hemizygous for roughly 30-50 deletions larger than 5 kb, totaling around 550-750 kb of euchromatic sequence across their genomes. The detected deletions span a total of 267 known and predicted genes. Overall, however, the deleted regions are relatively gene-poor, consistent with the action of purifying selection against deletions. Deletion polymorphisms may well have an important role in the genetics of complex traits; however, they are not directly observed in most current gene mapping studies. Our new method will permit the identification of deletion polymorphisms in high-density SNP surveys of trio or other family data.     

Recombination of mitochondrial DNA in skeletal muscle of individuals with multiple mitochondrial DNA heteroplasmy

Experimental evidence for human mitochondrial DNA (mtDNA) recombination was recently obtained in an individual with paternal inheritance of mtDNA and in an in vitro cell culture system. Whether mtDNA recombination is a common event in humans remained to be determined. To detect mtDNA recombination in human skeletal muscle, we analyzed the distribution of alleles in individuals with multiple mtDNA heteroplasmy using single-cell PCR and allele-specific PCR. In all ten individuals who carried a heteroplasmic D-loop mutation and a distantly located tRNA point mutation or a large deletion, we observed a mixture of four allelic combinations (tetraplasmy), a hallmark of recombination. Twelve of 14 individuals with closely located heteroplasmic D-loop mutation pairs contained a mixture of only three types of mitochondrial genomes (triplasmy), consistent with the absence of recombination between adjacent markers. These findings indicate that mtDNA recombination is common in human skeletal muscle.     

Birth of a metabolic gene cluster in yeast by adaptive gene relocation

Although most eukaryotic genomes lack operons, they contain some physical clusters of genes that are related in function despite being unrelated in sequence. How these clusters are formed during evolution is unknown. The DAL cluster is the largest metabolic gene cluster in yeast and consists of six adjacent genes encoding proteins that enable Saccharomyces cerevisiae to use allantoin as a nitrogen source. We show here that the DAL cluster was assembled, quite recently in evolutionary terms, through a set of genomic rearrangements that happened almost simultaneously. Six of the eight genes involved in allantoin degradation, which were previously scattered around the genome, became relocated to a single subtelomeric site in an ancestor of S. cerevisiae and Saccharomyces castellii. These genomic rearrangements coincided with a biochemical reorganization of the purine degradation pathway, which switched to importing allantoin instead of urate. This change eliminated urate oxidase, one of several oxygen-consuming enzymes that were lost by yeasts that can grow vigorously in anaerobic conditions. The DAL cluster is located in a domain of modified chromatin involving both H2A.Z histone exchange and Hst1-Sum1-mediated histone deacetylation, and it may be a coadapted gene complex formed by epistatic selection.     

Epistatic buffering of fitness loss in yeast double deletion strains

Interactions between deleterious mutations have been insufficiently studied, despite the fact that their strength and direction are critical for understanding the evolution of genetic recombination and the buildup of mutational load in populations. We compiled a list of 758 yeast gene deletions causing growth defects (from the Munich Information Center for Protein Sequences database and ref. 7). Using BY4741 and BY4742 single-deletion strains, we carried out 639 random crosses and assayed growth curves of the resulting progeny. We show that the maximum growth rate averaged over strains lacking deletions and those with double deletions is higher than that of strains with single deletions, indicating a positive epistatic effect. This tendency is shared by genes belonging to a variety of functional classes. Based on our data and former theoretical work, we suggest that epistasis is likely to diminish the negative effects of mutations when the ability to produce biomass at high rates contributes significantly to fitness.