Comprehensive analysis of NRG1 common and rare variants in Hirschsprung patients.

Hirschsprung disease (HSCR, OMIM 142623) is a developmental disorder characterized by the absence of ganglion cells along variable lengths of the distal gastrointestinal tract, which results in tonic contraction of the aganglionic gut segment and functional intestinal obstruction. The RET proto-oncogene is the major gene for HSCR with differential contributions of its rare and common, coding and noncoding mutations to the multifactorial nature of this pathology. Many other genes have been described to be associated with the pathology, as NRG1 gene (8p12), encoding neuregulin 1, which is implicated in the development of the enteric nervous system (ENS), and seems to contribute by both common and rare variants. Here we present the results of a comprehensive analysis of the NRG1 gene in the context of the disease in a series of 207 Spanish HSCR patients, by both mutational screening of its coding sequence and evaluation of 3 common tag SNPs as low penetrance susceptibility factors, finding some potentially damaging variants which we have functionally characterized. All of them were found to be associated with a significant reduction of the normal NRG1 protein levels. The fact that those mutations analyzed alter NRG1 protein would suggest that they would be related with HSCR disease not only in Chinese but also in a Caucasian population, which reinforces the implication of NRG1 gene in this pathology.

Genetic variations in drug-induced liver injury (DILI): resolving the puzzle.

Despite stringent requirements for drug development imposed by regulatory agencies, drug-induced liver injury (DILI) is an increasing health problem and a significant cause for failure to approve drugs, market withdrawal of commercialized medications, and adoption of regulatory measures. The pathogenesis is yet undefined, though the rare occurrence of idiosyncratic DILI (1/100,000–1/10,000) and the fact that hepatotoxicity often recurs after re-exposure to the culprit drug under different environmental conditions strongly points toward a major role for genetic variations in the underlying mechanism and susceptibility. Pharmacogenetic studies in DILI have to a large extent focused on genes involved in drug metabolism, as polymorphisms in these genes may generate increased plasma drug concentrations as well as lower clearance rates when treated with standard medication doses. A range of studies have identified a number of genetic variants in drug metabolism Phase I, II, and III genes, including cytochrome P450 (CYP) 2E1, N-acetyltransferase 2, UDP-glucuronosyltransferase 2B7, glutathione S-transferase M1/T1, ABCB11, and ABCC2, that enhance DILI susceptibility (Andrade et al., 2009; Agundez et al., 2011). Several metabolic gene variants, such as CYP2E1c1 and NAT2 slow, have been associated with DILI induced by specific drugs based on individual drug metabolism information. Others, such as GSTM1 and T1 null alleles have been associated with enhanced risk of DILI development induced by a large range of drugs. Hence, these variants appear to have a more general role in DILI susceptibility due to their role in reducing the cell's antioxidative capacity (Lucena et al., 2008). Mitochondrial superoxide dismutase (SOD2) and glutathione peroxidase 1 (GPX1) are two additional enzymes involved in combating oxidative stress, with specific genetic variants shown to enhance the risk of developing DILI

Rab18 dynamics in adipocytes in relation to lipogenesis, lipolysis and obesity.

Lipid droplets (LDs) are organelles that coordinate lipid storage and mobilization, both processes being especially important in cells specialized in managing fat, the adipocytes. Proteomic analyses of LDs have consistently identified the small GTPase Rab18 as a component of the LD coat. However, the specific contribution of Rab18 to adipocyte function remains to be elucidated. Herein, we have analyzed Rab18 expression, intracellular localization and function in relation to the metabolic status of adipocytes. We show that Rab18 production increases during adipogenic differentiation of 3T3-L1 cells. In addition, our data show that insulin induces, via phosphatidylinositol 3-kinase (PI3K), the recruitment of Rab18 to the surface of LDs. Furthermore, Rab18 overexpression increased basal lipogenesis and Rab18 silencing impaired the lipogenic response to insulin, thereby suggesting that this GTPase promotes fat accumulation in adipocytes. On the other hand, studies of the β-adrenergic receptor agonist isoproterenol confirmed and extended previous evidence for the participation of Rab18 in lipolysis. Together, our data support the view that Rab18 is a common mediator of lipolysis and lipogenesis and suggests that the endoplasmic reticulum (ER) is the link that enables Rab18 action on these two processes. Finally, we describe, for the first time, the presence of Rab18 in human adipose tissue, wherein the expression of this GTPase exhibits sex- and depot-specific differences and is correlated to obesity. Taken together, these findings indicate that Rab18 is involved in insulin-mediated lipogenesis, as well as in β-adrenergic-induced lipolysis, likely facilitating interaction of LDs with ER membranes and the exchange of lipids between these compartments. A role for Rab18 in the regulation of adipocyte biology under both normal and pathological conditions is proposed.

Zinc-alpha 2-glycoprotein gene expression in adipose tissue is related with insulin resistance and lipolytic genes in morbidly obese patients.

OBJECTIVE Zinc-α(2) glycoprotein (ZAG) stimulates lipid loss by adipocytes and may be involved in the regulation of adipose tissue metabolism. However, to date no studies have been made in the most extreme of obesity. The aims of this study are to analyze ZAG expression levels in adipose tissue from morbidly obese patients, and their relationship with lipogenic and lipolytic genes and with insulin resistance (IR). METHODS mRNA expression levels of PPARγ, IRS-1, IRS-2, lipogenic and lipolytic genes and ZAG were quantified in visceral (VAT) and subcutaneous adipose tissue (SAT) of 25 nondiabetic morbidly obese patients, 11 with low IR and 14 with high IR. Plasma ZAG was also analyzed. RESULTS The morbidly obese patients with low IR had a higher VAT ZAG expression as compared with the patients with high IR (p = 0.023). In the patients with low IR, the VAT ZAG expression was greater than that in SAT (p = 0.009). ZAG expression correlated between SAT and VAT (r = 0.709, p<0.001). VAT ZAG expression was mainly predicted by insulin, HOMA-IR, plasma adiponectin and expression of adiponectin and ACSS2. SAT ZAG expression was only predicted by expression of ATGL. CONCLUSIONS ZAG could be involved in modulating lipid metabolism in adipose tissue and is associated with insulin resistance. These findings suggest that ZAG may be a useful target in obesity and related disorders, such as diabetes.

Munc18c in adipose tissue is downregulated in obesity and is associated with insulin.

OBJECTIVE Munc18c is associated with glucose metabolism and could play a relevant role in obesity. However, little is known about the regulation of Munc18c expression. We analyzed Munc18c gene expression in human visceral (VAT) and subcutaneous (SAT) adipose tissue and its relationship with obesity and insulin. MATERIALS AND METHODS We evaluated 70 subjects distributed in 12 non-obese lean subjects, 23 overweight subjects, 12 obese subjects and 23 nondiabetic morbidly obese patients (11 with low insulin resistance and 12 with high insulin resistance). RESULTS The lean, overweight and obese persons had a greater Munc18c gene expression in adipose tissue than the morbidly obese patients (p<0.001). VAT Munc18c gene expression was predicted by the body mass index (B = -0.001, p = 0.009). In SAT, no associations were found by different multiple regression analysis models. SAT Munc18c gene expression was the main determinant of the improvement in the HOMA-IR index 15 days after bariatric surgery (B = -2148.4, p = 0.038). SAT explant cultures showed that insulin produced a significant down-regulation of Munc18c gene expression (p = 0.048). This decrease was also obtained when explants were incubated with liver X receptor alpha (LXRα) agonist, either without (p = 0.038) or with insulin (p = 0.050). However, Munc18c gene expression was not affected when explants were incubated with insulin plus a sterol regulatory element-binding protein-1c (SREBP-1c) inhibitor (p = 0.504). CONCLUSIONS Munc18c gene expression in human adipose tissue is down-regulated in morbid obesity. Insulin may have an effect on the Munc18c expression, probably through LXRα and SREBP-1c.

IL2RA/CD25 gene polymorphisms: uneven association with multiple sclerosis (MS) and type 1 diabetes (T1D).

BACKGROUND: IL-2 receptor (IL2R) alpha is the specific component of the high affinity IL2R system involved in the immune response and in the control of autoimmunity. METHODS AND RESULTS: Here we perform a replication and fine mapping of the IL2RA gene region analyzing 3 SNPs previously associated with multiple sclerosis (MS) and 5 SNPs associated with type 1 diabetes (T1D) in a collection of 798 MS patients and 927 matched Caucasian controls from the south of Spain. We observed association with MS in 6 of 8 SNPs. The rs1570538, at the 3'- UTR extreme of the gene, previously reported to have a weak association with MS, is replicated here (P = 0.032). The most associated T1D SNP (rs41295061) was not associated with MS in the present study. However, the rs35285258, belonging to another independent group of SNPs associated with T1D, showed the maximal association in this study but different risk allele. We replicated the association of only one (rs2104286) of the two IL2RA SNPs identified in the recently performed genome-wide association study of MS. CONCLUSIONS: These findings confirm and extend the association of this gene with MS and reveal a genetic heterogeneity of the associated polymorphisms and risk alleles between MS and T1D suggesting different immunopathological roles of IL2RA in these two diseases.

Human genetic selection on the MTHFR 677C>T polymorphism.

BACKGROUND The prevalence of genotypes of the 677C>T polymorphism for the MTHFR gene varies among humans. In previous studies, we found changes in the genotypic frequencies of this polymorphism in populations of different ages, suggesting that this could be caused by an increase in the intake of folate and multivitamins by women during the periconceptional period. The aim was to analyze changes in the allelic frequencies of this polymorphism in a Spanish population, including samples from spontaneous abortions (SA). METHODS A total of 1305 subjects born in the 20th century were genotyped for the 677C>T polymorphism using allele specific real-time PCR with Taqman probes. A section of our population (n = 276) born in 1980-1989 was compared with fetal samples (n = 344) from SA of unknown etiology from the same period. RESULTS An increase in the frequency of the T allele (0.38 vs 0.47; p < 0.001) and of the TT genotype (0.14 vs 0.24; p < 0.001) in subjects born in the last quarter of the century was observed. In the 1980-1989 period, the results show that the frequency of the wild type genotype (CC) is about tenfold lower in the SA samples than in the controls (0.03 vs 0.33; p < 0.001) and that the frequency of the TT genotype increases in the controls (0.19 to 0.27) and in the SA samples (0.20 to 0.33 (p < 0.01)); r = 0.98. CONCLUSION Selection in favor of the T allele has been detected. This selection could be due to the increased fetal viability in early stages of embryonic development, as is deduced by the increase of mutants in both living and SA populations.

The MRC1/CD68 Ratio Is Positively Associated with Adipose Tissue Lipogenesis and with Muscle Mitochondrial Gene Expression in Humans.

BACKGROUND Alternative macrophages (M2) express the cluster differentiation (CD) 206 (MCR1) at high levels. Decreased M2 in adipose tissue is known to be associated with obesity and inflammation-related metabolic disturbances. Here we aimed to investigate MCR1 relative to CD68 (total macrophages) gene expression in association with adipogenic and mitochondrial genes, which were measured in human visceral [VWAT, n = 147] and subcutaneous adipose tissue [SWAT, n = 76] and in rectus abdominis muscle (n = 23). The effects of surgery-induced weight loss were also longitudinally evaluated (n = 6). RESULTS MCR1 and CD68 gene expression levels were similar in VWAT and SWAT. A higher proportion of CD206 relative to total CD68 was present in subjects with less body fat and lower fasting glucose concentrations. The ratio MCR1/CD68was positively associated with IRS1gene expression and with the expression of lipogenic genes such as ACACA, FASN and THRSP, even after adjusting for BMI. The ratio MCR1/CD68 in SWAT increased significantly after the surgery-induced weight loss (+44.7%; p = 0.005) in parallel to the expression of adipogenic genes. In addition, SWAT MCR1/CD68ratio was significantly associated with muscle mitochondrial gene expression (PPARGC1A, TFAM and MT-CO3). AT CD206 was confirmed by immunohistochemistry to be specific of macrophages, especially abundant in crown-like structures. CONCLUSION A decreased ratio MCR1/CD68 is linked to adipose tissue and muscle mitochondrial dysfunction at least at the level of expression of adipogenic and mitochondrial genes.

Expression of the transcription factor NFAT2 in human neutrophils: IgE-dependent, Ca2+- and calcineurin-mediated NFAT2 activation.

NFAT (nuclear factors of activated T cells) proteins constitute a family of transcription factors involved in mediating signal transduction. The presence of NFAT isoforms has been described in all cell types of the immune system, with the exception of neutrophils. In the present work we report for the first time the expression in human neutrophils of NFAT2 mRNA and protein. We also report that specific antigens were able to promote NFAT2 protein translocation to the nucleus, an effect that was mimicked by the treatment of neutrophils with anti-immunoglobulin E (anti-IgE) or anti-Fcepsilon-receptor antibodies. Antigens, anti-IgE and anti-FcepsilonRs also increased Ca2+ release and the intracellular activity of calcineurin, which was able to interact physically with NFAT2, in parallel to eliciting an enhanced NFAT2 DNA-binding activity. In addition, specific chemical inhibitors of the NFAT pathway, such as cyclosporin A and VIVIT peptide, abolished antigen and anti-IgE-induced cyclooxygenase-2 (COX2) gene upregulation and prostaglandin (PGE(2)) release, suggesting that this process is through NFAT. Our results provide evidence that NFAT2 is constitutively expressed in human neutrophils, and after IgE-dependent activation operates as a transcription factor in the modulation of genes, such as COX2, during allergic inflammation.

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.

Novel (60%) and recurrent (40%) androgen receptor gene mutations in a series of 59 patients with a 46,XY disorder of sex development.

BACKGROUND Androgen receptor (AR) gene mutations are the most frequent cause of 46,XY disorders of sex development (DSD) and are associated with a variety of phenotypes, ranging from phenotypic women [complete androgen insensitivity syndrome (CAIS)] to milder degrees of undervirilization (partial form or PAIS) or men with only infertility (mild form or MAIS). OBJECTIVE The aim of the study was to characterize the contribution of the AR gene to the molecular cause of 46,XY DSD in a series of Spanish patients. SETTING We studied a series of 133 index patients with 46,XY DSD in whom gonads were differentiated as testes, with phenotypes including varying degrees of undervirilization, and in whom the AR gene was the first candidate for a molecular analysis. METHODS The AR gene was sequenced (exons 1 to 8 with intronic flanking regions) in all patients and in family members of 61% of AR-mutated gene patients. RESULTS AR gene mutations were found in 59 individuals (44.4% of index patients), of whom 46 (78%) were CAIS and 13 (22%) PAIS. Fifty-seven different mutations were found: 21.0% located in exon 1, 15.8% in exons 2 and 3, 57.9% in exons 4-8, and 5.3% intronic. Twenty-three mutations (40.4%) had been previously described and 34 (59.6%) were novel. CONCLUSIONS AR gene mutation is the most frequent cause of 46,XY DSD, with a clearly higher frequency in the complete phenotype. Mutations spread along the whole coding sequence, including exon 1. This series shows that 60% of mutations detected during the period 2002-2009 were novel.

Differences in osteogenic and apoptotic genes between osteoporotic and osteoarthritic patients.

BACKGROUND Osteoporosis is a metabolic disorder characterized by a reduction in bone mass and deterioration in the microarchitectural structure of the bone, leading to a higher risk for spontaneous and fragility fractures.The main aim was to study the differences between human bone from osteoporotic and osteoarthritic patients about gene expression (osteogenesis and apoptosis), bone mineral density, microstructural and biomechanic parameters. METHODS We analyzed data from 12 subjects: 6 with osteoporotic hip fracture (OP) and 6 with hip osteoarthritis (OA), as the control group. All subjects underwent medical history, analytical determinations, densitometry, histomorphometric and biochemical study. The expression of 86 genes of osteogenesis and 86 genes of apoptosis was studied in pool of bone samples from patients with OP and OA by PCR array. RESULTS We observed that most of the genes of apoptosis and osteogenesis show a decrease in gene expression in the osteoporotic group in comparison with the osteoarthritic group. The histomorphometric study shows a lower bone quality in the group of patients with hip fractures compared to the osteoarthritic group. CONCLUSIONS The bone tissue of osteoporotic fracture patients is more fragile than the bone of OA patients. Our results showed an osteoporotic bone with a lower capacities for differentiation and osteoblastic activity as well as a lower rate of apoptosis than osteoarthritic bone. These results are related with structural and biochemical parameters.

Fine mapping and functional analysis of the multiple sclerosis risk gene CD6.

CD6 has recently been identified and validated as risk gene for multiple sclerosis (MS), based on the association of a single nucleotide polymorphism (SNP), rs17824933, located in intron 1. CD6 is a cell surface scavenger receptor involved in T-cell activation and proliferation, as well as in thymocyte differentiation. In this study, we performed a haptag SNP screen of the CD6 gene locus using a total of thirteen tagging SNPs, of which three were non-synonymous SNPs, and replicated the recently reported GWAS SNP rs650258 in a Spanish-Basque collection of 814 controls and 823 cases. Validation of the six most strongly associated SNPs was performed in an independent collection of 2265 MS patients and 2600 healthy controls. We identified association of haplotypes composed of two non-synonymous SNPs [rs11230563 (R225W) and rs2074225 (A257V)] in the 2(nd) SRCR domain with susceptibility to MS (P max(T) permutation = 1×10(-4)). The effect of these haplotypes on CD6 surface expression and cytokine secretion was also tested. The analysis showed significantly different CD6 expression patterns in the distinct cell subsets, i.e. - CD4(+) naïve cells, P = 0.0001; CD8(+) naïve cells, P<0.0001; CD4(+) and CD8(+) central memory cells, P = 0.01 and 0.05, respectively; and natural killer T (NKT) cells, P = 0.02; with the protective haplotype (RA) showing higher expression of CD6. However, no significant changes were observed in natural killer (NK) cells, effector memory and terminally differentiated effector memory T cells. Our findings reveal that this new MS-associated CD6 risk haplotype significantly modifies expression of CD6 on CD4(+) and CD8(+) T cells.

Virulence profiles of bacteremic extended-spectrum β-lactamase-producing Escherichia coli: association with epidemiological and clinical features.

There is scarce data about the importance of phylogroups and virulence factors (VF) in bloodstream infections (BSI) caused by extended-spectrum β-lactamase-producing Escherichia coli (ESBLEC). A prospective multicenter Spanish cohort including 191 cases of BSI due to ESBLEC was studied. Phylogroups and 25 VF genes were investigated by PCR. ESBLEC were classified into clusters according to their virulence profiles. The association of phylogropus, VF, and clusters with epidemiological features were studied using multivariate analysis. Overall, 57.6%, 26.7%, and 15.7% of isolates belonged to A/B1, D and B2 phylogroups, respectively. By multivariate analysis (adjusted OR [95% CI]), virulence cluster C2 was independently associated with urinary tract source (5.05 [0.96-25.48]); cluster C4 with sources other than urinary of biliary tract (2.89 [1.05-7.93]), and cluster C5 with BSI in non-predisposed patients (2.80 [0.99-7.93]). Isolates producing CTX-M-9 group ESBLs and from phylogroup D predominated among cluster C2 and C5, while CTX-M-1 group of ESBL and phylogroup B2 predominantes among C4 isolates. These results suggest that host factors and previous antimicrobial use were more important than phylogroup or specific VF in the occurrence of BSI due to ESBLEC. However, some associations between virulence clusters and some specific epidemiological features were found.

Candidate gene study of TRAIL and TRAIL receptors: association with response to interferon beta therapy in multiple sclerosis patients.

TRAIL and TRAIL Receptor genes have been implicated in Multiple Sclerosis pathology as well as in the response to IFN beta therapy. The objective of our study was to evaluate the association of these genes in relation to the age at disease onset (AAO) and to the clinical response upon IFN beta treatment in Spanish MS patients. We carried out a candidate gene study of TRAIL, TRAILR-1, TRAILR-2, TRAILR-3 and TRAILR-4 genes. A total of 54 SNPs were analysed in 509 MS patients under IFN beta treatment, and an additional cohort of 226 MS patients was used to validate the results. Associations of rs1047275 in TRAILR-2 and rs7011559 in TRAILR-4 genes with AAO under an additive model did not withstand Bonferroni correction. In contrast, patients with the TRAILR-1 rs20576-CC genotype showed a better clinical response to IFN beta therapy compared with patients carrying the A-allele (recessive model: p = 8.88×10(-4), pc = 0.048, OR = 0.30). This SNP resulted in a non synonymous substitution of Glutamic acid to Alanine in position 228 (E228A), a change previously associated with susceptibility to different cancer types and risk of metastases, suggesting a lack of functionality of TRAILR-1. In order to unravel how this amino acid change in TRAILR-1 would affect to death signal, we performed a molecular modelling with both alleles. Neither TRAIL binding sites in the receptor nor the expression levels of TRAILR-1 in peripheral blood mononuclear cell subsets (monocytes, CD4+ and CD8+ T cells) were modified, suggesting that this SNP may be altering the death signal by some other mechanism. These findings show a role for TRAILR-1 gene variations in the clinical outcome of IFN beta therapy that might have relevance as a biomarker to predict the response to IFN beta in MS.