7E olfactory receptor gene clusters and evolutionary chromosome rearrangements

Emergence and influence of sequence bias in evolutionarily malleable, mammalian tandem arrays

BACKGROUND

The radiation of mammals at the extinction of the dinosaurs produced a plethora of new forms-as diverse as bats, dolphins, and elephants-in only 10-20 million years. Behind the scenes, adaptation to new niches is accompanied by extensive innovation in large families of genes that allow animals to contact the environment, including chemosensors, xenobiotic enzymes, and immune and barrier proteins. Genes in these "outward-looking" families are allelically diverse among humans and exhibit tissue-specific and sometimes stochastic expression.

RESULTS

Here, we show that these tandem arrays of outward-looking genes occupy AT-biased isochores and comprise the "tissue-specific" gene class that lack CpG islands in their promoters. Models of mammalian genome evolution have not incorporated the sharply different functions and transcriptional patterns of genes in AT- versus GC-biased regions. To examine the relationship between gene family expansion, sequence content, and allelic diversity, we use population genetic data and comparative analysis. First, we find that AT bias can emerge during evolutionary expansion of gene families in cis. Second, human genes in AT-biased isochores or with GC-poor promoters experience relatively low rates of de novo point mutation today but are enriched for non-synonymous variants. Finally, we find that isochores containing gene clusters exhibit low rates of recombination.

CONCLUSIONS

Our analyses suggest that tolerance of non-synonymous variation and low recombination are two forces that have produced the depletion of GC bases in outward-facing gene arrays. In turn, high AT content exerts a profound effect on their chromatin organization and transcriptional regulation.

Expression profile of ectopic olfactory receptors determined by deep sequencing

Olfactory receptors (ORs) provide the molecular basis for the detection of volatile odorant molecules by olfactory sensory neurons. The OR supergene family encodes G-protein coupled proteins that belong to the seven-transmembrane-domain receptor family. It was initially postulated that ORs are exclusively expressed in the olfactory epithelium. However, recent studies have demonstrated ectopic expression of some ORs in a variety of other tissues. In the present study, we conducted a comprehensive expression analysis of ORs using an extended panel of human tissues. This analysis made use of recent dramatic technical developments of the so-called Next Generation Sequencing (NGS) technique, which encouraged us to use open access data for the first comprehensive RNA-Seq expression analysis of ectopically expressed ORs in multiple human tissues. We analyzed mRNA-Seq data obtained by Illumina sequencing of 16 human tissues available from Illumina Body Map project 2.0 and from an additional study of OR expression in testis. At least some ORs were expressed in all the tissues analyzed. In several tissues, we could detect broadly expressed ORs such as OR2W3 and OR51E1. We also identified ORs that showed exclusive expression in one investigated tissue, such as OR4N4 in testis. For some ORs, the coding exon was found to be part of a transcript of upstream genes. In total, 111 of 400 OR genes were expressed with an FPKM (fragments per kilobase of exon per million fragments mapped) higher than 0.1 in at least one tissue. For several ORs, mRNA expression was verified by RT-PCR. Our results support the idea that ORs are broadly expressed in a variety of tissues and provide the basis for further functional studies.

Comparative genomics reveals birth and death of fragile regions in mammalian evolution

Background

An important question in genome evolution is whether there exist fragile regions (rearrangement hotspots) where chromosomal rearrangements are happening over and over again. Although nearly all recent studies supported the existence of fragile regions in mammalian genomes, the most comprehensive phylogenomic study of mammals raised some doubts about their existence.

Results

Here we demonstrate that fragile regions are subject to a birth and death process, implying that fragility has a limited evolutionary lifespan.

Conclusions

This finding implies that fragile regions migrate to different locations in different mammals, explaining why there exist only a few chromosomal breakpoints shared between different lineages. The birth and death of fragile regions as a phenomenon reinforces the hypothesis that rearrangements are promoted by matching segmental duplications and suggests putative locations of the currently active fragile regions in the human genome.

Breakpoint graphs and ancestral genome reconstructions

Recently completed whole-genome sequencing projects marked the transition from gene-based phylogenetic studies to phylogenomics analysis of entire genomes. We developed an algorithm MGRA for reconstructing ancestral genomes and used it to study the rearrangement history of seven mammalian genomes: human, chimpanzee, macaque, mouse, rat, dog, and opossum. MGRA relies on the notion of the multiple breakpoint graphs to overcome some limitations of the existing approaches to ancestral genome reconstructions. MGRA also generates the rearrangement-based characters guiding the phylogenetic tree reconstruction when the phylogeny is unknown.

Update on the olfactory receptor (OR) gene superfamily

The olfactory receptor gene (OR) superfamily is the largest in the human genome. The superfamily contains 390 putatively functional genes and 465 pseudogenes arranged into 18 gene families and 300 subfamilies. Even members within the same subfamily are often located on different chromosomes. OR genes are located on all autosomes except chromosome 20, plus the X chromosome but not the Y chromosome. The gene:pseudogene ratio is lowest in human, higher in chimpanzee and highest in rat and mouse -- most likely reflecting the greater need of olfaction for survival in the rodent than in the human. The OR genes undergo allelic exclusion, each sensory neurone expressing usually only one odourant receptor allele; the mechanism by which this phenomenon is regulated is not yet understood. The nomenclature system (based on evolutionary divergence of genes into families and subfamilies of the OR gene superfamily) has been designed similarly to that originally used for the CYP gene superfamily.

Physical map of two tammar wallaby chromosomes: a strategy for mapping in non-model mammals

Marsupials are especially valuable for comparative genomic studies of mammals. Two distantly related model marsupials have been sequenced: the South American opossum (Monodelphis domestica) and the tammar wallaby (Macropus eugenii), which last shared a common ancestor about 70 Mya. The six-fold opossum genome sequence has been assembled and assigned to chromosomes with the help of a cytogenetic map. A good cytogenetic map will be even more essential for assembly and anchoring of the two-fold wallaby genome. As a start to generating a physical map of gene locations on wallaby chromosomes, we focused on two chromosomes sharing homology with the human X, wallaby chromosomes X and 5. We devised an efficient strategy for mapping large conserved synteny blocks in non-model mammals, and applied this to generate dense maps of the X and 'neo-X' regions and to determine the arrangement of large conserved synteny blocks on chromosome 5. Comparisons between the wallaby and opossum chromosome maps revealed many rearrangements, highlighting the need for comparative gene mapping between South American and Australian marsupials. Frequent rearrangement of the X, along with the absence of a marsupial XIST gene, suggests that inactivation of the marsupial X chromosome does not depend on a whole-chromosome repression by a control locus.

DA and Xiao-two giant and composite LTR-retrotransposon-like elements identified in the human genome

We discovered two new complex elements while studying large genomic rearrangements and segmental duplications in the human genome. Both resemble bacterial composite DNA transposon Tn9, consisting of a core flanked by mobile elements, except that the flanking element is not a DNA transposon but instead is long terminal repeat retrotransposon-like with human endogenous retrovirus and satellite sequences. Based on the core size, we named them Xiao ( approximately 30 kb) and DA ( approximately 280 kb), meaning small and big, respectively, in Chinese. Xiao originated from a 19p region encoding olfactory receptor 7E members after the human/ape divergence from Old World monkeys, while DA likely evolved from a Xiao by inserting approximately 200 kb of chimeric sequence from 16p and 21q into the Xiao core, resulting in a target site duplication of 3.4 kb. DA/Xiao was identified in 30 loci on 12 chromosomes, and only DAs mediated intrachromosomal rearrangements, based on our reconstructed human-mouse-rat ancestral genome and the rhesus macaque genome.

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.

Segmental duplications and evolutionary plasticity at tumor chromosome break-prone regions

We have previously found that the borders of evolutionarily conserved chromosomal regions often coincide with tumor-associated deletion breakpoints within human 3p12-p22. Moreover, a detailed analysis of a frequently deleted region at 3p21.3 (CER1) showed associations between tumor breaks and gene duplications. We now report on the analysis of 54 chromosome 3 breaks by multipoint FISH (mpFISH) in 10 carcinoma-derived cell lines. The centromeric region was broken in five lines. In lines with highly complex karyotypes, breaks were clustered near known fragile sites, FRA3B, FRA3C, and FRA3D (three lines), and in two other regions: 3p12.3-p13 ( approximately 75 Mb position) and 3q21.3-q22.1 ( approximately 130 Mb position) (six lines). All locations are shown based on NCBI Build 36.1 human genome sequence. The last two regions participated in three of four chromosome 3 inversions during primate evolution. Regions at 75, 127, and 131 Mb positions carry a large ( approximately 250 kb) segmental duplication (tumor break-prone segmental duplication [TBSD]). TBSD homologous sequences were found at 15 sites on different chromosomes. They were located within bands frequently involved in carcinoma-associated breaks. Thirteen of them have been involved in inversions during primate evolution; 10 were reused by breaks during mammalian evolution; 14 showed copy number polymorphism in man. TBSD sites showed an increase in satellite repeats, retrotransposed sequences, and other segmental duplications. We propose that the instability of these sites stems from specific organization of the chromosomal region, associated with location at a boundary between different CG-content isochores and with the presence of TBSDs and "instability elements," including satellite repeats and retroviral sequences.

Autosomal-dominant microtia linked to five tandem copies of a copy-number-variable region at chromosome 4p16

Recently, large-scale benign copy-number variations (CNVs)--encompassing over 12% of the genome and containing genes considered to be dosage tolerant for human development--were uncovered in the human population. Here we present a family with a novel autosomal-dominantly inherited syndrome characterized by microtia, eye coloboma, and imperforation of the nasolacrimal duct. This phenotype is linked to a cytogenetically visible alteration at 4pter consisting of five copies of a copy-number-variable region, encompassing a low-copy repeat (LCR)-rich sequence. We demonstrate that the approximately 750 kb amplicon occurs in exact tandem copies. This is the first example of an amplified CNV associated with a Mendelian disorder, a discovery that implies that genome screens for genetic disorders should include the analysis of so-called benign CNVs and LCRs.

Are there rearrangement hotspots in the human genome?

In a landmark paper, Nadeau and Taylor [18] formulated the random breakage model (RBM) of chromosome evolution that postulates that there are no rearrangement hotspots in the human genome. In the next two decades, numerous studies with progressively increasing levels of resolution made RBM the de facto theory of chromosome evolution. Despite the fact that RBM had prophetic prediction power, it was recently refuted by Pevzner and Tesler [4], who introduced the fragile breakage model (FBM), postulating that the human genome is a mosaic of solid regions (with low propensity for rearrangements) and fragile regions (rearrangement hotspots). However, the rebuttal of RBM caused a controversy and led to a split among researchers studying genome evolution. In particular, it remains unclear whether some complex rearrangements (e.g., transpositions) can create an appearance of rearrangement hotspots. We contribute to the ongoing debate by analyzing multi-break rearrangements that break a genome into multiple fragments and further glue them together in a new order. In particular, we demonstrate that (1) even if transpositions were a dominant force in mammalian evolution, the arguments in favor of FBM still stand, and (2) the "gene deletion" argument against FBM is flawed.

Modeling non-random deletions in cancer

Chromosome deletions do abound in cancer and are detected in certain regions in a non-random manner. Although their relevance remains elusive, it is a general agreement that segmental losses provide the cell with selective growth advantage. Consequently these may contain genes and/or regulatory sequences that control normal growth and inhibit malignancy. We have developed a monochromosomal hybrid based experimental model for the generation and functional analysis of deletions, that is called "elimination test" (Et). Focused on human chromosome 3 - that was known to carry multiple 3p deletions - the Et was expected to restrict a 3p tumor suppressor region to a sufficiently small segment that permits the selection of a critically important candidate gene. Surprisingly, we detected three regions that were lost in all or majority of tumors: CER1 (3p21.3, Mb: 43.32-45.74), CER2 (3p22, Mb: 37.83-39.06) and FER (3p14.3-p21.2, Mb: 50.12-58.03). In contrast a 3q26-qter region (CRR) was regularly retained. CER1 - our main focus - contains multiple genes that may inhibit tumor growth, but 3 genes, RIS1, LF (LTF) and LIMD1 have already the necessary experimental support to be considered bona fide tumor suppressors. Tumor suppressor region borders display instability features including: (1) they break in evolution and in tumors, (2) they evolve horizontally, and (3) they are enriched with pseudogene insertions. The most remarkable features at the breakpoint cluster regions were segmental duplications that drive horizontal evolution and contribute to cancer associated instability.

High-resolution mapping identifies a commonly amplified 11q13.3 region containing multiple genes flanked by segmental duplications

DNA amplification of the 11q13 region is observed frequently in many carcinomas. Within the amplified region several candidate oncogenes have been mapped, including cyclin D1, TAOS1 and cortactin. Yet, it is unknown which gene(s) is/are responsible for the selective pressure enabling amplicon formation. This is probably due to the use of low-resolution detection methods. Furthermore, the size and structure of the amplified 11q13 region is complex and consists of multiple amplicon cores that differ between different tumor types. We set out to test whether the borders of the 11q13 amplicon are restricted to regions that enable DNA breakage and subsequent amplification. A high-resolution array of the 11q13 region was generated to study the structure of the 11q13 amplicon and analyzed 29 laryngeal and pharyngeal carcinomas and nine cell lines with 11q13 amplification. We found that boundaries of the commonly amplified region were restricted to four segments. Three boundaries coincided with a syntenic breakpoint. Such regions have been suggested to be putatively fragile. Sequence comparisons revealed that the amplicon was flanked by two large low copy repeats known as segmental duplications. These segmental duplications might be responsible for the typical structure and size of the 11q13 amplicon. We hypothesize that the selection for genes through amplification of the 11q13.3 region is determined by the ability to form DNA breaks within specific regions and, consequently, results in large amplicons containing multiple genes.