JNK controls the onset of mitosis of planarian stem cells and triggers apoptotic cell death required for regeneration and remodeling

Regeneration of lost tissues depends on the precise interpretation of molecular signals that control and coordinate the onset of proliferation, cellular differentiation and cell death. However, the nature of those molecular signals and the mechanisms that integrate the cellular responses remain largely unknown. The planarian flatworm is a unique model in which regeneration and tissue renewal can be comprehensively studied in vivo. The presence of a population of adult pluripotent stem cells combined with the ability to decode signaling after wounding enable planarians to regenerate a complete, correctly proportioned animal within a few days after any kind of amputation, and to adapt their size to nutritional changes without compromising functionality. Here, we demonstrate that the stress-activated c-jun-NH2-kinase (JNK) links wound-induced apoptosis to the stem cell response during planarian regeneration. We show that JNK modulates the expression of wound-related genes, triggers apoptosis and attenuates the onset of mitosis in stem cells specifically after tissue loss. Furthermore, in pre-existing body regions, JNK activity is required to establish a positive balance between cell death and stem cell proliferation to enable tissue renewal, remodeling and the maintenance of proportionality. During homeostatic degrowth, JNK RNAi blocks apoptosis, resulting in impaired organ remodeling and rescaling. Our findings indicate that JNK-dependent apoptotic cell death is crucial to coordinate tissue renewal and remodeling required to regenerate and to maintain a correctly proportioned animal. Hence, JNK might act as a hub, translating wound signals into apoptotic cell death, controlled stem cell proliferation and differentiation, all of which are required to coordinate regeneration and tissue renewal.

Consequences of lineage-specific gene loss on functional evolution of surviving paralogs: ALDH1A and retinoic acid signaling in vertebrate genomes

Genome duplications increase genetic diversity and may facilitate the evolution of gene subfunctions. Little attention, however, has focused on the evolutionary impact of lineage-specific gene loss. Here, we show that identifying lineage-specific gene loss after genome duplication is important for understanding the evolution of gene subfunctions in surviving paralogs and for improving functional connectivity among human and model organism genomes. We examine the general principles of gene loss following duplication, coupled with expression analysis of the retinaldehyde dehydrogenase Aldh1a gene family during retinoic acid signaling in eye development as a case study. Humans have three ALDH1A genes, but teleosts have just one or two. We used comparative genomics and conserved syntenies to identify loss of ohnologs (paralogs derived from genome duplication) and to clarify uncertain phylogenies. Analysis showed that Aldh1a1 and Aldh1a2 form a clade that is sister to Aldh1a3-related genes. Genome comparisons showed secondarily loss of aldh1a1 in teleosts, revealing that Aldh1a1 is not a tetrapod innovation and that aldh1a3 was recently lost in medaka, making it the first known vertebrate with a single aldh1a gene. Interestingly, results revealed asymmetric distribution of surviving ohnologs between co-orthologous teleost chromosome segments, suggesting that local genome architecture can influence ohnolog survival. We propose a model that reconstructs the chromosomal history of the Aldh1a family in the ancestral vertebrate genome, coupled with the evolution of gene functions in surviving Aldh1a ohnologs after R1, R2, and R3 genome duplications. Results provide evidence for early subfunctionalization and late subfunction-partitioning and suggest a mechanistic model based on altered regulation leading to heterochronic gene expression to explain the acquisition or modification of subfunctions by surviving ohnologs that preserve unaltered ancestral developmental programs in the face of gene loss.

Erwinia amylovora Novel Plasmid pEI70: Complete Sequence, Biogeography, and Role in Aggressiveness in the Fire Blight Phytopathogen

Comparative genomics of several strains of Erwinia amylovora, a plant pathogenic bacterium causal agent of fire blight disease, revealed that its diversity is primarily attributable to the flexible genome comprised of plasmids. We recently identified and sequenced in full a novel 65.8 kb plasmid, called pEI70. Annotation revealed a lack of known virulence-related genes, but found evidence for a unique integrative conjugative element related to that of other plant and human pathogens. Comparative analyses using BLASTN showed that pEI70 is almost entirely included in plasmid pEB102 from E. billingiae, an epiphytic Erwinia of pome fruits, with sequence identities superior to 98%. A duplex PCR assay was developed to survey the prevalence of plasmid pEI70 and also that of pEA29, which had previously been described in several E. amylovora strains. Plasmid pEI70 was found widely dispersed across Europe with frequencies of 5–92%, but it was absent in E. amylovora analyzed populations from outside of Europe. Restriction analysis and hybridization demonstrated that this plasmid was identical in at least 13 strains. Curing E. amylovora strains of pEI70 reduced their aggressiveness on pear, and introducing pEI70 into low-aggressiveness strains lacking this plasmid increased symptoms development in this host. Discovery of this novel plasmid offers new insights into the biogeography, evolution and virulence determinants in E. amylovora

Bioinformática: interfaz web para estudiar el efecto de diferentes condiciones sobre la expresión de los genes

El trabajo realizado se divide en dos bloques bien diferenciados, ambos relacionados con el análisis de microarrays. El primer bloque consiste en agrupar las condiciones muestrales de todos los genes en grupos o clústers. Estas agrupaciones se obtienen al aplicar directamente sobre la microarray los siguientes algoritmos de agrupación: SOM,PAM,SOTA,HC y al aplicar sobre la microarray escalada con PC y MDS los siguientes algoritmos: SOM,PAM,SOTA,HC y K-MEANS. El segundo bloque consiste en realizar una búsqueda de genes basada en los intervalos de confianza de cada clúster de la agrupación activa. Las condiciones de búsqueda ajustadas por el usuario se validan para cada clúster respecto el valor basal 0 y respecto el resto de clústers, para estas validaciones se usan los intervalos de confianza. Estos dos bloques se integran en una aplicación web ya existente, el applet PCOPGene, alojada en el servidor: http://revolutionresearch.uab.es.

Teleost fish larvae adapt to dietary arachidonic acid supply through modulation of the expression of lipid metabolism and stress response genes

Dietary fatty acid supply can affect stress response in fish during early development. Although knowledge on the mechanisms involved in fatty acid regulation of stress tolerance is scarce, it has often been hypothesised that eicosanoid profiles can influence cortisol production. Genomic cortisol actions are mediated by cytosolic receptors which may respond to cellular fatty acid signalling. An experiment was designed to test the effects of feeding gilthead sea-bream larvae with four microdiets, containing graded arachidonic acid (ARA) levels (0·4, 0·8, 1·5 and 3·0 %), on the expression of genes involved in stress response (steroidogenic acute regulatory protein, glucocorticoid receptor and phosphoenolpyruvate carboxykinase), lipid and, particularly, eicosanoid metabolism (hormone-sensitive lipase, PPARα, phospholipase A2, cyclo-oxygenase-2 and 5-lipoxygenase), as determined by real-time quantitative PCR. Fish fatty acid phenotypes reflected dietary fatty acid profiles. Growth performance, survival after acute stress and similar whole-body basal cortisol levels suggested that sea-bream larvae could tolerate a wide range of dietary ARA levels. Transcription of all genes analysed was significantly reduced at dietary ARA levels above 0·4 %. Nonetheless, despite practical suppression of phospholipase A2 transcription, higher leukotriene B4 levels were detected in larvae fed 3·0 % ARA, whereas a similar trend was observed regarding PGE2 production. The present study demonstrates that adaptation to a wide range of dietary ARA levels in gilthead sea-bream larvae involves the modulation of the expression of genes related to eicosanoid synthesis, lipid metabolism and stress response. The roles of ARA, other polyunsaturates and eicosanoids as signals in this process are discussed.

Genomics of wine yeast: functional analysis of new and variable genes

Projecte de recerca elaborat a partir d’una estada a l’Institut National de la Recherche Agronomique, França, entre 2007 i 2009. Saccharomyces cerevisiae ha estat el llevat utilitzat durant mil.lenis en l'elaboració de vins. Tot i així, es té poc coneixement sobre les pressions de selecció que han actuat en la modelització del genoma dels llevats vínics. S’ha seqüenciat el genoma d'una soca vínica comercial, EC1118, obtenint 31 supercontigs que cobreixen el 97% del genoma de la soca de referència, S288c. S’ha trobat que el genoma de la soca vínica es diferencia bàsicament en la possessió de 3 regions úniques que contenen 34 gens implicats en funcions claus per al procés fermentatiu. A banda, s’han dut a terme estudis de filogènia i synteny (ordre dels gens) que mostren que una d'aquestes tres regions és pròxima a una espècie relacionada amb el gènere Saccharomyces, mentre que les altres dos regions tenen un origen no-Saccharomyces. S’ha identificat mitjançant PCR i seqüenciació a Zygosaccharomyces bailii, una espècie contaminant de les fermentacions víniques, com a espècie donadora d'una de les dues regions. Les hibridacions naturals entre soques de diferents espècies dins del grup Saccharomyces sensu stricto ja han estat descrites. El treball és el primer que presenta hibridacions entre espècies Saccharomyces i no-Saccharomyces (Z. bailii, en aquest cas). També s’assenyala que les noves regions es troben freqüent i diferencialment presents entre els clades de S. cerevisiae, trobant-se de manera gairebé exclusiva en el grup de les soques víniques, suggerint que es tracta d'una adquisició recent de transferència gènica. En general, les dades demostren que el genoma de les soques víniques pateix una constant remodelació mitjançant l'adquisició de gens exògens. Els resultats suggereixen que aquests processos estan afavorits per la proximitat ecològica i estan implicats en l'adaptació molecular de les soques víniques a les condicions d'elevada concentració en sucres, poc nitrogen i elevades concentracions en etanol.

I. Identification and characterisation of epigenetically altered tumor supressor genes by proteomic and genomic approches / II. Studies of genes regualted by MeCP2 union to its promoter regions by proteomic and genmomic approches

DNA methylation has an important impact on normal cell physiology, thus any defects in this mechanism may be related to the development of various diseases In this project we are interested in identifying epigeneticaliy modified genes, in general controlled by processes related to the DNA methylation, by means of a new strategy combining protomic and genomic analyses. First, the two Dimensional-Difference Gel Electrophoresis (2-DIGE) protein analyses of extracts obtained from HCT-116 wt and double knockout for DNMT1 and DNMT3b (DKO) cells revealed 34 proteins overexpressed in the condition of DNMTs depletion. From five genes with higher transcript lavels in DKO cells, comparing with HCT-116 wt. oniy AKR1B1, UCHLl and VIM are melhylated in HCT-116. As expected. the DNA methvlation 1s lost in DKO cells. The rneth,vl ation of VIM and UCHLl promoters in some cancer samples has already been repaired, thus further studies has been focused on AKRlBI. AKR1B1 expression due lo DNA methyiaton of promoter region seems to occur specilfically in the colon cancer cell Iines. which was confirmed in the DNA rnethylation status and expression analyses. performed on 32 different cancer cell lines (including colon, breast, lymphoma, leukemia, neuroblastoma, glioma and lung cancer cell Iines) as well as normal colon and normal lymphocytes samples. AKRIBI expression after treatments with DNA demethvlating agent (AZA) was rescued in 5 coloncancer cell lines (including genetic regulation of the candidate gene. The methylation status of the rest of the genes identified in proteomic analysis was checked by methylation specific PCR (MSP) experiment and all appeared to be unmethylated. The similar research has been done also bv means of Mecp2-null mouse model For 14 selected candidate genes the analyses of expression leveis, methylation Status and MeCP2 interaction with promoters are currently being performed.

Decay of linkage disequilibrium within genes across HGDP-CEPH human samples: most population isolates do not show increased LD

Background: It is well known that the pattern of linkage disequilibrium varies between human populations, with remarkable geographical stratification. Indirect association studies routinely exploit linkage disequilibrium around genes, particularly in isolated populations where it is assumed to be higher. Here, we explore both the amount and the decay of linkage disequilibrium with physical distance along 211 gene regions, most of them related to complex diseases, across 39 HGDP-CEPH population samples, focusing particularly on the populations defined as isolates. Within each gene region and population we use r2 between all possible single nucleotide polymorphism (SNP) pairs as a measure of linkage disequilibrium and focus on the proportion of SNP pairs with r2 greater than 0.8.Results: Although the average r2 was found to be significantly different both between and within continental regions, a much higher proportion of r2 variance could be attributed to differences between continental regions (2.8% vs. 0.5%, respectively). Similarly, while the proportion of SNP pairs with r2 > 0.8 was significantly different across continents for all distance classes, it was generally much more homogenous within continents, except in the case of Africa and the Americas. The only isolated populations with consistently higher LD in all distance classes with respect to their continent are the Kalash (Central South Asia) and the Surui (America). Moreover, isolated populations showed only slightly higher proportions of SNP pairs with r2 > 0.8 per gene region than non-isolated populations in the same continent. Thus, the number of SNPs in isolated populations that need to be genotyped may be only slightly less than in non-isolates. Conclusion: The "isolated population" label by itself does not guarantee a greater genotyping efficiency in association studies, and properties other than increased linkage disequilibrium may make these populations interesting in genetic epidemiology.

OSIRISv1.2: a named entity recognition system for sequence variants of genes in biomedical literature

Background: Single Nucleotide Polymorphisms, among other type of sequence variants, constitute key elements in genetic epidemiology and pharmacogenomics. While sequence data about genetic variation is found at databases such as dbSNP, clues about the functional and phenotypic consequences of the variations are generally found in biomedical literature. The identification of the relevant documents and the extraction of the information from them are hampered by the large size of literature databases and the lack of widely accepted standard notation for biomedical entities. Thus, automatic systems for the identification of citations of allelic variants of genes in biomedical texts are required. Results: Our group has previously reported the development of OSIRIS, a system aimed at the retrieval of literature about allelic variants of genes http://ibi.imim.es/osirisform.html. Here we describe the development of a new version of OSIRIS (OSIRISv1.2, http://ibi.imim.es/OSIRISv1.2.html webcite) which incorporates a new entity recognition module and is built on top of a local mirror of the MEDLINE collection and HgenetInfoDB: a database that collects data on human gene sequence variations. The new entity recognition module is based on a pattern-based search algorithm for the identification of variation terms in the texts and their mapping to dbSNP identifiers. The performance of OSIRISv1.2 was evaluated on a manually annotated corpus, resulting in 99% precision, 82% recall, and an F-score of 0.89. As an example, the application of the system for collecting literature citations for the allelic variants of genes related to the diseases intracranial aneurysm and breast cancer is presented. Conclusion: OSIRISv1.2 can be used to link literature references to dbSNP database entries with high accuracy, and therefore is suitable for collecting current knowledge on gene sequence variations and supporting the functional annotation of variation databases. The application of OSIRISv1.2 in combination with controlled vocabularies like MeSH provides a way to identify associations of biomedical interest, such as those that relate SNPs with diseases.

A peroxisomal glutathione transferase of Saccharomyces cerevisiae is functionally related to sulfur amino acid metabolism

Saccharomyces cerevisiae cells contain three omega-class glutathione transferases with glutaredoxin activity (Gto1, Gto2, and Gto3), in addition to two glutathione transferases (Gtt1 and Gtt2) not classifiable into standard classes. Gto1 is located at the peroxisomes, where it is targeted through a PTS1-type sequence, whereas Gto2 and Gto3 are in the cytosol. Among the GTO genes, GTO2 shows the strongest induction of expression by agents such as diamide, 1-chloro-2,4-dinitrobenzene, tert-butyl hydroperoxide or cadmium, in a manner that is dependent on transcriptional factors Yap1 and/or Msn2/4. Diamide and 1-chloro-2,4-dinitrobenzene (causing depletion of reduced glutathione) also induce expression of GTO1 over basal levels. Phenotypic analyses with single and multiple mutants in the S. cerevisiae glutathione transferase genes show that, in the absence of Gto1 and the two Gtt proteins, cells display increased sensitivity to cadmium. A gto1-null mutant also shows growth defects on oleic acid-based medium, which is indicative of abnormal peroxisomal functions, and altered expression of genes related to sulfur amino acid metabolism. As a consequence, growth of the gto1 mutant is delayed in growth medium without lysine, serine, or threonine, and the mutant cells have low levels of reduced glutathione. The role of Gto1 at the S. cerevisiae peroxisomes could be related to the redox regulation of the Str3 cystathionine -lyase protein. This protein is also located at the peroxisomes in S. cerevisiae, where it is involved in transulfuration of cysteine into homocysteine, and requires a conserved cysteine residue for its biological activity.

Uncovering the functional constraints underlying the genomic organisation of the Odorant-Binding Protein genes

Animal olfactory systems have a critical role for the survival and reproduction of individuals. In insects, the odorant-binding proteins (OBPs) are encoded by a moderately sized gene family, and mediate the first steps of the olfactory processing. Most OBPs are organized in clusters of a few paralogs, which are conserved over time. Currently, the biological mechanism explaining the close physical proximity among OBPs is not yet established. Here, we conducted a comprehensive study aiming to gain insights into the mechanisms underlying the OBP genomic organization. We found that the OBP clusters are embedded within large conserved arrangements. These organizations also include other non-OBP genes, which often encode proteins integral to plasma membrane. Moreover, the conservation degree of such large clusters is related to the following: 1) the promoter architecture of the confined genes, 2) a characteristic transcriptional environment, and 3) the chromatin conformation of the chromosomal region. Our results suggest that chromatin domains may restrict the location of OBP genes to regions having the appropriate transcriptional environment, leading to the OBP cluster structure. However, the appropriate transcriptional environment for OBP and the other neighbor genes is not dominated by reduced levels of expression noise. Indeed, the stochastic fluctuations in the OBP transcript abundance may have a critical role in the combinatorial nature of the olfactory coding process.

Planarians as a model to assess in vivo the role of matrix metalloproteinase genes during homeostasis and regeneration

Matrix metalloproteinases (MMPs) are major executors of extracellular matrix remodeling and, consequently, play key roles in the response of cells to their microenvironment. The experimentally accessible stem cell population and the robust regenerative capabilities of planarians offer an ideal model to study how modulation of the proteolytic system in the extracellular environment affects cell behavior in vivo. Genome-wide identification of Schmidtea mediterranea MMPs reveals that planarians possess four mmp-like genes. Two of them (mmp1 and mmp2) are strongly expressed in a subset of secretory cells and encode putative matrilysins. The other genes (mt-mmpA and mt-mmpB) are widely expressed in postmitotic cells and appear structurally related to membrane-type MMPs. These genes are conserved in the planarian Dugesia japonica. Here we explore the role of the planarian mmp genes by RNA interference (RNAi) during tissue homeostasis and regeneration. Our analyses identify essential functions for two of them. Following inhibition of mmp1 planarians display dramatic disruption of tissues architecture and significant decrease in cell death. These results suggest that mmp1 controls tissue turnover, modulating survival of postmitotic cells. Unexpectedly, the ability to regenerate is unaffected by mmp1(RNAi). Silencing of mt-mmpA alters tissue integrity and delays blastema growth, without affecting proliferation of stem cells. Our data support the possibility that the activity of this protease modulates cell migration and regulates anoikis, with a consequent pivotal role in tissue homeostasis and regeneration. Our data provide evidence of the involvement of specific MMPs in tissue homeostasis and regeneration and demonstrate that the behavior of planarian stem cells is critically dependent on the microenvironment surrounding these cells. Studying MMPs function in the planarian model provides evidence on how individual proteases work in vivo in adult tissues. These results have high potential to generate significant information for development of regenerative and anti cancer therapies.

Decay of linkage disequilibrium within genes across HGDP-CEPH human samples: most population isolates do not show increased LD

Background It is well known that the pattern of linkage disequilibrium varies between human populations, with remarkable geographical stratification. Indirect association studies routinely exploit linkage disequilibrium around genes, particularly in isolated populations where it is assumed to be higher. Here, we explore both the amount and the decay of linkage disequilibrium with physical distance along 211 gene regions, most of them related to complex diseases, across 39 HGDP-CEPH population samples, focusing particularly on the populations defined as isolates. Within each gene region and population we use r2 between all possible single nucleotide polymorphism (SNP) pairs as a measure of linkage disequilibrium and focus on the proportion of SNP pairs with r2 greater than 0.8. Results Although the average r2 was found to be significantly different both between and within continental regions, a much higher proportion of r2 variance could be attributed to differences between continental regions (2.8% vs. 0.5%, respectively). Similarly, while the proportion of SNP pairs with r2 > 0.8 was significantly different across continents for all distance classes, it was generally much more homogenous within continents, except in the case of Africa and the Americas. The only isolated populations with consistently higher LD in all distance classes with respect to their continent are the Kalash (Central South Asia) and the Surui (America). Moreover, isolated populations showed only slightly higher proportions of SNP pairs with r2 > 0.8 per gene region than non-isolated populations in the same continent. Thus, the number of SNPs in isolated populations that need to be genotyped may be only slightly less than in non-isolates. Conclusion The 'isolated population' label by itself does not guarantee a greater genotyping efficiency in association studies, and properties other than increased linkage disequilibrium may make these populations interesting in genetic epidemiology.

Analysis of a Plant Complex Resistance Gene Locus Underlying Immune-Related Hybrid Incompatibility and Its Occurrence in Nature

Mechanisms underlying speciation in plants include detrimental (incompatible) genetic interactions between parental alleles that incur a fitness cost in hybrids. We reported on recessive hybrid incompatibility between an Arabidopsis thaliana strain from Poland, Landsberg erecta (Ler), and many Central Asian A. thaliana strains. The incompatible interaction is determined by a polymorphic cluster of Toll/interleukin-1 receptor-nucleotide binding-leucine rich repeat (TNL) RPP1 (Recognition of Peronospora parasitica1)-like genes in Ler and alleles of the receptor-like kinase Strubbelig Receptor Family 3 (SRF3) in Central Asian strains Kas-2 or Kond, causing temperature-dependent autoimmunity and loss of growth and reproductive fitness. Here, we genetically dissected the RPP1-like Ler locus to determine contributions of individual RPP1-like Ler (R1R8) genes to the incompatibility. In a neutral background, expression of most RPP1-like Ler genes, except R3, has no effect on growth or pathogen resistance. Incompatibility involves increased R3 expression and engineered R3 overexpression in a neutral background induces dwarfism and sterility. However, no individual RPP1-like Ler gene is sufficient for incompatibility between Ler and Kas-2 or Kond, suggesting that co-action of at least two RPP1-like members underlies this epistatic interaction. We find that the RPP1-like Ler haplotype is frequent and occurs with other Ler RPP1-like alleles in a local population in Gorzów Wielkopolski (Poland). Only Gorzów individuals carrying the RPP1-like Ler haplotype are incompatible with Kas-2 and Kond, whereas other RPP1-like alleles in the population are compatible. Therefore, the RPP1-like Ler haplotype has been maintained in genetically different individuals at a single site, allowing exploration of forces shaping the evolution of RPP1-like genes at local and regional population scales.