A transcriptomic analysis of striped catfish (Pangasianodon hypophthalmus) in response to salinity adaptation: De novo assembly, gene annotation and marker discovery

The striped catfish (Pangasianodon hypophthalmus) culture industry in the Mekong Delta in Vietnam has developed rapidly over the past decade. The culture industry now however, faces some significant challenges, especially related to climate change impacts notably from predicted extensive saltwater intrusion into many low topographical coastal provinces across the Mekong Delta. This problem highlights a need for development of culture stocks that can tolerate more saline culture environments as a response to expansion of saline water-intruded land. While a traditional artificial selection program can potentially address this need, understanding the genomic basis of salinity tolerance can assist development of more productive culture lines. The current study applied a transcriptomic approach using Ion PGM technology to generate expressed sequence tag (EST) resources from the intestine and swim bladder from striped catfish reared at a salinity level of 9 ppt which showed best growth performance. Total sequence data generated was 467.8 Mbp, consisting of 4,116,424 reads with an average length of 112 bp. De novo assembly was employed that generated 51,188 contigs, and allowed identification of 16,116 putative genes based on the GenBank non-redundant database. GO annotation, KEGG pathway mapping, and functional annotation of the EST sequences recovered with a wide diversity of biological functions and processes. In addition, more than 11,600 simple sequence repeats were also detected. This is the first comprehensive analysis of a striped catfish transcriptome, and provides a valuable genomic resource for future selective breeding programs and functional or evolutionary studies of genes that influence salinity tolerance in this important culture species.

A transcriptomic analysis of the kidney tissue of Tra catfish (pangasianodon hypophthalmus) reared in saline condition: De novo assembly, annotation, SNP discovery

Pangasianodon hypophthalmus is a commercially important freshwater fish used in inland aquaculture in the Mekong Delta, Vietnam. The current study using Ion Torrent technology generated EST resources from the kidney for Tra catfish reared at a salinity level of 9 ppt. We obtained 2,623,929 reads after trimming and processing with an average length of 104 bp. De novo assemblies were generated using CLC Genomic Workbench, Trinity and Velvet/Oases with the best overall contig performance resulting from the CLC assembly. De novo assembly using CLC yielded 29,940 contigs, and allowing identification of 5,710 putative genes when comppared with NCBI non-redundant database. A large number of single nucleotide polymorphisms (SNPs) were also detected. The sequence collection generated in our study represents the most comprehensive transcriptomic resource for P. hypophthalmus available to date.

Genomic and transcriptomic analysis of the streptomycin-dependent Mycobacterium tuberculosis strain 18b.

The ability of Mycobacterium tuberculosis to establish a latent infection (LTBI) in humans confounds the treatment of tuberculosis. Consequently, there is a need to discover new therapeutic agents that can kill M. tuberculosis both during active disease and LTBI. The streptomycin-dependent strain of M. tuberculosis, 18b, provides a useful tool for this purpose since upon removal of streptomycin (STR) it enters a non-replicating state that mimics latency both in vitro and in animal models. The 4.41 Mb genome sequence of M. tuberculosis 18b was determined and this revealed the strain to belong to clade 3 of the ancient ancestral lineage of the Beijing family. STR-dependence was attributable to insertion of a single cytosine in the 530 loop of the 16S rRNA and to a single amino acid insertion in the N-terminal domain of initiation factor 3. RNA-seq was used to understand the genetic programme activated upon STR-withdrawal and hence to gain insight into LTBI. This revealed reconfiguration of gene expression and metabolic pathways showing strong similarities between non-replicating 18b and M. tuberculosis residing within macrophages, and with the core stationary phase and microaerophilic responses. The findings of this investigation confirm the validity of 18b as a model for LTBI, and provide insight into both the evolution of tubercle bacilli and the functioning of the ribosome.

Transcriptomic analysis of rhizobium leguminosarum biovar viciae in symbiosis with host plants pisum sativum and Vicia cracca

Rhizobium leguminosarum bv. viciae forms nitrogen-fixing nodules on several legumes, including pea (Pisum sativum) and vetch (Vicia cracca), and has been widely used as a model to study nodule biochemistry. To understand the complex biochemical and developmental changes undergone by R. leguminosarum bv. viciae during bacteroid development, microarray experiments were first performed with cultured bacteria grown on a variety of carbon substrates (glucose, pyruvate, succinate, inositol, acetate, and acetoacetate) and then compared to bacteroids. Bacteroid metabolism is essentially that of dicarboxylate-grown cells (i.e., induction of dicarboxylate transport, gluconeogenesis and alanine synthesis, and repression of sugar utilization). The decarboxylating arm of the tricarboxylic acid cycle is highly induced, as is gamma-aminobutyrate metabolism, particularly in bacteroids from early (7-day) nodules. To investigate bacteroid development, gene expression in bacteroids was analyzed at 7, 15, and 21 days postinoculation of peas. This revealed that bacterial rRNA isolated from pea, but not vetch, is extensively processed in mature bacteroids. In early development (7 days), there were large changes in the expression of regulators, exported and cell surface molecules, multidrug exporters, and heat and cold shock proteins. fix genes were induced early but continued to increase in mature bacteroids, while nif genes were induced strongly in older bacteroids. Mutation of 37 genes that were strongly upregulated in mature bacteroids revealed that none were essential for nitrogen fixation. However, screening of 3,072 mini-Tn5 mutants on peas revealed previously uncharacterized genes essential for nitrogen fixation. These encoded a potential magnesium transporter, an AAA domain protein, and proteins involved in cytochrome synthesis.

Transcriptomic analysis of Enterohaemorrhagic Escherichia coli O157:H7 in response to plant extracts

Enterohaemorrhagic Escherichia coli (EHEC) are a group of food and contact-borne pathogens responsible for haemorrhagic colitis. The bacteria can be transmitted by contaminated meat, but importantly, also by plants. The bacteria can use plants as an alternative host, where they associate with both the leaves and the roots. Colonisation in the rhizosphere of plants is thought to be the main habitat for colonisation. Four different plant species, commonly associated with EHEC outbreaks, were infected with EHEC O157:H7 isolates Sakai and TUV 93-0 over ten days to assess the colonisation potential of the bacteria in both the phyllosphere and rhizosphere of plants. The rhizosphere was found to sustain a higher population level of bacteria over time in comparison to the phyllosphere, yet both strains were unable to utilize root exudates for growth. Global gene expression changes of EHEC O157:H7 strain Sakai were measured in response to plant extracts such as leaf lysates, root exudates and leaf cell wall polysaccharides from spinach cultivar Amazon and lettuce cultivar Salinas. Microarrays analysis showed a significant change in expression of 17 % of genes on exposure to leaf lysates of spinach. A more specific response was seen to spinach leaf cell wall polysaccharides with only a 1.5 % change. In contrast, when exposed to lettuce leaf cell wall polysaccharides a higher change of 4.8 % was seen. Genes that were differentially expressed belonged to multiple functional groups, including metabolism, indicating the utilization of plant-specific polysaccharides. Several areas of further investigation have been determined from this project, including the importance of culturing bacterial strains at a relevant temperature, the proposed lack of the type III secretion system in plant colonization by EHEC O157:H7 and the utilization of plant components for growth and persistence in the plant environment.

An Integrative Genomic and Transcriptomic Analysis Reveals Potential Targets Associated with Cell Proliferation in Uterine Leiomyomas

Background: Uterine Leiomyomas (ULs) are the most common benign tumours affecting women of reproductive age. ULs represent a major problem in public health, as they are the main indication for hysterectomy. Approximately 40-50% of ULs have non-random cytogenetic abnormalities, and half of ULs may have copy number alterations (CNAs). Gene expression microarrays studies have demonstrated that cell proliferation genes act in response to growth factors and steroids. However, only a few genes mapping to CNAs regions were found to be associated with ULs. Methodology: We applied an integrative analysis using genomic and transcriptomic data to identify the pathways and molecular markers associated with ULs. Fifty-one fresh frozen specimens were evaluated by array CGH (JISTIC) and gene expression microarrays (SAM). The CONEXIC algorithm was applied to integrate the data. Principal Findings: The integrated analysis identified the top 30 significant genes (P<0.01), which comprised genes associated with cancer, whereas the protein-protein interaction analysis indicated a strong association between FANCA and BRCA1. Functional in silico analysis revealed target molecules for drugs involved in cell proliferation, including FGFR1 and IGFBP5. Transcriptional and protein analyses showed that FGFR1 (P = 0.006 and P<0.01, respectively) and IGFBP5 (P = 0.0002 and P = 0.006, respectively) were up-regulated in the tumours when compared with the adjacent normal myometrium. Conclusions: The integrative genomic and transcriptomic approach indicated that FGFR1 and IGFBP5 amplification, as well as the consequent up-regulation of the protein products, plays an important role in the aetiology of ULs and thus provides data for potential drug therapies development to target genes associated with cellular proliferation in ULs. © 2013 Cirilo et al.

Combined snake venomics and venom gland transcriptomic analysis of Bothropoides pauloensis

Unraveling the repertoire of venom toxins of Bothropoides pauloensis was assessed by snake venomics and venom gland transcriptomic surveys. Both approaches yielded converging overall figures, pointing to metalloproteinases (similar to 37%), PLA(2)s (26-32%), and vasoactive (bradykinin-potentiating) peptides (12-17%) as the major toxin classes. The high occurrence of SVMPs, PLA(2) molecules, vasoactive peptides, along with serine proteinases, explains the local and systemic effects observed in envenomations by B. pauloensis. Minor (<3%) C-type lectin, serine proteinase, L-amino acid oxidase, nerve growth factor, and CRISP molecules were also identified in the transcriptome and the proteome. Low abundance (0.3%) EST singletons coding for vascular endothelial growth factor (svVEGF), ohanin, hyaluronidase, and 5' nucleotidase were found only in the venom gland cDNA library. At the molecular level, the transcriptomic and proteomic datasets display low compositional concordance. In particular, although there is good agreement between transcriptome and proteome in the identity of BPPs, PLA(2) molecules and L-amino acid oxidase, both datasets strongly depart in their C-type lectin and SVMP complements. These data support the view that venom composition is influenced by transcriptional and translational mechanisms and emphasize the value of combining proteomic and transcriptomic approaches to acquire a more complete understanding of the toxinological profile and natural history of the snake venom. (C) 2012 Elsevier B.V. All rights reserved.

Effects of nitrogen sources on the nitrate assimilation in Haloferax mediterranei: growth kinetics and transcriptomic analysis

The haloarchaeon Haloferax mediterranei is able to grow in a defined culture media not only in the presence of inorganic nitrogen salt but also with amino acid as the sole nitrogen source. Assimilatory nitrate and nitrite reductases, respectively, catalyze the first and second reactions. The genes involved in this process are nasA, which encodes nitrate reductase and is found within the operon nasABC, and nasD, which encodes nitrite reductase. These genes are subjected to transcriptional regulation, being repressed in the presence of ammonium and induced with either nitrate or nitrite. This type of regulation has also been described when the amino acids are used as nitrogen source in the minimal media. Furthermore, it has been observed that the microorganism growth depends on nitrogen source, obtaining the lowest growth rate in the presence of nitrate and aspartate. In this paper, we present the results of a comparative study of microorganism growth and transcriptomic analysis of the operon nasABC and gene nasD in different nitrogen sources. The results are the first ever produced in relation to amino acids as nitrogen sources within the Halobacteriaceae family.

Transcriptomic Analysis of the Host Response and Innate Resilience to Enterotoxigenic Escherichia coli Infection in Humans.

BACKGROUND: Enterotoxigenic Escherichia coli (ETEC) is a globally prevalent cause of diarrhea. Though usually self-limited, it can be severe and debilitating. Little is known about the host transcriptional response to infection. We report the first gene expression analysis of the human host response to experimental challenge with ETEC. METHODS: We challenged 30 healthy adults with an unattenuated ETEC strain, and collected serial blood samples shortly after inoculation and daily for 8 days. We performed gene expression analysis on whole peripheral blood RNA samples from subjects in whom severe symptoms developed (n = 6) and a subset of those who remained asymptomatic (n = 6) despite shedding. RESULTS: Compared with baseline, symptomatic subjects demonstrated significantly different expression of 406 genes highlighting increased immune response and decreased protein synthesis. Compared with asymptomatic subjects, symptomatic subjects differentially expressed 254 genes primarily associated with immune response. This comparison also revealed 29 genes differentially expressed between groups at baseline, suggesting innate resilience to infection. Drug repositioning analysis identified several drug classes with potential utility in augmenting immune response or mitigating symptoms. CONCLUSIONS: There are statistically significant and biologically plausible differences in host gene expression induced by ETEC infection. Differential baseline expression of some genes may indicate resilience to infection.

Transcriptomic analysis of genes involved in the lipid metabolism of the Alentejano pig

The Alentejano pig is an autochthonous breed scarcely selected, that due to its high trend for fat deposition present poorer meat yields than modern commercial breeds. However, its higher contents of intramuscular fat (IMF) increase pork sensory attributes and consumers’ acceptability. Animal cells can obtain fatty acids (FA) from three distinct pathways: diet ingested fats, lipolysis of stored lipids in cells and through de novo synthesis. Betaine has been used as a dietary supplement in pig nutrition to reduce fat deposition and increase lean muscle mass with inconsistent results so far. This study compares the expression of genes involved in lipid metabolism from pigs consuming a control diet, and the control diet supplemented with betaine (WB). The expression of two genes involved in lipogenesis and lipolysis were evaluated in L. lumborum and B. femoris: ACC, which mediates the carboxylation of acetyl CoA into malonyl CoA concluding the first step of de novo synthesis, and MCPT1 which is responsible for the transport of acyl groups into the mitochondria for the start of β-oxidation.

Different approaches to identify markers for meat quality and production traits in pigs: analysis of the transcriptome in post mortem skeletal muscle and investigation of candidate genes

Due to the growing attention of consumers towards their food, improvement of quality of animal products has become one of the main focus of research. To this aim, the application of modern molecular genetics approaches has been proved extremely useful and effective. This innovative drive includes all livestock species productions, including pork. The Italian pig breeding industry is unique because needs heavy pigs slaughtered at about 160 kg for the production of high quality processed products. For this reason, it requires precise meat quality and carcass characteristics. Two aspects have been considered in this thesis: the application of the transcriptome analysis in post mortem pig muscles as a possible method to evaluate meat quality parameters related to the pre mortem status of the animals, including health, nutrition, welfare, and with potential applications for product traceability (chapters 3 and 4); the study of candidate genes for obesity related traits in order to identify markers associated with fatness in pigs that could be applied to improve carcass quality (chapters 5, 6, and 7). Chapter three addresses the first issue from a methodological point of view. When we considered this issue, it was not obvious that post mortem skeletal muscle could be useful for transcriptomic analysis. Therefore we demonstrated that the quality of RNA extracted from skeletal muscle of pigs sampled at different post mortem intervals (20 minutes, 2 hours, 6 hours, and 24 hours) is good for downstream applications. Degradation occurred starting from 48 h post mortem even if at this time it is still possible to use some RNA products. In the fourth chapter, in order to demonstrate the potential use of RNA obtained up to 24 hours post mortem, we present the results of RNA analysis with the Affymetrix microarray platform that made it possible to assess the level of expression of more of 24000 mRNAs. We did not identify any significant differences between the different post mortem times suggesting that this technique could be applied to retrieve information coming from the transcriptome of skeletal muscle samples not collected just after slaughtering. This study represents the first contribution of this kind applied to pork. In the fifth chapter, we investigated as candidate for fat deposition the TBC1D1 [TBC1 (tre-2/USP6, BUB2, cdc16) gene. This gene is involved in mechanisms regulating energy homeostasis in skeletal muscle and is associated with predisposition to obesity in humans. By resequencing a fragment of the TBC1D1 gene we identified three synonymous mutations localized in exon 2 (g.40A>G, g.151C>T, and g.172T>C) and 2 polymorphisms localized in intron 2 (g.219G>A and g.252G>A). One of these polymorphisms (g.219G>A) was genotyped by high resolution melting (HRM) analysis and PCR-RFLP. Moreover, this gene sequence was mapped by radiation hybrid analysis on porcine chromosome 8. The association study was conducted in 756 performance tested pigs of Italian Large White and Italian Duroc breeds. Significant results were obtained for lean meat content, back fat thickness, visible intermuscular fat and ham weight. In chapter six, a second candidate gene (tribbles homolog 3, TRIB3) is analyzed in a study of association with carcass and meat quality traits. The TRIB3 gene is involved in energy metabolism of skeletal muscle and plays a role as suppressor of adipocyte differentiation. We identified two polymorphisms in the first coding exon of the porcine TRIB3 gene, one is a synonymous SNP (c.132T> C), a second is a missense mutation (c.146C> T, p.P49L). The two polymorphisms appear to be in complete linkage disequilibrium between and within breeds. The in silico analysis of the p.P49L substitution suggests that it might have a functional effect. The association study in about 650 pigs indicates that this marker is associated with back fat thickness in Italian Large White and Italian Duroc breeds in two different experimental designs. This polymorphisms is also associated with lactate content of muscle semimembranosus in Italian Large White pigs. Expression analysis indicated that this gene is transcribed in skeletal muscle and adipose tissue as well as in other tissues. In the seventh chapter, we reported the genotyping results for of 677 SNPs in extreme divergent groups of pigs chosen according to the extreme estimated breeding values for back fat thickness. SNPs were identified by resequencing, literature mining and in silico database mining. analysis, data reported in the literature of 60 candidates genes for obesity. Genotyping was carried out using the GoldenGate (Illumina) platform. Of the analyzed SNPs more that 300 were polymorphic in the genotyped population and had minor allele frequency (MAF) >0.05. Of these SNPs, 65 were associated (P<0.10) with back fat thickness. One of the most significant gene marker was the same TBC1D1 SNPs reported in chapter 5, confirming the role of this gene in fat deposition in pig. These results could be important to better define the pig as a model for human obesity other than for marker assisted selection to improve carcass characteristics.

Early gene expression analysis in 9L orthotopic tumor-bearing rats identifies immune modulation in molecular response to synchrotron microbeam radiation therapy

Synchrotron Microbeam Radiation Therapy (MRT) relies on the spatial fractionation of the synchrotron photon beam into parallel micro-beams applying several hundred of grays in their paths. Several works have reported the therapeutic interest of the radiotherapy modality at preclinical level, but biological mechanisms responsible for the described efficacy are not fully understood to date. The aim of this study was to identify the early transcriptomic responses of normal brain and glioma tissue in rats after MRT irradiation (400Gy). The transcriptomic analysis of similarly irradiated normal brain and tumor tissues was performed 6 hours after irradiation of 9 L orthotopically tumor-bearing rats. Pangenomic analysis revealed 1012 overexpressed and 497 repressed genes in the irradiated contralateral normal tissue and 344 induced and 210 repressed genes in tumor tissue. These genes were grouped in a total of 135 canonical pathways. More than half were common to both tissues with a predominance for immunity or inflammation (64 and 67% of genes for normal and tumor tissues, respectively). Several pathways involving HMGB1, toll-like receptors, C-type lectins and CD36 may serve as a link between biochemical changes triggered by irradiation and inflammation and immunological challenge. Most immune cell populations were involved: macrophages, dendritic cells, natural killer, T and B lymphocytes. Among them, our results highlighted the involvement of Th17 cell population, recently described in tumor. The immune response was regulated by a large network of mediators comprising growth factors, cytokines, lymphokines. In conclusion, early response to MRT is mainly based on inflammation and immunity which appear therefore as major contributors to MRT efficacy.

Endophytic colonization of barley (Hordeum vulgare) roots by the nematophagous fungus Pochonia chlamydosporia reveals plant growth promotion and a general defense and stress transcriptomic response

Plant crop yields are negatively conditioned by a large set of biotic and abiotic factors. An alternative to mitigate these adverse effects is the use of fungal biological control agents and endophytes. The egg-parasitic fungus Pochonia chlamydosporia has been traditionally studied because of its potential as a biological control agent of plant-parasitic nematodes. This fungus can also act as an endophyte in monocot and dicot plants, and has been shown to promote plant growth in different agronomic crops. An Affymetrix 22K Barley GeneChip was used in this work to analyze the barley root transcriptomic response to P. chlamydosporia root colonization. Functional gene ontology (GO) and gene set enrichment analyses showed that genes involved in stress response were enriched in the barley transcriptome under endophytism. An 87.5 % of the probesets identified within the abiotic stress response group encoded heat shock proteins. Additionally, we found in our transcriptomic analysis an up-regulation of genes implicated in the biosynthesis of plant hormones, such as auxin, ethylene and jasmonic acid. Along with these, we detected induction of brassinosteroid insensitive 1-associated receptor kinase 1 (BR1) and other genes related to effector-triggered immunity (ETI) and pattern-triggered immunity (PTI). Our study supports at the molecular level the growth-promoting effect observed in plants endophytically colonized by P. chlamydosporia, which opens the door to further studies addressing the capacity of this fungus to mitigate the negative effects of biotic and abiotic factors on plant crops.

The E2F2 Transcription Factor Sustains Hepatic Glycerophospholipid Homeostasis in Mice

Increasing evidence links metabolic signals to cell proliferation, but the molecular wiring that connects the two core machineries remains largely unknown. E2Fs are master regulators of cellular proliferation. We have recently shown that E2F2 activity facilitates the completion of liver regeneration after partial hepatectomy (PH) by regulating the expression of genes required for S-phase entry. Our study also revealed that E2F2 determines the duration of hepatectomy-induced hepatic steatosis. A transcriptomic analysis of normal adult liver identified "lipid metabolism regulation" as a major E2F2 functional target, suggesting that E2F2 has a role in lipid homeostasis. Here we use wild-type (E2F2(+/+)) and E2F2 deficient (E2F2(-/-)) mice to investigate the in vivo role of E2F2 in the composition of liver lipids and fatty acids in two metabolically different contexts: quiescence and 48-h post-PH, when cellular proliferation and anabolic demands are maximal. We show that liver regeneration is accompanied by large triglyceride and protein increases without changes in total phospholipids both in E2F2(+/+) and E2F2(-/-) mice. Remarkably, we found that the phenotype of quiescent liver tissue from E2F2(-/-) mice resembles the phenotype of proliferating E2F2(+/+) liver tissue, characterized by a decreased phosphatidylcholine to phosphatidylethanolamine ratio and a reprogramming of genes involved in generation of choline and ethanolamine derivatives. The diversity of fatty acids in total lipid, triglycerides and phospholipids was essentially preserved on E2F2 loss both in proliferating and non-proliferating liver tissue, although notable exceptions in inflammation-related fatty acids of defined phospholipid classes were detected. Overall, our results indicate that E2F2 activity sustains the hepatic homeostasis of major membrane glycerolipid components while it is dispensable for storage glycerolipid balance.