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

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

Young but not relatively old retrotransposons are preferentially located in gene‐rich euchromatic regions in tomato (Solanum lycopersicum) plants

Long terminal repeat (LTR) retrotransposons are the major DNA components of flowering plants. They are generally enriched in pericentromeric heterochromatin regions of their host genomes, which could result from the preferential insertion of LTR retrotransposons and the low effectiveness of purifying selection in these regions. To estimate the relative importance of the actions of these two factors on their distribution pattern, the LTR retrotransposons in Solanum lycopersicum (tomato) plants were characterized at the genome level, and then the distribution of young elements was compared with that of relatively old elements. The current data show that old elements are mainly located in recombination‐suppressed heterochromatin regions, and that young elements are preferentially located in the gene‐rich euchromatic regions. Further analysis showed a negative correlation between the insertion time of LTR retrotransposons and the recombination rate. The data also showed there to be more solo LTRs in genic regions than in intergenic regions or in regions close to genes. These observations indicate that, unlike in many other plant genomes, the current LTR retrotransposons in tomatoes have a tendency to be preferentially located into euchromatic regions, probably caused by their severe suppression of activities in heterochromatic regions. These elements are apt to be maintained in heterochromatin regions, probably as a consequence of the pericentromeric effect in tomatoes. These results also indicate that local recombination rates and intensities of purifying selection in different genomic regions are largely responsible for structural variation and non‐random distribution of LTR retrotransposons in tomato plants.     

Transposons but not retrotransposons are located preferentially in regions of high recombination rate in Caenorhabditis elegans

We analyzed the distribution of transposable elements (TEs: transposons, LTR retrotransposons, and non-LTR retrotransposons) in the chromosomes of the nematode Caenorhabditis elegans. The density of transposons (DNA-based elements) along the chromosomes was found to be positively correlated with recombination rate, but this relationship was not observed for LTR or non-LTR retrotransposons (RNA-based elements). Gene (coding region) density is higher in regions of low recombination rate. However, the lower TE density in these regions is not due to the counterselection of TE insertions within exons since the same positive correlation between TE density and recombination rate was found in noncoding regions (both in introns and intergenic DNA). These data are not compatible with a global model of selection acting against TE insertions, for which an accumulation of elements in regions of reduced recombination is expected. We also found no evidence for a stronger selection against TE insertions on the X chromosome compared to the autosomes. The difference in distribution of the DNA and RNA-based elements along the chromosomes in relation to recombination rate can be explained by differences in the transposition processes.     

Patterns of insertion and deletion in contrasting chromatin domains

Transposable elements (TEs) play a fundamental role in the evolution of genomes. In Drosophila they are disproportionately represented in regions of low recombination, such as in heterochromatin. This pattern has been attributed to selection against repeated elements in regions of normal recombination, owing to either (1) the slightly deleterious position effects of TE insertions near or into genes, or (2) strong selection against chromosomal abnormalities arising from ectopic exchange between TE repeats. We have used defective non-long-terminal repeat (LTR) TEs that are "dead-on-arrival" (DOA) and unable to transpose in order to estimate spontaneous deletion rates in different constituents of chromatin. These elements have previously provided evidence for an extremely high rate of spontaneous deletion in Drosophila as compared with mammals, potentially explaining at least part of the differences in the genome sizes in these organisms. However, rates of deletion could be overestimated due to positive selection for a smaller likelihood of ectopic exchange. In this article, we show that rates of spontaneous deletion in DOA repeats are as high in heterochromatin and regions of euchromatin with low recombination as they are in regions of euchromatin with normal recombination. We have also examined the age distribution of five non-LTR families throughout the genome. We show that there is substantial variation in the historical pattern of transposition of these TEs. The overrepresentation of TEs in the heterochromatin is primarily due to their longer retention time in heterochromatin, as evidenced by the average time since insertion. Fragments inserted recently are much more evenly distributed in the genome. This contrast demonstrates that the accumulation of TEs in heterochromatin and in euchromatic regions of low recombination is not due to biased transposition but by greater probabilities of fixation in these regions relative to regions of normal recombination.     

Size matters: non-LTR retrotransposable elements and ectopic recombination in Drosophila

The Drosophila melanogaster genome contains approximately 100 distinct families of transposable elements (TEs). In the euchromatic part of the genome, each family is present in a small number of copies (5-150 copies), with individual copies of TEs often present at very low frequencies in populations. This pattern is likely to reflect a balance between the inflow of TEs by transposition and the removal of TEs by natural selection. The nature of natural selection acting against TEs remains controversial. We provide evidence that selection against chromosome abnormalities caused by ectopic recombination limits the spread of some TEs. We also demonstrate for the first time that some TE families in the Drosophila euchromatin appear to be only marginally affected by purifying selection and contain many copies at high population frequencies. We argue that TEs in these families attain high population frequencies and even reach fixation as a result of low family-wide transposition rates leading to low TE copy numbers and consequently reduced strength of selection acting on individual TE copies. Fixation of TEs in these families should provide an upward pressure on the size of intergenic sequences counterbalancing rapid DNA loss through small deletions. Copy-number-dependent selection on TE families caused by ectopic recombination may also promote diversity among TEs in the Drosophila genome.     

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.     

Selection at linked sites in the partial selfer Caenorhabditis elegans

Natural selection can produce a correlation between local recombination rates and levels of neutral DNA polymorphism as a consequence of genetic hitchhiking and background selection. Theory suggests that selection at linked sites should affect patterns of neutral variation in partially selfing populations more dramatically than in outcrossing populations. However, empirical investigations of selection at linked sites have focused primarily on outcrossing species. To assess the potential role of selection as a determinant of neutral polymorphism in the context of partial self-fertilization, we conducted a multivariate analysis of single-nucleotide polymorphism (SNP) density throughout the genome of the nematode Caenorhabditis elegans. We based the analysis on a published SNP data set and partitioned the genome into windows to calculate SNP densities, recombination rates, and gene densities across all six chromosomes. Our analyses identify a strong, positive correlation between recombination rate and neutral polymorphism (as estimated by noncoding SNP density) across the genome of C. elegans. Furthermore, we find that levels of neutral polymorphism are lower in gene-dense regions than in gene-poor regions in some analyses. Analyses incorporating local estimates of divergence between C. elegans and C. briggsae indicate that a mutational explanation alone is unlikely to explain the observed patterns. Consequently, we interpret these findings as evidence that natural selection shapes genome-wide patterns of neutral polymorphism in C. elegans. Our study provides the first demonstration of such an effect in a partially selfing animal. Explicit models of genetic hitchhiking and background selection can each adequately describe the relationship between recombination rate and SNP density, but only when they incorporate selfing rate. Clarification of the relative roles of genetic hitchhiking and background selection in C. elegans awaits the development of specific theoretical predictions that account for partial self-fertilization and biased sex ratios.     

The Impact of Recombination Hotspots on Genome Evolution of a Fungal Plant Pathogen

Recombination has an impact on genome evolution by maintaining chromosomal integrity, affecting the efficacy of selection, and increasing genetic variability in populations. Recombination rates are a key determinant of the coevolutionary dynamics between hosts and their pathogens. Historic recombination events created devastating new pathogens, but the impact of ongoing recombination in sexual pathogens is poorly understood. Many fungal pathogens of plants undergo regular sexual cycles, and sex is considered to be a major factor contributing to virulence. We generated a recombination map at kilobase-scale resolution for the haploid plant pathogenic fungus Zymoseptoria tritici. To account for intraspecific variation in recombination rates, we constructed genetic maps from two independent crosses. We localized a total of 10,287 crossover events in 441 progeny and found that recombination rates were highly heterogeneous within and among chromosomes. Recombination rates on large chromosomes were inversely correlated with chromosome length. Short accessory chromosomes often lacked evidence for crossovers between parental chromosomes. Recombination was concentrated in narrow hotspots that were preferentially located close to telomeres. Hotspots were only partially conserved between the two crosses, suggesting that hotspots are short-lived and may vary according to genomic background. Genes located in hotspot regions were enriched in genes encoding secreted proteins. Population resequencing showed that chromosomal regions with high recombination rates were strongly correlated with regions of low linkage disequilibrium. Hence, genes in pathogen recombination hotspots are likely to evolve faster in natural populations and may represent a greater threat to the host.     

Evolution of the phosphoglycerate mutase processed gene in human and chimpanzee revealing the origin of a new primate gene

Processed genes are created by retroposition from messenger RNA of expressed genes. The estimated amount of processed copies of genes in the human genome is 10,000-14,000. Some of these might be pseudogenes with the expected pattern for nonfunctional sequences, but some others might be an important source of new genes. We have studied the evolution of a Phosphoglycerate mutase processed gene (PGAM3) described in humans and believed to be a pseudogene. We sequenced PGAM3 in chimpanzee and macaque and obtained polymorphism data for human coding region. We found evidence that PGAM3 likely produces a functional protein, as an example of addressing functionality for human processed pseudogenes. First, the open reading frame was intact despite many deletions that occurred in the 3' untranslated region. Second, it appears that the gene is expressed. Finally, interspecies and intraspecies variation for PGAM3 was not consistent with the neutral model proposed for pseudogenes, suggesting that a new functional primate gene has originated. Amino acid divergence was significantly higher than synonymous divergence in PGAM3 lineage, supporting positive selection acting in this gene. This role of selection was further supported by the excess of rare alleles in a population genetic analysis. PGAM3 is located in a region of very low recombination; therefore, it is conceivable that the rapid fixation events in this newly arising gene may have contributed to a selective sweep of variation in the region.     

人类基因组上的假基因

假基因是基因组上与编码基因序列非常相似的非功能性基因组DNA拷贝,一般情况都不被转录,且没有明确生理意义。假基因根据其来源可分为复制假基因和已加工假基因。迄今为止,明确鉴定的人类假基因多为已加工假基因,有8000个之多。在Swiss-Prot/TrEMBL收录的编码蛋白质的将近25500个基因序列中,约10％在基因组中有一个或多个近全长已加工假基因。其余的功能基因都没有已加工假基因。核糖体蛋白基因具有最多数量的已加工假基因,约有l700个(占已加工假基因数的22％),少数基因,如cyclophilinA、肌动蛋白(actin)、角蛋白(keratin)、GAPDH、细胞色素C(cytochromec)和nucleophosmin等则有很多份已加工假基因。总体上讲,假基因在人类染色体上的分布与染色体长度成比例,但已加工假基因在GC含量为41％～46％的染色体区域密度最高。已加工假基因的拷贝数和功能基因在生殖器官中的表达高度一致,说明许多假基因发生在胚胎阶段,另外也和基因中GC含量和基因大小密切相关。假基因的准确鉴定对基因组进化、分子医学研究和医学应用具有重要意义。     

Patterns of selection against transposons inferred from the distribution of Tc1, Tc3 and Tc5 insertions in the mut-7 line of the nematode Caenorhabditis elegans

To identify the factors (selective or mutational) that affect the distribution of transposable elements (TEs) within a genome, it is necessary to compare the pattern of newly arising element insertions to the pattern of element insertions that have been fixed in a population. To do this, we analyzed the distribution of recent mutant insertions of the Tc1, Tc3, and Tc5 elements in a mut-7 background of the nematode Caenorhabditis elegans and compared it to the distribution of element insertions (presumably fixed) within the sequenced genome. Tc1 elements preferentially insert in regions with high recombination rates, whereas Tc3 and Tc5 do not. Although Tc1 and Tc3 both insert in TA dinucleotides, there is no clear relationship between the frequency of insertions and the TA dinucleotide density. There is a strong selection against TE insertions within coding regions: the probability that a TE will be fixed is at least 31 times lower in coding regions than in noncoding regions. Contrary to the prediction of theoretical models, we found that the selective pressure against TE insertions does not increase with the recombination rate. These findings indicate that the distribution of these three transposon families in the genome of C. elegans is determined essentially by just two factors: the pattern of insertions, which is a characteristic of each family, and the selection against insertions within coding regions.     

Natural selection in gene-dense regions shapes the genomic pattern of polymorphism in wild and domesticated rice

Levels of nucleotide variability are frequently positively correlated with recombination rate and negatively associated with gene density due to the effects of selection on linked variation. These relationships are determined by properties that frequently differ among species, including the mating system, and aspects of genome organization such as how genes are distributed along chromosomes. In rice, genes are found at highest density in regions with frequent crossing-over. This association between gene density and recombination rate provides an opportunity to evaluate the effects of selection in a genomic context that differs from other model organisms. Using single-nucleotide polymorphism data from Asian domesticated rice Oryza sativa ssp. japonica and ssp. indica and their progenitor species O. rufipogon, we observe a significant negative association between levels of polymorphism and both gene and coding site density, but either no association, or a negative correlation, between nucleotide variability and recombination rate. We establish that these patterns are unlikely to be explained by neutral mutation rate biases and demonstrate that a model of background selection with variable rates of deleterious mutation is sufficient to account for the gene density effect in O. rufipogon. In O. sativa ssp. japonica, we report a strong negative correlation between polymorphism and recombination rate and greater losses of variation during domestication in the euchromatic chromosome arms than heterochromatin. This is consistent with Hill-Robertson interference in low-recombination regions, which may limit the efficacy of selection for domestication traits. Our results suggest that the physical distribution of selected mutations is a primary factor that determines the genomic pattern of polymorphism in wild and domesticated rice species.     

Effects of recombination rate on human endogenous retrovirus fixation and persistence

Endogenous retroviruses (ERVs) result from germ line infections by exogenous retroviruses. They can proliferate within the genome of their host species until they are either inactivated by mutation or removed by recombinational deletion. ERVs belong to a diverse group of mobile genetic elements collectively termed transposable elements (TEs). Numerous studies have attempted to elucidate the factors determining the genomic distribution and persistence of TEs. Here we show that, within humans, gene density and not recombination rate correlates with fixation of endogenous retroviruses, whereas the local recombination rate determines their persistence in a full-length state. Recombination does not appear to influence fixation either via the ectopic exchange model or by indirect models based on the efficacy of selection. We propose a model linking rates of meiotic recombination to the probability of recombinational deletion to explain the effect of recombination rate on persistence. Chromosomes 19 and Y are exceptions, possessing more elements than other regions, and we suggest this is due to low gene density and elevated rates of human ERV integration in males for chromosome Y and segmental duplication for chromosome 19.     

Antigen Diversity in the Parasitic Bacterium Anaplasma phagocytophilum Arises from Selectively-Represented,Spatially Clustered Functional Pseudogenes

Anaplasma phagocytophilum is a tick-transmitted bacterial pathogen of humans and other animals, and is an obligate intracellular parasite. Throughout the course of infection, hosts acquire temporary resistance to granulocytic anaplasmosis as they develop immunity specific for the major antigen, major surface protein 2 (Msp2). However, the bacterium then utilizes a novel recombination mechanism shuffling functional pseudogenes sequentially into an expression cassette with conserved 5′ and 3′ ends, bypassing host immunity. Approximately 100 pseudogenes are present in the only fully sequenced human-origin HZ genome, representing the possibility for almost unlimited antigenic diversity. In the present study, we identified a select group of 20% of the A. phagocytophilum HZ msp2 pseudogenes that have matched preferentially to human, canine, and equine expression cassettes. Pseudogenes cluster predominantly in one spatial run limited to a single genomic island in less than 50% of the genome but phylogenetically related pseudogenes are neither necessarily located in close proximity on the genome nor share similar percent identity with expression cassettes. Pseudogenes near the expression cassette (and the origin) are more likely to be expressed than those farther away. Taken together, these findings suggest that there may be natural selection pressure to retain pseudogenes in one cluster near the putative origin of replication, even though global recombination shuffles pseudogenes around the genome, separating pseudogenes that share genetic origins as well as those with similar identities.