Homologous Gene Pairs Research Articles

Two non-homologous brain diseases-related genes, SERPINI1 and PDCD10, are tightly linked by an asymmetric bidirectional promoter in an evolutionarily conserved manner

BackgroundDespite of the fact that mammalian genomes are far more spacious than prokaryotic genomes, recent nucleotide sequencing data have revealed that many mammalian genes are arranged in a head-to-head orientation and separated by a small intergenic sequence. Extensive studies on some of these neighboring genes, in particular homologous gene pairs, have shown that these genes are often co-expressed in a symmetric manner and regulated by a shared promoter region. Here we report the identification of two non-homologous brain disease-related genes, with one coding for a serine protease inhibitor (SERPINI1) and the other for a programmed cell death-related gene (PDCD10), being tightly linked together by an asymmetric bidirectional promoter in an evolutionarily conserved fashion. This asymmetric bidirectional promoter, in cooperation with some cis-acting elements, is responsible for the co-regulation of the gene expression pattern as well as the tissue specificity of SERPINI1 and PDCD10.ResultsWhile SERPINI1 is predominantly expressed in normal brain and down-regulated in brain tumors, PDCD10 is ubiquitously expressed in all normal tissues but its gene transcription becomes aberrant in different types of cancers. By measuring the luciferase activity in various cell lysates, their 851-bp intergenic sequence was shown to be capable of driving the reporter gene expression in either direction. A 175-bp fragment from nt 1 to 175 in the vicinity of PDCD10 was further determined to function as a minimal bidirectional promoter. A critical regulatory fragment, from nt 176-473 outside the minimal promoter in the intergenic region, was identified to contain a strong repressive element for SERPINI1 and an enhancer for PDCD10. These cis-acting elements may exist to help coordinate the expression and regulation of the two flanking genes.ConclusionFor all non-homologous genes that have been described to be closely adjacent in the mammalian genomes, the intergenic region of the head-to-head PDCD10-SERPINI1 gene pair provides an interesting and informative example of a complex regulatory system that governs the expression of both genes not only through an asymmetric bidirectional promoter, but also through fine-tuned regulations with some cis-acting elements.

Methylation of twoHomo sapiens‐specific X‐Y homologous genes in Klinefelter's syndrome (XXY)

An increased incidence of psychiatric and structural brain abnormalities in individuals with Klinefelter syndrome (KS, 47 XXY) could be due to the presence of extra copies of X-Y homologous genes that escape X inactivation. Of particular interest are the two brain-expressed genes Protocadherin11XY (PCDH11XY) and the Synaptobrevin-like gene (SYBL1) which have been duplicated from the X chromosome to the Y chromosome to give X-Y homologous gene pairs that are specific to modern humans. We examined the DNA of KS individuals reported recently by DeLisi et al. 2005 and determined the parental origin of the X alleles, the degree of skewed X inactivation and investigated the CpG island methylation status of PCDH11XY and SYBL1 by bisulphite sequencing and quantification of methylated HpaII sites. We used a novel method for quantification of unmethylated CpGs with the restriction enzyme McrBC which cuts methylated but not unmethylated CpGs. The results showed that KS individuals have two methylated and one unmethylated SYBL1 allele whereas PCDH11XY is unmethylated and escapes X inactivation on the extra X chromosome. Overexpression of PCDH11XY in KS is probable and variable escape from inactivation of this Homo sapiens-specific gene could account for some abnormalities in KS. The origin of the parental alleles or their preferential X inactivation was not associated with psychotic symptoms.

Positional Homology in Bacterial Genomes

In comparative genomic studies, syntenic groups of homologous sequence in the same order have been used as supplementary information that can be used in helping to determine the orthology of the compared sequences. The assumption is that orthologous gene copies are more likely to share the same genome positions and share the same gene neighbors. In this study we have defined positional homologs as those that also have homologous neighboring genes and we investigated the usefulness of this distinction for bacterial comparative genomics. We considered the identification of positionaly homologous gene pairs in bacterial genomes using protein and DNA sequence level alignments and found that the positional homologs had on average relatively lower rates of substitution at the DNA level (synonymous substitutions) than duplicate homologs in different genomic locations, regardless of the level of protein sequence divergence (measured with non-synonymous substitution rate). Since gene order conservation can indicate accuracy of orthology assignments, we also considered the effect of imposing certain alignment quality requirements on the sensitivity and specificity of identification of protein pairs by BLAST and FASTA when neighboring information is not available and in comparisons where gene order is not conserved. We found that the addition of a stringency filter based on the second best hits was an efficient way to remove dubious ortholog identifications in BLAST and FASTA analyses. Gene order conservation and DNA sequence homology are useful to consider in comparative genomic studies as they may indicate different orthology assignments than protein sequence homology alone.

Sex chromosome evolution: Molecular aspects of Y-chromosome degeneration inDrosophila

Ancient Y-chromosomes of various organisms contain few active genes and an abundance of repetitive DNA. The neo-Y chromosome of Drosophila miranda is in transition from an ordinary autosome to a genetically inert Y-chromosome, while its homolog, the neo-X chromosome, is evolving partial dosage compensation. Here, I compare four large genomic regions located on the neo-sex chromosomes that contain a total of 12 homologous genes. In addition, I investigate the partial coding sequence for 56 more homologous gene pairs from the neo-sex chromosomes. Little modification has occurred on the neo-X chromosome, and genes are highly constrained at the protein level. In contrast, a diverse array of molecular changes is contributing to the observed degeneration of the neo-Y chromosome. In particular, the four large regions surveyed on the neo-Y chromosome harbor several transposable element insertions, large deletions, and a large structural rearrangement. About one-third of all neo-Y-linked genes are nonfunctional, containing either premature stop codons and/or frameshift mutations. Intact genes on the neo-Y are accumulating amino acid and unpreferred codon changes. In addition, both 5'- and 3'-flanking regions of genes and intron sequences are less constrained on the neo-Y relative to the neo-X. Despite heterogeneity in levels of dosage compensation along the neo-X chromosome of D. miranda, the neo-Y chromosome shows surprisingly uniform signs of degeneration.

Compositional transitions between Oryza sativa and Arabidopsis thaliana genes are linked to the functional change of encoded proteins

Compositional transition was previously reported between homologous pairs of Oryza sativa and Arabidopsis thaliana. We analysed 2041 homologous gene pairs of O. sativa and A. thaliana and found that the genes that underwent compositional transition have significantly high frequency of amino acids encoded by GC rich codons. Amino acid substitution matrix calculated among hydrophobic, hydrophilic and amphipathic amino acid groups’ shows that amphipathic and hydrophilic amino acids are strongly substituted by hydrophobic amino acids, whereas the substitution from hydrophobic to hydrophilic and amphipathic acids are not so frequent. Moreover, O. sativa genes that experienced compositional transition are more hydrophobic than their A. thaliana counterpart. This indicates that nucleotide bias induces a directional changes in proteome composition in such a way that underwent strong changes in hydropathy values. Further analysis of amino acid composition and hydrophobicity on the three different protein secondary structures confirm that natural selection is the main driving force for compositional transition in O. sativa genome.

Indirect evidence from DNA sequence diversity for genetic degeneration of the Y‐chromosome in dioecious species of the plant Silene: the SlY4/SlX4 and DD44‐X/DD44‐Y gene pairs

The action of natural selection is expected to reduce the effective population size of a nonrecombining chromosome, and this is thought to be the chief factor leading to genetic degeneration of Y-chromosomes, which cease recombining during their evolution from ordinary chromosomes. Low effective population size of Y chromosomes can be tested by studying DNA sequence diversity of Y-linked genes. In the dioecious plant, Silene latifolia, which has sex chromosomes, one comparison (SlX1 vs. SlY1) indeed finds lower Y diversity compared with the homologous X-linked gene, and one Y-linked gene with no X-linked homologue has lower species-wide diversity than a homologous autosomal copy (SlAp3Y vs. SlAp3A). To test whether this is a general pattern for Y-linked genes, we studied two further recently described X and Y homologous gene pairs in samples from several populations of S. latifolia and S. dioica. Diversity is reduced for both Y-linked genes, compared with their X-linked homologues. Our new data are analysed to show that the low Y effective size cannot be explained by different levels of gene flow for the X vs. the Y chromosomes, either between populations or between these closely related species. Thus, all four Y-linked genes that have now been studied in these plants (the two studied here, and two previously studied genes, have low diversity). This supports other evidence for an ongoing degeneration process in these species.

Sex differences in structure and expression of the sex chromosome genes CHD1Z and CHD1W in zebra finches.

Genes on the sex chromosomes are unique because of their sex-specific inheritance. One question is whether homologous gene pairs on the sex chromosomes, which have diverged in their sequence, have acquired different functions. We have analyzed the first homologous pair of genes (CHD1Z and CHD1W) discovered on the avian Z and W sex chromosomes of the zebra finch (Taeniopygia guttata) to examine whether functional differences may have evolved. Sequence analysis revealed that the two genes maintained a high degree of similarity especially within the C, H, and D domains, but outside of these regions larger differences were observed. Expression studies showed that CHD1W was unique to females and has the potential to produce a protein that CHD1Z does not. CHD1Z mRNA was expressed at a higher level in the male brain than in the female brain at various post-hatch ages. Reporter constructs containing the 5' flanking regions of each gene showed they had the ability to drive reporter expression in primary cell cultures. The 5' flanking region sequence of CHD1Z and CHD1W exhibited little homology, and differences in putative promoter elements were apparent. These differences between CHD1Z and CHD1W suggest that the two proteins may have diverged in their function.

Expression pattern of yeast sporulation: a case study for regulatory changes after yeast genome duplication

While nucleotide substitutions in both regulatory and coding regions are the two facets of functional divergence after gene duplication, evidence is limited for our understanding the extent to which regulatory changes can occur after gene (genome) duplication. Fortunately, the complete yeast genome sequence and the whole-genome micro-array expression data provide a unique opportunity for resolving this issue. In the yeast genome, a genome-wide duplication may result in 54 homologous chromosome block pairs, including 452 gene pairs. On the other hand, Chu et al. [Science 282 (1998) 699] used DNA micro-arrays containing nearly every yeast gene to assay changes in gene expression during sporulation, which consists of meiosis and spore morphogenesis. Our results have shown that, for a subset of gene pairs that at least one duplicate gene is induced during sporulation, there is a significant association between expression and the presence of DNA binding site MSE, while no such association is observed for the rest of homologous gene pairs.

Sex differences in sex chromosome gene expression in mouse brain.

A major question is whether genes encoded on the sex chromosomes act directly in non-gonadal tissues to cause sex differences in development or function, or whether all sex differences in somatic tissues are induced by gonadal secretions. As part of this question we asked whether mouse X-Y homologous gene pairs are expressed in brain in a sex-specific fashion. Using RT-PCR and northern blot analysis, we assessed mRNA expression in brain of eight Y-linked genes as well as their X-linked homologues, at three ages: 13.5 days post coitum, the day of birth (P1) and adult. Transcripts of six Y genes were expressed at one or more ages: Usp9y, Ube1y, Smcy, Eif2s3y, Uty and Dby. Their expression also occurred in XY female brain, and therefore does not require testicular secretions. Six X-linked homologues (Usp9x, Ube1x, Smcx, Eif2s3x, Utx and Dbx) were also expressed in brain, and in adulthood all of these transcripts were expressed at significantly higher levels in brains of females than in brains of males, irrespective of their X-inactivation status. For five of these gene pairs, the expression of the Y-linked homologue in males was not sufficient to compensate for the female bias in X gene expression. Three X-Y gene pairs, Usp9x/y, Ube1x/y and Eif2s3x/y, appeared to be differentially regulated (expressed in brain in a different age- or tissue-dependent pattern), and hence may not be functionally equivalent. These sex differences in X-Y gene expression suggest several mechanisms by which these genes may participate in sex differences in brain development and function.

An exhaustive search for tandem repeats and long duplicated chromosomal regions within the mouse genome.

Traces of ancient duplications of extensive chromosomal regions are being discovered within the human genome. For example, the MHC (major histocompatibility complex) gene region on chromosome 6 (6p21.3) encodes a cluster of genes that are homologous to those on chromosome 9 (9q33-q34) . These regions appear to have arisen from regional duplication around the time of vertebrate emergence. To elucidate the evolution of mammalian genomes, it is crucial to search for duplicated chromosomal regions in the mouse genome. In this study, we searched for long duplicated chromosomal regions in the mouse genome using the map information derived from the Mouse Genome Database and the numerous homologous gene pairs from GenBank. To identify the regions, we searched for homologous gene pairs among the genomic sequences extracted from GenBank (r.118) . We defined the candidates of duplicated regions as those having more than two homologous gene pairs located within 5cM on each chromosome. We conducted a statistical test considering the distribution of homologous gene pairs and tandem repeats in the mouse genome. Twenty-seven pairs of duplicated regions were found in this study. Seven procollagen genes, three Hox gene clusters, six integrin genes, three fibroblast growth factor receptor genes and three eye absent homolog genes were located in the duplicated regions.

A human sex-chromosomal gene family expressed in male germ cells and encoding variably charged proteins.

Approximately 12 X-Y homologous gene pairs have been identified in the non-recombining portions of human sex chromosomes. These X-Y gene pairs fall into two categories. In the first category, both X and Y homologs are ubiquitously expressed. In the second category, the X homolog is ubiquitously expressed, whereas the Y homolog is expressed exclusively in the testis. Here we describe a family of human X-Y genes that cannot be assigned to either category. Designated VCX / Y ( Variable Charge X / Y; VCY previously known as BPY1 ), this gene family has multiple members on both X and Y, and all appear to be expressed exclusively in male germ cells. Members of the VCX / Y family share a high degree of sequence identity, with the exception that a 30 nucleotide unit is tandemly repeated in X-linked members but is present only once in Y-linked members. These atypical features suggest that the VCX / Y family has evolved in a manner previously unrecognized for mammalian X-Y genes. We also found that a copy of VCX is present in CRI-S232, a previously described genomic fragment derived from the X chromosome. Studies have shown that aberrant recombination between arrays of CRI-S232-homologous repeats flanking the steroid sulfatase ( STS ) gene results in STS deletion, which is manifested clinically as X-linked ichthyosis. The revelation that CRI-S232 contains VCX offers a more precise description of the genetic etiology of X-linked ichthyosis: it results from aberrant recombination between VCX gene arrays that flank the STS locus.

Genealogy and subpopulation differentiation under various models of population structure

The structured coalescent is used to calculate some quantities relating to the genealogy of a pair of homologous genes and to the degree of subpopulation differentiation, under a range of models of subdivided populations and assuming the infinite alleles model of neutral mutation. The classical island and stepping-stone models of population structure are considered, as well as two less symmetric models. For each model, we calculate the Laplace transform of the distribution of the coalescence time of a pair of genes from specified locations and the corresponding mean and variance. These results are then used to calculate the values of Wright’s coefficient F ST , its limit as the mutation rate tends to zero and the limit of its derivative with respect to the mutation rate as the mutation rate tends to zero. From this derivative it is seen that F ST can depend strongly on the mutation rate, for example in the case of an essentially one-dimensional habitat with many subpopulations where gene flow is restricted to neighbouring subpopulations.

Silent substitutions in mammalian genomes and their evolutionary implications.

An analysis of silent substitutions in pairwise comparisons of homologous genes from different mammals has shown that, in spite of individual fluctuations, their frequencies (which are very strongly correlated with the frequency of substitutions per synonymous site calculated according to Li et al. 1985) do not vary, on the average, with the GC levels of silent positions. This holds in the general case, in which silent positions of pairs of homologous genes share the same composition, namely in the human/other primates, human/artiodactyls, and in the mouse/rat pairs, as well as in the special cases in which the composition of silent positions are different, namely in the human/rabbit and the human/rat (or human/mouse) pairs. A slightly lower frequency found for low GC values in the human/bovine and human/pig pairs seems to be due to the specific gene samples used. These results contradict the previously claimed existence of differences in mutation rates and of mutational biases in third codon positions of coding sequences located in different isochores of mammalian genomes. They also imply that the variations in nucleotide precursor pools through the cell cycle and the differences in replication timing, or in repair efficiency, which were reported for different isochores, do not lead, as claimed, to differences in mutation rates, not in mutational biases in mammals. The differences claimed appear to be due to using small gene samples when individual fluctuations from gene to gene are relatively large.

Two novel transmembrane protein tyrosine kinases expressed during Caenorhabditis elegans hypodermal development.

We describe our characterization of kin-15 and kin-16, a tandem pair of homologous Caenorhabditis elegans genes encoding transmembrane protein tyrosine kinases (PTKs) with an unusual structure: the predicted extracellular domain of each putative gene product is only about 50 amino acids, and there are no potential autophosphorylation sites in the C-terminal domain. Using lacZ fusions, we found that kin-15 and kin-16 both appear to be expressed during postembryonic development in the large hypodermal syncytium (hyp7) around the time that specific hypodermal cells fuse with hyp7. kin-15 and kin-16 were positioned on the genetic and physical maps, but extrachromosomal arrays containing wild-type kin-15 and/or kin-16 genes were unable to complement candidate lethal mutations. The results suggest that kin-15 and kin-16 may be specifically involved in cell-cell interactions regulating cell fusions that generate the hypodermis during postembryonic development.

Two Novel Transmembrane Protein Tyrosine Kinases Expressed during Caenorhabditis elegans Hypodermal Development

We describe our characterization of kin-15 and kin-16, a tandem pair of homologous Caenorhabditis elegans genes encoding transmembrane protein tyrosine kinases (PTKs) with an unusual structure: the predicted extracellular domain of each putative gene product is only about 50 amino acids, and there are no potential autophosphorylation sites in the C-terminal domain. Using lacZ fusions, we found that kin-15 and kin-16 both appear to be expressed during postembryonic development in the large hypodermal syncytium (hyp7) around the time that specific hypodermal cells fuse with hyp7. kin-15 and kin-16 were positioned on the genetic and physical maps, but extrachromosomal arrays containing wild-type kin-15 and/or kin-16 genes were unable to complement candidate lethal mutations. The results suggest that kin-15 and kin-16 may be specifically involved in cell-cell interactions regulating cell fusions that generate the hypodermis during postembryonic development.

Genetic linkage analysis of the murine developmental mutant velvet coat (Ve) and the distal chromosome 15 developmental genes Hox-3.1, Rar-g, Wnt-1, and Krt-2.

We have identified restriction fragment length polymorphisms between Mus musculus and Mus spretus for the Chromosome 15 loci Hox-3, Wnt-1, Krt-2, Rar-g, and Ly-6. We followed the inheritance of these alleles in interspecific genetic test crosses between velvet coat (Ve) heterozygotes and M. spretus. The results suggest a gene order and recombination distances (in cM) of Ly-6-22-Wnt-1-2-Ve/Krt-2/Rar-g-3-Hox-3. No recombination was found between Ve, Krt-2, and Rar-g. The data also provide evidence for the hypothesis of a large-scale genomic duplication involving homologous gene pairs on mouse Chromosomes 15 and 11.

Codon usage changes and sequence dissimilarity between human and rat

This paper reports on the relationship between the number of silent differences and the codon usage changes in the lineages leading to human and rat. Examination of 102 pairs of homologous genes gives rise to four main conclusions: (1) We have previously demonstrated the existence of a codon usage change (called the minor shift) between human and rat; this was confirmed here with a larger sample. For genes with extreme C & G frequencies, the C & G level in the third codon position is less extreme in rat than in human. (2) Protein similarity and percentage of positive differences are the two main factors that discriminate homologous genes when characterized by differences between rat and human. By definition, positive differences result from silent changes between A or T and C or G with a direction implying a C & G content variation in the same direction as the overall gene variation. (3) For genes showing both codon usage change and low protein similarity, a majority of amino acid replacements contributes to C & G level variation in positions I and II in the same direction as the variation in position III. This is thus a new example of protein evolution due to constraints acting at the DNA level. (4) In heavy isochores (high C & G content) no direct correlation exists between codon usage change (measured by the dissymmetry of differences) and silent dissimilarity. In light isochores the opposite situation is observed: modification of codon usage is associated with a high synonymous dissimilarity. This result shows that, in some cases, modification of constrains acting at the DNA level could accelerate divergence between genomes.

DNA microenvironments and the molecular clock

A few years ago we presented a stationary Markov model of gene evolution according to which only homologous genes from not too divergent species obeying the condition of being stationary may behave as reliable molecular clocks. A compartmentalized model of the nuclear genome in which the genes are distributed in compartments, the isochores, defined by their G + C content has been proposed recently. We have found that only homologous gene pairs that are stationary, and belong to the same isochore, can be used consistently for the determination of phylogeny and base substitution rate. In particular, for the rodent-human couple, only about half of the homologous gene pairs are stationary. Stationary genes evolve at the third silent codon position with the same velocity independent of the genes and base composition. By contrast, nonstationary genes display apparent rate values (pseudovelocities) that are significantly higher. Our results cast doubt upon recent claims of a large acceleration in the rate of molecular evolution in rodents.

Homologous expressed genes in the human sex chromosome pairing region.

The human sex chromosomes share a pair of homologous genes which independently encode a cell-surface antigen defined by the monoclonal antibody 12E7 (refs 1, 2; see refs 3, 4 for review). The X-located homologue, MIC2X, escapes X-inactivation and the equivalent Y-located locus, MIC2Y, was one of the first genes shown to reside on a mammalian Y chromosome. By using a bacterial expression system we have previously cloned a complementary DNA sequence corresponding to a MIC2 gene and have used this probe to show that the MIC2X and MIC2Y loci are closely related, if not identical. Here we report the use of the cloned probe to confirm the localization of the MIC2X locus to the region Xpter-Xp22.32 (ref. 7) and demonstrate, for the first time, that the MIC2Y locus is located on the short arm of the Y chromosome in the distal region Ypter-Yp11.2. The MIC2 sequences and the sequences described in the accompanying papers by Cooke et al. and Simmler et al. are the first which have been shown to be shared by the sex chromosomes in the pairing region.

Selection against dam methylation sites in the genomes of DNA of enterobacteriophages.

Postreplicative methylation of adenine in Escherichia coli DNA to produce G6m ATC (where 6mA is 6-methyladenine) has been associated with preferential daughter-strand repair and possibly regulation of replication. An analysis was undertaken to determine if these, or other, as yet unknown roles of GATC, have had an effect on the frequency of GATC in E. coli or bacteriophage DNA. It was first ascertained that the most accurate predictions of GATC frequency were based on the observed frequencies of GAT and ATC, which would be expected since these predictors take into account preferences in codon usage. The predicted frequencies were compared with observed GATC frequencies in all available bacterial and phage nucleotide sequences. The frequency of GATC was close to the predicted frequency in most genes of E. coli and its RNA bacteriophages and in the genes of nonenteric bacteria and their bacteriophages. However, for DNA enterobacteriophages the observed frequency of GATC was generally significantly lower than predicted when assessed by the chi square test. No elevation in the rate of mutation of 6mA in GATC relative to other bases was found when pairs of DNA sequences from closely related phages or pairs of homologous genes from enterobacteria were compared, nor was any preferred pathway for mutation of 6mA evident in the E. coli DNA bacteriophages. This situation contrasts with that of 5-methylcytosine, which is hypermutable, with a preferred pathway to thymine. Thus, the low level of GATC in enterobacteriophages is probably due not to 6mA hypermutability, but no selection against GATC in order to bypass a GATC-mediated host function.(ABSTRACT TRUNCATED AT 250 WORDS)