Deep sequencing of the zebrafish transcriptome response to mycobacterium infection

Re-identification and characterization of grass carp Ctenopharyngodon idella TLR20

Toll-like receptors (TLRs) play a crucial role in the recognition of microbial-associated molecular patterns in the innate immune system. Fish TLRs have undergone significant gene expansion to adapt to complex aquatic environments. Among them, TLR20 from the TLR11 family actively responds to viral and bacterial invasions. Previous studies have reported two TLR20s in grass carp (Ctenopharyngodon idella), and in this study, we revised this conclusion. Based on the latest grass carp genome, we identified a new TLR20 member. These three TLR20s are arranged in tandem on chromosome 9, indicating that they are generated by gene duplication events. They were renamed CiTLR20.1 to CiTLR20.3 based on their chromosomal positions. The CiTLR20s in C. idella exhibit higher similarities with those in Danio rerio, Cyprinus carpio, and Megalobrama amblycephala, and lower similarities with those in other distantly related fish species. Selective pressure analysis revealed low conservation and negative evolution of TLR20s during evolution. The 3D structures of the three TLR20s showed significant differences, reflecting functional variations and different downstream adaptor molecule recruitment. Transcriptome data revealed tissue distribution differences of TLR20s, with TLR20.1 showing relatively low expression levels in all the tissues, while TLR20.2 and TLR20.3 showed higher expression in the head kidney, spleen, and gill. Additionally, TLR20.2 and TLR20.3 actively responded to GCRV-II infection, with higher upregulation of TLR20.2 in response to Aeromonas hydrophila challenge. In conclusion, this study corrected the number of grass carp TLR20 members and analyzed TLR20 from an evolutionary and structural perspective, exploring its role in antiviral and antibacterial defense. This study provides reference for future research on fish TLR20.

Biosynthetic constraints on amino acid synthesis at the base of the food chain may determine their use in higher-order consumer genomes

Dietary nutrient composition is essential for shaping important fitness traits and behaviours. Many organisms are protein limited, and for Drosophila melanogaster this limitation manifests at the level of the single most limiting essential Amino Acid (AA) in the diet. The identity of this AA and its effects on female fecundity is readily predictable by a procedure called exome matching in which the sum of AAs encoded by a consumer's exome is used to predict the relative proportion of AAs required in its diet. However, the exome matching calculation does not weight AA contributions to the overall profile by protein size or expression. Here, we update the exome matching calculation to include these weightings. Surprisingly, although nearly half of the transcriptome is differentially expressed when comparing male and female flies, we found that creating transcriptome-weighted exome matched diets for each sex did not enhance their fecundity over that supported by exome matching alone. These data indicate that while organisms may require different amounts of dietary protein across conditions, the relative proportion of the constituent AAs remains constant. Interestingly, we also found that exome matched AA profiles are generally conserved across taxa and that the composition of these profiles might be explained by energetic and elemental limitations on microbial AA synthesis. Thus, it appears that ecological constraints amongst autotrophs shape the relative proportion of AAs that are available across trophic levels and that this constrains biomass composition.

Transcriptome analysis of Macrobrachium rosenbergii hemocytes reveals in-depth insights into the immune response to Vibrio parahaemolyticus infection

Macrobrachium rosenbergii as one of the common freshwater prawn species in Southeast Asia, which breeding industry is seriously threatened by vibriosis and causes high mortality. In this study, the RNA-seq was employed for assessing the M. rosenbergii hemocytes transcriptomes following Vibrio parahaemolyticus challenge. After challenge for 6 h (h), there were overall 1849 DEGs or differentially expressed genes, including 1542 up-regulated and 307 down-regulated genes, and there was a total of 1048 DEGs, including 510 up-regulated genes and 538 down-regulated genes, after challenge for 12 h. Mitogen-activated protein kinase (MAPK) immune-related pathways, Toll, immune deficiency (IMD), and Janus kinase (JAK)/signal transducer and activator of transcription (STAT) were among the immune pathways where a lot of the DEGs were connected. The expression patterns of 18 chosen immune-related genes were examined utilizing qRT-PCR or quantitative real-time polymerase chain reaction, which revealed that the V. parahaemolyticus infection activated the M. rosenbergii's immune response. Permutational multivariate analysis of variance (PERMANOVA) showed that V. parahaemolyticus infection modulated immune regulation and apoptosis pathways. The gathered information provided new insight into M. rosenbergii's immunity and suggested a novel approach to fight against bacterial infection.

Comparative Transcriptome Analysis Revealed the Key Genes Regulating Ascorbic Acid Synthesis in Actinidia

Actinidia (kiwifruit) is known as ‘the king of vitamin C’ due to its rich ascorbic acid (AsA) concentration, which makes it an important model for studying the regulation of AsA metabolism. Herein, transcriptomic analysis was employed to identify candidate genes that regulate AsA synthesis in Actinidia species with 100-fold variations in fruit AsA content (A. latifolia and A. rufa). Approximately 1.16 billion high-quality reads were generated, and an average of 66.68% of the data was uniquely aligned against the reference genome. AsA-associated DEGs that predominately respond to abiotic signals, and secondary metabolic pathways were identified. The key candidate genes, for instance, GDP-L-galactose phosphorylase-3 (GGP3), were explored according to integrated analysis of the weighted gene co-expression network and L-galactose pathway. Transgenic kiwifruit plants were generated, and the leaves of GGP3 (OE-GGP3) overexpressing lines had AsA contents 2.0- to 6.4-fold higher than those of the wild type. Transcriptomic analysis of transgenic kiwifruit lines was further implemented to identify 20 potential downstream target genes and understand GGP3-regulated cellular processes. As a result, two transcription factors (AcESE3 and AcMYBR) were selected to carry out yeast two-hybrid and BiFC assays, which verified that there were obvious AcESE3–AcMYBR and AcESE3–AcGGP3 protein–protein interactions. This study provides insight into the mechanism of AsA synthesis and provides candidate factors and genes involved in AsA accumulation in kiwifruit.

De novo transcriptome assembly, polymorphic SSR markers development and population genetics analyses for southern corn rust (Puccinia polysora)

Southern corn rust is a destructive maize disease caused by Puccinia polysora Underw that can lead to severe yield losses. However, genomic information and microsatellite markers are currently unavailable for this disease. In this study, we generated a total of 27,295,216 high-quality cDNA sequence reads using Illumina sequencing technology. These reads were assembled into 17,496 unigenes with an average length of 1015 bp. The functional annotation indicated that 8113 (46.37%), 1933 (11.04%) and 5516 (31.52%) unigenes showed significant similarity to known proteins in the NCBI Nr, Nt and Swiss-Prot databases, respectively. In addition, 2921 (16.70%) unigenes were assigned to KEGG database categories; 4218 (24.11%), to KOG database categories; and 6,603 (37.74%), to GO database categories. Furthermore, we identified 8,798 potential SSRs among 6653 unigenes. A total of 9 polymorphic SSR markers were developed to evaluate the genetic diversity and population structure of 96 isolates collected from Guangdong Province in China. Clonal reproduction of P. polysora in Guangdong was dominant. The YJ (Yangjiang) population had the highest genotypic diversity and the greatest number of the multilocus genotypes, followed by the HY (Heyuan), HZ (Huizhou) and XY (Xinyi) populations. These results provide valuable information for the molecular genetic analysis of P. polysora and related species.

De novo transcriptome analysis and identification of reproduction control genes from the red palm weevil Rhynchophorus ferrugineus

Recent attacks by the red palm weevil, Rhynchophorus ferrugineus (Olivier), have become a severe problem for palm species. In present work, fat body transcriptome of adult female red palm weevil was analyzed, focusing on the identification of reproduction control genes. Transcriptome study was completed by means of next-generation sequencing (NGS) using Illumina Hiseq 2000 sequencing system. A total of 105,938,182 raw reads, 102,645,544 clean reads, and 9,238,098,960 clean nucleotides with a guanine-cytosine content of 40.31%, were produced. The processed transcriptome data resulted in 43,789 unique transcripts (with mean lengths of 1,172 bp). It was found that 20% of total unique transcripts shared up to 80%-100% sequence identity with homologous species, mainly the mountain pine beetle Dendroctonus ponderosae (59.9%) and red flour beetle Tribolium castaneum (26.9%). Nearly 25 annotated genes were predicted to be involved in red palm weevil reproduction, including five vitellogenin (Vg) transcripts. Among the five Vg gene transcripts, one was highly expressed compared with the other four (FPKM values of 1.963, 1.471, 1.028, and 1.017, respectively), and the five Vg gene transcripts were designated as RfVg, RfVg-equivalent1, RfVg-equivalent2, RfVg-equivalent3, and RfVg-equivalent4, respectively. The high expression level of RfVg verified by RT-polymerase chain reaction analysis suggested that RfVg is the primary functional Vg gene in red palm weevil. A high similarity of RfVg with other Coleopterans was also reflected in a phylogenetic tree, where RfVg was placed within the clade of the order Coleoptera. Awareness of the major genes that play critical roles in reproduction and proliferation of red palm weevil is valuable to understand their reproduction mechanism at a molecular level. In addition, for future molecular studies, the NGS dataset obtained will be useful and will promote the exploration of biotech-based control strategies against red palm weevil, a primary pest of palm trees.

The effects of different temperatures on the storage characteristics of lotus (Nelumbo nucifera G.) root

Lotus root (Nelumbo nucifera G.) is a high economic value crop in the world. In this study, the storage characteristics (color, sensory, texture, and fatty acids) of lotus root ("Elian No.5″) were evaluated at different harvest periods (September 2018, October 2018, November 2018, December 2018, and January 2019). Moreover, the storage characteristics were evaluated after the short- term and long-term storage of lotus root at 4 °C and 20 °C. The hardness of lotus root significantly decreased at both temperatures (4 °C and 20 °C) during the first 3 days of storage. In contrast, the decrease in hardness delayed at 4 °C (beyond 3 days of storage). Further, genes related to hardness at different storage temperatures were identified using the RNA-seq and qRT-PCR. The results of this study provide a reference for lotus root storage and a basis for the molecular breeding of longterm-storable lotus root.

RNA-Seq analysis of the guppy immune response against Gyrodactylus bullatarudis infection

Gyrodactylids are ubiquitous ectoparasites of teleost fish, but our understanding of the host immune response against them is fragmentary. Here, we used RNA-Seq to investigate genes involved in the primary response to infection with Gyrodactylus bullatarudis on the skin of guppies, Poecilia reticulata, an important evolutionary model, but also one of the most common fish in the global ornamental trade. Analysis of differentially expressed genes identified several immune-related categories, including IL-17 signalling pathway and Th17 cell differentiation, cytokine-cytokine receptor interaction, chemokine signalling pathway, NOD-like receptor signalling pathway, natural killer cell-mediated cytotoxicity and pathways involved in antigen recognition, processing and presentation. Components of both the innate and the adaptive immune responses play a role in response to gyrodactylid infection. Genes involved in IL-17/Th17 response were particularly enriched among differentially expressed genes, suggesting a significant role for this pathway in fish responses to ectoparasites. Our results revealed a sizable list of genes potentially involved in the teleost-gyrodactylid immune response.

Integrating fish models in tuberculosis vaccine development

Tuberculosis is a chronic infection by Mycobacterium tuberculosis that results in over 1.5 million deaths worldwide each year. Currently, there is only one vaccine against tuberculosis, the Bacillus Calmette–Guérin (BCG) vaccine. Despite widespread vaccination programmes, over 10 million new M. tuberculosis infections are diagnosed yearly, with almost half a million cases caused by antibiotic-resistant strains. Novel vaccination strategies concentrate mainly on replacing BCG or boosting its efficacy and depend on animal models that accurately recapitulate the human disease. However, efforts to produce new vaccines against an M. tuberculosis infection have encountered several challenges, including the complexity of M. tuberculosis pathogenesis and limited knowledge of the protective immune responses. The preclinical evaluation of novel tuberculosis vaccine candidates is also hampered by the lack of an appropriate animal model that could accurately predict the protective effect of vaccines in humans. Here, we review the role of zebrafish (Danio rerio) and other fish models in the development of novel vaccines against tuberculosis and discuss how these models complement the more traditional mammalian models of tuberculosis.

Summary: In this Review, we discuss how zebrafish (Danio rerio) and other fish models can complement the more traditional mammalian models in the development of novel vaccines against tuberculosis.

The combined analysis as the best strategy for Dual RNA-Seq mapping

In Dual RNA-Seq experiments the simultaneous extraction of RNA and analysis of gene expression data from both interacting organisms could be a challenge. One alternative is separating the reads during in silico data analysis. There are two main mapping methods used: sequential and combined. Here we present a combined approach in which the libraries were aligned to a concatenated genome to sort the reads before mapping them to the respective annotated genomes. A comparison of this method with the sequential analysis was performed. Two RNA-Seq libraries available in public databases consisting of a eukaryotic (Zea mays) and a prokaryotic (Herbaspirillum seropediceae) organisms were mixed to simulate a Dual RNA-Seq experiment. Libraries from real Dual RNA-Seq experiments were also used. The sequential analysis consistently attributed more reads to the first reference genome used in the analysis (due to cross-mapping) than the combined approach. More importantly, the combined analysis resulted in lower numbers of cross-mapped reads. Our results highlight the necessity of combining the reference genomes to sort reads previously to the counting step to avoid losing information in Dual RNA-Seq experiments. Since most studies first map the RNA-Seq libraries to the eukaryotic genome, much prokaryotic information has probably been lost.

Characterization of immune response against Mycobacterium marinum infection in the main hematopoietic organ of adult zebrafish (Danio rerio)

Tuberculosis remains a major global health challenge. To gain information about genes important for defense against tuberculosis, we used a well-established tuberculosis model; Mycobacterium marinum infection in adult zebrafish. To characterize the immunological response to mycobacterial infection at 14 days post infection, we performed a whole-genome level transcriptome analysis using cells from kidney, the main hematopoietic organ of adult zebrafish. Among the upregulated genes, those associated with immune signaling and regulation formed the largest category, whereas the largest group of downregulated genes had a metabolic role. We also performed a forward genetic screen in adult zebrafish and identified a fish line with severely impaired survival during chronic mycobacterial infection. Based on transcriptome analysis, these fish have decreased expression of several immunological genes. Taken together, these results give new information about the genes involved in the defense against mycobacterial infection in zebrafish.

Identification of drought response genes by digital gene expression (DGE) analysis in Caragana korshinskii Kom

Caragana korshinskii Kom. is a legume shrub that is widely distributed across desert habitats with gravely, sandy, and saline soils in Asia and Africa. C. korshinskii has highly developed roots and a strong tolerance to abiotic stress. At present, there are few genetic studies of C. korshinskii because of the limited availability of genomic resources. To understand the comprehensive mechanisms that are associated with drought tolerance, we used RNA-seq to survey the differentially expressed genes (DEGs) in comparisons of drought-treated and control plants. After analysing the sequencing results, we found 440 differentially expressed genes existing in drought-treated and control plants. Among the DEGs, 39 unigenes showed up-regulated expression after drought treatment, while 401 unigenes were down-regulated. We used the KEGG database to annotate these drought-induced genes; 126 unigenes were identified by KEGG pathway annotation, and approximately 28% of the unigenes with known function fell into categories related to fatty acid metabolism, starch, sucrose metabolism, and nitrogen metabolism, suggesting that these pathways or processes may be involved in the drought response. Finally, we confirmed that one gene has a potential function in drought tolerance. Our study is the first to provide transcriptomic resources for Caragana korshinskii and to determine its digital gene expression profile under conditions of drought stress using the assembled transcriptomic data for reference. These data provide a valuable resource for genetic and genomic studies of desert plants under abiotic stress conditions.

Piscirickettsia salmonis Cryptic Plasmids: Source of Mobile DNA and Virulence Factors

Four large cryptic plasmids were identified in the salmon pathogen Piscirickettsia salmonis reference strain LF-89. These plasmids appeared highly novel, with less than 7% nucleotidic identity to the nr plasmid database. Plasmid copy number analysis revealed that they are harbored in chromosome equivalent ratios. In addition to plasmid-related genes (plasmidial autonomous replication, partitioning, maintenance, and mobilization genes), mobile genetic elements such as transposases, integrases, and prophage sequences were also identified in P. salmonis plasmids. However, bacterial lysis was not observed upon the induction of prophages. A total of twelve putative virulence factors (VFs) were identified, in addition to two global transcriptional regulators, the widely conserved CsrA protein and the regulator Crp/Fnr. Eleven of the putative VFs were overexpressed during infection in two salmon-derived cellular infection models, supporting their role as VFs. The ubiquity of these plasmids was also confirmed by sequence similarity in the genomes of other P. salmonis strains. The ontology of P. salmonis plasmids suggests a role in bacterial fitness and adaptation to the environment as they encode proteins related to mobilization, nutrient transport and utilization, and bacterial virulence. Further functional characterization of P. salmonis plasmids may improve our knowledge regarding virulence and mobile elements in this intracellular pathogen.

Transcriptome analysis of the almond moth, Cadra cautella, female abdominal tissues and identification of reproduction control genes

Background

The almond moth, Cadra cautella is a destructive pest of stored food commodities including dates that causes severe economic losses for the farming community worldwide. To date, no genetic information related to the molecular mechanism/strategies of its reproduction is available. Thus, transcriptome analysis of C. cautella female abdominal tissues was performed via next-generation sequencing (NGS) to recognize the genes responsible for reproduction.

Results

The NGS was performed with an Illumina Hiseq 2000 sequencer (Beijing Genomics Institute: BGI). From the transcriptome data, 9,804,804,120 nucleotides were generated and their assemblage resulted in 62,687 unigenes. The functional annotation analyses done by different databases, annotated, 27,836 unigenes in total. The transcriptome data of C. cautella female abdominal tissue was submitted to the National Center for Biotechnology Information (accession no: PRJNA484692). The transcriptome analysis yielded several genes responsible for C. cautella reproduction including six Vg gene transcripts. Among the six Vg gene transcripts, only one was highly expressed with 3234.95 FPKM value (fragments per kilobase per million mapped reads) that was much higher than that of the other five transcripts. Higher differences in the expression level of the six Vg transcripts were confirmed by running the RT-PCR using gene specific primers, where the expression was observed only in one transcript it was named as the CcVg.

Conclusions

This is the first study to explore C. cautella reproduction control genes and it might be supportive to explore the reproduction mechanism in this pest at the molecular level. The NGS based transcriptome pool is valuable to study the functional genomics and will support to design biotech-based management strategies for C. cautella.

Transcriptome Analysis of Low-Temperature-Induced Breaking of Garlic Aerial Bulb Dormancy

The long history of asexual reproduction of garlic using garlic cloves has resulted in virus accumulation and genetic depression. Propagation of garlic seedlings by tissue culture can both eliminate viruses and improve breeding efficiency. Aerial bulbs are the first-choice materials for breeding virus-free garlic seedlings under external conditions, but they show dormancy just like garlic bulbs. However, low temperatures can quickly break dormancy. In this research, we used a high-throughput sequencing method to sequence aerial bulbs during dormancy and after low-temperature-induced breaking of dormancy to screen out the key differentially expressed genes (DEGs) associated with low temperature and to provide a theoretical basis for exploring the molecular mechanism of low-temperature-induced breaking of aerial bulb dormancy. The high-throughput transcriptome sequencing results showed that 6,675 DEGs were upregulated and 36,596 DEGs were downregulated in the aerial bulbs given low-temperature treatment. Then, 19,507 DEGs were assigned KEGG annotations, among which most DEGs were annotated to the metabolism pathway (11,817 genes, accounting for 60.58%), followed by the genetic information processing pathway (4,521 genes, accounting for 23.18%). The DEGs were mostly concentrated in pathways such as protein processing in the endoplasmic reticulum, plant-pathogen interaction, plant hormone signal transduction, and ribosome biogenesis in eukaryotes, with significant enrichment. The key DEGs related to calcium signaling, hormonal signaling, and transcription factors were screened out, including CaM, CDPK, and CML in accessory pathways of calcium signaling; GA20ox, GAI1, and GA2ox in accessory pathways of hormonal signaling; and transcription factor genes such as MYB, AP2/ERF, bHLH, MADS, and bZIP. qRT-PCR verification results were consistent with the sequencing results, indicating that the transcriptome sequencing data were accurate and reliable. Our results provide a theoretical basis for breaking the dormancy of aerial bulbs with low-temperature treatment to produce virus-free seedlings and increase the output and quality of garlic.

Transcription Profiles Reveal the Regulatory Synthesis of Phenols during the Development of Lotus Rhizome (Nelumbo nucifera Gaertn)

Lotus (Nelumbo nucifera Gaertn) is a wetland vegetable famous for its nutritional and medicinal value. Phenolic compounds are secondary metabolites that play important roles in the browning of fresh-cut fruits and vegetables, and chemical constituents are extracted from lotus for medicine due to their high antioxidant activity. Studies have explored in depth the changes in phenolic compounds during browning, while little is known about their synthesis during the formation of lotus rhizome. In this study, transcriptomic analyses of six samples were performed during lotus rhizome formation using a high-throughput tag sequencing technique. About 23 million high-quality reads were generated, and 92.14% of the data was mapped to the reference genome. The samples were divided into two stages, and we identified 23,475 genes in total, 689 of which were involved in the biosynthesis of secondary metabolites. A complex genetic crosstalk-regulated network involved in the biosynthesis of phenolic compounds was found during the development of lotus rhizome, and 25 genes in the phenylpropanoid biosynthesis pathway, 18 genes in the pentose phosphate pathway, and 30 genes in the flavonoid biosynthesis pathway were highly expressed. The expression patterns of key enzymes assigned to the synthesis of phenolic compounds were analyzed. Moreover, several differentially expressed genes required for phenolic compound biosynthesis detected by comparative transcriptomic analysis were verified through qRT-PCR. This work lays a foundation for future studies on the molecular mechanisms of phenolic compound biosynthesis during rhizome formation.

Comparative Transcriptome Analysis Reveals a More Complicated Adipogenic Process in Intramuscular Stem Cells than That of Subcutaneous Vascular Stem Cells

Fat-related traits have great influences on pork quality. As different fat tissues have different biochemical profiles depending on their location, intramuscular fat contributes to gustatory qualities, while subcutaneous fat is considered as a negative factor associated with growth performance. In this study, both primary intramuscular and subcutaneous vascular stem cells (IVSCs and SVSCs) could be differentiated into mature adipocytes, though the IVSC differentiation efficiency was lower. By comparative analysis of transcriptomes, 2524 differentially expressed genes (DEGs) were found between two VSCs before differentiation, while only 551 DEGs were found and enriched in two pathways including biosynthesis of unsaturated fatty acids after differentiation. This result indicated that differentiated VSCs were more similar. During differentiation, more DEGs existed in IVSCs than that in SVSCs, suggesting that adipogenesis of IVSCs might be more complex. Additionally, the expression level of DEGs involved in the adipogenic process helps to explain the difference of differentiation efficiency between IVSCs and SVSCs.

Occurrence of diseases in fish used for experimental research

The objective of the present study was to evaluate the occurrence of pathogens and diseases in laboratory fish over a 10-year period at the Centre for Experimental Fish Pathology of Sicily, University of Messina. This report also emphasizes the adverse effects of subclinical infections on research endpoints, as well as the importance of animal health with respect to welfare. Infections in fish used for research can alter experimental outcomes, increase the variability of data, and impede experimental reproducibility. For this purpose, 411 diseased fish of different species (out of a total of 2820 fish) that belonged to four marine species ( Dicentrarchus labrax, Sparus aurata, Argyrosomus regius and Mugil cephalus) and to four fresh water species ( Danio rerio, Carassius auratus, Xiphophorus variatus and Poecilia reticulata) were examined in this study. Our results showed that mycobacteriosis and myxosporidiosis were the most important diseases found in our research fish, and the results represent a useful tool to obtain wider knowledge on the incidence of various diseases in different fish species. Further studies in this field are necessary to improve knowledge on the state of the health of fish used for research.

Transcriptomic profiles of tumor-associated neutrophils reveal prominent roles in enhancing angiogenesis in liver tumorigenesis in zebrafish

We have previously demonstrated the pro-tumoral role of neutrophils using a kras-induced zebrafish hepatocarcinogenesis model. To further illustrate the molecular basis of the pro-tumoral role, Tumor-associated neutrophils (TANs) were isolated by fluorescence-activated cell sorting (FACS) and transcriptomic analyses were carried out by RNA-Seq. Differentially expressed gene profiles of TANs from larvae, male and female livers indicate great variations during liver tumorigenesis, but the common responsive canonical pathways included an immune pathway (Acute Phase Response Signaling), a liver metabolism-related pathway (LXR/RXR Activation) and Thrombin Signaling. Consistent with the pro-tumoral role of TANs, gene module analysis identified a consistent down-regulation of Cytotoxicity module, which may allow continued proliferation of malignant cells. Gene Set Enrichment Analysis indicated up-regulation of several genes promoting angiogenesis. Consistent with this, we found decreased density of blood vessels accompanied with decreased oncogenic liver sizes in neutrophil-depleted larvae. Collectively, our study has indicated some molecular mechanisms of the pro-tumoral roles of TANs in hepatocarcinogenesis, including weakened immune clearance against tumor cells and enhanced function in angiogenesis.

Intelectin 3 is dispensable for resistance against a mycobacterial infection in zebrafish (Danio rerio)

Tuberculosis is a multifactorial bacterial disease, which can be modeled in the zebrafish (Danio rerio). Abdominal cavity infection with Mycobacterium marinum, a close relative of Mycobacterium tuberculosis, leads to a granulomatous disease in adult zebrafish, which replicates the different phases of human tuberculosis, including primary infection, latency and spontaneous reactivation. Here, we have carried out a transcriptional analysis of zebrafish challenged with low-dose of M. marinum, and identified intelectin 3 (itln3) among the highly up-regulated genes. In order to clarify the in vivo significance of Itln3 in immunity, we created nonsense itln3 mutant zebrafish by CRISPR/Cas9 mutagenesis and analyzed the outcome of M. marinum infection in both zebrafish embryos and adult fish. The lack of functional itln3 did not affect survival or the mycobacterial burden in the zebrafish. Furthermore, embryonic survival was not affected when another mycobacterial challenge responsive intelectin, itln1, was silenced using morpholinos either in the WT or itln3 mutant fish. In addition, M. marinum infection in dexamethasone-treated adult zebrafish, which have lowered lymphocyte counts, resulted in similar bacterial burden in both WT fish and homozygous itln3 mutants. Collectively, although itln3 expression is induced upon M. marinum infection in zebrafish, it is dispensable for protective mycobacterial immune response.

De Novo Assembly and Discovery of Genes That Involved in Drought Tolerance in the Common Vetch

The common vetch (Vicia sativa) is often used as feed for livestock because of its high nutritional value. However, drought stress reduces forage production through plant damage. Here, we studied the transcriptional profiles of common vetch exposed to drought in order to understand the molecular mechanisms of drought tolerance in this species. The genome of the common vetch has not been sequenced, therefore we used Illumina sequencing to generate de novo transcriptomes. Nearly 500 million clean reads were used to generate 174,636 transcripts, including 122,299 unigenes. In addition, 5313 transcription factors were identified and these transcription factors were classified into 79 different gene families. We also identified 11,181 SSR loci from di- to hexa-nucleotides whose repeat number was greater than five. On the basis of differentially expressed genes, Gene Ontology analysis identified many drought-relevant categories, including “oxidation-reduction process”, “lipid metabolic process” and “oxidoreductase activity”. In addition to these, Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis identified pathways, such as “Plant hormone signal transduction”, “Glycolysis/Gluconeogenesis” and “Phenylpropanoid biosynthesis”, as differentially expressed in the plants exposed to drought. The expression results in this study will be useful for further extending our knowledge on the drought tolerance of common vetch.

Dual Transcriptomics of Host-Pathogen Interaction of Cystic Fibrosis Isolate Pseudomonas aeruginosa PASS1 With Zebrafish

Pseudomonas aeruginosa is a significant cause of mortality in patients with cystic fibrosis (CF). To explore the interaction of the CF isolate P. aeruginosa PASS1 with the innate immune response, we have used Danio rerio (zebrafish) as an infection model. Confocal laser scanning microscopy (CLSM) enabled visualization of direct interactions between zebrafish macrophages and P. aeruginosa PASS1. Dual RNA-sequencing of host-pathogen was undertaken to profile RNA expression simultaneously in the pathogen and the host during P. aeruginosa infection. Following establishment of infection in zebrafish embryos with PASS1, 3 days post infection (dpi), there were 6739 genes found to be significantly differentially expressed in zebrafish and 176 genes in PASS1. A range of virulence genes were upregulated in PASS1, including genes encoding pyoverdine biosynthesis, flagellin, non-hemolytic phospholipase C, proteases, superoxide dismutase and fimbrial subunits. Additionally, iron and phosphate acquisition genes were upregulated in PASS1 cells in the zebrafish. Transcriptional changes in the host immune response genes highlighted phagocytosis as a key response mechanism to PASS1 infection. Transcriptional regulators of neutrophil and macrophage phagocytosis were upregulated alongside transcriptional regulators governing response to tissue injury, infection, and inflammation. The zebrafish host showed significant downregulation of the ribosomal RNAs and other genes involved in translation, suggesting that protein translation in the host is affected by PASS1 infection.

RNA-Seq-based transcriptome analysis of reproduction- and growth-related genes in Lateolabrax japonicus ovaries at four different ages

Lateolabrax japonicus is an abundant marine aquatic fish species that is commonly cultured in East Asia due to its high commercial value. In this study, RNA-Seq analysis of L. japonicus was carried out to identify reproduction- and growth-related genes expressed in L. japonicus ovaries at different ages using Illumina sequencing technology. In total, 334,388,688 high-quality reads were obtained in four libraries, i.e., 4-year-old ovaries (4th_Ovary), 3-year-old ovaries (3rd_Ovary), 2-year-old ovaries (2nd_Ovary), and 1-year-old ovaries (1st_Ovary). The reads were then de novo assembled into 101,860 unigenes with an average unigene length of 879 bp. In total, 30,142 unigenes (29.59%) were annotated in public databases, including Nr database (Nr), Swiss-Prot, Kyoto Encyclopedia of Genes and Genomes (KEGG), and Cluster of Orthologous Groups of proteins (COG), based on sequence similarity. Comparative analysis revealed that there were 35,749, 43,657, and 36,819 differentially expressed genes (DEGs) in three comparisons (4th_Ovary versus 3rd_Ovary, 4th_Ovary versus 2rd_Ovary, and 4th_Ovary versus 1st_Ovary, respectively). In total, 24,295 DEGs were different expressed in 4th_Ovary. Enrichment and pathway analyses of the DEGs were also carried out to excavate the candidate genes related to reproduction and growth, and 402 genes that potential involved in the regulation of reproduction and growth were identified, e.g., GnRHR (GnRH receptor), GHR 2 (growth hormone receptor 2), I_LGF1R (insulin-like growth factor 1 receptor), etc. Our findings expanded the genomic resources of L. japonicus and provided fundamental information for further studies.

Differential transcriptome analysis of zebrafish (Danio rerio) larvae challenged by Vibrio parahaemolyticus

Zebrafish embryo and larva represent a useful in vivo model for identification of host innate immune responses to bacterial infection. Vibrio parahaemolyticus is a typical zoonotic pathogen worldwide that causes acute gastroenteritis in humans and vibriosis in fishes. However, the mechanism of the innate immune response in the zebrafish larvae infected by V. parahaemolyticus has not been clear. We analysed the transcriptomic profile of 3 days post-fertilization (dpf) zebrafish larvae immersed in V. parahaemolyticus 13 (Vp13) strain suspension for 2 hr. A total of 602 differentially expressed genes (DEGs) were identified in the infection group, of which 175 (29.07%) genes were upregulated and 427 (70.93%) genes were downregulated. These altered genes encoded complement and coagulation cascades, chemokine, TNF signalling pathway, NF-κB signalling pathway and JAK-STAT signalling pathway. Some significant DEGs, such as mmp13, cxcr4a, ccl20, hsp70, gngt, serpina1l, il8, cofilin and il11, were subjected to quantitative gene expression analysis, and the results were consistent with those of the transcriptome profile. These results clearly demonstrated that exposure to V. parahaemolyticus for 2 hr could activate innate immune response in 3dpf larvae by altered expression of downstream signalling pathway genes of pattern recognition receptors (PRRs). Our results also provide a useful reference for future analysis of signal transduction pathways and pathogenesis mechanisms underlying the systemic innate immune response to the external bacteria of V. parahaemolyticus.

Impaired terpenoid backbone biosynthesis reduces saponin accumulation in Panax notoginseng under Cd stress

Panax notoginseng saponins (PNS) are major secondary metabolite of Panax notoginseng (Burk.) F.H. Chen. Previous studies identified that P. notoginseng planting soil usually with high content of Cd. However, the effects of Cd stress on the accumulation of PNS and the corresponding regulation mechanisms have yet to be reported. In the present study, the impact of Cd stress on the PNS accumulation of P. notoginseng was studied in pot culture experiments. The effect of Cd stress on antioxidant enzyme activity was studied using hydroponics. In addition, transcriptase sequencing analysis was used to study the effect of Cd stress on the expression of PNS metabolism transcripts in hydroponic experiments. Cd treatments significantly decreased the accumulation of PNS in the rhizome and main root. The sensitive concentration of antioxidant enzyme activity for both leaf and stem was 2.5μM, whereas the sensitive concentration for the root was 5.0μM. Transcriptome analysis showed that 5132 genes (2930 up- and 2202 downregulated) were regulated by 5.0μM Cd stress in the root of P. notoginseng. Among them, six upregulated differentially expressed genes (DEGs) were related to the methylerythritol 4-phosphate (MEP) pathway, whereas three of the downregulated DEGs were mevalonate kinase (MVK), phosphomevalonate kinase (PMK), and geranylgeranyl diphosphate synthase (type II, GGPS). Of the 15 transcripts selected for real-time quantitative-PCR, 13 were expressed in the same manner as identified using RNA-seq. In conclusion, Cd stress inhibited the accumulation of PNS in the root of P. notoginseng by reducing the expression of MVK, PMK, and GGPS in the terpenoid backbone biosynthesis pathway, and also caused by the removal of reactive oxygen species.

Omics and cytokine discovery in fish: Presenting the Yellowtail kingfish (Seriola lalandi) as a case study

A continued programme of research is essential to overcome production bottlenecks in any aquacultured fish species. Since the introduction of genetic and molecular techniques, the quality of immune research undertaken in fish has greatly improved. Thousands of species specific cytokine genes have been discovered, which can be used to conduct more sensitive studies to understand how fish physiology is affected by aquaculture environments or disease. Newly available transcriptomic technologies, make it increasingly easier to study the immunogenetics of farmed species for which little data exists. This paper reviews how the application of transcriptomic procedures such as RNA Sequencing (RNA-Seq) can advance fish research. As a case study, we present some preliminary findings using RNA-Seq to identify cytokine related genes in Seriola lalandi. These will allow in-depth investigations to understand the immune responses of these fish in response to environmental change or disease and help in the development of therapeutic approaches.

Large-scale identification of Gossypium hirsutum genes associated with Verticillium dahliae by comparative transcriptomic and reverse genetics analysis

Verticillium wilt is a devastating disease of cotton, which is caused by the soil-borne fungus Verticillium dahliae (V. dahliae). Although previous studies have identified some genes or biological processes involved in the interaction between cotton and V. dahliae, its underlying molecular mechanism remains unclear, especially in G. hirsutum. In the present study, we obtained an overview of transcriptome characteristics of resistant upland cotton (G. hirsutum) after V. dahliae infection at 24 h post-inoculation (hpi) via a high-throughput RNA-sequencing technique. A total of 4,794 differentially expressed genes (DEGs) were identified, including 820 up-regulated genes and 3,974 down-regulated genes. The enrichment analysis showed that several important processes were induced upon V. dahliae infection, such as plant hormone signal transduction, plant-pathogen interaction, phenylpropanoid-related and ubiquitin-mediated signals. Moreover, we investigated some key regulatory gene families involved in the defense response, such as receptor-like protein kinases (RLKs), WRKY transcription factors and cytochrome P450 (CYPs), via virus-induced gene silencing (VIGS). GhSKIP35, a partner of SKP1 protein, was involved in ubiquitin-mediated signal. Over-expression of GhSKIP35 in Arabidopsis improved its tolerance to Verticillium wilt in transgenic plants. Collectively, global transcriptome analysis and functional gene characterization provided significant insights into the molecular mechanisms of G. hirsutum-V. dahliae interaction and offered a number of candidate genes as potential sources for breeding wilt-tolerance in cotton.

Characterization of Mycobacterium marinum infections in zebrafish wounds and sinus tracts

The impaired healing of Mycobacterium tuberculosis-infected wounds is a clinical challenge, and the mechanisms involved are still not clear. The zebrafish model of Mycobacterium marinum infection has provided surprising insights into the pathogenesis of tuberculosis in humans. Similarly, the major principles and phases of cutaneous wound healing are conserved among adult mammals and adult zebrafish. Here, we injected Mycobacterium marinum into the dorsal muscles of adult zebrafish and observed the development of chronic wound pathology. Deep sequencing showed that gene expression related to muscles was down-regulated, whereas expressions of the IL-1β, TNF-α, dram1 genes and the transcript of mir1-2 gene were up-regulated in infected wounds of zebrafish compared with control zebrafish. Muscles are immune-responsive tissues. Thus, muscles may play a role in the anti-Mycobacterium tuberculosis immunologic process, which leads to apoptosis of the infected muscle cell and formation of the subcutaneous sinus tract.

Reproductive toxicity of β-diketone antibiotic mixtures to zebrafish (Danio rerio)

So far, few data are available on the reproductive toxicological assessment of β-diketone antibiotics (DKAs), a class of ubiquitous pseudo-persistent pollutant, in zebrafish (Danio rerio). Herein, we reported the reproductive effects of DKAs by means of transcriptome analysis (F1-zebrafish), changes in a series of reproductive indices (F0-zebrafish) and histopathological observations. A total of 1170, 983 and 1399 genes were found to be differentially expressed when compared control vs. 6.25mg/L, control vs. 12.5mg/L and 6.25 vs. 12.5mg/L DKA-exposure treatments, respectively. Among three comparison groups, 670, 569 and 821 genes were respectively assigned for GO analyses based on matches with sequences of known functions. In 149 KEGG-noted metabolic pathways, the preferential one was the MAPK (mitogen-activated protein kinase) signaling pathway, followed by oxidative phosphorylation, neuroactive ligand-receptor interaction and so on. By qPCR verification, 6 genes (c6ast4, igfbp1b, mrpl42, tnnc2, emc4 and ddit4) showed consistent gene expression with those identified by transcriptome sequencing. Due to DKA-exposure, the concentrations of plasma estradiol and testosterone, and the gonado-somatic index were significantly dose-dependently declined. Also, DKA-exposure led to declining in zebrafish reproductive capacity, reflecting in fertilization, hatchability and egg production. Histopathological observations demonstrated that zebrafish ovary and testis suffered serious damage after DKA-exposure. The 4-oxo-TEMP signals increased obviously with increasing DKA-exposed concentrations, implying disruption of balance between generation and clearance of ¹O₂. In summary, DKAs not only produce reproductive toxicological effects on F0-zebrafish, but also result in adverse consequences for growth and development of F1-zebrafish.

Transcriptome profiling and digital gene expression analysis of the skin of Dybowski's frog (Rana dybowskii) exposed to Aeromonas hydrophila

Recently, populations of Rana dybowskii, an important amphibian species in Northeast China, have decreased, mainly owing to the disease caused by Aeromonas hydrophila. However, effective control methods have not yet been developed. In order to explore the immune responses of R. dybowskii upon exposure to A. hydrophila infection, Illumina high-throughput transcriptome sequencing and digital gene expression (DGE) technology were employed to investigate transcriptomic changes in the skin of R. dybowskii exposed to A. hydrophila. In this work, a total of 26,244,446 transcriptome sequencing reads were obtained and assembled into 109,089 unique unigenes using de novo assembly, and a total of 37,105 unigenes (34.0%) were functionally annotated against the non-redundant (Nr), Swiss-Prot, Cluster of Orthologous Groups of Proteins (COG), Kyoto Encyclopedia of Genes and Genomes (KEGG), and Gene Ontology (GO) databases. Gene expression changes in the skin tissue of R. dybowskii exposed to A. hydrophila were investigated by a tag-based DGE system, and a total of 1435 significantly differentially expressed genes were identified, including 460 that were up-regulated and 975 that were down-regulated, indicating a large change in the host transcriptome profile exposed to A. hydrophila. Among these, 478 genes were associated with immune-relevant pathways, metabolic pathways, cellular components, growth, migration, and muscle and hormone signaling pathways. We confirmed the differential expression of 106 immune-relevant genes associated with innate and adaptive immune responses. Our data provide a fairly comprehensive molecular biology background for the deeper understanding of the amphibian immune system following A. hydrophila infection.

Transcriptomic responses in the fish intestine

The intestine, being a multifunctional organ central to both nutrient uptake, pathogen recognition and regulating the intestinal microbiome, has been subjected to intense research. This review will focus on the recent studies carried out using high-throughput gene expression approaches, such as microarray and RNA sequencing (RNA-seq). These techniques have advanced greatly in recent years, mainly as a result of the massive changes in sequencing methodologies. At the time of writing, there is a transition between relatively well characterised microarray platforms and the developing RNA-seq, with the prediction that within a few years as costs decrease and computation power increase, RNA-seq related approaches will supersede the microarrays. Comparisons between the approaches are made and specific examples of how the techniques have been used to examine intestinal responses to pathogens, dietary manipulations and osmoregulatory challenges are given.

Chronic toxicological effects of β-diketone antibiotics on Zebrafish (Danio rerio) using transcriptome profiling of deep sequencing

Transcriptome analysis is important for interpreting the functional elements of the genome and revealing the molecular constituents of cells and tissues. Herein, differentially transcribed genes were identified by deep sequencing after zebrafish (Danio rerio) were exposed to β-diketone antibiotics (DKAs); 23,129 and 23,550 mapped genes were detected in control and treatment groups, a total of 3238 genes were differentially expressed between control and treatment groups. Of these genes, 328 genes (213 up- and 115 down-regulation) had significant differential expression (p < 0.05) and an expression ratio (control/treatment) of >2 or <0.5. Additionally, we performed Gene Ontology (GO) category and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analyses, and found 266 genes in the treatment group with annotation terms linked to the GO category. A total of 77 differentially expressed transcriptional genes were associated with 132 predicted KEGG metabolic pathways. Serious liver tissue damage was reflected and consistent with the differences in genetic classification and function from the transcriptome analysis. These results enhance our understanding of zebrafish developmental processes under exposure to DKA stress. © 2015 Wiley Periodicals, Inc. Environ Toxicol 31: 1357-1371, 2016.

De Novo Assembly and Transcriptome Analysis of Bulb Onion (Allium cepa L.) during Cold Acclimation Using Contrasting Genotypes

Bulb onion (Allium cepa) is the second most widely cultivated and consumed vegetable crop in the world. During winter, cold injury can limit the production of bulb onion. Genomic resources available for bulb onion are still very limited. To date, no studies on heritably durable cold and freezing tolerance have been carried out in bulb onion genotypes. We applied high-throughput sequencing technology to cold (2°C), freezing (-5 and -15°C), and control (25°C)-treated samples of cold tolerant (CT) and cold susceptible (CS) genotypes of A. cepa lines. A total of 452 million paired-end reads were de novo assembled into 54,047 genes with an average length of 1,331 bp. Based on similarity searches, these genes were aligned with entries in the public non-redundant (nr) database, as well as KEGG and COG database. Differentially expressed genes (DEGs) were identified using log₁₀ values with the FPKM method. Among 5,167DEGs, 491 genes were differentially expressed at freezing temperature compared to the control temperature in both CT and CS libraries. The DEG results were validated with qRT-PCR. We performed GO and KEGG pathway enrichment analyses of all DEGs and iPath interactive analysis found 31 pathways including those related to metabolism of carbohydrate, nucleotide, energy, cofactors and vitamins, other amino acids and xenobiotics biodegradation. Furthermore, a large number of molecular markers were identified from the assembled genes, including simple sequence repeats (SSRs) 4,437 and SNP substitutions of transition and transversion types of CT and CS. Our study is the first to provide a transcriptome sequence resource for Allium spp. with regard to cold and freezing stress. We identified a large set of genes and determined their DEG profiles under cold and freezing conditions using two different genotypes. These data represent a valuable resource for genetic and genomic studies of Allium spp.

De Novo Transcriptome Analysis of Differential Functional Gene Expression in Largemouth Bass (Micropterus salmoides) after Challenge with Nocardia seriolae

Largemouth bass (Micropterus salmoides) are common hosts of an epizootic bacterial infection by Nocardia seriolae. We conducted transcriptome profiling of M. salmoides to understand the host immune response to N. seriolae infection, using the Illumina sequencing platform. De novo assembly of paired-end reads yielded 47,881 unigenes, the total length, average length, N50, and GC content of which were 49,734,288, 1038, 1983 bp, and 45.94%, respectively. Annotation was performed by comparison against non-redundant protein sequence (NR), non-redundant nucleotide (NT), Swiss-Prot, Clusters of Orthologous Groups (COG), Kyoto Encyclopaedia of Genes and Genomes (KEGG), Gene Ontology (GO), and Interpro databases, yielding 28,964 (NR: 60.49%), 36,686 (NT: 76.62%), 24,830 (Swissprot: 51.86%), 8913 (COG: 18.61%), 20,329 (KEGG: 42.46%), 835 (GO: 1.74%), and 22,194 (Interpro: 46.35%) unigenes. Additionally, 8913 unigenes were classified into 25 Clusters of Orthologous Groups (KOGs) categories, and 20,329 unigenes were assigned to 244 specific signalling pathways. RNA-Seq by Expectation Maximization (RSEM) and PossionDis were used to determine significantly differentially expressed genes (False Discovery Rate (FDR) < 0.05) and we found that 1384 were upregulated genes and 1542 were downregulated genes, and further confirmed their regulations using reverse transcription quantitative PCR (RT-qPCR). Altogether, these results provide information on immune mechanisms induced during bacterial infection in largemouth bass, which may facilitate the prevention of nocardiosis.

Comparative Digital Gene Expression Analysis of the Arabidopsis Response to Volatiles Emitted by Bacillus amyloliquefaciens

Some plant growth-promoting rhizobacteria (PGPR) regulated plant growth and elicited plant basal immunity by volatiles. The response mechanism to the Bacillus amyloliquefaciens volatiles in plant has not been well studied. We conducted global gene expression profiling in Arabidopsis after treatment with Bacillus amyloliquefaciens FZB42 volatiles by Illumina Digital Gene Expression (DGE) profiling of different growth stages (seedling and mature) and tissues (leaves and roots). Compared with the control, 1,507 and 820 differentially expressed genes (DEGs) were identified in leaves and roots at the seedling stage, respectively, while 1,512 and 367 DEGs were identified in leaves and roots at the mature stage. Seventeen genes with different regulatory patterns were validated using quantitative RT-PCR. Numerous DEGs were enriched for plant hormones, cell wall modifications, and protection against stress situations, which suggests that volatiles have effects on plant growth and immunity. Moreover, analyzes of transcriptome difference in tissues and growth stage using DGE profiling showed that the plant response might be tissue-specific and/or growth stage-specific. Thus, genes encoding flavonoid biosynthesis were downregulated in leaves and upregulated in roots, thereby indicating tissue-specific responses to volatiles. Genes related to photosynthesis were downregulated at the seedling stage and upregulated at the mature stage, respectively, thereby suggesting growth period-specific responses. In addition, the emission of bacterial volatiles significantly induced killing of cells of other organism pathway with up-regulated genes in leaves and the other three pathways (defense response to nematode, cell morphogenesis involved in differentiation and trichoblast differentiation) with up-regulated genes were significantly enriched in roots. Interestingly, some important alterations in the expression of growth-related genes, metabolic pathways, defense response to biotic stress and hormone-related genes were firstly founded response to FZB42 volatiles.

Transcriptome analysis of genes responding to NNV infection in Asian seabass epithelial cells

Asian seabass is an important food fish in Southeast Asia. Viral nervous necrosis (VNN) disease, triggered by nervous necrosis virus (NNV) infection, has caused mass mortality of Asian seabass larvae, resulting in enormous economic losses in the Asian seabass industry. In order to better understand the complex molecular interaction between Asian seabass and NNV, we investigated the transcriptome profiles of Asian seabass epithelial cells, which play an essential role in immune regulation, after NNV infection. Using the next generation sequencing (NGS) technology, we sequenced mRNA from eight samples (6, 12, 24, 48 h post-inoculation) of mock and NNV-infected Asian seabass epithelial cell line, respectively. Clean reads were de novo assembled into a transcriptome consisting of 89026 transcripts with a N50 of 2617 bp. Furthermore, 251 differentially expressed genes (DEGs) in response to NNV infection were identified. Top DEGs include protein asteroid homolog 1-like (ASTE1), receptor-transporting protein 3 (RTP3), heat shock proteins 30 (HSP30) and 70 (HSP70), Viperin, interferon regulatory factor 3 (IRF3) and other genes related to innate immunity. Our data suggest that abundant and diverse genes corresponding to NNV infection. The results of this study could also offer vital information not only for identification of novel genes involved in Asian seabass-NNV interaction, but also for our understanding of the molecular mechanism of Asian seabass' response to viral infection. In addition, 24807 simple sequence repeats (SSRs) were detected in the assembled transcriptome, providing valuable resources for studying genetic variations and accelerating quantitative trait loci (QTL) mapping for disease resistance in Asian seabass in the future.

Transcriptomic Approaches in the Zebrafish Model for Tuberculosis-Insights Into Host- and Pathogen-specific Determinants of the Innate Immune Response

Mycobacterium marinum infection in zebrafish has become a well-established model of tuberculosis. Both embryonic and adult zebrafish infection studies have contributed to our knowledge of the development and function of tuberculous granulomas, which are typical of mycobacterial pathogenesis. In this review we discuss how transcriptome profiling studies have helped to characterize this infection process. We illustrate this using new RNA sequencing (RNA-Seq) data that reveals three main phases in the host response to M. marinum during the early stages of granuloma development in zebrafish embryos and larvae. The early phase shows induction of complement and transcription factors, followed by a relatively minor induction of pro-inflammatory cytokines within hours following phagocytosis of M. marinum. A minimal response is observed in the mid-phase, between 6 hours and 1day post infection, when the tissue dissemination of M. marinum begins. During subsequent larval development the granulomas expand and a late-phase response is apparent, which is characterized by progressively increasing induction of complement, transcription factors, pro-inflammatory cytokines, matrix metalloproteinases, and other defense and inflammation-related gene groups. This late-phase response shares common components with the strong and acute host transcriptome response that has previously been reported for Salmonella typhimurium infection in zebrafish embryos. In contrast, the early/mid-phase response to M. marinum infection, characterized by suppressed pro-inflammatory signaling, is strikingly different from the acute response to S. typhimurium infection. Furthermore, M. marinum infection shows a collective and strongly fluctuating regulation of lipoproteins, while S. typhimurium infection has pronounced effects on amino acid metabolism and glycolysis.

Different RNA splicing mechanisms contribute to diverse infective outcome of classical swine fever viruses of differing virulence: insights from the deep sequencing data in swine umbilical vein endothelial cells

Molecular mechanisms underlying RNA splicing regulation in response to viral infection are poorly understood. Classical swine fever (CSF), one of the most economically important and highly contagious swine diseases worldwide, is caused by classical swine fever virus (CSFV). Here, we used high-throughput sequencing to obtain the digital gene expression (DGE) profile in swine umbilical vein endothelial cells (SUVEC) to identify different response genes for CSFV by using both Shimen and C strains. The numbers of clean tags obtained from the libraries of the control and both CSFV-infected libraries were 3,473,370, 3,498,355, and 3,327,493 respectively. In the comparison among the control, CSFV-C, and CSFV-Shimen groups, 644, 158, and 677 differentially expressed genes (DEGs) were confirmed in the three groups. Pathway enrichment analysis showed that many of these DEGs were enriched in spliceosome, ribosome, proteasome, ubiquitin-mediated proteolysis, cell cycle, focal adhesion, Wnt signalling pathway, etc., where the processes differ between CSFV strains of differing virulence. To further elucidate important mechanisms related to the differential infection by the CSFV Shimen and C strains, we identified four possible profiles to assess the significantly expressed genes only by CSFV Shimen or CSFV C strain. GO analysis showed that infection with CSFV Shimen and C strains disturbed ‘RNA splicing’ of SUVEC, resulting in differential ‘gene expression’ in SUVEC. Mammalian target of rapamycin (mTOR) was identified as a significant response regulator contributed to impact on SUVEC function for CSFV Shimen. This computational study suggests that CSFV of differing virulence could induce alterations in RNA splicing regulation in the host cell to change cell metabolism, resulting in acute haemorrhage and pathological damage or infectious tolerance.

Immunotoxicity of β-Diketone Antibiotic Mixtures to Zebrafish (Danio rerio) by Transcriptome Analysis

Fluoroquinolones and tetracyclines are known as β-diketone antibiotics (DKAs) because of bearing a diketone group in their molecular structure. DKAs are the most widely used antibiotics to prevent generation of disease in humans and animals and to suppress bacterial growth in aquaculture. In recent years, overuse of DKAs has caused serious environmental risk due to their pseudo-persistence in the environment, even though their half-lives are not long. So far, no reports were concerned with the joint immunotoxicity of DKAs. Herein, we reported on the immunotoxicity of DKAs on zebrafish after a 3-month DKAs exposure using transcriptomic techniques. According to transcriptome sequencing, 10 differentially expressed genes were screened out among the genes related to KEGG pathways with high enrichment. The identified 7 genes showed to be consistent between RNA-seq and qRT-PCR. Due to DKAs exposure, the content or activity for a series of immune-related biomarkers (Complement 3, lysozyme, IgM and AKP) showed the inconsistent changing trends as compared with the control group. Histopathological observations showed that the number of goblet cells increased sharply, the columnar epithelial cells swelled, the nucleus became slender in intestinal villi, and numerous brown metachromatic granules occurred in spleens of DKAs-exposed groups. Overall, both detection of biomarkers and histopathological observation corroborated that chronic DKAs exposure could result in abnormal expression of immune genes and enzymes, and variable levels of damage to immune-related organs. These complex effects of DKAs may lead to zebrafish dysfunction and occurrence of diseases related to the immune system.

Genome-Wide Host-Pathogen Interaction Unveiled by Transcriptomic Response of Diamondback Moth to Fungal Infection

Genome-wide insight into insect pest response to the infection of Beauveria bassiana (fungal insect pathogen) is critical for genetic improvement of fungal insecticides but has been poorly explored. We constructed three pairs of transcriptomes of Plutella xylostella larvae at 24, 36 and 48 hours post treatment of infection (hpt_I) and of control (hpt_C) for insight into the host-pathogen interaction at genomic level. There were 2143, 3200 and 2967 host genes differentially expressed at 24, 36 and 48 hpt_I/hpt_C respectively. These infection-responsive genes (~15% of the host genome) were enriched in various immune processes, such as complement and coagulation cascades, protein digestion and absorption, and drug metabolism-cytochrome P450. Fungal penetration into cuticle and host defense reaction began at 24 hpt_I, followed by most intensive host immune response at 36 hpt_I and attenuated immunity at 48 hpt_I. Contrastingly, 44% of fungal genes were differentially expressed in the infection course and enriched in several biological processes, such as antioxidant activity, peroxidase activity and proteolysis. There were 1636 fungal genes co-expressed during 24–48 hpt_I, including 116 encoding putative secretion proteins. Our results provide novel insights into the insect-pathogen interaction and help to probe molecular mechanisms involved in the fungal infection to the global pest.

Transcriptome profiling of grass carp (Ctenopharyngodon idellus) infected with Aeromonas hydrophila

Aeromonas hydrophila is the causative pathogen of intestinal hemorrhage which has caused great economic loss in grass carp aquaculture. In order to understand the immunological response of grass carp to infection by A. hydrophila, the transcriptomic profiles of the spleens from infected and non-infected grass carp groups were obtained using HiSeq™ 2500 (Illumina). An average of 63 million clean reads per library was obtained, and approximately 80% of these genes were successfully mapped to the reference genome. A total of 1591 up-regulated and 530 down-regulated genes were identified. Eight immune-related categories involving 105 differently expressed genes were scrutinized. 16 of the differently expressed genes involving immune response were further validated by qRT-PCR. Our results provide valuable information for further analysis of the mechanisms of grass carp defense against A. hydrophila invasion.

Co-transcriptomic Analysis by RNA Sequencing to Simultaneously Measure Regulated Gene Expression in Host and Bacterial Pathogen

Intramacrophage pathogens subvert antimicrobial defence pathways using various mechanisms, including the targeting of host TLR-mediated transcriptional responses. Conversely, TLR-inducible host defence mechanisms subject intramacrophage pathogens to stress, thus altering pathogen gene expression programs. Important biological insights can thus be gained through the analysis of gene expression changes in both the host and the pathogen during an infection. Traditionally, research methods have involved the use of qPCR, microarrays and/or RNA sequencing to identify transcriptional changes in either the host or the pathogen. Here we describe the application of RNA sequencing using samples obtained from in vitro infection assays to simultaneously quantify both host and bacterial pathogen gene expression changes, as well as general approaches that can be undertaken to interpret the RNA sequencing data that is generated. These methods can be used to provide insights into host TLR-regulated transcriptional responses to microbial challenge, as well as pathogen subversion mechanisms against such responses.

Transcriptome and Gene Expression Analysis of Cylas formicarius (Coleoptera: Brentidae) During Different Development Stages

The sweet potato weevil, Cylas formicarius (F.) (Coleoptera: Brentidae), is an important pest of sweet potato worldwide. However, there is limited knowledge on the molecular mechanisms underlying growth and differentiation of C. formicarius The transcriptomes of the eggs, second instar larvae, third instar larvae (L3), pupae, females, and males of C. formicarius were sequenced using Illumina sequencing technology for obtaining global insights into developing transcriptome characteristics and elucidating the relative functional genes. A total of 54,255,544 high-quality reads were produced, trimmed, and de novo assembled into 115,281 contigs. 61,686 unigenes were obtained, with an average length of 1,009 nt. Among these unigenes, 17,348 were annotated into 59 Gene Ontology (GO) terms and 12,660 were assigned to 25 Cluster of Orthologous Groups classes, whereas 24,796 unigenes were mapped to 258 pathways. Differentially expressed unigenes between various developmental stages of C. formicarius were detected. Higher numbers of differentially expressed genes (DEGs) were recorded in the eggs versus L3 and eggs versus male samples (2,141 and 2,058 unigenes, respectively) than the others. Genes preferentially expressed in each stage were also identified. GO and pathway-based enrichment analysis were used to further investigate the functions of the DEGs. In addition, the expression profiles of ten DEGs were validated by quantitative real-time PCR. The transcriptome profiles presented in this study and these DEGs detected by comparative analysis of different developed stages of C. formicarius will facilitate the understanding of the molecular mechanism of various living process and will contribute to further genome-wide research.

Ixodes scapularis and Ixodes ricinus tick cell lines respond to infection with tick-borne encephalitis virus: transcriptomic and proteomic analysis

BACKGROUND

Ixodid ticks are important vectors of a wide variety of viral, bacterial and protozoan pathogens of medical and veterinary importance. Although several studies have elucidated tick responses to bacteria, little is known about the tick response to viruses. To gain insight into the response of tick cells to flavivirus infection, the transcriptomes and proteomes of two Ixodes spp cell lines infected with the flavivirus tick-borne encephalitis virus (TBEV) were analysed.

METHODS

RNA and proteins were isolated from the Ixodes scapularis-derived cell line IDE8 and the Ixodes ricinus-derived cell line IRE/CTVM19, mock-infected or infected with TBEV, on day 2 post-infection (p.i.) when virus production was increasing, and on day 6 p.i. when virus production was decreasing. RNA-Seq and mass spectrometric technologies were used to identify changes in abundance of, respectively, transcripts and proteins. Functional analyses were conducted on selected transcripts using RNA interference (RNAi) for gene knockdown in tick cells infected with the closely-related but less pathogenic flavivirus Langat virus (LGTV).

RESULTS

Differential expression analysis using DESeq resulted in totals of 43 and 83 statistically significantly differentially-expressed transcripts in IDE8 and IRE/CTVM19 cells, respectively. Mass spectrometry detected 76 and 129 statistically significantly differentially-represented proteins in IDE8 and IRE/CTVM19 cells, respectively. Differentially-expressed transcripts and differentially-represented proteins included some that may be involved in innate immune and cell stress responses. Knockdown of the heat-shock proteins HSP90, HSP70 and gp96, the complement-associated protein Factor H and the protease trypsin resulted in increased LGTV replication and production in at least one tick cell line, indicating a possible antiviral role for these proteins. Knockdown of RNAi-associated proteins Argonaute and Dicer, which were included as positive controls, also resulted in increased LGTV replication and production in both cell lines, confirming their role in the antiviral RNAi pathway.

CONCLUSIONS

This systems biology approach identified several molecules that may be involved in the tick cell innate immune response against flaviviruses and highlighted that ticks, in common with other invertebrate species, have other antiviral responses in addition to RNAi.

Differentially-Expressed Genes Associated with Faster Growth of the Pacific Abalone, Haliotis discus hannai

The Pacific abalone Haliotis discus hannai is used for commercial aquaculture in Korea. We examined the transcriptome of Pacific abalone Haliotis discus hannai siblings using NGS technology to identify genes associated with high growth rates. Pacific abalones grown for 200 days post-fertilization were divided into small-, medium-, and large-size groups with mean weights of 0.26 ± 0.09 g, 1.43 ± 0.405 g, and 5.24 ± 1.09 g, respectively. RNA isolated from the soft tissues of each group was subjected to RNA sequencing. Approximately 1%–3% of the transcripts were differentially expressed in abalones, depending on the growth rate. RT-PCR was carried out on thirty four genes selected to confirm the relative differences in expression detected by RNA sequencing. Six differentially-expressed genes were identified as associated with faster growth of the Pacific abalone. These include five up-regulated genes (including one specific to females) encoding transcripts homologous to incilarin A, perlucin, transforming growth factor-beta-induced protein immunoglobulin-heavy chain 3 (ig-h3), vitelline envelope zona pellucida domain 4, and defensin, and one down-regulated gene encoding tomoregulin in large abalones. Most of the transcripts were expressed predominantly in the hepatopancreas. The genes identified in this study will lead to development of markers for identification of high-growth-rate abalones and female abalones.

Disparate developmental patterns of immune responses to bacterial and viral infections in fish

During early stages of development vertebrates rely on an immature immune system to fight pathogens, but in non mammalian species few studies have taken an in-depth analysis of the transition from reliance on innate immune mechanisms to the appearance of adaptive immunity. Using rainbow trout as a model we characterized responses to two natural pathogens of this species, the Gram negative bacterium Aeromonas salmonicida and the virus VHSV, using microarray analysis at four early life history stages; eyed egg, post hatch, first feeding and three weeks post first feeding when adaptive immunity starts to be effective. All stages responded to both infections, but the complexity of the response increased with developmental stage. The response to virus showed a clear interferon response only from first feeding. In contrast, bacterial infection induced a marked response from early stages, with modulation of inflammatory, antimicrobial peptide and complement genes across all developmental stages. Whilst the viral and bacterial responses were distinct, there were modulated genes in common, mainly of general inflammatory molecules. This work provides a first platform to explore the development of fish immunity to infection, and to compare the age-dependent changes (from embryo to adults) across vertebrates.

Neutrophils in host defense: new insights from zebrafish

Neutrophils are highly motile phagocytic cells that play a critical role in the immune response to infection. Zebrafish (Danio rerio) are increasingly used to study neutrophil function and host-pathogen interactions. The generation of transgenic zebrafish lines with fluorescently labeled leukocytes has made it possible to visualize the neutrophil response to infection in real time by use of optically transparent zebrafish larvae. In addition, the genetic tractability of zebrafish has allowed for the generation of models of inherited neutrophil disorders. In this review, we discuss several zebrafish models of infectious disease, both in the context of immunocompetent, as well as neutrophil-deficient hosts and how these models have shed light on neutrophil behavior during infection.

Deep sequencing analysis of transcriptomes in Aspergillus flavus in response to resveratrol

Background

Resveratrol has been reported as a natural phytoalexin that inhibits infection or the growth of certain fungi including Aspergillus flavus. Our previous research revealed that aflatoxin production in A. flavus was reduced in medium with resveratrol. To understand the molecular mechanism of the A. flavus response to resveratrol treatment, the high-throughput paired-end RNA-Seq was applied to analyze the transcriptomic profiles of A. flavus.

Results

In total, 366 and 87 genes of A. flavus were significantly up- and down- regulated, respectively, when the fungus was treated with resveratrol. Gene Ontology (GO) functional enrichment analysis revealed that 48 significantly differentially expressed genes were involved in 6 different terms. Most genes in the aflatoxin biosynthetic pathway genes cluster (#54) did not show a significant change when A. flavus was treated with resveratrol, but 23 of the 30 genes in the #54 cluster were down-regulated. The transcription of aflA and aflB was significantly suppressed under resveratrol treatment, resulting in an insufficient amount of the starter unit hexanoate for aflatoxin biosynthesis. In addition, resveratrol significantly increased the activity of antioxidative enzymes that destroy radicals, leading to decreased aflatoxin production. Moreover, stuA, fluG, flbC, and others genes involved in mycelial and conidial development were down-regulated, which disrupted the cell’s orderly differentiation and blocked conidia formation and mycelia development. The transcripts of laeA and veA were slightly inhibited by resveratrol, which may partly decrease aflatoxin production and depress conidia formation.

Conclusions

Resveratrol can affect the expression of A. flavus genes that are related to developmental and secondary metabolic processes, resulting in decreased aflatoxin production and conidia formation and could also cause abnormal mycelia development. These results provide insight into the transcriptome of A. flavus in response to resveratrol and a new clew for further study in regulation of aflatoxin biosynthesis in A. flavus.

Electronic supplementary material

The online version of this article (doi:10.1186/s12866-015-0513-6) contains supplementary material, which is available to authorized users.

Transcriptome analysis of the pearl oyster (Pinctada fucata) hemocytes in response to Vibrio alginolyticus infection

The pearl oyster Pinctada fucata is cultured widely for production of marine pearls in China, while mass mortalities, likely related to pathogenic infections, have occurred frequently in juvenile, mother and operated oysters. To address this issue, understanding host defense mechanisms of P. fucata against pathogenic challenge is extremely important. In the present study, a comparative analysis of hemocyte transcriptomes of P. fucata before and after Vibrio alginolyticus infection was conducted using the Illumina/Hiseq-2000 RNA-Seq technology. A total of 56,345,139 clean reads were generated and then assembled into 74,007 unigenes with an average length of 680 bp and an N50 of 1197 bp. Unigenes were annotated by comparing against non-redundant protein sequence (nr), non-redundant nucleotide (nt), Swiss-Prot, Pfam, Gene Ontology database (GO), Clusters of Orthologous Groups (COG) and Kyoto Encyclopedia of Genes and Genomes (KEGG) databases, and 29,615 unigenes (40.01%) were annotated in at least one database. There were 636 genes (518 up-regulated and 118 down-regulated) that were significantly differentially expressed after bacterial challenge, and among which 369 were associated with 122 pathways, including classical immune-related pathways, such as 'MAPK signaling pathway', 'Chemokine signaling pathway', 'Apoptosis' and 'Wnt signaling pathway'. These findings provide information on the pearl oyster innate immunity and may contribute to developing strategies for management of diseases and long-term sustainability of P. fucata culture.