Survey of polymorphic sequence variation in the immediate 5' region of human DNA repair genes

Systematic screens have revealed extensive DNA sequence variation existing in the human population. Studies of the role of polymorphic genetic variants in explaining the association of family history with risk of common disease have generally focused on variants predicted to disrupt protein structure and activity. Recent studies have identified genetic variation in the level of expression of many genes, variation that is potentially biologically relevant in explaining individual variation in disease risk. In a survey of data available for 108 DNA repair genes that have been systematically screened for sequence variation, an average of 3.3 SNPs per gene were found to exist at a variant allele frequency of at least 0.02 in the region 2kb upstream from the 5'-untranslated region. One-third of the genes harbored a SNP with an allele frequency of at least 0.02 within a predicted promotor element. These variants are distributed among promoter elements that average 20 elements per gene. The frequency of polymorphic SNPs in CpG islands was 0.8 per gene, while the frequency of SNPs in the 5'-UTR was 0.7 per gene. The recognition of extensive genetic variation with potential to impact levels of gene expression, and thereby exacerbate the impact of amino acid substitution variants on the activity of proteins, increases the complexity of analyses required to explain the molecular genetic basis for the familial contribution to the sporadic incidence of common disease.     

Allelic variation in the intron 2 and 3 of the MICA gene.

We have determined by sequencing the allelic variation in intron 2 and 3 of hte MICA gene for a total of 22 different alleles. Sequencing of introns was performed in two directions, using DNA from homozygous cell lines from families and from unrelated individuals. Intron 2 is 273 bp long and did in the alleles investigated not reveal any length polymorphism. We found a total of eight polymorphic positions which exhibit a strict biallelism, as it is also found in the polymorphisms for exon 2, 3 and 4 of MICA. Intron 3 is 586 bp long an required an additional set of primers placed near the middle of this intron in order to allow a complete bidirectional sequence. In intron 3, a total of 10 polymorphic positions were identified. Interestingly, we found two variants of the allele MICA*002 which are distinguished only by one basepair difference in intron 3. The variant MICA*002A is associated with HLA-B35 and B58, while the allele MICA*002B is associated with B38 and B39.     

Structure of allelic variants of subtype 5 of histone H1 in pea Pisum sativum L

The pea genome contains seven histone H1 genes encoding different subtypes. Previously, the DNA sequence of only one gene, His1, coding for the subtype H1-1, had been identified. We isolated a histone H1 allele from a pea genomic DNA library. Data from the electrophoretic mobility of the pea H1 subtypes and their N-bromosuccinimide cleavage products indicated that the newly isolated gene corresponded to the H1-5 subtype encoded by His5. We confirmed this result by sequencing the gene from three pea lines with H1-5 allelic variants of altered electrophoretic mobility. The allele of the slow H1-5 variant differed from the standard allele by a nucleotide substitution that caused the replacement of the positively charged lysine with asparagine in the DNA-interacting domain of the histone molecule. A temperature-related occurrence had previously been demonstrated for this H1-5 variant in a study on a worldwide collection of pea germplasm. The variant tended to occur at higher frequencies in geographic regions with a cold climate. The fast allelic variant of H1-5 displayed a deletion resulting in the loss of a duplicated pentapeptide in the C-terminal domain.     

Mitochondrial base excision repair of uracil and AP sites takes place by single-nucleotide insertion and long-patch DNA synthesis

Base excision repair (BER) corrects a variety of small base lesions in DNA. The UNG gene encodes both the nuclear (UNG2) and the mitochondrial (UNG1) forms of the human uracil-DNA glycosylase (UDG). We prepared mitochondrial extracts free of nuclear BER proteins from human cell lines. Using these extracts we show that UNG is the only detectable UDG in mitochondria, and mitochondrial BER (mtBER) of uracil and AP sites occur by both single-nucleotide insertion and long-patch repair DNA synthesis. Importantly, extracts of mitochondria carry out repair of modified AP sites which in nuclei occurs through long-patch BER. Such lesions may be rather prevalent in mitochondrial DNA because of its proximity to the electron transport chain, the primary site of production of reactive oxygen species. Furthermore, mitochondrial extracts remove 5' protruding flaps from DNA which can be formed during long-patch BER, by a "flap endonuclease like" activity, although flap endonuclease (FEN1) is not present in mitochondria. In conclusion, combined short- and long-patch BER activities enable mitochondria to repair a broader range of lesions in mtDNA than previously known.     

Identification of base substitutions in ten types of rare variants of phosphoglucomutase-1 (PGM1) encountered in Japanese.

In a previous starch-gel electrophoresis study of erythrocyte phosphoglucomutase-1 (PGM1) in 23,095 Japanese from Hiroshima and Nagasaki, we detected 14 types of rare variant alleles. To determine sequence differences in these rare alleles, cell lines were established from peripheral B-lymphocytes from 24 unrelated individuals in whom nine types of rare variants are presumed to exist on the basis of earlier electrophoresis studies. cDNAs reverse transcribed from mRNAs extracted from these cell lines were amplified by polymerase chain reaction and sequences determined. Amino acid substitution types were deduced from each cDNA sequence. Although two individuals were reported to have an identical electromorph (PGM1 4HR3), sequence analysis revealed that alleles encoding these electromorphs possessed different base substitutions, and one was renamed PGM1 4HR4. As the amino acid substitution of ten different variants could be deduced by cDNA sequence in this study, the effect of each amino acid substitution on enzyme activity could be precisely simulated. The secondary structure of each variant predicted by computer simulations revealed that very decreased activity observed on PGM1 4HR2 protein was caused by significant secondary structure change introduced by the amino acid substitution. On the basis of the crystal structure, the amino acid substitutions of the ten types of rare variants seem to be outside the active center of this enzyme.     

Mitochondrial DNA Polymerase POLG1 Disease Mutations and Germline Variants Promote Tumorigenic Properties

Germline mutations in mitochondrial DNA polymerase gamma (POLG1) induce mitochondrial DNA (mtDNA) mutations, depletion, and decrease oxidative phosphorylation. Earlier, we identified somatic mutations in POLG1 and the contribution of these mutations in human cancer. However, a role for germline variations in POLG1 in human cancers is unknown. In this study, we examined a role for disease associated germline variants of POLG1, POLG1 gene expression, copy number variation and regulation in human cancers. We analyzed the mutations, expression and copy number variation in POLG1 in several cancer databases and validated the analyses in primary breast tumors and breast cancer cell lines. We discovered 5-aza-2''-deoxycytidine led epigenetic regulation of POLG1, mtDNA-encoded genes and increased mitochondrial respiration. We conducted comprehensive race based bioinformatics analyses of POLG1 gene in more than 33,000 European-Americans and 5,000 African-Americans. We identified a mitochondrial disease causing missense variation in polymerase domain of POLG1 protein at amino acid 1143 (E1143G) to be 25 times more prevalent in European-Americans (allele frequency 0.03777) when compared to African-American (allele frequency 0.00151) population. We identified T251I and P587L missense variations in exonuclease and linker region of POLG1 also to be more prevalent in European-Americans. Expression of these variants increased glucose consumption, decreased ATP production and increased matrigel invasion. Interestingly, conditional expression of these variants revealed that matrigel invasion properties conferred by these germline variants were reversible suggesting a role of epigenetic regulators. Indeed, we identified a set of miRNA whose expression was reversible after variant expression was turned off. Together, our studies demonstrate altered genetic and epigenetic regulation of POLG1 in human cancers and suggest a role for POLG1 germline variants in promoting tumorigenic properties.     

ABO Genetic Variation in Neanderthals and Denisovans

Variation at the ABO locus was one of the earliest sources of data in the study of human population identity and history, and to this day remains widely genotyped due to its importance in blood and tissue transfusions. Here, we look at ABO blood type variants in our archaic relatives: Neanderthals and Denisovans. Our goal is to understand the genetic landscape of the ABO gene in archaic humans, and how it relates to modern human ABO variation. We found two Neanderthal variants of the O allele in the Siberian Neanderthals (O1 and O2), one of these variants is shared with an European Neanderthal, who is a heterozygote for this O1 variant and a rare cis-AB variant. The Denisovan individual is heterozygous for two variants of the O1 allele, functionally similar to variants found widely in modern humans. Perhaps more surprisingly, the O2 allele variant found in Siberian Neanderthals can be found at low frequencies in modern Europeans and Southeast Asians, and the O1 allele variant found in Siberian and European Neanderthal is also found at very low frequency in modern East Asians. Our genetic distance analyses suggest both alleles survive in modern humans due to inbreeding with Neanderthals. We find that the sequence backgrounds of the surviving Neanderthal-like O alleles in modern humans retain a higher sequence divergence than other surviving Neanderthal genome fragments, supporting a view of balancing selection operating in the Neanderthal ABO alleles by retaining highly diverse haplotypes compared with portions of the genome evolving neutrally.     

Nuclear and mitochondrial forms of human uracil-DNA glycosylase are encoded by the same gene.

Recent cloning of a cDNA (UNG15) encoding human uracil-DNA glycosylase (UDG), indicated that the gene product of M(r) = 33,800 contains an N-terminal sequence of 77 amino acids not present in the presumed mature form of M(r) = 25,800. This led to the hypothesis that the N-terminal sequence might be involved in intracellular targeting. To examine this hypothesis, we analysed UDG from nuclei, mitochondria and cytosol by western blotting and high resolution gel filtration. An antibody that recognises a sequence in the mature form of the UNG protein detected all three forms, indicating that they are products of the same gene. The nuclear and mitochondrial form had an apparent M(r) = 27,500 and the cytosolic form an apparent M(r) = 38,000 by western blotting. Gel filtration gave essentially similar estimates. An antibody with specificity towards the presequence recognised the cytosolic form of M(r) = 38,000 only, indicating that the difference in size is due to the presequence. Immunofluorescence studies of HeLa cells clearly demonstrated that the major part of the UDG activity was localised in the nuclei. Transfection experiments with plasmids carrying full-length UNG15 cDNA or a truncated form of UNG15 encoding the presumed mature UNG protein demonstrated that the UNG presequence mediated sorting to the mitochondria, whereas UNG lacking the presequence was translocated to the nuclei. We conclude that the same gene encodes nuclear and mitochondrial uracil-DNA glycosylase and that the signals for mitochondrial translocation resides in the presequence, whereas signals for nuclear import are within the mature protein.     

Sequence variants in the FcepsilonRI alpha chain gene.

There is a relationship between IgE levels and expression of high-affinity IgE receptors (FcepsilonRI). Because the alpha chain is the only portion of the receptor that binds directly to IgE, we reasoned that sequence variants in the FcepsilonRI alpha gene may exist that alter these binding events. We screened all of the exons and the promoter region of the FcepsilonRI alpha chain gene with genomic DNA from 389 asthmatic and 341 normal control subjects for mutations by using single-stranded conformational polymorphism analysis. No nonsynonomous single nucleotide polymorphisms (SNPs) were identified in the coding region. Three SNPs were found in the promoter region: an A/C transversion at -770 from the translation start site; a G/A transition at -664; and a T/C transition at -335. No differences in allele frequencies were detected between asthmatic subjects and controls. Homozygosity for the C variant at locus -335 was more common in Caucasian asthmatic patients with IgE levels in the lower quartile than in the upper quartile (P = 0.032). An analysis of highly polymorphic SNPs indicated that this association is unlikely to be due to population substructure. We conclude that homozygosity for the C allele of FcepsilonRI alpha chain variant is associated with lower IgE levels.     

Characterization of a single nucleotide polymorphism in the coding sequence of the bovine transferrin gene.

A single nucleotide polymorphism was identified in the coding sequence of the bovine transferrin gene. Two alleles (SSCP1 and SSCP2) were detected by SSCP analysis and the mutation point was identified and confirmed by direct sequencing of the PCR products. The relationship between protein and DNA polymorphism was established. Protein variants A, D1 and E correspond to SSCP allele 1 and variant D2 corresponds to SSCP allele 2. DNA sequences from genotypes AA, AE, AD2, D1E, D2E and D2D2 reveal an A/G substitution at position 1455 of the cDNA which causes a Gly/Glu substitution which could be responsible for the mobility difference between D1 and D2 variants. Because of the number of variants, this suggests that other SNPs exist in the bovine transferrin gene. A linkage analysis between the SSCPs and two microsatellites (UWCA46 and CSSM019) mapped the transferrin gene to BTA1. Two-point analysis revealed a tight linkage within the transferrin protein variants and the SSCPs.     

A genome-wide comparison of the functional properties of rare and common genetic variants in humans

One of the longest running debates in evolutionary biology concerns the kind of genetic variation that is primarily responsible for phenotypic variation in species. Here, we address this question for humans specifically from the perspective of population allele frequency of variants across the complete genome, including both coding and noncoding regions. We establish simple criteria to assess the likelihood that variants are functional based on their genomic locations and then use whole-genome sequence data from 29 subjects of European origin to assess the relationship between the functional properties of variants and their population allele frequencies. We find that for all criteria used to assess the likelihood that a variant is functional, the rarer variants are significantly more likely to be functional than the more common variants. Strikingly, these patterns disappear when we focus on only those variants in which the major alleles are derived. These analyses indicate that the majority of the genetic variation in terms of phenotypic consequence may result from a mutation-selection balance, as opposed to balancing selection, and have direct relevance to the study of human disease.     

Phenylalanine hydroxylase gene: novel missense mutation in exon 7 causing severe phenylketonuria.

By direct sequence analysis of 94 mutant phenylalanine hydroxylase alleles using polymerase chain reaction-based techniques, we identified a C to T transition in exon 7 of the human phenylalanine hydroxylase gene that is associated with RFLP haplotypes 1 and 4. A leucine for proline substitution at position 281 can be predicted from the nucleotide sequence of the mutant codon. Expression analysis in cultured mammalian cells after site-directed mutagenesis proved that the base substitution is a disease causing gene lesion. Dot-blot hybridization analysis using allele-specific oligonucleotides revealed that 25% of all mutant haplotype 1 alleles in the German population bear this mutation. In addition, this mutation could be detected on one mutant haplotype 4 allele. The fact that this mutation is associated with only 25% of all mutant haplotype 1 alleles suggests that multiple mutations may be associated with this haplotype. The occurrence of several different mutations would be in agreement with the clinical heterogeneity observed in the group of patients whose PKU alleles belong to haplotype 1.     

Accurate prediction of the functional significance of single nucleotide polymorphisms and mutations in the ABCA1 gene

The human genome contains an estimated 100,000 to 300,000 DNA variants that alter an amino acid in an encoded protein. However, our ability to predict which of these variants are functionally significant is limited. We used a bioinformatics approach to define the functional significance of genetic variation in the ABCA1 gene, a cholesterol transporter crucial for the metabolism of high density lipoprotein cholesterol. To predict the functional consequence of each coding single nucleotide polymorphism and mutation in this gene, we calculated a substitution position-specific evolutionary conservation score for each variant, which considers site-specific variation among evolutionarily related proteins. To test the bioinformatics predictions experimentally, we evaluated the biochemical consequence of these sequence variants by examining the ability of cell lines stably transfected with the ABCA1 alleles to elicit cholesterol efflux. Our bioinformatics approach correctly predicted the functional impact of greater than 94% of the naturally occurring variants we assessed. The bioinformatics predictions were significantly correlated with the degree of functional impairment of ABCA1 mutations (r² = 0.62, p = 0.0008). These results have allowed us to define the impact of genetic variation on ABCA1 function and to suggest that the in silico evolutionary approach we used may be a useful tool in general for predicting the effects of DNA variation on gene function. In addition, our data suggest that considering patterns of positive selection, along with patterns of negative selection such as evolutionary conservation, may improve our ability to predict the functional effects of amino acid variation.     

Variants of the mannose-binding lectin gene in the Benin population: heterozygosity for the p.G57E allele may confer a selective advantage

Human mannose-binding lectin (MBL) plays an important role in innate immunity. MBL deficiency is associated with mutations in the promoter region and in exon 1 of the MBL2 gene. Such deficiency has been correlated with elevated incidence of infections in infancy and in immunocompromised adults. We determined the distribution profile of the MBL2 gene variants in the general population of Benin (West Africa) and in a vulnerable subset of children with sickle cell disease (SCD) (SS homozygotes). Five hundred forty-two healthy individuals (274 newborns, 268 adults) and 128 patients with SCD (35 newborns, 93 children) were screened for the common variant alleles in the MBL2 secretor haplotype region (exon 1 and promoter). The p.G57E variant allele was the most frequent allele compared to p.G54D (27.5% vs. 1.6%, respectively). The p.R52C allele was not found in this population. There was no difference in allele or genotype frequencies between healthy newborns and newborns with SCD. Alleles associated with MBL deficiency were more frequent in adults than in newborns (69.8% vs. 57.3%, respectively; p = 0.002). This enrichment was exclusively due to an elevated proportion of heterozygotes for the p.G57E allele (47.0% vs. 35.3%, respectively; p = 0.004), supporting a potential selective advantage of this genotype. Our results, compared to those reported in other African countries, support the implication of the MBL2 gene in various major infections in Africa, such as meningitis and tuberculosis in HIV-positive patients.     

Germline EPHB2 receptor variants in familial colorectal cancer

Familial clustering of colorectal cancer occurs in 15-20% of cases, however recognized cancer syndromes explain only a small fraction of this disease. Thus, the genetic basis for the majority of hereditary colorectal cancer remains unknown. EPHB2 has recently been implicated as a candidate tumor suppressor gene in colorectal cancer. The aim of this study was to evaluate the contribution of EPHB2 to hereditary colorectal cancer. We screened for germline EPHB2 sequence variants in 116 population-based familial colorectal cancer cases by DNA sequencing. We then estimated the population frequencies and characterized the biological activities of the EPHB2 variants identified. Three novel nonsynonymous missense alterations were detected. Two of these variants (A438T and G787R) result in significant residue changes, while the third leads to a conservative substitution in the carboxy-terminal SAM domain (V945I). The former two variants were found once in the 116 cases, while the V945I variant was present in 2 cases. Genotyping of additional patients with colorectal cancer and control subjects revealed that A438T and G787R represent rare EPHB2 alleles. In vitro functional studies show that the G787R substitution, located in the kinase domain, causes impaired receptor kinase activity and is therefore pathogenic, whereas the A438T variant retains its receptor function and likely represents a neutral polymorphism. Tumor tissue from the G787R variant case manifested loss of heterozygosity, with loss of the wild-type allele, supporting a tumor suppressor role for EPHB2 in rare colorectal cancer cases. Rare germline EPHB2 variants may contribute to a small fraction of hereditary colorectal cancer.     

Development of molecular approaches to estimating germinal mutation rates. I. Detection of insertion/deletion/rearrangement variants in the human genome

DNA from 130 individuals was studied with up to 18 (primarily cDNA) probes for the frequency of variants in this initial experiment to determine the feasibility of this approach to screening for germinal gene mutations. This approach, a modification of the usual restriction enzyme mapping strategy, focuses on the detection of insertion/deletion/rearrangement (I/D/R) variants, because the DNA is digested with only two restriction enzymes before transfer to membranes and hybridization with an extensive series of unrelated probes. Some 4000 noncontiguous, independent DNA fragments ("loci"), functional loci, pseudogenes or anonymous fragments, (a total of approximately 77,400 kb) were screened. 19 different classes and 31 copies of presumably I/D/R variants were detected while 4 different classes and 24 individuals exhibiting base substitution variants were observed. 18 of the 19 I/D/R classes were rare variants, that is, each were observed at a frequency, within this population, of less than 0.01; 3 of the base substitution classes existed at polymorphic frequencies and only 1 was a rare variant. 10 of the I/D/R classes, occurring in a total of 18 individuals, were detected with probes which are not known to be associated with repetitive elements. This is a variant frequency for I/D/R variants without known repetitive elements of 0.15 classes and 0.23 copies for each 1000 kb screened; this would extrapolate to 1600 such variant sites in the genome of each individual. Within the context of a mutation screening program, the rare variants, either with or without repetitive elements, would have a higher probability of being de novo mutations than would polymorphic variants; this former group would be the focus of family studies to test for the heritability of the allele (fragment pattern). Sufficient DNA probes are available to screen a significant portion of the human genome for genetic variation and de novo mutations of this type.     

Identification of genetic polymorphisms of CYP2S1 in a Finnish Caucasian population

CYP2S1 is a recently discovered member of the cytochrome P450 (CYP) gene superfamily. Interestingly, even though the DNA sequence identifies it as the sole member of the new CYP2S family, CYP2S1 exhibits many features typical to CYP1 family members, e.g. dioxin-inducibility mediated by the aryl hydrocarbon receptor (AHR) and the aryl hydrocarbon receptor nuclear translocator (ARNT). In addition, CYP2S1 metabolises some aromatic hydrocarbons as well as cellular substances. These characteristics, together with a wide extrahepatic tissue distribution, suggest that CYP2S1 may have an important role in both exogenous and endogenous metabolism. This is the first study characterising CYP2S1 alleles and naming them with the recommended CYP allele nomenclature. We used denaturing gradient gel electrophoresis (DGGE) and direct sequencing to investigate genetic variation of CYP2S1 in 100 male Finnish Caucasians. Those exons in which variation was found were examined in subsequent 100 subjects. The coding region of all of the nine exons, as well as a 449 bp fragment of the proximal promoter region, was analysed. This systematic investigation revealed eight single nucleotide polymorphisms (SNPs), which comprise nine different variant alleles (haplotypes), in addition to the wild-type allele. Seven of the SNPs occurred in the protein-coding areas and one in the proximal 3' untranslated region (3'UTR). Two of these sequence variations (10347C > T and 13106C > T) result in non-conservative amino acid substitutions, i.e. Arg380Cys and Pro466Leu, respectively. The respective allelic variants, CYP2S1*2 ([10347C > T]) and CYP2S1*3 (13106C > T; 13255A > G]), occurred in our study population at frequencies of 0.50 and 3.75%, respectively. The most common of the variant alleles was CYP2S1*1H (23.8%), harbouring a 13255A > G substitution located in the 3'UTR.     

Widespread pyrethroid resistance in Australian diamondback moth, Plutella xylostella (L.), is related to multiple mutations in the para sodium channel gene

Populations of Plutella xylostella, extending over 3800 km in southern Australia, show no genetic structure as assessed by microsatellite markers; yet outbreaks of pyrethroid resistance occur sporadically in cropping areas. Since mutations in the para voltage-gated sodium channel gene have been implicated in pyrethroid resistance, we looked for DNA sequence variation at this target among Australian moths. We found two resistance mutations previously reported for this species (L1014F and T929I), as well as a novel substitution (F1020S). Of the eight possible haplotypes formed by combinations of these three biallelic polymorphisms, only four were found in Australian populations: the wild-type allele (w), the kdr mutation allele (kdr) with only L1014F, the super-kdr-like combination of L1014F and T929I (skdrl), and the crashdown allele with only F1020S (cdr). Comparison of genotype frequencies among survivors of permethrin assays with those from untreated controls identified three resistant genotypes: skdrl homozygotes, cdr homozygotes and the corresponding heterozygote, cdr/skrdl - the heterozygote being at least as resistant as either homozygote. Spatial heterogeneity of allele frequencies was conspicuous, both across the continent and among local collections, consistent with reported spatial heterogeneity of pyrethroid resistance. Further, high resistance samples were sometimes associated with high frequency of cdr, sometimes high frequency of skdrl, or sometimes with a high combined cdr+skdrl frequency. The skdrl and cdr alleles explain a high proportion of the Australia-wide resistance variation. These data add to evidence that nerve insensitivity by mutations in the para-sodium channel gene is a common pyrethroid resistance mechanism in P. xylostella.