Replication timing of extremely large genes on human chromosomes 11q and 21q

Many human genes have been mapped precisely in the genome. These genes vary from a few kb to more than 1 Mb in length. Previously, we measured replication timing along the entire lengths of human chromosomes 11q and 21q at the sequence level. In the present study, we used the newest information for human chromosomes 11q and 21q to analyze the replication timing of 30 extremely large genes (>250 kb) in two human cell lines (THP-1 and Jurkat). The timing of replication differed between the 5'- and 3'-ends of each of extremely large genes on 11q and 21q, and the time interval between their replication varied among genes of different lengths. The large genes analyzed here included several tissue-specific genes associated with neural diseases and genes encoding cell adhesion molecules: some of these genes had different patterns of replication timing between the two cell lines. The amyloid precursor protein gene (APP), which is associated with familial Alzheimer's disease (AD1), showed the largest difference in timing of replication between its 5'- and 3'-ends in relation to gene length of all the large genes studied on 11q and 21q. These extremely large genes were concentrated in and around genomic regions in which replication timing switches from early to late on both 11q and 21q. The differences of replication timing between the 5'- and 3'-terminal regions of large genes may be related to the molecular mechanisms that underlie tissue-specific expression.     

Replication timing in a single human chromosome 11 transferred into the Chinese hamster ovary (CHO) cell line

DNA replication in eukaryotes initiates from discrete genomic regions, termed origins, according to a strict and often tissue-specific temporal program. However, the genetic program that controls activation of replication origins has still not been fully elucidated in mammalian cells. Previously, we measured replication timing at the sequence level along human chromosomes 11q and 21q. In the present study, we sought to obtain a greater understanding of the relationship between replication timing programs and human chromosomes by analysis of the timing of replication of a single human chromosome 11 that had been transferred into the Chinese hamster ovary (CHO) cell line by chromosome engineering. Timing of replication was compared for three 11q chromosomal regions in the transformed CHO cell line (CHO(h11)) and the original human fibroblast cell line, namely, the R/G-band boundary at 11q13.5/q14.1, the centromere and the distal telomere. We found that the pattern of replication timing in and around the R/G band boundary at 11q13.5/q14.1 was similar in CHO(h11) cells and fibroblasts. The 11q centromeric region, which replicates late in human fibroblasts, replicated in the second half of S phase in CHO(h11) cells. By contrast, however, the telomeric region at 11q25, which is late replicating in fibroblasts (and in several other human cell lines), replicated in the first half of S phase or in very early S phase in CHO(h11) cells. Our observations suggest that the replication timing programs of the R/G-band boundary and the centromeric region of human chromosome 11q are maintained in CHO(h11) cells, whereas that for the telomeric region is altered. The replication timing program of telomeric regions on human chromosomes might be regulated by specific mechanisms that differ from those for other chromosomal regions.     

Cell line differences in replication timing of human glutamate receptor genes and other large genes associated with neural disease

There is considerable current interest in the function of epigenetic mechanisms in neuroplasticity with regard to learning and memory formation and to a range of neural diseases. Previously, we described replication timing on human chromosome 21q in the THP-1 human cell line (2n = 46, XY) and showed that several genes associated with neural diseases, such as the neuronal glutamate receptor subunit GluR-5 (GRIK1) and amyloid precursor protein (APP), were located in regions where replication timing transitioned from early to late S phase. Here, we compared replication timing of all known human glutamate receptor genes (26 genes in total) and APP in 6 different human cell lines including human neuron-related cell lines. Replication timings were obtained by integrating our previously reported data with new data generated here and information from the online database ReplicationDomain. We found that many of the glutamate receptor genes were clearly located in replication timing transition zones in neural precursor cells, but this relationship was less clear in embryonic stem cells before neural differentiation; in the latter, the genes were often located in later replication timing zones that displayed DNA hypermethylation. Analysis of selected large glutamate receptor genes (>200 kb), and of APP, showed that their precise replication timing patterns differed among the cell lines. We propose that the transition zones of DNA replication timing are altered by epigenetic mechanisms, and that these changes may affect the neuroplasticity that is important to memory and learning, and may also have a role in the development of neural diseases.     

Cell line differences in replication timing of human glutamate receptor genes and other large genes associated with neural disease

There is considerable current interest in the function of epigenetic mechanisms in neuroplasticity with regard to learning and memory formation and to a range of neural diseases. Previously, we described replication timing on human chromosome 21q in the THP-1 human cell line (2n = 46, XY) and showed that several genes associated with neural diseases, such as the neuronal glutamate receptor subunit GluR-5 (GRIK1) and amyloid precursor protein (APP), were located in regions where replication timing transitioned from early to late S phase. Here, we compared replication timing of all known human glutamate receptor genes (26 genes in total) and APP in 6 different human cell lines including human neuron-related cell lines. Replication timings were obtained by integrating our previously reported data with new data generated here and information from the online database ReplicationDomain. We found that many of the glutamate receptor genes were clearly located in replication timing transition zones in neural precursor cells, but this relationship was less clear in embryonic stem cells before neural differentiation; in the latter, the genes were often located in later replication timing zones that displayed DNA hypermethylation. Analysis of selected large glutamate receptor genes (>200 kb), and of APP, showed that their precise replication timing patterns differed among the cell lines. We propose that the transition zones of DNA replication timing are altered by epigenetic mechanisms, and that these changes may affect the neuroplasticity that is important to memory and learning, and may also have a role in the development of neural diseases.     

Evolutionarily plastic regions at human 3p21.3 coincide with tumor breakpoints identified by the "elimination test"

We have previously found with the microcell hybrid-based "elimination test" that human chromosome 3 transferred into murine or human tumor cells regularly lost certain 3p regions during tumor growth in SCID mice. The most common eliminated region, CER1, is approximately 2.4 Mb at 3p21.3. CER1 breakpoints were clustered in approximately 200-kb regions at both telomeric and centromeric borders. We have also shown, earlier, that tumor-related deletions often coincide with human/mouse synteny breakpoints on 3p12-p22. Here we describe the results of a comparative genomic analysis on the CER1 region in Caenorhabditis elegans, Drosophila melanogaster, Fugu rubripes, Gallus gallus, Mus musculus, Rattus norvegicus, and Canis familiaris. First, four independent synteny breaks were found within the CER1 telomeric breakpoint cluster region, comparing human, dog, and chicken genomes, and two independent synteny breaks within the CER1 centromeric breakpoint cluster region, comparing human, mouse, and chicken genomes, suggesting a nonrandom involvement of tumor breakpoint regions in chromosome evolution. Second, both CER1 breakpoint cluster regions show recent tandem duplications (seven Zn finger protein family genes at the telomeric and eight chemokine receptor genes at the centromeric side). Finally, all genes from these regions underwent horizontal evolution in mammals, with formation of new genes and expansion of gene families, which were displayed in the human genome as tandem gene duplications and pseudogene insertions. In contrast the CER1 middle region contained evolutionarily well-conserved solitary genes and a minimal amount of retroposed genes. The coincidence of evolutionary plasticity with CER1 breakpoints may suggest that regional structural instability is expressed in both evolutionary and cancer-associated chromosome rearrangements.     

Comparative sequence analysis of a single-gene conserved segment in mouse and human

Comparative mapping and sequencing of the mouse and human genomes have defined large, conserved chromosomal segments in which gene content and order are highly conserved. These regions span megabase-sized intervals and together comprise the vast majority of both genomes. However, the evolutionary relationships among the small remaining portions of these genomes are not as well characterized. Here we describe the sequencing and annotation of a 341-kb region of mouse Chr 2 containing nine genes, including biliverdin reductase A (Blvra), and its comparison with the orthologous regions of the human and rat genomes. These analyses reveal that the known conserved synteny between mouse Chromosome (Chr) 2 and human Chr 7 reflects an interval containing one gene (Blvra/BLVRA) that is, at most, just 34 kb in the mouse genome. In the mouse, this segment is flanked proximally by genes orthologous to human chromosome 15q21 and distally by genes orthologous to human Chr 2q11. The observed differences between the human and mouse genomes likely resulted from one or more rearrangements in the rodent lineage. In addition to the resulting changes in gene order and location, these rearrangements also appear to have included genomic deletions that led to the loss of at least one gene in the rodent lineage. Finally, we also have identified a recent mouse-specific segmental duplication. These finding illustrate that small genomic regions outside the large mouse–human conserved segments can contain a single gene as well as sequences that are apparently unique to one genome. The nucleotide sequence data reported in this paper have been submitted to GenBank and assigned the accession numbers AC074224 and AC074041.     

Fine mapping of the muscular dystrophy (AM) gene on chicken chromosome 2q

Our previous studies revealed that the genetic locus for chicken muscular dystrophy of abnormal muscle (AM) mapped to chromosome 2q, and that the region showed conserved synteny with human chromosome 8q11-24.3. In the current study, we mapped the chicken orthologues of genes from human chromosome 8q11-24 in order to identify the responsible gene. Polymorphisms in the chicken orthologues were identified in the parents of the resource family. Twenty-three genes and expressed sequence tags (ESTs) were mapped to chicken chromosome 2 by linkage analysis. The detailed comparative map shows a high conservation of synteny between chicken chromosome 2q and human chromosome 8q. The AM locus was mapped between [inositol(myo)-1(or4)-monophosphatase 1] (IMPA1) gene and [core-binding factor, runt domain, alpha-subunit 2; translocated to 1; cyclin D-related] (CBFA2T1) gene. The genes located between IMPA1 and CBFA2T1 are the most likely candidates for chicken muscular dystrophy.     

Conserved synteny of mammalian imprinted genes in chicken, frog, and fish genomes

Conservation of synteny of mammalian imprinted genes between chicken and human suggested that highly conserved gene clusters were selected long before these genes were recruited for genomic imprinting in mammals. Here we have applied in silico mapping of orthologous genes in pipid frog, zebrafish, spotted green and Japanese pufferfish to show considerable conservation of synteny in lower vertebrates. More than 400 million years ago in a common ancestor of teleost fish and tetrapods, 'preimprinted' chromosome regions homologous to human 6q25, 7q21, 7q32, 11p15, and 15q11-->q12 already contained most present-day mammalian imprinted genes. Interestingly, some imprinted gene orthologues which are isolated from imprinted clusters in mouse and human could be linked to preimprinted regions in lower vertebrates, indicating that separation occurred during mammalian evolution. On the contrary, newly arisen genes by segmental duplication in the mammalian lineage, i.e. SNRPN and FRAT3, were transposed or translocated to imprinted clusters and recruited for parent-specific activity. By analysis of currently available sequences of non-mammalian vertebrates, the imprinted gene clusters homologous to human chromosomes 14q32 and 19q12 are only poorly conserved in chicken, frog, and fish and, therefore, may not have evolved from ancestral preimprinted gene arrays. Evidently, evolution of imprinted gene clusters is an ongoing and dynamic process in mammals. In general, imprinted gene orthologues do not show a higher degree of synteny conservation in vertebrates than non-imprinted genes interspersed with or adjacent to an imprinted cluster.     

The zebrafish fth1, slc3a2, men1, pc, fgf3 and cycd1 genes define two regions of conserved synteny between linkage group 7 and human chromosome 11q13

In addition to being an excellent model system for studying vertebrate development, the zebrafish has become a great tool for gene discovery by mutational analysis. The recent availability of the zebrafish EST database and radiation hybrid mapping panels has dramatically expanded the framework for genomic research in this species. Developing comparative maps of the zebrafish and human genomes is of particular importance for zebrafish mutagenesis studies in which human orthologs are sought for zebrafish genes. However, only partial cDNA sequences are determined routinely for mapped ESTs, leaving the identity of the EST in question. It previously had been reported that zebrafish linkage group 7 shares conserved synteny with human chromosome 11q13. In an effort to further define this relationship, five full-length zebrafish cDNAs, fth1, slc3a2, prkri, cd81, and pc, as well as one putative human gene, DBX were identified and their map positions ascertained. These six genes, along with men1, fgf3 and cycd1 define two regions of conserved synteny between linkage group 7 and 11q13.     

New gene in the homologous human 11q13–q14 and mouse 7F chromosomal regions

Alterations in the chromosomal region 11q13–11q14 are involved in several pathologies in which most of the key genes remain to be identified. In an effort to isolate as many candidates as possible, we are cloning genes from this region. We report here the mapping for a new sequence from 11q13.5–11q14. This sequence, designated D11S833E, putatively encodes a new gene, provisionally named GARP. We cloned its homologous sequence in the mouse and located it on Chromosome (Chr) 7, region F. The human and mouse genes belong to a conserved group of synteny. This, together with the similar conservation of the FGF and TYR genes, indicates that the human 11q13–q14 and mouse 7E-7F regions share homology.     

Chromosomal localization of the chicken and mammalian orthologues of the orphan phosphatase PHOSPHO1 gene

PHOSPHO1 is a recently identified phosphatase expressed at high levels in the chicken growth plate and which may be involved in generating inorganic phosphate for skeletal matrix mineralization. Using a degenerate RT-PCR approach a fragment of human PHOSPHO1 was cloned. This enabled the identification of the human orthologue on HSA17q21, and the mouse orthologue on a region of MMU11 that exhibits conservation of synteny with HSA17q21. Chicken PHOSPHO1 was mapped by SSCP analysis to position 44 cM on GGA27, adjacent to the HOXB@ (44 cM) and COL1A1 (36 cM) loci. Comparison of genes on GGA27 with their orthologues on the preliminary draft of the human genome identifies regions of conserved synteny equivalent to 25 Mb on HSA17q21.2-23.3 and approximately 20 Mb on GGA27 in which the gene order appears to be conserved. Mapping of the PHOSPHO1 genes to regions of HSA17q21.3, MMU11 and GGA27 that exhibit conservation of synteny provides strong evidence that they are orthologous.     

Mapping of FASN and ACACA on two chicken microchromosomes disrupts the human 17q syntenic group well conserved in mammals

Fatty acid synthase and Acetyl-CoA carboxylase are both key enzymes of lipogenesis and may play a crucial role in the weight variability of abdominal adipose tissue in the growing chicken. They are encoded by the FASN and ACACA genes, located on human Chromosome (Chr) 17q25 and on Chr 17q12 or 17q21 respectively, a large region of conserved synteny among mammals. We have localized the homologous chicken genes FASN and ACACA coding for these enzymes, by single-strand conformation polymorphism analysis on different linkage groups of the Compton and East Lansing consensus genetic maps and by FISH on two different chicken microchromosomes. Although synteny is not conserved between these two genes, our results revealed linkage in chicken between FASN and NDPK (nucleoside diphosphate kinase), a homolog to the human NME1 and NME2 genes (non-metastatic cell proteins 1 and 2), both located on human Chr 17q21.3, and also between FASN and H3F3B (H3 histone family 3B), located on human Chr 17q25. The analysis of mapping data from the literature for other chicken and mammalian genes indicates rearrangements have occurred in this region in the mammalian lineage since the mammalian and avian radiation. Received: 8 August 1997 / Accepted: 24 November 1997     

HSA4 and GGA4: remarkable conservation despite 300-Myr divergence

The chicken (GGA) and human (HSA) genomes diverged around 300-350 Myr ago. Due to this large phylogenetic distance, significant synteny conservation has not been anticipated between the genomes of the two species. However, Zoo-FISH with HSA4 chromosome-specific paint on chicken metaphase chromosomes shows that the human chromosome corresponds largely to the GGA4cen-->q26 region. Comparative gene mapping data in the two species, though limited, provide strong support for these observations. The findings, together with the very recently published data on HSA9-GGAZ and HSA12-GGA1, show that some large chromosomal segments share conserved synteny in the two species. These syntenies are considerably disrupted in the mouse. This makes us believe that despite very early divergence, parts of the human and chicken genomes are more conserved than those of human and mouse, which radiated only 100-120 Myr ago. Moreover, the HSA4-GGA4q correspondence points to a "candidate" chromosome from the karyotype of a mammal-bird ancestor. The findings are thus a small but important step toward understanding the evolution of the two genomes.     

The radial arrangement of the human chromosome 7 in the lymphocyte cell nucleus is associated with chromosomal band gene density

In the nuclei of human lymphocytes, chromosome territories are distributed according to the average gene density of each chromosome. However, chromosomes are very heterogeneous in size and base composition, and can contain both very gene-dense and very gene-poor regions. Thus, a precise analysis of chromosome organisation in the nuclei should consider also the distribution of DNA belonging to the chromosomal bands in each chromosome. To improve our understanding of the chromatin organisation, we localised chromosome 7 DNA regions, endowed with different gene densities, in the nuclei of human lymphocytes. Our results showed that this chromosome in cell nuclei is arranged radially with the gene-dense/GC-richest regions exposed towards the nuclear interior and the gene-poorest/GC-poorest ones located at the nuclear periphery. Moreover, we found that chromatin fibres from the 7p22.3 and the 7q22.1 bands are not confined to the territory of the bulk of this chromosome, protruding towards the inner part of the nucleus. Overall, our work demonstrates the radial arrangement of the territory of chromosome 7 in the lymphocyte nucleus and confirms that human genes occupy specific radial positions, presumably to enhance intra- and inter-chromosomal interaction among loci displaying a similar expression pattern, and/or similar replication timing.     

Identification of Cis-regulatory elements in the mouse Pax9/Nkx2-9 genomic region: implication for evolutionary conserved synteny

We previously reported close physical linkage between Pax9 and Nkx2-9 in the human, mouse, and pufferfish (Fugu rubripes) genomes. In this study, we analyzed cis-regulatory elements of the two genes by comparative sequencing in the three species and by transgenesis in the mouse. We identified two regions including conserved noncoding sequences that possessed specific enhancer activities for expression of Pax9 in the medial nasal process and of Nkx2-9 in the ventral neural tube. Remarkably, the latter contained the consensus Gli-binding motif. Interestingly, the identified Pax9 cis-regulatory sequences were located in an intron of the neighboring gene Slc25a21. Close examination of an extended genomic interval around Pax9 revealed the presence of strong synteny conservation in the human, mouse, and Fugu genomes. We propose such an intersecting organization of cis-regulatory sequences in multigenic regions as a possible mechanism that maintains evolutionary conserved synteny.     

Mapping of theARIXHomeodomain Gene to Mouse Chromosome 7 and Human Chromosome 11q13

The recently described homeodomain protein ARIX is expressed specifically in noradrenergic cell types of the sympathetic nervous system, brain, and adrenal medulla. ARIX interacts with regulatory elements of the genes encoding the noradrenergic biosynthetic enzymes tyrosine hydroxylase and dopamine β-hydroxylase, suggesting a role for ARIX in expression of the noradrenergic phenotype. In the study described here, the mouse and human ARIX genes are mapped. Using segregation analysis of two panels of mouse backcross DNA, mouseArixwas positioned approximately 50 cM distal to the centromere of chromosome 7, nearHbb.HumanARIXwas positioned through analysis of somatic cell hybrids and fluorescencein situhybridization of human metaphase chromosomes to chromosome 11q13.3–q13.4. These map locations extend and further define regions of conserved synteny between mouse and human genomes and identify a new candidate gene for inherited developmental disorders linked to human 11q13.     

Murine segmental duplications are hot spots for chromosome and gene evolution

Mouse and rat genomic sequences permit us to obtain a global view of evolutionary rearrangements that have occurred between the two species and to define hallmarks that might underlie these events. We present a comparative study of the sequence assemblies of mouse and rat genomes and report an enrichment of rodent-specific segmental duplications in regions where synteny is not preserved. We show that segmental duplications present higher rates of molecular evolution and that genes in rearranged regions have evolved faster than those located elsewhere. Previous studies have shown that synteny breakpoints between the mouse and the human genomes are enriched in human segmental duplications, suggesting a causative connection between such structures and evolutionary rearrangements. Our work provides further evidence to support the role of segmental duplications in chromosomal rearrangements in the evolution of the architecture of mammalian chromosomes and in the speciation processes that separate the mouse and the rat.     

Genomic arrangement of salinity tolerance QTLs in salmonids: A comparative analysis of Atlantic salmon (Salmo salar) with Arctic charr (Salvelinus alpinus) and rainbow trout (

ABSTRACT: BACKGROUND: Quantitative trait locus (QTL) studies show that variation in salinity tolerance in Arctic charr and rainbow trout has a genetic basis, even though both these species have low to moderate salinity tolerance capacities. QTL were observed to localize to homologous linkage group segments within putative chromosomal regions possessing multiple candidate genes. We compared salinity tolerance QTL in rainbow trout and Arctic charr to those detected in a higher salinity tolerant species, Atlantic salmon. The highly derived karyotype of Atlantic salmon allows for the assessment of whether disparity in salinity tolerance in salmonids is associated with differences in genetic architecture. To facilitate these comparisons, we examined the genomic synteny patterns of key candidate genes in the other model teleost fishes that have experienced three whole-genome duplication (3R) events which preceded a fourth (4R) whole genome duplication event common to all salmonid species. RESULTS: Nine linkage groups contained chromosome-wide significant QTL (AS-2, -4p, -4q, -5, -9, -12p, -12q, -14q -17q, -22, and [MINUS SIGN]23), while a single genome-wide significant QTL was located on AS-4q. Salmonid genomes shared the greatest marker homology with the genome of three-spined stickleback. All linkage group arms in Atlantic salmon were syntenic with at least one stickleback chromosome, while 18 arms had multiple affinities. Arm fusions in Atlantic salmon were often between multiple regions bearing salinity tolerance QTL. Nine linkage groups in Arctic charr and six linkage group arms in rainbow trout currently have no synteny alignments with stickleback chromosomes, while eight rainbow trout linkage group arms were syntenic with multiple stickleback chromosomes. Rearrangements in the stickleback lineage involving fusions of ancestral arm segments could account for the 21 chromosome pairs observed in the stickleback karyotype. CONCLUSIONS: Salinity tolerance in salmonids from three genera is to some extent controlled by the same loci. Synteny between QTL in salmonids and candidate genes in stickleback suggests genetic variation at candidate gene loci could affect salinity tolerance in all three salmonids investigated. Candidate genes often occur in pairs on chromosomes, and synteny patterns indicate these pairs are generally conserved in 2R, 3R, and 4R genomes. Synteny maps also suggest that the Atlantic salmon genome contains three larger syntenic combinations of candidate genes that are not evident in any of the other 2R, 3R, or 4R genomes examined. These larger synteny tracts appear to have resulted from ancestral arm fusions that occurred in the Atlantic salmon ancestor. We hypothesize that the superior hypo-osmoregulatory efficiency that is characteristic of Atlantic salmon may be related to these clusters.