Comparative transcriptome profiling of human foreskin fibroblasts infected with the Sylvio and Y strains of Trypanosoma cruzi

Trypanosoma cruzi, the causative agent of Chagas Disease, is phylogenetically distributed into nearly identical genetic strains which show divergent clinical presentations including differences in rates of cardiomyopathy in humans, different vector species and transmission cycles, and differential congenital transmission in a mouse model. The population structure of these strains divides into two groups, which are geographically and clinically distinct. The aim of this study was to compare the transcriptome of two strains of T. cruzi, Sylvio vs. Y to identify differences in expression that could account for clinical and biochemical differences. We collected and sequenced RNA from T. cruzi-infected and control Human Foreskin Fibroblasts at three timepoints. Differential expression analysis identified gene expression profiles at different timepoints in Sylvio infections, and between Sylvio and Y infections in both parasite and host. The Sylvio strain parasite and the host response to Sylvio infection largely mirrored the host-pathogen interaction seen in our previous Y strain work. IL-8 was more highly expressed in Sylvio-infected HFFs than in Y-infected HFFs.

Transcriptome, carbohydrate, and phytohormone analysis of Petunia hybrida reveals a complex disturbance of plant functional integrity under mild chilling stress

Cultivation of chilling-tolerant ornamental crops at lower temperature could reduce the energy demands of heated greenhouses. To provide a better understanding of how sub-optimal temperatures (12 degrees C vs. 16 degrees C) affect growth of the sensitive Petunia hybrida cultivar 'SweetSunshine Williams', the transcriptome, carbohydrate metabolism, and phytohormone homeostasis were monitored in aerial plant parts over 4 weeks by use of a microarray, enzymatic assays and GC-MS/MS. The data revealed three consecutive phases of chilling response. The first days were marked by a strong accumulation of sugars, particularly in source leaves, preferential up-regulation of genes in the same tissue and down-regulation of several genes in the shoot apex, especially those involved in the abiotic stress response. The midterm phase featured a partial normalization of carbohydrate levels and gene expression. After 3 weeks of chilling exposure, a new stabilized balance was established. Reduced hexose levels in the shoot apex, reduced ratios of sugar levels between the apex and source leaves and a higher apical sucrose/hexose ratio, associated with decreased activity and expression of cell wall invertase, indicate that prolonged chilling induced sugar accumulation in source leaves at the expense of reduced sugar transport to and reduced sucrose utilization in the shoot. This was associated with reduced levels of indole-3-acetic acid and abscisic acid in the apex and high numbers of differentially, particularly up-regulated genes, especially in the source leaves, including those regulating histones, ethylene action, transcription factors, and a jasmonate-ZIM-domain protein. Transcripts of one Jumonji C domain containing protein and one expansin accumulated in source leaves throughout the chilling period. The results reveal a dynamic and complex disturbance of plant function in response to mild chilling, opening new perspectives for the comparative analysis of differently tolerant cultivars.

A transcriptome approach towards understanding the development of ripening capacity in 'Bartlett' pears (Pyrus communis L.)

Background: The capacity of European pear fruit (Pyrus communis L.) to ripen after harvest develops during the final stages of growth on the tree. The objective of this study was to characterize changes in 'Bartlett' pear fruit physico-chemical properties and transcription profiles during fruit maturation leading to attainment of ripening capacity. Results: The softening response of pear fruit held for 14days at 20°C after harvest depended on their maturity. We identified four maturity stages: S1-failed to soften and S2- displayed partial softening (with or without ET-ethylene treatment); S3 - able to soften following ET; and S4 - able to soften without ET. Illumina sequencing and Trinity assembly generated 68,010 unigenes (mean length of 911bp), of which 32.8% were annotated to the RefSeq plant database. Higher numbers of differentially expressed transcripts were recorded in the S3-S4 and S1-S2 transitions (2805 and 2505 unigenes, respectively) than in the S2-S3 transition (2037 unigenes). High expression of genes putatively encoding pectin degradation enzymes in the S1-S2 transition suggests pectic oligomers may be involved as early signals triggering the transition to responsiveness to ethylene in pear fruit. Moreover, the co-expression of these genes with Exps (Expansins) suggests their collaboration in modifying cell wall polysaccharide networks that are required for fruit growth. K-means cluster analysis revealed that auxin signaling associated transcripts were enriched in cluster K6 that showed the highest gene expression at S3. AP2/EREBP (APETALA 2/ethylene response element binding protein) and bHLH (basic helix-loop-helix) transcripts were enriched in all three transition S1-S2, S2-S3, and S3-S4. Several members of Aux/IAA (Auxin/indole-3-acetic acid), ARF (Auxin response factors), and WRKY appeared to play an important role in orchestrating the S2-S3 transition. Conclusions: We identified maturity stages associated with the development of ripening capacity in 'Bartlett' pear, and described the transcription profile of fruit at these stages. Our findings suggest that auxin is essential in regulating the transition of pear fruit from being ethylene-unresponsive (S2) to ethylene-responsive (S3), resulting in fruit softening. The transcriptome will be helpful for future studies about specific developmental pathways regulating the transition to ripening. © 2015 Nham et al.

Development of single nucleotide polymorphism (SNP) markers from the mango (Mangifera indica) transcriptome for mapping and estimation of genetic diversity

The development of molecular markers for genomic studies in Mangifera indica (mango) will allow marker-assisted selection and identification of genetically diverse germplasm, greatly aiding mango breeding programs. We report here our identification of thousands of unambiguous molecular markers that can be easily assayed across genotypes of the species. With origin centered in Southeast Asia, mangos are grown throughout the tropics and subtropics as a nutritious fruit that exhibits remarkable intraspecific phenotypic diversity. With the goal of building a high density genetic map, we have undertaken discovery of sequence variation in expressed genes across a broad range of mango cultivars. A transcriptome sequence reference was built de novo from extensive sequencing and assembly of RNA from cultivar 'Tommy Atkins'. Single nucleotide polymorphisms (SNPs) in protein coding transcripts were determined from alignment of RNA reads from 24 mango cultivars of diverse origins: 'Amin Abrahimpur' (India), 'Aroemanis' (Indonesia), 'Burma' (Burma), 'CAC' (Hawaii), 'Duncan' (Florida), 'Edward' (Florida), 'Everbearing' (Florida), 'Gary' (Florida), 'Hodson' (Florida), 'Itamaraca' (Brazil), 'Jakarata' (Florida), 'Long' (Jamaica), 'M. Casturi Purple' (Borneo), 'Malindi' (Kenya), 'Mulgoba' (India), 'Neelum' (India), 'Peach' (unknown), 'Prieto' (Cuba), 'Sandersha' (India), 'Tete Nene' (Puerto Rico), 'Thai Everbearing' (Thailand), 'Toledo' (Cuba), 'Tommy Atkins' (Florida) and 'Turpentine' (West Indies). SNPs in a selected subset of protein coding transcripts are currently being converted into Fluidigm assays for genotyping of mapping populations and germplasm collections. Using an alternate approach, SNPs (144) discovered by sequencing of candidate genes in 'Kensington Pride' have already been converted and used for genotyping.

Transcriptome Analysis of Capsicum Chlorosis Virus-Induced Hypersensitive Resistance Response in Bell Capsicum

Background Capsicum chlorosis virus (CaCV) is an emerging pathogen of capsicum, tomato and peanut crops in Australia and South-East Asia. Commercial capsicum cultivars with CaCV resistance are not yet available, but CaCV resistance identified in Capsicum chinense is being introgressed into commercial Bell capsicum. However, our knowledge of the molecular mechanisms leading to the resistance response to CaCV infection is limited. Therefore, transcriptome and expression profiling data provide an important resource to better understand CaCV resistance mechanisms. Methodology/Principal Findings We assembled capsicum transcriptomes and analysed gene expression using Illumina HiSeq platform combined with a tag-based digital gene expression system. Total RNA extracted from CaCV/mock inoculated CaCV resistant (R) and susceptible (S) capsicum at the time point when R line showed a strong hypersensitive response to CaCV infection was used in transcriptome assembly. Gene expression profiles of R and S capsicum in CaCV- and buffer-inoculated conditions were compared. None of the genes were differentially expressed (DE) between R and S cultivars when mock-inoculated, while 2484 genes were DE when inoculated with CaCV. Functional classification revealed that the most highly up-regulated DE genes in R capsicum included pathogenesis-related genes, cell death-associated genes, genes associated with hormone-mediated signalling pathways and genes encoding enzymes involved in synthesis of defense-related secondary metabolites. We selected 15 genes to confirm DE expression levels by real-time quantitative PCR. Conclusion/Significance DE transcript profiling data provided comprehensive gene expression information to gain an understanding of the underlying CaCV resistance mechanisms. Further, we identified candidate CaCV resistance genes in the CaCV-resistant C. annuum x C. chinense breeding line. This knowledge will be useful in future for fine mapping of the CaCV resistance locus and potential genetic engineering of resistance into CaCV-susceptible crops.

Microarrays as a functional approach to the transcriptome

Knowing a cell’s transcriptome is a fundamental requisite in order to analyze its response to the environment. Microarrays have supposed a revolution on this field as they are able to yield an overview of gene expression at any environmental condition on a genome-wide scale. This technique consists in the hybridisation of a nucleic acid sample, previously marked, with a probe (which might be made up of cDNA, oligonucleotides or PCR products) anchored to a solid surface (made of glass, plastic, silicon...) giving as a result a dot grid which reveals, after image analysis, which genes are being expressed. Nevertheless, this only can be achieved if information on the species genome has been generated. Different kinds of expression microarrays exist attending to the probe’s nature and the method used in its synthesis. In this poster two of these will be treated: Spotted Microarrays, for which the probe is synthesised prior to its fixation to the array and allow the analysis of two targets simultaneously. They can be easily customized, but lack high reproducibility and sensitivity. Oligonucleotide Microarrays, which are characterized by the direct printing of the probe on the array. In this case the probes consist on, invariably, oligonucleotides that are complementary to a small fraction of the gene it is representing at the microarray. Their application is somewhat restricted. This fact, however, makes them more reproducible. Currently, the approach towards the transcriptome studies from the Next Generation Sequencing technologies offers a large volume of information in a short amount of time needing less previous information on the target organism than that needed by microarrays, but their expensive price limits their use. The versatility of the latter, together with their reduced costs in comparison to other techniques, makes them an interesting resource in applications that may need less complexity.

The transcriptome of metamorphosing flatfish

Background Flatfish metamorphosis denotes the extraordinary transformation of a symmetric pelagic larva into an asymmetric benthic juvenile. Metamorphosis in vertebrates is driven by thyroid hormones (THs), but how they orchestrate the cellular, morphological and functional modifications associated with maturation to juvenile/adult states in flatfish is an enigma. Since THs act via thyroid receptors that are ligand activated transcription factors, we hypothesized that the maturation of tissues during metamorphosis should be preceded by significant modifications in the transcriptome. Targeting the unique metamorphosis of flatfish and taking advantage of the large size of Atlantic halibut (Hippoglossus hippoglossus) larvae, we determined the molecular basis of TH action using RNA sequencing. Results De novo assembly of sequences for larval head, skin and gastrointestinal tract (GI-tract) yielded 90,676, 65,530 and 38,426 contigs, respectively. More than 57 % of the assembled sequences were successfully annotated using a multi-step Blast approach. A unique set of biological processes and candidate genes were identified specifically associated with changes in morphology and function of the head, skin and GI-tract. Transcriptome dynamics during metamorphosis were mapped with SOLiD sequencing of whole larvae and revealed greater than 8,000 differentially expressed (DE) genes significantly (p < 0.05) up- or down-regulated in comparison with the juvenile stage. Candidate transcripts quantified by SOLiD and qPCR analysis were significantly (r = 0.843; p < 0.05) correlated. The majority (98 %) of DE genes during metamorphosis were not TH-responsive. TH-responsive transcripts clustered into 6 groups based on their expression pattern during metamorphosis and the majority of the 145 DE TH-responsive genes were down-regulated. Conclusions A transcriptome resource has been generated for metamorphosing Atlantic halibut and over 8,000 DE transcripts per stage were identified. Unique sets of biological processes and candidate genes were associated with changes in the head, skin and GI-tract during metamorphosis. A small proportion of DE transcripts were TH-responsive, suggesting that they trigger gene networks, signalling cascades and transcription factors, leading to the overt changes in tissue occurring during metamorphosis.

Mineral stress response of Vitis cell wall transcriptome: identification, characterization and expression of genes from key families involved in wall biosynthesis and modification

Doutoramento em Engenharia Agronómica - Instituto Superior de Agronomia - UL

Next-generation transcriptome assembly of an Amazon palm (Euterpe precatoria).

2016

Transcriptome analysis of Euterpe edulis and identification of microsatellite markers.

2016

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.

A Transcriptome Atlas of Physcomitrella patens Provides Insights into the Evolution and Development of Land Plants

Post-print version of the article.

Régulation du transcriptome codant et non-codant chez Schizosaccharomyces pombe: facteurs et mécanismes impliqués dans la maturation 3’ des ARNs et la terminaison de la transcription

La synthèse d’un ARNm eucaryotique dépend d’une suite d’étapes qui inclut notamment l’ajout d’une queue poly(A) à son extrémité 3’. Au noyau, la queue poly(A) des ARNms est liée par PABPN1 (poly(A)-binding protein nuclear 1). PABPN1 fut notamment caractérisée, d’après des études in vitro, pour stimuler la réaction de polyadénylation en plus de contrôler la taille ultime des queues poly(A). Cela dit, la ou les fonction(s) biologique(s) de PABPN1 est/sont cependant largement méconnue(s). Chez Schizosaccharomyces pombe (S. pombe), Pab2 est l’orthologue présumé de PABPN1. Or, mes travaux indiquent que Pab2 est fonctionnellement différente de PABPN1 à l’égard de son rôle sur le processus général de polyadénylation. Ainsi, in vivo, l’absence de Pab2 entraîne l’expression et l’accumulation d’un groupe limité d’ARNs hyperadénylés parmi lesquels se trouvent de nombreux petits ARNs nucléolaires non-codants (snoRNAs) lesquels constituent normalement un groupe abondant d’ARN poly(A)-. Mes résultats supportent ainsi un mécanisme par lequel des snoRNAs immatures poly(A)+, sont convertis en une forme mature poly(A)- par le biais de Pab2 et de l’activité 3’-->5’ exoribonucléase de l’exosome à ARN. Ces observations sont inusitées dans la mesure où elles associent une fonction pour une PABP dans la maturation d'ARNs non-codants, contrairement à la notion que les PABPs travaillent exclusivement au niveau des ARNms, en plus de procurer une nouvelle perspective face au mécanisme de recrutement de l'exosome à ARN à des substrats poly(A)+. La formation de l’extrémité 3’ d’un ARN est un processus étroitement lié à la terminaison de sa transcription. Pour les gènes codants, la terminaison transcriptionnelle est initiée par le clivage endonucléolytique du pré-ARNm. Ce clivage génère une extrémité d’ARN 5’ libre laquelle sera ciblée par une exoribonucléase 5'-->3’ afin de mener à bien l’éviction de l’ARNPII de la matrice d’ADN (terminaison transcriptionnelle de type torpedo). Au contraire, chez Saccharomyces cerevisiae (S. cerevisiae), la majorité des gènes non-codants, incluant les snoRNAs, dépendent plutôt du complexe NNS (Nrd1/Nab3/Sen1) pour la terminaison de leur transcription. Cela dit, il est incertain si le complexe NNS est conservé chez d’autres espèces. À cet égard, mes travaux indiquent que S. pombe est dépourvu d’un mécanisme de terminaison de la transcription de type NNS. Seb1, l’orthologue présumé de Nrd1 chez S. pombe, s’associe plutôt à la machinerie de clivage et de polyadénylation et influence la sélection de site de polyadénylation à l’échelle du génome. Mes résultats supportent ainsi l’utilisation de la machinerie de maturation 3’ des ARNms comme principal vecteur de terminaison transcriptionnelle chez S. pombe et identifient Seb1 comme un facteur clé de ce processus. L’évènement transcriptionnel étant hautement complexe, des erreurs peuvent arriver de manière stochastique menant à l’accumulation d’ARNs aberrants potentiellement néfastes pour la cellule. Or, mes travaux ont mis en lumière un mécanisme de surveillance co-transcriptionnel des ARNs impliquant l’exosome à ARN et lié à la terminaison de la transcription. Pour ce faire, l’exosome à ARN promeut la terminaison transcriptionnelle via la dégradation d’une extrémité 3’ libre d’ARN devenue émergente suite au recul de l’ARNPII le long de la matrice d’ADN (phénomène de backtracking). Mes résultats supportent ainsi une terminaison de la transcription de type torpedo inversé (3'-->5’) réévaluant par la même occasion le concept voulant que la terminaison de la transcription s’effectue uniquement selon une orientation 5’-->3’. Somme toute, mes travaux de doctorat auront permis d’identifier et de caractériser plus en détail les facteurs et mécanismes impliqués dans la maturation 3’ et la terminaison de la transcription des gènes codants et non-codants chez l’organisme modèle S. pombe.

Transcriptome analysis of wheat under drought.

2013

Functional relationship of colorectal cancer-associated glycans with malignant phenotype and transcriptome changes

Glycosyltransferases ST6GAL1 and B4GALNT2 (and their cognate antigens Sia6LacNAc and Sda, respectively) are associated with colorectal cancer (CRC) but it is not fully clear their biological and clinical significance. We explored the clinical relevance of both glycosyltransferases by interrogating The Cancer Genome Atlas (TCGA) database while the phenotypic/transcriptomic effects of ST6GAL1/B4GALNT2 overexpression were studied in genetically modified CRC cell lines. Transcriptomic data from CRC patients in TCGA database suggested a moderate impact of ST6GAL1 on CRC progression, although it was not possible to define a clear role for this glycosyltransferase. Transcriptomic analysis of ST6GAL1-transduced cell lines revealed a much deeper effect of ST6GAL1 on gene expression in SW948 than in SW48. The overexpression of ST6GAL1 induced opposite effects on soft agar growth and wound healing in both cell lines. These results indicate that the impact of a cancer-associated glycosyltransferase change on phenotype/transcriptome can be extremely variable, depending on the molecular context of the tumor cell. On the contrary, transcriptomic analysis of B4GALNT2-modified cell lines together with TCGA database survey demonstrated a strong impact of B4GALNT2 on the transcriptional activity of CRC cells, in particular its association with a better prognosis. We suggest an anti-tumoral role of B4GALNT2 in CRC. We also investigated the glycan changes related to ST6GAL1/B4GALNT2 expression in a small cohort of tissues/plasma as well as the N-glycomic profile of CRC, normal and polyp tissues. We found an increase of ST6GAL1 activity in CRC and inflammatory bowel disease plasma samples comparing with plasma from healthy donors. A different Sda protein carrier pattern was observed between healthy donors and CRC plasma samples. β-arrestin 1 is a possible candidate as Sda carrier protein in plasma samples although future validation studies are needed. The alterations found in the N-glycan pattern highlight the importance of N-glycome as a molecular signature in cancer.