The multiple sclerosis whole blood mRNA transcriptome and genetic associations indicate dysregulation of specific T cell pathways in pathogenesis

Multiple sclerosis (MS) is an autoimmune disease with a genetic component, caused at least in part by aberrant lymphocyte activity. The whole blood mRNA transcriptome was measured for 99 untreated MS patients: 43 primary progressive MS, 20 secondary progressive MS, 36 relapsing remitting MS and 45 age-matched healthy controls. The ANZgene Multiple Sclerosis Genetics Consortium genotyped more than 300 000 SNPs for 115 of these samples. Transcription from genes on translational regulation, oxidative phosphorylation, immune synapse and antigen presentation pathways was markedly increased in all forms of MS. Expression of genes tagging T cells was also upregulated (P < 10-12) in MS. A T cell gene signature predicts disease state with a concordance index of 0.79 with age and gender as co-variables, but the signature is not associated with clinical course or disability. The ANZgene genome wide association screen identified two novel regions with genome wide significance: one encoding the T cell co-stimulatory molecule, CD40; the other a region on chromosome 12q13-14. The CD40 haplotype associated with increased MS susceptibility has decreased gene expression in MS (P < 0.0007). The second MS susceptibility region includes 17 genes on 12q13-14 in tight linkage disequilibrium. Of these, only 13 are expressed in leukocytes, and of these the expression of one, FAM119B, is much lower in the susceptibility haplotype (P tdthomlt; 10-14). Overall, these data indicate dysregulation of T cells can be detected in the whole blood of untreated MS patients, and supports targeting of activated T cells in therapy for all forms of MS.

Regulation of strawberry growth and development

Strawberries (Fragaria sp.) are adapted to diverse environmental conditions from the tropics to about 70ºN, so different responses to environmental conditions can be found. Most genotypes of garden strawberry (F. x ananassa Duch.) and woodland strawberry (F. vesca L.) are short-day (SD) plants that are induced to flowering by photoperiods under a critical limit, but also various photoperiod x temperature interactions can be found. In addition, continuously flowering everbearing (EB) genotypes are found. In addition to flowering, axillary bud differentiation in strawberry is regulated by photoperiod. In SD conditions, axillary buds differentiate to rosette-like structures called "branch crowns", whereas in long-day conditions (LD) they form runners, branches with 2 long internodes followed by a daughter plant (leaf rosette). The number of crown branches determines the yield of the plant, since inflorescences are formed from the apical meristems of the crown. Although axillary bud differentiation is an important developmental process in strawberries, its environmental and hormonal regulation has not been characterized in detail. Moreover, the genetic mechanisms underlying axillary bud differentiation and regulation of flowering time in these species are almost completely unresolved. These topics have been studied in this thesis in order to enhance strawberry research, cultivation and breeding. The results showed that 8-12 SD cycles suppressed runner initiation from the axillary buds of the garden strawberry cv. Korona with the concomitant induction of crown branching, and 3 weeks of SD was sufficient for the induction of flowering in the main crown. Furthermore, a second SD treatment given a few weeks after the first SD period can be used to induce flowering in the primary branch crowns and to induce the formation of secondary branches. Thus, artificial SD treatments effectively stimulate crown branching, providing one means for the increase of cropping (yield) potential in strawberry. It was also shown by growth regulation applications, quantitave hormone analysis and gene expression analysis that gibberellin (GA) is one of the key signals involved in the photoperiod control of shoot differentiation. The results indicate that photoperiod controls GA activity specifically in axillary buds, thereby determining bud fate. It was further shown that chemical control of GA biosynthesis by prohexadione-calcium can be utilized to prevent excessive runner formation and induce crown branching in strawberry fields. Moreover, ProCa increased berry yield up to 50%, showing that it is an easier and more applicable alternative to artificial SD treatments for controlling strawberry crown development and yield. Finally, flowering gene pathways in Fragaria were explored by searching for homologs of 118 Arabidopsis thaliana flowering-time genes. In total, 66 gene homologs were identified, and they distributed to all known flowering pathways, suggesting the presence of these pathways also in strawberry. Expression analysis of selected genes revealed that the mRNA of putative floral identity gene APETALA1 accumulated in the shoot apex of the EB genotype after the induction of flowering, whereas it was absent in vegetative SD genotype, indicating the usefulness of this gene product as the marker of floral initiation. The present data enables the further exploration of strawberry flowering pathways with genetic transformation, gene mapping and transcriptomics methods.

Determination of the kallikrein-related peptidase 7 degradome and Transciptome in ovarian cancer

Ovarian cancer is the leading cause of cancer-related death in women, and the need for curative treatments is urgent. This study characterised an enzyme associated with the most lethal form of ovarian cancer, showing this enzyme to be a promising therapeutic target. Fifteen novel protein targets and key signalling pathways were determined to be regulated by this enzyme, kallikrein-related peptidase 7, in the ovarian tumour microenvironment, highlighting its involvement in cancer progression. Inhibition of this enzyme may be a useful therapeutic option to improve the life expectancy of women suffering from this cancer.

Evaluation of potential candidate genes involved in salinity tolerance in striped catfish (Pangasianodon hypophthalmus) using an RNA-Seq approach

Increasing salinity levels in freshwater and coastal environments caused by sea level rise linked to climate change is now recognized to be a major factor that can impact fish growth negatively, especially for freshwater teleost species. Striped catfish (Pangasianodon hypophthalmus) is an important freshwater teleost that is now widely farmed across the Mekong River Delta in Vietnam. Understanding the basis for tolerance and adaptation to raised environmental salinity conditions can assist the regional culture industry to mitigate predicted impacts of climate change across this region. Attempt of next generation sequencing using the ion proton platform results in more than 174 million raw reads from three tissue libraries (gill, kidney and intestine). Reads were filtered and de novo assembled using a variety of assemblers and then clustered together to generate a combined reference transcriptome. Downstream analysis resulted in a final reference transcriptome that contained 60,585 transcripts with an N50 of 683 bp. This resource was further annotated using a variety of bioinformatics databases, followed by differential gene expression analysis that resulted in 3062 transcripts that were differentially expressed in catfish samples raised under two experimental conditions (0 and 15 ppt). A number of transcripts with a potential role in salinity tolerance were then classified into six different functional gene categories based on their gene ontology assignments. These included; energy metabolism, ion transportation, detoxification, signal transduction, structural organization and detoxification. Finally, we combined the data on functional salinity tolerance genes into a hypothetical schematic model that attempted to describe potential relationships and interactions among target genes to explain the molecular pathways that control adaptive salinity responses in P. hypophthalmus. Our results indicate that P. hypophthalmus exhibit predictable plastic regulatory responses to elevated salinity by means of characteristic gene expression patterns, providing numerous candidate genes for future investigations.

A note on sampling chironomids for RNA-based studies of natural populations that retains critical morphological vouchers

The rapid uptake of transcriptomic approaches in freshwater ecology has seen a wealth of data produced concerning the ways in which organisms interact with their environment on a molecular level. Typically, such studies focus either at the community level and so don’t require species identifications, or on laboratory strains of known species identity or natural populations of large, easily identifiable taxa. For chironomids, impediments still exist for applying these technologies to natural populations because they are small-bodied and often require time-consuming secondary sorting of stream material and morphological voucher preparation to confirm species diagnosis. These procedures limit the ability to maintain RNA quantity and quality in such organisms because RNA degrades rapidly and gene expression can be altered rapidly in organisms; thereby limiting the inclusion of such taxa in transcriptomic studies. Here, we demonstrate that these limitations can be overcome and outline an optimised protocol for collecting, sorting and preserving chironomid larvae that enables retention of both morphological vouchers and RNA for subsequent transcriptomics purposes. By ensuring that sorting and voucher preparation are completed within <4 hours after collection and that samples are kept cold at all times, we successfully retained both RNA and morphological vouchers from all specimens. Although not prescriptive in specific methodology, we anticipate that this paper will assist in promoting transcriptomic investigations of the sublethal impact on chironomid gene expression of changes to aquatic environments.

Expression profiling of lymph nodes in tuberculosis patients reveal inflammatory milieu at site of infection

Extrapulmonary manifestations constitute 15 to 20% of tuberculosis cases, with lymph node tuberculosis (LNTB) as the most common form of infection. However, diagnosis and treatment advances are hindered by lack of understanding of LNTB biology. To identify host response, Mycobacterium tuberculosis infected lymph nodes from LNTB patients were studied by means of transcriptomics and quantitative proteomics analyses. The selected targets obtained by comparative analyses were validated by quantitative PCR and immunohistochemistry. This approach provided expression data for 8,728 transcripts and 102 proteins, differentially regulated in the infected human lymph node. Enhanced inflammation with upregulation of T-helper1-related genes, combined with marked dysregulation of matrix metalloproteinases, indicates tissue damage due to high immunoactivity at infected niche. This expression signature was accompanied by significant upregulation of an immunoregulatory gene, leukotriene A4 hydrolase, at both transcript and protein levels. Comparative transcriptional analyses revealed LNTB-specific perturbations. In contrast to pulmonary TB-associated increase in lipid metabolism, genes involved in fatty-acid metabolism were found to be downregulated in LNTB suggesting differential lipid metabolic signature. This study investigates the tissue molecular signature of LNTB patients for the first time and presents findings that indicate the possible mechanism of disease pathology through dysregulation of inflammatory and tissue-repair processes.

Brca1 is expressed in human microglia and is dysregulated in human and animal model of ALS

Background: There is growing evidence that microglia are key players in the pathological process of amyotrophic lateral sclerosis (ALS). It is suggested that microglia have a dual role in motoneurone degeneration through the release of both neuroprotective and neurotoxic factors. Results: To identify candidate genes that may be involved in ALS pathology we have analysed at early symptomatic age (P90), the molecular signature of microglia from the lumbar region of the spinal cord of hSOD1(G93A) mice, the most widely used animal model of ALS. We first identified unique hSOD1(G93A) microglia transcriptomic profile that, in addition to more classical processes such as chemotaxis and immune response, pointed toward the potential involvement of the tumour suppressor gene breast cancer susceptibility gene 1 (Brca1). Secondly, comparison with our previous data on hSOD1(G93A) motoneurone gene profile substantiated the putative contribution of Brca1 in ALS. Finally, we established that Brca1 protein is specifically expressed in human spinal microglia and is up-regulated in ALS patients. Conclusions: Overall, our data provide new insights into the pathogenic concept of a non-cell-autonomous disease and the involvement of microglia in ALS. Importantly, the identification of Brca1 as a novel microglial marker and as possible contributor in both human and animal model of ALS may represent a valid therapeutic target. Moreover, our data points toward novel research strategies such as investigating the role of oncogenic proteins in neurodegenerative diseases.

Nutrigenómica

Projeto de Pós-Graduação/Dissertação apresentado à Universidade Fernando Pessoa como parte dos requisitos para obtenção do grau de Mestre em Ciências Farmacêuticas

Functional genomics of motile commensal intestinal bacteria

Flagella confer upon bacteria the ability to move and are therefore organelles of significant bacteriological importance. The innate immune system has evolved to recognise flagellin, (the major protein component of the bacterial flagellar filament). Flagellate microbes can potentially stimulate the immune systems of mammals, and thus have significant immunomodulatory potential. The flagellum-biogenesis genotype and phenotype of Lactobacillus ruminis, an autochthonous intestinal commensal, was studied. The flagellum-biogenesis genotypes of motile enteric Eubacterium and Roseburia species were also investigated. Flagellin proteins were recovered from these commensal species, their amino-termini were sequenced and the proteins were found to be pro-inflammatory, as assessed by measurement of interleukin-8 (IL-8) secretion from human intestinal epithelial cell lines. For L. ruminis, this IL-8 secretion required signalling through Toll Like Receptor 5. A model for the regulation of flagellum-biogenesis in L. ruminis was inferred from transcriptomics data and bioinformatics analyses. Motility gene expression in this species may be under the control of a novel regulator, LRC_15730. Potential promoters for genes encoding flagellin proteins in the Eubacterium and Roseburia genomes analysed were inferred in silico. Relative abundances of the target Eubacterium and Roseburia species in the intestinal microbiota of 25 elderly individuals were determined. These species were found to be variably abundant in these individuals. Motility genes from these species were variably detected in the shotgun metagenome databases generated by the ELDERMET project. This suggested that a greater depth of sequencing, or improved evenness of sequencing, would be required to capture the full diversity of microbial functions for specific target or low abundance species in microbial communities by metagenomics. In summary, this thesis used a functional genomics approach to describe flagellum-mediated motility in selected Gram-positive commensal bacteria. The regulation of flagellum biosynthesis in these species, and the consequences of flagella expression from a host-interaction perspective were also considered.

Investigation of toll-like receptor tolerance in the regulation of inflammation

Inflammation is a complex and highly organised immune response to microbes and tissue injury. Recognition of noxious stimuli by pathogen recognition receptor families including Toll-like receptors results in the expression of hundreds of genes that encode cytokines, chemokines, antimicrobials and regulators of inflammation. Regulation of TLR activation responses is controlled by TLR tolerance which induces a global change in the cellular transcriptional expression profile resulting in gene specific suppression and induction of transcription. In this thesis the plasticity of TLR receptor tolerance is investigated using an in vivo, transcriptomics and functional approach to determine the plasticity of TLR tolerance in the regulation of inflammation. Firstly, using mice deficient in the negative regulator of TLR gene transcription, Bcl-3 (Bcl-3-/-) in a model of intestinal inflammation, we investigated the role of Bcl-3 in the regulation of intestinal inflammatory responses. Our data revealed a novel role for Bcl-3 in the regulation of epithelial cell proliferation and regeneration during intestinal inflammation. Furthermore this data revealed that increased Bcl-3 expression contributes to the development of inflammatory bowel disease (IBD). Secondly, we demonstrate that lipopolysaccharide tolerance is transient and recovery from LPS tolerance results in polarisation of macrophages to a previously un-described hybrid state (RM). In addition, we identified that RM cells have a unique transcriptional profile with suppression and induction of genes specific to this polarisation state. Furthermore, using a functional approach to characterise the outcomes of TLR tolerance plasticity, we demonstrate that cytokine transcription is uncoupled from cytokine secretion in macrophages following recovery from LPS tolerance. Here we demonstrate a novel mechanism of regulation of TLR tolerance through suppression of cytokine secretion in macrophages. We show that TNF-α is alternatively trafficked towards a degradative intracellular compartment. These studies demonstrate that TLR tolerance is a complex immunological response with the plasticity of this state playing an important role in the regulation of inflammation.

Sea Urchin Body Plan Development and Evolution: An Integrative Transcriptomic Approach

My dissertation work integrates comparative transcriptomics and functional analyses to investigate gene expression changes underlying two significant aspects of sea urchin evolution and development: the dramatic developmental changes associated with an ecologically significant shift in life history strategy and the development of the unusual radial body plan of adult sea urchins.

In Chapter 2, I investigate evolutionary changes in gene expression underlying the switch from feeding (planktotrophic) to nonfeeding (lecithotrophic) development in sea urchins. In order to identify these changes, I used Illumina RNA-seq to measure expression dynamics across 7 developmental stages in three sea urchin species: the lecithotroph Heliocidaris erythrogramma, the closely related planktotroph Heliocidaris tuberculata, and an outgroup planktotroph Lytechinus variegatus. My analyses draw on a well-characterized developmental gene regulatory network (GRN) in sea urchins to understand how the ancestral planktotrophic developmental program was altered during the evolution of lecithotrophic development. My results suggest that changes in gene expression profiles occurred more frequently across the transcriptome during the evolution of lecithotrophy than during the persistence of planktotrophy. These changes were even more pronounced within the GRN than across the transcriptome as a whole, and occurred in each network territory (skeletogenic, endomesoderm and ectoderm). I found evidence for both conservation and divergence of regulatory interactions in the network, as well as significant changes in the expression of genes with known roles in larval skeletogenesis, which is dramatically altered in lecithotrophs. I further explored network dynamics between species using coexpression analyses, which allowed me to identify novel players likely involved in sea urchin neurogenesis and endoderm patterning.

In Chapter 3, I investigate developmental changes in gene expression underlying radial body plan development and metamorphosis in H. erythrogramma. Using Illumina RNA-seq, I measured gene expression profiles across larval, metamorphic, and post-metamorphic life cycle phases. My results present a high-resolution view of gene expression dynamics during the complex transition from pre- to post-metamorphic development and suggest that distinct sets of regulatory and effector proteins are used during different life history phases.

Collectively, my investigations provide an important foundation for future, empirical studies to investigate the functional role of gene expression change in the evolution of developmental differences between species and also for the generation of the unusual radial body plan of sea urchins.

Bordetella pertussis from functional genomics to intranasal vaccination.

Whooping cough still represents a major health problem, despite the use of effective vaccines for several decades. Being classically a typical childhood disease, whooping cough in young adults is now more common than it used to be, suggesting that protection after vaccination wanes during adolescence. As an alternative to the current vaccines, we wish to develop live attenuated vaccines to be delivered by the nasal route, such as to mimic the natural route of infection and to induce long lasting immunity. Bordetella pertussis, the etiological agent of whooping cough, produces a number of virulence factors, including toxins. Its recently determined genome sequence makes it now possible to apply functional genomics, such as transcriptomics and systematic knock-out mutagenesis. The expression of most known B. pertussis virulence genes is controlled by the two-component system BvgA/S. DNA microarray analyses have led to the identification of novel genes in the BvgA/S regulon, some of which are activated by BvgA/S and others are repressed by BvgA/S. In addition, some genes appear to be differentially modulated by nicotinic acid and MgSO4, both known to modulate the expression of BvgA/S-regulated genes. Among others, the functional genomics approach has uncovered two strongly BvgA/S-activated genes, named hotA and hotB (for 'homolog of toxin'), the products of which show high sequence similarities to pertussis toxin subunits. The identification of the full array of virulence factors, as well as an integrated understanding of the bacterial physiology should allow us to design attenuated B. pertussis strains useful for intranasal vaccination. A first generation of attenuated strains has already shown full protection in mice after a single intranasal administration. Such strains may also serve as vaccine carriers for heterologous antigens, in order to vaccinate against several different pathogens simultaneously.

Transcriptome analysis of sulfate deficiency response in marine microalga Emiliania huxleyi

Summary The response to sulfate deficiency of plants and freshwater green algae has been extensively analysed by system biology approaches. By contrast, seawater sulfate concentration is high and very little is known about the sulfur metabolism of marine organisms. Here, we used a combination of metabolite analysis and transcriptomics to analyse the response of the marine microalga Emiliania huxleyi as it acclimated to sulfate limitation. Lowering sulfate availability in artificial seawater from 25 to 5 mM resulted in significant reduction in growth and intracellular concentrations of dimethylsulfoniopropionate and glutathione. Sulfate-limited E. huxleyi cells showed increased sulfate uptake but sulfate reduction to sulfite did not seem to be regulated. Sulfate limitation in E. huxleyi affected expression of 1718 genes. The vast majority of these genes were upregulated, including genes involved in carbohydrate and lipid metabolism, and genes involved in the general stress response. The acclimation response of E. huxleyi to sulfate deficiency shows several similarities to the well-described responses of Arabidopsis and Chlamydomonas, but also has many unique features. This dataset shows that even though E. huxleyi is adapted to constitutively high sulfate concentration, it retains the ability to re-program its gene expression in response to reduced sulfate availability.

Rice-arsenate interactions in hydroponics: a three-gene model for tolerance

In this study, the genetic mapping of the tolerance of root growth to 13.3 muM arsenate [As(V)] using the BalaxAzucena population is improved, and candidate genes for further study are identified. A remarkable three-gene model of tolerance is advanced, which appears to involve epistatic interaction between three major genes, two on chromosome 6 and one on chromosome 10. Any combination of two of these genes inherited from the tolerant parent leads to the plant having tolerance. Lists of potential positional candidate genes are presented. These are then refined using whole genome transcriptomics data and bioinformatics. Physiological evidence is also provided that genes related to phosphate transport are unlikely to be behind the genetic loci conferring tolerance. These results offer testable hypotheses for genes related to As(V) tolerance that might offer strategies for mitigating arsenic (As) accumulation in consumed rice.

Functional genomics to identify therapeutic prophylactic targets

Infectious diseases are a leading cause of global human mortality. The use of antimicrobials remains the most common strategy for treatment. However, the isolation of pathogens resistant to virtually all antimicrobials makes it urgent to develop effective therapeutics based on new targets. Here we review a new drug discovery paradigm focusing on identifying and targeting host factors important for infection as well as pathogen determinants involved in disease progression. We summarize innovative strategies which by combining bioinformatics with transcriptomics and chemical genetics have already identified host factors essential for pathogen entry, survival and replication. We describe how the discovery of RNA interference which allows loss-of-function studies has facilitated functional genomic studies in human cells. It is expected that these studies will identify targets to be used as host-directed drug therapy which, together with antimicrobials targeting microbial virulence factors, will efficiently eliminate the invading pathogen.