Genetic polymorphisms of manganese-superoxide dismutase and glutathione-S-transferase in chronic alcoholic pancreatitis

Chronic alcohol consumption is a major risk factor for the development of chronic pancreatitis. However, chronic pancreatitis occurs only in a minority of heavy drinkers. This variability may be due to yet unidentified genetic factors. Several enzymes involved in the degradation of reactive oxidants and xenobiotics, such as glutathione-S-transferase P1 (GSTP1) and manganese-superoxide dismutase (MnSOD) reveal functional polymorphisms that affect the antioxidative capacity and may therefore modulate the development of chronic pancreatitis and long-term complications like endocrine and exocrine pancreatic insufficiency. Two functional polymorphisms of the MnSOD and the GSTP1 gene were assessed by polymerase chain reaction and restriction fragment length polymorphism in 165 patients with chronic alcoholic pancreatitis, 140 alcoholics without evidence of pancreatic disease and 160 healthy control subjects. The distribution of GSTP1 and MnSOD genotypes were in Hardy-Weinberg equilibrium in the total cohort. Genotype and allele frequencies for both genes were not statistically different between the three groups. Although genotype MnSOD Ala/Val was seemingly associated with the presence of exocrine pancreatic insufficiency, this subgroup was too small and the association statistically underpowered. None of the tested genotypes affected the development of endocrine pancreatic insufficiency. Polymorphisms of MnSOD and GSTP1 are not associated with chronic alcoholic pancreatitis. The present data emphasize the need for stringently designed candidate gene association studies with well-characterized cases and controls and sufficient statistical power to exclude chance observations.

Study of a Swiss dopa-responsive dystonia family with a deletion in GCH1: redefining DYT14 as DYT5

OBJECTIVE: To report the study of a multigenerational Swiss family with dopa-responsive dystonia (DRD). METHODS: Clinical investigation was made of available family members, including historical and chart reviews. Subject examinations were video recorded. Genetic analysis included a genome-wide linkage study with microsatellite markers (STR), GTP cyclohydrolase I (GCH1) gene sequencing, and dosage analysis. RESULTS: We evaluated 32 individuals, of whom 6 were clinically diagnosed with DRD, with childhood-onset progressive foot dystonia, later generalizing, followed by parkinsonism in the two older patients. The response to levodopa was very good. Two additional patients had late onset dopa-responsive parkinsonism. Three other subjects had DRD symptoms on historical grounds. We found suggestive linkage to the previously reported DYT14 locus, which excluded GCH1. However, further study with more stringent criteria for disease status attribution showed linkage to a larger region, which included GCH1. No mutation was found in GCH1 by gene sequencing but dosage methods identified a novel heterozygous deletion of exons 3 to 6 of GCH1. The mutation was found in seven subjects. One of the patients with dystonia represented a phenocopy. CONCLUSIONS: This study rules out the previously reported DYT14 locus as a cause of disease, as a novel multiexonic deletion was identified in GCH1. This work highlights the necessity of an accurate clinical diagnosis in linkage studies as well as the need for appropriate allele frequencies, penetrance, and phenocopy estimates. Comprehensive sequencing and dosage analysis of known genes is recommended prior to genome-wide linkage analysis.

No role of matrixmetalloproteinase-3 genetic promoter polymorphism 1171 as a risk factor for cirrhosis in alcoholic liver disease

BACKGROUND: As only a minority of alcoholics develop cirrhosis, polymorphic genes, whose products are involved in fibrosis development were suggested to confer individual susceptibility. We tested whether a functional promoter polymorphism in the gene encoding matrix metalloproteinase-3 (MMP-3; 1171 5A/6A) was associated liver cirrhosis in alcoholics. METHODS: Independent cohorts from the UK and Germany were studied. (i) UK cohort: 320 alcoholic cirrhotics and 183 heavy drinkers without liver damage and (ii) German cohort: 149 alcoholic cirrhotics, 220 alcoholic cirrhotics who underwent liver transplantation and 151 alcoholics without liver disease. Patients were genotyped for MMP-3 variants by restriction fragment length polymorphism, single strand confirmation polymorphism, and direct sequencing. In addition, MMP-3 transcript levels were correlated with MMP-3 genotype in normal liver tissues. RESULTS: Matrix metalloproteinase-3 genotype and allele distribution in all 1023 alcoholic patients were in Hardy-Weinberg equilibrium. No significant differences in MMP-3 genotype and allele frequencies were observed either between alcoholics with or without cirrhosis. There were no differences in hepatic mRNA transcription levels according to MMP-3 genotype. CONCLUSIONS: Matrix metalloproteinase-3 1171 promoter polymorphism plays no role in the genetic predisposition for liver cirrhosis in alcoholics. Stringently designed candidate gene association studies are required to exclude chance observations.

Simulating the evolution of a clonal trait in plants with sexual and vegetative reproduction

Aims Phenotypic optimality models neglect genetics. However, especially when heterozygous genotypes ire fittest, evolving allele, genotype and phenotype frequencies may not correspond to predicted optima. This was not previously addressed for organisms with complex life histories. Methods Therefore, we modelled the evolution of a fitness-relevant trait of clonal plants, stolon internode length. We explored the likely case of air asymmetric unimodal fitness profile with three model types. In constant selection models (CSMs), which are gametic, but not spatially explicit, evolving allele frequencies in the one-locus and five-loci cases did not correspond to optimum stolon internode length predicted by the spatially explicit, but not gametic, phenotypic model. This deviation was due to the asymmetry of the fitness profile. Gametic, spatially explicit individual-based (SEIB) modeling allowed us relaxing the CSM assumptions of constant selection with exclusively sexual reproduction. Important findings For entirely vegetative or sexual reproduction, predictions. of the gametic SEIB model were close to the ones of spatially explicit CSMs gametic phenotypic models, hut for mixed modes of reproduction they appoximated those of gametic, not spatially explicit CSMs. Thus, in contrast to gametic SEIB models, phenotypic models and, especially for few loci, also CSMs can be very misleading. We conclude that the evolution of trails governed by few quantitative trait loci appears hardly predictable by simple models, that genetic algorithms aiming at technical optimization may actually, miss the optimum and that selection may lead to loci with smaller effects, in derived compared with ancestral lines.

A MIF haplotype is associated with the outcome of patients with severe sepsis: a case control study

BACKGROUND: Macrophage migration inhibitory factor (MIF) plays an important regulatory role in sepsis. In the promoter region a C/G single nucleotide polymorphism (SNP) at position -173 (rs755622) and a CATT5-8 microsatellite at position -794 are related to modified promoter activity. The purpose of the study was to analyze their association with the incidence and outcome of severe sepsis. METHODS: Genotype distributions and allele frequencies in 169 patients with severe sepsis, 94 healthy blood donors and 183 postoperative patients without signs of infection or inflammation were analyzed by real time PCR and Sequence analysis. All included individuals were Caucasians. RESULTS: Genotype distribution and allele frequencies of severe sepsis patients were comparable to both control groups. However, the genotype and allele frequencies of both polymorphisms were associated significantly with the outcome of severe sepsis. The highest risk of dying from severe sepsis was detectable in patients carrying a haplotype with the alleles -173 C and CATT7 (p = 0.0005, fisher exact test, RR = 1,806, CI: 1.337 to 2.439). CONCLUSION: The haplotype with the combination of the -173 C allele and the -794 CATT7 allele may not serve as a marker for susceptibility to sepsis, but may help identify septic patients at risk of dying.

Quantitative sequence analysis of FBN1 premature termination codons provides evidence for incomplete NMD in leukocytes

We improved, evaluated, and used Sanger sequencing for quantification of single nucleotide polymorphism (SNP) variants in transcripts and gDNA samples. This improved assay resulted in highly reproducible relative allele frequencies (e.g., for a heterozygous gDNA 50.0+/-1.4%, and for a missense mutation-bearing transcript 46.9+/-3.7%) with a lower detection limit of 3-9%. It provided excellent accuracy and linear correlation between expected and observed relative allele frequencies. This sequencing assay, which can also be used for the quantification of copy number variations (CNVs), methylations, mosaicisms, and DNA pools, enabled us to analyze transcripts of the FBN1 gene in fibroblasts and blood samples of patients with suspected Marfan syndrome not only qualitatively but also quantitatively. We report a total of 18 novel and 19 known FBN1 sequence variants leading to a premature termination codon (PTC), 26 of which we analyzed by quantitative sequencing both at gDNA and cDNA levels. The relative amounts of PTC-containing FBN1 transcripts in fresh and PAXgene-stabilized blood samples were significantly higher (33.0+/-3.9% to 80.0+/-7.2%) than those detected in affected fibroblasts with inhibition of nonsense-mediated mRNA decay (NMD) (11.0+/-2.1% to 25.0+/-1.8%), whereas in fibroblasts without NMD inhibition no mutant alleles could be detected. These results provide evidence for incomplete NMD in leukocytes and have particular importance for RNA-based analyses not only in FBN1 but also in other genes.

Sequencing of TGF-beta pathway genes in familial cases of intracranial aneurysm.

BACKGROUND AND PURPOSE: Familial aggregation of intracranial aneurysms (IA) strongly suggests a genetic contribution to pathogenesis. However, genetic risk factors have yet to be defined. For families affected by aortic aneurysms, specific gene variants have been identified, many affecting the receptors to transforming growth factor-beta (TGF-beta). In recent work, we found that aortic and intracranial aneurysms may share a common genetic basis in some families. We hypothesized, therefore, that mutations in TGF-beta receptors might also play a role in IA pathogenesis. METHODS: To identify genetic variants in TGF-beta and its receptors, TGFB1, TGFBR1, TGFBR2, ACVR1, TGFBR3, and ENG were directly sequenced in 44 unrelated patients with familial IA. Novel variants were confirmed by restriction digestion analyses, and allele frequencies were analyzed in cases versus individuals without known intracranial disease. Similarly, allele frequencies of a subset of known SNPs in each gene were also analyzed for association with IA. RESULTS: No mutations were found in TGFB1, TGFBR1, TGFBR2, or ACVR1. Novel variants identified in ENG (p.A60E) and TGFBR3 (p.W112R) were not detected in at least 892 reference chromosomes. ENG p.A60E showed significant association with familial IA in case-control studies (P=0.0080). No association with IA could be found for any of the known polymorphisms tested. CONCLUSIONS: Mutations in TGF-beta receptor genes are not a major cause of IA. However, we identified rare variants in ENG and TGFBR3 that may be important for IA pathogenesis in a subset of families.

Gene flow and the maintenance of species boundaries

Hybrid zones are regions where individuals from genetically differentiated populations meet and mate, resulting in at least some offspring of mixed ancestry. Patterns of gene flow (introgression) in hybrid zones vary across the genome, allowing assessment of the role of individual genes or genome regions in reproductive isolation. Here, we document patterns of introgression between two recently diverged species of field crickets. We sampled at a very fine spatial scale and genotyped crickets for 110 highly differentiated single nucleotide polymorphisms (SNPs) identified through transcriptome scans. Using both genomic and geographic cline analysis, we document remarkably abrupt transitions (<100 m) in allele frequencies for 50 loci, despite high levels of gene flow at other loci. These are among the steepest clines documented for any hybridizing taxa. Furthermore, the cricket hybrid zone provides one of the clearest examples of the semi-permeability of species boundaries. Comparisons between data from the fine-scale transect and data (for the same set of markers) from sampling a much larger area in a different region of the cricket hybrid zone reveal consistent patterns of introgression for individual loci. The consistency in patterns of introgression between these two distant and distinct regions of the hybrid zone suggests that strong selection is acting to maintain abrupt discontinuities within the hybrid zone and that genomic regions with restricted introgression likely include genes that contribute to nonecological prezygotic barriers.

Distance from Sub-Saharan Africa Predicts Mutational Load in Diverse Human Genomes

The Out-of-Africa (OOA) dispersal ∼50,000 y ago is characterized by a series of founder events as modern humans expanded into multiple continents. Population genetics theory predicts an increase of mutational load in populations undergoing serial founder effects during range expansions. To test this hypothesis, we have sequenced full genomes and high-coverage exomes from seven geographically divergent human populations from Namibia, Congo, Algeria, Pakistan, Cambodia, Siberia, and Mexico. We find that individual genomes vary modestly in the overall number of predicted deleterious alleles. We show via spatially explicit simulations that the observed distribution of deleterious allele frequencies is consistent with the OOA dispersal, particularly under a model where deleterious mutations are recessive. We conclude that there is a strong signal of purifying selection at conserved genomic positions within Africa, but that many predicted deleterious mutations have evolved as if they were neutral during the expansion out of Africa. Under a model where selection is inversely related to dominance, we show that OOA populations are likely to have a higher mutation load due to increased allele frequencies of nearly neutral variants that are recessive or partially recessive.

Implications of Longterm Diameter-Limit Harvesting: Effects on Radial Growth of Red Spruce (Picea rubens) and Genetic Diversity of White Pine (Pinus strobus)

For over 3 centuries, diameter-limit harvesting has been a predominant logging method in the northeastern United States. Silvicultural theory asserts that such intensively selective harvesting can lead to genetic degradation. A decrease in softwood productivity has recently been reported in Maine - has a long history of dysgenic selection degraded the genetic resources of Maine softwoods, contributing to a decrease in growth and productivity? This study examines two aspects of potential implications of diameter-limit harvesting: effects on residual phenotypes of red spruce and impacts on genetic diversity of white pine. Radial growth of residual red spruce trees in stands experiencing 50 years of fixed diameter-limit harvesting was measured using annual increment rings and compared with residual red spruce trees in positive selection stands. Trees remaiaing after several rounds of diameter-limit harvesting exhibited sigdicantl y smaller radial sizes throughout their lives, and displayed significantly slower growth rates for the first 80 years of measured growth. These results strongly suggest that the largest and fastest-growing genotypes and their respective gene complexes determining good radial growth have been removed from the diameter-limit stand. Dysgenic selection can be observed in fixed diarneter-limit stands, resulting in a diminished genetic resource and decreased residual stand value. To examine more direct genetic implications of long-term diameter-limit harvesting, microsatellite DNA markers were implemented to study genetic diversity of eastern white pine in Maine. Three age groups of trees were studied: mature trees older than 200 years, juvenile trees 5-30 years old, and embryos. Trees were genotyped at 10 microsatellite loci. Overall genetic diversity levels of eastern white pine in Maine were extremely high, with an average observed heterozygosity of 0.762. Genetic differentiation was minimal among and between all three age groups, although an excess of heterozygotes was shown in the mature and juvenile groups that was not reflected in the embryo group, which actually had a slight heterozygote deficiency. Allele frequencies did not differ significantly between age groups, but did reveal more rare and low frequency alleles in the embryo groups than in the mature group. Overall, low frequency alleles comprise the largest portion of alleles in the sample population, with no common alleles evident overall. These results suggest that significant genetic degradation has either not occurred for white pine, or that the results of dysgenic selection have not yet emerged. It is clear, however, that selective harvesting could result in a loss of low frequency alleles, which are a primary reserve of evolutionary potential in a species. Implications of these studies affect industrial forestry, regional economics, and ecological concerns for the northeast. Long-term diameter-limit harvesting can lead to a degradation of residual phenotypes, and an overall decrease in stand quality. Potentially, a loss of low frequency, locally adapted alleles could result in a decrease of allelic richness and degradation of the regidnal genetic resource. Decreased genetic variation can lead to seriously limited evolutionary potential of species and ecosystems, particularly in rapidly changing environments. Based on these findings, I recommend a reassessment of any harvesting prescription that includes fixed diameter-limit removals, particularly for species that have low natural genetic diversity levels or a limited natural range, such as red spruce. Maintenance of a healthy genetic reserve can avoid effects of dysgenic harvesting.

Exact moment calculations for genetic models with migration, mutation, and drift

Using properties of moment stationarity we develop exact expressions for the mean and covariance of allele frequencies at a single locus for a set of populations subject to drift, mutation, and migration. Some general results can be obtained even for arbitrary mutation and migration matrices, for example: (1) Under quite general conditions, the mean vector depends only on mutation rates, not on migration rates or the number of populations. (2) Allele frequencies covary among all pairs of populations connected by migration. As a result, the drift, mutation, migration process is not ergodic when any finite number of populations is exchanging genes. in addition, we provide closed form expressions for the mean and covariance of allele frequencies in Wright's finite-island model of migration under several simple models of mutation, and we show that the correlation in allele frequencies among populations can be very large for realistic rates of mutation unless an enormous number of populations are exchanging genes. As a result, the traditional diffusion approximation provides a poor approximation of the stationary distribution of allele frequencies among populations. Finally, we discuss some implications of our results for measures of population structure based on Wright's F-statistics.

Bayesian models for the analysis of genetic structure when populations are correlated

Motivation: Population allele frequencies are correlated when populations have a shared history or when they exchange genes. Unfortunately, most models for allele frequency and inference about population structure ignore this correlation. Recent analytical results show that among populations, correlations can be very high, which could affect estimates of population genetic structure. In this study, we propose a mixture beta model to characterize the allele frequency distribution among populations. This formulation incorporates the correlation among populations as well as extending the model to data with different clusters of populations. Results: Using simulated data, we show that in general, the mixture model provides a good approximation of the among-population allele frequency distribution and a good estimate of correlation among populations. Results from fitting the mixture model to a dataset of genotypes at 377 autosomal microsatellite loci from human populations indicate high correlation among populations, which may not be appropriate to neglect. Traditional measures of population structure tend to over-estimate the amount of genetic differentiation when correlation is neglected. Inference is performed in a Bayesian framework.

Genetic studies of the human plasma orosomucoid (alpha-1 acid glycoprotein)

Orosomucoid (ORM) or alpha-1 acid glycoprotein is an acute phase protein of human plasma whose function is suggested to be the competitive inhibition of cellular recognition by infective agents. Isoelectric focusing (IEF) and immunoblotting have been combined and optimum conditions have been determined for reliable classification of different ORM phenotypes. Addition of 6 M urea in an IEF gel revealed additional microheterogeneity in the ORM system which has not been previously reported. 1,667 individuals from different native ethnic groups of North and South America, Africa and New Guinea have been screened to determine the distribution of ORM alleles. Two common alleles, ORM1*1 and ORM1*2 have been observed and their frequencies were determined. Genetically independent variation consistent with expression of the ORM2 locus was observed in American and African blacks but was not observed in other sampled populations. The population allele frequencies for this new locus were 0.958, 0.025, 0.006, 0.011, for alleles ORM2*1, ORM2*2, ORM2*3, ORM2*4, respectively. Family studies confirm the autosomal codominant inheritance of the phenotypes observed at both ORM loci. ^

THE NATURAL AND ORTHOGONAL INTERACTION (NOIA) MODELS FOR QUANTITATIVE TRAITS (QTs) AND COMPLEX DISEASES

My dissertation focuses on developing methods for gene-gene/environment interactions and imprinting effect detections for human complex diseases and quantitative traits. It includes three sections: (1) generalizing the Natural and Orthogonal interaction (NOIA) model for the coding technique originally developed for gene-gene (GxG) interaction and also to reduced models; (2) developing a novel statistical approach that allows for modeling gene-environment (GxE) interactions influencing disease risk, and (3) developing a statistical approach for modeling genetic variants displaying parent-of-origin effects (POEs), such as imprinting. In the past decade, genetic researchers have identified a large number of causal variants for human genetic diseases and traits by single-locus analysis, and interaction has now become a hot topic in the effort to search for the complex network between multiple genes or environmental exposures contributing to the outcome. Epistasis, also known as gene-gene interaction is the departure from additive genetic effects from several genes to a trait, which means that the same alleles of one gene could display different genetic effects under different genetic backgrounds. In this study, we propose to implement the NOIA model for association studies along with interaction for human complex traits and diseases. We compare the performance of the new statistical models we developed and the usual functional model by both simulation study and real data analysis. Both simulation and real data analysis revealed higher power of the NOIA GxG interaction model for detecting both main genetic effects and interaction effects. Through application on a melanoma dataset, we confirmed the previously identified significant regions for melanoma risk at 15q13.1, 16q24.3 and 9p21.3. We also identified potential interactions with these significant regions that contribute to melanoma risk. Based on the NOIA model, we developed a novel statistical approach that allows us to model effects from a genetic factor and binary environmental exposure that are jointly influencing disease risk. Both simulation and real data analyses revealed higher power of the NOIA model for detecting both main genetic effects and interaction effects for both quantitative and binary traits. We also found that estimates of the parameters from logistic regression for binary traits are no longer statistically uncorrelated under the alternative model when there is an association. Applying our novel approach to a lung cancer dataset, we confirmed four SNPs in 5p15 and 15q25 region to be significantly associated with lung cancer risk in Caucasians population: rs2736100, rs402710, rs16969968 and rs8034191. We also validated that rs16969968 and rs8034191 in 15q25 region are significantly interacting with smoking in Caucasian population. Our approach identified the potential interactions of SNP rs2256543 in 6p21 with smoking on contributing to lung cancer risk. Genetic imprinting is the most well-known cause for parent-of-origin effect (POE) whereby a gene is differentially expressed depending on the parental origin of the same alleles. Genetic imprinting affects several human disorders, including diabetes, breast cancer, alcoholism, and obesity. This phenomenon has been shown to be important for normal embryonic development in mammals. Traditional association approaches ignore this important genetic phenomenon. In this study, we propose a NOIA framework for a single locus association study that estimates both main allelic effects and POEs. We develop statistical (Stat-POE) and functional (Func-POE) models, and demonstrate conditions for orthogonality of the Stat-POE model. We conducted simulations for both quantitative and qualitative traits to evaluate the performance of the statistical and functional models with different levels of POEs. Our results showed that the newly proposed Stat-POE model, which ensures orthogonality of variance components if Hardy-Weinberg Equilibrium (HWE) or equal minor and major allele frequencies is satisfied, had greater power for detecting the main allelic additive effect than a Func-POE model, which codes according to allelic substitutions, for both quantitative and qualitative traits. The power for detecting the POE was the same for the Stat-POE and Func-POE models under HWE for quantitative traits.

Distancias genéticas entre perfiles moleculares obtenidos desde marcadores multilocus multialélicos

Para expresar la magnitud de la identidad genética (similaridad) o su complemento (distancia) entre dos individuos caracterizados molecularmente a través de marcadores del tipo microsatélites (SSR), que son multilocusmultialélicos, es necesario elegir una métrica acorde con la naturaleza multivariada de los datos. Comúnmente, las métricas de distancias genéticas son diseñadas para expresar, en un único número, la diferencia genética entre dos poblaciones y son expresadas como función de la frecuencia alélica poblacional. Dichas métricas pueden también ser utilizadas para calcular la distancia entre perfiles individuales, pero las frecuencias alélicas no son continuas en este caso. Alternativamente, se pueden usar distancias geométricas obtenidas como el complemento del índice de similaridad para datos binarios que indican la presencia/ ausencia de cada alelo en un individuo. El objetivo de este trabajo fue evaluar simultáneamente el desempeño de ambos tipos de métricas para ordenar y clasificar individuos en una base de datos generadas a partir de loci de marcadores microsatélites SSR. Se calcularon 11 métricas de distancias a partir de 17 loci SSR obtenidos desde 17 introducciones de un banco de germoplasma de soja [Glycine max (L.) Merr.]. Se evaluó el consenso de los resultados obtenidos para la clasificación de los 17 perfiles moleculares desde varias métricas. Los resultados sugieren que los diferentes tipos de métricas producen información similar para comparar individuos. No obstante, se realizó una clasificación de las métricas que responden a diferencias entre los núcleos de las expresiones de cálculo.