Comparative analysis of sequence characteristics of imprinted genes in human, mouse, and cattle

Genomic imprinting is an epigenetic mechanism that results in monoallelic expression of genes depending on parent-of-origin of the allele. Although the conservation of genomic imprinting among mammalian species has been widely reported for many genes, there is accumulating evidence that some genes escape this conservation. Most known imprinted genes have been identified in the mouse and human, with few imprinted genes reported in cattle. Comparative analysis of genomic imprinting across mammalian species would provide a powerful tool for elucidating the mechanisms regulating the unique expression of imprinted genes. In this study we analyzed the imprinting of 22 genes in human, mouse, and cattle and found that in only 11 was imprinting conserved across the three species. In addition, we analyzed the occurrence of the sequence elements CpG islands, C + G content, tandem repeats, and retrotransposable elements in imprinted and in nonimprinted (control) cattle genes. We found that imprinted genes have a higher G + C content and more CpG islands and tandem repeats. Short interspersed nuclear elements (SINEs) were notably fewer in number in imprinted cattle genes compared to control genes, which is in agreement with previous reports for human and mouse imprinted regions. Long interspersed nuclear elements (LINEs) and long terminal repeats (LTRs) were found to be significantly underrepresented in imprinted genes compared to control genes, contrary to reports on human and mouse. Of considerable significance was the finding of highly conserved tandem repeats in nine of the genes imprinted in all three species. Electronic supplementary material The online version of this article (doi: ) contains supplementary material, which is available to authorized users.     

The distinguishing sequence characteristics of mouse imprinted genes

Sequence comparison analysis has been carried out for 31 imprinted mouse genes and a set of 150 control genes. The imprinted genes were found to be associated with significantly reduced numbers of short interspersed transposable elements (SINEs), in particular SINE Alu repeats. This is similar to recent analyses of human imprinted genes and supports the suggestion that there is either active selection against SINE elements in imprinted regions or a reduced rate of insertion of these elements. The reduction in numbers of SINEs was more consistent in paternally expressed genes, whereas for maternally expressed genes significantly reduced numbers of SINE-B2 elements were coupled with increased numbers of SINE-B4 and SINE-ID elements. Paternally expressed genes were also found to be associated with a lower GC content. Discriminant analysis revealed that the two sub-groups of imprinted genes can be cleanly separated from each other on the basis of their genomic sequence characteristics and that they tend to localize to different genomic compartments. The differences between the sequence characteristics of imprinted and control genes have also enabled us to develop a discriminant function that can be used in a genome-wide screen to identify candidate imprinted genes.     

A novel approach for identifying candidate imprinted genes through sequence analysis of imprinted and control genes

Through the sequence analysis of 27 imprinted human genes and a set of 100 control genes we have developed a novel approach for identifying candidate imprinted genes based on the differences in sequence composition observed. The imprinted genes were found to be associated with significantly reduced numbers of short interspersed transposable element (SINE) Alus and mammalian-wide interspersed repeat (MIR) repeat elements, as previously reported. In addition, a significant association between imprinted genes and increased numbers of low-complexity repeats was also evident. Numbers of the Alu classes AluJ and AluS were found to be significantly depleted in some parts of the flanking regions of imprinted genes. A recent study has proposed that there is active selection against SINE elements in imprinted regions. Alternatively, there may be differences in the rates of insertion of Alu elements. Our study indicates that this difference extends both upstream and downstream of the coding region. This and other consistent differences between the sequence characteristics of imprinted and control genes has enabled us to develop discriminant analysis, which can be used to screen the genome for candidate imprinted genes. We have applied this function to a number of genes whose imprinting status is disputed or uncertain.     

Comparative sequence analysis of equine and human MHC class II DQB genes

The MHC class II DQB gene of horse was isolated and characterized. No obvious mutations causing frame shifts, or destruction of putative protein structure and splicing machinery were detected. Nucleotide sequence of exon 2 was consistent with an allelic sequence of the W23 haplotype. The cytoplasmic region of the equine DQB gene comprised two exons and an intron. A novel fragment of the gene was identified at the 3' intergenic region proximal to the ELA-DQB gene by sequence comparison between the human and horse DQB genes. This sequence showed the highest identity to exon 3 region of the DQB gene, however the 5' half of this exon was truncated as compared with the intact exon. This gene fragment was also identified in the same site of the HLA-DQB gene.     

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.     

Comparative sequence analysis of the mouse and human Lgn1/SMA interval.

Human chromosome 5q11.2-q13.3 and its ortholog on mouse chromosome 13 contain candidate genes for an inherited human neurodegenerative disorder called spinal muscular atrophy (SMA) and for an inherited mouse susceptibility to infection with Legionella pneumophila (Lgn1). These homologous genomic regions also have unusual repetitive organizations that create practical difficulties in mapping and raise interesting issues about the evolutionary origin of the repeats. In an attempt to analyze this region in detail, and as a way to identify additional candidate genes for these diseases, we have determined the sequence of 179 kb of the mouse Lgn1/SMA interval. We have analyzed this sequence using BLAST searches and various exon prediction programs to identify potential genes. Since these methods can generate false-positive exon declarations, our alignments of the mouse sequence with available human orthologous sequence allowed us to discriminate rapidly among this collection of potential coding regions by indicating which regions were well conserved and were more likely to represent actual coding sequence. As a result of our analysis, we accurately mapped two additional genes in the SMA interval that can be tested for involvement in the pathogenesis of SMA. While no new Lgn1 candidates emerged, we have identified new genetic markers that exclude Smn as an Lgn1 candidate. In addition to providing important resources for studying SMA and Lgn1, our data provide further evidence of the value of sequencing the mouse genome as a means to help with the annotation of the human genomic sequence and vice versa.     

Cloning and sequence analysis reveal structural variation among related zein genes in maize

We have isolated a gene encoding one of the 19,000 dalton zein proteins from a maize genomic library constructed in Charon 4A. This gene occurs on a 7.7 kb Eco RI fragment, and based on Southern hybridization analysis, represents one of several homologous sequences present in the maize genome. The nucleotide sequence of the gene predicts a protein composed of 235 amino acids, including a signal peptide of 21 amino acids. There are no intervening sequences in the gene. By comparing the nucleotide sequence of this gene with that of a homologous cDNA clone, we have identified a basis for microheterogeneity within the gene family. The 5′ nucleotide sequences of the genomic and cDNA clones are identical, but they differ in the center of the protein, where repeated amino acid sequences occur. A nucleotide sequence encoding a conserved peptide of 20 amino acids is repeated nine times in the center of both of these clones.     

Cloning and sequence comparison of the mouse, human, and chicken engrailed genes reveal potential functional domains and regulatory regions.

We have isolated and characterized genomic DNA clones for the human and chicken homologues of the mouse En-1 and En-2 genes and determined the genomic structure and predicted protein sequences of both En genes in all three species. Comparison of these vertebrate En sequences with the Xenopus En-2 [Hemmati-Brivanlou et al., 1991) and invertebrate engrailed-like genes showed that the two previously identified highly conserved regions within the En protein ]reviewed in Joyner and Hanks, 1991] can be divided into five distinct subregions, designated EH1 to EH5. Sequences 5' and 3' to the predicted coding regions of the vertebrate En genes were also analyzed in an attempt to identify cis-acting DNA sequences important for the regulation of En gene expression. Considerable sequence similarity was found between the mouse and human homologues both within the putative 5' and 3' untranslated as well as 5' flanking regions. Between the mouse and Xenopus En-2 genes, shorter stretches of sequence similarity were found within the 3' untranslated region. The 5' untranslated regions of the mouse, chicken and Xenopus En-2 genes, however, showed no similarly conserved stretches. In a preliminary analysis of the expression pattern of the human En genes, En-2 protein and RNA were detected in the embryonic and adult cerebellum respectively and not in other tissues tested. These patterns are analogous to those seen in other vertebrates. Taken together these results further strengthen the suggestion that En gene function and regulation has been conserved throughout vertebrate evolution and, along with the five highly conserved regions within the En protein, raise an interesting question about the presence of conserved genetic pathways.     

Comparative analysis of genome maintenance genes in naked mole rat,mouse, and human

Genome maintenance (GM) is an essential defense system against aging and cancer, as both are characterized by increased genome instability. Here, we compared the copy number variation and mutation rate of 518 GM‐associated genes in the naked mole rat (NMR), mouse, and human genomes. GM genes appeared to be strongly conserved, with copy number variation in only four genes. Interestingly, we found NMR to have a higher copy number of CEBPG, a regulator of DNA repair, and TINF2, a protector of telomere integrity. NMR, as well as human, was also found to have a lower rate of germline nucleotide substitution than the mouse. Together, the data suggest that the long‐lived NMR, as well as human, has more robust GM than mouse and identifies new targets for the analysis of the exceptional longevity of the NMR.     

Genomic organization and comparative sequence analysis of the mouse and human FRS2, FRS3 genes

The signaling adapter proteins FRS2 and FRS3 are implicated in the transmission of extracellular signals from nerve growth factor (NGF) or fibroblast growth factor (FGF) receptors to the Ras/mitogen-activated protein kinase signaling cascade. This study presents the genomic sequence and exon-intron organization of the mouse FRS2 and FRS3 loci as well as their evolutionary conservation with their human counterparts. Both FRS2 and FRS3 contain 5 coding exons spanning over 7 kb of genomic sequence with similar exon sizes and organization. Comparative genomic sequence analyses show a highly conserved genomic organization between mouse and human in both FRS2 and FRS3 genes. Non-coding sequences, highly conserved between mouse and human, were identified in the FRS3 introns that may potentially function as regulatory elements. To assay potential differences in their patterns of expression, RT-PCR analysis was used to assay FRS2 and FRS3 expression in the developing embryo and neural tube (NT) during the time of neurogenesis.     

Identification and characterisation of imprinted genes in the mouse.

Imprinted genes are expressed specifically from one or other parental allele. Over 70 are now known, and about one-half of these are expressed from the paternal allele and one-half from the maternal allele. Most imprinted genes are clustered within imprinting regions of the mouse genome, regions which are associated with abnormal phenotypes when inherited uniparentally. Imprinted genes have been identified from surveys based on differential expression or differential methylation according to parental origin, as well as analyses of candidate genes, mutants and imprinted gene clusters. Many imprinted genes affect growth and development, and more than 25 per cent determine non-coding RNAs that may have a function in controlling imprinted gene expression.     

Detecting correlation between sequence and expression divergences in a comparative analysis of human serpin genes

Physiological functions and characteristic structures of the serpin gene superfamily have been studied extensively, yet the evolution of the serpin genes remains unclear. Gene duplication in this superfamily may shed light on this issue. Two models are used to predict the preservation of duplicated genes: the classical model and the duplication–degeneration–complementation (DDC) model. In this study, we analyzed the phylogenetic relationships of 33 human serpin genes and the expression data of some members of the serpin superfamily from a DNA microarray of human leukemia U937 cells with stably inducible expression of the leukemia-related AML1-ETO gene. We then determined the utility of the DDC model by mapping serpin superfamily expression data to the phylogenetic tree. The correlation between sequence and expression divergences as measured by the Pearson correlation coefficient indicated that human serpin genes evolved under the DDC model. Our study provides a new strategy for comparative analysis of gene sequences and microarray data.     

Detecting correlation between sequence and expression divergences in a comparative analysis of human serpin genes

Physiological functions and characteristic structures of the serpin gene superfamily have been studied extensively, yet the evolution of the serpin genes remains unclear. Gene duplication in this superfamily may shed light on this issue. Two models are used to predict the preservation of duplicated genes: the classical model and the duplication-degeneration-complementation (DDC) model. In this study, we analyzed the phylogenetic relationships of 33 human serpin genes and the expression data of some members of the serpin superfamily from a DNA microarray of human leukemia U937 cells with stably inducible expression of the leukemia-related AML1-ETO gene. We then determined the utility of the DDC model by mapping serpin superfamily expression data to the phylogenetic tree. The correlation between sequence and expression divergences as measured by the Pearson correlation coefficient indicated that human serpin genes evolved under the DDC model. Our study provides a new strategy for comparative analysis of gene sequences and microarray data.