Recombination drives the evolution of GC-content in the human genome

Unraveling the evolutionary forces responsible for variations of neutral substitution patterns among taxa or along genomes is a major issue in the identification of functional sequence features. Mammalian genomes show large-scale regional variations of GC-content (the isochores), but the substitution processes at the origin of this structure are poorly understood. We have analyzed the pattern of neutral substitutions in 14.3 Mb of primate noncoding regions. We show that the GC-content toward which sequences are evolving is strongly correlated (r(2) = 0.61, P 相似文献   

Context dependent substitution biases vary within the human genome

Background  

Models of sequence evolution typically assume that different nucleotide positions evolve independently. This assumption is widely appreciated to be an over-simplification. The best known violations involve biases due to adjacent nucleotides. There have also been suggestions that biases exist at larger scales, however this possibility has not been systematically explored.     

Phase-dependent nucleotide substitution in protein-coding sequences

It is well known that due to the degeneracy of genetic code, most of the silent substitutions appear in the third codon position, so the mutation frequency of the third codon position is much higher than that of the first two positions. However, it remains unknown whether the directionality of point mutation in three codon positions is similar or not. In this paper, through analyzing 15 sets of orthologous genes, it is revealed that most of the substitution types are significantly different between any two codon positions, especially between the 2nd and the 3rd phases. Furthermore, the average frequencies of each type of substitution calculated from the fifteen sets of orthologous genes are similar to those identified in single nucleotide polymorphisms (SNPs) of human and mouse genome. The present analyses suggest that the nucleotide substitution in protein-coding sequences is not only context-dependent (so called neighboring-nucleotide effects), but also phase-dependent, which is of significance to improving the prevalent nucleotide-evolution models.     

Compositional evolution of noncoding DNA in the human and chimpanzee genomes

We have examined the compositional evolution of noncoding DNA in the primate genome by comparison of lineage-specific substitutions observed in 1.8 Mb of genomic alignments of human, chimpanzee, and baboon with 6542 human single-nucleotide polymorphisms (SNPs) rooted using chimpanzee sequence. The pattern of compositional evolution, measured in terms of the numbers of GC-->AT and AT-->GC changes, differs significantly between fixed and polymorphic sites, and indicates that there is a bias toward fixation of AT-->GC mutations, which could result from weak directional selection or biased gene conversion in favor of high GC content. Comparison of the frequency distributions of a subset of the SNPs revealed no significant difference between GC-->AT and AT-->GC polymorphisms, although AT-->GC polymorphisms in regions of high GC segregate at slightly higher frequencies on average than GC-->AT polymorphisms, which is consistent with a fixation bias favoring high GC in these regions. However, the substitution data suggest that this fixation bias is relatively weak, because the compositional structure of the human and chimpanzee genomes is becoming homogenized, with regions of high GC decreasing in GC content and regions of low GC increasing in GC content. The rate and pattern of nucleotide substitution in 333 Alu repeats within the human-chimpanzee-baboon alignments are not significantly affected by the GC content of the region in which they are inserted, providing further evidence that, since the time of the human-chimpanzee ancestor, there has been little or no regional variation in mutation bias.     

Organization of nucleotide sequences in the chicken genome

The four major components of chicken DNA were prepared by density gradient centrifugation and characterized in several basic properties: relative amounts, dG + dC content, buoyant densities, compositional heterogeneity, and reassociation kinetics. While the relative amounts and the compositions of the major components of chicken DNA were similar to those found in mammalian genomes, their compositional heterogeneities were found to be narrower. The relative amounts of interspersed repeated and unique sequences were strikingly different in different components and also different from those found in the corresponding major components of mouse and human DNAs. If one takes into consideration that major DNA components (a) account for practically all of main-band DNA and (b) derive by preparative breakage from very long DNA segments of fairly homogeneous composition, the isochores, our findings indicate that the distribution of interspersed repeats is different in different chromosomal regions and is species-specific.     

Functional genes mapped on the chicken genome

Microsatellite polymorphisms are finding increasing use in genetics. In addition to the random isolation of microsatellite markers, such markers can also be developed from sequences already present in public domain databases. An advantage of public domain databases is that these microsatellites are known to be located within or close to identified functional genes. In this study the GenBank and EMBL databases were screened for microsatellite markers and primers were defined for amplification. Subsequently, these markers were tested on a panel of five different birds from layer and broiler stocks and on the international reference families: the East Lansing reference family and the Compton reference family. Of the 33 loci tested, 25 were polymorphic on the test panel and from these 25, 14 were polymorphic in one or both reference families. Twelve of the 14 loci that could be mapped fell into previously defined linkage groups. The other two markers were not linked. Because three of the loci had previously been mapped to specific chromosomes by in situ hybridization, linkage groups E6 and C3 could be assigned to chromosome 6, E5 and C17 to chromosome 4 and E21 to one of the microchromosomes.     

Evolution of the isochore structure in the scale of chromosome: insight from the mutation bias and fixation bias

Following the development of reliable methods for inferring the direction of mutations of the single nucleotide polymorphism (SNP), and the revealing of the human isochore map, it has become possible to investigate the evolution of the isochore structure in a continuous region. In this study, the recent evolution of the isochore structure on human chromosome 18, as inferred from the SNP, was examined. A remarkable mutation bias was found, which was destroying the present isochore structure. However, a fixation bias contributed by the biased gene conversion (BGC) effect and a rising fixation probability of derived alleles with increasing GC content was extending the present isochore structure. Combining the two opposing processes, the old isochore structure was declining and a more homogenous isochore structure with higher GC content was being formed on the chromosome. During this process, both the CpG and genic sites, which were present in the isochore but were paid little attention to before, played an important role. In addition, the recombination was confirmed to promote the GC alleles fixed in the genome because of the BGC effect. For the first time, it was observed that with the occurrence of little recombination, AT alleles had the identical fixation probability with GC alleles in the recombination cold spots.     

A high recombination frequency within the chicken major histocompatibility (B) complex

Chickens of a commercial pure White Leghorn line were typed for B-F and B-G by serological, biochemical and molecular biological methods. Amongst 287 typed animals of one particular line, three animals with recombinant haplotypes were identified. Compared to earlier reports this revealed a statistically significant (P < 0 –05), tenfold higher recombination frequency in this chicken line.     

Minisatellite DNA markers in the chicken genome

This paper reports the detailed characterization of multilocus minisatellite DNA fingerprints in the chicken. Results are presented of DNA fingerprint segregation analyses carried out in three chicken pedigrees, calculating the number of detected loci, testing for Mendelian inheritance, and cosegregation among fingerprint bands. Two pedigrees (families 1 and 2) were analysed using the Jeffreys probes 33.6 and 33.15 only, and one pedigree (family 3) was analysed using 33.6, 33.15. 3′α-globin HVR and M13 protein III gene repeat. Mean band transmission frequencies in families 1 and 2 were near to the Mendelian expectation of 0.5 and no mutations were observed. Family 3 showed transmission frequencies slightly exceeding 0.5. Linkage among bands was higher than observed in some other avian species, with each allele represented by a mean of 1.48 HaeIII fragments. Cosegregation of heterozygous parental fragments representing distinguishable loci followed the expected binomial distribution. The number of minisatellites detectable by the four probes was estimated to be 217. The pattern of cosegregation among those minisatellite loci was tested against that expected for different levels of recombination through the use of a simulation model. We conclude that most minisatellites are unlinked and probably widely dispersed in the chicken genome.     

Parameters of the chicken genome (Gallus gallus)

As more information on the chicken genome is gathered, it is becoming increasingly more important to be able to correlate genetic and physical maps. Quantitation of the chicken karyotype is important in establishing parameters which define the genome. Here we report on the physical lengths of the chicken macrochromosomes and establish the DNA content of each, thus identifying implicitly how much of the genome is represented by the microchromosomal component. For the first time, genetic and physical data on the chicken karyotype are presented in relation to one another.     

Distinct changes of genomic biases in nucleotide substitution at the time of Mammalian radiation

Differences in the regional substitution patterns in the human genome created patterns of large-scale variation of base composition known as genomic isochores. To gain insight into the origin of the genomic isochores, we develop a maximum-likelihood approach to determine the history of substitution patterns in the human genome. This approach utilizes the vast amount of repetitive sequence deposited in the human genome over the past approximately 250 Myr. Using this approach, we estimate the frequencies of seven types of substitutions: the four transversions, two transitions, and the methyl-assisted transition of cytosine in CpG. Comparing substitutional patterns in repetitive elements of various ages, we reconstruct the history of the base-substitutional process in the different isochores for the past 250 Myr. At around 90 MYA (around the time of the mammalian radiation), we find an abrupt fourfold to eightfold increase of the cytosine transition rate in CpG pairs compared with that of the reptilian ancestor. Further analysis of nucleotide substitutions in regions with different GC content reveals concurrent changes in the substitutional patterns. Although the substitutional pattern was dependent on the regional GC content in such ways that it preserved the regional GC content before the mammalian radiation, it lost this dependence afterward. The substitutional pattern changed from an isochore-preserving to an isochore-degrading one. We conclude that isochores have been established before the radiation of the eutherian mammals and have been subject to the process of homogenization since then.     

Heterogeneity in regional GC content and differential usage of codons and amino acids in GC-poor and GC-rich regions of the genome of Apis mellifera

The honeybee (Apis mellifera) has a genome with a wide variation in GC content showing 2 clear modal GC values, in some ways reminiscent of an isochore-like structure. To gain insight into causes and consequences of this pattern, we used a comparative approach to study the genome-wide alignment of primarily coding sequence of A. mellifera with Drosophila melanogaster and Anopheles gambiae. The latter 2 species show a higher average GC content than A. mellifera and no indications of bimodality, suggesting that the GC-poor mode is a derived condition in honeybee. In A. mellifera, synonymous sites of genes generally adopt the GC content of the region in which they reside. A large proportion of genes in GC-poor regions have not been assigned to the honeybee assembly because of the low sequence complexity of their genome neighborhood. The synonymous substitution rate between A. mellifera and the other species is very close to saturation, but analyses of nonsynonymous substitutions as well as amino acid substitutions indicate that the GC-poor regions are not evolving faster than the GC-rich regions. We describe the codon usage and amino acid usage and show that they are remarkably heterogeneous within the honeybee genome between the 2 different GC regions. Specifically, the genes located in GC-poor regions show a much larger deviation in both codon usage bias and amino acid usage from the Dipterans than the genes located in the GC-rich regions.     

Environmental effects on gene expression phenotype have regional biases in the human genome

Phenotypic discordance between monozygotic twins, such as a difference in disease susceptibility, implicates the role of the environment in determining phenotype. To assess genomewide environmental effects on "gene expression phenotype," we employed a published microarray data set for twins. We found that variations in expression phenotypes between monozygotic twins have biases in their chromosomal locations. They also showed a strong inverse correlation with gene density. Genomic regions of low gene density were environmentally sensitive, containing genes involved in response to external signals, cell differentiation, and development, etc. Genetic factors were found to make no contribution to the observed regional biases, stressing the role of epigenetics. We propose that epigenetic modifications might occur more frequently in heterochromatic, gene-poor regions in response to environmental signals while gene-rich regions tend to remain in an active chromatin configuration for the constitutive expression of underlying genes.     

Autosomal genetic control of the activity of a new variant ornithine transcarbamylase in chicken kidney

The mode of inheritance of the gene for chick kidney ornithine transcarbamylase (OTC), found previously as a genetic variant, was investigated. White Leghorn B line males homozygous for the allele for the variant OTC gene were selected using the California Gray breed, having a near-absolute deficiency of the enzyme. Then further crosses of the two breeds were made. The mean value of the OTC level of F₁ progeny was about 170 units. Chicks from the backcross generation were divided into two groups, of high activity and low activity, in a ratio of 1:1. F₂ chicks were divided into three groups: one-fourth of the chicks were classified as a super high group, one-half were high, and the remaining one-fourth were low the mean values for OTC level were 356.7, 196.4, and 15.6 units, respectively. From these results, it was suggested that the variant OTC represents a simple autosomal incompletely dominant trait.This work was supported in part by a grant-in-aid (No. 012207) for the scientific research from the Ministry of Education, Science, and Culture of Japan.     

A gradient of silent substitution rate in the human pseudoautosomal region

It has been demonstrated that recombination in the human p-arm pseudoautosomal region (p-PAR) is at least twenty times more frequent than the genomic average of approximately 1 cM/Mb, which may affect substitution patterns and rates in this region. Here I report the analysis of substitution patterns and rates in 10 human, chimpanzee, gorilla, and orangutan genes across the p-PAR. Between species silent divergence in the p-PAR forms a gradient, increasing toward the telomere. The correlation of silent divergence with distance from the p-PAR boundary is highly significant (rho = 0.911, P < 0.001). After exclusion of the CpG dinucleotides this correlation is still significant (rho = 0.89, P < 0.01), thus the substitution rate gradient cannot be explained solely by the differences in the extent of methylation across the p-PAR. Frequent recombination in the PAR may result in a relatively strong effect of biased gene conversion (BGC), which, because of the increased probability of fixation of the G or C nucleotides at (A or T)/(G or C) segregating sites, may affect substitution rates. BGC, however, does not seem to be the factor creating the substitution rate gradient in the p-PAR, because the only gradient is still detactable if only A<-->T and G<-->C substitutions are taken into account (rho = 0.82, P < 0.01). I hypothesize that the substitution rate gradient in the p-PAR is due to the mutagenic effect of recombination, which is very frequent in the distal human p-PAR and might be lower near the p-PAR boundary.     

Patterns of nucleotide substitution in pseudogenes and functional genes

Summary The pattern of point mutations is inferred from nucleotide substitutions in pseudogenes. The pattern obtained suggests that transition mutations occur somewhat more frequently than transversion mutations and that mutations result more often in A or T than in G or C. Our results are discussed with respect to the predictions from Topal and Fresco's model for the molecular basis of point (substitution) mutations (Nature 263:285–289, 1976). The pattern of nucleotide substitution at the first and second positions of codons in functional genes is quite similar to that in pseudogenes, but the relative frequency of the transition CT in the sense strand is drastically reduced and those of the transversions CG and GC are doubled. The differences between the two patterns can be explained by the observation that in the protein evolution amino acid substitutions occur mainly between amino acids with similar biochemical properties (Grantham, Science 185:862–864, 1974). Our results for the patterns of nucleotide substitutions in pseudogenes and in functional genes lead to the prediction that both the coding and non-coding regions of protein coding genes should have high frequencies of A and T. Available data show that the non-coding regions are indeed high in A and T but the coding regions are low in T, though high in A.     

Single nucleotide polymorphisms in chicken lmbr1 gene were associated with chicken growth and carcass traits

Lmbr1 is the key candidate gene controlling vertebrate limb development, but its effects on animal growth and carcass traits have never been reported. In this experiment, lmbr1 was taken as the candi-date gene affecting chicken growth and carcass traits. T/C and G/A mutations located in exon 16 and one A/C mutation located in intron 5 of chicken lmbr1 were detected from Silky, White Plymouth Rock broilers and their F2 crossing chickens by PCR-SSCP and sequencing methods. The analysis of vari-ance (ANOVA) results suggests that T/C polymorphism of exon 16 had significant association with eviscerated yield rate (EYR), gizzard rate (GR), shank and claw rate (SCR) and shank girth (SG); A/C polymorphism of intron 5 was significantly associated with SCR, liver rate and head-neck weight (HNW), while both sites had no significant association with other growth and carcass traits. These results demonstrate that lmbr1 gene could be a genetic locus or linked to a major gene significantly affecting these growth and carcass traits in chicken.     

Enhancer analysis by chicken embryo electroporation with aid of genome comparison

The identification of the enhancers associated with each developmentally regulated gene is a first step to clarify the regulatory mechanisms underlying embryogenesis. The electroporation technique using chicken embryo is a powerful tool to identify such enhancers. The technique enables us to survey a large genomic region and to analyze the enhancers in great detail. Comparison of the genomic sequences of the chicken and other vertebrate species identifies conserved non-coding sequence blocks to which the functionally identified enhancers often correspond. In this review, I describe in detail the methods to analyze the enhancers using the chicken embryo electroporation and genome comparison.     

Using the chicken genome sequence in the development and mapping of genetic markers in the turkey (Meleagris gallopavo)

The efficacy of employing the chicken genome sequence in developing genetic markers and in mapping the turkey genome was studied. Eighty previously uncharacterized microsatellite markers were identified for the turkey using BLAST alignment to the chicken genome. The chicken sequence was then used to develop primers for polymerase chain reaction where the turkey sequence was either unavailable or insufficient. A total of 78 primer sets were tested for amplification and polymorphism in the turkey, and informative markers were genetically mapped. Sixty-five (83%) amplified turkey genomic DNA, and 33 (42%) were polymorphic in the University of Minnesota/Nicholas Turkey Breeding Farms mapping families. All but one marker genetically mapped to the position predicted from the chicken genome sequence. These results demonstrate the usefulness of the chicken sequence for the development of genomic resources in other avian species.