High-throughput SNP genotyping in the highly heterozygous genome of Eucalyptus: assay success, polymorphism and transferability across species

Background: High-throughput SNP genotyping has become an essential requirement for molecular breeding and population genomics studies in plant species. Large scale SNP developments have been reported for several mainstream crops. A growing interest now exists to expand the speed and resolution of genetic analysis to outbred species with highly heterozygous genomes. When nucleotide diversity is high, a refined diagnosis of the target SNP sequence context is needed to convert queried SNPs into high-quality genotypes using the Golden Gate Genotyping Technology (GGGT). This issue becomes exacerbated when attempting to transfer SNPs across species, a scarcely explored topic in plants, and likely to become significant for population genomics and inter specific breeding applications in less domesticated and less funded plant genera. Results: We have successfully developed the first set of 768 SNPs assayed by the GGGT for the highly heterozygous genome of Eucalyptus from a mixed Sanger/454 database with 1,164,695 ESTs and the preliminary 4.5X draft genome sequence for E. grandis. A systematic assessment of in silico SNP filtering requirements showed that stringent constraints on the SNP surrounding sequences have a significant impact on SNP genotyping performance and polymorphism. SNP assay success was high for the 288 SNPs selected with more rigorous in silico constraints; 93% of them provided high quality genotype calls and 71% of them were polymorphic in a diverse panel of 96 individuals of five different species. SNP reliability was high across nine Eucalyptus species belonging to three sections within subgenus Symphomyrtus and still satisfactory across species of two additional subgenera, although polymorphism declined as phylogenetic distance increased. Conclusions: This study indicates that the GGGT performs well both within and across species of Eucalyptus notwithstanding its nucleotide diversity >= 2%. The development of a much larger array of informative SNPs across multiple Eucalyptus species is feasible, although strongly dependent on having a representative and sufficiently deep collection of sequences from many individuals of each target species. A higher density SNP platform will be instrumental to undertake genome-wide phylogenetic and population genomics studies and to implement molecular breeding by Genomic Selection in Eucalyptus.

Pharmacogenetics of OATP Transporters Reveals That SLCO1B1 c.388A > G Variant Is Determinant of Increased Atorvastatin Response

Aims: The relationship between variants in SLCO1B1 and SLCO2B1 genes and lipid-lowering response to atorvastatin was investigated. Material and Methods: One-hundred-thirty-six unrelated individuals with hypercholesterolemia were selected and treated with atorvastatin (10 mg/day/4 weeks). They were genotyped with a panel of ancestry informative markers for individual African component of ancestry (ACA) estimation by SNaPshot (R) and SLCO1B1 (c.388A>G, c.463C>A and c.521T>C) and SLCO2B1 (-71T>C) gene polymorphisms were identified by TaqMan (R) Real-time PCR. Results: Subjects carrying SLCO1B1 c.388GG genotype exhibited significantly high low-density lipoprotein (LDL) cholesterol reduction relative to c.388AA+c.388AG carriers (41 vs. 37%, p = 0.034). Haplotype analysis revealed that homozygous of SLCO1B1*15 (c.521C and c.388G) variant had similar response to statin relative to heterozygous and non-carriers. A multivariate logistic regression analysis confirmed that c.388GG genotype was associated with higher LDL cholesterol reduction in the study population (OR: 3.2, CI95%: 1.3-8.0, p < 0.05). Conclusion: SLCO1B1 c.388A>G polymorphism causes significant increase in atorvastatin response and may be an important marker for predicting efficacy of lipid-lowering therapy.

TCF7L2 Polymorphism rs7903146 Is Associated with Coronary Artery Disease Severity and Mortality

Background: TCF7L2 polymorphisms have been consistently associated with type 2 diabetes mellitus in different populations and type 2 diabetes mellitus is a major risk factor for cardiovascular disease, especially coronary artery disease. This study aimed to evaluate the association between TCF7L2 polymorphism rs7903146 and coronary artery disease in diabetic and non-diabetic subjects. Methods and Results: two populations were studied in order to assess severity of coronary artery disease and cardiovascular events incidence. Eight-hundred and eighty nine subjects who were referred for cardiac catheterization for coronary artery disease diagnosis were cross-sectionally evaluated for coronary lesions (atherosclerotic burden) and 559 subjects from the MASS-II Trial were prospectively followed-up for 5 years and assessed for major cardiovascular events incidence. As expected, rs7903146 T allele was associated with diabetes. Although diabetic patients had a higher prevalence of coronary lesions, no association between TCF7L2 genotype and coronary lesions was found in this subgroup. However, non-diabetic individuals carrying the T allele were associated with a significantly higher frequency of coronary lesions than non-diabetic non-carriers of the risk allele (adjusted OR = 2.32 95% CI 1.27-4.24, p = 0.006). Moreover, presence of multi-vessel coronary artery disease was also associated with the CT or TT genotypes in non-diabetics. Similarly, from the prospective sample analysis, non-diabetics carrying the CT/TT genotypes had significantly more composite cardiovascular end-points events than CC carriers (p = 0.049), mainly due to an increased incidence of death (p = 0.004). Conclusions: rs7903146 T allele is associated with diabetes and, in non-diabetic individuals, with a higher prevalence and severity of coronary artery disease and cardiovascular events. name of registry site (see list below), registration number, trial registration URL in brackets.

TCF7L2 variant genotypes and type 2 diabetes risk in Brazil: significant association, but not a significant tool for risk stratification in the general population

Background: Genetic polymorphisms of the TCF7L2 gene are strongly associated with large increments in type 2 diabetes risk in different populations worldwide. In this study, we aimed to confirm the effect of the TCF7L2 polymorphism rs7903146 on diabetes risk in a Brazilian population and to assess the use of this genetic marker in improving diabetes risk prediction in the general population. Methods: We genotyped the single nucleotide polymorphisms (SNP) rs7903146 of the TCF7L2 gene in 560 patients with known coronary disease enrolled in the MASS II (Medicine, Angioplasty, or Surgery Study) Trial and in 1,449 residents of Vitoria, in Southeast Brazil. The associations of this gene variant to diabetes risk and metabolic characteristics in these two different populations were analyzed. To access the potential benefit of using this marker for diabetes risk prediction in the general population we analyzed the impact of this genetic variant on a validated diabetes risk prediction tool based on clinical characteristics developed for the Brazilian general population. Results: SNP rs7903146 of the TCF7L2 gene was significantly associated with type 2 diabetes in the MASS-II population (OR = 1.57 per T allele, p = 0.0032), confirming, in the Brazilian population, previous reports of the literature. Addition of this polymorphism to an established clinical risk prediction score did not increased model accuracy (both area under ROC curve equal to 0.776). Conclusion: TCF7L2 rs7903146 T allele is associated with a 1.57 increased risk for type 2 diabetes in a Brazilian cohort of patients with known coronary heart disease. However, the inclusion of this polymorphism in a risk prediction tool developed for the general population resulted in no improvement of performance. This is the first study, to our knowledge, that has confirmed this recent association in a South American population and adds to the great consistency of this finding in studies around the world. Finally, confirming the biological association of a genetic marker does not guarantee improvement on already established screening tools based solely on demographic variables.

Genetic Variants of Diabetes Risk and Incident Cardiovascular Events in Chronic Coronary Artery Disease

Objective: To determine whether information from genetic risk variants for diabetes is associated with cardiovascular events incidence. Methods: From the about 30 known genes associated with diabetes, we genotyped single-nucleotide polymorphisms at the 10 loci most associated with type-2 diabetes in 425 subjects from the MASS-II Study, a randomized study in patients with multi-vessel coronary artery disease. The combined genetic information was evaluated by number of risk alleles for diabetes. Performance of genetic models relative to major cardiovascular events incidence was analyzed through Kaplan-Meier curve comparison and Cox Hazard Models and the discriminatory ability of models was assessed for cardiovascular events by calculating the area under the ROC curve. Results: Genetic information was able to predict 5-year incidence of major cardiovascular events and overall-mortality in non-diabetic individuals, even after adjustment for potential confounders including fasting glycemia. Non-diabetic individuals with high genetic risk had a similar incidence of events then diabetic individuals (cumulative hazard of 33.0 versus 35.1% of diabetic subjects). The addition of combined genetic information to clinical predictors significantly improved the AUC for cardiovascular events incidence (AUC = 0.641 versus 0.610). Conclusions: Combined information of genetic variants for diabetes risk is associated to major cardiovascular events incidence, including overall mortality, in non-diabetic individuals with coronary artery disease.

Poly (A)(+) Transcriptome Assessment of ERBB2-Induced Alterations in Breast Cell Lines

We report the first quantitative and qualitative analysis of the poly (A)(+) transcriptome of two human mammary cell lines, differentially expressing (human epidermal growth factor receptor) an oncogene over-expressed in approximately 25% of human breast tumors. Full-length cDNA populations from the two cell lines were digested enzymatically, individually tagged according to a customized method for library construction, and simultaneously sequenced by the use of the Titanium 454-Roche-platform. Comprehensive bioinformatics analysis followed by experimental validation confirmed novel genes, splicing variants, single nucleotide polymorphisms, and gene fusions indicated by RNA-seq data from both samples. Moreover, comparative analysis showed enrichment in alternative events, especially in the exon usage category, in ERBB2 over-expressing cells, data indicating regulation of alternative splicing mediated by the oncogene. Alterations in expression levels of genes, such as LOX, ATP5L, GALNT3, and MME revealed by large-scale sequencing were confirmed between cell lines as well as in tumor specimens with different ERBB2 backgrounds. This approach was shown to be suitable for structural, quantitative, and qualitative assessment of complex transcriptomes and revealed new events mediated by ERBB2 overexpression, in addition to potential molecular targets for breast cancer that are driven by this oncogene.

The Genomic Ancestry of Individuals from Different Geographical Regions of Brazil Is More Uniform Than Expected

Based on pre-DNA racial/color methodology, clinical and pharmacological trials have traditionally considered the different geographical regions of Brazil as being very heterogeneous. We wished to ascertain how such diversity of regional color categories correlated with ancestry. Using a panel of 40 validated ancestry-informative insertion-deletion DNA polymorphisms we estimated individually the European, African and Amerindian ancestry components of 934 self-categorized White, Brown or Black Brazilians from the four most populous regions of the Country. We unraveled great ancestral diversity between and within the different regions. Especially, color categories in the northern part of Brazil diverged significantly in their ancestry proportions from their counterparts in the southern part of the Country, indicating that diverse regional semantics were being used in the self-classification as White, Brown or Black. To circumvent these regional subjective differences in color perception, we estimated the general ancestry proportions of each of the four regions in a form independent of color considerations. For that, we multiplied the proportions of a given ancestry in a given color category by the official census information about the proportion of that color category in the specific region, to arrive at a ""total ancestry"" estimate. Once such a calculation was performed, there emerged a much higher level of uniformity than previously expected. In all regions studied, the European ancestry was predominant, with proportions ranging from 60.6% in the Northeast to 77.7% in the South. We propose that the immigration of six million Europeans to Brazil in the 19(th) and 20(th) centuries - a phenomenon described and intended as the ""whitening of Brazil"" -is in large part responsible for dissipating previous ancestry dissimilarities that reflected region-specific population histories. These findings, of both clinical and sociological importance for Brazil, should also be relevant to other countries with ancestrally admixed populations.

Disequilibrium Coefficient: A Bayesian Perspective

Hardy-Weinberg Equilibrium (HWE) is an important genetic property that populations should have whenever they are not observing adverse situations as complete lack of panmixia, excess of mutations, excess of selection pressure, etc. HWE for decades has been evaluated; both frequentist and Bayesian methods are in use today. While historically the HWE formula was developed to examine the transmission of alleles in a population from one generation to the next, use of HWE concepts has expanded in human diseases studies to detect genotyping error and disease susceptibility (association); Ryckman and Williams (2008). Most analyses focus on trying to answer the question of whether a population is in HWE. They do not try to quantify how far from the equilibrium the population is. In this paper, we propose the use of a simple disequilibrium coefficient to a locus with two alleles. Based on the posterior density of this disequilibrium coefficient, we show how one can conduct a Bayesian analysis to verify how far from HWE a population is. There are other coefficients introduced in the literature and the advantage of the one introduced in this paper is the fact that, just like the standard correlation coefficients, its range is bounded and it is symmetric around zero (equilibrium) when comparing the positive and the negative values. To test the hypothesis of equilibrium, we use a simple Bayesian significance test, the Full Bayesian Significance Test (FBST); see Pereira, Stern andWechsler (2008) for a complete review. The disequilibrium coefficient proposed provides an easy and efficient way to make the analyses, especially if one uses Bayesian statistics. A routine in R programs (R Development Core Team, 2009) that implements the calculations is provided for the readers.

Genetic Diversity on the Integrase Region of the pol Gene among HIV Type 1-Infected Patients Naive for Integrase Inhibitors in Sao Paulo City, Brazil

The presence of mutations associated with integrase inhibitor (INI) resistance among INI-naive patients may play an important clinical role in the use of those drugs Samples from 76 HIV-1-infected subjects naive to INIs were submitted to direct sequencing. No differences were found between naive (25%) subjects and subjects on HAART (75%). No primary mutation associated with raltegravir or elvitegravir resistance was found. However, 78% of sequences showed at least one accessory mutation associated with resistance. The analysis of the 76 IN sequences showed a high polymorphic level on this region among Brazilian HIV-1-infected subjects, including a high prevalence of aa substitutions related to INI resistance. The impact of these findings remains unclear and further studies are necessary to address these questions.

An assessment of population structure in eight breeds of cattle using a whole genome SNP panel

Background: Analyses of population structure and breed diversity have provided insight into the origin and evolution of cattle. Previously, these studies have used a low density of microsatellite markers, however, with the large number of single nucleotide polymorphism markers that are now available, it is possible to perform genome wide population genetic analyses in cattle. In this study, we used a high-density panel of SNP markers to examine population structure and diversity among eight cattle breeds sampled from Bos indicus and Bos taurus. Results: Two thousand six hundred and forty one single nucleotide polymorphisms ( SNPs) spanning all of the bovine autosomal genome were genotyped in Angus, Brahman, Charolais, Dutch Black and White Dairy, Holstein, Japanese Black, Limousin and Nelore cattle. Population structure was examined using the linkage model in the program STRUCTURE and Fst estimates were used to construct a neighbor-joining tree to represent the phylogenetic relationship among these breeds. Conclusion: The whole-genome SNP panel identified several levels of population substructure in the set of examined cattle breeds. The greatest level of genetic differentiation was detected between the Bos taurus and Bos indicus breeds. When the Bos indicus breeds were excluded from the analysis, genetic differences among beef versus dairy and European versus Asian breeds were detected among the Bos taurus breeds. Exploration of the number of SNP loci required to differentiate between breeds showed that for 100 SNP loci, individuals could only be correctly clustered into breeds 50% of the time, thus a large number of SNP markers are required to replace the 30 microsatellite markers that are currently commonly used in genetic diversity studies.

Relatedness and HLA-DRB1 typing may discriminate the magnitude of the genetic susceptibility to tuberculosis using a household contact model

Background Tuberculosis clusters in families may be due to increased household exposure, shared genetic factors, or both. Household contact studies are useful to control exposure because socioeconomic and environmental conditions are similar to all subjects, allowing the evaluation of the contribution of relatedness to disease development. Methods In this study, the familial aggregation of tuberculosis using relatedness and a specific inherited marker (HLA-DRB1) was evaluated. Fifty families, which had at least two cases of tuberculosis diagnosed within the past 5 years, were selected from a cohort of tuberculosis carried out in Recife, Brazil. The first case diagnosed was considered to be a primary case. The secondary attack rate of tuberculosis in household contacts was estimated according to the degree of relatedness. The relative risk of having tuberculosis based on the degree of relatedness household and the population attributable fraction to relatedness were also estimated. HLA-DRB1 typing and attributable etiologic/preventive fractions were calculated among sick and healthy household contacts. Results Compared to unrelated contacts, the relative risk for tuberculosis adjusted for age was 1.38 (95% CI 0.86 to 2.21). Relatedness contributed 23% to the development of tuberculosis at the population levels. The HLA-DRB1*04 allele group (OR = 2.44; p =0.0324; etiologic fraction =0.15) was overrepresented and the DRB1*15 allele group (OR=0.48; p=0.0488; protective fraction=0.19) was underrepresented among household contacts exhibiting tuberculosis. The presence of DRB1 shared alleles between primary cases and their contacts was a risk factor for tuberculosis (p=0.0281). Conclusion This household contact model together with the utilisation of two genetic variables permitted the evaluation of genetic factors contributing towards tuberculosis development.

Vascular Endothelial Growth Factor Haplotypes Associated with Childhood Obesity

Expansion of adipose tissue in obesity is associated with angiogenesis and adipose tissue mass depends on neovascularization. Vascular endothelial growth factor (VEGF) is the main angiogenic factor in the adipose tissue, and VEGF expression is tightly regulated at both transcriptional and translational levels. However, no previous study has tested the hypothesis that genetic polymorphisms in the VEGF gene could affect susceptibility to obesity. To test this hypothesis, we compared the distribution of genotypes and haplotypes including three VEGF genetic polymorphisms in obese children and adolescents with those found in healthy controls. We studied 172 healthy children and adolescents and 113 obese children and adolescents. Genotypes of three clinically relevant VEGF polymorphisms in the promoter region (C-2578A, G-1154A, and G-634C) of the VEGF gene were determined by TaqMan allele discrimination assay and real-time polymerase chain reaction. VEGF haplotypes were inferred using Haplo. stats and PHASE 2.1 programs. We found no differences in the distributions of VEGF genotypes and alleles (p > 0.05). However, the CAG haplotype was more frequent in the obese group than in the control group (4% versus 0%, respectively, in white subjects; p = 0.008; odds ratio 10.148 (95% confidence interval: 1.098-93.788). Our findings suggest that VEGF haplotypes affect susceptibility to obesity in children and adolescents.

Endothelial Nitric Oxide Genotypes and Haplotypes Are Not Associated with End-Stage Renal Disease

The identification of genetic markers associated with chronic kidney disease (CKD) may help to predict its development. Because reduced nitric oxide (NO) bioavailability and endothelial dysfunction are involved in CKD, genetic polymorphisms in the gene encoding the enzyme involved in NO synthesis (endothelial NO synthase [eNos]) may affect the susceptibility to CKD and the development of end-stage renal disease (ESRD). We compared genotype and haplotype distributions of three relevant eNOS polymorphisms (T(-786) C in the promoter region, Glu298Asp in exon 7, and 4b/4a in intron 4) in 110 healthy control subjects and 127 ESRD patients. Genotypes for the T(-786) C and Glu298Asp polymorphisms were determined by TaqMan (R) Allele Discrimination assay and real-time polymerase chain reaction. Genotypes for the intron 4 polymorphism were determined by polymerase chain reaction and fragment separation by electrophoresis. The software program PHASE 2.1 was used to estimate the haplotypes frequencies. We considered significant a probability value of p < 0.05/number of haplotypes (p < 0.05/8 = 0.0063). We found no significant differences between groups with respect to age, ethnicity, and gender. CKD patients had higher blood pressure, total cholesterol, and creatinine levels than healthy control subjects (all p < 0.05). Genotype and allele distributions for the three eNOS polymorphisms were similar in both groups (p > 0.05). We found no significant differences in haplotype distribution between groups (p > 0.05). The lack of significant associations between eNOS polymorphisms and ESRD suggests that eNOS polymorphisms may not be relevant to the genetic component of CKD that leads to ESRD.

Endothelial Nitric Oxide Synthase Haplotypes Associated with Hypertension Do Not Predispose to Cardiac Hypertrophy

Left ventricular hypertrophy (LVH) is a complication that may result from chronic hypertension. While nitric oxide (NO) deficiency has been associated with LVH, inconsistent results have been reported with regards to the association of endothelial NO synthase (eNOS) polymorphisms and LVH in hypertensive patients. This study aims to assess whether eNOS haplotypes are associated with LVH in hypertensive patients. This study included 101 healthy controls and 173 hypertensive patients submitted to echocardiography examination. Genotypes for three eNOS polymorphisms were determined: a single-nucleotide polymorphism in the promoter region (T-786C) and in exon 7 (Glu298Asp), and variable number of tandem repeats in intron 4. We found no significant association between eNOS genotypes and hypertension or with LVH (all p>0.05). However, while we found two eNOS haplotypes associated with variable risk of hypertension (all p<0.05), we found no significant associations between eNOS haplotypes and LVH (all p>0.05), even after adjustment in multiple linear regression analysis. These findings suggest that eNOS haplotypes that have been associated with variable susceptibility to hypertension were not associated with LVH in hypertensive patients. Further studies are necessary to examine whether other genes downstream may interact with eNOS polymorphisms and predispose to LVH in hypertensive patients.

GSTP1 Ile105Val polymorphism in astrocytomas and glioblastomas

Glutathione S-transferases (GSTs) constitute a superfamily of ubiquitous multifunctional enzymes that are involved in the cellular detoxification of a large number of endogenous and exogenous chemical agents that have electrophilic functional groups. People who have deficiencies in this family of genes are at increased risk of developing some types of tumors. We examined GSTP1 Ile105Val polymorphism using PCR-RFLP in 80 astrocytoma and glioblastoma samples. Patients who had the Val allele of the GSTP1 Ile105Val polymorphism had an increased risk of tumor development (odds ratio = 8.60; 95% confidence interval = 4.74-17.87; P < 0.001). Overall survival of patients did not differ significantly. We suggest that GSTP1 Ile105Val polymorphisms are involved in susceptibility to developing astrocytomas and glioblastomas.