Unusual evolutionary conservation of 5S rRNA pseudogenes inAspergillus nidulans: Similarity of the DNA sequence associated with the pseudogenes with the mouse immunoglobulin switch region

Summary AllAspergillus nidulans 5S rRNA pseudogenes known so far are the result of integration of an approx. 0.2-kbp-long DNA sequence into the 5S rRNA genes. This sequence, called block C, is present in at least five copies in theA. nidulans genome and seems to be associated either with 5S rRNA genes or pseudogenes. In contrast to the 78% sequence conservation of the C-block in pseudogenes, the truncated 5 halves of the pseudogenes are very highly conserved (96.9–100%). We postulate that the 5S rRNA pseudogenes are still a subject of concerted evolution. The C-block sequence shows similarity to the switch region of the mouse heavy chain immunoglobulin gene. A characteristic motif GGGTGAG is repeated several times in both sequences; the sequence conservation is 63%.     

The size distribution of insertions and deletions in human and rodent pseudogenes suggests the logarithmic gap penalty for sequence alignment

The size distributions of deletions, insertions, and indels (i.e., insertions or deletions) were studied, using 78 human processed pseudogenes and other published data sets. The following results were obtained: (1) Deletions occur more frequently than do insertions in sequence evolution; none of the pseudogenes studied shows significantly more insertions than deletions. (2) Empirically, the size distributions of deletions, insertions, and indels can be described well by a power law, i.e., f _k = Ck ^–b , where f _k is the frequency of deletion, insertion, or indel with gap length k, b is the power parameter, and C is the normalization factor. (3) The estimates of b for deletions and insertions from the same data set are approximately equal to each other, indicating that the size distributions for deletions and insertions are approximately identical. (4) The variation in the estimates of b among various data sets is small, indicating that the effect of local structure exists but only plays a secondary role in the size distribution of deletions and insertions. (5) The linear gap penalty, which is most commonly used in sequence alignment, is not supported by our analysis; rather, the power law for the size distribution of indels suggests that an appropriate gap penalty is w _k = a + b ln k, where a is the gap creation cost and blnk is the gap extension cost. (6) The higher frequency of deletion over insertion suggests that the gap creation cost of insertion (a _i ) should be larger than that of deletion (a _d ); that is, a _i – a _d = In R, where R is the frequency ratio of deletions to insertions. Correspondence to: W.-H. Li     

Mutation Pattern Variation Among Regions of the Primate Genome

We sequenced three argininosuccinate-synthetase-processed pseudogenes (ΨAS-A1, ΨAS-A3, ΨAS-3) and their noncoding flanking sequences in human, orangutan, baboon, and colobus. Our data showed that these pseudogenes were incorporated into the genome of the Old World monkeys after the divergence of the Old World and New World monkey lineages. These pseudogene flanking regions show variable mutation rates and patterns. The variation in the G/C to A/T mutation rate (u) can account for the unequal GC contents at equilibrium: 34.9, 36.9, and 41.7% in the pseudogene ΨAS-A1, ΨAS-A3, and ΨAS-3 flanking regions, respectively. The A/T to G/C mutation rate (v) seems stable and the u/v ratios equal 1.9, 1.7, and 1.4 in the flanking regions of ΨAS-A1, ΨAS-A3, and ΨAS-3, respectively. These ``regional' variations of the mutation rate affect the evolution of the pseudogenes, too. The ratio u/v being greater than 1.0 in each case, the overall mutation rate in the GC-rich pseudogenes is, as expected, higher than in their GC-poor flanking regions. Moreover, a ``sequence effect' has been found. In the three cases examined u and v are higher (at least 20%) in the pseudogene than in its flanking region—i.e., the pseudogene appears as mutation ``hot' spots embedded in ``cold' regions. This observation could be partly linked to the fact that the pseudogene flanking regions are long-standing unconstrained DNA sequences, whereas the pseudogenes were relieved of selection on their coding functions only around 30–40 million years ago. We suspect that relatively more mutable sites maintained unchanged during the evolution of the argininosuccinate gene are able to change in the pseudogenes, such sites being eliminated or rare in the flanking regions which have been void of strong selective constraints over a much longer period. Our results shed light on (1) the multiplicity of factors that tune the spontaneous mutation rate and (2) the impact of the genomic position of a sequence on its evolution. Received: 10 February 1997 / Accepted: 21 April 1997     

The influence of nearest neighbors on the rate and pattern of spontaneous point mutations

Summary The numbers and local sequence environments of the two types of substitution mutation plus additions and deletions have been obtained directly in this study from differences between a large number of extant primate gene and pseudogene sequences. A total of 3786 mutations were scored in regions where similarities between pseudogene and corresponding gene sequences is 85%, comprising 30% of the pseudogene database of 80,584 bp. The pattern of mutations obtained in this fashion is almost identical to that obtained by Li et al. (1984) using a slightly different, more direct approach and with a smaller database. When mutations were scored, the neighbor pairs on the 5 and 3 sides were also noted, leading to a large 16 × 12 matrix of transitions and transversions. Biases of varying magnitude are found in the rates of substitution of the same base pair in different local sequence environments. The overall order for the effect of the 5 neighbor on the rates of substitution mutation of a pyrimidine is A > C T > G, and G > A > T > C for the 3 neighbor; where these results represent the average of substitution rates for the complement purine with complement neighbors of bases ordered above. The order for the 3 neighbor is essentially the same for the two transitions and most of the four transversions as well; however, the order for the 5 neighbor is more variable. The overall rate for the C · G T · A transition is not unusual, however the presence of a 3 neighboring G · C pair boosts the rate substantially, presumably due to specific cytosine methylation of the CG doublet in primate DNAs. The rate of the T · A C · G transition is also well above average when the 3 neighbor is an A · T, and to a lesser extent a G · C, pair. The latter bias is typical in that it reflects the association of alternating pyrimidine-purine sequences with increasing mutation rates. The substitution of the pyrimidine in a 5 purine-pyrimi-dine-purine3 sequence generally occurs much faster than in a pyrimidine tract and points to the local conformation as a major determining factor of the substitution rate. An apparent inverse relationship is found between starting and product doublet frequencies of base pairs undergoing mutations with specific 3 neighbors, indicating that differences in intrinsic substitution rates of base pairs with specific neighbors are a key factor in producing the familiar biases of nearest-neighbor frequencies.Offprint requests to: R. D. Blake     

The Muridae glyceraldehyde-3-phosphate dehydrogenase family

Summary Although only one gene is known to be functional, numerous glyceraldehyde-3-phosphate dehydrogenase (GAPDH) related sequences are scattered throughoutMus musculus andRattus rattus genomes. In this report we show that: (1) GAPDH pseudogenes are repeated to comparable extents, at least 400 copies, in 12 other Muridae species; (2) the complete, or nearly so, sequence of GAPDH messenger RNA is amplified, and a high proportion, if not all of these copies, are intronless; (3) GAPDH pseudogenes are preferentially located in heavily methylated and DNAse I-insensitive regions of chromatin; and (4) the presence of atypical GAPDH-related mRNAs in different cellular contexts raises the possibility that more than one GAPDH gene is transcribed.     

Active MHC class Ib genes in rat are pseudogenes in the mouse

The most telomeric class I region of the MHC in rat and mouse is the M region, which contains about 20 class I genes or gene fragments. The central part carries three class I genes—M4, M5, and M6—which are orthologous between the two species. M4 and M6 are pseudogenes in the mouse but transcribed, intact genes in the rat. To analyze the pseudogene status for the mouse genes in more detail, we have sequenced the respective exons in multiple representative haplotypes. The stop codons are conserved in all mouse strains analyzed, and, consistent with the pseudogene status, all strains show additional insertions and deletions, taking the genes further away from functionality. Thus, M4 and M6 indeed have a split status. They are silent in the mouse but intact in the closely related rodent, the rat.GenBank accession numbers: AF057065 to AF057072 (exon 3 of H2-M4 of reported mouse strains), AF057976 to AF057985 (exon 3 of RT1.M4 of reported rat strains), AF058923 and AF058924 (exon 2 of RT1.M4 of strains PVG and BN), AY286080 to AY286092 (exon 4 of H2-M6 of reported mouse stains), and AY303772 (full-length genomic sequence of RT1.M6-1^l)     

Leprosy – clues about the biochemistry of Mycobacterium leprae and its host-dependency from the genome

Deletions and the appearance of pseudogenes in pathways of carbon source utilisation and energy metabolism best explain the host-dependency and failure to culture Mycobacterium leprae axenically. From the genome sequence it is possible to predict that acetate and galactose cannot be used as carbon sources, while pyruvate can only be catabolised. Glycerol, glucose, and fatty acids could be used for glycolysis, the pentose cycle and -oxidation which are complete. Retrospective functional genomics – interpreting work before the completion of the genome project – supports the failure of M. leprae to use acetate as well as another prediction that metabolic flux from pyruvate to acetyl-CoA would be very low. However, the loss of a second icd gene (compared with M. tuberculosis), predicted to encode isocitrate dehydrogenase, did not diminish the specific activity of the enzyme. The genes for respiratory pathways are extremely limited, being present for oxidative phosphorylation as a result of electron transport only using FADH as an electron donor. In contrast, all the major biosynthetic pathways are complete except that M. leprae is a natural methionine auxotroph: this is predicted not to be attenuating, or explain host-dependency since methionine would be present in rich culture media.     

Local hopping mobile DNA implicated in pseudogene formation and reductive evolution in an obligate cyanobacteria-plant symbiosis

Background

Insertion sequences (ISs) are approximately 1 kbp long “jumping” genes found in prokaryotes. ISs encode the protein Transposase, which facilitates the excision and reinsertion of ISs in genomes, making these sequences a type of class I (“cut-and-paste”) Mobile Genetic Elements. ISs are proposed to be involved in the reductive evolution of symbiotic prokaryotes. Our previous sequencing of the genome of the cyanobacterium ‘Nostoc azollae’ 0708, living in a tight perpetual symbiotic association with a plant (the water fern Azolla), revealed the presence of an eroding genome, with a high number of insertion sequences (ISs) together with an unprecedented large proportion of pseudogenes. To investigate the role of ISs in the reductive evolution of ‘Nostoc azollae’ 0708, and potentially in the formation of pseudogenes, a bioinformatic investigation of the IS identities and positions in 47 cyanobacterial genomes was conducted. To widen the scope, the IS contents were analysed qualitatively and quantitatively in 20 other genomes representing both free-living and symbiotic bacteria.

Results

Insertion Sequences were not randomly distributed in the bacterial genomes and were found to transpose short distances from their original location (“local hopping”) and pseudogenes were enriched in the vicinity of IS elements. In general, symbiotic organisms showed higher densities of IS elements and pseudogenes than non-symbiotic bacteria. A total of 1108 distinct repeated sequences over 500 bp were identified in the 67 genomes investigated. In the genome of ‘Nostoc azollae’ 0708, IS elements were apparent at 970 locations (14.3%), with 428 being full-length. Morphologically complex cyanobacteria with large genomes showed higher frequencies of IS elements, irrespective of life style.

Conclusions

The apparent co-location of IS elements and pseudogenes found in prokaryotic genomes implies earlier IS transpositions into genes. As transpositions tend to be local rather than genome wide this likely explains the proximity between IS elements and pseudogenes. These findings suggest that ISs facilitate the reductive evolution in for instance in the symbiotic cyanobacterium ‘Nostoc azollae’ 0708 and in other obligate prokaryotic symbionts.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1386-7) contains supplementary material, which is available to authorized users.     

Study of molecular genetic diversity and evolutionary history of medicinally important endangered genus Chlorophytum using DNA barcodes

This study was aimed to authenticate and present phylogenetic relationship among 19 species of genus Chlorophytum using DNA barcoding. In all, 107 accessions were analyzed with eight plastid (matK, rbcL, trnH-psbA, rpoC1, ycf5, rpoB, atp and psbK-psbI) and six nuclear (ITS) markers. The matK and rbcL were found to be ideal markers for identification and discrimination of Chlorophytum species. Phylogenetic analysis based on matK and rbcL sequences resolved the species in two major clades. All markers, except matK and rbcL, showed ambiguous reads and paralogy in analysis. DGGE analysis showed the presence of pseudogenes and/or co-amplification in these markers, which caused poor sequence quality. Phylogeny and probable evolution of genus Chlorophytum was proposed on the basis of cytological, morphological and genetic information.