Evidence of en bloc duplication in vertebrate genomes

It has been 30 years since it was first proposed that the vertebrate genome evolved through several rounds of genome-wide duplications (polyploidizations). Despite rapid advances in genetics, including sequencing of the complete genomes of several divergent species, this hypothesis has not been tested rigorously and is still a matter of debate. If polyploidizations occurred during chordate evolution, there should be a network of paralogous regions in the present-day jawed vertebrate (Gnathostomata) genomes. Here we present an investigation of the major histocompatibility complex (MHC) paralogous regions, which we accomplished by characterizing the corresponding region in amphioxus by identifying nine anchor genes and sequencing both the anchor genes and the regions that flank them (a total of 400 kb). Phylogenetic analysis of 31 genes (including the anchor genes) in these regions shows that duplications occurred after the divergence of cephalochordates and vertebrates but before the Gnathostomata radiation. The distribution of human and amphioxus orthologs in their respective genomes and the relationship between these distributions support the en bloc duplication events. Our analysis represents the first step towards demonstrating that the human ancestral genome has undergone polyploidization. Moreover, reconstruction of the pre-duplicated region indicates that one of the duplicated regions retains the ancestral organization.     

Discovery of previously unidentified genomic disorders from the duplication architecture of the human genome

Genomic disorders are characterized by the presence of flanking segmental duplications that predispose these regions to recurrent rearrangement. Based on the duplication architecture of the genome, we investigated 130 regions that we hypothesized as candidates for previously undescribed genomic disorders. We tested 290 individuals with mental retardation by BAC array comparative genomic hybridization and identified 16 pathogenic rearrangements, including de novo microdeletions of 17q21.31 found in four individuals. Using oligonucleotide arrays, we refined the breakpoints of this microdeletion, defining a 478-kb critical region containing six genes that were deleted in all four individuals. We mapped the breakpoints of this deletion and of four other pathogenic rearrangements in 1q21.1, 15q13, 15q24 and 17q12 to flanking segmental duplications, suggesting that these are also sites of recurrent rearrangement. In common with the 17q21.31 deletion, these breakpoint regions are sites of copy number polymorphism in controls, indicating that these may be inherently unstable genomic regions.     

Conservation of hotspots for recombination in low-copy repeats associated with the NF1 microdeletion

Several large-scale studies of human genetic variation have provided insights into processes such as recombination that have shaped human diversity. However, regions such as low-copy repeats (LCRs) have proven difficult to characterize, hindering efforts to understand the processes operating in these regions. We present a detailed study of genetic variation and underlying recombination processes in two copies of an LCR (NF1REPa and NF1REPc) on chromosome 17 involved in the generation of NF1 microdeletions and in a third copy (REP19) on chromosome 19 from which the others originated over 6.7 million years ago. We find evidence for shared hotspots of recombination among the LCRs. REP19 seems to contain hotspots in the same place as the nonallelic recombination hotspots in NF1REPa and NF1REPc. This apparent conservation of patterns of recombination hotspots in moderately diverged paralogous regions contrasts with recent evidence that these patterns are not conserved in less-diverged orthologous regions of chimpanzees.     

Conserved noncoding sequences are selectively constrained and not mutation cold spots

Noncoding genetic variants are likely to influence human biology and disease, but recognizing functional noncoding variants is difficult. Approximately 3% of noncoding sequence is conserved among distantly related mammals, suggesting that these evolutionarily conserved noncoding regions (CNCs) are selectively constrained and contain functional variation. However, CNCs could also merely represent regions with lower local mutation rates. Here we address this issue and show that CNCs are selectively constrained in humans by analyzing HapMap genotype data. Specifically, new (derived) alleles of SNPs within CNCs are rarer than new alleles in nonconserved regions (P = 3 x 10(-18)), indicating that evolutionary pressure has suppressed CNC-derived allele frequencies. Intronic CNCs and CNCs near genes show greater allele frequency shifts, with magnitudes comparable to those for missense variants. Thus, conserved noncoding variants are more likely to be functional. Allele frequency distributions highlight selectively constrained genomic regions that should be intensively surveyed for functionally important variation.     

Extensive genomic duplication during early chordate evolution

Opinions on the hypothesis that ancient genome duplications contributed to the vertebrate genome range from strong skepticism to strong credence. Previous studies concentrated on small numbers of gene families or chromosomal regions that might not have been representative of the whole genome, or used subjective methods to identify paralogous genes and regions. Here we report a systematic and objective analysis of the draft human genome sequence to identify paralogous chromosomal regions (paralogons) formed during chordate evolution and to estimate the ages of duplicate genes. We found that the human genome contains many more paralogons than would be expected by chance. Molecular clock analysis of all protein families in humans that have orthologs in the fly and nematode indicated that a burst of gene duplication activity took place in the period 350 650 Myr ago and that many of the duplicate genes formed at this time are located within paralogons. Our results support the contention that many of the gene families in vertebrates were formed or expanded by large-scale DNA duplications in an early chordate. Considering the incompleteness of the sequence data and the antiquity of the event, the results are compatible with at least one round of polyploidy.     

Juxtaposed regions of extensive and minimal linkage disequilibrium in human Xq25 and Xq28

Linkage disequilibrium (LD), or the non-random association of alleles, is poorly understood in the human genome. Population genetic theory suggests that LD is determined by the age of the markers, population history, recombination rate, selection and genetic drift. Despite the uncertainties in determining the relative contributions of these factors, some groups have argued that LD is a simple function of distance between markers. Disease-gene mapping studies and a simulation study gave differing predictions on the degree of LD in isolated and general populations. In view of the discrepancies between theory and experimental observations, we constructed a high-density SNP map of the Xq25-Xq28 region and analysed the male genotypes and haplotypes across this region for LD in three populations. The populations included an outbred European sample (CEPH males) and isolated population samples from Finland and Sardinia. We found two extended regions of strong LD bracketed by regions with no evidence for LD in all three samples. Haplotype analysis showed a paucity of haplotypes in regions of strong LD. Our results suggest that, in this region of the X chromosome, LD is not a monotonic function of the distance between markers, but is more a property of the particular location in the human genome.     

Dichotomy of single-nucleotide polymorphism haplotypes in olfactory receptor genes and pseudogenes

Substantial efforts are focused on identifying single-nucleotide polymorphisms (SNPs) throughout the human genome, particularly in coding regions (cSNPs), for both linkage disequilibrium and association studies. Less attention, however, has been directed to the clarification of evolutionary processes that are responsible for the variability in nucleotide diversity among different regions of the genome. We report here the population sequence diversity of genomic segments within a 450-kb cluster of olfactory receptor (OR) genes on human chromosome 17. We found a dichotomy in the pattern of nucleotide diversity between OR pseudogenes and introns on the one hand and the closely interspersed intact genes on the other. We suggest that weak positive selection is responsible for the observed patterns of genetic variation. This is inferred from a lower ratio of polymorphism to divergence in genes compared with pseudogenes or introns, high non-synonymous substitution rates in OR genes, and a small but significant overall reduction in variability in the entire OR gene cluster compared with other genomic regions. The dichotomy among functionally different segments within a short genomic distance requires high recombination rates within this OR cluster. Our work demonstrates the impact of weak positive selection on human nucleotide diversity, and has implications for the evolution of the olfactory repertoire.     

Microsatellites are preferentially associated with nonrepetitive DNA in plant genomes.

Microsatellites are a ubiquitous class of simple repetitive DNA sequence. An excess of such repetitive tracts has been described in all eukaryotes analyzed and is thought to result from the mutational effects of replication slippage. Large-scale genomic and EST sequencing provides the opportunity to evaluate the abundance and relative distribution of microsatellites between transcribed and nontranscribed regions and the relationship of these features to haploid genome size. Although this has been studied in microbial and animal genomes, information in plants is limited. We assessed microsatellite frequency in plant species with a 50-fold range in genome size that is mostly attributable to the recent amplification of repetitive DNA. Among species, the overall frequency of microsatellites was inversely related to genome size and to the proportion of repetitive DNA but remained constant in the transcribed portion of the genome. This indicates that most microsatellites reside in regions pre-dating the recent genome expansion in many plants. The microsatellite frequency was higher in transcribed regions, especially in the untranslated portions, than in genomic DNA. Contrary to previous reports suggesting a preferential mechanism for the origin of microsatellites from repetitive DNA in both animals and plants, our findings show a significant association with the low-copy fraction of plant genomes.     

A complete BAC-based physical map of the Arabidopsis thaliana genome.

Arabidopsis thaliana is a small flowering plant that serves as the major model system in plant molecular genetics. The efforts of many scientists have produced genetic maps that provide extensive coverage of the genome (http://genome-www. stanford.edu/Arabidopsis/maps.html). Recently, detailed YAC, BAC, P1 and cosmid-based physical maps (that is, representations of genomic regions as sets of overlapping clones of corresponding libraries) have been established that extend over wide genomic areas ranging from several hundreds of kilobases to entire chromosomes. These maps provide an entry to gain deeper insight into the A. thaliana genome structure. A. thaliana has been chosen as the subject of the first large-scale project intended to determine the full genome sequence of a plant. This sequencing project, together with the increasing interest in map-based gene cloning, has highlighted the requirement for a complete and accurate physical map of this plant species. To supply the scientific community with a high-quality resource, we present here a complete physical map of A. thaliana using essentially the IGF BAC library. The map consists of 27 contigs that cover the entire genome, except for the presumptive centromeric regions, nucleolar organization regions (NOR) and telomeric areas. This is the first reported map of a complex organism based entirely on BAC clones and it represents the most homogeneous and complete physical map established to date for any plant genome. Furthermore, the analysis performed here serves as a model for an efficient physical mapping procedure using BAC clones that can be applied to other complex genomes.     

Mosaicism for a specific somatic mitochondrial DNA mutation in adult human brain.

The levels of a specific mitochondrial DNA deletion (mtDNA4977) measured in 12 brain regions of 6 normal adults 39 to 82 years old exhibited striking variation among anatomical locations. Comparisons of the same region among individuals showed an increase of mtDNA4977 with age. The three regions with the highest levels, caudate, putamen and substantia nigra, are characterized by a high dopamine metabolism. The breakdown of dopamine by mitochondrial MAO produces H2O2 which can lead to oxygen radical formation. We suggest that mtDNA4977 may be the "tip of the iceberg" of the spectrum of somatic mutations produced by oxidative damage.     

On the subspecific origin of the laboratory mouse

The genome of the laboratory mouse is thought to be a mosaic of regions with distinct subspecific origins. We have developed a high-resolution map of the origin of the laboratory mouse by generating 25,400 phylogenetic trees at 100-kb intervals spanning the genome. On average, 92% of the genome is of Mus musculus domesticus origin, and the distribution of diversity is markedly nonrandom among the chromosomes. There are large regions of extremely low diversity, which represent blind spots for studies of natural variation and complex traits, and hot spots of diversity. In contrast with the mosaic model, we found that most of the genome has intermediate levels of variation of intrasubspecific origin. Finally, mouse strains derived from the wild that are supposed to represent different mouse subspecies show substantial intersubspecific introgression, which has strong implications for evolutionary studies that assume these are pure representatives of a given subspecies.     

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.     

Genome-wide association study identifies new susceptibility loci for Crohn disease and implicates autophagy in disease pathogenesis

We present a genome-wide association study of ileal Crohn disease and two independent replication studies that identify several new regions of association to Crohn disease. Specifically, in addition to the previously established CARD15 and IL23R associations, we identified strong and significantly replicated associations (combined P < 10(-10)) with an intergenic region on 10q21.1 and a coding variant in ATG16L1, the latter of which was also recently reported by another group. We also report strong associations with independent replication to variation in the genomic regions encoding PHOX2B, NCF4 and a predicted gene on 16q24.1 (FAM92B). Finally, we demonstrate that ATG16L1 is expressed in intestinal epithelial cell lines and that functional knockdown of this gene abrogates autophagy of Salmonella typhimurium. Together, these findings suggest that autophagy and host cell responses to intracellular microbes are involved in the pathogenesis of Crohn disease.     

Exclusion of DNA changes in the beta-subunit of the c-GMP phosphodiesterase gene as the cause for Huntington's disease.

To identify expressed sequences within candidate regions for the Huntington's disease (HD) gene in 4p16.3, we isolated the gene encoding the beta subunit of the human cGMP phosphodiesterase (PDEB). We formally assessed this as a candidate gene for HD based on it's expression in brain, the demonstration of linkage disequilibrium between intragenic DNA markers and HD, and the demonstration that mice with a mutation in this gene have a reduction of neurons in particular brain regions. We investigated all 22 exons of PDEB and 5'-flanking region for point mutations in 16 HD patients of different ethnic origins using single strand conformational polymorphism analysis. The underlying DNA changes found initially exclusively in HD patients were excluded as the cause for HD.