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.

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.

Genome-wide SNP genotyping highlights the role of natural selection in Plasmodium falciparum population divergence

Background: The malaria parasite Plasmodium falciparum exhibits abundant genetic diversity, and this diversity is key to its success as a pathogen. Previous efforts to study genetic diversity in P. falciparum have begun to elucidate the demographic history of the species, as well as patterns of population structure and patterns of linkage disequilibrium within its genome. Such studies will be greatly enhanced by new genomic tools and recent large-scale efforts to map genomic variation. To that end, we have developed a high throughput single nucleotide polymorphism (SNP) genotyping platform for P. falciparum. Results: Using an Affymetrix 3,000 SNP assay array, we found roughly half the assays (1,638) yielded high quality, 100% accurate genotyping calls for both major and minor SNP alleles. Genotype data from 76 global isolates confirm significant genetic differentiation among continental populations and varying levels of SNP diversity and linkage disequilibrium according to geographic location and local epidemiological factors. We further discovered that nonsynonymous and silent (synonymous or noncoding) SNPs differ with respect to within-population diversity, interpopulation differentiation, and the degree to which allele frequencies are correlated between populations. Conclusions: The distinct population profile of nonsynonymous variants indicates that natural selection has a significant influence on genomic diversity in P. falciparum, and that many of these changes may reflect functional variants deserving of follow-up study. Our analysis demonstrates the potential for new high-throughput genotyping technologies to enhance studies of population structure, natural selection, and ultimately enable genome-wide association studies in P. falciparum to find genes underlying key phenotypic traits.

Single-nucleotide polymorphism, linkage disequilibrium and geographic structure in the malaria parasite Plasmodium vivax: prospects for genome-wide association studies

Background: The ideal malaria parasite populations for initial mapping of genomic regions contributing to phenotypes such as drug resistance and virulence, through genome-wide association studies, are those with high genetic diversity, allowing for numerous informative markers, and rare meiotic recombination, allowing for strong linkage disequilibrium (LD) between markers and phenotype-determining loci. However, levels of genetic diversity and LD in field populations of the major human malaria parasite P. vivax remain little characterized. Results: We examined single-nucleotide polymorphisms (SNPs) and LD patterns across a 100-kb chromosome segment of P. vivax in 238 field isolates from areas of low to moderate malaria endemicity in South America and Asia, where LD tends to be more extensive than in holoendemic populations, and in two monkey-adapted strains (Salvador-I, from El Salvador, and Belem, from Brazil). We found varying levels of SNP diversity and LD across populations, with the highest diversity and strongest LD in the area of lowest malaria transmission. We found several clusters of contiguous markers with rare meiotic recombination and characterized a relatively conserved haplotype structure among populations, suggesting the existence of recombination hotspots in the genome region analyzed. Both silent and nonsynonymous SNPs revealed substantial between-population differentiation, which accounted for similar to 40% of the overall genetic diversity observed. Although parasites clustered according to their continental origin, we found evidence for substructure within the Brazilian population of P. vivax. We also explored between-population differentiation patterns revealed by loci putatively affected by natural selection and found marked geographic variation in frequencies of nucleotide substitutions at the pvmdr-1 locus, putatively associated with drug resistance. Conclusion: These findings support the feasibility of genome-wide association studies in carefully selected populations of P. vivax, using relatively low densities of markers, but underscore the risk of false positives caused by population structure at both local and regional levels.

Detection of the TLR4 1196C > T Polymorphism by Mismatched-Polymerase Chain Reaction Using Plasmid DNA as Internal Control in Restriction Fragment Length Polymorphism Assays

Background: Restriction fragment length polymorphism (RFLP) is a common molecular assay used for genotyping, and it requires validated quality control procedures to prevent mistyping caused by impaired endonuclease activity. We have evaluated the usefulness of a plasmid-based internal control in RFLP assays. Results: Blood samples were collected from 102 individuals with acute myocardial infarction (AMI) and 108 non-AMI individuals (controls) for DNA extraction and laboratory analyses. The 1196C> T polymorphism in the toll-like receptor 4 (TLR4) gene was amplified by mismatched-polymerase chain reaction (PCR). Amplicons and pBluescript II SK-plasmid were simultaneously digested with endonuclease HincII. Fragments were separated on 2% agarose gels. Plasmid was completely digested using up to 55.2 nmL/L DNA solutions and 1 mu L PCR product. Nevertheless, plasmid DNA with 41.4 nM or higher concentrations was incompletely digested in the presence of 7 mL PCR product. In standardized conditions, TLR4 1196C> T variant was accurately genotyped. TLR4 1196T allele frequency was similar between AMI (3.1%) and controls (2.0%, p = 0.948). TLR4 SNP was not associated with AMI in this sample population. In conclusion, the plasmid-based control is a useful approach to prevent mistyping in RFLP assays, and it is validate for genetic association studies such as TLR4 1196C> T.

Matrix metalloproteinase-2 polymorphism is associated with prognosis in prostate cancer

Objective: Prostate cancer (PCa) is the most frequent tumor in males in Brazil. Single nucleotide polymorphisms (SNP) have been demonstrated in the promoter region of matrix metalloproteinases (MMPs) genes and have been associated with development and progression of some cancers. In this study, our aim was to investigate a possible relation between polymorphism of the promoter region of the MMP2 gene and classical prognostic parameters in prostate cancer. Materials and methods: Genomic DNA was extracted using conventional protocols. The DNA sequence containing the polymorphic site was amplified by real-time polymerase chain reaction, using fluorescent probes (TaqMan). Results: In patients with tumors of a higher stage (pT3), a polymorphic allele in the MMP2 gene was more frequent (P = 0.026) than in patients with lower tumor stage. A polymorphic allele in the MMP2 gene was more frequent in Gleason >= 7 than in Gleason <= 6 (P = 0.042). Conclusions: We conclude that MMP2 polymorphism can be used together with pathological stage and Gleason score to identify patients with worse prognosis. Our results illustrate the potential use of MMP2 SNP as a molecular marker for prostate cancer. (C) 2010 Elsevier Inc. All rights reserved.

CXCL12 rs1801157 polymorphism and expression in peripheral blood from breast cancer patients

The role of chemokines has been extensively analyzed both in cancer risk and tumor progression. Among different cytokines, CXCR4 and its ligand CXCL12 have been recently subjected to a closer examination. The single-nucleotide polymorphism (SNP) rs1801157 (previously known as CXCL12-A/SDF1-3`A) in the CXCL12 gene and the relative expression of mRNA CXCL12 in peripheral blood were assessed in breast cancer patients, since the chemokine CXCL12 and its receptor CXCR4 regulate leukocyte trafficking and many essential biological processes, including tumor growth, angiogenesis and metastasis of different types of tumors. Genotyping was performed by PCR-RFLP (polymerase chain reaction followed by restriction fragment length polymorphism) using MspI restriction enzyme and the expression analyses by quantitative RT-PCR. No difference in GG genotype and allele A carrier frequencies were observed between breast cancer patients and healthy blood donors and nor when CXCL12 mRNA expression was assessed among patients with different tumor stages. However a significant difference was observed when CXCL12 mRNA relative expression was analyzed in breast cancer patients in accordance to the presence or absence of the CXCL12 rs1801157 allele A. Allele A breast cancer patients presented a mRNA CXCL12 expression about 2.1-fold smaller than GG breast cancer patients. Estrogen positive patients presenting CXCL12 allele A presented a significantly lower expression of CXCL12 in peripheral blood (p = 0.039) than GG hormone positive patients. Our findings demonstrated that allele A is associated with low expression of CXCL12 in the peripheral blood from ER-positive breast cancer patients, which suggests implications on breast cancer clinical outcome. (C) 2011 Elsevier Ltd. All rights reserved.

CXCL12 rs1801157 Polymorphism in Patients with Breast Cancer, Hodgkin`s Lymphoma, and Non-Hodgkin`s Lymphoma

Chemokines and their receptors regulate the trafficking of immune cells during their development, inflammation, and tissue repair. The single-nucleotide polymorphism (SNP) rs1801157 (previously known as CXCL12-A/ stromal cell-derived factor-1 (SDF1)-3`A) in CXCL12/SDF1 gene was assessed in breast cancer, Hodgkin`s lymphoma (HL), and non-Hodgkin`s lymphoma (NHL), since the chemokine CXCL12, previously known as SDF1, and its receptor CXCR4 regulate leukocyte trafficking and many essential biological processes, including tumor growth, angiogenesis, and metastasis of different types of tumors. Genotyping was performed by PCR-RFLP (polymerase chain reaction followed by restriction fragment length polymorphism) using a restriction enzyme Hpall cleavage. No significant difference was observed in genotype distribution between breast cancer patients (GG: 57.3%; GA: 39.8%; AA: 2.9%) and healthy female controls (GG: 62.9%; GA: 33%; AA: 4.1%) nor between HL patients (GG: 61.1%; GA:27.8%; AA: 11.1%) and healthy controls (GG: 65.6%; GA: 28.9%; AA: 5.5%), whereas a significant difference was observed in genotype distribution between NHL patients (GG: 51.4%; GA: 47.1%; AA: 1.5%) and healthy controls (GG: 65.6%; GA: 28.9%; AA: 5.5%). Further studies will be necessary to elucidate the cancer chemokine network. However, this study suggests that CXCL12 rs1801157 polymorphism may have important implications in the pathogenesis of NHL. J. Clin. Lab. Anal. 23:387-393, 2009. (C) 2009 Wiley-Liss, Inc.

Val/Leu(247) Polymorphism of beta 2-glycoprotein I in Brazilian Patients with Antiphospholipid Syndrome-A Genetic Risk Factor?

A genetic polymorphism of the beta 2-glycoprotein I (beta 2-GPI) is recognized by antiphospholipid antibodies (aPL) and may even play a role in the development of antiphospholipid syndrome (APS). The objectives of this study were to determine a Val/Leu SNP at position 247 of the beta 2-GPI gene in Brazilian patients with APS and to compare these data with clinical and laboratory manifestations. Polymorphism assignment was performed by PCR followed by Rsa I restriction endonuclease. The titration of anti-beta 2-GPI antibodies was detected by ELISA. The results showed significantly higher frequencies of the V-encoding allele and the homozygous VV genotype in patients with APS than in control subjects (OR = 1.781, P = 0.0068; and OR = 6.413, P < 0.0001, respectively). The frequency of this genotype was also significantly higher in patients with arterial and venous thrombosis than in the control group (52% and 44%, respectively, versus 13%). Anti-beta 2-GPI-positive patients had significantly higher frequencies of the VV genotype than the controls subjects (OR = 8.179, P < 0.0001). These results suggest that the V-encoding allele and the homozygous VV genotype at position 247 of the beta 2-GPI gene may play a role in the generation of anomalous beta 2-GPI, with consequent auto-antibody production, and in phenotype expression of arterial and venous thrombosis in APS patients.

Correlation between beta-2-glycoprotein I gene polymorphism and anti-beta-2 glycoprotein I antibodies in patients with multibacillary leprosy

Antiphospholipid antibodies, such as anti-beta 2-glycoprotein I (beta 2GPI), are present in multibacillary leprosy (MB) patients; however, MB patients do not usually present with antiphospholipid antibody syndrome (APS), which is characterized by thromboembolic phenomena (TEP). Rare cases of TEP occur in leprosy patients, but the physiopathology of this condition remains unclear. In this case-control study, we examined whether single-nucleotide polymorphisms (SNPs) of the beta 2GPI gene contributed to the risk of leprosy and APS co-morbidity. SNPs Ser88Asn, Leu247Val, Cys306Gly and Trp316Ser were identified in 113 Brazilian leprosy patients. Additionally, anti-beta 2GPI antibodies and plasma concentrations of beta 2GPI were quantified. The Ser88Asn, Cys306Gly and Trp316Ser SNPs were not risk factors for APS in leprosy. A higher frequency of Val/Val homozygosity was observed in leprosy patients compared to controls (36 vs. 5%; P < 0.001). Forty-two percent of MB and 17% of paucibacillary leprosy patients were positive for anti-beta 2GPI IgM (P = 0.014). There was no correlation between SNP Ser88Asn or Cys306Gly and anti-beta 2GPI antibody levels. In MB patients with positive anti-beta 2GPI IgM, the frequency of Val/Val homozygosity was higher than in controls (32 vs. 15%; P = 0.042). The frequency of the mutant allele Ser316 was higher in MB patients with positive rather than negative anti-beta 2GPI IgM levels (6 vs. 0%; P = 0.040) and was greater than in the control group (6 vs. 1%; P = 0.034). The studied polymorphisms did not influence the plasma concentrations of beta 2GPI. These results suggest that Leu247Val and Trp316Ser SNPs may represent genetic risk factors for anti-beta 2GPI antibody production in MB patients.

Association study of an epidermal growth factor gene functional polymorphism with the risk and prognosis of gliomas in Brazil

Epidermal growth factor (EGF) plays an important role in cancer. A functional single nucleotide polymorphism (SNP) in the 5`-untranslated region of the EGF gene (+61 A>G) may influence its expression and contribute to cancer predisposition and aggressiveness. Aiming to investigate the role of EGF +61 A>G in the susceptibility to glioma and its prognosis, we performed a case-control study with 165 patients and 200 healthy controls from Brazil. Comparisons of genotype distributions and allele frequencies did not reveal any significant differences between the groups. The mean overall survival was 9.2 months for A/A, 8.2 months for A/G, and 7.7 months for GIG. When survival curves were plotted we found that the +61G allele is associated with poor overall survival (p=0.023) but not with disease-free survival (p=0.527). Our data suggest that, although there is no association between the EGF +61 A>G genotype and glioma susceptibility, this SNP is associated with shorter overall survival of glioma patients in the Brazilian population. Nevertheless, future studies utilizing a larger series are essential for a definitive conclusion. (Int J Biol Markers 2009; 24: 277-81)

WRN Cys1367Arg SNP is not associated with risk and prognosis of gliomas in Southeast Brazil

Werner syndrome (WS) is a premature aging disorder characterized by early onset of symptoms related to normal aging and by a high predisposition to various types of cancer, including gliomas. WS is caused by inherited recessive mutations in the WRN gene, which encodes a helicase considered a caretaker of the genome. Aiming to study the role of WRN Cys1367Arg in glioma susceptibility and oncologic prognosis of patients, we investigated the genotype distribution of this single nucleotide polymorphism in 94 glioma patients and 100 healthy subjects. Comparisons of genotype distributions and allele frequencies did not reveal any significant difference between the groups. Overall and disease-free survival rates were calculated, but no statistically significant difference was observed. Our data suggest that WRN Cys1367Arg SNP is not involved either in susceptibility to developing gliomas or in patient survival, at least in the Brazilian population.

Polymorphisms at GLUT1 Gene Are Not Associated With the Development of TSP/HAM in Brazilian HTLV-1 Infected Individuals and the Discovery of a New Polymorphism at GLUT1 Gene

The development of HTLV-1 associated clinical manifestations, such as TSP/HAM and ATLL, occur in 2-4% of the infected population and it is still unclear why this infection remains asymptomatic in most infected carriers. Recently, it has been demonstrated that HTLV uses the Glucose transporter type 1 (GLUT1) to infect T-CD4(+) lymphocytes and that single nucleotide polymorphisms (SNP) in the GLUT1 gene are associated with diabetic nephropathy in patients with diabetes mellitus in different populations. These polymorphisms could contribute to a higher GLUT1 protein expression on cellular membrane, facilitating the entry of HTLV and its transmission cell by cell. This could result in a higher provirus load and consequently in the development of TSP/HAM. To evaluate the role of GLUT1 gene polymorphisms in the development of TSP/HAM in HTLV-1 infected individuals, the g.22999G > T, g.15339T > C and c.-2841A > T sites were analyzed by PCR/RFLP or sequencing in 244 infected individuals and 102 normal controls. The proviral load of the HTLV-1 infected patients was also analyzed using Real Time Quantitative PCR. Genotypic and allelic frequencies of the three sites did not differ significantly between controls and HTLV-1 infected individuals. There was no difference in genotypic and allelic distributions among patients as to the presence or absence of HTLV-1 associated clinic manifestations. As regards the quantification of the provirus load, we observed a significant reduction in the asymptomatic individuals compared with the oligosymptomatic and TSP/HAM individuals. These results suggest that g.22999G > T, g.15339T > C, and c.-2841A > T SNP do not contribute to HTLV-1 infection nor to the genetic susceptibility of TSP/HAM in Brazilian HTLV-1 infected individuals. J. Med. Virol. 81:552557, 2009. (C) 2009 Wiley-Liss, Inc.

The broad effects of the functional IL-10 promoter-592 polymorphism: modulation of IL-10, TIMP-3, and OPG expression and their association with periodontal disease outcome

Periodontal diseases are infectious diseases, in which periodontopathogens trigger chronic inflammatory and immune responses that lead to tissue destruction. It occurs through the generation of metalloproteinases and the activation of bone resorption mechanisms. Anti-inflammatory cytokines such as IL-10 seem to attenuate periodontal tissue destruction through the induction of tissue inhibitors of metalloproteinases (TIMPs) and the inhibitor of osteoclastogenesis osteoprotegerin (OPG). A high individual variation in levels of IL-10 mRNA is verified in periodontitis patients, which is possibly determined by genetic polymorphisms. In this study, the IL-10 promoter -592C/A single nucleotide polymorphism ( SNP), which is associated with a decrease in IL-10 production, was analyzed by RFLP in 116 chronic periodontitis (CP) patients and 173 control (C) subjects, and the IL-10, TIMPs, and OPG mRNA expression levels in diseased gingival tissues were determined by real-time-PCR. The IL-10-592 SNP CA (P=0.0012/OR=2.4/CI:1.4-4.1), AA (P=0.0458/OR=2.3/CI:1.1-4.9), and CA+AA (P=0.0006/OR=2.4/CI: 1.4-3.4) genotypes and the allele A (P=0.0036/OR=1.7/CI:1.2-2.4) were found to be significantly more prevalent in the CP group when compared with control subjects. Both CA and AA genotypes were associated with lower levels of IL-10, TIMP-3, and OPG mRNA expression in diseased periodontal tissues and were also associated with disease severity as mean pocket depth. Taken together, the results presented here demonstrate that IL10-592 SNP is functional in CP, being associated with lower levels of IL-10 mRNA expression, which is supposed to consequently decrease the expression of the downstream genes TIMP-3 and OPG, and influence periodontal disease outcome. J. Leukoc. Biol. 84: 1565-1573; 2008.

An interleukin-1 beta (IL-1 beta) single-nucleotide polymorphism at position 3954 and red complex periodontopathogens independently and additively modulate the levels of IL-1 beta in diseased periodontal tissues

Inflammatory cytokines such as interieukin-1 beta (IL-1 beta) are involved in the pathogenesis of periodontal diseases. A high individual variation in the levels of IL-10 mRNA has been verified, which is possibly determined by genetic polymorphisms and/or by the presence of periodontopathogens such as Porphyromonas gingivalis, Tannerella forsythia, Treponema denticola, and Aggregatibacter actinomycetemcomitans. In this study, we investigated the role of an IL-10 promoter single-nucleotide polymorphism at position 3954 [IL-1 beta(3954) SNP] and the presence of the periodontopathogens in the determination of the IL-1 beta levels in the periodontal tissues of nonsmoking chronic periodontitis (CP) patients (n = 117) and control (C) subjects in = 175) and the possible correlations with the clinical parameters of the disease. IL-1 beta(3954) SNP was investigated by restriction fragment length polymorphism, while the IL-1 beta levels and the presence of the periodontopathogens were determined by real-time PCR. Similar frequencies of IL-1 beta(3954) SNP were found in the C and CP groups, in spite of a trend toward a higher incidence of T alleles in the CP group. The IL-1 beta (3954) SNP CT and TT genotypes, as well as P. gingivalis, T. forsythia, and T. denticola, were associated with higher IL-1 beta levels and with higher values of the clinical parameters of disease severity. Concomitant analyses demonstrate that IL-1 beta(3954) and the red complex periodontopathogens were found to independently and additively modulate the levels of IL-1 beta in periodontal tissues. Similarly, the concurrent presence of both factors was associated with increased scores of disease severity. IL-1 beta(3954) genotypes and red complex periodontopathogens, individually and additively, modulate the levels of IL-1 beta in the diseased tissues of nonsmoking CP patients and, consequently, are potentially involved in the determination of the disease outcome.