A systematic evaluation of normalization methods and probe replicability using infinium EPIC methylation data.

BACKGROUND: The Infinium EPIC array measures the methylation status of > 850,000 CpG sites. The EPIC BeadChip uses a two-array design: Infinium Type I and Type II probes. These probe types exhibit different technical characteristics which may confound analyses. Numerous normalization and pre-processing methods have been developed to reduce probe type bias as well as other issues such as background and dye bias. METHODS: This study evaluates the performance of various normalization methods using 16 replicated samples and three metrics: absolute beta-value difference, overlap of non-replicated CpGs between replicate pairs, and effect on beta-value distributions. Additionally, we carried out Pearson's correlation and intraclass correlation coefficient (ICC) analyses using both raw and SeSAMe 2 normalized data. RESULTS: The method we define as SeSAMe 2, which consists of the application of the regular SeSAMe pipeline with an additional round of QC, pOOBAH masking, was found to be the best performing normalization method, while quantile-based methods were found to be the worst performing methods. Whole-array Pearson's correlations were found to be high. However, in agreement with previous studies, a substantial proportion of the probes on the EPIC array showed poor reproducibility (ICC < 0.50). The majority of poor performing probes have beta values close to either 0 or 1, and relatively low standard deviations. These results suggest that probe reliability is largely the result of limited biological variation rather than technical measurement variation. Importantly, normalizing the data with SeSAMe 2 dramatically improved ICC estimates, with the proportion of probes with ICC values > 0.50 increasing from 45.18% (raw data) to 61.35% (SeSAMe 2).

Genome-wide meta-analysis associates GPSM1 with type 2 diabetes, a plausible gene involved in skeletal muscle function.

Genome-wide association studies (GWASs) have identified many genetic variations associated with type 2 diabetes mellitus (T2DM) in Asians, but understanding the functional genetic variants that influence traits is often a complex process. In this study, fine mapping and other analytical strategies were performed to investigate the effects of G protein signaling modulator 1 (GPSM1) on insulin resistance in skeletal muscle. A total of 128 single-nucleotide polymorphisms (SNPs) within GPSM1 were analysed in 21,897 T2DM cases and 32,710 healthy controls from seven GWASs. The SNP rs28539249 in intron 9 of GPSM1 showed a nominally significant association with T2DM in Asians (OR = 1.07, 95% CI = 1.04-1.10, P < 10-4). The GPSM1 mRNA was increased in skeletal muscle and correlated with T2DM traits across obese mice model. An eQTL for the cis-acting regulation of GPSM1 expression in human skeletal muscle was identified for rs28539249, and the increased GPSM1 expression related with T2DM traits within GEO datasets. Another independent Asian cohort showed that rs28539249 is associated with the skeletal muscle expression of CACFD1, GTF3C5, SARDH, and FAM163B genes, which are functionally enriched for endoplasmic reticulum stress (ERS) and unfolded protein response (UPR) pathways. Moreover, rs28539249 locus was predicted to disrupt regulatory regions in human skeletal muscle with enriched epigenetic marks and binding affinity for CTCF. Supershift EMSA assays followed luciferase assays demonstrated the CTCF specifically binding to rs28539249-C allele leading to decreased transcriptional activity. Thus, the post-GWAS annotation confirmed the Asian-specific association of genetic variant in GPSM1 with T2DM, suggesting a role for the variant in the regulation in skeletal muscle.

Genome-wide association and meta-analysis in populations from Starr County, Texas, and Mexico City identify type 2 diabetes susceptibility loci and enrichment for expression quantitative trait loci in top signals.

AIMS/HYPOTHESIS: We conducted genome-wide association studies (GWASs) and expression quantitative trait loci (eQTL) analyses to identify and characterise risk loci for type 2 diabetes in Mexican-Americans from Starr County, TX, USA. METHOD: Using 1.8 million directly interrogated and imputed genotypes in 837 unrelated type 2 diabetes cases and 436 normoglycaemic controls, we conducted Armitage trend tests. To improve power in this population with high disease rates, we also performed ordinal regression including an intermediate class with impaired fasting glucose and/or glucose tolerance. These analyses were followed by meta-analysis with a study of 967 type 2 diabetes cases and 343 normoglycaemic controls from Mexico City, Mexico. RESULT: The top signals (unadjusted p value <1 × 10(-5)) included 49 single nucleotide polymorphisms (SNPs) in eight gene regions (PER3, PARD3B, EPHA4, TOMM7, PTPRD, HNT [also known as RREB1], LOC729993 and IL34) and six intergenic regions. Among these was a missense polymorphism (rs10462020; Gly639Val) in the clock gene PER3, a system recently implicated in diabetes. We also report a second signal (minimum p value 1.52 × 10(-6)) within PTPRD, independent of the previously implicated SNP, in a population of Han Chinese. Top meta-analysis signals included known regions HNF1A and KCNQ1. Annotation of top association signals in both studies revealed a marked excess of trans-acting eQTL in both adipose and muscle tissues. CONCLUSIONS/INTERPRETATION: In the largest study of type 2 diabetes in Mexican populations to date, we identified modest associations of novel and previously reported SNPs. In addition, in our top signals we report significant excess of SNPs that predict transcript levels in muscle and adipose tissues.

Genome-wide association study of type 2 diabetes in a sample from Mexico City and a meta-analysis of a Mexican-American sample from Starr County, Texas.

AIMS/HYPOTHESIS: We report a genome-wide association study of type 2 diabetes in an admixed sample from Mexico City and describe the results of a meta-analysis of this study and another genome-wide scan in a Mexican-American sample from Starr County, TX, USA. The top signals observed in this meta-analysis were followed up in the Diabetes Genetics Replication and Meta-analysis Consortium (DIAGRAM) and DIAGRAM+ datasets. METHODS: We analysed 967 cases and 343 normoglycaemic controls. The samples were genotyped with the Affymetrix Genome-wide Human SNP array 5.0. Associations of genotyped and imputed markers with type 2 diabetes were tested using a missing data likelihood score test. A fixed-effects meta-analysis including 1,804 cases and 780 normoglycaemic controls was carried out by weighting the effect estimates by their inverse variances. RESULTS: In the meta-analysis of the two Hispanic studies, markers showing suggestive associations (p < 10(-5)) were identified in two known diabetes genes, HNF1A and KCNQ1, as well as in several additional regions. Meta-analysis of the two Hispanic studies and the recent DIAGRAM+ dataset identified genome-wide significant signals (p < 5 × 10(-8)) within or near the genes HNF1A and CDKN2A/CDKN2B, as well as suggestive associations in three additional regions, IGF2BP2, KCNQ1 and the previously unreported C14orf70. CONCLUSIONS/INTERPRETATION: We observed numerous regions with suggestive associations with type 2 diabetes. Some of these signals correspond to regions described in previous studies. However, many of these regions could not be replicated in the DIAGRAM datasets. It is critical to carry out additional studies in Hispanic and American Indian populations, which have a high prevalence of type 2 diabetes.

MGEA5-14 polymorphism and type 2 diabetes in Mexico City.

A family-based study has recently reported that a variant located in intron 10 of the gene MGEA5 increases susceptibility to Type 2 Diabetes (T2D). We evaluated the distribution of this SNP in a sample of T2D patients (N = 271) and controls (N = 244) from Mexico City. The frequency of the T allele was higher in the cases (2.6%) than in the controls (1.8%). After adjusting for age, sex, BMI, education, and individual ancestry the odds ratio was 1.60 but the 95% confidence interval was wide and overlapped 1 (0.52-4.86, P-value : 0.404). In order to characterize the distribution of the MGEA5-14 polymorphism in the relevant parental populations, we genotyped this variant in European (and European Americans), West African, and Native American samples. The T-allele was present at a frequency of 2.3% in Spain, 4.2% in European Americans, and 13% in Western Africans, but was absent in two Native American samples from Mexico and Peru. Given the low frequency of the T-allele, further studies using large sample sizes will be required to confirm the role of this variant in T2D.

Association of TCF7L2 polymorphisms with type 2 diabetes in Mexico City.

Polymorphisms within the transcription factor 7-like 2 gene (TCF7L2) have been associated with type 2 diabetes (T2D) in several recent studies. We characterized three of these polymorphisms (rs12255372, rs7903146 and the microsatellite DG10S478) in an admixed sample of 286 patients with T2D and 275 controls from Mexico City. We also analyzed three samples representative of the relevant parental populations: Native Americans from the state of Guerrero (Mexico), Spanish from Valencia and Nigerians (Bini from the Edo region). In order to minimize potential confounding because of the presence of population stratification in the sample, we evaluated the association of the three TCF7L2 polymorphisms with T2D by using the program admixmap to fit a logistic regression model incorporating individual ancestry, sex, age, body mass index and education. The markers rs12255372, rs7903146 and DG10S478 are in tight disequilibrium in the Mexican sample. We observed a significant association between the single-nucleotide polymorphism (SNP) rs12255372 and the microsatellite DG10S478 with T2D in the Mexican sample [rs12255372, odds ratio (OR) = 1.78, p = 0.017; DG10S478, OR = 1.62, p = 0.041]. The SNP rs7903146 shows similar trends, but its association with T2D is not as strong (OR = 1.39, p = 0.152). Analysis of the parental samples, as well as other available data, indicates that there are substantial population frequency differences for these polymorphisms: The frequencies of the T2D risk factors are more than 20% higher in European and West African populations than in East Asian and Native American populations.

Implications of correlations between skin color and genetic ancestry for biomedical research.

Skin pigmentation is a central element of most discussions on 'race' and genetics. Research on the genetic basis of population variation in this phenotype, which is important in mediating both social experiences and environmental exposures, is sparse. We studied the relationship between pigmentation and ancestry in five populations of mixed ancestry with a wide range of pigmentation and ancestral proportions (African Americans from Washington, DC; African Caribbeans living in England; Puerto Ricans from New York; Mexicans from Guerrero; and Hispanics from San Luis Valley). The strength of the relationship between skin color and ancestry was quite variable, with the correlations ranging in intensity from moderately strong (Puerto Rico, rho = 0.633) to weak (Mexico, rho = 0.212). These results demonstrate the utility of ancestry-informative genetic markers and admixture methods and emphasize the need to be cautious when using pigmentation as a proxy of ancestry or when extrapolating the results from one admixed population to another.

Admixture in the Hispanics of the San Luis Valley, Colorado, and its implications for complex trait gene mapping.

Hispanic populations are a valuable resource that can and should facilitate the identification of complex trait genes by means of admixture mapping (AM). In this paper we focus on a particular Hispanic population living in the San Luis Valley (SLV) in Southern Colorado. We used a set of 22 Ancestry Informative Markers (AIMs) to describe the admixture process and dynamics in this population. AIMs are defined as genetic markers that exhibit allele frequency differences between parental populations >or=30%, and are more informative for studying admixed populations than random markers. The ancestral proportions of the SLV Hispanic population are estimated as 62.7 +/- 2.1% European, 34.1 +/- 1.9% Native American and 3.2 +/- 1.5% West African. We also estimated the ancestral proportions of individuals using these AIMs. Population structure was demonstrated by the excess association of unlinked markers, the correlation between estimates of admixture based on unlinked marker sets, and by a highly significant correlation between individual Native American ancestry and skin pigmentation (R2= 0.082, p < 0.001). We discuss the implications of these findings in disease gene mapping efforts.

Ancestral proportions and admixture dynamics in geographically defined African Americans living in South Carolina.

We analyzed admixture in samples of six different African-American populations from South Carolina: Gullah-speaking Sea Islanders in coastal South Carolina, residents of four counties in the "Low Country" (Berkeley, Charleston, Colleton, and Dorchester), and persons living in the city of Columbia, located in central South Carolina. We used a battery of highly informative autosomal, mtDNA, and Y-chromosome markers. Two of the autosomal markers (FY and AT3) are linked and lie 22 cM apart on chromosome 1. The results of this study indicate, in accordance with previous historical, cultural, and anthropological evidence, a very low level of European admixture in the Gullah Sea Islanders (m = 3.5 +/- 0.9%). The proportion of European admixture is higher in the Low Country (m ranging between 9. 9 +/- 1.8% and 14.0 +/- 1.9%), and is highest in Columbia (m = 17.7 +/- 3.1%). A sex-biased European gene flow and a small Native American contribution to the African-American gene pool are also evident in these data. We studied the pattern of pairwise allelic associations between the FY locus and the nine other autosomal markers in our samples. In the combined sample from the Low Country (N = 548), a high level of linkage disequilibrium was observed between the linked markers, FY and AT3. Additionally, significant associations were also detected between FY and 4 of the 8 unlinked markers, suggesting the existence of significant genetic structure in this population. A continuous gene flow model of admixture could explain the observed pattern of genetic structure. A test conditioning on the overall admixture of each individual showed association of ancestry between the two linked markers (FY and AT3), but not between any of the unlinked markers, as theory predicts. Thus, even in the presence of genetic structure due to continuous gene flow or some other factor, it is possible to differentiate associations due to linkage from spurious associations due to genetic structure.

Population structure in admixed populations: effect of admixture dynamics on the pattern of linkage disequilibrium.

Gene flow between genetically distinct populations creates linkage disequilibrium (admixture linkage disequilibrium [ALD]) among all loci (linked and unlinked) that have different allele frequencies in the founding populations. We have explored the distribution of ALD by using computer simulation of two extreme models of admixture: the hybrid-isolation (HI) model, in which admixture occurs in a single generation, and the continuous-gene-flow (CGF) model, in which admixture occurs at a steady rate in every generation. Linkage disequilibrium patterns in African American population samples from Jackson, MS, and from coastal South Carolina resemble patterns observed in the simulated CGF populations, in two respects. First, significant association between two loci (FY and AT3) separated by 22 cM was detected in both samples. The retention of ALD over relatively large (>10 cM) chromosomal segments is characteristic of a CGF pattern of admixture but not of an HI pattern. Second, significant associations were also detected between many pairs of unlinked loci, as observed in the CGF simulation results but not in the simulated HI populations. Such a high rate of association between unlinked markers in these populations could result in false-positive linkage signals in an admixture-mapping study. However, we demonstrate that by conditioning on parental admixture, we can distinguish between true linkage and association resulting from shared ancestry. Therefore, populations with a CGF history of admixture not only are appropriate for admixture mapping but also have greater power for detection of linkage disequilibrium over large chromosomal regions than do populations that have experienced a pattern of admixture more similar to the HI model, if methods are employed that detect and adjust for disequilibrium caused by continuous admixture.

Comparison of narrow-band reflectance spectroscopy and tristimulus colorimetry for measurements of skin and hair color in persons of different biological ancestry.

We have used two modern computerized handheld reflectometers, the Photovolt ColorWalk colorimeter (a tristimulus colorimeter; Photovolt, UMM Electronics, Indianapolis, IN) and the DermaSpectrometer (a specialized narrow-band reflectometer; Cortex Technology, Hadsund, Denmark), to compare two methods for the objective determination of skin and hair color. These instruments both determine color by measuring the intensity of reflected light of particular wavelengths. The Photovolt ColorWalk instrument does so by shining a white light and sensing the intensity of the reflected light with a linear photodiode array. The ColorWalk results can then be expressed in terms of several standard color systems, most importantly, the Commission International d'Eclairage (CIE) Lab system, in which any color can be described by three values: L*, the lightness; a*, the amount of green or red; and b*, the amount of yellow or blue. Instead of a white light and photodiodes, the DermaSpectrometer uses two light-emitting diodes (LEDs), one green and one red, to illuminate a surface, and then it records the intensity of the reflected light. The results of these readings are expressed in terms of erythema (E) and melanin (M) indices. We measured the unexposed skin of the inner upper arm, the exposed skin of the forehead, and the hair, of 80 persons using these two instruments. Since it is important for the application of these measures in anthropology that we understand their relationship across a number of different pigmentation levels, we sampled persons from several different groups, namely, European Americans (n = 55), African Americans (n = 9), South Asians (n = 7), and East Asians (n = 9). In these subjects, there is a very high correlation between L* and the M index for the inner arm (R(2) = 0.928, P < 0.001), the forehead (R(2) = 0.822, P < 0.001), and the hair (R(2) = 0.827, P < 0.001). The relationship between a* and the E index is complex and dependent on the pigmentation level. We conclude that while both types of instruments provide accurate estimates of pigment level in skin and hair, measurements using narrow-band instruments may be less affected by the greater redness of certain body sites due to increased vascularization.

Estimation of admixture and detection of linkage in admixed populations by a Bayesian approach: application to African-American populations.

We describe a novel method for analysis of marker genotype data from admixed populations, based on a hybrid of Bayesian and frequentist approaches in which the posterior distribution is generated by Markov chain simulation and score tests are obtained from the missing-data likelihood. We analysed data on unrelated individuals from eight African-American populations, genotyped at ten marker loci of which two (FY and AT3) are linked (22 cM apart). Linkage between these two loci was detected by testing for association of ancestry conditional on parental admixture. The strength of this association was consistent with European gene flow into the African-American population between five and nine generations ago. To mimic the mapping of an unknown gene in an 'affecteds- only' analysis, a binary trait was constructed from the genotype at the AT3 locus and a score test was shown to detect linkage of this 'trait' with the FY locus. Mis-specification of the ancestry-specific allele frequencies - the probabilities of each allelic state given the ancestry of the allele - was detected at three of the ten marker loci. The methods described here have wide application to the analysis of data from admixed populations, allowing the effects of linkage and population structure (variation of admixture between individuals) to be distinguished. With more markers and a more complex statistical model, genes underlying ethnic differences in disease risk could be mapped by this approach.

Estimating African American admixture proportions by use of population-specific alleles.

We analyzed the European genetic contribution to 10 populations of African descent in the United States (Maywood, Illinois; Detroit; New York; Philadelphia; Pittsburgh; Baltimore; Charleston, South Carolina; New Orleans; and Houston) and in Jamaica, using nine autosomal DNA markers. These markers either are population-specific or show frequency differences >45% between the parental populations and are thus especially informative for admixture. European genetic ancestry ranged from 6.8% (Jamaica) to 22.5% (New Orleans). The unique utility of these markers is reflected in the low variance associated with these admixture estimates (SEM 1.3%-2.7%). We also estimated the male and female European contribution to African Americans, on the basis of informative mtDNA (haplogroups H and L) and Y Alu polymorphic markers. Results indicate a sex-biased gene flow from Europeans, the male contribution being substantially greater than the female contribution. mtDNA haplogroups analysis shows no evidence of a significant maternal Amerindian contribution to any of the 10 populations. We detected significant nonrandom association between two markers located 22 cM apart (FY-null and AT3), most likely due to admixture linkage disequilibrium created in the interbreeding of the two parental populations. The strength of this association and the substantial genetic distance between FY and AT3 emphasize the importance of admixed populations as a useful resource for mapping traits with different prevalence in two parental populations.

Genetic structure of the population of Cabo Verde (west Africa): evidence of substantial European admixture.

The population of Cabo Verde was founded in the fifteenth century (1462), on the basis of slaves brought from the West African coast and a few Europeans, mainly from Portugal. The polymorphism of six red cell enzymes (ADA, AK1, ALAD, ESD, GLO1, and PGD) and ten plasma proteins (AHSG, BF, F13A, F13B, GC, HP, ORM, PLG, TBG, and TF) was studied in a sample of 268 individuals from Cabo Verde (West Africa). There is no statistical evidence of genetic heterogeneity between the two groups of islands which constitute the archipelago, Barlavento and Sotavento. The gene frequency distribution observed in Cabo Verde differs, in many markers, from that of West African populations, suggesting an important European influence. The proportion of Caucasian genes in the population of Cabo Verde has been calculated to be M = 0.3634 +/- 0.0510, and the considerable dispersion of the locus-specific admixture estimates seems to indicate random drift has also played a role in the evolution of the allele frequencies in the archipelago. Partition of the variance of the mean estimate in evolutionary and sampling variance shows the evolutionary variance is more than ten times higher than the sampling variance. When dendrograms are constructed on the basis of different genetic distances, the population of Cabo Verde clusters with Afro-Americans, forming a different group from the populations of the African continent. This is interpreted as a consequence of the importance of Caucasian admixture both in Afro-Americans and in the population of Cabo Verde.

Red cell enzymes and plasma protein polymorphisms in the population of Cabo Verde (West Africa).

The genetic polymorphism of six red cell enzymes (ADA, AK1, ALAD, ESD, GLO1, and PGD) and 10 plasma proteins (AHSG, BF, F13A, F13B, GC, HP, ORM, PLG, TBG, and TF) is analyzed in a sample of 268 unrelated individuals from Cabo Verde (West Africa). The population of Cabo Verde was founded in the 15th century (1462), on the basis of a great number of slaves brought from the West African coast and a few Europeans, mainly from Portugal. The frequencies found in Cabo Verde for the majority of the markers are intermediate between those reported for Africans and Europeans. Further, the presence of alleles which are rarely or never seen in Blacks, but are common in Caucasians, suggest a substantial contribution of Europeans to the gene pool of the population of Cabo Verde. © 1995 Wiley-Liss, Inc.

Distribution of alpha-2-HS-glycoprotein (AHSG) phenotypes in Cabo Verde (west Africa): description of a new allele, AHSG*32.

The genetic polymorphism of alpha-2-HS-glycoprotein (AHSG) was studied in the population of Cabo Verde (West Africa), using isoelectric focusing in polyacrylamide gels followed by immunofixation-silver stain. AHSG frequencies are reported for the first time in a subsaharan African population. In addition to the common variants, AHSG 1 and AHSG 2, five AHSG variants were observed, including a new variant, tentatively designated AHSG 32. The allele frequencies were, AHSG*1: 0.7289, AHSG*2: 0.2111, AHSG*10: 0.0276, AHSG*3: 0.0162, AHSG*11: 0.0081, AHSG*22: 0.0065, AHSG*32:0.0016.

[Osteoarticular infections associated with catheterization of the subclavian vein]. / Infecciones osteoarticulares asociadas a cateterización de la vena subclavia.

Access to a central way through catheterization of the subclavian vein is a widely used technique. Not uncommonly, the procedure is followed by infective complications among which clavicular osteomyelitis and septic sternoclavicular arthritis represent a rare eventuality. We report two cases of staphylococcic bacteremia produced after subclavian vein catheterization. Both patients presented septic sternoclavicular arthritis and osteomyelitis of the sternal manubrium. The isolated microorganisms were Staphylococcus aureus and Staphylococcus epidermidis methicillin-resistant, respectively. The clinical course under antibiotic therapy was satisfactory in both cases. Septic metastases appear to be the most likely pathogenic mechanism for the osteoarticular complications.