Diversity of community acquired MRSA carrying the PVL gene in Queensland and New South Wales, Australia

Emergence and dissemination of community acquired methicillin resistant Staphylococcus aureus (CA-MRSA) strains are being reported with increasing frequency in Australia and worldwide. These strains of CA-MRSA are genetically diverse and distinct in Australia. Genotyping of CA-MRSA using eight highly-discriminatory single nucleotide polymorphisms (SNPs) is a rapid and robust method for monitoring the dissemination of these strains in the community. In this study, a SNP genotyping method was used to investigate the molecular epidemiology of 249 community acquired non-multiresistant MRSA (nm-MRSA) isolates over a 12-month period from routine diagnostic specimens. A real-time PCR for the presence of Panton-Valentine leukocidin (PVL) was also performed on these isolates. The CA-MRSA isolates were sourced from a large private laboratory in Brisbane, Australia that serves a wide geographic region encompassing Queensland and Northern New South Wales. This study identified 16 different STs and 98% of the CA-MRSA isolates were positive for the PVL gene. The most common ST was ST93 with 41% of isolates testing positive for this clone.

Phylogeographic reconstruction of a bacterial species with high levels of lateral gene transfer

Background Phylogeographic reconstruction of some bacterial populations is hindered by low diversity coupled with high levels of lateral gene transfer. A comparison of recombination levels and diversity at seven housekeeping genes for eleven bacterial species, most of which are commonly cited as having high levels of lateral gene transfer shows that the relative contributions of homologous recombination versus mutation for Burkholderia pseudomallei is over two times higher than for Streptococcus pneumoniae and is thus the highest value yet reported in bacteria. Despite the potential for homologous recombination to increase diversity, B. pseudomallei exhibits a relative lack of diversity at these loci. In these situations, whole genome genotyping of orthologous shared single nucleotide polymorphism loci, discovered using next generation sequencing technologies, can provide very large data sets capable of estimating core phylogenetic relationships. We compared and searched 43 whole genome sequences of B. pseudomallei and its closest relatives for single nucleotide polymorphisms in orthologous shared regions to use in phylogenetic reconstruction. Results Bayesian phylogenetic analyses of >14,000 single nucleotide polymorphisms yielded completely resolved trees for these 43 strains with high levels of statistical support. These results enable a better understanding of a separate analysis of population differentiation among >1,700 B. pseudomallei isolates as defined by sequence data from seven housekeeping genes. We analyzed this larger data set for population structure and allele sharing that can be attributed to lateral gene transfer. Our results suggest that despite an almost panmictic population, we can detect two distinct populations of B. pseudomallei that conform to biogeographic patterns found in many plant and animal species. That is, separation along Wallace's Line, a biogeographic boundary between Southeast Asia and Australia. Conclusion We describe an Australian origin for B. pseudomallei, characterized by a single introduction event into Southeast Asia during a recent glacial period, and variable levels of lateral gene transfer within populations. These patterns provide insights into mechanisms of genetic diversification in B. pseudomallei and its closest relatives, and provide a framework for integrating the traditionally separate fields of population genetics and phylogenetics for other bacterial species with high levels of lateral gene transfer.

Genomic evidence for the parallel evolution of coastal forms in the Senecio lautus complex

Instances of parallel ecotypic divergence where adaptation to similar conditions repeatedly cause similar phenotypic changes in closely related organisms are useful for studying the role of ecological selection in speciation. Here we used a combination of traditional and next generation genotyping techniques to test for the parallel divergence of plants from the Senecio lautus complex, a phenotypically variable groundsel that has adapted to disparate environments in the South Pacific. Phylogenetic analysis of a broad selection of Senecio species showed that members of the S. lautus complex form a distinct lineage that has diversified recently in Australasia. An inspection of thousands of polymorphisms in the genome of 27 natural populations from the S. lautus complex in Australia revealed a signal of strong genetic structure independent of habitat and phenotype. Additionally, genetic differentiation between populations was correlated with the geographical distance separating them, and the genetic diversity of populations strongly depended on geographical location. Importantly, coastal forms appeared in several independent phylogenetic clades, a pattern that is consistent with the parallel evolution of these forms. Analyses of the patterns of genomic differentiation between populations further revealed that adjacent populations displayed greater genomic heterogeneity than allopatric populations and are differentiated according to variation in soil composition. These results are consistent with a process of parallel ecotypic divergence in face of gene flow.

No association between MTHFR A1298C and MTRR A66G polymorphisms, and MS in an Australian cohort

Multiple sclerosis (MS) is a complex neurological disease that affects the central nervous system (CNS) resulting in debilitating neuropathology. Pathogenesis is primarily defined by CNS inflammation and demyelination of nerve axons. Methionine synthase reductase (MTRR) is an enzyme that catalyzes the remethylation of homocysteine (Hcy) to methionine via cobalamin and folate dependant reactions. Cobalamin acts as an intermediate methyl carrier between methylenetetrahydrofolate reductase (MTHFR) and Hcy. MTRR plays a critical role in maintaining cobalamin in an active form and is consequently an important determinant of total plasma Hcy (pHcy) concentrations. Elevated intracellular pHcy levels have been suggested to play a role in CNS dysfunction, neurodegenerative, and cerebrovascular diseases. Our investigation entailed the genotyping of a cohort of 140 cases and matched controls for MTRR and MTHFR, by restriction length polymorphism (RFLP) techniques. Two polymorphisms: MTRR A66G and MTHFR A1298C were investigated in an Australian age and gender matched case-control study. No significant allelic frequency difference was observed between cases and controls at the α = 0.05 level (MTRR χ2 = 0.005, P = 0.95, MTHFR χ2 = 1.15, P = 0.28). Our preliminary findings suggest no association between the MTRR A66G and MTHFR A1298C polymorphisms and MS

Polymorphisms in inflammation pathway genes and endometrial cancer risk

BACKGROUND Experimental and epidemiologic evidence have suggested that chronic inflammation may play a critical role in endometrial carcinogenesis. METHODS To investigate this hypothesis, a two-stage study was carried out to evaluate single-nucleotide polymorphisms (SNP) in inflammatory pathway genes in association with endometrial cancer risk. In stage I, 64 candidate pathway genes were identified and 4,542 directly genotyped or imputed SNPs were analyzed among 832 endometrial cancer cases and 2,049 controls, using data from the Shanghai Endometrial Cancer Genetics Study. Linkage disequilibrium of stage I SNPs significantly associated with endometrial cancer (P < 0.05) indicated that the majority of associations could be linked to one of 24 distinct loci. One SNP from each of the 24 loci was then selected for follow-up genotyping. Of these, 21 SNPs were successfully designed and genotyped in stage II, which consisted of 10 additional studies including 6,604 endometrial cancer cases and 8,511 controls. RESULTS Five of the 21 SNPs had significant allelic odds ratios (ORs) and 95% confidence intervals (CI) as follows: FABP1, 0.92 (0.85-0.99); CXCL3, 1.16 (1.05-1.29); IL6, 1.08 (1.00-1.17); MSR1, 0.90 (0.82-0.98); and MMP9, 0.91 (0.87-0.97). Two of these polymorphisms were independently significant in the replication sample (rs352038 in CXCL3 and rs3918249 in MMP9). The association for the MMP9 polymorphism remained significant after Bonferroni correction and showed a significant association with endometrial cancer in both Asian- and European-ancestry samples. CONCLUSIONS These findings lend support to the hypothesis that genetic polymorphisms in genes involved in the inflammatory pathway may contribute to genetic susceptibility to endometrial cancer. Impact statement: This study adds to the growing evidence that inflammation plays an important role in endometrial carcinogenesis.

Single nucleotide polymorphisms and mRNA expression of VEGF-A in papillary thyroid carcinoma: Potential markers for aggressive phenotypes

BACKGROUND AND OBJECTIVES Polymorphisms of the VEGF gene are known to affect the biological behaviour of cancers but have seldom been studied in thyroid cancer. The aim of the current study is to evaluate the prevalence and relevance of VEGF-A polymorphisms and mRNA expression in papillary thyroid carcinoma (PTC). MATERIALS AND METHODS Genomic DNA and total RNA were isolated from paraffin-embedded tissue from 91 PTC (51 conventional PTC and 40 follicular variant) and 78 control thyroid tissues. Three DNA polymorphisms (+936C > T, +405C > G and -141A > C) in the 3' and 5' untranslated region (3'-UTR, 5'-UTR) of VEGF-A were studied using PCR and RFLP. Also, the mRNA expression of VEGF-A in these tissues was studied by real-time PCR. RESULTS Distribution of polymorphisms in the 5'-UTR (VEGF-VEGF -141A > C and +405C > G) and 3'-UTR (VEGF +936C > T) were all significantly different in PTC and benign thyroid tissue (p = 0.0001, 0.001 and 0.028 respectively). The VEGF -141 C allele was more common in PTC with lymph node metastases (p = 0.026). VEGF + 405 Galleles andVEGF +936 CC genotype were more common in PTC of advanced pathological staging (p = 0.018 and 0.017 respectively). Also, increased expression of VEGF-A mRNA was noted in PTC compared to control (p = 0.009). Within the group of patients with conventional PTC, those with lymph nodal metastases had a higher level of VEGF-A mRNA expression than other patients (p = 0.0003). CONCLUSION These findings suggest that VEGF polymorphisms and mRNA expression may predict the aggressiveness behaviour of thyroid cancer.

Association of IL23R and ERAP1 genes with ankylosing spondylitis in a Portuguese population

Objective: Association between ankylosing spondylitis (AS) and two genes, ERAP1 and IL23R, has recently been reported in North American and British populations. The population attributable risk fraction for ERAP1 in this study was 25%, and for IL23R, 9%. Confirmation of these findings to ERAP1 in other ethnic groups has not yet been demonstrated. We sought to test the association between single nucleotide polymorphisms (SNPs) in these genes and susceptibility to AS among a Portuguese population. We also investigated the role of these genes in clinical manifestations of AS, including age of symptom onset, the Bath Ankylosing Spondylitis Disease Activity, Metrology and Functional Indices, and the modified Stoke Ankylosing Spondylitis Spinal Score. Methods: The study was conducted on 358 AS cases and 285 ethnically matched Portuguese healthy controls. AS was defined according to the modified New York Criteria. Genotyping of IL23R and ERAP1 allelic variants was carried out with TaqMan allelic discrimination assays. Association analysis was performed using the Cochrane-Armitage and linear regression tests of genotypes as implemented in PLINK for dichotomous and quantitative variables respectively. A meta-analysis for Portuguese and previously published Spanish IL23R data was performed using the StatsDirect® Statistical tools, by fixed and random effects models. Results: A total of 14 nsSNPs markers (8 for IL23R, 5 for ERAPl, 1 for LN-PEP) were analysed. Three markers (2 for IL23R and 1 for ERAP1) showed significant single-locus disease associations, confirming that the association of these genes with AS in the Portuguese population. The strongest associated SNP in IL23R was rs1004819 (OR=1.4, p=0.0049), and in ERAP1 was rs30187 (OR=1.26, p=0.035). The population attributable risk fractions in the Portuguese population for these SNPs are 11% and 9.7% respectively. No association was seen with any SNP in LN-PEP, which flanks ERAP1 and was associated with AS in the British population. No association was seen with clinical manifestations of AS. Conclusions: These results show that IL23R and ERAP1 genes are also associated with susceptibility to AS in the Portuguese population, and that they contribute a significant proportion of the population risk for this disease.

Effect of an estrogen receptor-α intron 4 polymorphism on fat mass in 11-year-old children

Context: in the ESR1 gene encoding estrogen receptor (ER)-α may be associated with fat mass in adults. Objectives: The objective of the study was to establish whether ESR1 polymorphisms influence fat mass in childhood. Design: This was a cross-sectional analysis after genotyping of rs9340799, rs2234693, and rs7757956 ESR1 polymorphisms. Setting: The Avon Longitudinal Study of Parents and Children (ALSPAC) was a population-based prospective study. Participants: Participants included 3097 11-yr-old children with results for ESR1 genotyping, puberty measures, and dual-energy x-ray absorptiometry results. Outcomes: Relationships between ESR1 polymorphisms and indices of body composition were measured. Results: The rs7757956 polymorphism was associated with fat mass (P = 0.002). Total body fat mass (adjusted for height) was reduced by 6% in children with TA/AA genotypes, and risk of being overweight (≥85th centile of fat mass) was decreased by 20%. This genetic effect appeared to interact with puberty in girls (P = 0.05 for interaction): in those with the TT genotype, total body fat mass (adjusted for height) was 18% higher in Tanner stages 3-5 vs. stages 1-2; the equivalent difference was 7% in those with TA/AA genotypes. Furthermore, the risk of being overweight was 36% lower in girls with TA/AA genotypes in Tanner stages 3-5, but no reduction was seen in those in stages 1-2. Neither rs9340799 nor rs2234693 polymorphisms were associated with body composition measures. Conclusions: Fat mass in 11-yr-old children was related to the rs7757956 ESR1 polymorphism. This association was strongest in girls in more advanced puberty, in whom the risk of being overweight was reduced by 36% in those with the TA/AA genotype.

Genetic association analysis of miRNA SNPs implicates MIR145 in breast cancer susceptibility

Background MicroRNAs (miRNAs) are important small non-coding RNA molecules that regulate gene expression in cellular processes related to the pathogenesis of cancer. Genetic variation in miRNA genes could impact their synthesis and cellular effects and single nucleotide polymorphisms (SNPs) are one example of genetic variants studied in relation to breast cancer. Studies aimed at identifying miRNA SNPs (miR-SNPs) associated with breast malignancies could lead towards further understanding of the disease and to develop clinical applications for early diagnosis and treatment. Methods We genotyped a panel of 24 miR-SNPs using multiplex PCR and chip-based matrix assisted laser desorption ionization time-of-flight (MALDI-TOF) mass spectrometry (MS) analysis in two Caucasian breast cancer case control populations (Primary population: 173 cases and 187 controls and secondary population: 679 cases and 301 controls). Association to breast cancer susceptibility was determined using chi-square (X 2 ) and odds ratio (OR) analysis. Results Statistical analysis showed six miR-SNPs to be non-polymorphic and twelve of our selected miR-SNPs to have no association with breast cancer risk. However, we were able to show association between rs353291 (located in MIR145) and the risk of developing breast cancer in two independent case control cohorts (p = 0.041 and p = 0.023). Conclusions Our study is the first to report an association between a miR-SNP in MIR145 and breast cancer risk in individuals of Caucasian background. This finding requires further validation through genotyping of larger cohorts or in individuals of different ethnicities to determine the potential significance of this finding as well as studies aimed to determine functional significance. Keywords: Association analysis; Breast cancer; microRNA; miR-SNPs; MIR145

Genotyping for genetic association studies: Methods and applications

In this thesis, two separate single nucleotide polymorphism (SNP) genotyping techniques were set up at the Finnish Genome Center, pooled genotyping was evaluated as a screening method for large-scale association studies, and finally, the former approaches were used to identify genetic factors predisposing to two distinct complex diseases by utilizing large epidemiological cohorts and also taking environmental factors into account. The first genotyping platform was based on traditional but improved restriction-fragment-length-polymorphism (RFLP) utilizing 384-microtiter well plates, multiplexing, small reaction volumes (5 µl), and automated genotype calling. We participated in the development of the second genotyping method, based on single nucleotide primer extension (SNuPeTM by Amersham Biosciences), by carrying out the alpha- and beta tests for the chemistry and the allele-calling software. Both techniques proved to be accurate, reliable, and suitable for projects with thousands of samples and tens of markers. Pooled genotyping (genotyping of pooled instead of individual DNA samples) was evaluated with Sequenom s MassArray MALDI-TOF, in addition to SNuPeTM and PCR-RFLP techniques. We used MassArray mainly as a point of comparison, because it is known to be well suited for pooled genotyping. All three methods were shown to be accurate, the standard deviations between measurements being 0.017 for the MassArray, 0.022 for the PCR-RFLP, and 0.026 for the SNuPeTM. The largest source of error in the process of pooled genotyping was shown to be the volumetric error, i.e., the preparation of pools. We also demonstrated that it would have been possible to narrow down the genetic locus underlying congenital chloride diarrhea (CLD), an autosomal recessive disorder, by using the pooling technique instead of genotyping individual samples. Although the approach seems to be well suited for traditional case-control studies, it is difficult to apply if any kind of stratification based on environmental factors is needed. Therefore we chose to continue with individual genotyping in the following association studies. Samples in the two separate large epidemiological cohorts were genotyped with the PCR-RFLP and SNuPeTM techniques. The first of these association studies concerned various pregnancy complications among 100,000 consecutive pregnancies in Finland, of which we genotyped 2292 patients and controls, in addition to a population sample of 644 blood donors, with 7 polymorphisms in the potentially thrombotic genes. In this thesis, the analysis of a sub-study of pregnancy-related venous thromboses was included. We showed that the impact of factor V Leiden polymorphism on pregnancy-related venous thrombosis, but not the other tested polymorphisms, was fairly large (odds ratio 11.6; 95% CI 3.6-33.6), and increased multiplicatively when combined with other risk factors such as obesity or advanced age. Owing to our study design, we were also able to estimate the risks at the population level. The second epidemiological cohort was the Helsinki Birth Cohort of men and women who were born during 1924-1933 in Helsinki. The aim was to identify genetic factors that might modify the well known link between small birth size and adult metabolic diseases, such as type 2 diabetes and impaired glucose tolerance. Among ~500 individuals with detailed birth measurements and current metabolic profile, we found that an insertion/deletion polymorphism of the angiotensin converting enzyme (ACE) gene was associated with the duration of gestation, and weight and length at birth. Interestingly, the ACE insertion allele was also associated with higher indices of insulin secretion (p=0.0004) in adult life, but only among individuals who were born small (those among the lowest third of birth weight). Likewise, low birth weight was associated with higher indices of insulin secretion (p=0.003), but only among carriers of the ACE insertion allele. The association with birth measurements was also found with a common haplotype of the glucocorticoid receptor (GR) gene. Furthermore, the association between short length at birth and adult impaired glucose tolerance was confined to carriers of this haplotype (p=0.007). These associations exemplify the interaction between environmental factors and genotype, which, possibly due to altered gene expression, predisposes to complex metabolic diseases. Indeed, we showed that the common GR gene haplotype associated with reduced mRNA expression in thymus of three individuals (p=0.0002).

Microsatellite polymorphisms of Sichuan golden monkeys

Previous study using protein electrophoresis shows no polymorphism in 44 nuclear loci of Sichuan golden monkey (Rhinopithecus roxellana), which limits our understandings of its population genetic patterns in the nuclear genome. In order to obtain sufficient information, we scanned 14 microsatellite loci in a sample of 32 individuals from its three major habitats (Minshan, Qinling and Shennongiia). A considerable amount of polymorphisms were detected. The average heterozygosities in the local populations were all above 0.5. The differentiations among local populations were significant. There was evidence of geneflow among subpopullations, but geneflow between Qinling and Shennongjia local populations was the weakest. Minshan and Qinling populations might have gone through recent bottlenecks. The estimation of the ratio of the effective population sizes among local populations was close to that from census sizes. Comparisons to available mitochondria data suggested that R. roxellana's social structures played an important role in shaping its population genetic patterns. Our study showed that the polymorphism level of R. roxellana was no higher than other endangered species; therefore, measures should be taken to preserve genetic diversity of this species.

Replication of association between schizophrenia and ZNF804A in the Irish Case-Control Study of Schizophrenia sample

A recent genome-wide association study reported association between schizophrenia and the ZNF804A gene on chromosome 2q32.1. We attempted to replicate these findings in our Irish Case-Control Study of Schizophrenia (ICCSS) sample (N=1021 cases, 626 controls). Following consultation with the original investigators, we genotyped three of the most promising single-nucleotide polymorphisms (SNPs) from the Cardiff study. We replicate association with rs1344706 (trend test one-tailed P=0.0113 with the previously associated A allele) in ZNF804A. We detect no evidence of association with rs6490121 in NOS1 (one-tailed P=0.21), and only a trend with rs9922369 in RGRIP1L (one-tailed P=0.0515). On the basis of these results, we completed genotyping of 11 additional linkage disequilibrium-tagging SNPs in ZNF804A. Of 12 SNPs genotyped, 11 pass quality control criteria and 4 are nominally associated, with our most significant evidence of association at rs7597593 (P=0.0013) followed by rs1344706. We observe no evidence of differential association in ZNF804A on the basis of family history or sex of case. The associated SNP rs1344706 lies in approximately 30 bp of conserved mammalian sequence, and the associated A allele is predicted to maintain binding sites for the brain-expressed transcription factors MYT1l and POU3F1/OCT-6. In controls, expression is significantly increased from the A allele of rs1344706 compared with the C allele. Expression is increased in schizophrenic cases compared with controls, but this difference does not achieve statistical significance. This study replicates the original reported association of ZNF804A with schizophrenia and suggests that there is a consistent link between the A allele of rs1344706, increased expression of ZNF804A and risk for schizophrenia.

Identification of domestication-related loci associated with flowering time and seed size in soybean with the RAD-seq genotyping method

Flowering time and seed size are traits related to domestication. However, identification of domestication-related loci/genes of controlling the traits in soybean is rarely reported. In this study, we identified a total of 48 domestication-related loci based on RAD-seq genotyping of a natural population comprising 286 accessions. Among these, four on chromosome 12 and additional two on chromosomes 11 and 15 were associated with flowering time, and four on chromosomes 11 and 16 were associated with seed size. Of the five genes associated with flowering time and the three genes associated with seed size, three genes Glyma11g18720, Glyma11g15480 and Glyma15g35080 were homologous to Arabidopsis genes, additional five genes were found for the first time to be associated with these two traits. Glyma11g18720 and Glyma05g28130 were co-expressed with five genes homologous to flowering time genes in Arabidopsis, and Glyma11g15480 was co-expressed with 24 genes homologous to seed development genes in Arabidopsis. This study indicates that integration of population divergence analysis, genome-wide association study and expression analysis is an efficient approach to identify candidate domestication-related genes.

Genes of cathepsin L-like proteases in Trypanosoma rangeli isolates: Markers for diagnosis, genotyping and phylogenetic relationships

We have sequenced genes encoding cathepsin L-like (CatL-like) cysteine proteases from isolates of Trypanosoma rangeli from humans, wild mammals and Rhodnius species of Central and South America. Phylogenetic trees of sequences encoding mature CatL-like enzymes of T rangeli and homologous genes from other trypanosomes, Leishmania spp. and bodonids positioned sequences of T rangeli (rangelipain) closest to T cruzi (cruzipain). Phylogenetic tree of kinetoplastids based on sequences of CatL-like was totally congruent with those derived from SSU rRNA and gGAPDH genes. Analysis of sequences from the CatL-like catalytic domains of 17 isolates representative of the overall phylogenetic diversity and geographical range of T rangeli supported all the lineages (A-D) previously defined using ribosomal and spliced leader genes. Comparison of the proteolytic activities of T rangeli isolates revealed heterogeneous banding profiles of cysteine proteases in gelatin gels, with differences even among isolates of the same lineage. CatL-like sequences proved to be excellent targets for diagnosis and genotyping of T rangeli by PCR. Data from CatL-like encoding genes agreed with results from previous studies of kDNA markers, and ribosomal and spliced leader genes, thereby corroborating clonal evolution, independent transmission cycles and the divergence of T rangeli lineages associated with sympatric species of Rhodnius. (c) 2009 Elsevier B.V. All rights reserved.

Asthma Polymorphic Loci genotyping in Madeira population: identification of potential genetic markers of disease susceptibility, severity and clinical relevance

Asthma is a complex disease, influenced by both environmental and genetic factors. In this study, the analysis of multiple environmental factos assessed by questionnaire and the genotyping of SNPs IL131c.144 G/A, IL41590 C/T, IL41RP2 253183, ADRB21c.16 A/G, ADAM331V4 C/G, ADAM331S1 c.710 G/A, GSDML1236 C/T and STAT6121 C/T were performed in a sample of Madeiran asthmatic patients and their families, and their association to asthma susceptibility and severity was assessed. Family, environmental, social and individual factos such as the presence of rhinitis in one of the parents,the habitation conditions, the family smoking habits, individual food habits and allergen sensitivity, were found to account for asthma severity. IL41590*T and IL41RP2*183$ alleles as well as the combined genotypes IL41590*CT/IL41590*TT and IL41 RP2*253183/IL41RP2*253183 were associated to both asthma susceptibility and severity.GSDML1236*TT was found associated only to asthma severity.Allele ADAM331 V4*C was significantly overM transmitted to asthmatic offspring being linked with the disease by TDT. These findings suggest that in addition to environmental influences, IL41 590 C/T, IL41RP2 253183, ADAM331V4 C/G and GSDML1236 C/T SNPs may constitute important genetic factos contributing to asthmasusceptibility and/or severity in Madeira population.