Large multi-chromosomal duplications encompass many members of the olfactory receptor gene family in the human genome

Olfactory receptor genes and chromosome 11 structural aberrations: Players or spectators?

The largest multi-gene family in metazoans is the family of olfactory receptor (OR) genes. Human ORs are organized in clusters over most chromosomes and seem to include >0.1% the human genome. Because 369 out of 856 OR genes are mapped on chromosome 11 (HSA11), we sought to determine whether they mediate structural rearrangements involving this chromosome. To this aim, we analyzed 220 specimens collected during diagnostic procedures involving structural rearrangements of chromosome 11. A total of 222 chromosomal abnormalities were included, consisting of inversions, deletions, translocations, duplications, and one insertion, detected by conventional chromosome analysis and/or fluorescence in situ hybridization (FISH) and array comparative genomic hybridization (array-CGH). We verified by bioinformatics and statistical approaches the occurrence of breakpoints in cytobands with or without OR genes. We found that OR genes are not involved in chromosome 11 reciprocal translocations, suggesting that different DNA motifs and mechanisms based on homology or non-homology recombination can cause chromosome 11 structural alterations. We also considered the proximity between the chromosomal territories of chromosome 11 and its partner chromosomes involved in the translocations by using the deposited Hi-C data concerning the possible occurrence of chromosome interactions. Interestingly, most of the breakpoints are located in regions highly involved in chromosome interactions. Further studies should be carried out to confirm the potential role of chromosome territories' proximity in promoting genome structural variation, so fundamental in our understanding of the molecular basis of medical genetics and evolutionary genetics.

Ethnic and functional differentiation of copy number polymorphisms in Tunisian and HapMap population unveils insights on genome organizational plasticity

Admixture mapping has been useful in identifying genetic variations linked to phenotypes, adaptation and diseases. Copy number variations (CNVs) represents genomic structural variants spanning large regions of chromosomes reaching several megabases. In this investigation, the "Canary" algorithm was applied to 102 Tunisian samples and 991 individuals from eleven HapMap III populations to genotype 1279 copy number polymorphisms (CNPs). In this present work, we investigate the Tunisian population structure using the CNP makers previously identified among Tunisian. The study revealed that Sub-Saharan African populations exhibited the highest diversity with the highest proportions of allelic CNPs. Among all the African populations, Tunisia showed the least diversity. Individual ancestry proportions computed using STRUCTURE analysis revealed a major European component among Tunisians with lesser contribution from Sub-Saharan Africa and Asia. Population structure analysis indicated the genetic proximity with Europeans and noticeable distance from the Sub-Saharan African and East Asian clusters. Seven genes harbouring Tunisian high-frequent CNPs were identified known to be associated with 9 Mendelian diseases and/or phenotypes. Functional annotation of genes under selection highlighted a noteworthy enrichment of biological processes to receptor pathway and activity as well as glutathione metabolism. Additionally, pathways of potential concern for health such as drug metabolism, infectious diseases and cancers exhibited significant enrichment. The distinctive genetic makeup of the Tunisians might have been influenced by various factors including natural selection and genetic drift, resulting in the development of distinct genetic variations playing roles in specific biological processes. Our research provides a justification for focusing on the exclusive genome organization of this population and uncovers previously overlooked elements of the genome.

A map of copy number variations in the Tunisian population: a valuable tool for medical genomics in North Africa

Copy number variation (CNV) is considered as the most frequent type of structural variation in the human genome. Some CNVs can act on human phenotype diversity, encompassing rare Mendelian diseases and genomic disorders. The North African populations remain underrepresented in public genetic databases in terms of single-nucleotide variants as well as for larger genomic mutations. In this study, we present the first CNV map for a North African population using the Affymetrix Genome-Wide SNP (single-nucleotide polymorphism) array 6.0 array genotyping intensity data to call CNVs in 102 Tunisian healthy individuals. Two softwares, PennCNV and Birdsuite, were used to call CNVs in order to provide reliable data. Subsequent bioinformatic analyses were performed to explore their features and patterns. The CNV map of the Tunisian population includes 1083 CNVs spanning 61.443 Mb of the genome. The CNV length ranged from 1.017 kb to 2.074 Mb with an average of 56.734 kb. Deletions represent 57.43% of the identified CNVs, while duplications and the mixed loci are less represented. One hundred and three genes disrupted by CNVs are reported to cause 155 Mendelian diseases/phenotypes. Drug response genes were also reported to be affected by CNVs. Data on genes overlapped by deletions and duplications segments and the sequence properties in and around them also provided insights into the functional and health impacts of CNVs. These findings represent valuable clues to genetic diversity and personalized medicine in the Tunisian population as well as in the ethnically similar populations from North Africa.

High-throughput single-molecule mapping links subtelomeric variants and long-range haplotypes with specific telomeres

Accurate maps and DNA sequences for human subtelomere regions, along with detailed knowledge of subtelomere variation and long-range telomere-terminal haplotypes in individuals, are critical for understanding telomere function and its roles in human biology. Here, we use a highly automated whole genome mapping technology in nano-channel arrays to analyze large terminal human chromosome segments extending from chromosome-specific subtelomere sequences through subtelomeric repeat regions to terminal (TTAGGG)n repeat tracts. We establish detailed maps for subtelomere gap regions in the human reference sequence, detect many new large subtelomeric variants and demonstrate the feasibility of long-range haplotyping through segmentally duplicated subtelomere regions. These features make the method a uniquely valuable new tool for improving the quality of genome assemblies in complex DNA regions. Based on single molecule mapping of telomere-terminal DNA fragments, we provide proof of principle for a novel method to estimate telomere lengths linked to distinguishable telomeric haplotypes; this single-telomere genotyping method may ultimately enable delineation of human cis elements involved in telomere length regulation.

Clinical detection of deletion structural variants in whole-genome sequences

Optimal management of acutely ill infants with monogenetic diseases requires rapid identification of causative haplotypes. Whole-genome sequencing (WGS) has been shown to identify pathogenic nucleotide variants in such infants. Deletion structural variants (DSVs, >50 nt) are implicated in many genetic diseases, and tools have been designed to identify DSVs using short-read WGS. Optimisation and integration of these tools into a WGS pipeline could improve diagnostic sensitivity and specificity of WGS. In addition, it may improve turnaround time when compared with current CNV assays, enhancing utility in acute settings. Here we describe DSV detection methods for use in WGS for rapid diagnosis in acutely ill infants: SKALD (Screening Konsensus and Annotation of Large Deletions) combines calls from two tools (Breakdancer and GenomeStrip) with calibrated filters and clinical interpretation rules. In four WGS runs, the average analytic precision (positive predictive value) of SKALD was 78%, and recall (sensitivity) was 27%, when compared with validated reference DSV calls. When retrospectively applied to a cohort of 36 families with acutely ill infants SKALD identified causative DSVs in two. The first was heterozygous deletion of exons 1–3 of MMP21 in trans with a heterozygous frame-shift deletion in two siblings with transposition of the great arteries and heterotaxy. In a newborn female with dysmorphic features, ventricular septal defect and persistent pulmonary hypertension, SKALD identified the breakpoints of a heterozygous, de novo 1p36.32p36.13 deletion. In summary, consensus DSV calling, implemented in an 8-h computational pipeline with parameterised filtering, has the potential to increase the diagnostic yield of WGS in acutely ill neonates and discover novel disease genes.

Chromosomal Rearrangements as Barriers to Genetic Homogenization between Archaic and Modern Humans

Chromosomal rearrangements, which shuffle DNA throughout the genome, are an important source of divergence across taxa. Using a paired-end read approach with Illumina sequence data for archaic humans, I identify changes in genome structure that occurred recently in human evolution. Hundreds of rearrangements indicate genomic trafficking between the sex chromosomes and autosomes, raising the possibility of sex-specific changes. Additionally, genes adjacent to genome structure changes in Neanderthals are associated with testis-specific expression, consistent with evolutionary theory that new genes commonly form with expression in the testes. I identify one case of new-gene creation through transposition from the Y chromosome to chromosome 10 that combines the 5'-end of the testis-specific gene Fank1 with previously untranscribed sequence. This new transcript experienced copy number expansion in archaic genomes, indicating rapid genomic change. Among rearrangements identified in Neanderthals, 13% are transposition of selfish genetic elements, whereas 32% appear to be ectopic exchange between repeats. In Denisovan, the pattern is similar but numbers are significantly higher with 18% of rearrangements reflecting transposition and 40% ectopic exchange between distantly related repeats. There is an excess of divergent rearrangements relative to polymorphism in Denisovan, which might result from nonuniform rates of mutation, possibly reflecting a burst of transposable element activity in the lineage that led to Denisovan. Finally, loci containing genome structure changes show diminished rates of introgression from Neanderthals into modern humans, consistent with the hypothesis that rearrangements serve as barriers to gene flow during hybridization. Together, these results suggest that this previously unidentified source of genomic variation has important biological consequences in human evolution.

Improved genome inference in the MHC using a population reference graph

Although much is known about human genetic variation, such information is typically ignored in assembling new genomes. Instead, reads are mapped to a single reference, which can lead to poor characterization of regions of high sequence or structural diversity. We introduce a population reference graph, which combines multiple reference sequences and catalogs of variation. The genomes of new samples are reconstructed as paths through the graph using an efficient hidden Markov model, allowing for recombination between different haplotypes and additional variants. By applying the method to the 4.5-Mb extended MHC region on human chromosome 6, combining 8 assembled haplotypes, the sequences of known classical HLA alleles and 87,640 SNP variants from the 1000 Genomes Project, we demonstrate using simulations, SNP genotyping, and short-read and long-read data how the method improves the accuracy of genome inference and identified regions where the current set of reference sequences is substantially incomplete.

Mammalian olfactory receptors: molecular mechanisms of odorant detection, 3D-modeling, and structure-activity relationships

This chapter describes the main characteristics of olfactory receptor (OR) genes of vertebrates, including generation of this large multigenic family and pseudogenization. OR genes are compared in relation to evolution and among species. OR gene structure and selection of a given gene for expression in an olfactory sensory neuron (OSN) are tackled. The specificities of OR proteins, their expression, and their function are presented. The expression of OR proteins in locations other than the nasal cavity is regulated by different mechanisms, and ORs display various additional functions. A conventional olfactory signal transduction cascade is observed in OSNs, but individual ORs can also mediate different signaling pathways, through the involvement of other molecular partners and depending on the odorant ligand encountered. ORs are engaged in constitutive dimers. Ligand binding induces conformational changes in the ORs that regulate their level of activity depending on odorant dose. When present, odorant binding proteins induce an allosteric modulation of OR activity. Since no 3D structure of an OR has been yet resolved, modeling has to be performed using the closest G-protein-coupled receptor 3D structures available, to facilitate virtual ligand screening using the models. The study of odorant binding modes and affinities may infer best-bet OR ligands, to be subsequently checked experimentally. The relationship between spatial and steric features of odorants and their activity in terms of perceived odor quality are also fields of research that development of computing tools may enhance.

Genome structural variation discovery and genotyping

Comparisons of human genomes show that more base pairs are altered as a result of structural variation - including copy number variation - than as a result of point mutations. Here we review advances and challenges in the discovery and genotyping of structural variation. The recent application of massively parallel sequencing methods has complemented microarray-based methods and has led to an exponential increase in the discovery of smaller structural-variation events. Some global discovery biases remain, but the integration of experimental and computational approaches is proving fruitful for accurate characterization of the copy, content and structure of variable regions. We argue that the long-term goal should be routine, cost-effective and high quality de novo assembly of human genomes to comprehensively assess all classes of structural variation.

Evolutionary adaptations to dietary changes

Through cultural innovation and changes in habitat and ecology, there have been a number of major dietary shifts in human evolution, including meat eating, cooking, and those associated with plant and animal domestication. The identification of signatures of adaptations to such dietary changes in the genome of extant primates (including humans) may shed light not only on the evolutionary history of our species, but also on the mechanisms that underlie common metabolic diseases in modern human populations. In this review, we provide a brief overview of the major dietary shifts that occurred during hominin evolution, and we discuss the methods and approaches used to identify signals of natural selection in patterns of sequence variation. We then review the results of studies aimed at detecting the genetic loci that played a major role in dietary adaptations and conclude by outlining the potential of future studies in this area.

Global features of neural activity in the olfactory system form a parallel code that predicts olfactory behavior and perception

Odor identity is coded in spatiotemporal patterns of neural activity in the olfactory bulb. Here we asked whether meaningful olfactory information could also be read from the global olfactory neural population response. We applied standard statistical methods of dimensionality-reduction to neural activity from 12 previously published studies using seven different species. Four studies reported olfactory receptor activity, seven reported glomerulus activity, and one reported the activity of projection-neurons. We found two linear axes of neural population activity that accounted for more than half of the variance in neural response across species. The first axis was correlated with the total sum of odor-induced neural activity, and reflected the behavior of approach or withdrawal in animals, and odorant pleasantness in humans. The second and orthogonal axis reflected odorant toxicity across species. We conclude that in parallel with spatiotemporal pattern coding, the olfactory system can use simple global computations to read vital olfactory information from the neural population response.

Population analysis of large copy number variants and hotspots of human genetic disease

Copy number variants (CNVs) contribute to human genetic and phenotypic diversity. However, the distribution of larger CNVs in the general population remains largely unexplored. We identify large variants in approximately 2500 individuals by using Illumina SNP data, with an emphasis on "hotspots" prone to recurrent mutations. We find variants larger than 500 kb in 5%-10% of individuals and variants greater than 1 Mb in 1%-2%. In contrast to previous studies, we find limited evidence for stratification of CNVs in geographically distinct human populations. Importantly, our sample size permits a robust distinction between truly rare and polymorphic but low-frequency copy number variation. We find that a significant fraction of individual CNVs larger than 100 kb are rare and that both gene density and size are strongly anticorrelated with allele frequency. Thus, although large CNVs commonly exist in normal individuals, which suggests that size alone can not be used as a predictor of pathogenicity, such variation is generally deleterious. Considering these observations, we combine our data with published CNVs from more than 12,000 individuals contrasting control and neurological disease collections. This analysis identifies known disease loci and highlights additional CNVs (e.g., 3q29, 16p12, and 15q25.2) for further investigation. This study provides one of the first analyses of large, rare (0.1%-1%) CNVs in the general population, with insights relevant to future analyses of genetic disease.

A highly conserved segmental duplication in the subtelomeres of Plasmodium falciparum chromosomes varies in copy number

Background

Segmental duplications (SD) have been found in genomes of various organisms, often accumulated at the ends of chromosomes. It has been assumed that the sequence homology in-between the SDs allow for ectopic interactions that may contribute to the emergence of new genes or gene variants through recombinatorial events.

Methods

In silico analysis of the 3D7 Plasmodium falciparum genome, conducted to investigate the subtelomeric compartments, led to the identification of subtelomeric SDs. Sequence variation and copy number polymorphisms of the SDs were studied by DNA sequencing, real-time quantitative PCR (qPCR) and fluorescent in situ hybridization (FISH). The levels of transcription and the developmental expression of copy number variant genes were investigated by qPCR.

Results

A block of six genes of >10 kilobases in size, including var, rif, pfmc-2tm and three hypothetical genes (n-, o- and q-gene), was found duplicated in the subtelomeric regions of chromosomes 1, 2, 3, 6, 7, 10 and 11 (SD1). The number of SD1 per genome was found to vary from 4 to 8 copies in between different parasites. The intragenic regions of SD1 were found to be highly conserved across ten distinct fresh and long-term cultivated P. falciparum. Sequence variation was detected in a ≈ 23 amino-acid long hypervariable region of a surface-exposed loop of PFMC-2TM. A hypothetical gene within SD1, the n-gene, encoding a PEXEL/VTS-containing two-transmembrane protein was found expressed in ring stage parasites. The n-gene transcription levels were found to correlate to the number of n-gene copies. Fragments of SD1 harbouring two or three of the SD1-genes (o-gene, pfmc-2tm, q-gene) were also found in the 3D7 genome. In addition a related second SD, SD2, of ≈ 55% sequence identity to SD1 was found duplicated in a fresh clinical isolate but was only present in a single copy in 3D7 and in other P. falciparum lines or clones.

Conclusion

Plasmodium falciparum carries multiple sequence conserved SDs in the otherwise highly variable subtelomeres of its chromosomes. The uniqueness of the SDs amongst plasmodium species, and the conserved nature of the genes within, is intriguing and suggests an important role of the SD to P. falciparum.

Co-regulation of a large and rapidly evolving repertoire of odorant receptor genes

The olfactory system meets niche- and species-specific demands by an accelerated evolution of its odorant receptor repertoires. In this review, we describe evolutionary processes that have shaped olfactory and vomeronasal receptor gene families in vertebrate genomes. We emphasize three important periods in the evolution of the olfactory system evident by comparative genomics: the adaptation to land in amphibian ancestors, the decline of olfaction in primates, and the delineation of putative pheromone receptors concurrent with rodent speciation. The rapid evolution of odorant receptor genes, the sheer size of the repertoire, as well as their wide distribution in the genome, presents a developmental challenge: how are these ever-changing odorant receptor repertoires coordinated within the olfactory system? A central organizing principle in olfaction is the specialization of sensory neurons resulting from each sensory neuron expressing only ~one odorant receptor allele. In this review, we also discuss this mutually exclusive expression of odorant receptor genes. We have considered several models to account for co-regulation of odorant receptor repertoires, as well as discussed a new hypothesis that invokes important epigenetic properties of the system.

The Facioscapulohumeral muscular dystrophy region on 4qter and the homologous locus on 10qter evolved independently under different evolutionary pressure

Background

The homologous 4q and 10q subtelomeric regions include two distinctive polymorphic arrays of 3.3 kb repeats, named D4Z4. An additional BlnI restriction site on the 10q-type sequence allows to distinguish the chromosomal origin of the repeats. Reduction in the number of D4Z4 repeats below a threshold of 10 at the 4q locus is tightly linked to Facioscapulohumeral Muscular Dystrophy (FSHD), while similar contractions at 10q locus, are not pathogenic. Sequence variations due to the presence of BlnI-sensitive repeats (10q-type) on chromosome 4 or viceversa of BlnI-resistant repeats (4q-type) on chromosome 10 are observed in both alleles.

Results

We analysed DNA samples from 116 healthy subiects and 114 FSHD patients and determined the size distributions of polymorphic 4q and 10q alleles, the frequency and the D4Z4 repeat assortment of variant alleles, and finally the telomeric sequences both in standard and variant alleles.

We observed the same frequency and types of variant alleles in FSHD patients and controls, but we found marked differences between the repeat arrays of the 4q and 10q chromosomes. In particular we detected 10q alleles completely replaced by the 4q subtelomeric region, consisting in the whole set of 4q-type repeats and the distal telomeric markers. However the reciprocal event, 10q-type subtelomeric region on chromosome 4, was never observed. At 4q locus we always identified hybrid alleles containing a mixture of 4q and 10q-type repeats.

Conclusion

The different size distribution and different structure of 10q variant alleles as compared with 4q suggests that these loci evolved in a different manner, since the 4q locus is linked to FSHD, while no inheritable disease is associated with mutations in 10qter genomic region. Hybrid alleles on chromosome 4 always retain a minimum number of 4q type repeats, as they are probably essential for maintaining the structural and functional properties of this subtelomeric region.

In addition we found: i) several instances of variant alleles that could be misinterpreted and interfere with a correct diagnosis of FSHD; ii) the presence of borderline alleles in the range of 30–40 kb that carried a qA type telomere and were not associated with the disease.

Analysis of segmental duplications reveals a distinct pattern of continuation-of-synteny between human and mouse genomes

About 5% of the human genome consists of large-scale duplicated segments of almost identical sequences. Segmental duplications (SDs) have been proposed to be involved in non-allelic homologous recombination leading to recurrent genomic variation and disease. It has also been suggested that these SDs are associated with syntenic rearrangements that have shaped the human genome. We have analyzed 14 members of a single family of closely related SDs in the human genome, some of which are associated with common inversion polymorphisms at chromosomes 8p23 and 4p16. Comparative analysis with the mouse genome revealed syntenic inversions for these two human polymorphic loci. In addition, 12 of the 14 SDs, while absent in the mouse genome, occur at the breaks of synteny; suggesting a non-random involvement of these sequences in genome evolution. Furthermore, we observed a syntenic familial relationship between 8 and 12 breakpoint-loci, where broken synteny that ends at one family member resumes at another, even across different chromosomes. Subsequent genome-wide assessment revealed that this relationship, which we named continuation-of-synteny, is not limited to the 8p23 family and occurs 46 times in the human genome with high frequency at specific chromosomes. Our analysis supports a non-random breakage model of genomic evolution with an active involvement of segmental duplications for specific regions of the human genome.

Novel deletion variants of 9q13–q21.12 and classical euchromatic variants of 9q12/qh involve deletion, duplication and triplication of large tracts of segmentally duplicated pericentromeric euchromatin

Large-scale copy number variation that is cytogenetically visible in normal individuals has been described as euchromatic variation but needs to be distinguished from pathogenic euchromatic deletion or duplication. Here, we report eight patients (three families and two individuals) with interstitial deletions of 9q13-q21.12. Fluorescence in situ hybridisation with a large panel of BACs showed that all the deleted clones were from extensive tracts of segmentally duplicated euchromatin, copies of which map to both the long and short arms of chromosome 9. The variety of reasons for which these patients were ascertained, and the phenotypically normal parents, indicates that this is a novel euchromatic variant with no phenotypic effect. Further, four patients with classical euchromatic variants of 9q12/qh or 9p12 were also shown to have duplications or triplications of this segmentally duplicated material common to both 9p and 9q. The cytogenetic boundaries between the segmentally duplicated regions and flanking unique sequences were mapped to 9p13.1 in the short arm (BAC RP11-402N8 at 38.7 Mb) and to 9q21.12 in the long arm (BAC RP11-88I18 at 70.3 Mb). The BACs identified in this study should in future make it possible to differentiate between clinically significant deletions or duplications and euchromatic variants with no established phenotypic consequences.

7E olfactory receptor gene clusters and evolutionary chromosome rearrangements

Olfactory receptor (OR) genes of the 7E subfamily have been duplicated to multiple regions throughout the human genome. Segmental duplications containing 7E OR genes have been associated with both pathological and evolutionary chromosome rearrangements. Many of these breakpoint regions coincide with breaks of chromosomal synteny in the mouse, rat and/or chicken genomes. Collectively, these data suggest that 7E OR-containing regions represent hot spots of genomic instability.

Duplications and copy number variants of 8p23.1 are cytogenetically indistinguishable but distinct at the molecular level

It has been proposed that duplications of 8p23.1 are either euchromatic variants of the 8p23.1 defensin domain with no phenotypic consequences or true duplications associated with developmental delay and heart defects. Here, we provide evidence for both alternatives in two new families. A duplication of most of band 8p23.1 (circa 5 Mb) was found in a girl of 8 years with pulmonary stenosis and mild language delay. BAC fluorescence in situ hybridisation (FISH) and multiplex amplifiable probe hybridisation (MAPH) showed that the two copies of the duplicated segment were sited, in an alternating fashion, between three copies of a circa 300-450 kb segment from 8p23.1 distal to REPD. Copy number of the variable 8p23.1 defensin domain was consistent with duplication but within the normal range. Duplication of the GATA-binding protein 4 gene (GATA4) in this patient and others with and without heart defects, suggests it is a dosage-sensitive gene with variable penetrance. A cytogenetically similar duplication of 8p23.1 was found at prenatal diagnosis in a fetus, father and grandmother. There was no duplication using BAC FISH but MAPH showed 11 copies of the 360 kb variable defensin domain which is within the expanded range found in previous euchromatic variant carriers. Semiquantitative FISH (SQ-FISH) was consistent with a simultaneous expansion of the adjacent olfactory receptor repeats. These results distinguish duplications of 8p23.1 with clinically significant consequences from benign copy number variants, which have not yet been associated with qualitative or quantitative traits.

Characterization of novel GPCR gene coding locus in amphioxus genome: gene structure, expression, and phylogenetic analysis with implications for its involvement in chemoreception

Chemosensation is the primary sensory modality in almost all metazoans. The vertebrate olfactory receptor genes exist as tandem clusters in the genome, so that identifying their evolutionary origin would be useful for understanding the expansion of the sensory world in relation to a large-scale genomic duplication event in a lineage leading to the vertebrates. In this study, I characterized a novel GPCR (G-protein-coupled receptor) gene-coding locus from the amphioxus genome. The genomic DNA contains an intronless ORF whose deduced amino acid sequence encodes a seven-transmembrane protein with some amino acid residues characteristic of vertebrate olfactory receptors (ORs). Surveying counterparts in the Ciona intestinalis (Asidiacea, Urochordata) genome by querying BLAST programs against the Ciona genomic DNA sequence database resulted in the identification of a remotely related gene. In situ hybridization analysis labeled primary sensory neurons in the rostral epithelium of amphioxus adults. Based on these findings, together with comparison of the developmental gene expression between amphioxus and vertebrates, I postulate that chemoreceptive primary sensory neurons in the rostrum are an ancient cell population traceable at least as far back in phylogeny as the common ancestor of amphioxus and vertebrates.

An adaptive radiation model for the origin of new gene functions

Current models for the evolution of new gene functions after gene duplication presume that the duplication events themselves have no effect on fitness. But those duplications that result in new gene functions are likely to be positively selected from their inception. The evolution of new function may start with the amplification of an existing gene with some level of preadaptation for that function, followed by a period of competitive evolution among the gene copies, resulting in the preservation of the most effective variant and the 'pseudogenization' and eventual loss of the rest.

A comparison of the human and chimpanzee olfactory receptor gene repertoires

Olfactory receptor (OR) genes constitute the basis of the sense of smell and are encoded by the largest mammalian gene superfamily, with >1000 members. In humans, but not in mice or dogs, the majority of OR genes have become pseudogenes, suggesting that OR genes in humans evolve under different selection pressures than in other mammals. To explore this further, we compare the OR gene repertoire of human with its closest living evolutionary relative, by taking advantage of the recently sequenced genome of the chimpanzee. In agreement with previous reports based on a small number of ORs, we find that humans have a significantly higher proportion of OR pseudogenes than chimpanzees. Moreover, we can reject the possibility that humans have been accumulating OR pseudogenes at a constant neutral rate since the divergence of human and chimpanzee. The comparison of the two repertoires reveals two chimpanzee-specific OR subfamily expansions and three expansions specific to humans. It also suggests that a subset of OR genes are under positive selection in either the human or the chimpanzee lineage. Thus, although overall there is relaxed constraint on human olfaction relative to chimpanzee, species-specific sensory requirements appear to have shaped the evolution of the functional OR gene repertoires in both species.

Characterization of the segmental duplication LCR7-20 in the human genome

Our previous study described the amplification of a genomic sequence containing exon 9 of CFTR in the human genome. Here we report that this CFTR sequence is part of a large duplicated sequence unit, provisionally named LCR7-20. Through successive screening of two human chromosome 7-specific cosmid libraries to construct a cosmid contig, we assembled two sequenced BAC clones into a single contig containing a prototypic LCR7-20 unit. Subsequent searches of existing human genome sequences identified additional six copies of LCR7-20-like sequences with more than 90% sequence homology. Additional genomic clones containing LCR7-20-like sequences were then isolated from total genomic BAC and PAC libraries. Restriction fragment analysis and limited sequencing data indicated that there could be around 30 copies of LCR7-20-like sequences in the human genome and that the average region of homology could extend over 120 kb. As indicated by fluorescence in situ hybridization analysis, LCR7-20-like sequences are dispersed on different chromosomes, mainly in the centromeric and pericentromeric regions, and some may exist in tandem copies. Our study also indicates that many genomic regions containing LCR7-20's either have been misassembled or are missing in current versions of the human genome sequence.

Reconstruction of the ancestral karyotype of eutherian mammals

Applying the parsimony principle, i.e. that chromosomes identical in species belonging to different taxa were likely to be present in their common ancestor, the ancestral karyotype of eutherian mammals (about 100 million years old) was tentatively reconstructed. Comparing chromosome banding with all ZOO-FISH data from literature or studied by us, this reconstruction can be proposed with only limited uncertainties. This karyotype comprised 50 chromosomes of which 40-42 were acrocentrics. Ten ancestral pairs of chromosomes were homologous to a single human chromosome: 5, 6, 9, 11, 13, 17, 18, 20, X and Y (human nomenclature). Nine others were homologous to a part of a human chromosome: 1p + q (proximal), 1q, 2p + q (proximal), 2q, part of 7, 8q, 10p, 10q and 19p (human nomenclature). Finally, seven pairs of chromosomes, homologs to human chromosomes 3 + 21, 4 + 8p, part of 7 + 16p, part of 12 + part of 22 (twice), 14 + 15, 16q + 19q, formed syntenies disrupted in man.

The human olfactory receptor gene family

Humans perceive an immense variety of chemicals as having distinct odors. Odor perception initiates in the nose, where odorants are detected by a large family of olfactory receptors (ORs). ORs have diverse protein sequences but can be assigned to subfamilies on the basis of sequence relationships. Members of the same subfamily have related sequences and are likely to recognize structurally related odorants. To gain insight into the mechanisms underlying odor perception, we analyzed the human OR gene family. By searching the human genome database, we identified 339 intact OR genes and 297 OR pseudogenes. Determination of their genomic locations showed that OR genes are unevenly distributed among 51 different loci on 21 human chromosomes. Sequence comparisons showed that the human OR family is composed of 172 subfamilies. Types of odorant structures that may be recognized by some subfamilies were predicted by identifying subfamilies that contain ORs with known odor ligands or human homologs of such ORs. Analysis of the chromosomal locations of members of each OR subfamily revealed that most subfamilies are encoded by a single chromosomal locus. Moreover, many loci encode only one or a few subfamilies, suggesting that different parts of the genome may, to some extent, be involved in the detection of different types of odorant structural motifs.

Loss of olfactory receptor genes coincides with the acquisition of full trichromatic vision in primates

Olfactory receptor (OR) genes constitute the molecular basis for the sense of smell and are encoded by the largest gene family in mammalian genomes. Previous studies suggested that the proportion of pseudogenes in the OR gene family is significantly larger in humans than in other apes and significantly larger in apes than in the mouse. To investigate the process of degeneration of the olfactory repertoire in primates, we estimated the proportion of OR pseudogenes in 19 primate species by surveying randomly chosen subsets of 100 OR genes from each species. We find that apes, Old World monkeys and one New World monkey, the howler monkey, have a significantly higher proportion of OR pseudogenes than do other New World monkeys or the lemur (a prosimian). Strikingly, the howler monkey is also the only New World monkey to possess full trichromatic vision, along with Old World monkeys and apes. Our findings suggest that the deterioration of the olfactory repertoire occurred concomitant with the acquisition of full trichromatic color vision in primates.

Examination of olfactory receptor genes in 19 primate species suggests that the olfactory repertoire lost complexity as our ancestors acquired full-color vision

Co-duplication of olfactory receptor and MHC class I genes in the mouse major histocompatibility complex

We report the 897 kb sequence of a cluster of olfactory receptor (OR) genes located at the distal end of the major histocompatibility complex (MHC) class I region on mouse chromosome 17 of strain 129/SvJ (H2bc). With additional information from the mouse genome draft sequence, we identified 59 OR loci (approximately 20% pseudogenes) in contrast to only 25 OR loci (approximately 50% pseudogenes) in the corresponding centromeric OR cluster that is part of the 'extended MHC class I region' on human chromosome 6. Comparative analysis leads to three major observations: (i) most of the OR subfamilies have evolved independently in the two species, expanding more in the mouse, and resulting in co-orthologs--subfamilies of highly similar paralogs that keep orthologous relationships with their human counterparts; (ii) three of the mouse OR subfamilies have no orthologs in humans; and (iii) MHC class I loci are interspersed in the OR cluster in mouse but not in human, and were subjected to co-duplication with OR genes. Screening of our sequence against the available sequences of other strains/haplotypes revealed that most of the OR loci are polymorphic and that the number of OR loci may vary among strains/haplotypes. Our findings that MHC-linked OR loci share duplication with MHC class I loci, have duplicated extensively and are polymorphic revives questions about potential reciprocal influences acting on the dynamics and evolution of the H2 region and the H2-linked OR loci.

The major shifts of human duplicated genes

Since many gene duplications in the human genome are ancient duplications going back to the origin of vertebrates, the question may be asked about the fate of such duplicated genes at the compositional genome transitions that occurred between cold- and warm-blooded vertebrates. Indeed, at that transition, about half of the (GC-poor) genes of cold-blooded vertebrates (the genes of the gene-dense "ancestral genome core") underwent a GC enrichment to become the genes of the "genome core" of warm-blooded vertebrates. Since the compositional distribution of the human duplicated genes investigated (1111 pairs) mimics the general distribution of human genes (about 50% GC(3)-poor and 50% GC(3)-rich genes, the border being at 60% GC(3)), we considered two possibilities, namely that the compositional transition affected either (i) about half of the copies on a random basis, or (ii) preferentially only one copy of the duplicated genes. The two possibilities could be distinguished if each copy is put into one of two subsets according to its GC(3) level. Indeed, in the first case, the two distributions would be similar, whereas in the second case, the two distributions would be different, one copy having maintained the ancestral GC-poor composition, and one copy having undergone the compositional change. Using this approach, we could show that, by far and large, one copy of the duplicated genes preferentially underwent the GC enrichment. This result implies that this copy, which had possibly acquired a different function and/or regulation, was preferentially translocated into the gene-dense compartment of the genome, the "ancestral genome core", namely the "gene space" which underwent the compositional transition at the emergence of warm-blooded vertebrates.

Natural selection on the olfactory receptor gene family in humans and chimpanzees

The olfactory receptor (OR) genes constitute the largest gene family in mammalian genomes. Humans have >1,000 OR genes, of which only approximately 40% have an intact coding region and are therefore putatively functional. In contrast, the fraction of intact OR genes in the genomes of the great apes is significantly greater (68%-72%), suggesting that selective pressures on the OR repertoire vary among these species. We have examined the evolutionary forces that shaped the OR gene family in humans and chimpanzees by resequencing 20 OR genes in 16 humans, 16 chimpanzees, and one orangutan. We compared the variation at the OR genes with that at intergenic regions. In both humans and chimpanzees, OR pseudogenes seem to evolve neutrally. In chimpanzees, patterns of variability are consistent with purifying selection acting on intact OR genes, whereas, in humans, there is suggestive evidence for positive selection acting on intact OR genes. These observations are likely due to differences in lifestyle, between humans and great apes, that have led to distinct sensory needs.

Complex low-copy repeats associated with a common polymorphic inversion at human chromosome 8p23

To characterize a submicroscopic, common 8p23 polymorphic inversion, we constructed a complete BAC/PAC-based physical map covering the entire 4.7-Mb inversion and its flanking regions. Two low-copy repeats (LCRs), REPD (approximately 1.3 Mb) and REPP (approximately 0.4 Mb), were identified at each of the inversion breakpoints. Comparison of the REPD and REPP sequences revealed that REPD showed high homology to REPP, with complex direct and inverted orientations. REPD and REPP contain six and five olfactory receptor gene-related sequences, respectively. LCRs at 8p23 showed multiple FISH signals from an Old World monkey to the human. Thus, multiplication of the LCR may have occurred at least 21-25 million years ago. We also investigated the frequency of the 4.7-Mb inversion in the general Japanese population and found that the allele frequency for the 8p23 inversion was estimated to be 27%.

The 200-kb segmental duplication on human chromosome 21 originates from a pericentromeric dissemination involving human chromosomes 2, 18 and 13

Regions close to human centromeres contain DNA fragments spanning hundreds of kilobases that exhibit a high degree of sequence identity (>95%). Here we report the genomic structure and evolution of a family of four paralogous regions related to a 220-kb genomic fragment present on the long arm of human chromosome 21 (21q22.1). Phylogenetic classification of the paralogous sequences obtained from the draft of the Human Genome Project are in agreement with results from comparative fluorescence in situ hybridization on metaphase chromosomes from human and great apes. The original copy present in 21q22.1 in human was duplicated in great apes after the divergence of the orang-utan and inserted in a pericentromeric region, most likely the ancestor of HSA2q, then disseminated by transposition of a larger fragment to other pericentromeric locations: HSA18p11, HSA13q11 and HSA21q11.1. The degree of dissemination varies among species.

From subgenome analysis to protein structure

Groups of related genes abound in large eukaryotic genomes. In such 'subgenomes', homology modeling carried out for a few genes will probably have relevance to the entire group. Subgenomes also afford unique ways of determining protein structural information. In addition to analyses based on the quantification of residue variability in paralogs, two-way comparisons, both within and among species, help to disclose functional amino acids. Comparative studies of gene families throughout the mammalian genome will also help elucidate the functional significance of single nucleotide polymorphisms in coding regions.

Complex evolution of 7E olfactory receptor genes in segmental duplications

Large segmental duplications (SDs) constitute at least 3.6% of the human genome and have increased its size, complexity, and diversity. SDs can mediate ectopic sequence exchange resulting in gross chromosomal rearrangements that could contribute to speciation and disease. We have identified and evaluated a subset of human SDs that harbor an 88-member subfamily of olfactory receptor (OR)-like genes called the 7Es. At least 92% of these genes appear to be pseudogenes when compared to other OR genes. The 7E-containing SDs (7E SDs) have duplicated to at least 35 regions of the genome via intra- and interchromosomal duplication events. In contrast to many human SDs, the 7E SDs are not biased towards pericentromeric or subtelomeric regions. We find evidence for gene conversion among 7E genes and larger sequence exchange between 7E SDs, supporting the hypothesis that long, highly similar stretches of DNA facilitate ectopic interactions. The complex structure and history of the 7E SDs necessitates extension of the current model of large-scale DNA duplication. Despite their appearance as pseudogenes, some 7E genes exhibit a signature of purifying selection, and at least one 7E gene is expressed.

Large-scale variation among human and great ape genomes determined by array comparative genomic hybridization

Large-scale genomic rearrangements are a major force of evolutionary change and the ascertainment of such events between the human and great ape genomes is fundamental to a complete understanding of the genetic history and evolution of our species. Here, we present the results of an evolutionary analysis utilizing array comparative genomic hybridization (array CGH), measuring copy-number gains and losses among these species. Using an array of 2460 human bacterial artificial chromosomes (BACs) (12% of the genome), we identified a total of 63 sites of putative DNA copy-number variation between humans and the great apes (chimpanzee, bonobo, gorilla, and orangutan). Detailed molecular characterization of a subset of these sites confirmed rearrangements ranging from 40 to at least 175 kb in size. Surprisingly, the majority of variant sites differentiating great ape and human genomes were found within interstitial euchromatin. These data suggest that such large-scale events are not restricted solely to subtelomeric or pericentromeric regions, but also occur within genic regions. In addition, 5/9 of the verified variant sites localized to areas of intrachromosomal segmental duplication within the human genome. On the basis of the frequency of duplication in humans, this represents a 14-fold positional bias. In contrast to previous cytogenetic and comparative mapping studies, these results indicate extensive local repatterning of hominoid chromosomes in euchromatic regions through a duplication-driven mechanism of genome evolution.

Genomic sequence and transcriptional profile of the boundary between pericentromeric satellites and genes on human chromosome arm 10p

Contiguous finished sequence from highly duplicated pericentromeric regions of human chromosomes is needed if we are to understand the role of pericentromeric instability in disease, and in gene and karyotype evolution. Here, we have constructed a BAC contig spanning the transition from pericentromeric satellites to genes on the short arm of human chromosome 10, and used this to generate 1.4 Mb of finished genomic sequence. Combining RT-PCR, in silico gene prediction, and paralogy analysis, we can identify two domains within the sequence. The proximal 600 kb consists of satellite-rich pericentromerically duplicated DNA which is transcript poor, containing only three unspliced transcripts. In contrast, the distal 850 kb contains four known genes (ZNF248, ZNF25, ZNF33A, and ZNF37A) and up to 32 additional transcripts of unknown function. This distal region also contains seven out of the eight intrachromosomal duplications within the sequence, including the p arm copy of the approximately 250-kb duplication which gave rise to ZNF33A and ZNF33B. By sequencing orthologs of the duplicated ZNF33 genes we have established that ZNF33A has diverged significantly at residues critical for DNA binding but ZNF33B has not, indicating that ZNF33B has remained constrained by selection for ancestral gene function. These results provide further evidence of gene formation within intrachromosomal duplications, but indicate that recent interchromosomal duplications at this centromere have involved transcriptionally inert, satellite rich DNA, which is likely to be heterochromatic. This suggests that any novel gene structures formed by these interchromosomal events would require relocation to a more open chromatin environment to be expressed.

Human Gene-Centric Databases at the Weizmann Institute of Science: GeneCards, UDB, CroW 21 and HORDE

Recent enhancements and current research in the GeneCards (GC) (http://bioinfo.weizmann.ac.il/cards/) project are described, including the addition of gene expression profiles and integrated gene locations. Also highlighted are the contributions of specialized associated human gene-centric databases developed at the Weizmann Institute. These include the Unified Database (UDB) (http://bioinfo.weizmann.ac.il/udb) for human genome mapping, the human Chromosome 21 database at the Weizmann Insti-tute (CroW 21) (http://bioinfo.weizmann.ac.il/crow21), and the Human Olfactory Receptor Data Explora-torium (HORDE) (http://bioinfo.weizmann.ac.il/HORDE). The synergistic relationships amongst these efforts have positively impacted the quality, quantity and usefulness of the GeneCards gene compendium.

Molecular mechanisms for genomic disorders

Genomic rearrangements play a major role in the pathogenesis of human genetic diseases. Nonallelic homologous recombination (NAHR) between low-copy repeats (LCRs) that flank unique genomic segments results in changes of genome organization and can cause a loss or gain of genomic segments. These LCRs appear to have arisen recently during primate speciation via paralogous segmental duplication, thus making the human species particularly susceptible to genomic rearrangements. Genomic disorders are defined as a group of diseases that result from genomic rearrangements, mostly mediated by NAHR. Molecular investigations of genomic disorders have revealed genome architectural features associated with susceptibility to rearrangements and the recombination mechanisms responsible for such rearrangements. The human genome sequence project reveals that LCRs may account for 5% of the genome, suggesting that many novel genomic disorders might still remain to be recognized.

Use of fluorescent sequence-specific polyamides to discriminate human chromosomes by microscopy and flow cytometry

In this paper, we demonstrate the use of synthetic polyamide probes to fluorescently label heterochromatic regions on human chromosomes for discrimination in cytogenetic preparations and by flow cytometry. Polyamides bind to the minor groove of DNA in a sequence-specific manner. Unlike conventional sequence-specific DNA or RNA probes, polyamides can recognize their target sequence without the need to subject chromosomes to harsh denaturing conditions. For this study, we designed and synthesized a polyamide to target the TTCCA-motif repeated in the heterochromatic regions of chromosome 9, Y and 1. We demonstrate that the fluorescently labeled polyamide binds to its target sequence in both conventional cytogenetic preparations of metaphase chromosomes and suspended chromosomes without denaturation. Chromosomes 9 and Y can be discriminated and purified by flow sorting on the basis of polyamide binding and Hoechst 33258 staining. We generate chromosome 9- and Y-specific 'paints' from the sorted fractions. We demonstrate the utility of this technology by characterizing the sequence of an olfactory receptor gene that is duplicated on multiple chromosomes. By separating chromosome 9 from chromosomes 10-12 on the basis of polyamide fluorescence, we determine and differentiate the haplotypes of the highly similar copies of this gene on chromosomes 9 and 11.

Haplotype exclusion: the unique case presented by multiple immunoglobulin gene loci in cartilaginous fish

Cartilaginous fish represent the most phylogenetically distant species from man in which immunoglobulin and T cell antigen receptor genes have been identified. Immunoglobulin genes in cartilaginous fish are organized in hundreds of clusters, located on different chromosomes and presumably are under independent regulation; large numbers of immunoglobulin gene clusters are germline-joined and thus their expression is not directly dependent on somatic rearrangement. Despite the unusual nature of immunoglobulin gene genetics in these species, preliminary characterization of the transcription products of immunoglobulin loci in single cell isolates is consistent with haplotype exclusion. Certain features of immunoglobulin gene organization and expression in cartilaginous fish are remarkably similar to that of odorant receptors and suggest that at the level of transcriptional regulation, at least two different mechanisms could exist that relate to haplotype exclusion.

G proteins and olfactory signal transduction

The olfactory system sits at the interface of the environment and the nervous system and is responsible for correctly coding sensory information from thousands of odorous stimuli. Many theories existed regarding the signal transduction mechanism that mediates this difficult task. The discovery that odorant transduction utilizes a unique variation (a novel family of G protein-coupled receptors) based upon a very common theme (the G protein-coupled adenylyl cyclase cascade) to accomplish its vital task emphasized the power and versatility of this motif. We now must understand the downstream consequences of this cascade that regulates multiple second messengers and perhaps even gene transcription in response to the initial interaction of ligand with G protein-coupled receptor.

Genome architecture, rearrangements and genomic disorders

An increasing number of human diseases are recognized to result from recurrent DNA rearrangements involving unstable genomic regions. These are termed genomic disorders, in which the clinical phenotype is a consequence of abnormal dosage of gene(s) located within the rearranged genomic fragments. Both inter- and intrachromosomal rearrangements are facilitated by the presence of region-specific low-copy repeats (LCRs) and result from nonallelic homologous recombination (NAHR) between paralogous genomic segments. LCRs usually span approximately 10-400 kb of genomic DNA, share >or= 97% sequence identity, and provide the substrates for homologous recombination, thus predisposing the region to rearrangements. Moreover, it has been suggested that higher order genomic architecture involving LCRs plays a significant role in karyotypic evolution accompanying primate speciation.

Mammalian chemosensory receptors

Chemosensory receptors are critical for the survival of many mammalian species, and their genes can comprise up to 1% of mammalian genomes. Odorant, taste, and vomeronasal receptors are being discovered and functionally characterized at a rapid pace which has been further accelerated by the availability of the human genome sequence. Five multigene families, consisting of >1,000 genes in the mouse, have been proposed to encode functional chemoreceptors. Although all of the chemoreceptor gene families encode G-protein coupled receptors, they are largely unrelated and uniquely specialized for the processing of different chemosensory modalities. Using members of the families as molecular probes, great insights are being gained into the different organizational strategies used by these sensory systems to encode information in both the periphery and the brain.

The human repertoire of odorant receptor genes and pseudogenes

The nose of Homo sapiens is a sophisticated chemical sensor. It is able to smell almost any type of volatile molecule, often at extraordinarily low concentrations, and can make fine perceptual discriminations between structurally related molecules. The diversity of odor recognition is mediated by odorant receptor (OR) genes, discovered in 1991 by Buck & Axel. OR genes form the largest gene families in mammalian genomes. A decade after their discovery, advances in the sequencing of the human genome have provided a first draft of the human OR repertoire: It consists of approximately 1000 sequences, residing in multiple clusters spread throughout the genome, with more than half being pseudogenes. Allelic variants are beginning to be recognized and may provide an opportunity for genotype-phenotype correlations. Here, I review the current knowledge of the human OR repertoire and summarize the limited information available regarding putative pheromone and taste receptors in humans.

Human paralogs of KIAA0187 were created through independent pericentromeric-directed and chromosome-specific duplication mechanisms

KIAA0187 is a gene of unknown function that maps to 10q11 and has been subject to recent duplication events. Here we analyze 18 human paralogs of this gene and show that paralogs of exons 14-23 were formed through satellite-associated pericentromeric-directed duplication, whereas paralogs of exons 1-9 were created via chromosome-specific satellite-independent duplications. In silico, Northern, and RT-PCR analyses indicate that nine paralogs are transcribed, including four in which KIAA0187 exons are spliced onto novel sequences. Despite this, no new genes appear to have been created by these events. The chromosome 10 paralogs map to 10q11, 10q22, 10q23.1, and 10q23.3, forming part of a complex family of chromosome-specific repeats that includes GLUD1, Cathepsin L, and KIAA1099 pseudogenes. Phylogenetic analyses and comparative FISH indicates that the 10q23.1 and 10q23.3 repeats were created in 10q11 and relocated by a paracentric inversion 13 to 27 Myr ago. Furthermore, the most recent duplications, involving the KIAA1099 pseudogenes, have largely been confined to 10q11. These results indicate a simple model for the evolution of this repeat family, involving multiple rounds of centromere-proximal duplication and dispersal through intrachromosomal rearrangement. However, more complex events must be invoked to account for high sequence identity between some paralogs.

Genomic analysis of the olfactory receptor region of the mouse and human T-cell receptor alpha/delta loci

We have conducted a comparative genomic analysis of several olfactory receptor (OR) genes that lie immediately 5' to the V-alpha gene segments at the mouse and human T-cell receptor (TCR) alpha/delta loci. Five OR genes are identified in the human cluster. The murine cluster has at least six OR genes; the first five are orthologous to the human genes. The sixth mouse gene has arisen since mouse-human divergence by a duplication of a approximately 10-kb block. One pair of OR paralogs found at the mouse and human loci are more similar to each other than to their corresponding orthologs. This paralogous "twinning" appears to be under selection, perhaps to increase sensitivity to particular odorants or to resolve structurally-similar odorants. The promoter regions of the mouse OR genes were identified by RACE-PCR. Orthologs share extensive 5' UTR homology, but we find no significant similarity among paralogs. These findings extend previous observations that suggest that OR genes do not share local significant regulatory homology despite having a common regulatory agenda. We also identified a diverged TCR-alpha gene segment that uses a divergent recombination signal sequence (RSS) to initiate recombination in T-cells from within the OR region. We explored the hypothesis that OR genes may use DNA recombination in expressing neurons, e.g., to recombine ORs into a transcriptionally active locus. We searched the mouse sequence for OR-flanking RSS motifs, but did not find evidence to suggest that these OR genes use TCR-like recombination target sequences.

Human-specific duplication and mosaic transcripts: the recent paralogous structure of chromosome 22

In recent decades, comparative chromosomal banding, chromosome painting, and gene-order studies have shown strong conservation of gross chromosome structure and gene order in mammals. However, findings from the human genome sequence suggest an unprecedented degree of recent (<35 million years ago) segmental duplication. This dynamism of segmental duplications has important implications in disease and evolution. Here we present a chromosome-wide view of the structure and evolution of the most highly homologous duplications (> or = 1 kb and > or = 90%) on chromosome 22. Overall, 10.8% (3.7/33.8 Mb) of chromosome 22 is duplicated, with an average sequence identity of 95.4%. To organize the duplications into tractable units, intron-exon structure and well-defined duplication boundaries were used to define 78 duplicated modules (minimally shared evolutionary segments) with 157 copies on chromosome 22. Analysis of these modules provides evidence for the creation or modification of 11 novel transcripts. Comparative FISH analyses of human, chimpanzee, gorilla, orangutan, and macaque reveal qualitative and quantitative differences in the distribution of these duplications--consistent with their recent origin. Several duplications appear to be human specific, including a approximately 400-kb duplication (99.4%-99.8% sequence identity) that transposed from chromosome 14 to the most proximal pericentromeric region of chromosome 22. Experimental and in silico data further support a pericentromeric gradient of duplications where the most recent duplications transpose adjacent to the centromere. Taken together, these data suggest that segmental duplications have been an ongoing process of primate genome evolution, contributing to recent gene innovation and the dynamic transformation of genome architecture within and among closely related species.