Decay of linkage disequilibrium within genes across HGDP-CEPH human samples: most population isolates do not show increased LD

Background It is well known that the pattern of linkage disequilibrium varies between human populations, with remarkable geographical stratification. Indirect association studies routinely exploit linkage disequilibrium around genes, particularly in isolated populations where it is assumed to be higher. Here, we explore both the amount and the decay of linkage disequilibrium with physical distance along 211 gene regions, most of them related to complex diseases, across 39 HGDP-CEPH population samples, focusing particularly on the populations defined as isolates. Within each gene region and population we use r2 between all possible single nucleotide polymorphism (SNP) pairs as a measure of linkage disequilibrium and focus on the proportion of SNP pairs with r2 greater than 0.8. Results Although the average r2 was found to be significantly different both between and within continental regions, a much higher proportion of r2 variance could be attributed to differences between continental regions (2.8% vs. 0.5%, respectively). Similarly, while the proportion of SNP pairs with r2 > 0.8 was significantly different across continents for all distance classes, it was generally much more homogenous within continents, except in the case of Africa and the Americas. The only isolated populations with consistently higher LD in all distance classes with respect to their continent are the Kalash (Central South Asia) and the Surui (America). Moreover, isolated populations showed only slightly higher proportions of SNP pairs with r2 > 0.8 per gene region than non-isolated populations in the same continent. Thus, the number of SNPs in isolated populations that need to be genotyped may be only slightly less than in non-isolates. Conclusion The 'isolated population' label by itself does not guarantee a greater genotyping efficiency in association studies, and properties other than increased linkage disequilibrium may make these populations interesting in genetic epidemiology.

Rarity of DNA sequence alterations in the promoter region of the human androgen receptor gene

The human androgen receptor (AR) gene promoter lies in a GC-rich region containing two principal sites of transcription initiation and a putative Sp1 protein-binding site, without typical "TATA" and "CAAT" boxes. It has been suggested that mutations within the 5'untranslated region (5'UTR) may contribute to the development of prostate cancer by changing the rates of gene transcription and/or translation. In order to investigate this question, the aim of the present study was to search for the presence of mutations or polymorphisms at the AR-5'UTR in 92 prostate cancer patients, where histological diagnosis of adenocarcinoma was established in specimens obtained from transurethral resection or after prostatectomy. The AR-5'UTR was amplified by PCR from genomic DNA samples of the patients and of 100 healthy male blood donors, included as controls. Conformation-sensitive gel electrophoresis was used for DNA sequence alteration screening. Only one band shift was detected in one individual from the blood donor group. Sequencing revealed a new single nucleotide deletion (T) in the most conserved portion of the promoter region at position +36 downstream from the transcription initiation site I. Although the effect of this specific mutation remains unknown, its rarity reveals the high degree of sequence conservation of the human androgen promoter region. Moreover, the absence of detectable variation within the critical 5'UTR in prostate cancer patients indicates a low probability of its involvement in prostate cancer etiology.

High baseline serum total and LDL cholesterol levels are associated with MDR1 haplotypes in Brazilian hypercholesterolemic individuals of European descent

The MDR1 gene encodes the P-glycoprotein, an efflux transporter with broad substrate specificity. P-glycoprotein has raised great interest in pharmacogenetics because it transports a variety of structurally divergent drugs, including lipid-lowering drugs. The synonymous single-nucleotide polymorphism C3435T and the nonsynonymous single-nucleotide polymorphism G2677T/A in MDR1 have been indicated as potential determinants of variability in drug disposition and efficacy. In order to evaluate the effect of G2677T/A and C3435T MDR1 polymorphisms on serum levels of lipids before and after atorvastatin administration, 69 unrelated hypercholesterolemic individuals from São Paulo city, Brazil, were selected and treated with 10 mg atorvastatin orally once daily for four weeks. MDR1 polymorphisms were analyzed by PCR-RFLP. C3435T and G2677T polymorphisms were found to be linked. The allelic frequencies for C3435T polymorphism were 0.536 and 0.464 for the 3435C and 3435T alleles, respectively, while for G2677T/A polymorphism allele frequencies were 0.580 for the 2677G allele, 0.384 for the 2677T allele and 0.036 for the 2677A allele. There was no significant relation between atorvastatin response and MDR1 polymorphisms (repeated measures ANOVA; P > 0.05). However, haplotype analysis revealed an association between T/T carriers and higher basal serum total (TC) and LDL cholesterol levels (TC: 303 ± 56, LDL-C: 216 ± 57 mg/dl, respectively) compared with non-T/T carriers (TC: 278 ± 28, LDL-C: 189 ± 24 mg/dl; repeated measures ANOVA/Tukey test; P < 0.05). These data indicate that MDR1 polymorphism may have an important contribution to the control of basal serum cholesterol levels in Brazilian hypercholesterolemic individuals of European descent.

The search for circadian clock components in humans: new perspectives for association studies

Individual circadian clocks entrain differently to environmental cycles (zeitgebers, e.g., light and darkness), earlier or later within the day, leading to different chronotypes. In human populations, the distribution of chronotypes forms a bell-shaped curve, with the extreme early and late types _ larks and owls, respectively _ at its ends. Human chronotype, which can be assessed by the timing of an individual's sleep-wake cycle, is partly influenced by genetic factors - known from animal experimentation. Here, we review population genetic studies which have used a questionnaire probing individual daily timing preference for associations with polymorphisms in clock genes. We discuss their inherent limitations and suggest an alternative approach combining a short questionnaire (Munich ChronoType Questionnaire, MCTQ), which assesses chronotype in a quantitative manner, with a genome-wide analysis (GWA). The advantages of these methods in comparison to assessing time-of-day preferences and single nucleotide polymorphism genotyping are discussed. In the future, global studies of chronotype using the MCTQ and GWA may also contribute to understanding the influence of seasons, latitude (e.g., different photoperiods), and climate on allele frequencies and chronotype distribution in different populations.

Simultaneous detection of the C282Y, H63D and S65C mutations in the hemochromatosis gene using quenched-FRET real-time PCR

Hereditary hemochromatosis (HH) is a common autosomal disorder of iron metabolism mainly affecting Caucasian populations. Three recurrent disease-associated mutations have been detected in the hemochromatosis gene (HFE): C282Y, H63D, and S65C. Although HH phenotype has been associated with all three mutations, C282Y is considered the most relevant mutation responsible for hemochromatosis. Clinical complications of HH include cirrhosis of the liver, congestive cardiac failure and cardiac arrhythmias, endocrine pancreatic disease, which can be prevented by early diagnosis and treatment. Therefore, a reliable genotyping method is required for presymptomatic diagnosis. We describe the simultaneous detection of the C282Y, H63D and S65C mutations in the hemochromatosis gene by real-time PCR followed by melting curve analysis using fluorescence resonance energy transfer (FRET) probes. The acceptor fluorophore may be replaced by a quencher, increasing multiplex possibilities. Real-time PCR results were compared to the results of sequencing and conventional PCR followed by restriction digestion and detection by agarose gel electrophoresis (PCR-RFLP). Genotypes from 80 individuals obtained both by the conventional PCR-RFLP method and quenched-FRET real-time PCR were in full agreement. Sequencing also confirmed the results obtained by the new method, which proved to be an accurate, rapid and cost-effective diagnostic assay. Our findings demonstrate the usefulness of real-time PCR for the simultaneous detection of mutations in the HFE gene, which allows a reduction of a significant amount of time in sample processing compared to the PCR-RFLP method, eliminates the use of toxic reagents, reduces the risk of contamination in the laboratory, and enables full process automation.

Do Caucasian and Asian clocks tick differently?

The Period 3 and Clock genes are important components of the mammalian molecular circadian system. Studies have shown association between polymorphisms in these clock genes and circadian phenotypes in different populations. Nevertheless, differences in the pattern of allele frequency and genotyping distribution are systematically observed in studies with different ethnic groups. To investigate and compare the pattern of distribution in a sample of Asian and Caucasian populations living in Brazil, we evaluated two well-studied polymorphisms in the clock genes: a variable number of tandem repeats (VNTR) in PER3 and a single nucleotide polymorphism (SNP) in CLOCK. The aim of this investigation was to search for clues about human evolutionary processes related to circadian rhythms. We selected 109 Asian and 135 Caucasian descendants. The frequencies of the shorter allele (4 repeats) in the PER3 gene and the T allele in the CLOCK gene among Asians (0.86 and 0.84, respectively) were significantly higher than among Caucasians (0.69 and 0.71, respectively). Our results directly confirmed the different distribution of these polymorphisms between the Asian and Caucasian ethnic groups. Given the genetic differences found between groups, two points became evident: first, ethnic variations may have implications for the interpretation of results in circadian rhythm association studies, and second, the question may be raised about which evolutionary conditions shaped these genetic clock variations.

Étude génétique et fonctionnelle de variantes de la région chromosomique 3p21 associée aux maladies inflammatoires de l'intestin

Des études de liaison et d’association génétiques ont permis d’identifier certains des facteurs de risque génétiques aux maladies inflammatoires de l’intestin (MII) dans la région chromosomique 3p21. Dans cette région, le polymorphisme nucléotidique simple (SNP) codant non-synonyme du gène MST1, rs3197999, encodant pour la mutation R689C, a été associé et répliqué à la fois à la colite ulcéreuse (CU) et à la maladie de Crohn (MC). Un autre SNP, corrélé à des SNP codants non-synonymes du gène MST1R, a également été associé à la MC. Afin de déterminer si d’autres variantes des gènes MST1 et MST1R sont associés à la CU, nous avons testé pour association des SNP de ces gènes. Seul un proxy de R689C a montré un signal d’association significatif aux MII, ce qui suggère que R689C est la variante causale aux MII dans le gène MST1. En cherchant à déterminer si la région 3p21 contenait plusieurs signaux d’association mutuellement indépendants, trois SNP ont été identifiés comme possible facteurs de risque indépendants, et ont été génotypés dans des cas de CU et de MC et des témoins, puis nos résultats d’association ont été combinés à ceux provenant de trois autres cohortes indépendantes. Les trois SNP, R689C (MST1), rs6802890 et rs7629936 (CDHR4), sont associés aux MII, mais une étude d’association conditionnelle suggère qu’il existe en fait deux signaux d’association mutuellement indépendants dans la région 3p21. Le signal principal provient de R689C, une mutation de la protéine MSP. Cette protéine a un rôle dans l’inflammation chez les macrophages murins, et la migration, la cicatrisation et la survie chez les cellules épithéliales. Dans cette étude, le rôle de la MSP a été investigué dans des modèles de macrophages humains et de cellules épithéliales de côlon, et seule la phosphorylation d’AKT, un acteur dans la voie de signalisation de la survie cellulaire, a été modulée par la MSP dans nos modèles. Ce projet a donc permis d’apporter des connaissances sur les facteurs de risques génétiques aux MII dans la région 3p21, en identifiant 2 signaux d’association indépendants, et en nous informant sur le rôle de MST1, duquel provient le signal d’association principal, chez les cellules humaines.

Analyse génotypique des cellules initiatrices de tumeurs exprimant CD133 dans le neuroblastome

Le neuroblastome (NB) est la tumeur solide extracranienne la plus fréquente et mortelle chez les jeunes enfants. Il se caractérise par une résistance à la chimiothérapie possiblement en partie dû à la présence de cellules initiatrices de tumeurs (TICs). Des études ont mis en évidence le rôle de CD133 comme un marqueur des TICs dans divers types de cancers. Les buts de notre travail étaient d’abord de démontrer les vertus de TICs des cellules exprimant CD133 et ensuite, en utilisant une analyse globale du génome avec des polymorphismes nucléotidiques simples (SNPs), d’effectuer une analyse différentielle entre les TICs et les autres cellules du NB afin d’en identifier les anomalies génétiques spécifiques. Des lignées cellulaires de NB ont été triées par cytométrie de flux afin d’obtenir deux populations: une enrichie en CD133 (CD133high), l’autre faible en CD133 (CD133low). Afin de déterminer si ces populations cellulaires présentent des propriétés de TICs, des essais sur les neurosphères, les colonies en agar mou et les injections orthotopiques de 500 cellules sélectionnées dans 11 souris ont été réalisées. Après une isolation de l’ADN des populations sélectionnées, nous avons effectué une analyse génotypique par SNP utilisant les puces « Affymetrix Genome-Wide Human SNP Array 6.0 ». Pour vérifier l’expression des gènes identifiés, des Western Blots ont été réalisés. Nos résultats ont démontré que la population CD133 avait des propriétés de TICs in vitro et in vivo. L’analyse génotypique différentielle a permis d’identifier deux régions communes (16p13.3 and 19p13.3) dans la population CD133high ayant des gains et deux autres régions (16q12.1 and 21q21.3) dans la population CD133low possédant des pertes d’hétérozygoties (LOH). Aucune perte n’a été observée. Parmi les gènes étudiés, l’expression protéique d’éphrine-A2 était corrélée à celle de CD133 dans 6 tumeurs et 2 lignées cellulaires de NB. De plus, l’augmentation de la concentration d’anticorps anti-éphrine-A2 dans le milieu diminue la taille des neurosphères. Ainsi, la population CD133high, qui a des vertus de TICs, possède des caractéristiques génotypiques différentes par rapport à celle CD133low. La présence d’éphrine-A2 dans les cellules exprimant CD133 souligne son importance dans le développement des TICs. Ces résultats suggèrent la présence de potentielle cible pour de nouvelles thérapeutiques ciblant les TICs mise en évidence par l’étude génomique.

Further characterization and determination of the single amino acid change in the tsI138 reovirus thermosensitive mutant

Many temperature-sensitive mutants have been isolated in early studies of mammalian reovirus. However, the bio- logical properties and nature of the genetic alterations remain incompletely explored for most of these mutants. The mutation harbored by the tsI138 mutant was already assigned to the L3 gene encoding the l1 protein. In the present study, this mu- tant was further studied as a possible tool to establish the role of the putative l1 enzymatic activities in viral multiplication. It was observed that synthesis of viral proteins is only marginally reduced, while it was difficult to recover viral particles at the nonpermissive temperature. A single nucleotide substitution resulting in an amino acid change was found; the position of this amino acid is consistent with a probable defect in assembly of the inner capsid at the nonpermissive temperature.

Preliminary population study of Methylenetetrahydrofolate Reductase (MTHFR) C677T polymorphism determination in a pilot group of students from the University of Rosario

 Introduction: the 5, 10-methylenetetrahydrofolate reductase (MTHFR) is an essential enzyme in folate metabolism; their polymorphisms have been associated with heart disease risk increase, obstetric problems, neural tube defects in fetuses and cancer susceptibility. This gene has a single nucleotide polymorphism, a C-T change at nucleotide 677, which affects significantly its enzymatic activity. Objective: because of the biological importance of this enzyme and the Colombian population genetic heterogeneity characteristic, a study was performed to determine allele and genotype frequencies of MTHFR C677T polymorphism in healthy individuals, taking into account that in Colombia there are only studies that have involved case-control methodology. Methods: we analyzed this polymorphism trough the amplification of the DNA of a 206 students sample population. Additionally, Colombian overall frequencies were calculated, using data from healthy controls reported in other studies. Results: a Hardy-Weinberg disequilibri m was found in the sample tested. For the Colombian data, we found that the global population was in equilibrium. Conclusion: T allele population frequency seems to be under positive selection pressure, which is reflected in the population allele increase, despite its deleterious effect. A Spanish study reported similar results and identified folic acid supplementation on expectant mothers as a probably cause of this change. 

Enfermedad tiroidea autoinmune en pacientes colombianos con lupus eritematoso sistémico

Objetivos: Determinar la prevalencia y los factores asociados con el desarrollo de hipotiroidismo autoinmune (HA) en una cohorte de pacientes con lupus eritematoso sistémico (LES), y analizar la información actual en cuanto a la prevalencia e impacto de la enfermedad tiroidea autoinmune y la autoinmunidad tiroidea en pacientes con LES. Métodos: Este fue un estudio realizado en dos pasos. Primero, un total de 376 pacientes con LES fueron evaluados sistemáticamente por la presencia de: 1) HA confirmado, 2) positividad para anticuerpos tiroperoxidasa/tiroglobulina (TPOAb/TgAb) sin hipotiroidismo, 3) hipotiroidismo no autoinmune, y 4) pacientes con LES sin hipotiroidismo ni positividad para TPOAb/TgAb. Se construyeron modelos multivariados y árboles de regresión y clasificación para analizar los datos. Segundo, la información actual fue evaluada a través de una revisión sistemática de la literatura (RLS). Se siguieron las guías PRISMA para la búsqueda en las bases de datos PubMed, Scopus, SciELO y Librería Virtual en Salud. Resultados: En nuestra cohorte, la prevalencia de HA confirmado fue de 12% (Grupo 1). Sin embargo, la frecuencia de positividad para TPOAb y TgAb fue de 21% y 10%, respectivamente (Grupo 2). Los pacientes con LES sin HA, hipotiroidismo no autoinmune ni positividad para TPOAb/TgAb constituyeron el 40% de la corhorte. Los pacientes con HA confirmada fueron estadísticamente significativo de mayor edad y tuvieron un inicio tardío de la enfermedad. El tabaquismo (ORA 6.93, IC 95% 1.98-28.54, p= 0.004), la presencia de Síndrome de Sjögren (SS) (ORA 23.2, IC 95% 1.89-359.53, p= 0.015) y la positividad para anticuerpos anti-péptido cíclico citrulinado (anti-CCP) (ORA 10.35, IC 95% 1.04-121.26, p= 0.047) se asociaron con la coexistencia de LES-HA, ajustado por género y duración de la enfermedad. El tabaquismo y el SS fueron confirmados como factores predictivos para LES-HA (AUC del modelo CART = 0.72). En la RSL, la prevalencia de ETA en LES varío entre 1% al 60%. Los factores asociados con esta poliautoinmunidad fueron el género femenino, edad avanzada, tabaquismo, positividad para algunos anticuerpos, SS y el compromiso articular y cutáneo. Conclusiones: La ETA es frecuente en pacientes con LES, y no afecta la severidad del LES. Los factores de riesgo identificados ayudarán a los clínicos en la búsqueda de ETA. Nuestros resultados deben estimular políticas para la suspensión del tabaquismo en pacientes con LES.

Sequential genome-wide association studies for monitoring adverse events in the clinical evaluation of new drugs

Pharmacovigilance, the monitoring of adverse events (AEs), is an integral part in the clinical evaluation of a new drug. Until recently, attempts to relate the incidence of AEs to putative causes have been restricted to the evaluation of simple demographic and environmental factors. The advent of large-scale genotyping, however, provides an opportunity to look for associations between AEs and genetic markers, such as single nucleotides polymorphisms (SNPs). It is envisaged that a very large number of SNPs, possibly over 500 000, will be used in pharmacovigilance in an attempt to identify any genetic difference between patients who have experienced an AE and those who have not. We propose a sequential genome-wide association test for analysing AEs as they arise, allowing evidence-based decision-making at the earliest opportunity. This gives us the capability of quickly establishing whether there is a group of patients at high-risk of an AE based upon their DNA. Our method provides a valid test which takes account of linkage disequilibrium and allows for the sequential nature of the procedure. The method is more powerful than using a correction, such as idák, that assumes that the tests are independent. Copyright © 2006 John Wiley & Sons, Ltd.

Sequentially testing for a gene-drug interaction in a genomewide analysis

Assaying a large number of genetic markers from patients in clinical trials is now possible in order to tailor drugs with respect to efficacy. The statistical methodology for analysing such massive data sets is challenging. The most popular type of statistical analysis is to use a univariate test for each genetic marker, once all the data from a clinical study have been collected. This paper presents a sequential method for conducting an omnibus test for detecting gene-drug interactions across the genome, thus allowing informed decisions at the earliest opportunity and overcoming the multiple testing problems from conducting many univariate tests. We first propose an omnibus test for a fixed sample size. This test is based on combining F-statistics that test for an interaction between treatment and the individual single nucleotide polymorphism (SNP). As SNPs tend to be correlated, we use permutations to calculate a global p-value. We extend our omnibus test to the sequential case. In order to control the type I error rate, we propose a sequential method that uses permutations to obtain the stopping boundaries. The results of a simulation study show that the sequential permutation method is more powerful than alternative sequential methods that control the type I error rate, such as the inverse-normal method. The proposed method is flexible as we do not need to assume a mode of inheritance and can also adjust for confounding factors. An application to real clinical data illustrates that the method is computationally feasible for a large number of SNPs. Copyright (c) 2007 John Wiley & Sons, Ltd.

Influence of apoA-V gene variants on postprandial triglyceride metabolism: impact of gender

Although apolipoprotein AN (apoA-V) polymorphisms have been consistently associated with fasting triglyceride (TG) levels, their impact on postprandial lipemia remains relatively unknown. In this study, we investigate the impact of two common apoA-V polymorphisms (-1131 T>C and S19W) and apoA-V haplotypes on fasting and postprandial lipid metabolism in adults in the United Kingdom (n = 259). Compared with the wild-type TT, apoA-V -1131 TC heterozygotes had 15% (P = 0.057) and 21% (P = 0.002) higher fasting TG and postprandial TG area under the curve (AUC), respectively. Significant (P = 0.038) and nearly significant (P = 0.057) gender X genotype interactions were observed for fasting TG and TG AUC, with a greater impact of genotype in males. Lower HDL-cholesterol was associated with the rare TC genotype (P = 0.047). Significant linkage disequilibrium was found between the apoA-V -1131 T>C and the apoC-III 3238 C>G variants, with univariate analysis indicating an impact of this apoC-III single nucleotide polymorphism (SNP) on TG AUC (P = 0.015). However, in linear regression analysis, a significant independent association with TG AUC (P = 0.007) was only evident for the apoA-V -1131 T>C SNP, indicating a greater relative importance of the apoA-V genotype.

Detection of somaclonal variation during cocoa somatic embryogenesis characterised using cleaved amplified polymorphic sequence and the new freeware Artbio

The scarcity and stochastic nature of genetic mutations presents a significant challenge for scientists seeking to characterise de novo mutation frequency at specific loci. Such mutations can be particularly numerous during regeneration of plants from in vitro culture and can undermine the value of germplasm conservation efforts. We used cleaved amplified polymorphic sequence (CAPS) analysis to characterise new mutations amongst a clonal population of cocoa plants regenerated via a somatic embryogenesis protocol used previously for cocoa cryopreservation. Efficacy of the CAPS system for mutation detection was greatly improved after an ‘a priori’ in silico screen of reference target sequences for actual and potential restriction enzyme recognition sites using a new freely available software called Artbio. Artbio surveys known sequences for existing restriction enzyme recognition sites but also identifies all single nucleotide polymorphism (SNP) deviations from such motifs. Using this software, we performed an in silico screen of seven loci for restriction sites and their potential mutant SNP variants that were possible from 21 restriction enzymes. The four most informative locus-enzyme combinations were then used to survey the regenerant populations for de novo mutants. We characterised the pattern of point mutations and, using the outputs of Artbio, calculated the ratio of base substitution in 114 somatic embryo-derived cocoa regenerants originating from two explant genotypes. We found 49 polymorphisms, comprising 26.3% of the samples screened, with an inferred rate of 2.8 × 10−3 substitutions/screened base. This elevated rate is of a similar order of magnitude to previous reports of de novo microsatellite length mutations arising in the crop and suggests caution should be exercised when applying somatic embryogenesis for the conservation of plant germplasm.