The effect of relaxed functional constraints on the photosynthetic gene rbcL in photosynthetic and nonphotosynthetic parasitic plants

The photosynthetic gene rbcL has been lost or dramatically altered in some lineages of nonphotosynthetic parasitic plants, but the dynamics of these events following loss of photosynthesis and whether rbcL has sustained functionally significant changes in photosynthetic parasitic plants are unknown. To assess the changes to rbcL associated with the loss of functional constraints for photosynthesis, nucleotide sequences from nonparasitic and parasitic plants of Scrophulariales were used for phylogeny reconstruction and character analysis. Plants in this group display a broad range of parasitic abilities, from photosynthetic ("hemiparasites") to nonphotosynthetic ("holoparasites"). With the exception of Conopholis (Orobanchaceae), the rbcL locus is present in all parasitic plants of Scrophulariales examined. Several holoparasitic genera included in this study, including Boschniakia, Epifagus, Orobanche, and Hyobanche, have rbcL pseudogenes. However, the holoparasites Alectra orobanchoides, Harveya capensis, Harveya purpurea, Lathraea clandestina, Orobanche corymbosa, O. fasciculata, and Striga gesnerioides have intact open reading frames (ORFs) for the rbcL gene. Phylogenetic hypotheses based on rbcL are largely in agreement with those based on sequences of the nonphotosynthetic genes rps2 and matK and show a single origin of parasitism, and loss of photosynthesis and pseudogene formation have been independently derived several times in Scrophulariales. The mutations in rbcL in nonparasitic and hemiparasitic plants would result in largely conservative amino acid substitutions, supporting the hypothesis that functional proteins can experience only a limited range of changes, even in minimally photosynthetic plants. In contrast, ORFs in some holoparasites had many previously unobserved missense substitutions at functionally important amino acid residues, suggesting that rbcL genes in these plants have evolved under relaxed or altered functional constraints.      

Molecular evolution of rbcL in the mycoheterotrophic coralroot orchids (Corallorhiza Gagnebin, Orchidaceae)

The RuBisCO large subunit gene (rbcL) has been the focus of numerous plant phylogenetic studies and studies on molecular evolution in parasitic plants. However, there has been a lack of investigation of photosynthesis gene molecular evolution in fully mycoheterotrophic plants. These plants invade pre-existing mutualistic associations between ectomycorrhizal trees and fungi, from which they obtain fixed carbon and nutrients. The mycoheterotrophic orchid Corallorhiza contains both green (photosynthetic) and non-green (putatively nonphotosynthetic) species. We sequenced rbcL from 31 accessions of eight species of Corallorhiza and hypothesized that some lineages would have pseudogenes resulting from relaxation of purifying selection on RuBisCO's carboxylase function. Phylogenetic analysis of rbcL+ITS gave high jackknife support for relationships among species. We found evidence of pseudogene formation in all lineages of the Corallorhiza striata complex and in some lineages of the C. maculata complex. Evidence includes: stop codons, frameshifts, decreased d(S)/d(N) ratios, replacements not observed in photosynthetic species, rate heterogeneity, and high likelihood of neutral evolution. The evolution of rbcL in Corallorhiza may serve as an exemplary system in which to study the effects of relaxed evolutionary constraints on photosynthesis genes for >400 documented fully mycoheterotrophic plant species.     

Alternate paths of evolution for the photosynthetic gene rbcL in four nonphotosynthetic species of Orobanche

We have determined the nucleotide sequence for the Rubisco large subunit from four holoparasitic species of Orobanche. Intact open reading frames are present in two species (O. corymbosa and O. fasciculata), whereas the remaining species (O. cernua and O. ramosa) have rbcL pseudogenes. Sequences for rbcL 5'-UTRs from species of Orobanche have few changes in the promoter and ribosome binding sites compared to photosynthetic higher plants. Comparison of rbcL 3'-UTR sequences for Nicotiana, Ipomoea, Cuscuta, and Orobanche reveal that nucleotide sequences from parasitic plants have regions capable of forming stem-loop structures, but 56–69 nt are deleted upstream of the stem-loop in the parasitic plants compared to their photosynthetic relatives. Although rbcL pseudogenes of O. cernua and O ramosa have many large and small deletions, few indels are shared in common, implying that their common ancestor probably had an intact rbcL reading frame. Intact rbcL reading frames in O. corymbosa and O. fasciculata retain a bias of synonymous over nonsynonymous substitutions and deduced protein sequences are consistent with potentially functional Rubisco large subunit proteins. A conservative model of random substitution processes in pseudogene sequences estimates that the probability is low (P<0.028) that these sequences would retain an open reading frame by chance. Species of Orobanche have either had recent photosynthetic ancestors, implying multiple independent losses of photosynthesis in this genus, or the rbcL gene may serve an unknown function in some nonphotosynthetic plants.     

Evolutionary Aspects of Functional and Pseudogene Members of the Phytochrome Gene Family in Scots Pine

According to the neutral theory of evolution, mutation and genetic drift are the only forces that shape unconstrained, neutral, gene evolution. Thus, pseudogenes (which often evolve neutrally) provide opportunities to obtain direct estimates of mutation rates that are not biased by selection, and gene families comprising functional and pseudogene members provide useful material for both estimating neutral mutation rates and identifying sites that appear to be under positive or negative selection pressures. Conifers could be very useful for such analyses since they have large and complex genomes. There is evidence that pseudogenes make significant contributions to the size and complexity of gene families in pines, although few studies have examined the composition and evolution of gene families in conifers. In this work, I examine the complexity and rates of mutation of the phytochrome gene family in Pinus sylvestris and show that it includes not only functional genes but also pseudogenes. As expected, the functional PHYO does not appear to have evolved neutrally, while phytochrome pseudogenes show signs of unconstrained evolution.     

Pseudogene evolution and natural selection for a compact genome

Pseudogenes are nonfunctional copies of protein-coding genes that are presumed to evolve without selective constraints on their coding function. They are of considerable utility in evolutionary genetics because, in the absence of selection, different types of mutations in pseudogenes should have equal probabilities of fixation. This theoretical inference justifies the estimation of patterns of spontaneous mutation from the analysis of patterns of substitutions in pseudogenes. Although it is possible to test whether pseudogene sequences evolve without constraints for their protein-coding function, it is much more difficult to ascertain whether pseudogenes may affect fitness in ways unrelated to their nucleotide sequence. Consider the possibility that a pseudogene affects fitness merely by increasing genome size. If a larger genome is deleterious--for example, because of increased energetic costs associated with genome replication and maintenance--then deletions, which decrease the length of a pseudogene, should be selectively advantageous relative to insertions or nucleotide substitutions. In this article we examine the implications of selection for genome size relative to small (1-400 bp) deletions, in light of empirical evidence pertaining to the size distribution of deletions observed in Drosophila and mammalian pseudogenes. There is a large difference in the deletion spectra between these organisms. We argue that this difference cannot easily be attributed to selection for overall genome size, since the magnitude of selection is unlikely to be strong enough to significantly affect the probability of fixation of small deletions in Drosophila.     

Decay of mutualistic potential in aphid endosymbionts through silencing of biosynthetic loci: Buchnera of Diuraphis

Buchnera, the primary bacterial endosymbiont of aphids, is known to provision essential amino acids lacking in the hosts' diet of plant sap. The recent discovery of silenced copies of genes for tryptophan biosynthesis (trpEG) in certain Buchnera lineages suggests a decay in symbiotic functions in some aphid species. However, neither the distribution of pseudogenes among lineages nor the impact of this gene silencing on amino-acid availability in hosts has been assessed. In Buchnera of the aphid Diuraphis noxia, tandem repeats of these pseudogenes have persisted in diverse lineages, and thpEG pseudogenes have originated at least twice within this aphid genus. Measures of amino-acid concentrations in Diuraphis species have shown that the presence of the pseudogene is associated with a decreased availability of tryptophan, indicating that gene silencing decreases nutrient provisioning by symbionts. In Buchnera of Diuraphis, rates of nonsynonymous substitutions are elevated in functional trpE copies, supporting the hypothesis that pseudogene origin and persistence reflect a reduced selection for symbiont biosynthetic contributions. The parallel evolution of trpEG pseudogenes in Buchnera of Diuraphis and certain other aphid hosts suggests that either selection at the host level is not effective or that fitness in these aphids is not limited by tryptophan availability.     

A maximum likelihood method for analyzing pseudogene evolution: implications for silent site evolution in humans and rodents.

We present a new likelihood method for detecting constrained evolution at synonymous sites and other forms of nonneutral evolution in putative pseudogenes. The model is applicable whenever the DNA sequence is available from a protein-coding functional gene, a pseudogene derived from the protein-coding gene, and an orthologous functional copy of the gene. Two nested likelihood ratio tests are developed to test the hypotheses that (1) the putative pseudogene has equal rates of silent and replacement substitutions; and (2) the rate of synonymous substitution in the functional gene equals the rate of substitution in the pseudogene. The method is applied to a data set containing 74 human processed-pseudogene loci, 25 mouse processed-pseudogene loci, and 22 rat processed-pseudogene loci. Using the informatics resources of the Human Genome Project, we localized 67 of the human-pseudogene pairs in the genome and estimated the GC content of a large surrounding genomic region for each. We find that, for pseudogenes deposited in GC regions similar to those of their paralogs, the assumption of equal rates of silent and replacement site evolution in the pseudogene is upheld; in these cases, the rate of silent site evolution in the functional genes is approximately 70% the rate of evolution in the pseudogene. On the other hand, for pseudogenes located in genomic regions of much lower GC than their functional gene, we see a sharp increase in the rate of silent site substitutions, leading to a large rate of rejection for the pseudogene equality likelihood ratio test.     

RT-PCR检测HUTP14A mRNA时排除假基因HUTP14C干扰的方法

人类U3蛋白14C基因(HUTP14C)是人类U3蛋白14A基因(HUTP14A)的假基因。两者转录本序列同源性高达95%。常规RT-qPCR技术在检测HUTP14A mRNA丰度时,HUTP14C的存在会影响检测结果。本研究旨在建立检测HUTP14A mRNA时排除HUTP14C干扰的RT PCR方法。本研究设计出能分别从多种肿瘤细胞DNA和RNA中特异性扩增HUTP14A和HUTP14C的引物,避免假基因HUTP14C对其同源基因HUTP14A检测的干扰。在检测细胞系HUTP14A mRNA时,通过DNaseⅠ消除RNA中污染的HUTP14C DNA,用靶向HUTP14C 3′-UTR的siRNA沉默HUTP14C mRNA后,再用RT PCR检测HUTP14A mRNA丰度,使结果更加准确。在18对肝癌及癌旁组织中,利用特异性引物进行RT PCR检测,HUTP14A和HUTP14C mRNA的表达略高于癌旁组织。本研究提示,针对有假基因存在的功能基因,对其mRNA丰度进行检测时,在提取细胞或组织总RNA后,用DNaseⅠ处理,再用RNA直接进行PCR扩增,排除DNA污染后,再进行RT-PCR或RT-qPCR扩增。大多假基因具有较长的3′-UTR区,在该区域设计siRNA特异性沉默假基因的mRNA后,用RT-qPCR检测功能基因的mRNA丰度,可以排除假基因mRNA的影响。在病理组织中检测功能基因的mRNA丰度时,可以根据假基因和其功能基因的序列差异设计出特异扩增功能基因的引物,从假基因的3′-UTR区设计特异扩增假基因的引物,通过RT-qPCR技术分别检测二者的mRNA。     

人类基因组中加工假基因分布与重组率和基因密度的关系

分析了人类加工假基因在染色体上的分布,发现加工假基因密度与重组率负相关,而与基因密度正相关。加工假基因在低重组区的积累与插入有害模型和异位重组模型相吻合：在插入有害模型下,低重组区的选择强度由于Hill．Robertson干涉而变弱,所以加工假基因较多地插入到低重组区;在异位重组模型下,同源加工假基因家族（包括同源祖先基因）之内可能发生异位重组而对机体造成危害,所以加工假基因在高重组区的插入受到较强的负选择,导致加工假基因较多地分布在低重组区。除以上两种模型以外,加工假基因还可能通过降低重组率的方式对加工假基因密度与重组率的负相关有所贡献。加工假基因偏好分布在基因密区,这可能与异位重组在该区较少发生有关。     

Identification of three human pseudogenes for subunit VIb of cytochrome c oxidase: a molecular record of gene evolution.

Three pseudogenes for the nuclear-encoded subunit VIb of cytochrome c oxidase (COX) were isolated by screening a human genomic library with cloned human cDNA coding for COX subunit VIb. The nucleotide sequences of the pseudogenes, designated psi COX6b-1, psi COX6b-2 and psi COX6b-3, were determined. Pseudogene psi COX6b-1 bears all the hallmarks of a processed pseudogene and diverged from the parental gene after the divergence of man and cow. Alu repetitive elements were integrated into the structural sequences of the other two pseudogenes. Comparison with the human and bovine cDNA sequences encoding COX subunit VIb suggests that psi COX6b-2 and psi COX6b-3 were formed earlier in evolution than psi COX6b-1. Genomic Southern analysis indicated that a few more pseudogenes for COX subunit VIb are likely to be present in the human genome. Identical nt differences with respect to the human cDNA sequence in the pseudogenes provide some clues on the evolution of the ancestral gene coding for COX subunit VIb.     

The recent evolution of a pseudogene: diversity and divergence of a mitochondria-localized small heat shock protein in Arabidopsis thaliana

In this study we examined the evolution of the genes for three organelle-localized small heat shock proteins in Arabidopsis thaliana: the chloroplast-localized (CP) protein HSP21 and two mitochondria-localized (MT) proteins, HSP23.5 and HSP23.6. We found that the CP protein and one of the MT proteins, HSP23.6, are evolving under purifying selection to maintain function. In contrast, the gene for HSP23.5, the other MT protein, is highly variable within A. thaliana, and in some accessions or ecotypes this gene may be a pseudogene. HSP23.5 and HSP23.6 are related via a segmental duplication event, and the presence of orthologs of each gene in other species within the Brassicaceae indicates that the duplication generating HSP23.5 and HSP23.6 may have occurred as much as 20 million years ago. This is considerably longer than the 4 million year half-life of gene duplicates (functional genes as well as pseudogenes) reported by some studies. Our results are consistent with the prediction that after gene duplication one gene duplicate can be maintained for some time under relaxed selection while it accumulates random mutations. By capturing a pseudogene in the making our study provides important information on how pseudogenes are formed.     

Differences in pyrenoid morphology are correlated with differences in the rbcL genes of members of the Chloromonas lineage (volvocales,chlorophyceae)

Chloromonas is distinguished from Chlamydomonas primarily by the absence of pyrenoids, which are structures that are present in the chloroplasts of most algae and are composed primarily of the CO2-fixing enzyme Rubisco. In this study we compared sequences of the rbcL (Rubisco large subunit-encoding) genes of pyrenoid-less Chloromonas species with those of closely related pyrenoid-containing Chlamydomonas species in the "Chloromonas lineage" and with those of 45 other green algae. We found that the proteins encoded by the rbcL genes had a much higher level of amino acid substitution in members of the Chloromonas lineage than they did in other algae. This kind of elevated substitution rate was not observed, however, in the deduced proteins encoded by two other chloroplast genes that we analyzed: atpB and psaB. The rates of synonymous and nonsynonymous nucleotide substitutions in the rbcL genes indicate that the rapid evolution of these genes in members of the Chloromonas lineage is not due to relaxed selection (as it preasumably is in parasitic land plants). A phylogenetic tree based on rbcL nucleotide sequences nested two Chlamydomonas species as a "pyrenoid-regained" clade within a monophyletic Chloromonas "pyrenoid-lost" clade. Character-state optimization with this tree suggested that the loss and the regain of pyrenoids were accompanied by eight synapomorphic amino acid replacements in the Rubisco large subunit, four of which are positioned in the region involved in its dimerization. However, both the atpB and the psaB sequence data gave robust support for a rather different set of phylogenetic relationships in which neither the "pyrenoid-lost" nor the "pyrenoid-regained" clade was resolved. The appearance of such clades in the rbcL-based tree may be an artifact of convergent evolutionary changes that have occurred in a region of the large subunit that determines whether Rubisco molecules will aggregate to form a visible pyrenoid.     

A possible role of genetic translation ambiguity in evolution

A hypothesis has been suggested that the variability of translation machinery is one of the key factors of evolutionary transformations of genetic material. It considers a module principle of the evolution theory based on the concepts of duplication and divergence of genetic material, which is required for origination of new genes and proteins with new functions. The duplication results in the appearance of pseudogenes, functionally inactive, but serving a material for creating new functions. The possible mechanisms changing the translation machinery have been considered, which may lead to the sporadic pseudogene activation "supplying" natural selection with mutational changes accumulated by pseudogenes to assess their adaptive value. This takes into account not only potential possibilities of mutational variability of the translation machinery, but also the possibility of protein prioritization: a prion mechanism of inheritance is also considered which is intensively studied nowadays.     

Comparative Analysis of Evolution in a Rodent Histone H2a Pseudogene

Sequences were obtained from five species of rodents that are orthologous to an H2a histone pseudogene from Mus musculus. The pseudogene is part of the cluster of replication-dependent histone genes found on Mus musculus chromosome 13. Comparative analysis of these five sequences together with the previously published sequence from M. musculus shows that this gene has likely been a pseudogene throughout the evolution of the genus Mus, while the gene from Rattus norvegicus is likely functional. Three large (>20 bp) deletions were found among the Mus pseudogenes, a feature that is very unusual compared to surveys of processed pseudogenes. In addition, there are two single-base deletions and one 4-bp insertion among the Mus pseudogenes. The species distributions of one of the large deletions and the 4-bp insertion require either independent insertions of an identical sequence, independent deletions with identical boundaries, or a deletion followed by precise reintegration of the original sequence. The evidence favors the hypothesis of multiple deletions with identical boundaries. The ``coding' regions of the Mus pseudogenes show a much reduced level of among-species variability in the 3′ half of the pseudogene, compared both to the 5′ half and to flanking sequences. This supports a hypothesis that the 3′ end of the pseudogene is the target of frequent gene conversion by functional H2a genes. Received: 1 April 1997 / Accepted: 12 June 1997