Tracing the origin of functional and conserved domains in the human proteome: implications for protein evolution at the modular level

Background

The highly homologous PE_PGRS (Proline-glutamic acid_polymorphic GC-rich repetitive sequence) genes are members of the PE multigene family which is found only in mycobacteria. PE genes are particularly abundant within the genomes of pathogenic mycobacteria where they seem to have expanded as a result of gene duplication events. PE_PGRS genes are characterized by their high GC content and extensive repetitive sequences, making them prone to recombination events and genetic variability.

Results

Comparative sequence analysis of Mycobacterium tuberculosis genes PE_PGRS17 (Rv0978c) and PE_PGRS18 (Rv0980c) revealed a striking genetic variation associated with this typical tandem duplicate. In comparison to the M. tuberculosis reference strain H37Rv, the variation (named the 12/40 polymorphism) consists of an in-frame 12-bp insertion invariably accompanied by a set of 40 single nucleotide polymorphisms (SNPs) that occurs either in PE_PGRS17 or in both genes. Sequence analysis of the paralogous genes in a representative set of worldwide distributed tubercle bacilli isolates revealed data which supported previously proposed evolutionary scenarios for the M. tuberculosis complex (MTBC) and confirmed the very ancient origin of " M. canettii " and other smooth tubercle bacilli. Strikingly, the identified polymorphism appears to be coincident with the emergence of the post-bottleneck successful clone from which the MTBC expanded. Furthermore, the findings provide direct and clear evidence for the natural occurrence of gene conversion in mycobacteria, which appears to be restricted to modern M. tuberculosis strains.

Conclusion

This study provides a new perspective to explore the molecular events that accompanied the evolution, clonal expansion, and recent diversification of tubercle bacilli.     

First-generation SNP/InDel markers tagging loci for pathogen resistance in the potato genome

A panel of 17 tetraploid and 11 diploid potato genotypes was screened by comparative sequence analysis of polymerase chain reaction (PCR) products for single nucleotide polymorphisms (SNPs) and insertion-deletion polymorphisms (InDels), in regions of the potato genome where genes for qualitative and/or quantitative resistance to different pathogens have been localized. Most SNP and InDel markers were derived from bacterial artificial chromosome (BAC) insertions that contain sequences similar to the family of plant genes for pathogen resistance having nucleotide-binding-site and leucine-rich-repeat domains (NBS-LRR-type genes). Forty-four such NBS-LRR-type genes containing BAC-insertions were mapped to 14 loci, which tag most known resistance quantitative trait loci (QTL) in potato. Resistance QTL not linked to known resistance-gene-like (RGL) sequences were tagged with other markers. In total, 78 genomic DNA fragments with an overall length of 31 kb were comparatively sequenced in the panel of 28 genotypes. 1498 SNPs and 127 InDels were identified, which corresponded, on average, to one SNP every 21 base pairs and one InDel every 243 base pairs. The nucleotide diversity of the tetraploid genotypes (pi = 0.72 x 10(-3)) was lower when compared with diploid genotypes (pi = 2.31 x 10(-3)). RGL sequences showed higher nucleotide diversity when compared with other sequences, suggesting evolution by divergent selection. Information on sequences, sequence similarities, SNPs and InDels is provided in a database that can be queried via the Internet.     

Characterization of the cow stomach lysozyme genes: Repetitive DNA and concerted evolution

Cow stomach lysozyme genes have evolved in a mosaic pattern. The majority of the intronic and flanking sequences show an amount of sequence difference consistent with divergent evolution since duplication of the genes 40–50 million years ago. In contrast, exons 1, 2, and 4 and immediately adjacent intronic sequences differ little between genes and show evidence of recent concerted evolution. Exon 3 appears to be evolving divergently. The three characterized genes vary from 5.6 to 7.9 kilobases in length. Different distributions of repetitive DNA are found in each gene, which accounts for the majority of length differences between genes. The different distributions of repetitive DNA in each gene suggest the repetitive elements were inserted into each gene after the duplications that give rise to these three genes and provide additional support for divergent evolution for the majority of each gene. The observation that intronic and flanking sequences are evolving divergently suggests that the concerted evolution events involved in homogenizing the coding regions of lysozyme genes involve only one exon at a time. This model of concerted evolution would allow the shuffling of exon-sized pieces of information between genes, a phenomenon that may have aided in the early adaptive evolution of stomach lysozyme.Deceased July 21, 1991 Correspondence to: D.M. Irwin     

Repeated sequences in CASPASE-5 and FANCD2 but not NF1 are targets for mutation in microsatellite-unstable acute leukemia/myelodysplastic syndrome

Microsatellite instability (MSI) in tumors is diagnostic for inactive DNA mismatch repair. It is widespread among some tumor types, such as colorectal or endometrial carcinoma, but is rarely found in leukemia. Therapy-related acute myeloid leukemia/myelodysplastic syndrome (tAML/MDS) is an exception, and MSI is frequent in tAML/MDS following cancer chemotherapy or organ transplantation. The development of MSI+ tumors is associated with an accumulation of insertion/deletion mutations in repetitive sequences. These events can cause inactivating frameshifts or loss of expression of key growth control proteins. We examined established MSI+ cell lines and tAML/MDS cases for frameshift-like mutations of repetitive sequences in several genes that have known, or suspected, relevance to leukemia. CASPASE-5, an acknowledged frameshift target in MSI+ gastrointestinal tract tumors, was frequently mutated in MSI+ cell lines (67%) and in tAML/MDS (29%). Frameshift-like mutations were also observed in the NF1 and FANCD2 genes that are associated with genetic conditions conferring a predisposition to leukemia. Both genes were frequent targets for mutation in MSI+ cell lines and colorectal carcinomas. FANCD2 mutations were also common in MSI+ tAML/MDS, although NF1 mutations were not observed. A novel FANCD2 polymorphism was also identified.     

Frequent homologous recombination events in Mycobacterium tuberculosis PE/PPE multigene families: potential role in antigenic variability

The PE and PPE (PE/PPE) multigene families of Mycobacterium tuberculosis are particularly GC-rich and share extensive homologous repetitive sequences. We hypothesized that they may undergo homologous recombination events, a mechanism rarely described in the natural evolution of mycobacteria. To test our hypothesis, we developed a specific oligonucleotide-based microarray targeting nearly all of the PE/PPE genes, aimed at detecting signals for homologous recombination. Such a microarray has never before been reported due to the multiplicity and highly repetitive and homologous nature of these sequences. Application of the microarray to a collection of M. tuberculosis clinical isolates (n = 33) representing prevalent spoligotype strain families in Tunisia allowed successful detection of six deleted genomic regions involving a total of two PE and seven PPE genes. Some of these deleted genes are known to be immunodominant or involved in virulence. The four precisely determined deletions were flanked by 400- to 500-bp stretches of nearly identical sequences lying mainly at the conserved N-terminal region of the PE/PPE genes. These highly homologous sequences thus serve as substrates to mediate both intergenic and intragenic homologous recombination events, indicating an important function in generating strain variation. Importantly, all recombination events yielded a new in-frame fusion chimeric gene. Hence, homologous recombination within and between PE/PPE genes likely increased their antigenic variability, which may have profound implications in pathogenicity and/or host adaptation. The finding of high prevalence (approximately 45% and approximately 58%) for at least two of the genomic deletions suggests that they likely confer advantageous biological attributes.     

Rapid genome divergence at orthologous low molecular weight glutenin loci of the A and Am genomes of wheat

To study genome evolution in wheat, we have sequenced and compared two large physical contigs of 285 and 142 kb covering orthologous low molecular weight (LMW) glutenin loci on chromosome 1AS of a diploid wheat species (Triticum monococcum subsp monococcum) and a tetraploid wheat species (Triticum turgidum subsp durum). Sequence conservation between the two species was restricted to small regions containing the orthologous LMW glutenin genes, whereas >90% of the compared sequences were not conserved. Dramatic sequence rearrangements occurred in the regions rich in repetitive elements. Dating of long terminal repeat retrotransposon insertions revealed different insertion events occurring during the last 5.5 million years in both species. These insertions are partially responsible for the lack of homology between the intergenic regions. In addition, the gene space was conserved only partially, because different predicted genes were identified on both contigs. Duplications and deletions of large fragments that might be attributable to illegitimate recombination also have contributed to the differentiation of this region in both species. The striking differences in the intergenic landscape between the A and A(m) genomes that diverged 1 to 3 million years ago provide evidence for a dynamic and rapid genome evolution in wheat species.     

The characterization and comparative analysis of high-molecular-weight glutenin genes from genomes A and B of a hexaploid bread wheat

Summary Two high-molecular-weight subunit (HMWS) glutenin genes from the A and B genomes of the hexaploid bread wheat Triticum aestivum L. cv Cheyenne have been isolated and sequenced. Both of these genes are of the high M_r class (x-type) of HMW glutenins, and have not been previously reported. The entire set of six HMW genes from cultivar Cheyenne have now been isolated and characterized. An analysis of the Ax and Bx sequences shows that the Ax sequence is similar to the homoeologous gene from the D genome, while the Bx repeat structure is significantly different. The repetitive region of these proteins can be modelled as a series of interspersed copies of repeat modifs of 6, 9, and 15 amino acid residues. The evolution of these genes includes single-base substitutions over the entire coding region, plus insertion/deletions of single or blocks of repeats in the central repetitive domain.     

Evolutionary Process of the Genomic Sequence Around the 100 Map Unit of Chromosome 1 in <Emphasis Type="Italic">Arabidopsis thaliana</Emphasis>

Comparative analysis of nucleotide sequences of the genomic region located around 100 map unit of chromosome 1 using two accessions, Columbia (Col) and Landsberg erecta (Ler), of Arabidopsis thaliana was performed. High divergence was detected between them, and the length of the Ler sequence was half of corresponding sequence of Col. This divergence occurred by tandem duplication, deletion of large regions, and insertion of unrelated sequences. These events led to the high polymorphism of plant disease resistant genes, which are located in the analyzed region. It is highly probable that two-round duplication occurred, and the insertion sequences are transposable elements. The data suggest that the analyzed region had been evolving until quite recently.     

Microbial genome evolution: sources of variability

Comparative genome analyses of close relatives have yielded exciting insight into the sources of microbial genome variability with respect to gene content, gene order and evolution of genes with unknown functions. The genomes of free-living bacteria often carry phages and repetitive sequences that mediate genomic rearrangements in contrast to the small genomes of obligate host-associated bacteria. This suggests that genomic stability correlates with the genomic content of repeated sequences and movable genetic elements, and thereby with bacterial lifestyle. Genes with unknown functions present in a single species tend to be shorter than conserved, functional genes, indicating that the fraction of unique genes in microbial genomes has been overestimated.     

A unique structure of a mouse gamma-actin processed pseudogene

We have isolated several gamma-actin-related genes from a mouse genomic library. One of these has been shown to be a gamma-actin processed pseudogene (Tokunga, K., Yoda, K. and Sakiyama, S. (1985) Nucleic Acids Res. 13, 3031-3042). Here, we report the structure of another pseudogene (pMA131). pMA131 contained the sequences corresponding to the carboxyl half of a cytoskeletal actin in which random point mutations as well as insertion and deletion events took place. This region was flanked at its 5' end by the sequences related to mouse repetitive sequences, including the MIF-1 family, and was interrupted by the sequence homologous to the R family which is also a mouse repetitive sequence. The coding region was followed by the sequence corresponding to 3' untranslated region of gamma-actin mRNA.     

Exon junction sequences as cryptic splice sites: implications for intron origin

Introns are flanked by a partially conserved coding sequence that forms the immediate exon junction sequence following intron removal from pre-mRNA. Phylogenetic evidence indicates that these sequences have been targeted by numerous intron insertions during evolution, but little is known about this process. Here, we test the prediction that exon junction sequences were functional splice sites that existed in the coding sequence of genes prior to the insertion of introns. To do this, we experimentally identified nine cryptic splice sites within the coding sequence of actin genes from humans, Arabidopsis, and Physarum by inactivating their normal intron splice sites. We found that seven of these cryptic splice sites correspond exactly to the positions of exon junctions in actin genes from other species. Because actin genes are highly conserved, we could conclude that at least seven actin introns are flanked by cryptic splice sites, and from the phylogenetic evidence, we could also conclude that actin introns were inserted into these cryptic splice sites during evolution. Furthermore, our results indicate that these insertion events were dependent upon the splicing machinery. Because most introns are flanked by similar sequences, our results are likely to be of general relevance.     

Partial reversion at the bobbed locus of Drosophila melanogaster

In Drosophila melanogaster the tandemly arranged repetitive sequences coding for 18S and 28S rRNA are heterogenous at the level of the spacers between units and insertions that interrupt many 28S rRNA genes. This heterogeneity contrasts with the homogeneity of the regions transcribed into 18S and 28S rRNA. Homogenization and evolution of repetitive genes are usually explained by conversion, amplification events or unequal crossovers. In this paper we studied the change in rDNA patterns associated with partial reversion of bobbed mutations. In most cases, no increase in rDNA gene number, but a new repartition of gene types were found.     

A ribosomal orphon sequence from Xenopus laevis flanked by novel low copy number repetitive elements

We have used a differential cloning approach to isolate ribosomal/non-ribosomal frontier sequences from Xenopus laevis. A ribosomal intergenic spacer sequence (IGS) was cloned and shown not to be physically linked with the ribosomal locus. This ribosomal orphon contained the IGS sequences found immediately downstream of the 28S gene and included an array of enhancer repetitions and a non-functional spacer promoter. The orphon sequence was flanked by a member of the novel 'Frt' low copy repetitive element family. Three individual Frt repeats were sequenced and all members of this family were shown to lie clustered at two chromosomal sites, one of which contained the ribosomal orphon. One of the Frt elements contained an insertion of 297 bp that showed extensive homology to sequences within at least three other Xenopus genes. Each homology region was flanked by members of the T2 family of short interspersed repetitive elements, (SINEs), and by its target insertion sequence, suggesting multiple translocation events. The data are discussed in terms of the evolution of the ribosomal gene locus.     

Comparative analysis of the Hrp pathogenicity island of Rubus- and Spiraeoideae-infecting Erwinia amylovora strains identifies the IT region as a remnant of an integrative conjugative element

The Hrp pathogenicity island (hrpPAI) of Erwinia amylovora not only encodes a type III secretion system (T3SS) and other genes required for pathogenesis on host plants, but also includes the so-called island transfer (IT) region, a region that originates from an integrative conjugative element (ICE). Comparative genomic analysis of the IT regions of two Spiraeoideae- and three Rubus-infecting strains revealed that the regions in Spiraeoideae-infecting strains were syntenic and highly conserved in length and genetic information, but that the IT regions of the Rubus-infecting strains varied in gene content and length, showing a mosaic structure. None of the ICEs in E. amylovora strains were complete, as conserved ICE genes and the left border were missing, probably due to reductive genome evolution. Comparison of the hrpPAI region of E. amylovora strains to syntenic regions from other Erwinia spp. indicates that the hrpPAI and the IT regions are the result of several insertion and deletion events that have occurred within the ICE. It also suggests that the T3SS was present in a common ancestor of the pathoadapted Erwinia spp. and that insertion and deletion events in the IT region occurred during speciation.     

Circularization of transduced fragments: a mechanism for adding segments to the bacterial chromosome

Generalized transducing fragments that have redundant sequences in direct order can circularize during transduction events. The length of the required redundant sequences can be at least as short as IS10 (1.4 kb) (Kleckner 1977). The circular transduced fragment is able to recombine with homologous sequences in the chromosome. Circularization and insertion of transduced fragments allow addition of segments to the bacterial chromosome rather than replacement of recipient segments as in a normal transductional cross. It also provides a method for translocation of bacterial genes to a variety of specific sites on the chromosome in either orientation. The significance of these events to bacterial evolution is discussed.     

Partial reversion at the bobbed locus of Drosophila melanogaster

In Drosophila melanogaster the tandemly arranged repetitive sequences coding for 18S and 28S rRNA are heterogenous at the level of the spacers between units and insertions that interrupt many 28S rRNA genes. This heterogeneity contrasts with the homogeneity of the regions transcribed into 18S and 28S rRNA. Homogenization and evolution of repetitive genes are usually explained by conversion, amplification events or unequal crossovers. In this paper we studied the change in rDNA patterns associated with partial reversion of bobbed mutations. In most cases, no increase in rDNA gene number, but a new repartition of gene types were found.     

Integration site preferences of the Alu family and similar repetitive DNA sequences.

Numerous flanking nucleotide sequences from two primate interspersed repetitive DNA families have been aligned to determine the integration site preferences of each repetitive family. This analysis indicates that both the human Alu and galago Monomer families were preferentially inserted into short d(A+T)-rich regions. Moreover, both primate repeat families demonstrated an orientation specific integration with respect to dA-rich sequences within the flanking direct repeats. These observations suggest that a common mechanism exists for the insertion of many repetitive DNA families into new genomic sites. A modified mechanism for site-specific integration of primate repetitive DNA sequences is provided which requires insertion into dA-rich sequences in the genome. This model is consistent with the observed relationship between galago Type II subfamilies suggesting that they have arisen not by mere mutation but by independent integration events.     

Multiple repetitive elements and organization of the lux operons of luminescent terrestrial bacteria.

The complete nucleotide sequences of the luxA to luxE genes, as well as the flanking regions, were determined for the lux operons of two Xenorhabdus luminescens strains isolated from insects and humans. The nucleotide sequences of the corresponding lux genes (luxCDABE) were 85 to 90% identical but completely diverged 350 bp upstream of the first lux gene (luxC) and immediately downstream of the last lux gene (luxE). These results show that the luxG gene found immediately downstream of luxE in luminescent marine bacteria is missing at this location in terrestrial bacteria and raise the possibility that the lux operons are at different positions in the genomes of the X. luminescens strains. Four enteric repetitive intergenic consensus (ERIC) or intergenic repetitive unit (IRU) sequences of 126 bp were identified in the 7.7-kbp DNA fragment from the X.luminescens strain isolated from humans, providing the first example of multiple ERIC structures in the same operon including two ERIC structures at the same site. Only a single ERIC structure between luxB and luxE is present in the 7-kbp lux DNA from insects. Analysis of the genomic DNAs from five X. luminescens strains or isolates by polymerase chain reaction has demonstrated that an ERIC structure is between luxB and luxE in all of the strains, whereas only the strains isolated from humans had an ERIC structure between luxD and luxA. The results indicate that there has been insertion and/or deletion of multiple 126-bp repetitive elements in the lux operons of X.luminescens during evolution.     

A rapidly evolving region in the immunoglobulin heavy chain loci of rat and mouse: postulated role of (dC-dA)n.(dG-dT)n sequences

The nucleotide sequences of the introns that are located between the C4 exon and the first membrane exon of mouse and rat immunoglobulin epsilon-chain genes have been determined. The rat intron sequence was found to contain four separate clusters of repetitive sequences all of which consisted of (dC-dA)n.(dG-dT)n dinucleotide repeats. A comparison between this chromosomal region in mouse and rat revealed four deletions or duplications, three of which have occurred inside or at the borders of the CA clusters. Rearrangements have occurred inside or at the borders of all four repeats after the evolutionary separation of mouse and rat. The sequence comparison reveals in addition a duplication, connected to the CA repeats, which has occurred early in evolution, before the evolutionary divergence of mouse and rat. These findings suggest that (dC-dA)n.(dG-dT)n sequences are potential targets for recombination events.