Evaluation of the new advanced 15-loci MIRU-VNTR genotyping tool in Mycobacterium tuberculosis molecular epidemiology studies

Background. During the last few years, PCR-based methods have been developed to simplify and reduce the time required for genotyping Mycobacterium tuberculosis (MTB) by standard approaches based on IS6110-Restriction Fragment Length Polymorphism (RFLP). Of these, MIRU-12-VNTR (Mycobacterial interspersed repetitive units- variable number of tandem repeats) (MIRU-12) has been considered a good alternative. Nevertheless, some limitations and discrepancies with RFLP, which are minimized if the technique is complemented with spoligotyping, have been found. Recently, a new version of MIRU-VNTR targeting 15 loci (MIRU-15) has been proposed to improve the MIRU-12 format. Results. We evaluated the new MIRU-15 tool in two different samples. First, we analyzed the same convenience sample that had been used to evaluate MIRU-12 in a previous study, and the new 15-loci version offered higher discriminatory power (Hunter-Gaston discriminatory index [HGDI]: 0.995 vs 0.978; 34.4% of clustered cases vs 57.5%) and better correlation (full or high correlation with RFLP for 82% of the clusters vs 47%). Second, we evaluated MIRU-15 on a population-based sample and, once again, good correlation with the RFLP clustering data was observed (for 83% of the RFLP clusters). To understand the meaning of the discrepancies still found between MIRU-15 and RFLP, we analyzed the epidemiological data for the clustered patients. In most cases, splitting of RFLP-clustered patients by MIRU-15 occurred for those without epidemiological links, and RFLP-clustered patients with epidemiological links were also clustered by MIRU-15, suggesting a good epidemiological background for clustering defined by MIRU-15. Conclusion. The data obtained by MIRU-15 suggest that the new design is very efficient at assigning clusters confirmed by epidemiological data. If we add this to the speed with which it provides results, MIRU-15 could be considered a suitable tool for real-time genotyping.

Multiple sclerosis risk variant HLA-DRB1*1501 associates with high expression of DRB1 gene in different human populations.

The human leukocyte antigen (HLA) DRB1*1501 has been consistently associated with multiple sclerosis (MS) in nearly all populations tested. This points to a specific antigen presentation as the pathogenic mechanism though this does not fully explain the disease association. The identification of expression quantitative trait loci (eQTL) for genes in the HLA locus poses the question of the role of gene expression in MS susceptibility. We analyzed the eQTLs in the HLA region with respect to MS-associated HLA-variants obtained from genome-wide association studies (GWAS). We found that the Tag of DRB1*1501, rs3135388 A allele, correlated with high expression of DRB1, DRB5 and DQB1 genes in a Caucasian population. In quantitative terms, the MS-risk AA genotype carriers of rs3135388 were associated with 15.7-, 5.2- and 8.3-fold higher expression of DQB1, DRB5 and DRB1, respectively, than the non-risk GG carriers. The haplotype analysis of expression-associated variants in a Spanish MS cohort revealed that high expression of DRB1 and DQB1 alone did not contribute to the disease. However, in Caucasian, Asian and African American populations, the DRB1*1501 allele was always highly expressed. In other immune related diseases such as type 1 diabetes, inflammatory bowel disease, ulcerative colitis, asthma and IgA deficiency, the best GWAS-associated HLA SNPs were also eQTLs for different HLA Class II genes. Our data suggest that the DR/DQ expression levels, together with specific structural properties of alleles, seem to be the causal effect in MS and in other immunopathologies rather than specific antigen presentation alone.

Analysis of TNFAIP3, a feedback inhibitor of nuclear factor-κB and the neighbor intergenic 6q23 region in rheumatoid arthritis susceptibility

INTRODUCTION Genome-wide association studies of rheumatoid arthritis (RA) have identified an association of the disease with a 6q23 region devoid of genes. TNFAIP3, an RA candidate gene, flanks this region, and polymorphisms in both the TNFAIP3 gene and the intergenic region are associated with systemic lupus erythematosus. We hypothesized that there is a similar association with RA, including polymorphisms in TNFAIP3 and the intergenic region. METHODS To test this hypothesis, we selected tag-single nucleotide polymorphisms (SNPs) in both loci. They were analyzed in 1,651 patients with RA and 1,619 control individuals of Spanish ancestry. RESULTS Weak evidence of association was found both in the 6q23 intergenic region and in the TNFAIP3 locus. The rs582757 SNP and a common haplotype in the TNFAIP3 locus exhibited association with RA. In the intergenic region, two SNPs were associated, namely rs609438 and rs13207033. The latter was only associated in patients with anti-citrullinated peptide antibodies. Overall, statistical association was best explained by the interdependent contribution of SNPs from the two loci TNFAIP3 and the 6q23 intergenic region. CONCLUSIONS Our data are consistent with the hypothesis that several RA genetic factors exist in the 6q23 region, including polymorphisms in the TNFAIP3 gene, like that previously described for systemic lupus erythematosus.

Prospective universal application of mycobacterial interspersed repetitive-unit-variable-number tandem-repeat genotyping to characterize Mycobacterium tuberculosis isolates for fast identification of clustered and orphan cases.

The use of molecular tools for genotyping Mycobacterium tuberculosis isolates in epidemiological surveys in order to identify clustered and orphan strains requires faster response times than those offered by the reference method, IS6110 restriction fragment length polymorphism (RFLP) genotyping. A method based on PCR, the mycobacterial interspersed repetitive-unit-variable-number tandem-repeat (MIRU-VNTR) genotyping technique, is an option for fast fingerprinting of M. tuberculosis, although precise evaluations of correlation between MIRU-VNTR and RFLP findings in population-based studies in different contexts are required before the methods are switched. In this study, we evaluated MIRU-VNTR genotyping (with a set of 15 loci [MIRU-15]) in parallel to RFLP genotyping in a 39-month universal population-based study in a challenging setting with a high proportion of immigrants. For 81.9% (281/343) of the M. tuberculosis isolates, both RFLP and MIRU-VNTR types were obtained. The percentages of clustered cases were 39.9% (112/281) and 43.1% (121/281) for RFLP and MIRU-15 analyses, and the numbers of clusters identified were 42 and 45, respectively. For 85.4% of the cases, the RFLP and MIRU-15 results were concordant, identifying the same cases as clustered and orphan (kappa, 0.7). However, for the remaining 14.6% of the cases, discrepancies were observed: 16 of the cases clustered by RFLP analysis were identified as orphan by MIRU-15 analysis, and 25 cases identified as orphan by RFLP analysis were clustered by MIRU-15 analysis. When discrepant cases showing subtle genotypic differences were tolerated, the discrepancies fell from 14.6% to 8.6%. Epidemiological links were found for 83.8% of the cases clustered by both RFLP and MIRU-15 analyses, whereas for the cases clustered by RFLP or MIRU-VNTR analysis alone, links were identified for only 30.8% or 38.9% of the cases, respectively. The latter group of cases mainly comprised isolates that could also have been clustered, if subtle genotypic differences had been tolerated. MIRU-15 genotyping seems to be a good alternative to RFLP genotyping for real-time interventional schemes. The correlation between MIRU-15 and IS6110 RFLP findings was reasonable, although some uncertainties as to the assignation of clusters by MIRU-15 analysis were identified.

Characterization of Mycobacterium tuberculosis Beijing isolates from the Mediterranean area

BACKGROUND. The Beijing lineage of Mycobacterium tuberculosis is causing concern due to its global distribution and its involvement in severe outbreaks. Studies focused on this lineage are mainly restricted to geographical settings where its prevalence is high, whereas those in other areas are scarce. In this study, we analyze Beijing isolates in the Mediterranean area, where this lineage is not prevalent and is mainly associated with immigrant cases. RESULTS. Only 1% (N = 26) of the isolates from two population-based studies in Spain corresponded to Beijing strains, most of which were pan-susceptible and from Peruvian and Ecuadorian patients. Restriction fragment length polymorphism typing with the insertion sequence IS6110 identified three small clusters (2-3 cases). Mycobacterial interspersed repetitive unit-variable number tandem repeat typing (MIRU-15) offered low discriminatory power, requiring the introduction of five additional loci. A selection of the Beijing isolates identified in the Spanish sample, together with a sample of Beijing strains from Italy, to broaden the analysis context in the Mediterranean area, were assayed in an infection model with THP-1 cells. A wide range of intracellular growth rates was observed with only two isolates showing an increased intracellular replication, in both cases associated with contained production of TNF-alpha. No correlation was observed between virulence and the Beijing phylogenetic group, clustered/orphan status, or resistance. The Beijing strain responsible for extensive spread on Gran Canaria Island was also identified in Madrid, but did not lead to secondary cases and did not show high infectivity in the infection model. CONCLUSIONS. The Beijing lineage in our area is a non-homogeneous family, with only certain highly virulent representatives. The specific characterization of Beijing isolates in different settings could help us to accurately identify the virulent representatives before making general assumptions about this lineage.

Genome-wide association study of multiple sclerosis confirms a novel locus at 5p13.1.

Multiple Sclerosis (MS) is the most common progressive and disabling neurological condition affecting young adults in the world today. From a genetic point of view, MS is a complex disorder resulting from the combination of genetic and non-genetic factors. We aimed to identify previously unidentified loci conducting a new GWAS of Multiple Sclerosis (MS) in a sample of 296 MS cases and 801 controls from the Spanish population. Meta-analysis of our data in combination with previous GWAS was done. A total of 17 GWAS-significant SNPs, corresponding to three different loci were identified:HLA, IL2RA, and 5p13.1. All three have been previously reported as GWAS-significant. We confirmed our observation in 5p13.1 for rs9292777 using two additional independent Spanish samples to make a total of 4912 MS cases and 7498 controls (ORpooled = 0.84; 95%CI: 0.80-0.89; p = 1.36 × 10-9). This SNP differs from the one reported within this locus in a recent GWAS. Although it is unclear whether both signals are tapping the same genetic association, it seems clear that this locus plays an important role in the pathogenesis of MS.

Real-time molecular epidemiology of tuberculosis by direct genotyping of smear-positive clinical specimens.

We applied MIRU-VNTR (mycobacterial interspersed repetitive-unit-variable-number tandem-repeat typing) to directly analyze the bacilli present in 61 stain-positive specimens from tuberculosis patients. A complete MIRU type (24 loci) was obtained for all but one (no amplification in one locus) of the specimens (98.4%), and the allelic values fully correlated with those obtained from the corresponding cultures. Our study is the first to demonstrate that real-time genotyping of Mycobacterium tuberculosis can be achieved, fully transforming the way in which molecular epidemiology techniques can be integrated into control programs.

MicroRNA-200 family modulation in distinct breast cancer phenotypes.

The epithelial to mesenchymal transition (EMT) contributes to tumor invasion and metastasis in a variety of cancer types. In human breast cancer, gene expression studies have determined that basal-B/claudin-low and metaplastic cancers exhibit EMT-related characteristics, but the molecular mechanisms underlying this observation are unknown. As the family of miR-200 microRNAs has been shown to regulate EMT in normal tissues and cancer, here we evaluated whether the expression of the miR-200 family (miR-200f) and their epigenetic state correlate with EMT features in human breast carcinomas. We analyzed by qRT-PCR the expression of miR-200f members and various EMT-transcriptional inducers in a series of 70 breast cancers comprising an array of phenotypic subtypes: estrogen receptor positive (ER+), HER2 positive (HER2+), and triple negative (TN), including a subset of metaplastic breast carcinomas (MBCs) with sarcomatous (homologous or heterologous) differentiation. No MBCs with squamous differentiation were included. The DNA methylation status of miR-200f loci in tumor samples were inspected using Sequenom MassArray® MALDI-TOF platform. We also used two non-tumorigenic breast basal cell lines that spontaneously undergo EMT to study the modulation of miR-200f expression during EMT in vitro. We demonstrate that miR-200f is strongly decreased in MBCs compared with other cancer types. TN and HER2+ breast cancers also exhibited lower miR-200f expression than ER+ tumors. Significantly, the decreased miR-200f expression found in MBCs is accompanied by an increase in the expression levels of EMT-transcriptional inducers, and hypermethylation of the miR-200c-141 locus. Similar to tumor samples, we demonstrated that downregulation of miR-200f and hypermethylation of the miR-200c-141 locus, together with upregulation of EMT-transcriptional inducers also occur in an in vitro cellular model of spontaneous EMT. Thus, the expression and methylation status of miR-200f could be used as hypothetical biomarkers to assess the occurrence of EMT in breast cancer.

HLA Alleles Influence the Clinical Signature of Amoxicillin-Clavulanate Hepatotoxicity.

BACKGROUND AND AIM The genotype-phenotype interaction in drug-induced liver injury (DILI) is a subject of growing interest. Previous studies have linked amoxicillin-clavulanate (AC) hepatotoxicity susceptibility to specific HLA alleles. In this study we aimed to examine potential associations between HLA class I and II alleles and AC DILI with regards to phenotypic characteristics, severity and time to onset in Spanish AC hepatotoxicity cases. METHODS High resolution genotyping of HLA loci A, B, C, DRB1 and DQB1 was performed in 75 AC DILI cases and 885 controls. RESULTS The distributions of class I alleles A*3002 (P/Pc = 2.6E-6/5E-5, OR 6.7) and B*1801 (P/Pc = 0.008/0.22, OR 2.9) were more frequently found in hepatocellular injury cases compared to controls. In addition, the presence of the class II allele combination DRB1*1501-DQB1*0602 (P/Pc = 5.1E-4/0.014, OR 3.0) was significantly increased in cholestatic/mixed cases. The A*3002 and/or B*1801 carriers were found to be younger (54 vs 65 years, P = 0.019) and were more frequently hospitalized than the DRB1*1501-DQB1*0602 carriers. No additional alleles outside those associated with liver injury patterns were found to affect potential severity as measured by Hy's Law criteria. The phenotype frequencies of B*1801 (P/Pc = 0.015/0.42, OR 5.2) and DRB1*0301-DQB1*0201 (P/Pc = 0.0026/0.07, OR 15) were increased in AC DILI cases with delayed onset compared to those corresponding to patients without delayed onset, while the opposite applied to DRB1*1302-DQB1*0604 (P/Pc = 0.005/0.13, OR 0.07). CONCLUSIONS HLA class I and II alleles influence the AC DILI signature with regards to phenotypic expression, latency presentation and severity in Spanish patients.

Mutational spectrum of the SPG4 (SPAST) and SPG3A (ATL1) genes in Spanish patients with hereditary spastic paraplegia.

BACKGROUND Hereditary Spastic Paraplegias (HSP) are characterized by progressive spasticity and weakness of the lower limbs. At least 45 loci have been identified in families with autosomal dominant (AD), autosomal recessive (AR), or X-linked hereditary patterns. Mutations in the SPAST (SPG4) and ATL1 (SPG3A) genes would account for about 50% of the ADHSP cases. METHODS We defined the SPAST and ATL1 mutational spectrum in a total of 370 unrelated HSP index cases from Spain (83% with a pure phenotype). RESULTS We found 50 SPAST mutations (including two large deletions) in 54 patients and 7 ATL1 mutations in 11 patients. A total of 33 of the SPAST and 3 of the ATL1 were new mutations. A total of 141 (31%) were familial cases, and we found a higher frequency of mutation carriers among these compared to apparently sporadic cases (38% vs. 5%). Five of the SPAST mutations were predicted to affect the pre-mRNA splicing, and in 4 of them we demonstrated this effect at the cDNA level. In addition to large deletions, splicing, frameshifting, and missense mutations, we also found a nucleotide change in the stop codon that would result in a larger ORF. CONCLUSIONS In a large cohort of Spanish patients with spastic paraplegia, SPAST and ATL1 mutations were found in 15% of the cases. These mutations were more frequent in familial cases (compared to sporadic), and were associated with heterogeneous clinical manifestations.

A multicenter study confirms CD226 gene association with systemic sclerosis-related pulmonary fibrosis.

INTRODUCTION CD226 genetic variants have been associated with a number of autoimmune diseases and recently with systemic sclerosis (SSc). The aim of this study was to test the influence of CD226 loci in SSc susceptibility, clinical phenotypes and autoantibody status in a large multicenter European population. METHODS A total of seven European populations of Caucasian ancestry were included, comprising 2,131 patients with SSc and 3,966 healthy controls. Three CD226 single nucleotide polymorphisms (SNPs), rs763361, rs3479968 and rs727088, were genotyped using Taqman 5'allelic discrimination assays. RESULTS Pooled analyses showed no evidence of association of the three SNPs, neither with the global disease nor with the analyzed subphenotypes. However, haplotype block analysis revealed a significant association for the TCG haplotype (SNP order: rs763361, rs34794968, rs727088) with lung fibrosis positive patients (PBonf = 3.18E-02 OR 1.27 (1.05 to 1.54)). CONCLUSION Our data suggest that the tested genetic variants do not individually influence SSc susceptibility but a CD226 three-variant haplotype is related with genetic predisposition to SSc-related pulmonary fibrosis.

Evidence of new risk genetic factor to systemic lupus erythematosus: the UBASH3A gene

The ubiquitin associated and Src-homology 3 (SH3) domain containing A (UBASH3a) is a suppressor of T-cell receptor signaling, underscoring antigen presentation to T-cells as a critical shared mechanism of diseases pathogenesis. The aim of the present study was to determine whether the UBASH3a gene influence the susceptibility to systemic lupus erythematosus (SLE) in Caucasian populations. We evaluated five UBASH3a polymorphisms (rs2277798, rs2277800, rs9976767, rs13048049 and rs17114930), using TaqMan® allelic discrimination assays, in a discovery cohort that included 906 SLE patients and 1165 healthy controls from Spain. The SNPs that exhibit statistical significance difference were evaluated in a German replication cohort of 360 SLE patients and 379 healthy controls. The case-control analysis in the Spanish population showed a significant association between the rs9976767 and SLE (Pc = 9.9E-03 OR = 1.21 95%CI = 1.07-1.37) and a trend of association for the rs2277798 analysis (P = 0.09 OR = 0.9 95%CI = 0.79-1.02). The replication in a German cohort and the meta-analysis confirmed that the rs9976767 (Pc = 0.02; Pc = 2.4E-04, for German cohort and meta-analysis, respectively) and rs2277798 (Pc = 0.013; Pc = 4.7E-03, for German cohort and meta-analysis, respectively) UBASH3a variants are susceptibility factors for SLE. Finally, a conditional regression analysis suggested that the most likely genetic variation responsible for the association was the rs9976767 polymorphism. Our results suggest that UBASH3a gene plays a role in the susceptibility to SLE. Moreover, our study indicates that UBASH3a can be considered as a common genetic factor in autoimmune diseases.

No evidence of association between common autoimmunity STAT4 and IL23R risk polymorphisms and non-anterior uveitis.

OBJECTIVE: STAT4 and IL23R loci represent common susceptibility genetic factors in autoimmunity. We decided to investigate for the first time the possible role of different STAT4/IL23R autoimmune disease-associated polymorphisms on the susceptibility to develop non-anterior uveitis and its main clinical phenotypes. METHODS Four functional polymorphisms (rs3821236, rs7574865, rs7574070, and rs897200) located within STAT4 gene as well as three independent polymorphisms (rs7517847, rs11209026, and rs1495965) located within IL23R were genotyped using TaqMan® allelic discrimination in a total of 206 patients with non-anterior uveitis and 1553 healthy controls from Spain. RESULTS No statistically significant differences were found when allele and genotype distributions were compared between non-anterior uveitis patients and controls for any STAT4 (rs3821236: P=0.39, OR=1.12, CI 95%=0.87-1.43; rs7574865: P=0.59 OR=1.07, CI 95%=0.84-1.37; rs7574070: P=0.26, OR=0.89, CI 95%=0.72-1.10; rs897200: P=0.22, OR=0.88, CI 95%=0.71-1.08;) or IL23R polymorphisms (rs7517847: P=0.49, OR=1.08, CI 95%=0.87-1.33; rs11209026: P=0.26, OR=0.78, CI 95%=0.51-1.21; rs1495965: P=0.51, OR=0.93, CI 95%=0.76-1.15). CONCLUSION Our results do not support a relevant role, similar to that described for other autoimmune diseases, of IL23R and STAT4 polymorphisms in the non-anterior uveitis genetic predisposition. Further studies are needed to discard a possible weak effect of the studied variant.

The systemic lupus erythematosus IRF5 risk haplotype is associated with systemic sclerosis

Systemic sclerosis (SSc) is a fibrotic autoimmune disease in which the genetic component plays an important role. One of the strongest SSc association signals outside the human leukocyte antigen (HLA) region corresponds to interferon (IFN) regulatory factor 5 (IRF5), a major regulator of the type I IFN pathway. In this study we aimed to evaluate whether three different haplotypic blocks within this locus, which have been shown to alter the protein function influencing systemic lupus erythematosus (SLE) susceptibility, are involved in SSc susceptibility and clinical phenotypes. For that purpose, we genotyped one representative single-nucleotide polymorphism (SNP) of each block (rs10488631, rs2004640, and rs4728142) in a total of 3,361 SSc patients and 4,012 unaffected controls of Caucasian origin from Spain, Germany, The Netherlands, Italy and United Kingdom. A meta-analysis of the allele frequencies was performed to analyse the overall effect of these IRF5 genetic variants on SSc. Allelic combination and dependency tests were also carried out. The three SNPs showed strong associations with the global disease (rs4728142: P  = 1.34×10(-8), OR  = 1.22, CI 95%  = 1.14-1.30; rs2004640: P  = 4.60×10(-7), OR  = 0.84, CI 95%  = 0.78-0.90; rs10488631: P  = 7.53×10(-20), OR  = 1.63, CI 95%  = 1.47-1.81). However, the association of rs2004640 with SSc was not independent of rs4728142 (conditioned P  = 0.598). The haplotype containing the risk alleles (rs4728142*A-rs2004640*T-rs10488631*C: P  = 9.04×10(-22), OR  = 1.75, CI 95%  = 1.56-1.97) better explained the observed association (likelihood P-value  = 1.48×10(-4)), suggesting an additive effect of the three haplotypic blocks. No statistical significance was observed in the comparisons amongst SSc patients with and without the main clinical characteristics. Our data clearly indicate that the SLE risk haplotype also influences SSc predisposition, and that this association is not sub-phenotype-specific.

Evaluation of the IL2/IL21, IL2RA and IL2RB genetic variants influence on the endogenous non-anterior uveitis genetic predisposition.

BACKGROUND Recently, different genetic variants located within the IL2/IL21 genetic region as well as within both IL2RA and IL2RB loci have been associated to multiple autoimmune disorders. We aimed to investigate for the first time the potential influence of the IL2/IL21, IL2RA and IL2RB most associated polymorphisms with autoimmunity on the endogenous non-anterior uveitis genetic predisposition. METHODS A total of 196 patients with endogenous non-anterior uveitis and 760 healthy controls, all of them from Caucasian population, were included in the current study. The IL2/IL21 (rs2069762, rs6822844 and rs907715), IL2RA (2104286, rs11594656 and rs12722495) and IL2RB (rs743777) genetic variants were genotyped using TaqMan® allelic discrimination assays. RESULTS A statistically significant difference was found for the rs6822844 (IL2/IL21 region) minor allele frequency in the group of uveitis patients compared with controls (P(-value)=0.02, OR=0.64 CI 95%=0.43-0.94) although the significance was lost after multiple testing correction. Furthermore, no evidence of association with uveitis was detected for the analyzed genetic variants of the IL2RA or IL2RB loci. CONCLUSION Our results indicate that analyzed IL2/IL21, IL2RA and IL2RB polymorphisms do not seem to play a significant role on the non-anterior uveitis genetic predisposition although further studies are needed in order to clear up the influence of these loci on the non-anterior uveitis susceptibility.