Cancer-Associated Fibroblasts Influence Survival in Pleural Mesothelioma: Digital Gene Expression Analysis and Supervised Machine Learning Model

The exact mechanism of desmoplastic stromal reaction (DSR) formation is still unclear. The interaction between cancer cells and cancer-associated fibroblasts (CAFs) has an important role in tumor progression, while stromal changes are a poor prognostic factor in pleural mesothelioma (PM). We aimed to assess the impact of CAFs paracrine signaling within the tumor microenvironment and the DSR presence on survival, in a cohort of 77 PM patients. DSR formation was evaluated morphologically and by immunohistochemistry for Fibroblast activation protein alpha (FAP). Digital gene expression was analyzed using a custom-designed CodeSet (NanoString). Decision-tree-based analysis using the "conditional inference tree" (CIT) machine learning algorithm was performed on the obtained results. A significant association between FAP gene expression levels and the appearance of DSR was found (p = 0.025). DSR-high samples demonstrated a statistically significant prolonged median survival time. The elevated expression of MYT1, KDR, PIK3R1, PIK3R4, and SOS1 was associated with shortened OS, whereas the upregulation of VEGFC, FAP, and CDK4 was associated with prolonged OS. CIT revealed a three-tier system based on FAP, NF1, and RPTOR expressions. We could outline the prognostic value of CAFs-induced PI3K signaling pathway activation together with FAP-dependent CDK4 mediated cell cycle progression in PM, where prognostic and predictive biomarkers are urgently needed to introduce new therapeutic strategies.

Digital gene expression profiling and validation study highlight Cyclin F as an important regulator for sperm motility of chickens

In poultry industry, around 5 to 12% roosters were eliminated from the breeding program because of low sperm motility. Relatively few studies have been directed toward understanding and explaining the genetics mechanisms involved in sperm motility regulation in chickens. In the present study, digital gene expression (DGE) profiling and bioinformation analysis were used to explore the globally differentially expressed genes (DEG) in the testis of low sperm motility and high sperm motility roosters. Further validation study of key candidate genes was also performed. The DGE identified 652 DEGs, including 473 up-regulated and 179 down-regulated genes in the low sperm motility testis. Those DEGs were enriched on 21 terms of biological process category, 10 terms of cellular component category, including motile cilium, and 13 terms of molecular function category including microtubule motor activity and ATP binding. The kyoto encyclopedia of genes and genomes (KEGG) enrichment analysis indicated that these DEGs were involved in the FoxO signaling pathway and insulin resistance pathway. Quantitative real time PCR (qRT-PCR) studies of 8 DEGs were used to validate the DGE results. A key candidate gene Cyclin F (CCNF) was extremely low expressed in the low sperm motility testis (log2 ratio (low sperm motility/high sperm motility) = -5.23). The CCNF gene silencing in the chicken DF-1 cell line induced the reduced cell activity and proliferation. In summary, the present study provides insight into the potential genetic regulation of sperm motility and highlighted the underlying pathways (Insulin resistance and FoxO signaling pathways) and important candidate genes such as CCNF.

Digital gene expression analysis in mice lung with coinfection of influenza and streptococcus pneumoniae

Influenza A virus (IAV) and Streptococcus pneumoniae (SP) are two major upper respiratory tract pathogens that can also cause infection in polarized bronchial epithelial cells to exacerbate disease in coinfected individuals which may result in significant morbidity. However, the underlying molecular mechanism is poorly understood. Here, we employed BALB/c ByJ mice inflected with SP, IAV, IAV followed by SP (IAV+SP) and PBS (Control) as models to survey the global gene expression using digital gene expression (DGE) profiling. We attempt to gain insights into the underlying genetic basis of this synergy at the expression level. Gene expression profiles were obtain using the Illimina/Hisseq sequencing technique, and further analyzed by enrichment analysis of Gene Ontology (GO) and Pathway function. The hematoxylin-eosin (HE) staining revealed different tissue changes in groups during which IAV+SP group showed the most severe cell apoptosis. Compared with Control, a total of 2731, 3221 and 3946 differentially expressed genes (DEGs) were detected in SP, IAV and IAV+SP respectively. Besides, sixty-two GO terms were identified by Gene Ontology functional enrichment analysis, such as cell killing, biological regulation, response to stimulus, signaling, biological adhesion, enzyme regulator activity, receptor regulator activity and translation regulator activity. Pathway significant enrichment analysis indicated the dysregulation of multiple pathways, including apoptosis pathway. Among these, five selected genes were further verified by quantitative reverse transcription-polymerase chain reaction (qRT-PCR). This study shows that infection with SP, IAV or IAV+SP induces apoptosis with different degrees which might provide insights into the molecular mechanisms to facilitate further research.

Molecular characterisation of the broad-spectrum resistance to powdery mildew conferred by the Stpk-V gene from the wild species Haynaldia villosa

A key member of the Pm21 resistance gene locus, Stpk-V, derived from Haynaldia villosa, was shown to confer broad-spectrum resistance to wheat powdery mildew. The present study was planned to investigate the resistance mechanism mediated by Stpk-V. Transcriptome analysis was performed in Stpk-V transgenic plants and recipient Yangmai158 upon Bgt infection, and detailed histochemical observations were conducted. Chromosome location of Stpk-V orthologous genes in Triticeae species was conducted for evolutionary study and over-expression of Stpk-V both in barley and Arabidopsis was performed for functional study. The transcriptome results indicate, at the early infection stage, the ROS pathway, JA pathway and some PR proteins associated with the SA pathway were activated in both the resistant Stpk-V transgenic plants and susceptible Yangmai158. However, at the later infection stage, the genes up-regulated at the early stage were continuously held only in the transgenic plants, and a large number of new genes were also activated in the transgenic plants but not in Yangmai158. Results indicate that sustained activation of the early response genes combined with later-activated genes mediated by Stpk-V is critical for resistance in Stpk-V transgenic plants. Stpk-V orthologous genes in the representative grass species are all located on homologous group six chromosomes, indicating that Stpk-V is an ancient gene in the grasses. Over-expression of Stpk-V enhanced host resistance to powdery mildew in barley but not in Arabidopsis. Our results enable a better understanding of the resistance mechanism mediated by Stpk-V, and establish a solid foundation for its use in cereal breeding as a gene resource.

Profiling of differentially expressed genes in adipose tissues of multiple symmetric lipomatosis

Multiple symmetric lipomatosis (MSL) is a rare disorder characterized by aberrant multiple and symmetric subcutaneous adipose tissue accumulation in the face, neck, shoulders, back, chest and abdomen, severely affecting the quality of life of patients. At present, precise MSL etiology and pathogenesis remain to be elucidated. The present study first utilized a digital gene expression technique with a next‑generation sequencing platform to profile differentially expressed genes in three cases of MSL vs. normal control tissue. cDNA libraries from these tissue specimens were constructed and DNA sequenced for identification of differentially expressed genes, which underwent bioinformatic analysis using the Gene Ontology (GO) enrichment, Kyoto Encyclopedia of Genes and Genomes (KEGG) and protein‑protein interaction (PPI) network analyses. As a result, a total of 859 differentially expressed genes were identified, including 308 upregulated genes (C19orf80, Apelin, C21orf33, FAM166B and HSD11B2 were mostly upregulated 6.984‑, 4.670‑, 4.412‑, 3.693‑ and 3.561‑fold, respectively) and 551 downregulated genes [FosB proto‑oncogene, AP‑1 transcription factor subunit (FOSB), selectin (SEL) E, RAR related orphan receptor (ROR) B, salt inducible kinase (SIK)1 and epidermal growth factor‑like protein (EGFL)6 were mostly downregulated ‑9.845, ‑8.243, ‑8.123, ‑7.702 and ‑7.664 fold, respectively). The GO functional enrichment analysis demonstrated these differentially expressed genes were predominantly involved in biological processes and cellular components, while the KEGG pathway enrichment analysis demonstrated that ribosome, non‑alcoholic fatty liver disease, human T‑lymphotropic virus type 1 (HTLV‑I) infection and Alzheimer's disease pathways were altered in MSL. The PPI network data demonstrated ubiquitin C (UBC), translocator protein (TSPO), Jun Proto‑Oncogene, AP‑1 Transcription Factor (JUN) and FOS were among these differentially expressed genes that participated in regulation of adipocyte differentiation, although no previous study has linked them to MSL. In conclusion, the present study profiled differentially expressed genes in MSL and identified gene pathways that may be associated with MSL development and progression.

Phenotypic and Comparative Transcriptome Analysis of Different Ploidy Plants in Dendrocalamus latiflorus Munro

Elucidating the differences in gene expression profiles of plants with different ploidy levels and how they affect phenotypic traits is vital to allow genetic improvement of plants such as Ma bamboo (Dendrocalamus latiflorus Munro). We previously obtained triploid (2n = 3X = 36), hexaploid (2n = 6X = 72), and dodecaploid (2n = 12X = 144) Ma bamboo plants from embryogenic callus by anther culturing. Phenotypic differences between these plants appeared to be correlated with differences in ploidy. Here, we performed transcriptome profiling and sequencing of anther-regenerated plants and F1 seedlings of different ploidy levels using RNA-Seq technology. Pair-wise comparisons of the four resulting libraries revealed 8,396 differentially expressed genes. These differentially expressed genes were annotated, functionally classified, and partially validated. We found that the chromosome doubling led to substantially up- or down-regulation of genes that were involved in cell growth and differentiation; the polyploidy levels altered the anatomical, physiological and growth characteristics, such as leaf thickness, fusoid cell and stomatal size, shoot number, photosynthesis and respiration rate and so on. Additionally, two candidate genes, EXPB3 and TCP with potenitial regulatory roles in cell division and differentiation, were identified through gene coexpresseion network analysis. These results highlight the significance of potential applications of polyploidy, and provide valuable information for the genetic breeding of bamboo species.

Molecular and Cellular Mechanisms of Apoptosis during Dissociated Spermatogenesis

Apoptosis is a tightly controlled process by which tissues eliminate unwanted cells. Spontaneous germ cell apoptosis in testis has been broadly investigated in mammals that have an associated spermatogenesis pattern. However, the mechanism of germ cell apoptosis in seasonally breeding reptiles following a dissociated spermatogenesis has remained enigmatic. In the present study, morphological evidence has clearly confirmed the dissociated spermatogenesis pattern in Pelodiscus sinensis. TUNEL and TEM analyses presented dynamic changes and ultrastructural characteristics of apoptotic germ cells during seasonal spermatogenesis, implying that apoptosis might be one of the key mechanisms to clear degraded germ cells. Furthermore, using RNA-Seq and digital gene expression (DGE) profiling, a large number of apoptosis-related differentially expressed genes (DEGs) at different phases of spermatogenesis were identified and characterized in the testis. DGE and RT-qPCR analysis revealed that the critical anti-apoptosis genes, such as Bcl-2, BAG1, and BAG5, showed up-regulated patterns during intermediate and late spermatogenesis. Moreover, the increases in mitochondrial transmembrane potential in July and October were detected by JC-1 staining. Notably, the low protein levels of pro-apoptotic cleaved caspase-3 and CytC in cytoplasm were detected by immunohistochemistry and western blot analyses, indicating that the CytC-Caspase model might be responsible for the effects of germ cell apoptosis on seasonal spermatogenesis. These results facilitate understanding the regulatory mechanisms of apoptosis during spermatogenesis and uncovering the biological process of the dissociated spermatogenesis system in reptiles.

Knockdown of two trehalose-6-phosphate synthases severely affects chitin metabolism gene expression in the brown planthopper Nilaparvata lugens

BACKGROUND

RNA interference combined with digital gene expression (DGE) analysis can be used to study gene function. Trehalose-6-phosphate synthase (TPS) plays a key role in the synthesis of trehalose and insect development.

RESULTS

DGE analysis revealed that the expression of nine or four chitinase genes was reduced significantly 48 h after NlTPS1 and NlTPS2 knockdown by RNAi, respectively. Additionally, abnormal phenotypes were noted, and approximately 30% of insects died. HK and G6PI2 expression decreased significantly whereas GFAT, GNPNA and UAP expression increased significantly 72 h after NlTPS1 and NlTPS2 knockdown. PGM1 expression decreased significantly after TPS2 knockdown, whereas PGM2 expression increased significantly and the expression of three CHS genes decreased 48 h after TPS1 knockdown. The mRNA expression of all 12 chitin degradation genes decreased 48 h after NlTPS1 and NlTPS2 treatment, and Cht2, Cht3, Cht6, Cht7, Cht10 and ENGase levels remained significantly decreased up to 72 h after NlTPS1 and NlTPS2 knockdown.

CONCLUSIONS

These results demonstrate that silencing of TPS genes can lead to increased moulting deformities and mortality rates owing to the misregulation of genes involved in chitin metabolism, and TPS genes are potential pest control targets in the future. © 2016 Society of Chemical Industry.

Digital gene expression profiling analysis of DNA repair pathways in colon cancer stem population of HT29 cells

Cancer stem cells (CSCs) contribute to the relapse and development of new neoplasm lesions. While most available clinical approaches, such as chemical and radiation therapies, will kill the majority of cancer cells, they do not kill them all. Some resisting cells, like CSCs, are able to survive due to their excellent self-maintaining capabilities, even in challenging environments. In the present study, we investigated the mRNA level of DNA repair genes of colon CSCs from the HT29 cell line in response to single-strand damage and double-strand breaks, as well as the evident upregulation of key genes in base excision repair, mismatch repair, non-homologous end-joining, and homologous recombination pathways in these cells. Digital gene expression analysis identified upregulated genes in CD44⁺ HT29 cells that may play important roles in DNA repair. Our results reveal that colon CSCs bear efficient DNA repair abilities, which might explain the survival of colon CSCs after repeated chemical and radiation therapy.

Effect of oxidative phosphorylation signaling pathway on silkworm midgut following exposure to phoxim

Organophosphate pesticides are applied widely in the world for agricultural purposes, and their exposures often resulted in non-cocooning of Bombyx mori in China. Silkworm midgut is the major organ for digestion and nutrient absorption, importantly it is also a barrier against foreign substances and chemical pesticides. The purpose of this study was to determine the mechanism of oxidative injury in silkworm midgut with phoxim induction. The results showed that the transcription level of oxidative phosphorylation signaling pathway genes of midgut under phoxim stress. Digital gene expression (DGE) analysis revealed that 24 electron transport chain (ETC)-related genes were upregulated. Quantitative real time polymerase chain reaction results indicated that the ETC the genes encoding NADH-CoQ1, Succinic-Q, cyt c reductase-S, cyt c oxidase-S, cytochrome c oxidase polypeptide IV, ATP synthase, and vacuolar H+ ATP synthase were all significantly up-regulated by 1.50-, 1.31-, 1.42-, 1.44-, 1.70-, 2.03- and 1.43-fold, respectively. Phoxim induction enhanced the activity of ETC complex in mitochondria, and induced the accumulation of ROS in midgut. These results indicated that trace phoxim enhanced respiration in midgut, and the imbalance between the activity changes of ETC may led to reactive oxygen species accumulation. The ETC of mitochondria may be potential biomarkers of midgut toxicity in B. mori caused by phoxim exposure. © 2015 Wiley Periodicals, Inc. Environ Toxicol 32: 167-175, 2017.

Transcriptome Sequencing and De Novo Assembly of Golden Cuttlefish Sepia esculenta Hoyle

Golden cuttlefish Sepia esculenta Hoyle is an economically important cephalopod species. However, artificial hatching is currently challenged by low survival rate of larvae due to abnormal embryonic development. Dissecting the genetic foundation and regulatory mechanisms in embryonic development requires genomic background knowledge. Therefore, we carried out a transcriptome sequencing on Sepia embryos and larvae via mRNA-Seq. 32,597,241 raw reads were filtered and assembled into 98,615 unigenes (N50 length at 911 bp) which were annotated in NR database, GO and KEGG databases respectively. Digital gene expression analysis was carried out on cleavage stage embryos, healthy larvae and malformed larvae. Unigenes functioning in cell proliferation exhibited higher transcriptional levels at cleavage stage while those related to animal disease and organ development showed increased transcription in malformed larvae. Homologs of key genes in regulatory pathways related to early development of animals were identified in Sepia. Most of them exhibit higher transcriptional levels in cleavage stage than larvae, suggesting their potential roles in embryonic development of Sepia. The de novo assembly of Sepia transcriptome is fundamental genetic background for further exploration in Sepia research. Our demonstration on the transcriptional variations of genes in three developmental stages will provide new perspectives in understanding the molecular mechanisms in early embryonic development of cuttlefish.

Digital Gene Expression Profiling to Explore Differentially Expressed Genes Associated with Terpenoid Biosynthesis during Fruit Development in Litsea cubeba

Mountain pepper (Litseacubeba (Lour.) Pers.) (Lauraceae) is an important industrial crop as an ingredient in cosmetics, pesticides, food additives and potential biofuels. These properties are attributed to monoterpenes and sesquiterpenes. However, there is still no integrated model describing differentially expressed genes (DEGs) involved in terpenoid biosynthesis during the fruit development of L. cubeba. Here, we performed digital gene expression (DGE) using the Illumina NGS platform to evaluated changes in gene expression during fruit development in L. cubeba. DGE generated expression data for approximately 19354 genes. Fruit at 60 days after flowering (DAF) served as the control, and a total of 415, 1255, 449 and 811 up-regulated genes and 505, 1351, 1823 and 1850 down-regulated genes were identified at 75, 90, 105 and 135 DAF, respectively. Pathway analysis revealed 26 genes involved in terpenoid biosynthesis pathways. Three DEGs had continued increasing or declining trends during the fruit development. The quantitative real-time PCR (qRT-PCR) results of five differentially expressed genes were consistent with those obtained from Illumina sequencing. These results provide a comprehensive molecular biology background for research on fruit development, and information that should aid in metabolic engineering to increase the yields of L. cubeba essential oil.

Transcriptome Profiling of Taproot Reveals Complex Regulatory Networks during Taproot Thickening in Radish (Raphanus sativus L.)

Radish (Raphanus sativus L.) is one of the most important vegetable crops worldwide. Taproot thickening represents a critical developmental period that determines yield and quality in radish life cycle. To isolate differentially expressed genes (DGEs) involved in radish taproot thickening process and explore the molecular mechanism underlying taproot development, three cDNA libraries from radish taproot collected at pre-cortex splitting stage (L1), cortex splitting stage (L2), and expanding stage (L3) were constructed and sequenced by RNA-Seq technology. More than seven million clean reads were obtained from the three libraries, from which 4,717,617 (L1, 65.35%), 4,809,588 (L2, 68.24%) and 4,973,745 (L3, 69.45%) reads were matched to the radish reference genes, respectively. A total of 85,939 transcripts were generated from three libraries, from which 10,450, 12,325, and 7392 differentially expressed transcripts (DETs) were detected in L1 vs. L2, L1 vs. L3, and L2 vs. L3 comparisons, respectively. Gene Ontology and pathway analysis showed that many DEGs, including EXPA9, Cyclin, CaM, Syntaxin, MADS-box, SAUR, and CalS were involved in cell events, cell wall modification, regulation of plant hormone levels, signal transduction and metabolisms, which may relate to taproot thickening. Furthermore, the integrated analysis of mRNA-miRNA revealed that 43 miRNAs and 92 genes formed 114 miRNA-target mRNA pairs were co-expressed, and three miRNA-target regulatory networks of taproot were constructed from different libraries. Finally, the expression patterns of 16 selected genes were confirmed using RT-qPCR analysis. A hypothetical model of genetic regulatory network associated with taproot thickening in radish was put forward. The taproot formation of radish is mainly attributed to cell differentiation, division and expansion, which are regulated and promoted by certain specific signal transduction pathways and metabolism processes. These results could provide new insights into the complex molecular mechanism underlying taproot thickening and facilitate genetic improvement of taproot in radish.

Transcriptomic Analysis of Differentially Expressed Genes During Larval Development of Rapana venosa by Digital Gene Expression Profiling

During the life cycle of shellfish, larval development, especially metamorphosis, has a vital influence on the dynamics, distribution, and recruitment of natural populations, as well as seed breeding. Rapana venosa, a carnivorous gastropod, is an important commercial shellfish in China, and is an ecological invader in the United States, Argentina, and France. However, information about the mechanism of its early development is still limited, because research in this area has long suffered from a lack of genomic resources. In this study, 15 digital gene expression (DGE) libraries from five developmental stages of R. venosa were constructed and sequenced on the IIIumina Hi-Sequation 2500 platform. Bioinformaticsanalysis identified numerous differentially and specifically expressed genes, which revealed that genes associated with growth, nervous system, digestive system, immune system, and apoptosis participate in important developmental processes. The functional analysis of differentially expressed genes was further implemented by gene ontology, and Kyoto encyclopedia of genes and genomes enrichment. DGE profiling provided a general picture of the transcriptomic activities during the early development of R. venosa, which may provide interesting hints for further study. Our data represent the first comparative transcriptomic information available for the early development of R. venosa, which is a prerequisite for a better understanding of the physiological traits controlling development.

Comprehensive and Quantitative Proteomic Analysis of Metamorphosis-Related Proteins in the Veined Rapa Whelk, Rapana venosa

Larval metamorphosis of the veined rapa whelk (Rapana venosa) is a pelagic to benthic transition that involves considerable structural and physiological changes. Because metamorphosis plays a pivotal role in R. venosa commercial breeding and natural populations, the endogenous proteins that drive this transition attract considerable interest. This study is the first to perform a comprehensive and quantitative proteomic analysis related to metamorphosis in a marine gastropod. We analyzed the proteomes of competent R. venosa larvae and post-larvae, resulting in the identification of 5312 proteins, including 470 that were downregulated and 668 that were upregulated after metamorphosis. The differentially expressed proteins reflected multiple processes involved in metamorphosis, including cytoskeleton and cell adhesion, ingestion and digestion, stress response and immunity, as well as specific tissue development. Our data improve understanding of the physiological traits controlling R. venosa metamorphosis and provide a solid basis for further study.

Insight into the possible mechanism of the summer diapause of Delia antiqua (Diptera: Anthomyiidae) through digital gene expression analysis

The onion fly, Delia antiqua, is a major underground agricultural pest that can enter pupal diapause in the summer and winter seasons. However, little is known about its molecular regulation due to the lack of genomic resources. To gain insight into the possible mechanism of summer diapause (SD), high-throughput RNA-Seq data were generated from non-diapause (ND) and SD (initial, maintenance and quiescence phase) pupae. Three pair-wise comparisons were performed and identified, 1380, 1471 and 435, and were significantly regulated transcripts. Further analysis revealed that the enrichment of several functional terms related to juvenile hormone regulation, cell cycle, carbon hydrate and lipid metabolism, innate immune and stress responses, various signalling transductions, ubiquitin-dependent proteosome, and variation in cuticular and cytoskeleton components were found between ND and SD and between different phases of SD. Global characterization of transcriptome profiling between SD and ND contributes to the in-depth elucidation of the molecular mechanism of SD. Our results also offer insights into the evolution of insect diapause and support the importance of using the onion fly as a model to compare the molecular regulation events of summer and winter diapauses.

GeneGazer: A Toolkit Integrating Two Pipelines for Personalized Profiling and Biosignature Identification

BACKGROUND

Next-generation sequencing provides useful information about gene mutations, gene expression, epigenetic modification, microRNA expression, and copy number variations. More and more computing tools have been developed to analyze this large quantity of information. However, to test and find suitable analytical tools and integrate their results is tedious and challenging for users with little bioinformatics training. In the present study, we assembled the computing tools into a convenient toolkit to simplify the analysis and integration of data between bioinformatics tools.

MATERIALS AND METHODS

The toolkit, GeneGazer, comprises of two parts: the first, named Gaze_Profiler, was designed for personalized molecular profiling from next-generation sequencing data of paired samples; the other, named Gaze_BioSigner, was designed for the discovery of disease-associated biosignatures from expressional and mutational profiles of a cohort study.

RESULTS

To demonstrate the capabilities of Gaze_Profiler, we analyzed a pair (colon cancer and adjacent normal tissues) of RNA-sequencing data from one patient downloaded from the Sequencing Read Archive database and used them to profile somatic mutations and digital gene expression. In this case, alterations in the RAS/RAF/MEK/ERK signaling pathway (activated by KRAS G13D mutation) and canonical WNT signaling pathway (activated by truncated APC) were identified; no EGFR mutation or overexpression was found. These data suggested a limited efficacy of cetuximab in the patient. To demonstrate the ability of Gazer_BioSigner, we analyzed gene-expression data from 192 cancer tissues downloaded from The Cancer Genome Atlas and found that the activation of cAMP/PKA signaling, OCT-3/4 and SRF were associated with colon cancer progression and could be potential therapeutic targets.

CONCLUSION

GeneGazer is a reliable and robust toolkit for the analysis of data from high-throughput platforms and has potential for clinical application and biomedical research.

Transcriptome Analysis and Identification of Genes Associated with Floral Transition and Flower Development in Sugar Apple (Annona squamosa L.)

Sugar apple (Annona squamosa L.) is a semi-deciduous subtropical tree that progressively sheds its leaves in the spring. However, little information is available on the mechanism involved in flower developmental pattern. To gain a global perspective on the floral transition and flower development of sugar apple, cDNA libraries were prepared independently from inflorescent meristem and three flowering stages. Illumina sequencing generated 107,197,488 high quality reads that were assembled into 71,948 unigenes, with an average sequence length of 825.40 bp. Among the unigenes, various transcription factor families involved in floral transition and flower development were elucidated. Furthermore, a Kyoto Encyclopedia of Genes and Genomes pathway enrichment analysis revealed that unigenes exhibiting differential expressions were involved in various phytohormone signal transduction events and circadian rhythms. In addition, 147 unigenes exhibiting sequence similarities to known flowering-related genes from other plants were differentially expressed during flower development. The expression patterns of 20 selected genes were validated using quantitative-PCR. The expression data presented in our study is the most comprehensive dataset available for sugar apple so far and will serve as a resource for investigating the genetics of the flowering process in sugar apple and other Annona species.

Transcriptomic comparison of the self-pollinated and cross-pollinated flowers of Erigeron breviscapus to analyze candidate self-incompatibility-associated genes

BACKGROUND

Self-incompatibility (SI) is a widespread and important mating system that promotes outcrossing in plants. Erigeron breviscapus, a medicinal herb used widely in traditional Chinese medicine, is a self-incompatible species of Asteraceae. However, the genetic characteristics of SI responses in E. breviscapus remain largely unknown. To understand the possible mechanisms of E. breviscapus in response to SI, we performed a comparative transcriptomic analysis with capitulum of E. breviscapus after self- and cross-pollination, which may provide valuable information for analyzing the candidate SI-associated genes of E. breviscapus.

METHODS

Using a high-throughput next-generation sequencing (Illumina) approach, the transcriptionexpression profiling of the different genes of E. breviscapus were obtained, some results were verified by quantitative real time PCR (qRT-PCR).

RESULTS

After assembly, 63,485 gene models were obtained (average gene size 882 bp; N50 = 1485 bp), among which 38,540 unigenes (60.70% of total gene models) were annotated by comparisons with four public databases (Nr, Swiss-Prot, KEGG and COG): 38,338 unigenes (60.38% of total gene models) showed high homology with sequences in the Nr database. Differentially expressed genes were identified among the three cDNA libraries (non-, self- and cross-pollinated capitulum of E. breviscapus), and approximately 230 genes might be associated with SI responses. Several these genes were upregulated in self-pollinated capitulum but downregulated in cross-pollinated capitulum, such as SRLK (SRK-like) and its downstream signal factor, MLPK. qRT-PCR confirmed that the expression patterns of EbSRLK1 and EbSRLK3 genes were not closely related to SI of E. breviscapus.

CONCLUSIONS

This work represents the first large-scale analysis of gene expression in the self-pollinated and cross-pollinated flowers of E. breviscapus. A larger number of notable genes potentially involved in SI responses showed differential expression, including genes playing crucial roles in cell-cell communication, signal transduction and the pollination process. We thus hypothesized that those genes showing differential expression and encoding critical regulators of SI responses, such as MLPK, ARC1, CaM, Exo70A1, MAP, SF21 and Nod, might affect SI responses in E. breviscapus. Taken together, our study provides a pool of SI-related genes in E. breviscapus and offers a valuable resource for elucidating the mechanisms of SI in Asteraceae.

Identification of Ramie Genes in Response to Pratylenchus coffeae Infection Challenge by Digital Gene Expression Analysis

Root lesion disease, caused by Pratylenchus coffeae, seriously impairs the growth and yield of ramie, an important natural fiber crop. The ramie defense mechanism against P. coffeae infection is poorly understood, which hinders efforts to improve resistance via breeding programs. In this study, the transcriptome of the resistant ramie cultivar Qingdaye was characterized using Illumina sequence technology. About 46.3 million clean pair end (PE) reads were generated and assembled into 40,826 unigenes with a mean length of 830 bp. Digital gene expression (DGE) analysis was performed on both the control roots (CK) and P. coffeae-challenged roots (CH), and the differentially expressed genes (DEGs) were identified. Approximately 10.16 and 8.07 million cDNA reads in the CK and CH cDNA libraries were sequenced, respectively. A total of 137 genes exhibited different transcript abundances between the two libraries. Among them, the expressions of 117 and 20 DEGs were up- and down-regulated in P. coffeae-challenged ramie, respectively. The expression patterns of 15 candidate genes determined by qRT-PCR confirmed the results of DGE analysis. Time-course expression profiles of eight defense-related genes in susceptible and resistant ramie cultivars were different after P. coffeae inoculation. The differential expression of protease inhibitors, pathogenesis-related proteins (PRs), and transcription factors in resistant and susceptible ramie during P. coffeae infection indicated that cystatin likely plays an important role in nematode resistance.

c-Fos is involved in inhibition of human bladder carcinoma T24 cells by brazilin

Crude brazilin extract from Sappan wood has demonstrated strong anti tumor activity in the mouse model of human bladder carcinoma and clinical trial for intravesical therapy. Purified brazilin was confirmed the most active molecule in inhibition of bladder carcinoma T24 cells. Brazilin decreased proliferation and viability of T24 cells in a dose- and time-dependent manner, with a calculated LC50 of 32 µg/mL. More than 1,000 of genes were found upregulated and down regulated by brazilin treatment in digital gene expression profiling. Gene ontology analysis indicated that stress response, apoptosis, and cell cycle regulatory pathways were highly enriched. Among the regulated genes, c-Fos was the most and specifically upregulated. Overexpression of c-Fos in T24 cells resulted in tumor cell specific changes in cell morphology and viability. Over expression of stress-responsive gene, HSP70, and other highly upregulated genes did not have any effect on cell growth. Brazilin may inhibit T24 cell growth and trigger cell death through a c-Fos-mediated and tumor cell specific signaling pathway. Further studies of its down stream mediators may help to identify better tumor cell type specific drug targets.

Gene expression profiling in winged and wingless cotton aphids, Aphis gossypii (Hemiptera: Aphididae)

While trade-offs between flight capability and reproduction is a common phenomenon in wing dimorphic insects, the molecular basis is largely unknown. In this study, we examined the transcriptomic differences between winged and wingless morphs of cotton aphids, Aphis gossypii, using a tag-based digital gene expression (DGE) approach. Ultra high-throughput Illumina sequencing generated 5.30 and 5.39 million raw tags, respectively, from winged and wingless A. gossypii DGE libraries. We identified 1,663 differentially expressed transcripts, among which 58 were highly expressed in the winged A. gossypii, whereas 1,605 expressed significantly higher in the wingless morphs. Bioinformatics tools, including Gene Ontology, Cluster of Orthologous Groups, euKaryotic Orthologous Groups and Kyoto Encyclopedia of Genes and Genomes pathways, were used to functionally annotate these transcripts. In addition, 20 differentially expressed transcripts detected by DGE were validated by the quantitative real-time PCR. Comparative transcriptomic analysis of sedentary (wingless) and migratory (winged) A. gossyii not only advances our understanding of the trade-offs in wing dimorphic insects, but also provides a candidate molecular target for the genetic control of this agricultural insect pest.

Why mosaic? Gene expression profiling of African cassava mosaic virus-infected cassava reveals the effect of chlorophyll degradation on symptom development

Cassava mosaic disease, caused by cassava begomoviruses, is the most serious disease for cassava in Africa. However, the pathogenesis of this disease is poorly understood. We employed high throughput digital gene expression profiling based on the Illumina Solexa sequencing technology to investigate the global transcriptional response of cassava to African cassava mosaic virus infection. We found that 3,210 genes were differentially expressed in virus-infected cassava leaves. Gene ontology term and Kyoto Encyclopedia of Genes and Genomes pathway analysis indicated that genes implicated in photosynthesis were most affected, consistent with the chlorotic symptoms observed in infected leaves. The upregulation of chlorophyll degradation genes, including the genes encoding chlorophyllase, pheophytinase, and pheophorbide a oxygenase, and downregulation of genes encoding the major apoproteins in light-harvesting complex II were confirmed by qRT-PCR. These findings, together with the reduction of chlorophyll b content and fewer grana stacks in the infected leaf cells, reveal that the degradation of chlorophyll plays an important role in African cassava mosaic virus symptom development. This study will provide a road map for future investigations into viral pathogenesis.

The adverse effects of phoxim exposure in the midgut of silkworm, Bombyx mori

The silkworm is an important economic insect. Poisoning of silkworms by organophosphate pesticides causes tremendous loss to the sericulture. In this study, Solexa sequencing technology was performed to profile the gene expression changes in the midgut of silkworms in response to 24h of phoxim exposure and the impact on detoxification, apoptosis and immune defense were addressed. The results showed that 254 genes displayed at least 2.0-fold changes in expression levels, with 148 genes up-regulated and 106 genes down-regulated. Cytochrome P450 played an important role in detoxification. Histopathology examination and transmission electron microscope revealed swollen mitochondria and disappearance of the cristae of mitochondria, which are the important features in insect apoptotic cells. Cytochrome C release from mitochondria into the cytoplasm was confirmed. In addition, the Toll and immune deficiency (IMD) signal pathways were all inhibited using qRT-PCR. Our results could help better understand the impact of phoxim exposure on silkworm.

De novo characterization of transcriptome and gene expression dynamics in epidermis during the larval-pupal metamorphosis of common cutworm

Larval cuticle is degraded and replaced by the pupal counterpart during larval-pupal metamorphosis in the holometabolous insects. In addition to the extrinsic transformation, the epidermis goes through significant changes at molecular levels. To elucidate the intrinsic mechanism of epidermal metamorphosis, the dynamics of chitin content in the cuticle was examined in an important agricultural lepidopteran, the common cutworm, and the transcriptome was analyzed using Illumina sequencing technology. Gene expression profiles during the metamorphosis were further studied by both the digital gene expression (DGE) system and real-time quantitative PCR. The results showed that the chitin content decreased in prepupae and then increased in pupae. A total of 58 million sequencing reads were obtained and assembled into 70,346 unigenes. Over 9000 unigenes were identified to express differentially during the transformation process. As compared with the 6th instar feeding larvae, the most significant changes took place in the proteasome and metabolic pathways in prepupae and pupae, respectively. The cytochrome P450s, VHDLs, chitinase, serine protease and genes involved in sex pheromone biosynthesis changed their mRNA levels remarkably. Three chitinolytic enzymes (chitinase, β-N-acetylglucosaminidase and chitin deacetylase) showed distinct mRNA expression patterns, the former two enzymes revealed the highest expression in prepupae, however the latter one showed its climax mRNA level in pupae. The gene expression patterns suggest that chitinase and β-N-acetylglucosaminidase may be responsible for the degradation of larval cuticles, whereas chitin deacetylase may help to degrade the pupal counterparts. Gene expression dynamics also implied that the chitin of pupal cuticle might be formed by recycling of the degraded chitin of larval cuticle rather than through de novo synthesis. The 20E-induced nuclear receptors seem to be important factors regulating chitin metabolic enzymes during the cuticle remodeling. Our data provide a comprehensive resource for exploring the molecular mechanism of epidermal metamorphosis in insects.

Transcriptome analysis reveals novel genes potentially involved in photoperiodic tuberization in potato

Potato microtuber produced in vitro provides a model system to investigate photoperiod-dependent tuberization. However, the genes associated with potato tuberization remain to be elucidated. The present research involved three potato clones with distinct tuberization response to changes of photoperiod. Digital Gene Expression (DGE) Tag Profiling analysis of the short-day-sensitive clone identified 2218 genes that were regulated by day length. Both GO and KEGG pathway analysis provided insights into predominant biological processes and pathways, and enabled the selection of 56 genes associated with circadian rhythmicity, signal transduction, and development. Quantitative transcriptional analysis in the selected clones revealed 5 genes potentially associated with photoperiodic tuberization, which were predicted to encode a DOF protein, a blue light receptor, a lectin, a syntaxin-like protein, and a protein with unknown function. Our results strongly suggest that potato tuberization may be largely controlled by the homologs of genes shown to regulate flowering time in other plants.

Identification and expression analysis of a new glycoside hydrolase family 55 exo-β-1,3-glucanase-encoding gene in Volvariella volvacea suggests a role in fruiting body development

The edible straw mushroom Volvariella volvacea is an important crop in South East Asia and is predominantly harvested in the egg stage. Rapid stipe elongation and cap expansion result in a swift transition from the egg to elongation and maturation stage, which are subjected to fast senescence and deterioration. In other mushrooms, β-1,3-glucanases have been associated with degradation (softening) of the cell wall during stipe elongation and senescence. We present a new glycoside hydrolase family 55 (GH55) exo-β-1,3-glucanase gene, exg2, and highly conserved deduced EXG2 protein. The 3D model and presumed catalytic residues of V. volvacea EXG2 are identical to Lentinula edodes EXG2 and Phanerochaete chrysosporium Lam55A, supporting similar enzymatic functions. In addition to previous association to stipe elongation and senescence, our data clearly indicates a role for cap (pileus) expansion. Digital gene expression, quantitative PCR and isobaric tags for relative and absolute quantification analysis showed low exg2 and EXG2 levels in primordia, button, egg and elongation stages and significantly increased levels in the maturation stage. Subsequent relative quantitative PCR analysis designated expression of exg2 to the stipe in the elongation stage and to the pileus and stipe in the maturation stage. EXG2 cell wall softening activity, close correlation of exg2 expression with the principal expanding mushroom tissues and a strong conservation of expression patterns and protein sequences in other mushrooms, make V. volvacea exg2 an important candidate for future studies on mechanisms of fruiting body expansion and senescence causing commodity value loss.