Degradación bacteriana de cianuro y compuestos nitrogenados tóxicos

La ruta de asimilación de cianuro en P. pseudoalcaligenes CECT5344 transcurre a través de un nitrilo formado por la reacción química del cianuro con el oxalacetato, siendo este último acumulado como consecuencia de la acción conjunta de una malato:quinona oxidoreductasa (MQO) y la oxidasa terminal resistente a cianuro (CioAB) (Luque-Almagro et al., 2011b). Los nitrilos pueden ser convertidos en amonio por la acción de una nitrilasa o un sistema nitrilo hidratasa/amidasa. Con el objetivo de elucidar la ruta de asimilación de cianuro en P. pseudoalcalígenes CECT5344, se ha analizado el proteoma de este microorganismo en condiciones cianotróficas frente a nitrato como fuente de nitrógeno como control. En este estudio se identificaron proteínas relacionadas con la ruta de asimilación de cianuro en la estirpe CECT5344, que aparecían inducidas por cianuro, como NitB y NitG, cuyos genes se encuentran localizados en la agrupación génica nit1C. Además de NitB y NitG, de función desconocida, la agrupación génica nit1C codifica un regulador transcripcional del tipo Fis dependiente de σ54 (NitA), una nitrilasa (NitC), una proteína que pertenece a la superfamilia S-adenosilmetionina (NitD), un miembro de la superfamilia N-aciltransferasa (NitE), un polipéptido de la familia AIRS/GARS (NitF) y una oxidorreductasa dependiente de NADH (NitH). Un análisis transcripcional mediante RT-PCR determinó que los genes nitBCDEFGH se cotranscriben, mientras que el gen regulador nitA se transcribe de forma divergente. Además, resultados obtenidos por RT-PCR confirman que la expresión de los genes nitBCDEFGH está inducida por cianuro y reprimida por amonio. La relación entre el cianuro y el grupo de genes nit1C queda patente por el fenotipo de los mutantes deficientes nitA, nitB y nitC, incapaces de usar complejos cianuro-metálicos o 2-hidroxinitrilos como única fuente de nitrógeno. Todos estos datos indican que la nitrilasa NitC, junto con la proteína NitB, utilizan de forma específica determinados nitrilos alifáticos como sustrato, entre los que se encuentran el formado durante la asimilación de cianuro (Estepa et al., 2012). Además, entre las proteínas inducidas por cianuro se identificaron una dihidropicolinato sintasa (DapA), una fosfoserina transaminasa (SerC) y una proteína de función desconocida (Orf1), las tres codificadas por genes del operón cio, una cianasa (CynS), la proteína S6 de la subunidad ribosomal 30S (RpsF), una superóxido dismutasa (SodB), la ferritina (Dps), una oxidorreductasa (Fpr) y un factor de elongación P (EF-P). Una vez identificadas, estas proteínas se han analizado funcionalmente y se han localizado en el genoma de P. pseudoalcaligenes CECT5344 los genes correspondientes, así como los genes adyacentes. La inducción de estas proteínas en condiciones cianotróficas sugiere que el metabolismo del cianuro incluye, además de la resistencia y asimilación de este tóxico, otros procesos biológicos relacionados con el metabolismo del cianato y de algunos aminoácidos, el estrés oxidativo y la homeostasis de hierro, entre otros. Por otra parte, el conocimiento en profundidad y la interpretación de la secuencia génica de P. pseudoalcaligenes CECT5344, así como el análisis comparativo frente a organismos no cianotrofos ha permitido entender algunos de los mecanismos implicados en la resistencia y asimilación de cianuro, lo que permitiría conducir a la posterior mejora del proceso de biodegradación de cianuro. Además, el estudio del genoma de la estirpe CECT5344 permitirá explorar la capacidad de este organismo para ser utilizado en procesos de biorremediación de residuos cianurados en los que se encuentran metales y otros tóxicos (Luque-Almagro et al., 2013; Wibberg et al., 2014). En este trabajo se muestran y discuten los resultados de la secuenciación del genoma de P. pseudoalcaligenes, así como el estudio del análisis filogenético y evolutivo de la cepa, estableciéndose de esta manera relaciones con otras especies en base a los genomas secuenciados de las mismas, entre las que destaca P. mendocina ymp relacionada con P. pseudoalcaligenes CECT5344. El estudio de las características del genoma de P. pseudoalcaligenes CECT5344 ha sido completado con un análisis comparativo frente a los genomas de otras especies de Pseudomonas, encontrándose así semejanzas y diferencias en cuanto a la distribución génica funcional. Por último, se muestra un análisis del genoma de P. pseudoalcaligenes CECT5344 en relación con los genes implicados probablemente en los procesos de asimilación de cianuro y residuos cianurados, tales como los codificantes de nitrilasas y aquellos implicados en la resistencia a cianuro como los constituyentes del operón cio que codifican la oxidasa terminal insensible a cianuro. Finalmente, se discute la presencia de genes implicados posiblemente en otros procesos con una alto potencial biotecnológico, tales como la producción de bioplásticos y la biodegradación de diversos contaminantes.

Insights into Neural Crest Evolution

Neural crest cells are unique to vertebrates and essential to the development and evolution of the craniofacial skeleton. Using a combination of DiI cell lineage tracing, transcriptomics, and analysis of key transcription factors of the Sox Family, I examined neural crest development in the sea lamprey, Petromyzon marinus, as the most basal extant vertebrate from which it is possible to get embryos. The results have uncovered distinct cranial and trunk neural crest subpopulations along the anterior-posterior axis of the lamprey embryo, with a clear separation between the two. However, no evidence of the presence of an intermediate vagal neural crest population was uncovered. Comparing cranial neural crest genes between lamprey and chick, either by examining individual candidate genes or whole genome transcriptome analysis, reveals significant changes in the cranial neural crest gene regulatory network of lamprey compared with chick. In particular, the lamprey cranial neural crest is "missing" several gnathostome cranial crest genes. We speculate that these may underlie the evolutionary divergence of craniofacial development between jawed and jawless vertebrates. Despite the absence of vagal neural crest, DiI-labeling shows that trunk neural crest-derived cells, likely homologous to mammalian Schwann cell precursors, contribute to the lamprey enteric nervous system, potentially representing the most primitive form of neural crest cells contribution to the ENS. Finally, I characterized key members of the Sox Family (Sox B-F) due to their importance in neural crest specification in other species. In comparative studies of the SoxC genes (Sox4, Sox11, and Sox12) in both lamprey and Xenopus, I found similar expression patterns and a novel key role in early neural crest specification, suggesting a conserved role of the SoxC genes amongst vertebrates. Taken together, this work represents important progress in characterizing the early evolution of the neural crest in vertebrates and its role in the transition from jawless to jawed vertebrates.

Hormone-Mediated Gene-Specific Translation Regulation

The central role of translation regulation in the control of critical cellular processes has long been recognized. Yet the systematic exploration of quantitative changes in translation at a genome-wide scale in response to specific stimuli has only recently become technically feasible. Using a genetic approach, we have identified new Arabidopsis weak-ethylene insensitive mutants that also display defects in translation, which suggested the existence of a previously unknown molecular module involved in ethylene-mediated translation regulation of components of this signaling pathway. To explore this link in detail, we implemented for Arabidopsis the ribosome-footprinting technology, which enables the study of translation at a whole-genome level at single codon resolution[1]. Using ribosome-footprinting we examined the effects of short exposure to ethylene on the Arabidopsis translatome looking for ethylene-triggered changes in translation rates that could not be explained by changes in transcript levels. The results of this research, in combination with the characterization of a subset of the aforementioned weak-ethylene insensitive mutants that are defective in the UPF genes (core-components of the nonsense-mediated mRNA decay machinery), uncovered a translation-based branch of the ethylene signaling pathway[2]. In the presence of ethylene, translation of a negative regulator of ethylene signaling EBF2 is repressed, despite induced transcription of this gene. These translational effects of ethylene require the long 3´UTR of EBF2 (3´EBF2), which is recognized by the C-terminal end of the key ethylene-signaling protein EIN2 (EIN2C) in the cytoplasm once EIN2C is released from the ER-membrane by proteolytic cleavage. EIN2C binds the 3´EBF2, recruits the UPF proteins and moves to P-bodies, where the translation of EBF2 in inhibited despite its mRNA accumulation. Once the ethylene signal is withdrawn, the translation of the stored EBF2 mRNAs is resumed, thus rapidly dampening the ethylene response. These findings represent a mechanistic paradigm of gene-specific regulation of translation in response to a key growth regulator. Translation regulatory elements can be located in both 3′ and 5′ UTRs. We are now focusing on the ead1 and ead2 mutants, another set of ethylene-signaling mutants defective in translational regulation. Ribosome-footprinting on the ead1 mutant revealed an accumulation of translating ribosomes in the 5´UTRs of uORF-containing genes and reduction in the levels of ribosomes in the main ORF. The mutant is also impaired in the translation of GFP when this reporter is fused to WT 5´UTR of potential EAD1 targets but not when GFP is fused to the uORF-less versions of the same 5´UTRs. Our hypothesis is that EAD1/2 work as a complex that is required for the efficient translation of mRNAs that have common structural (complex 5´UTR with uORFs) and functional (regulation of key cellular processes) features. We are working towards the identification of the conditions where the EAD1 regulation of translation is required. [1] Ingolia, N. et al. (2009) Genome-Wide Analysis in Vivo of Translation with Nucleotide Resolution Using Ribosome Profiling. Science, 324; 218-222 [2] Merchante, C. et al. (2015) Gene-Specific Translation Regulation Mediated by the Hormone-Signaling Molecule EIN2. Cell, 163(3): 684-697

Comparative Genomics of Field Isolates of Mycobacterium bovis and M. caprae Provides Evidence for Possible Correlates with Bacterial Viability and Virulence

Mycobacteria of the Mycobacterium tuberculosis complex (MTBC) greatly affect humans and animals worldwide. The life cycle of mycobacteria is complex and the mechanisms resulting in pathogen infection and survival in host cells are not fully understood. Recently, comparative genomics analyses have provided new insights into the evolution and adaptation of the MTBC to survive inside the host. However, most of this information has been obtained using M. tuberculosis but not other members of the MTBC such as M. bovis and M. caprae. In this study, the genome of three M. bovis (MB1, MB3, MB4) and one M. caprae (MB2) field isolates with different lesion score, prevalence and host distribution phenotypes were sequenced. Genome sequence information was used for whole-genome and protein-targeted comparative genomics analysis with the aim of finding correlates with phenotypic variation with potential implications for tuberculosis (TB) disease risk assessment and control. At the whole-genome level the results of the first comparative genomics study of field isolates of M. bovis including M. caprae showed that as previously reported for M. tuberculosis, sequential chromosomal nucleotide substitutions were the main driver of the M. bovis genome evolution. The phylogenetic analysis provided a strong support for the M. bovis/M. caprae clade, but supported M. caprae as a separate species. The comparison of the MB1 and MB4 isolates revealed differences in genome sequence, including gene families that are important for bacterial infection and transmission, thus highlighting differences with functional implications between isolates otherwise classified with the same spoligotype. Strategic protein-targeted analysis using the ESX or type VII secretion system, proteins linking stress response with lipid metabolism, host T cell epitopes of mycobacteria, antigens and peptidoglycan assembly protein identified new genetic markers and candidate vaccine antigens that warrant further study to develop tools to evaluate risks for TB disease caused by M. bovis/M.caprae and for TB control in humans and animals.

Global transcriptome analysis of Lactococcus garvieae strains in response to temperature

Lactococcus garvieae is an important fish and an opportunistic human pathogen. The genomic sequences of several L. garvieae strains have been recently published, opening the possibility of global studies on the biology of this pathogen. In this study, a whole genome DNA microarray of two strains of L. garvieae was designed and validated. This DNA microarray was used to investigate the effects of growth temperature (18°C and 37°C) on the transcriptome of two clinical strains of L. garvieae that were isolated from fish (Lg8831) and from a human case of septicemia (Lg21881). The transcriptome profiles evidenced a strain-specific response to temperature, which was more evident at 18°C. Among the most significant findings, Lg8831 was found to up-regulate at 18°C several genes encoding different cold-shock and cold-induced proteins involved in an efficient adaptive response of this strain to low-temperature conditions. Another relevant result was the description, for the first time, of respiratory metabolism in L. garvieae, whose gene expression regulation was temperature-dependent in Lg21881. This study provides new insights about how environmental factors such as temperature can affect L. garvieae gene expression. These data could improve our understanding of the regulatory networks and adaptive biology of this important pathogen.

Staphylococcus aureus CC398: host adaptation and emergence of methicillin resistance in livestock

UNLABELLED Since its discovery in the early 2000s, methicillin-resistant Staphylococcus aureus (MRSA) clonal complex 398 (CC398) has become a rapidly emerging cause of human infections, most often associated with livestock exposure. We applied whole-genome sequence typing to characterize a diverse collection of CC398 isolates (n = 89), including MRSA and methicillin-susceptible S. aureus (MSSA) from animals and humans spanning 19 countries and four continents. We identified 4,238 single nucleotide polymorphisms (SNPs) among the 89 core genomes. Minimal homoplasy (consistency index = 0.9591) was detected among parsimony-informative SNPs, allowing for the generation of a highly accurate phylogenetic reconstruction of the CC398 clonal lineage. Phylogenetic analyses revealed that MSSA from humans formed the most ancestral clades. The most derived lineages were composed predominantly of livestock-associated MRSA possessing three different staphylococcal cassette chromosome mec element (SCCmec) types (IV, V, and VII-like) including nine subtypes. The human-associated isolates from the basal clades carried phages encoding human innate immune modulators that were largely missing among the livestock-associated isolates. Our results strongly suggest that livestock-associated MRSA CC398 originated in humans as MSSA. The lineage appears to have undergone a rapid radiation in conjunction with the jump from humans to livestock, where it subsequently acquired tetracycline and methicillin resistance. Further analyses are required to estimate the number of independent genetic events leading to the methicillin-resistant sublineages, but the diversity of SCCmec subtypes is suggestive of strong and diverse antimicrobial selection associated with food animal production. IMPORTANCE Modern food animal production is characterized by densely concentrated animals and routine antibiotic use, which may facilitate the emergence of novel antibiotic-resistant zoonotic pathogens. Our findings strongly support the idea that livestock-associated MRSA CC398 originated as MSSA in humans. The jump of CC398 from humans to livestock was accompanied by the loss of phage-carried human virulence genes, which likely attenuated its zoonotic potential, but it was also accompanied by the acquisition of tetracycline and methicillin resistance. Our findings exemplify a bidirectional zoonotic exchange and underscore the potential public health risks of widespread antibiotic use in food animal production.

FROM GENES TO GENOMES: LOCAL ADAPTATION AND ADAPTIVE POTENTIAL IN TWO ARCTIC SEABIRDS

By investigating the mechanisms underlying the evolution and the maintenance of local adaptations we can help predict how species will adapt to future environmental change. In this thesis I investigate local adaptation and adaptive potential in thick-billed and common murres (Uria lomvia and U. aalge), two arctic seabirds of international conservation concern. Thanks to the recent development of new genomic methods, I address three major themes that are relevant for both the development of evolutionary theory and conservation: 1) the role of gene flow in the origin and maintenance of adaptation; 2) levels and distribution of standing genetic variation, and their contribution to adaptive potential; and 3) the genomic mechanisms maintaining an adaptive dimorphism within a single interbreeding population. First, I review the literature on genomics of local adaptation with gene flow and find that adaptation can be maintained despite gene flow, that gene flow itself can promote adaptation, and that genetic architecture is important in the origin and maintenance of local adaptations. Second, I genotype genome-wide markers and toll-like receptor genes (TLRs) to investigate local adaptation and adaptive potential in thick-billed murres. Thick-billed murres do not show signatures of local adaptation to their breeding grounds, but outlier loci group birds according to their non-breeding distributions, suggesting that selection and/or demographic connectivity in the winter may explain patterns of differentiation in this species. Genetic variation at TLRs does not decrease with increasing latitude as predicted, but tests of selection and measures of genetic diversity suggest differences in local selective regimes at most genes. Thick-billed murres show high levels of standing genetic variation and their adaptive potential will mostly depend on rate and magnitude of environmental change. Finally, I improve and annotate the assembly of the highly heterozygous genome of the thick-billed murre. Using this assembly as a reference, I perform whole genome analyses to investigate the genomic basis of an adaptive dimorphism in Atlantic common murres. I show for the first time that a 60 kb complex copy number variant in a non-coding region maintains differences in plumage and cold adaptation despite high gene flow.

From field to fermentation: characterisation and application of non-dairy cultures in dairy foods

This thesis investigates the phenotypic and genotypic diversity of non-dairy L. lactis strains and their application to dairy fermentations. A bank of non-dairy lactococci were isolated from grass, vegetables and the bovine rumen. Subsequent analysis of these L. lactis strains revealed seven strains to possess cremoris genotypes which did not correlate with their observed phenotypes. Multi-locus sequence typing (MLST) and average nucleotide identity (ANI) highlighted the genetic diversity of lactis and cremoris subspecies. The application of these non-dairy lactococci to cheese production was also assessed. In milk, non-dairy strains formed diverse volatile profiles and selected strains were used as adjuncts in a mini Gouda-type cheese system. Sensory analysis showed non-dairy strains to be strongly associated with the development of off-flavours and bitterness. However, microfluidisation appeared to reduce bitterness. A novel bacteriophage, ɸL47, was isolated using the grass isolate L. lactis ssp. cremoris DPC6860 as a host. The phage, a member of the Siphoviridae, possessed a long tail fiber, previously unseen in dairy lactococcal phages. Genome sequencing revealed ɸL47 to be the largest sequenced lactococcal phage to date and owing to the high % similarity with ɸ949, a second member of the 949 group. Finally, to identify and characterise specific genes which may be important in niche adaptation and for applications to dairy fermentations, comparative genome sequence analysis was performed on L. lactis from corn (DPC6853), the bovine rumen (DPC6853) and grass (DPC6860). This study highlights the contribution of niche specialisation to the intra-species diversity of L. lactis and the adaptation of this organism to different environments. In summary this thesis describes the genetic diversity of L. lactis strains from outside the dairy environment and their potential application in dairy fermentations.

A transcription factor map as revealed by a genome-wide gene expression analysis of whole-blood mRNA transcriptome in multiple sclerosis

Background Several lines of evidence suggests that transcription factors are involved in the pathogenesis of Multiple Sclerosis (MS) but a complete mapping the whole network has been elusive. One of the reasons is that there are several clinical subtypes of MS and transcription factors which may be involved in one subtype may not be in others. We investigated the possibility that this network could be mapped using microarray technologies and modern bioinformatics methods on a dataset from whole blood in 99 untreated MS patients (36 Relapse Remitting MS, 43 Primary Progressive MS, and 20 Secondary Progressive MS) and 45 age-matched healthy controls, Methodology/Principal Findings We have used two different analytical methodologies: a differential expression analysis and a differential co-expression analysis, which have converged on a significant number of regulatory motifs that seem to be statistically overrepresented in genes which are either differentially expressed (or differentially co-expressed) in cases and controls (e.g. V$KROX_Q6, p-value < 3.31E-6; V$CREBP1_Q2, p-value < 9.93E-6, V$YY1_02, p-value < 1.65E-5). Conclusions/significance: Our analysis uncovered a network of transcription factors that potentially dysregulate several genes in MS or one or more of its disease subtypes. Analysing the published literature we have found that these transcription factors are involved in the early T-lymphocyte specification and commitment as well as in oligodendrocytes dedifferentiation and development. The most significant transcription factors motifs were for the Early Growth response EGR/KROX family, ATF2, YY1 (Yin and Yang 1), E2F-1/DP-1 and E2F-4/DP-2 heterodimers, SOX5, and CREB and ATF families.

Optimization of sample preparation for next-generation sequencing with Illumina Genome Analyzer II

The growing interest for sequencing with higher throughput in the last decade has led to the development of new sequencing applications. This thesis concentrates on optimizing DNA library preparation for Illumina Genome Analyzer II sequencer. The library preparation steps that were optimized include fragmentation, PCR purification and quantification. DNA fragmentation was performed with focused sonication in different concentrations and durations. Two column based PCR purification method, gel matrix method and magnetic bead based method were compared. Quantitative PCR and gel electrophoresis in a chip were compared for DNA quantification. The magnetic bead purification was found to be the most efficient and flexible purification method. The fragmentation protocol was changed to produce longer fragments to be compatible with longer sequencing reads. Quantitative PCR correlates better with the cluster number and should thus be considered to be the default quantification method for sequencing. As a result of this study more data have been acquired from sequencing with lower costs and troubleshooting has become easier as qualification steps have been added to the protocol. New sequencing instruments and applications will create a demand for further optimizations in future.

Elucidating the cancer-specific genetic alteration spectrum of glioblastoma derived cell lines from whole exome and RNA sequencing

Cell lines derived from tumor tissues have been used as a valuable system to study gene regulation and cancer development. Comprehensive characterization of the genetic background of cell lines could provide clues on novel genes responsible for carcinogenesis and help in choosing cell lines for particular studies. Here, we have carried out whole exome and RNA sequencing of commonly used glioblastoma (GBM) cell lines (U87, T98G, LN229, U343, U373 and LN18) to unearth single nucleotide variations (SNVs), indels, differential gene expression, gene fusions and RNA editing events. We obtained an average of 41,071 SNVs out of which 1,594 (3.88%) were potentially cancer-specific. The cell lines showed frequent SNVs and indels in some of the genes that are known to be altered in GBM-EGFR, TP53, PTEN, SPTA1 and NF1. Chromatin modifying genes-ATRX, MLL3, MLL4, SETD2 and SRCAP also showed alterations. While no cell line carried IDH1 mutations, five cell lines showed hTERT promoter activating mutations with a concomitant increase in hTERT transcript levels. Five significant gene fusions were found of which NUP93-CYB5B was validated. An average of 18,949 RNA editing events was also obtained. Thus we have generated a comprehensive catalogue of genetic alterations for six GBM cell lines.

Elucidating the cancer-specific genetic alteration spectrum of glioblastoma derived cell lines from whole exome and RNA sequencing

Cell lines derived from tumor tissues have been used as a valuable system to study gene regulation and cancer development. Comprehensive characterization of the genetic background of cell lines could provide clues on novel genes responsible for carcinogenesis and help in choosing cell lines for particular studies. Here, we have carried out whole exome and RNA sequencing of commonly used glioblastoma (GBM) cell lines (U87, T98G, LN229, U343, U373 and LN18) to unearth single nucleotide variations (SNVs), indels, differential gene expression, gene fusions and RNA editing events. We obtained an average of 41,071 SNVs out of which 1,594 (3.88%) were potentially cancer-specific. The cell lines showed frequent SNVs and indels in some of the genes that are known to be altered in GBM-EGFR, TP53, PTEN, SPTA1 and NF1. Chromatin modifying genes-ATRX, MLL3, MLL4, SETD2 and SRCAP also showed alterations. While no cell line carried IDH1 mutations, five cell lines showed hTERT promoter activating mutations with a concomitant increase in hTERT transcript levels. Five significant gene fusions were found of which NUP93-CYB5B was validated. An average of 18,949 RNA editing events was also obtained. Thus we have generated a comprehensive catalogue of genetic alterations for six GBM cell lines.

Hybrid error correction and de novo assembly of single-molecule sequencing reads.

Single-molecule sequencing instruments can generate multikilobase sequences with the potential to greatly improve genome and transcriptome assembly. However, the error rates of single-molecule reads are high, which has limited their use thus far to resequencing bacteria. To address this limitation, we introduce a correction algorithm and assembly strategy that uses short, high-fidelity sequences to correct the error in single-molecule sequences. We demonstrate the utility of this approach on reads generated by a PacBio RS instrument from phage, prokaryotic and eukaryotic whole genomes, including the previously unsequenced genome of the parrot Melopsittacus undulatus, as well as for RNA-Seq reads of the corn (Zea mays) transcriptome. Our long-read correction achieves >99.9% base-call accuracy, leading to substantially better assemblies than current sequencing strategies: in the best example, the median contig size was quintupled relative to high-coverage, second-generation assemblies. Greater gains are predicted if read lengths continue to increase, including the prospect of single-contig bacterial chromosome assembly.

Next-generation sequencing of the mitochondrial genome and association with IgA nephropathy in a renal transplant population

Kidneys are highly aerobic organs that are critically dependent on the normal functioning of mitochondria. Genetic variations disrupting mitochondrial function are associated with multifactorial disorders including kidney disease. This study sequenced the entire mitochondrial genome in a renal transplant cohort of 64 individuals, using next-generation sequencing, to evaluate the association of genetic variants with IgA nephropathy and end-stage renal disease (ESRD, n = 100).