Rare variants of blood proteins in human populations

Rare variants of blood proteins occur, due to mutations (mutant alleles) in monomorphic loci encoding various proteins. A number of authors studied the distribution of these variants in human populations using the method of electrophoresis. The population of USA, South America, Japan, Europe was analysed. 1334 rare variants (1.0.10(-3)) were discovered out of 1,329,558 alleles (test locus in 664,779 individuals). 7 mutant alleles (3.6.10(-6)) were found among 1,957,305 alleles. The low frequency of occurrence of mutations in the loci encoding rare blood protein variants, when testing the speed of mutagenicity and its alteration, necessitates electrophoresis of blood proteins to be done in large scales. A method was proposed, based on accounting rare variants in children with congenital disorders, which are supposed to have a heavy load of mutations. The data collected demonstrated that the majority of rare variants in a given generation were obtained from parents. Accumulation of rare protein variants at low concentrations, as neutral alleles, in conditions of low mutation frequency in monomorphic loci takes place in the population. Comparison of frequencies of rare variants among healthy newborns and the children with congenital disorders revealed their identity (1.0.10(-3)), as compared to 1.05.10(-3)). Simplification of the method for scoring mutations judging by rare blood protein variants, which is necessary for monitoring for gene mutations in human populations, stimulates development of novel approaches.     

Frequency of rare electrophoretic protein variants in normal human beings and in congenital pathology

17 blood proteins of infants with rough and multiple congenital malformations (CM), prematurely born infants and sick newborns without developmental anomalies were studied electrophoretically in polyacrylamide and starch gels (62422 locus tests). The control included blood samples of healthy newborns from ordinary maternity hospitals (60234 locus tests). The frequency of rare protein variants in all the cases was higher in sick children than in healthy ones. The frequency of rare genes (corrected for electrophoretically "silent" alleles) was 2.16 X 10(-3) in infants with CM and 0.99 X 10(-4) in the control. Examination of parents of 11 congenitally malformed infants with rare protein variants showed that at least in 5 cases such variants were absent in the parents and might be attributed to "fresh" mutations. However, only 3 variants (1 for serum albumin and 2 for red cell esterase) represented rare heterozygotes with codominant expression. This corresponds to the frequency of 0.59 X 10(-3). In the total population of newborns the proportion of infants with CM was 0.02, which means that the population mutation rate is 1.18 X 10(-5) per gene per generation. The data obtained support the conclusion about strong pressure of stabilizing selection against de novo mutations which change electrophoretic mobility of the protein molecule. The reasons for discrepancy between our data and the recent results of Neel and Mohrenweiser (1984) are discussed.     

Population genetics study of differential fertility in humans (based on an example of habitual abortion). II. Distribution of genotypes and heterozygosity levels for a variety of biochemical markers

Hereditary variation of 21 blood proteins coded by 22 monomorphic and 9 polymorphic loci was compared in 171 couples and 120 women with repeated spontaneous abortions in anamnesis (experimental group) and 183 couples with normal fertility (control group). Significant elevation of frequency of rare protein electrophoretic variants, marked deviation in distribution of genotypes for polymorphic loci, alteration in observed heterozygosity level were not found in experimental group. Statistically significant decrease in average value of D index (D = (hobs-hexp)/hexp; hobs, hexp--observed and expected heterozygosity per locus) was found in women of experimental group.     

A revised indirect estimate of mutation rates in Amerindians

We have previously raised the possibility that the mutation rate resulting in rare electrophoretic variants is higher in tribal/tropical-dwelling/nonindustrialized societies than in civilized/temperate-dwelling/industrialized societies. Here, we report the results of examining 11 additional proteins for the occurrence of rare electrophoretic variants in 10 Amerindian tribes, for a total of 8,968 determinations and a total of 17,648 locus tests. When these data are combined with the results of all our previous similar studies of Amerindians, a total of 272,298 polypeptides, the products of 43 different loci, have been examined for the occurrence of rare electrophoretic variants. On the assumption that these variants are maintained by mutation pressure and are essentially neutral in their phenotypic effects, we have calculated by three different approaches that it requires an average mutation rate of 1.3 X 10(-5)/locus per generation to maintain the observed variant frequency. Concurrently, we are reporting elsewhere that a direct estimate of the mutation rate resulting in electromorphs in various studies of civilized industrialized populations is 0.3 X 10(-5)/locus per generation. Although this difference appears to have statistical significance, the nonquantifiable uncertainties to both approaches are such that our enthusiasm for a true difference in mutation rates between the two types of populations has diminished. However, even the lower of these estimates, when corrected for all the types of genetic variation that electrophoresis does not detect, implies total locus and gametic mutation rates well above those which in the past have dominated genetic thinking.     

Genetic variants of the lipoprotein lipase gene and myocardial infarction in the Central Valley of Costa Rica

To assess common variants of the LPL gene that could influence susceptibility to myocardial infarction (MI), we assessed three functional single-nucleotide polymorphisms (SNPs), D9N, N291S, and S447X, in 1,321 survivors of a first acute MI and 1,321 population-based controls, matched for age, gender, and area of residence, all living in the Central Valley of Costa Rica. Conditional logistic regression was used to estimate odds ratio (OR) and 95% confidence interval (CI). The frequency of the X447 mutant allele was significantly lower in cases than in controls (6.2% vs. 7.6%; P < 0.01), whereas no association with MI was found for D9N or N291S. The OR (95% CI) for carriers vs. noncarriers of the X447 allele was 0.80 (0.63-1.01); when considering the haplotype that contained X447 and normal alleles of D9N and N291S, the OR (95% CI) was 0.66 (0.48-0.91). Twelve other SNPs were assessed in a subgroup of the population, of which the four functional SNPs were found to be monomorphic, and no correlation with MI was observed for the other eight neutral SNPs. The X447 mutant allele of the LPL gene may protect from MI risk, although this effect is small.     

Construction and behavior of circularly permuted and telocentric chromosomes in Saccharomyces cerevisiae.

We developed techniques that allow us to construct novel variants of Saccharomyces cerevisiae chromosomes. These modified chromosomes have precisely determined structures. A metacentric derivative of chromosome III which lacks the telomere-associated X and Y' elements, which are found at the telomeres of most yeast chromosomes, behaves normally in both mitosis and meiosis. We made a circularly permuted telocentric version of yeast chromosome III whose closest telomere was 33 kilobases from the centromere. This telocentric chromosome was lost at a frequency of 1.6 X 10(-5) per cell compared with a frequency of 4.0 X 10(-6) for the natural metacentric version of chromosome III. An extremely telocentric chromosome whose closet telomere was only 3.5 kilobases from the centromere was lost at a frequency of 6.0 X 10(-5). The mitotic stability of telocentric chromosomes shows that the very high frequency of nondisjunction observed for short linear artificial chromosomes is not due to inadequate centromere-telomere separation.     

Distribution of single‐nucleotide variants on protein–protein interaction sites and its relationship with minor allele frequency

Recent advances in DNA sequencing techniques have identified rare single‐nucleotide variants with less than 1% minor allele frequency. Despite the growing interest and physiological importance of rare variants in genome sciences, less attention has been paid to the allele frequency of variants in protein sciences. To elucidate the characteristics of genetic variants on protein interaction sites, from the viewpoints of the allele frequency and the structural position of variants, we mapped about 20,000 human SNVs onto protein complexes. We found that variants are less abundant in protein interfaces, and specifically the core regions of interfaces. The tendency to “avoid” the interfacial core is stronger among common variants than rare variants. As amino acid substitutions, the trend of mutating amino acids among rare variants is consistent in different interfacial regions, reflecting the fact that rare variants result from random mutations in DNA sequences, whereas amino acid changes of common variants vary between the interfacial core and rim regions, possibly due to functional constraints on proteins. This study illustrated how the allele frequency of variants relates to the protein structural regions and the functional sites in general and will lead to deeper understanding of the potential deleteriousness of rare variants at the structural level. Exceptional cases of the observed trends will shed light on the limitations of structural approaches to evaluate the functional impacts of variants.     

Electroporation-induced transformation of intact cells of Clostridium perfringens.

Electroporation-induced transformation of intact cells of Clostridium perfringens 3624A with plasmids pAMB1 and pHR106 resulted in 3.8 X 10(-5) and 4.2 X 10(-4) transformants per viable cell, respectively. With respect to shuttle plasmid pHR106, these values represent a greater than 100-fold increase in transformation frequency when compared with the values reported with polyethylene glycol-induced L-phase variants.     

Electroporation-induced transformation of intact cells of Clostridium perfringens

Electroporation-induced transformation of intact cells of Clostridium perfringens 3624A with plasmids pAMB1 and pHR106 resulted in 3.8 X 10(-5) and 4.2 X 10(-4) transformants per viable cell, respectively. With respect to shuttle plasmid pHR106, these values represent a greater than 100-fold increase in transformation frequency when compared with the values reported with polyethylene glycol-induced L-phase variants.     

Characterization of ouabain-resistant mutants of a canine kidney cell line, MDCK

Madin-Darby canine kidney (MDCK) cells were mutagenized and variants resistant to 10, 160, and 2000 times the ouabain lethal dose for wild type cells selected. The phenotypes were stable in the absence of selection. The frequencies with which variants were recovered were consistent with genetic alterations being responsible for drug resistance. It was shown that 50% of the (Na+, K+)-ATPase activity present in mutant cells had a higher Kd for ouabain than normal while 50% remained wild type for ouabain binding. Wild type MDCK cells were measured to have 2 X 10(6) ouabain binding sites per cell with a Kd for the drug of 0.6-1.0 X 10(-7) M. The novel (Na+, K+)-ATPase activities in the mutants demonstrated Kd values for ouabain of 10(-5) M, 3 X 10(-4) M, or 3 X 10(-3) M for the different mutant classes tested. The rate of synthesis of the (Na+, K+)-ATPase as well as the total amount of enzyme per unit of cell protein was unaltered in the mutants. Comparison of the alpha subunit of the enzyme, known to contain the ouabain-binding site, by sodium dodecyl sulfate-gel electrophoresis did not reveal any difference in the size of this subunit in mutant versus wild type cells.     

Alternative splice variants of alpha 7 beta 1 integrin selectively recognize different laminin isoforms.

The integrin alpha(7)beta(1) occurs in several cytoplasmic (alpha(7A), alpha(7B)) and extracellular splice variants (alpha(7X1), alpha(7X2)), which are differentially expressed during development of skeletal and heart muscle. The extracellular variants result from the alternative splicing of exons X1 and X2, corresponding to a segment within the putative ligand binding domain. To study the specificity and affinity of the X1/X2 variants to different laminin isoforms, soluble alpha(7)beta(1) complexes were prepared by recombinant coexpression of the extracellular domains of the alpha- and beta-subunits. The binding of these complexes to purified ligands was measured by solid phase binding assays. Surprisingly, the alternative splice variants revealed different and specific affinities to different laminin isoforms. While the alpha(7X2) variant bound much more strongly to laminin-1 than the alpha(7X1) variant, the latter showed a high affinity binding to laminins-8 and -10/11. Laminin-2, the major laminin isoform in skeletal muscle, was recognized by both variants, whereas none of the two variants were able to interact with laminin-5. A specific blocking antibody inhibited the binding of both variants to all laminins tested, indicating the involvement of common epitopes in alpha(7X1)beta(1) and alpha(7X2)beta(1). Because laminin-8 and -10/11 as well as alpha(7X1) are expressed in developing skeletal and cardiac muscle, these findings suggest that alpha(7X1)beta(1) may represent a physiological receptor with novel specificities for laminin-8 and -10.     

NextGen sequencing reveals short double crossovers contribute disproportionately to genetic diversity in Toxoplasma gondii

Background

Toxoplasma gondii is a widespread protozoan parasite of animals that causes zoonotic disease in humans. Three clonal variants predominate in North America and Europe, while South American strains are genetically diverse, and undergo more frequent recombination. All three northern clonal variants share a monomorphic version of chromosome Ia (ChrIa), which is also found in unrelated, but successful southern lineages. Although this pattern could reflect a selective advantage, it might also arise from non-Mendelian segregation during meiosis. To understand the inheritance of ChrIa, we performed a genetic cross between the northern clonal type 2 ME49 strain and a divergent southern type 10 strain called VAND, which harbors a divergent ChrIa.

Results

NextGen sequencing of haploid F1 progeny was used to generate a genetic map revealing a low level of conventional recombination, with an unexpectedly high frequency of short, double crossovers. Notably, both the monomorphic and divergent versions of ChrIa were isolated with equal frequency. As well, ChrIa showed no evidence of being a sex chromosome, of harboring an inversion, or distorting patterns of segregation. Although VAND was unable to self fertilize in the cat, it underwent successful out-crossing with ME49 and hybrid survival was strongly associated with inheritance of ChrIII from ME49 and ChrIb from VAND.

Conclusions

Our findings suggest that the successful spread of the monomorphic ChrIa in the wild has not been driven by meiotic drive or related processes, but rather is due to a fitness advantage. As well, the high frequency of short double crossovers is expected to greatly increase genetic diversity among progeny from genetic crosses, thereby providing an unexpected and likely important source of diversity.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-1168) contains supplementary material, which is available to authorized users.     

The dose-response relationship for ethyl methanesulfonate-induced mutations at the hypoxanthine-guanine phosphoribosyl transferase locus in Chinese hamster ovary cells.

The frequency of ethyl methanesulfonate (EMS)-induced mutations to 6-thioguanine resistance in a Chinese hamster ovary cells clone K1-BH4 was studied at many EMS doses including the minimally lethal range (0-100 microng/ml) as well as the exponential killing portion (100-800 microng/ml) of the survival curve. The mutation frequency increases approximately in proportion with increasing EMS concentration at a fixed treatment time. The pooled data for the observed mutation frequency, f(X), as a function of EMS dose X, is adequately described by a linear function f(X)=10(-6)(8.73+3.45 X), where 0 less than or equal to X less than or equal to 800 microng/ml. One interpretation of the linear dose-response is that, as a result of EMS treatment, ethylation of cellular constituents occurs, which is directly responsible for the mutation. Biochemical analyses demonstrate that most of the randomly isolated 6-thioguanine-resistant variants possess a highly reduced or undetectable level of HGPRT activity suggesting that the EMS-induced mutations to 6-thioguanine resistance affect primarily, if not exclusively, the HGPRT locus.     

BRCA1 variants in a family study of African-American and Latina women

We sequenced the entire coding region of BRCA1 to improve our understanding of the frequency and nature of BRCA1 variants in African-American and Latina women identified from a multiethnic cohort in Los Angeles, California. The study included 109 African-American and 140 Latina sibships from families with two or more cases of breast or ovarian cancer among first-degree relatives. BRCA1 was sequenced in 278 breast or ovarian cancer cases and 229 unaffected sisters. The proportion of cases with known disease-causing mutations was low (0.72, 95% confidence interval: 0–1.7%). In total, 33 sequence variants were identified, including two protein truncation mutations, one deletion, and six silent and 24 missense variants. Two novel rare variants were identified that appeared to act as benign polymorphisms. Four rare variants may be unique to women of African descent based on existing literature, and three have been described exclusively in Latina women. The frequency of common variants was similar for cases and controls, but the frequency of common variants for African-American women significantly differed from those previously described for Caucasian women. We believe this to be the largest study of high-risk African-American and Latina women sequenced for variants in the BRCA1 gene to date.     

Leveraging Identity-by-Descent for Accurate Genotype Inference in Family Sequencing Data

Sequencing family DNA samples provides an attractive alternative to population based designs to identify rare variants associated with human disease due to the enrichment of causal variants in pedigrees. Previous studies showed that genotype calling accuracy can be improved by modeling family relatedness compared to standard calling algorithms. Current family-based variant calling methods use sequencing data on single variants and ignore the identity-by-descent (IBD) sharing along the genome. In this study we describe a new computational framework to accurately estimate the IBD sharing from the sequencing data, and to utilize the inferred IBD among family members to jointly call genotypes in pedigrees. Through simulations and application to real data, we showed that IBD can be reliably estimated across the genome, even at very low coverage (e.g. 2X), and genotype accuracy can be dramatically improved. Moreover, the improvement is more pronounced for variants with low frequencies, especially at low to intermediate coverage (e.g. 10X to 20X), making our approach effective in studying rare variants in cost-effective whole genome sequencing in pedigrees. We hope that our tool is useful to the research community for identifying rare variants for human disease through family-based sequencing.     

Alz 50, a monoclonal antibody to Alzheimer's disease antigen, cross-reacts with tau proteins from bovine and normal human brain

Microtubule-associated protein tau from bovine brain reacted on immunoblots and on enzyme-linked immunosorbent assay with a monoclonal antibody, Alz 50, which has previously been found to bind to an Alzheimer disease-specific antigen. The apparent affinity of binding of Alz 50 to tau was 2.1 X 10(-9) M on competitive enzyme-linked immunosorbent assay, and it was in the same range as for Tau-1 (0.5 X 10(-9) M), an antibody raised against purified bovine tau proteins. Immunoblotting of trypsin-digested tau revealed differences between Alz 50 and Tau-1 binding sites. The binding of both antibodies to tau was not affected by prior treatment with phosphatase, indicating that the cross-reactivity of Alz 50 with tau is due to the presence of phosphate-independent epitope. This epitope then differs from phosphate-dependent tau epitopes often shared with other cytoskeletal proteins. Alz 50 and Tau-1 binding sites were present in all isoelectric (pI 6-8) and molecular weight variants of tau. In contrast, phosphate-dependent epitopes recognized by another tau-reactive antibody (NP14) were found mostly in acidic tau variants. Similarly to tau proteins from bovine brain, tau-enriched preparations from normal human brain contained Alz 50 and Tau-1 reactive sites in all isoelectric (pI 6.5-8.5) and molecular weight variants. Our observation of Alz 50 cross-reactivity with tau suggests a relationship between tau and the novel protein identified recently in Alzheimer brains.     

Identification of ABCG2 dysfunction as a major factor contributing to gout

The ATP-binding cassette, subfamily G, member 2 gene ABCG2/BCRP locates in a gout-susceptibility locus (MIM 138900) on chromosome 4q. Recent genome-wide association studies also showed that the ABCG2 gene relates to serum uric acid levels and gout. Since ABCG2 is also known as a transporter of nucleotide analogs that are structurally similar to urate, and is an exporter that has common polymorphic reduced functionality variants, ABCG2 could be a urate secretion transporter and a gene causing gout. To find candidate mutations in ABCG2, we performed a mutation analysis of the ABCG2 gene in 90 Japanese patients with hyperuricemia and found six non-synonymous mutations. Among the variants, ATP-dependent urate transport was reduced or eliminated in five variants, and two out of the five variants (Q126X and Q141K) were frequently detected in patients. Haplotype frequency analysis revealed that there is no simultaneous presence of Q126X and Q141K in one haplotype. As Q126X and Q141K are a nonfunctional and half-functional haplotype, respectively, their genotype combinations are divided into four estimated functional groups. The association study with 161 male gout patients and 865 male controls showed that all of those who had dysfunctional ABCG2 had an increased risk of gout, and that a remarkable risk was observed in those with ≤1/4 function (OR, 25.8; 95% CI, 10.3-64.6; p = 3.39 × 10(-21)). In 2,150 Japanese individuals, the frequency of those with dysfunctional ABCG2 was more than 50%. Our function-based clinicogenetic analysis identified the combinations of dysfunctional variants of ABCG2 as a major contributing factor in Japanese patients with gout.     

Statistical Studies on Protein Polymorphism in Natural Populations. III. Distribution of Allele Frequencies and the Number of Alleles per Locus

With the aim of understanding the mechanism of maintenance of protein polymorphism, we have studied the properties of allele frequency distribution and the number of alleles per locus, using gene-frequency data from a wide range of organisms (mammals, birds, reptiles, amphibians, Drosophila and non-Drosophila invertebrates) in which 20 or more loci with at least 100 genes were sampled. The observed distribution of allele frequencies was U-shaped in all of the 138 populations (mostly species or subspecies) examined and generally agreed with the theoretical distribution expected under the mutation-drift hypothesis, though there was a significant excess of rare alleles (gene frequency, 0 ~ 0.05) in about a quarter of the populations. The agreement between the mutation-drift theory and observed data was quite satisfactory for the numbers of polymorphic (gene frequency, 0.05 ~ 0.95) and monomorphic (0.95 ~ 1.0) alleles.—The observed pattern of allele-frequency distribution was incompatible with the prediction from the overdominance hypothesis. The observed correlations of the numbers of rare alleles, polymorphic alleles and monomorphic alleles with heterozygosity were of the order of magnitude that was expected under the mutation-drift hypothesis. Our results did not support the view that intracistronic recombination is an important source of genetic variation. The total number of alleles per locus was positively correlated with molecular weight in most of the species examined, and the magnitude of the correlation was consistent with the theoretical prediction from mutation-drift hypothesis. The correlation between molecular weight and the number of alleles was generally higher than the correlation between molecular weight and heterozygosity, as expected.     

The structure of the amino terminus of tropomyosin is critical for binding to actin in the absence and presence of troponin

Analysis of two recombinant variants of chicken striated muscle alpha-tropomyosin has shown that the structure of the amino terminus is crucial for most aspects of tropomyosin function: affinity to actin, promotion of binding to actin by troponin, and regulation of the actomyosin MgATPase. Initial characterization of variants expressed and isolated from Escherichia coli has been published (Hitchcock-DeGregori, S. E., and Heald, R. W. (1987) J. Biol. Chem. 262, 9730-9735). Fusion tropomyosin contains 80 amino acids of a nonstructural influenza virus protein (NS1) on the amino terminus. Nonfusion tropomyosin is a variant because the amino-terminal methionine is not acetylated (unacetylated tropomyosin). The affinity of tropomyosin labeled at Cys190 with N-[14C]ethylmaleimide for actin was measured by cosedimentation in a Beckman Airfuge. Fusion tropomyosin binds to actin with an affinity slightly greater than that of chicken striated muscle alpha-tropomyosin (Kapp = 1-2 X 10(7) versus 0.5-1 X 10(7) M-1) and more strongly than unacetylated tropomyosin (Kapp = 3 X 10(5) M-1). Both variants bind cooperatively to actin. Troponin increases the affinity of unacetylated tropomyosin for actin (+Ca2+, Kapp = 6 X 10(6) M-1; +EGTA, Kapp = 2 X 10(7) M-1), but the affinity is still lower than that of muscle tropomyosin for actin in the presence of troponin (Kapp much greater than 10(8) M-1). Troponin has no effect on the affinity of fusion tropomyosin for actin indicating that binding of troponin T to the over-lap region of the adjacent tropomyosin, presumably sterically prevented by the fusion peptide in fusion tropomyosin, is required for troponin to promote the binding of tropomyosin to actin. The role of troponin T in regulation and the mechanisms of cooperative binding of tropomyosin to actin have been discussed in relation to this work.     

Genetic studies of blood markers in Przewalski's horses

Ninety-six Przewalski's horses (Equus przewalskii) were blood typed using systems of inherited blood variants known to be highly effective for parentage testing of domestic horses (E. caballus). Sixteen red cell antigenic factors detected using sera prepared by alloimmunization of domestic horses were shown to be inherited in six systems (A, C, D, P, Q, and U) and in the same patterns as domestic horses. Family data confirmed autosomal, codominant inheritance at five loci of serum protein variants (Al, Tf, Xk, Pi, and Es) and three loci of red cell proteins (PGM, PHI, and Hb). One serum protein locus (Gc) and two red cell protein loci (PGD and CA) appeared to be monomorphic. Despite the narrow genetic base and high inbreeding coefficients of captive Przewalski's horses, average heterozygosity calculated over 18 loci was estimated to be 0.320 +/- 0.05, which was similar to that found in five breeds of domestic horses.