Microarray analysis of autoimmune diseases by machine learning procedures

—Microarray-based global gene expression profiling, with the use of sophisticated statistical algorithms is providing new insights into the pathogenesis of autoimmune diseases. We have applied a novel statistical technique for gene selection based on machine learning approaches to analyze microarray expression data gathered from patients with systemic lupus erythematosus (SLE) and primary antiphospholipid syndrome (PAPS), two autoimmune diseases of unknown genetic origin that share many common features. The methodology included a combination of three data discretization policies, a consensus gene selection method, and a multivariate correlation measurement. A set of 150 genes was found to discriminate SLE and PAPS patients from healthy individuals. Statistical validations demonstrate the relevance of this gene set from an univariate and multivariate perspective. Moreover, functional characterization of these genes identified an interferon-regulated gene signature, consistent with previous reports. It also revealed the existence of other regulatory pathways, including those regulated by PTEN, TNF, and BCL-2, which are altered in SLE and PAPS. Remarkably, a significant number of these genes carry E2F binding motifs in their promoters, projecting a role for E2F in the regulation of autoimmunity.

Molecular tools to improve chestnut management: El Bierzo as a case study

The European chestnut (Castanea sativa Mill.) is a multipurpose species that has been widely cultivated around the Mediterranean basin since ancient times. New varieties were brought to the Iberian Peninsula during the Roman Empire, which coexist since then with native populations that survived the last glaciation. The relevance of chestnut cultivation has being steadily growing since the Middle Ages, until the rural decline of the past century put a stop to this trend. Forest fires and diseases were also major factors. Chestnut cultivation is gaining momentum again due to its economic (wood, fruits) and ecologic relevance, and represents currently an important asset in many rural areas of Europe. In this Thesis we apply different molecular tools to help improve current management strategies. For this study we have chosen El Bierzo (Castile and Leon, NW Spain), which has a centenary tradition of chestnut cultivation and management, and also presents several unique features from a genetic perspective (next paragraph). Moreover, its nuts are widely appreciated in Spain and abroad for their organoleptic properties. We have focused our experimental work on two major problems faced by breeders and the industry: the lack of a fine-grained genetic characterization and the need for new strategies to control blight disease. To characterize with sufficient detail the genetic diversity and structure of El Bierzo orchards, we analyzed DNA from 169 trees grafted for nut production covering the entire region. We also analyzed 62 nuts from all traditional varieties. El Bierzo constitutes an outstanding scenario to study chestnut genetics and the influence of human management because: (i) it is located at one extreme of the distribution area; (ii) it is a major glacial refuge for the native species; (iii) it has a long tradition of human management (since Roman times, at least); and (iv) its geographical setting ensures an unusual degree of genetic isolation. Thirteen microsatellite markers provided enough informativeness and discrimination power to genotype at the individual level. Together with an unexpected level of genetic variability, we found evidence of genetic structure, with three major gene pools giving rise to the current population. High levels of genetic differentiation between groups supported this organization. Interestingly, genetic structure does not match with spatial boundaries, suggesting that the exchange of material and cultivation practices have strongly influenced natural gene flow. The microsatellite markers selected for this study were also used to classify a set of 62 samples belonging to all traditional varieties. We identified several cases of synonymies and homonymies, evidencing the need to substitute traditional classification systems with new tools for genetic profiling. Management and conservation strategies should also benefit from these tools. The avenue of high-throughput sequencing technologies, combined with the development of bioinformatics tools, have paved the way to study transcriptomes without the need for a reference genome. We took advantage of RNA sequencing and de novo assembly tools to determine the transcriptional landscape of chestnut in response to blight disease. In addition, we have selected a set of candidate genes with high potential for developing resistant varieties via genetic engineering. Our results evidenced a deep transcriptional reprogramming upon fungal infection. The plant hormones ET and JA appear to orchestrate the defensive response. Interestingly, our results also suggest a role for auxins in modulating such response. Many transcription factors were identified in this work that interact with promoters of genes involved in disease resistance. Among these genes, we have conducted a functional characterization of a two major thaumatin-like proteins (TLP) that belongs to the PR5 family. Two genes encoding chestnut cotyledon TLPs have been previously characterized, termed CsTL1 and CsTL2. We substantiate here their protective role against blight disease for the first time, including in silico, in vitro and in vivo evidence. The synergy between TLPs and other antifungal proteins, particularly endo-p-1,3-glucanases, bolsters their interest for future control strategies based on biotechnological approaches.

Adenoviral-mediated gene transfer in lymphocytes

Although adenovirus can infect a wide range of cell types, lymphocytes are not generally susceptible to adenovirus infection, in part because of the absence of the expression of the cellular receptor for the adenoviral fiber protein. The cellular receptor for adenovirus and coxsackievirus (CAR) recently was cloned and shown to mediate adenoviral entry by interaction with the viral fiber protein. We show that the ectopic expression of CAR in various lymphocyte cell lines, which are almost completely resistant to adenovirus infection, is sufficient to facilitate the efficient transduction of these cells by recombinant adenoviruses. Furthermore, this property of CAR does not require its cytoplasmic domain, consistent with the idea that CAR primarily serves as a high affinity binding site for the adenoviral fiber protein, and that viral entry is mediated by interaction of the viral penton base proteins with cellular integrins. As a demonstration of their functional utility, we used CAR-expressing lymphocytes transduced with an adenovirus expressing Fas ligand to efficiently kill Fas receptor-expressing tumor cells. The ability to efficiently manipulate gene expression in lymphocyte cells by using adenovirus vectors should facilitate the functional characterization of pathways affecting lymphocyte physiology.

The neuronal RNA-binding protein Nova-2 is implicated as the autoantigen targeted in POMA patients with dementia

Paraneoplastic opsoclonus myoclonus ataxia (POMA) is a neurologic disorder thought to be mediated by an autoimmune attack against onconeural disease antigens that are expressed by gynecologic or lung tumors and by neurons. One POMA disease antigen, termed Nova-1, has been identified as a neuron-specific KH-type RNA-binding protein. Nova-1 expression is restricted to specific regions of the central nervous system, primarily the hindbrain and ventral spinal cord, which correlate with the predominantly motor symptoms in POMA. However, POMA antisera recognize antigens that are widely expressed in both caudal and rostral regions of the central nervous system, and some patients develop cognitive symptoms. We have used POMA antisera to clone a cDNA encoding a second POMA disease antigen termed Nova-2. Nova-2 is closely related to Nova-1, and is expressed at high levels in neurons during development and in adulthood, and at lower levels in the adult lung. In the postnatal mouse brain, Nova-2 is expressed in a pattern that is largely reciprocal with Nova-1, including high levels of Nova-2 expression in the neocortex and hippocampus. Functional characterization of Nova-2 in RNA selection and nitrocellulose filter-binding assays reveals that Nova-2 binds RNA with high affinity and with sequence specificity that differs from Nova-1. Our results demonstrate that the immune response in POMA targets a family of highly related sequence-specific neuronal RNA-binding proteins. The expression pattern of the Nova-2 protein is likely to underlie the development of cognitive deficits in some POMA patients.

Human γ-aminobutyric acid type B receptors are differentially expressed and regulate inwardly rectifying K+ channels

γ-Aminobutyric acid type B receptors (GABABRs) are involved in the fine tuning of inhibitory synaptic transmission. Presynaptic GABABRs inhibit neurotransmitter release by down-regulating high-voltage activated Ca2+ channels, whereas postsynaptic GABABRs decrease neuronal excitability by activating a prominent inwardly rectifying K+ (Kir) conductance that underlies the late inhibitory postsynaptic potentials. Here we report the cloning and functional characterization of two human GABABRs, hGABABR1a (hR1a) and hGABABR1b (hR1b). These receptors closely match the pharmacological properties and molecular weights of the most abundant native GABABRs. We show that in transfected mammalian cells hR1a and hR1b can modulate heteromeric Kir3.1/3.2 and Kir3.1/3.4 channels. Heterologous expression therefore supports the notion that Kir3 channels are the postsynaptic effectors of GABABRs. Our data further demonstrate that in principle either of the cloned receptors could mediate inhibitory postsynaptic potentials. We find that in the cerebellum hR1a and hR1b transcripts are largely confined to granule and Purkinje cells, respectively. This finding supports a selective association of hR1b, and not hR1a, with postsynaptic Kir3 channels. The mapping of the GABABR1 gene to human chromosome 6p21.3, in the vicinity of a susceptibility locus (EJM1) for idiopathic generalized epilepsies, identifies a candidate gene for inherited forms of epilepsy.

Wild-type Escherichia coli grows on the chitin disaccharide, N,N′-diacetylchitobiose, by expressing the cel operon

We report here that wild-type Escherichia coli can grow on the chitin disaccharide, N,N′-diacetylchitobiose (GlcNAc)2, as the sole source of carbon. Transposon mutants were isolated that were unable to ferment (GlcNAc)2 but grew normally on the monosaccharide GlcNAc. One such mutant was used to screen a wild-type E. coli genomic cosmid library for restoration of (GlcNAc)2 fermentation. A partial sequence analysis of the isolated fragment mapped the clone to the (previously sequenced) E. coli genome between 39.0 and 39.2 min. The nucleotide ORFs at this region had been previously assigned to code for a “cryptic” cellobiose utilization (cel) operon. We report here, however, that functional analysis of the operon, including growth and chemotaxis, reveal that it encodes a set of proteins that are not cryptic, but are induced by (GlcNAc)2 and catabolize the disaccharide. We therefore propose to rename the cel operon as the chb (N,N′-diacetylchitobiose) operon, with the letter designation of the genes of the operon to be reassigned consistent with the nomenclature based on functional characterization of the gene products as follows: celA to chbB, celB to chbC, celC to chbA, celD to chbR, and celF to chbF. Furthermore, sequencing evidence indicates that the operon contains an additional gene of unknown function to be designated as chbG. Thus, the overall gene sequence is to be named chbBCARFG.

Phototactic Migration of Dictyostelium Cells Is Linked to a New Type of Gelsolin-related Protein

The molecular and functional characterization of a 125-kDa Ca2+-extractable protein of the Triton X-100–insoluble fraction of Dictyostelium cells identified a new type of a gelsolin-related molecule. In addition to its five gelsolin segments, this gelsolin-related protein of 125 kDa (GRP125) reveals a number of unique domains, two of which are predicted to form coiled-coil regions. Another distinct attribute of GRP125 concerns the lack of sequence elements known to be essential for characteristic activities of gelsolin-like proteins, i.e. the severing, capping, or nucleation of actin filaments. The subcellular distribution of GRP125 to vesicular compartments suggests an activity of GRP125 different from actin-binding, gelsolin-related proteins. GRP125 expression is tightly regulated and peaks at the transition to the multicellular pseudoplasmodial stage of Dictyostelium development. GRP125 was found indispensable for slug phototaxis, because slugs fail to correctly readjust their orientation in the absence of GRP125. Analysis of the GRP125-deficient mutant showed that GRP125 is required for coupling photodetection to the locomotory machinery of slugs. We propose that GRP125 is essential in the natural environment for the propagation of Dictyostelium spores. We also present evidence for further representatives of the GRP125 type in Dictyostelium, as well as in heterologous cells from lower to higher eukaryotes.

Disruption of a Dynamin Homologue Affects Endocytosis, Organelle Morphology, and Cytokinesis in Dictyostelium discoideum

The identification and functional characterization of Dictyostelium discoideum dynamin A, a protein composed of 853 amino acids that shares up to 44% sequence identity with other dynamin-related proteins, is described. Dynamin A is present during all stages of D. discoideum development and is found predominantly in the cytosolic fraction and in association with endosomal and postlysosomal vacuoles. Overexpression of the protein has no adverse effect on the cells, whereas depletion of dynamin A by gene-targeting techniques leads to multiple and complex phenotypic changes. Cells lacking a functional copy of dymA show alterations of mitochondrial, nuclear, and endosomal morphology and a defect in fluid-phase uptake. They also become multinucleated due to a failure to complete normal cytokinesis. These pleiotropic effects of dynamin A depletion can be rescued by complementation with the cloned gene. Morphological studies using cells producing green fluorescent protein-dynamin A revealed that dynamin A associates with punctate cytoplasmic vesicles. Double labeling with vacuolin, a marker of a postlysosomal compartment in D. discoideum, showed an almost complete colocalization of vacuolin and dynamin A. Our results suggest that that dynamin A is likely to function in membrane trafficking processes along the endo-lysosomal pathway of D. discoideum but not at the plasma membrane.

A Mutation in a Novel Yeast Proteasomal Gene, RPN11/MPR1, Produces a Cell Cycle Arrest, Overreplication of Nuclear and Mitochondrial DNA, and an Altered Mitochondrial Morphology

We report here the functional characterization of an essential Saccharomyces cerevisiae gene, MPR1, coding for a regulatory proteasomal subunit for which the name Rpn11p has been proposed. For this study we made use of the mpr1-1 mutation that causes the following pleiotropic defects. At 24°C growth is delayed on glucose and impaired on glycerol, whereas no growth is seen at 36°C on either carbon source. Microscopic observation of cells growing on glucose at 24°C shows that most of them bear a large bud, whereas mitochondrial morphology is profoundly altered. A shift to the nonpermissive temperature produces aberrant elongated cell morphologies, whereas the nucleus fails to divide. Flow cytometry profiles after the shift to the nonpermissive temperature indicate overreplication of both nuclear and mitochondrial DNA. Consistently with the identification of Mpr1p with a proteasomal subunit, the mutation is complemented by the human POH1 proteasomal gene. Moreover, the mpr1-1 mutant grown to stationary phase accumulates ubiquitinated proteins. Localization of the Rpn11p/Mpr1p protein has been studied by green fluorescent protein fusion, and the fusion protein has been found to be mainly associated to cytoplasmic structures. For the first time, a proteasomal mutation has also revealed an associated mitochondrial phenotype. We actually showed, by the use of [rho°] cells derived from the mutant, that the increase in DNA content per cell is due in part to an increase in the amount of mitochondrial DNA. Moreover, microscopy of mpr1-1 cells grown on glucose showed that multiple punctate mitochondrial structures were present in place of the tubular network found in the wild-type strain. These data strongly suggest that mpr1-1 is a valuable tool with which to study the possible roles of proteasomal function in mitochondrial biogenesis.

The RD114/simian type D retrovirus receptor is a neutral amino acid transporter

The RD114/simian type D retroviruses, which include the feline endogenous retrovirus RD114, all strains of simian immunosuppressive type D retroviruses, the avian reticuloendotheliosis group including spleen necrosis virus, and baboon endogenous virus, use a common cell-surface receptor for cell entry. We have used a retroviral cDNA library approach, involving transfer and expression of cDNAs from highly infectable HeLa cells to nonpermissive NIH 3T3 mouse cells, to clone and identify this receptor. The cloned cDNA, denoted RDR, is an allele of the previously cloned neutral amino acid transporter ATB0 (SLC1A5). Both RDR and ATB0 serve as retrovirus receptors and both show specific transport of neutral amino acids. We have localized the receptor by radiation hybrid mapping to a region of about 500-kb pairs on the long arm of human chromosome 19 at q13.3. Infection of cells with RD114/type D retroviruses results in impaired amino acid transport, suggesting a mechanism for virus toxicity and immunosuppression. The identification and functional characterization of this retrovirus receptor provide insight into the retrovirus life cycle and pathogenesis and will be an important tool for optimization of gene therapy using vectors derived from RD114/type D retroviruses.

Does the colonic H,K-ATPase also act as an Na,K-ATPase?

We previously have demonstrated that the colonic P-ATPase α subunit cDNA encodes an H,K-ATPase when expressed in Xenopus laevis oocytes. Besides its high level of amino acid homology (75%) with the Na,K-ATPase, the colonic H,K-ATPase also shares a common pharmacological profile with Na,K-ATPase, because both are ouabain-sensitive and Sch 28080-insensitive. These features raise the possibility that an unrecognized property of the colonic H,K-ATPase would be Na+ translocation. To test this hypothesis, ion-selective microelectrodes were used to measure the intracellular Na+ activity of X. laevis oocytes expressing various combinations of P-ATPase subunits. The results show that expression in oocytes of the colonic H,K-ATPase affects intracellular Na+ homeostasis in a way similar to the expression of the Bufo marinus Na,K-ATPase; intracellular Na+ activity is lower in oocytes expressing the colonic H,K-ATPase or the B. marinus Na,K-ATPase than in oocytes expressing the gastric H,K-ATPase or a β subunit alone. In oocytes expressing the colonic H,K-ATPase, the decrease in intracellular Na+ activity persists when diffusive Na+ influx is enhanced by functional expression of the amiloride-sensitive epithelial Na+ channel, suggesting that the decrease is related to increased active Na+ efflux. The Na+ decrease depends on the presence of K+ in the external medium and is inhibited by 2 mM ouabain, a concentration that inhibits the colonic H,K-ATPase. These data are consistent with the hypothesis that the colonic H,K-ATPase may transport Na+, acting as an (Na,H),K-ATPase. Despite its molecular and functional characterization, the physiological role of the colonic (Na,H),K-ATPase in colonic and renal ion homeostasis remains to be elucidated.

Complementation of the Arabidopsis pds1 Mutation with the Gene Encoding p-Hydroxyphenylpyruvate Dioxygenase

Plastoquinone and tocopherols are the two major quinone compounds in higher plant chloroplasts and are synthesized by a common pathway. In previous studies we characterized two loci in Arabidopsis defining key steps of this biosynthetic pathway. Mutation of the PDS1 locus disrupts the activity of p-hydroxyphenylpyruvate dioxygenase (HPPDase), the first committed step in the synthesis of both plastoquinone and tocopherols in plants. Although plants homozygous for the pds1 mutation could be rescued by growth in the presence of homogentisic acid, the product of HPPDase, we were unable to determine if the mutation directly or indirectly disrupted HPPDase activity. This paper reports the isolation of a cDNA, pHPPD, encoding Arabidopsis HPPDase and its functional characterization by expression in both plants and Escherichia coli. pHPPD encodes a 50-kD polypeptide with homology to previously identified HPPDases, including 37 highly conserved amino acid residues clustered in the carboxyl region of the protein. Expression of pHPPD in E. coli catalyzes the accumulation of homogentisic acid, indicating that it encodes a functional HPPDase enzyme. Mapping of pHPPD and co-segregation analysis of the pds1 mutation and the HPPD gene indicate tight linkage. Constitutive expression of pHPPD in a pds1 mutant background complements this mutation. Finally, comparison of the HPPD genomic sequences from wild type and pds1 identified a 17-bp deletion in the pds1 allele that results in deletion of the carboxyterminal 26 amino acids of the HPPDase protein. Together, these data conclusively demonstrate that pds1 is a mutation in the HPPDase structural gene.

Multidrug resistance mediated by a bacterial homolog of the human multidrug transporter MDR1.

Resistance of Lactococcus lactis to cytotoxic compounds shares features with the multidrug resistance phenotype of mammalian tumor cells. Here, we report the gene cloning and functional characterization in Escherichia coli of LmrA, a lactococcal structural and functional homolog of the human multidrug resistance P-glycoprotein MDR1. LmrA is a 590-aa polypeptide that has a putative topology of six alpha-helical transmembrane segments in the N-terminal hydrophobic domain, followed by a hydrophilic domain containing the ATP-binding site. LmrA is similar to each of the two halves of MDR1 and may function as a homodimer. The sequence conservation between LmrA and MDR1 includes particular regions in the transmembrane domains and connecting loops, which, in MDR1 and the MDR1 homologs in other mammalian species, have been implicated as determinants of drug recognition and binding. LmrA and MDR1 extrude a similar spectrum of amphiphilic cationic compounds, and the activity of both systems is reversed by reserpine and verapamil. As LmrA can be functionally expressed in E. coli, it offers a useful prokaryotic model for future studies on the molecular mechanism of MDR1-like multidrug transporters.

Comparison of the defensive elicitation activity of cerato-populin and cerato-platanin: a structural model

Plants can defend themselves from potential pathogenic microorganisms relying on a complex interplay of signaling pathways: activation of the MAPK cascade, transcription of defense related genes, production of reactive oxygen species, nitric oxide and synthesis of other defensive compounds such as phytoalexins. These events are triggered by the recognition of pathogen’s effectors (effector-triggered immunity) or PAMPs (PAMP-triggered immunity). The Cerato Platanin Family (CPF) members are Cys-rich proteins secreted and localized on fungal cell walls, involved in several aspects of fungal development and pathogen-host interactions. Although more than hundred genes of the CPF have been identified and analyzed, the structural and functional characterization of the expressed proteins has been restricted only to few members of the family. Interestingly, those proteins have been shown to bind chitin with diverse affinity and after foliar treatment they elicit defensive mechanisms in host and non-host plants. This property turns cerato platanins into interesting candidates, worth to be studied to develop new fungal elicitors with applications in sustainable agriculture. This study focus on cerato-platanin (CP), core member of the family and on the orthologous cerato-populin (Pop1). The latter shows an identity of 62% and an overall homology of 73% with respect to CP. Both proteins are able to induce MAPKs phosphorylation, production of reactive oxygen species and nitric oxide, overexpression of defense’s related genes, programmed cell death and synthesis of phytoalexins. CP, however, when compared to Pop1, induces a faster response and, in some cases, a stronger activity on plane leaves. Aim of the present research is to verify if the dissimilarities observed in the defense elicitation activity of these proteins can be associated to their structural and dynamic features. Taking advantage of the available CP NMR structure, Pop1’s 3D one was obtained by homology modeling. Experimental residual dipolar couplings and 1H, 15N, 13C resonance assignments were used to validate the model. Previous works on CPF members, addressed the highly conserved random coil regions (loops b1-b2 and b2-b3) as sufficient and necessary to induce necrosis in plants’ leaves: that region was investigated in both Pop1 and CP. In the two proteins the loops differ, in their primary sequence, for few mutations and an insertion with a consequent diversification of the proteins’ electrostatic surface. A set of 2D and 3D NMR experiments was performed to characterize both the spatial arrangement and the dynamic features of the loops. NOE data revealed a more extended network of interactions between the loops in Pop1 than in CP. In addition, in Pop1 we identified a salt bridge Lys25/Asp52 and a strong hydrophobic interaction between Phe26/Trp53. These structural features were expected not only to affect the loops’ spatial arrangement, but also to reduce the degree of their conformational freedom. Relaxation data and the order parameter S2 indeed highlighted reduced flexibility, in particular for loop b1-b2 of Pop1. In vitro NMR experiments, where Pop1 and CP were titrated with oligosaccharides, supported the hypothesis that the loops structural and dynamic differences may be responsible for the different chitin-binding properties of the two proteins: CP selectively binds tetramers of chitin in a shallow groove on one side of the barrel defined by loops b1-b2, b2-b3 and b4-b5, Pop1, instead, interacts in a non-specific fashion with oligosaccharides. Because the region involved in chitin-binding is also responsible for the defense elicitation activity, possibly being recognized by plant's receptors, it is reasonable to expect that those structural and dynamic modifications may also justify the different extent of defense elicitation. To test that hypothesis, the initial steps of a protocol aimed to the identify a receptor for CP, in silico, are presented.