Whole-genome sequencing suggests a chemokine gene cluster that modifies age at onset in familial Alzheimer's disease.

We have sequenced the complete genomes of 72 individuals affected with early-onset familial Alzheimer's disease caused by an autosomal dominant, highly penetrant mutation in the presenilin-1 (PSEN1) gene, and performed genome-wide association testing to identify variants that modify age at onset (AAO) of Alzheimer's disease. Our analysis identified a haplotype of single-nucleotide polymorphisms (SNPs) on chromosome 17 within a chemokine gene cluster associated with delayed onset of mild-cognitive impairment and dementia. Individuals carrying this haplotype had a mean AAO of mild-cognitive impairment at 51.0 ± 5.2 years compared with 41.1 ± 7.4 years for those without these SNPs. This haplotype thus appears to modify Alzheimer's AAO, conferring a large (~10 years) protective effect. The associated locus harbors several chemokines including eotaxin-1 encoded by CCL11, and the haplotype includes a missense polymorphism in this gene. Validating this association, we found plasma eotaxin-1 levels were correlated with disease AAO in an independent cohort from the University of California San Francisco Memory and Aging Center. In this second cohort, the associated haplotype disrupted the typical age-associated increase of eotaxin-1 levels, suggesting a complex regulatory role for this haplotype in the general population. Altogether, these results suggest eotaxin-1 as a novel modifier of Alzheimer's disease AAO and open potential avenues for therapy.

Meta-analysis of gene expression in human pancreatic islets after in vitro expansion.

Pancreatic islet transplantation as a potential cure for type 1 diabetes (T1D) cannot be scaled up due to a scarcity of human pancreas donors. In vitro expansion of beta-cells from mature human pancreatic islets provides an alternative source of insulin-producing cells. The exact nature of the expanded cells produced by diverse expansion protocols and their potential for differentiation into functional beta-cells remain elusive. We performed a large-scale meta-analysis of gene expression in human pancreatic islet cells, which were processed using three different previously described protocols for expansion and for which redifferentiation was attempted. All three expansion protocols induced dramatic changes in the expression profiles of pancreatic islets; many of these changes are shared among the three protocols. Attempts at redifferentiation of expanded cells induce a limited number of gene expression changes. Nevertheless, these fail to restore a pancreatic islet-like gene expression pattern. Comparison with a collection of public microarray datasets confirmed that expanded cells are highly comparable to mesenchymal stem cells. Genes induced in expanded cells are also enriched for targets of transcription factors important for pluripotency induction. The present data increase our understanding of the active pathways in expanded and redifferentiated islets. Knowledge of the mesenchymal stem cell potential may help development of drug therapeutics to restore beta-cell mass in T1D patients.

The RUNX3 gene--sequence, structure and regulated expression.

The RUNX3 gene belongs to the runt domain family of transcription factors that act as master regulators of gene expression in major developmental pathways. In mammals the family includes three genes, RUNX1, RUNX2 and RUNX3. Here, we describe a comparative analysis of the human chromosome 1p36.1 encoded RUNX3 and mouse chromosome 4 encoded Runx3 genomic regions. The analysis revealed high similarities between the two genes in the overall size and organization and showed that RUNX3/Runx3 is the smallest in the family, but nevertheless exhibits all the structural elements characterizing the RUNX family. It also revealed that RUNX3/Runx3 bears a high content of the ancient mammalian repeat MIR. Together, these data delineate RUNX3/Runx3 as the evolutionary founder of the mammalian RUNX family. Detailed sequence analysis placed the two genes at a GC-rich H3 isochore with a sharp transition of GC content between the gene sequence and the downstream intergenic region. Two large conserved CpG islands were found within both genes, one around exon 2 and the other at the beginning of exon 6. RUNX1, RUNX2 and RUNX3 gene products bind to the same DNA motif, hence their temporal and spatial expression during development should be tightly regulated. Structure/function analysis showed that two promoter regions, designated P1 and P2, regulate RUNX3 expression in a cell type-specific manner. Transfection experiments demonstrated that both promoters were highly active in the GM1500 B-cell line, which endogenously expresses RUNX3, but were inactive in the K562 myeloid cell line, which does not express RUNX3.

Comparative genomics of the human and mouse T cell receptor loci.

The availability of the complete genomic sequences of the human and mouse T cell receptor loci opens up new opportunities for understanding T cell receptors (TCRs) and their genes. The full complement of TCR gene segments is finally known and should prove a valuable resource for supporting functional studies. A rational nomenclature system has been implemented and is widely available through IMGT and other public databases. Systematic comparisons of the genomic sequences within each locus, between loci, and across species enable precise analyses of the various diversification mechanisms and some regulatory signals. The genomic landscape of the TCR loci provides fundamental insights into TCR evolution as highly localized and tightly regulated gene families.

The complete human olfactory subgenome.

Olfactory receptors likely constitute the largest gene superfamily in the vertebrate genome. Here we present the nearly complete human olfactory subgenome elucidated by mining the genome draft with gene discovery algorithms. Over 900 olfactory receptor genes and pseudogenes (ORs) were identified, two-thirds of which were not annotated previously. The number of extrapolated ORs is in good agreement with previous theoretical predictions. The sequence of at least 63% of the ORs is disrupted by what appears to be a random process of pseudogene formation. ORs constitute 17 gene families, 4 of which contain more than 100 members each. "Fish-like" Class I ORs, previously considered a relic in higher tetrapods, constitute as much as 10% of the human repertoire, all in one large cluster on chromosome 11. Their lower pseudogene fraction suggests a functional significance. ORs are disposed on all human chromosomes except 20 and Y, and nearly 80% are found in clusters of 6-138 genes. A novel comparative cluster analysis was used to trace the evolutionary path that may have led to OR proliferation and diversification throughout the genome. The results of this analysis suggest the following genome expansion history: first, the generation of a "tetrapod-specific" Class II OR cluster on chromosome 11 by local duplication, then a single-step duplication of this cluster to chromosome 1, and finally an avalanche of duplication events out of chromosome 1 to most other chromosomes. The results of the data mining and characterization of ORs can be accessed at the Human Olfactory Receptor Data Exploratorium Web site (http://bioinfo.weizmann.ac.il/HORDE).

Architecture and anatomy of the genomic locus encoding the human leukemia-associated transcription factor RUNX1/AML1.

The RUNX1 gene on human chromosome 21q22.12 belongs to the 'runt domain' gene family of transcription factors (also known as AML/CBFA/PEBP2alpha). RUNX1 is a key regulator of hematopoiesis and a frequent target of leukemia associated chromosomal translocations. Here we present a detailed analysis of the RUNX1 locus based on its complete genomic sequence. RUNX1 spans 260 kb and its expression is regulated through two distinct promoter regions, that are 160 kb apart. A very large CpG island complex marks the proximal promoter (promoter-2), and an additional CpG island is located at the 3' end of the gene. Hitherto, 12 different alternatively spliced RUNX1 cDNAs have been identified. Genomic sequence analysis of intron/exon boundaries of these cDNAs has shown that all consist of properly spliced authentic coding regions. This indicates that the large repertoire of RUNX1 proteins, ranging in size between 20-52 kDa, are generated through usage of alternatively spliced exons some of which contain in frame stop codons. The gene's introns are largely depleted of repetitive sequences, especially of the LINE1 family. The RUNX1 locus marks the transition from a ~1 Mb of gene-poor region containing only pseudogenes, to a gene-rich region containing several functional genes. A search for RUNX1 sequences that may be involved in the high frequency of chromosomal translocations revealed that a 555 bp long segment originating in chromosome 11 FLI1 gene was transposed into RUNX1 intron 4.1. This intron harbors the t(8;21) and t(3;21) chromosomal breakpoints involved in acute myeloid leukemia. Interestingly, the FLI1 homologous sequence contains a breakpoint of the t(11;22) translocation associated with Ewing's tumors, and may have a similar function in RUNX1.

The human olfactory subgenome: from sequence to structure and evolution.

Olfactory receptors (ORs) constitute the largest multigene family in multicellular organisms. Their evolutionary proliferation has been driven by the need to provide recognition capacity for millions of potential odorants with arbitrary chemical configurations. Human genome sequencing has provided a highly informative picture of the "olfactory subgenome", the repertoire of OR genes. We describe here an analysis of 224 human OR genes, a much larger number than hitherto systematically analyzed. These are derived by literature survey, data mining at 14 genomic clusters, and by an OR-targeted experimental sequencing strategy. The presented set contains at least 53% pseudogenes and is minimally divided into 11 gene families. One of these (no. 7) has undergone a particularly extensive expansion in primates. The analysis of this collection leads to insight into the origin of OR genes, suggesting a graded expansion through mammalian evolution. It also allows us to delineate a structural map of the respective proteins. A sequence database and analysis package is provided (http://bioinformatics.weizmann.ac.il/HORDE), which will be useful for analyzing human OR sequences genome-wide.

Harvesting the human genome: the Israeli perspective.

The post-genome era is at our door, and soon the complete human genome sequence will be available for the next set of goals. Israel is well equipped and skilled to join the worldwide harvest of the human genome, but additional massive government investment is required. This will affect various domains of activity, including the fields of diagnostics and therapeutics. The technologies and know-how described above constitute the basis for future human genome applications in Israel.

The olfactory receptor gene superfamily: data mining, classification, and nomenclature.

The vertebrate olfactory receptor (OR) subgenome harbors the largest known gene family, which has been expanded by the need to provide recognition capacity for millions of potential odorants. We implemented an automated procedure to identify all OR coding regions from published sequences. This led us to the identification of 831 OR coding regions (including pseudogenes) from 24 vertebrate species. The resulting dataset was subjected to neighbor-joining phylogenetic analysis and classified into 32 distinct families, 14 of which include only genes from tetrapodan species (Class II ORs). We also report here the first identification of OR sequences from a marsupial (koala) and a monotreme (platypus). Analysis of these OR sequences suggests that the ancestral mammal had a small OR repertoire, which expanded independently in all three mammalian subclasses. Classification of "fish-like" (Class I) ORs indicates that some of these ancient ORs were maintained and even expanded in mammals. A nomenclature system for the OR gene superfamily is proposed, based on a divergence evolutionary model. The nomenclature consists of the root symbol 'OR', followed by a family numeral, subfamily letter(s), and a numeral representing the individual gene within the subfamily. For example, OR3A1 is an OR gene of family 3, subfamily A, and OR7E12P is an OR pseudogene of family 7, subfamily E. The symbol is to be preceded by a species indicator. We have assigned the proposed nomenclature symbols for all 330 human OR genes in the database. A WWW tool for automated name assignment is provided.

The genomic structure of human olfactory receptor genes.

The genomic and cDNA structures were studied for eight human olfactory receptor (OR) genes within the chromosome 17p13.3 cluster. A common gene structure was revealed, which included an approximately 1-kb intronless coding region terminated by a signal for polyadenylation and a variable number of upstream noncoding exons. The latter were found to be alternatively spliced, giving rise to different isoforms of OR mRNA. While the initial exons mostly agreed with previous computer predictions and were conserved within OR subfamilies, other upstream exons were novel and idiosyncratic. In some cases, repetitive sequences were involved in the generation of splice sites and putative transcription control elements. Such gene structure is consistent with early repertoire enhancement by retrogene generation, which was likely followed by extensive genomic duplication. Each OR gene had a unique signature of transcription factor elements, consistent with a combinatorial expression control mechanism.

Identification of the gene causing mucolipidosis type IV.

Mucolipidosis type IV (MLIV) is an autosomal recessive, neurodegenerative, lysosomal storage disorder characterized by psychomotor retardation and ophthalmological abnormalities including corneal opacities, retinal degeneration and strabismus. Most patients reach a maximal developmental level of 12?15 months. The disease was classified as a mucolipidosis following observations by electron microscopy indicating the lysosomal storage of lipids together with water-soluble, granulated substances. Over 80% of the MLIV patients diagnosed are Ashkenazi Jews, including severely affected and mildly affected patients. The gene causing MLIV was previously mapped to human chromosome 19p13.2-13.3 in a region of approximately 1 cM (ref. 7). Haplotype analysis in the MLIV gene region of over 70 MLIV Ashkenazi chromosomes indicated the existence of two founder chromosomes among 95% of the Ashkenazi MLIV families: a major haplotype in 72% and a minor haplotype in 23% of the MLIV chromosomes (ref. 7, and G.B., unpublished data). The remaining 5% are distinct haplotypes found only in single patients. The basic metabolic defect causing the lysosomal storage in MLIV has not yet been identified. Thus, positional cloning was an alternative to identify the MLIV gene. We report here the identification of a new gene in this human chromosomal region in which MLIV-specific mutations were identified.

GESTALT: a workbench for automatic integration and visualization of large-scale genomic sequence analyses.

SUMMARY: The GESTALT Workbench is a WWW-based tool for genomic sequence analysis, comparison and annotation, with strong emphasis on visualization. GESTALT integrates graphically the output of diverse sequence analysis algorithms producing an information-rich, interactive genomic map. AVAILABILITY: The GESTALT Workbench, as well as a more detailed description, are available at http://bioinfo. weizmann.ac.il/GESTALT/.

Sequence, structure, and evolution of a complete human olfactory receptor gene cluster.

The olfactory receptor (OR) gene cluster on human chromosome 17p13.3 was subjected to mixed shotgun automated DNA sequencing. The resulting 412 kb of genomic sequence include 17 OR coding regions, 6 of which are pseudogenes. Six of the coding regions were discovered only upon genomic sequencing, while the others were previously reported as partial sequences. A comparison of DNA sequences in the vicinity of the OR coding regions revealed a common gene structure with an intronless coding region and at least one upstream noncoding exon. Potential gene control regions including specific pyrimidine:purine tracts and Olf-1 sites have been identified. One of the pseudogenes apparently has evolved into a CpG island. Four extensive CpG islands can be discerned within the cluster, not coupled to specific OR genes. The cluster is flanked at its telomeric end by an unidentified open reading frame (C17orf2) with no significant similarity to any known protein. A high proportion of the cluster sequence (about 60%) belongs to various families of interspersed repetitive elements, with a clear predominance of LINE repeats. The OR genes in the cluster belong to two families and seven subfamilies, which show a relatively high degree of intermixing along the cluster, in seemingly random orientations. This genomic organization may be best accounted for by a complex series of evolutionary events.

Identification and characterization of coding single-nucleotide polymorphisms within a human olfactory receptor gene cluster.

Single-nucleotide polymorphisms (SNPs) were studied in 15 olfactory receptor (OR) coding regions, one control region and two noncoding sequences all residing within a 412 kb OR gene cluster on human chromosome 17p13.3, as well as in other G-protein coupled receptors (GPCRs). A total of 26 SNPs were identified in ORs, 21 of which are coding SNPs (cSNPs). The mean nucleotide diversity of OR coding regions was 0.078% (ranging from 0 to 0.16%), which is about twice higher than that of other GPCRs, and similar to the nucleotide diversity levels of noncoding regions along the human genome. The high polymorphism level in the OR coding regions might be due to a weak positive selection pressure acting on the OR genes. In two cases, OR genes have been found to share the same cSNP. This could be explained by recent gene conversion events, which might be a part of a concerted evolution mechanism acting on the OR superfamily. Using the genotype data of 85 unrelated individuals in 15 SNPs, we found linkage disequilibrium (LD) between pairs of SNPs located on the centromeric part of the cluster. On the other hand, no LD was found between SNPs located on the telomeric part of the cluster, suggesting the presence of several hot-spots for recombination within this cluster. Thus, different regions of this gene cluster may have been subject to different recombination rates.

Primate evolution of an olfactory receptor cluster: diversification by gene conversion and recent emergence of pseudogenes.

The olfactory receptor (OR) subgenome harbors the largest known gene family in mammals, disposed in clusters on numerous chromosomes. We have carried out a comparative evolutionary analysis of the best characterized genomic OR gene cluster, on human chromosome 17p13. Fifteen orthologs from chimpanzee (localized to chromosome 19p15), as well as key OR counterparts from other primates, have been identified and sequenced. Comparison among orthologs and paralogs revealed a multiplicity of gene conversion events, which occurred exclusively within OR subfamilies. These appear to lead to segment shuffling in the odorant binding site, an evolutionary process reminiscent of somatic combinatorial diversification in the immune system. We also demonstrate that the functional mammalian OR repertoire has undergone a rapid decline in the past 10 million years: while for the common ancestor of all great apes an intact OR cluster is inferred, in present-day humans and great apes the cluster includes nearly 40% pseudogenes.

Organization and evolution of olfactory receptor genes on human chromosome 11.

Olfactory receptors (OR) are encoded by a large multigene family including hundreds of members dispersed throughout the human genome. Cloning and mapping studies have determined that a large proportion of the olfactory receptor genes are located on human chromosomes 6, 11, and 17, as well as distributed on other chromosomes. In this paper, we describe and characterize the organization of olfactory receptor genes on human chromosome 11 by using degenerate PCR-based probes to screen chromosome 11-specific and whole genome clone libraries for members of the OR gene family. OR genes were identified by DNA sequencing and then localized to regions of chromosome 11. Physical maps of several gene clusters were constructed to determine the chromosomal relationships between various members of the family. This work identified 25 new OR genes located on chromosome 11 in at least seven distinct regions. Three of these regions contain gene clusters that include additional members of this gene family not yet identified by sequencing. Phylogenetic analysis of the newly described OR genes suggests a mechanism for the generation of genetic diversity.

Genome dynamics, evolution, and protein modeling in the olfactory receptor gene superfamily.

The human olfactory subgenome represents several hundred olfactory receptor (OR) genes in a dozen or more clusters on several chromosomes. One OR gene cluster on human chromosome 17 has been characterized by us in detail. Based on a large-scale DNA sequence analysis, we have identified events of gene duplication and fusion as well as the generation of pseudogenes. The latter instances of 'gene death' could underlie the widespread phenomenon of human specific anosmias. Sixteen OR coding regions were found on this cluster, and six of them are pseudogenes. One of these pseudogenes, OR17-23, was found to be an intact open reading frame in an old world monkey. This may be a reflection of an OR repertoire diminution in man. A homology model of the OR protein was constructed by utilizing the rich information available on approximately 200 OR sequences. The putative odorant complementarity determining regions (CDR) was found to consist of 20 hypervariable residues facing an interior caving defined by transmembrane helices 3, 4 and 5. Such a model could be useful in analyzing additional OR gene sequences in the human genome in terms of odorant binding.