Identifying DNA polymorphisms in human TCRA/D variable genes by direct sequencing of PCR products

 The T-cell receptor (TCR) is a highly variable molecule composed of two polypeptide chains that recognize antigenic peptides in the context of major histocompatibility complex (MHC) molecules. In this study, we describe a sequence-based search for germline polymorphisms in the variable (V) gene segments of the human TCRA/D locus. Thirty different V gene segments were amplified from six to eight unrelated individuals and sequenced from low melting point agarose. Twenty-seven polymorphisms were identified in 15 V gene segments. These polymorphisms are mainly single nucleotide substitutions, but an insertion/deletion polymorphism and a single dinucleotide repeat with variable length were also seen. Of the 15 sequence variations found in the coding regions, six are silent and nine encode amino acid changes. All of the amino acid changes are found at non-conserved residues, frequently in the hypervariable regions, where they may influence MHC and/or peptide recognition. Therefore, it is possible that germline variations in TCR genes could influence an individual’s immune response, and may also contribute to susceptibility to diseases such as autoimmunity. Received: 9 January 1996 / Revised: 22 February 1996     

High level of functional polymorphism indicates a unique role of natural selection at human immune system loci

Several studies have shown that immune system proteins have on average a higher rate of amino acid evolution between different species of mammals than do most other proteins. To test whether immune-system-expressed loci show a correspondingly elevated rate of within-species nonsynonymous (amino acid altering) polymorphism, we examined gene diversity (heterozygosity) at 4,911 single nucleotide polymorphism (SNP) sites at 481 protein-coding loci. At loci with nonimmune functions, gene diversity at nonsynonymous SNP sites was typically lower than that at silent SNP sites (those not altering the amino acid sequence) in the same gene, a pattern that is an evidence of purifying selection acting to eliminate slightly deleterious variants. However, this pattern was not seen at nonsynonymous SNPs causing conservative amino acid replacements in immune system proteins, indicating that the latter are subject to a reduced level of functional constraint. Similarly, immune system genes showed higher gene diversities in their 5′ noncoding regions than did other proteins. These results identified certain immune system loci that are likely to be subject to balancing selection that acts to maintain polymorphism in either coding or regulatory regions. Electronic Supplementary Material Supplementary material is available for this article at .     

The relationship between allozyme and chromosomal polymorphism inferred from nucleotide variation at the Acph-1 gene region of Drosophila subobscura.

The Acph-1 gene region was sequenced in 51 lines of Drosophila subobscura. Lines differ in their chromosomal arrangement for segment I of the O chromosome (O(st) and O(3+4)) and in the Acph-1 electrophoretic allele (Acph-1(100), Acph-1(054), and Acph-1(>100)). The ACPH-1 protein exhibits much more variation than previously detected by electrophoresis. The amino acid replacements responsible for the Acph-1(054) and Acph-1(>100) electrophoretic variants are different within O(st) and within O(3+4), which invalidates all previous studies on linkage disequilibrium between chromosomal and allozyme polymorphisms at this locus. The Acph-1(>100) allele within O(3+4) has a recent origin, while both Acph-1(054) alleles are rather old. Levels of nucleotide variation are higher within the O(3+4) than within the O(st) arrangement except for nonsynonymous sites. The McDonald and Kreitman test shows a significant excess of nonsynonymous polymorphisms within O(st) when D. guanche is used as the outgroup. According to the nearly neutral model of molecular evolution, this excess is consistent with a smaller effective size of O(st) relative to O(3+4) arrangements. A smaller population size, a lower recombination, and a more recent bottleneck might be contributing to the smaller effective size of O(st).     

Hotspots of Biased Nucleotide Substitutions in Human Genes

Genes that have experienced accelerated evolutionary rates on the human lineage during recent evolution are candidates for involvement in human-specific adaptations. To determine the forces that cause increased evolutionary rates in certain genes, we analyzed alignments of 10,238 human genes to their orthologues in chimpanzee and macaque. Using a likelihood ratio test, we identified protein-coding sequences with an accelerated rate of base substitutions along the human lineage. Exons evolving at a fast rate in humans have a significant tendency to contain clusters of AT-to-GC (weak-to-strong) biased substitutions. This pattern is also observed in noncoding sequence flanking rapidly evolving exons. Accelerated exons occur in regions with elevated male recombination rates and exhibit an excess of nonsynonymous substitutions relative to the genomic average. We next analyzed genes with significantly elevated ratios of nonsynonymous to synonymous rates of base substitution (d_N/d_S) along the human lineage, and those with an excess of amino acid replacement substitutions relative to human polymorphism. These genes also show evidence of clusters of weak-to-strong biased substitutions. These findings indicate that a recombination-associated process, such as biased gene conversion (BGC), is driving fixation of GC alleles in the human genome. This process can lead to accelerated evolution in coding sequences and excess amino acid replacement substitutions, thereby generating significant results for tests of positive selection.     

Genetic Diversity in Exon 2 of the Major Histocompatibility Complex Class II DQB1 Locus in Nigerian Goats

The DQB1 locus is located in the major histocompatibility complex (MHC) class II region and involved in immune response. We identified 20 polymorphic sites in a 228 bp fragment of exon 2, one of the most critical regions of the MHC DQB1 gene, in 60 Nigerian goats. Four sites are located in the peptide binding region, and 10 amino acid substitutions are peculiar to Nigerian goats, compared with published sequences. A significantly higher ratio of nonsynonymous/synonymous substitutions (d _N/d _S) suggests that allelic sequence evolution is driven by balancing selection (P < 0.01). In silico functional analysis using PANTHER predicted that substitution P56R, with a subPSEC score of ?4.00629 (P_deleterious = 0.73229), is harmful to protein function. The phylogenetic tree from consensus sequences placed the two northern breeds closer to each other than either was to the southern goats. This first report of sequence diversity at the DQB1 locus for any African goat breed may be useful in the search for disease-resistant genotypes.     

Bayesian analysis suggests that most amino acid replacements in Drosophila are driven by positive selection

One of the principal goals of population genetics is to understand the processes by which genetic variation within species (polymorphism) becomes converted into genetic differences between species (divergence). In this transformation, selective neutrality, near neutrality, and positive selection may each play a role, differing from one gene to the next. Synonymous nucleotide sites are often used as a uniform standard of comparison across genes on the grounds that synonymous sites are subject to relatively weak selective constraints and so may, to a first approximation, be regarded as neutral. Synonymous sites are also interdigitated with nonsynonymous sites and so are affected equally by genomic context and demographic factors. Hence a comparison of levels of polymorphism and divergence between synonymous sites and amino acid replacement sites in a gene is potentially informative about the magnitude of selective forces associated with amino acid replacements. We have analyzed 56 genes in which polymorphism data from D. simulans are compared with divergence from a reference strain of D. melanogaster. The framework of the analysis is Bayesian and assumes that the distribution of selective effects (Malthusian fitnesses) is Gaussian with a mean that differs for each gene. In such a model, the average scaled selection intensity (gamma = N(e)s) of amino acid replacements eligible to become polymorphic or fixed is -7.31, and the standard deviation of selective effects within each locus is 6.79 (assuming homoscedasticity across loci). For newly arising mutations of this type that occur in autosomal or X-linked genes, the average proportion of beneficial mutations is 19.7%. Among the amino acid polymorphisms in the sample, the expected average proportion of beneficial mutations is 47.7%, and among amino acid replacements that become fixed the average proportion of beneficial mutations is 94.3%. The average scaled selection intensity of fixed mutations is +5.1. The presence of positive selection is pervasive with the single exception of kl-5, a Y-linked fertility gene. We find no evidence that a significant fraction of fixed amino acid replacements is neutral or nearly neutral or that positive selection drives amino acid replacements at only a subset of the loci. These results are model dependent and we discuss possible modifications of the model that might allow more neutral and nearly neutral amino acid replacements to be fixed.     

Testing for the Footprint of Sexually Antagonistic Polymorphisms in the Pseudoautosomal Region of a Plant Sex Chromosome Pair

The existence of sexually antagonistic (SA) polymorphism is widely considered the most likely explanation for the evolution of suppressed recombination of sex chromosome pairs. This explanation is largely untested empirically, and no such polymorphisms have been identified, other than in fish, where no evidence directly implicates these genes in events causing loss of recombination. We tested for the presence of loci with SA polymorphism in the plant Silene latifolia, which is dioecious (with separate male and female individuals) and has a pair of highly heteromorphic sex chromosomes, with XY males. Suppressed recombination between much of the Y and X sex chromosomes evolved in several steps, and the results in Bergero et al. (2013) show that it is still ongoing in the recombining or pseudoautosomal, regions (PARs) of these chromosomes. We used molecular evolutionary approaches to test for the footprints of SA polymorphisms, based on sequence diversity levels in S. latifolia PAR genes identified by genetic mapping. Nucleotide diversity is high for at least four of six PAR genes identified, and our data suggest the existence of polymorphisms maintained by balancing selection in this genome region, since molecular evolutionary (HKA) tests exclude an elevated mutation rate, and other tests also suggest balancing selection. The presence of sexually antagonistic alleles at a locus or loci in the PAR is suggested by the very different X and Y chromosome allele frequencies for at least one PAR gene.     

Sequence polymorphisms in <Emphasis Type="Italic">Zmisa2</Emphasis> gene are significantly associated with starch pasting and gelatinization properties in maize (<Emphasis Type="Italic">Zea mays</Emphasis> L.)

Researches on the genetic basis of starch pasting and gelatinization properties will provide foundation of maize improvement for quality, feed and industrial applications. Maize gene Zmisa2, encoding an isoamylase-type starch debranching enzyme, plays important roles in starch biosynthesis. In this study, the genomic sequences of the gene Zmisa2 in 72 elite maize inbred lines were obtained, and the nucleotide polymorphisms and haplotype diversity were detected. In addition, seven pasting and four gelatinization properties of maize were measured for the tested inbred lines using rapid visco analyzer and differential scanning calorimeter, respectively. A total of 99 sequence variants, including 91 SNPs and 8 indels, were identified at the promoter and coding regions of this gene. Although the frequency of polymorphism in promoter region is much higher than that of coding region, the SNPs in the coding region of maize gene Zmisa2 classified this gene into 21 haplotypes, which encode 11 different ISA2 proteins. Furthermore, the association of the variants of Zmisa2 gene with maize starch pasting and gelatinization properties was estimated, and the results revealed that seven SNPs in coding region, including four nonsynonymous sites, were significantly associated with phenotypic variations of pasting time and enthalpy of transition (ΔH). These results suggested that the polymorphism in maize Zmisa2 locus could be used in molecular marker-assisted selection for improvement of quality in maize breeding programs.     

Analyses of sequence polymorphism and haplotype diversity of LEAFY genes revealed post-domestication selection in the Chinese elite maize inbred lines

Post-domestication selection refers to the artificial selection on the loci controlling important agronomic traits during the process of genetic improvement in a population. The maize genes Zfl1 and Zfl2, duplicate orthologs of Arabidopsis LEAFY, are key regulators in plant branching, inflorescence and flower development, and reproduction. In this study, the full gene sequences of Zfl1 and Zfl2 from 62 Chinese elite inbred lines were amplified to evaluate their nucleotide polymorphisms and haplotype diversities. A total of 254 and 192 variants that included SNPs and indels were identified from the full sequences of Zfl1 and Zfl2, respectively. Although most of the variants were found to be located in the non-coding regions, the polymorphisms of CDS sequences classified Zfl1 into 16 haplotypes encoding 16 different proteins and Zfl2 into 18 haplotypes encoding eight different proteins. The population of Huangzaosi and its derived lines showed statistically significant signals of post-domestication selection on the Zfl1 CDS sequences, as well as lower nucleotide polymorphism and haplotype diversity than the whole set. However, the Zfl2 locus was only selected for in the heterotic group Reid. Further evidence revealed that at least 17 recombination events contributed to the genetic and haplotype diversities at the Zfl1 locus and 16 recombination events at the Zfl2 locus.     

Molecular Evolution of the TAC1 Gene from Rice （Oryza sativa L.）

Tiller angle is a key feature of the architecture of cultivated rice(Oryza sativa),since it determines planting density and influences rice yield.Our previous work identified Tiller Angle Control 1(TACl) as a major quantitative trait locus that controls rice tiller angle.To further clarify the evolutionary characterization of the TACl gene,we compared a TACl-containing 3164-bp genomic region among 113 cultivated varieties and 48 accessions of wild rice,including 43 accessions of O.rufipogon and five accessions of O.nivara.Only one single nucleotide polymorphism(SNP),a synonymous substitution,was detected in TACl coding regions of the cultivated rice varieties, whereas one synonymous and one nonsynonymous SNP were detected among the TACl coding regions of wild rice accessions.These data indicate that little natural mutation and modification in the TACl coding region occurred within the cultivated rice and its progenitor during evolution.Nucleotide diversities in the TACl gene regions of O.sativa and O.rufipogon of 0.00116 and 0.00112,respectively, further indicate that TACl has been highly conserved during the course of rice domestication.A functional nucleotide polymorphism (FNP) of TACl was only found in the japonica rice group.A neutrality test revealed strong selection,especially in the 3’-flanking region of the TACl coding region containing the FNP in the japonica rice group.However,no selection occurred in the indica and wild-rice groups.A phylogenetic tree derived from TACl sequence analysis suggests that the indica and japonica subspecies arose independently during the domestication of wild rice.     

Characterization of Mhc-DRB allelic diversity in white-tailed deer (Odocoileus virginianus) provides insight into Mhc-DRB allelic evolution within Cervidae

 Although white-tailed deer (Odocoileus virginianus) are one of North America's best studied mammals, no information is available concerning allelic diversity at any locus of the major histocompatibility complex in this taxon. Using the polymerase chain reaction, single-stranded conformation polymorphism analysis, and DNA sequencing techniques, 15 DRB exon 2 alleles were identified among 150 white-tailed deer from a single population in southeastern Oklahoma. These alleles represent a single locus and exhibit a high degree of nucleotide and amino acid polymorphism, with most amino acid variation occurring at positions forming the peptide binding sites. Furthermore, twenty-seven amino acid residues unique to white-tailed deer DRB alleles were detected, with 19 of these occurring at residues forming contact points of the peptide binding region. Significantly higher rates of nonsynonymous than synonymous substitutions were detected among these DRB alleles. In contrast to other studies of Artiodactyla DRB sequences, interallelic recombination does not appear to be playing a significant role in the generation of allelic diversity at this locus in white-tailed deer. To examine evolution of white-tailed deer (Odvi-DRB) alleles within Cervidae, we performed a phylogenetic analysis of all published red deer (Ceel-DRB), roe deer (Caca-DRB), and moose (Alal-DRB) DRB alleles. The phylogenetic tree clearly shows a trans-species persistence of DRB lineages among these taxa. Moreover, this phylogenetic tree provides insight into evolution of DRB allelic lineages within Cervidae and may aid in assignment of red deer DRB alleles to specific loci. Received: 25 June 1998 / Revised: 2 September 1998     

Characterization of chimpanzee TCRV gene polymorphism: how old are human TCRV alleles?

 The functional relevance of the majority of human T-cell receptor A and B variable region gene polymorphisms is controversial. Studies of human and nonhuman primate major histocompatibility complex (MHC) class I and II polymorphisms show that allelic lineages predate human speciation and indicate that selection favors the long-term maintenance of these advantageous mutations. We investigated at the DNA level whether 15 human TCRA and B polymorphisms exist in contemporary chimpanzee populations. Polymorphisms consisted of variable region replacements, a recombination signal sequence base change, and silent mutations. With one exception, none of these human TCR polymorphisms were observed in contemporary chimpanzees. Investigation of the same polymorphisms in a range of other nonhuman primates showed little evidence of the existence of human polymorphism prespeciation. Chimpanzee TCRAV and BV regions were however polymorphic for variation so far not observed in human groups. Levels of mitochondrial and nuclear DNA sequence variation in contemporary chimpanzees suggest that population bottlenecks have not been a feature of chimpanzee evolution and it is therefore probable that most human TCR polymorphisms have evolved in the estimated five million years since the speciation of human and chimpanzees. Thus, over the evolutionary time period studied, ancient TCRA and B polymorphisms have not been maintained by selection to the same degree as postulated for MHC polymorphisms. Received: 13 June 1997 / Received: 25 July 1997     

LEC-2, a highly variable lectin in the lichen Peltigera membranacea

Lectins are a diverse group of carbohydrate binding proteins often involved in cellular interactions. A lectin gene, lec-2, was identified in the mycobiont of the lichen Peltigera membranacea. Sequencing of lec-2 open reading frames from 21 individual samples showed an unexpectedly high level of polymorphism in the deduced protein (LEC-2), which was sorted into nine haplotypes based on amino acid sequence. Calculations showed that the rates of nonsynonymous versus synonymous nucleotide substitutions deviated significantly from the null hypothesis of neutrality, indicating strong positive selection. Molecular modeling revealed that most amino acid replacements were around the putative carbohydrate-binding pocket, indicating changes in ligand binding. Lectins have been thought to be involved in the recognition of photobiont partners in lichen symbioses, and the hypothesis that positive selection of LEC-2 is driven by variation in the Nostoc photobiont partner was tested by comparing mycobiont LEC-2 haplotypes and photobiont genotypes, as represented by the rbcLX region. It was not possible to pair up the two types of marker sequences without conflicts, suggesting that positive selection of LEC-2 was not due to variation in photobiont partners.     

Molecular analysis of intragenic recombination at the tryptophan synthetase locus in Neurospora crassa

Fifteen different classically generated and mapped mutations at the tryptophan synthetase locus in Neurospora crassa have been characterized to the level of the primary sequence of the gene. This sequence analysis has demonstrated that intragenic recombination is accurate to order mutations within one open reading frame. While classic genetic analysis correctly ordered the mutations, the position of mutations characterized by gene sequence analysis was more accurate. A leaky mutation was found to have a wild-type primary sequence. The presence of unique polymorphisms in the primary sequence of the trp-3 gene from strain 861 confirms that it has a unique history relative to the other strains studied. Most strains that were previously shown to be immunologically nonreactive with antibody preparations raised against tryptophan synthetase protein were shown to have nonsense mutations. This work defines 14 alleles of the N. crassa trp-3 gene.     

Caffeic acid O-methyltransferase allelic polymorphism characterization and analysis in different maize inbred lines

The caffeic acid O-methyltransferase (COMT) gene is a gene candidate forimproving silage maize digestibility. Before assessing possible functionalconsequences of COMT gene polymorphism, the nucleotide sequence variability atthe COMT locus in 6 normal maize inbred lines was characterized and analyzed.Inaddition, the COMT bm3-3 mutation characterization was specified. Even if theCOMT gene seemed to evolve according to the predictions of the neutral model,estimates of nucleotide diversity indicated that polymorphism distributionwithin the sequence was not uniform and that nonsynonymous sites were the leastvariable. Linkage disequilibrium calculation and estimates of recombinationshowed that polymorphism was a dimorphic segregating distribution caused by alot of recombination events.     

Nucleotide polymorphism of the Srlk gene that determines salt stress tolerance in alfalfa (Medicago sativa L.)

Based on legume genome syntheny, the nucleotide sequence of Srlk gene, key role of which in response to salt stress was demonstrated for the model species Medicago truncatula, was identified in the major forage and siderate crop alfalfa (Medicago sativa). In twelve alfalfa samples originating from regions with contrasting growing conditions, 19 SNPs were revealed in the Srlk gene. For two nonsynonymous SNPs, molecular markers were designed that could be further used to analyze the association between Srlk gene nucleotide polymorphism and the variability in salt stress tolerance among alfalfa cultivars.     

Rare McArdle disease locus polymorphic site on 11q13 contains CpG sequence

Summary When probes throughout the McArdle disease (myophosphorylase) gene region were used to search for DNA polymorphisms, only an MspI polymorphism was found in 94 enzyme-probe combinations. Along with an insertion/deletion polymorphism more 3 to the gene locus, these polymorphisms will be informative in 75% of at-risk patients. These results contrast strikingly to the six polymorphic sites detected in 15 enzyme-probe combinations in the homologous Her's disease (liver phosphorylase) gene region. This single MspI polymorphic site includes a CpG sequence of known increased mutability suggesting that DNA regions with rare polymorphisms will have most polymorphic sites at sequences with enhanced mutability. Fluorescence in situ hybridization sublocalized this gene to proximal band 11q13, establishing a point of cross-reference between the physical and genetic maps.     

An Efficient Method of Selectable Marker Gene Excision by Xer Recombination for Gene Replacement in Bacterial Chromosomes

A simple, effective method of unlabeled, stable gene insertion into bacterial chromosomes has been developed. This utilizes an insertion cassette consisting of an antibiotic resistance gene flanked by dif sites and regions homologous to the chromosomal target locus. dif is the recognition sequence for the native Xer site-specific recombinases responsible for chromosome and plasmid dimer resolution: XerC/XerD in Escherichia coli and RipX/CodV in Bacillus subtilis. Following integration of the insertion cassette into the chromosomal target locus by homologous recombination, these recombinases act to resolve the two directly repeated dif sites to a single site, thus excising the antibiotic resistance gene. Previous approaches have required the inclusion of exogenous site-specific recombinases or transposases in trans; our strategy demonstrates that this is unnecessary, since an effective recombination system is already present in bacteria. The high recombination frequency makes the inclusion of a counter-selectable marker gene unnecessary.