iAssembler: a package for de novo assembly of Roche-454/Sanger transcriptome sequences

Evolutionary analyses and expression patterns of TCP genes in Ranunculales

TCP transcription factors play a role in a large number of developmental processes and are at the crossroads of numerous hormonal biosynthetic and signaling pathways. The complete repertoire of TCP genes has already been characterized in several plant species, but not in any species of early diverging eudicots. We focused on the order Ranunculales because of its phylogenetic position as sister group to all other eudicots and its important morphological diversity. Results show that all the TCP genes expressed in the floral transcriptome of Nigella damascena (Ranunculaceae) are the orthologs of the TCP genes previously identified from the fully sequenced genome of Aquilegia coerulea. Phylogenetic analyses combined with the identification of conserved amino acid motifs suggest that six paralogous genes of class I TCP transcription factors were present in the common ancestor of angiosperms. We highlight independent duplications in core eudicots and Ranunculales within the class I and class II subfamilies, resulting in different numbers of paralogs within the main subclasses of TCP genes. This has most probably major consequences on the functional diversification of these genes in different plant clades. The expression patterns of TCP genes in Nigella damascena were consistent with the general suggestion that CIN and class I TCP genes may have redundant roles or take part in same pathways, while CYC/TB1 genes have more specific actions. Our findings open the way for future studies at the tissue level, and for investigating redundancy and subfunctionalisation in TCP genes and their role in the evolution of morphological novelties.

In-Depth Analysis of an Obligate Anaerobe Paraclostridium bifermentans Isolated from Uterus of Bubalus bubalis

Simple Summary

Non-specific uterine infections in bovine (uterine line inflammation) are a significant issue for the dairy industry. These infections are responsible for significant financial losses all over the world. Paraclostridium bifermentans is an obligate anaerobe, gram-positive rod-shaped bacteria belonging to the Clostridia class and the Peptostreptococcaceae family produces endospores. This bacterium has the ability to infiltrate bovine uterine endometrial epithelial cells and cause infection in the endometrium epithelial cells. Our study found that an examination of a buffalo uterus with yellowish purulent discharge reported the presence of pathogenic bacteria Paraclostridium bifermentans, where its genomic characterization, substrate utilization, and antibiotic susceptibility potentiality was studied. This discovery indicates the presence of virulence genes as well as pathogenic features. This is the first report of P. bifermentans from the bovine uterus environment.

Abstract

Chronic non-specific contamination of the reproductive tract in animals is a major issue during early postpartum, natural coitus, or artificial insemination. Uterine infection is one of the major concerns reducing fertility, production loss, and early culling of the animals. Therefore, the aim of this study was to identify any novel bacterium if present in the uterine environment of Bubalus bubalis causing infections. A strictly anaerobic bacterial strain designated as Paraclostridium bifermentans GBRC was isolated and characterized. Bacterium was found to be Gram positive moderate rod with motility. The optimum growth was observed at 40 ± 2 °C. The pathogenic characteristics of the GBRC strain, such as hemolysis, gelatin hydrolysis, and the production of volatile sulfur compounds, were similar to those seen in the epithelial layer invading pathogenic strains. Assembled genome size was 3.6 MB, with 78 contigs, and a G + C content of 28.10%. Furthermore, the whole genome sequence analysis confirmed the presence of genes encoding virulence factors and provided genomic insights on adaptation of the strain in the uterine environment. Based on the phenotypic and genetic differences with phylogenetic relatives, strain GBRC is proposed to represent a first reported species of the genus Paraclostridium with potential pathogenic character, from the buffalo uterine environment. This study analysis of the GBRC strain serves as a key reference point for the investigation of potential pathogenic strains that may cause endometritis and metritis in bovine.

Transcriptome Profiling Reveals a Petunia Transcription Factor, PhCOL4, Contributing to Antiviral RNA Silencing

RNA silencing is a common antiviral mechanism in eukaryotic organisms. However, the transcriptional regulatory mechanism that controls the RNA silencing process remains elusive. Here, we performed high-depth transcriptome analysis on petunia (Petunia hybrida) leaves infected with tobacco rattle virus (TRV) strain PPK20. A total of 7,402 differentially expressed genes (DEGs) were identified. Of them, some RNA silencing-related transcripts, such as RNA-dependent RNA polymerases (RDRs), Dicer-like RNase III enzymes (DCLs), and Argonautes (AGOs), were induced by viral attack. Furthermore, we performed TRV-based virus-induced gene silencing (VIGS) assay on 39 DEGs encoding putative transcription factors (TFs), using green fluorescent protein (GFP) and phytoene desaturase (PhPDS) as reporters. Results showed that the down-regulation of PhbHLH41, PhbHLH93, PhZPT4-3, PhCOL4, PhHSF-B3A, PhNAC90, and PhWRKY75 led to enhanced TRV accumulation and inhibited PhPDS-silenced photobleaching phenotype. In contrast, silencing of PhERF22 repressed virus accumulation and promoted photobleaching development. Thus, these TFs were identified as potential positive and negative regulators of antiviral RNA silencing, respectively. One positive regulator PhCOL4, belonging to the B-box zinc finger family, was selected for further functional characterization. Silencing and transient overexpression of PhCOL4 resulted in decreased and increased expression of several RNA silencing-related genes. DNA affinity purification sequencing analysis revealed that PhCOL4 targeted PhRDR6 and PhAGO4. Dual luciferase and yeast one-hybrid assays determined the binding of PhCOL4 to the PhRDR6 and PhAGO4 promoters. Our findings suggest that TRV-GFP-PhPDS-based VIGS could be helpful to identify transcriptional regulators of antiviral RNA silencing.

Cellular and molecular characterizations of the irregular internode division zone formation of a slow-growing bamboo variant

The key molecular mechanisms underlying the sectionalized growth within bamboo or other grass internodes remain largely unknown. Here, we genetically and morphologically compared the culm and rhizome internode division zones (DZs) of a slow-growing bamboo variant (sgv) having dwarf internodes, with those of the corresponding wild type (WT). Histological analysis discovers that the sgv has an irregular internode DZ. However, the shoot apical meristems in height, width, outside shape, cell number and cell width of the sgv and the WT were all similar. The DZ irregularities first appeared post apical meristem development, in 1-mm sgv rhizome internodes. Thus, the sgv is a DZ irregularity bamboo variant, which has been first reported in bamboo according to our investigation. Transcriptome sequencing analysis finds that a number of cell wall biogenesis and cell division-related genes are dramatically downregulated in the sgv DZ. Interestingly, both transcriptomic and brassinosteroid (BR) contents detection, as well as quantitative real-time PCR analyses show that these irregularities have resulted from the BR signaling pathway defects. Brassinosteroid defect might also cause the erect leaves and branches as well as the irregular epidermis of the sgv. These results suggest that BR signaling pathway plays critical roles in bamboo internode DZ and leaf development from a mutant perspective and also explain the upstream mechanisms causing the dwarf internode of the sgv bamboo.

C-CorA: A Cluster-Based Method for Correlation Analysis of RNA-Seq Data

Correlation analysis is a routine method of biological data analysis. In the process of RNA-Seq analysis, differentially expressed genes could be identified by calculating the correlation coefficients in the comparison of gene expression vs. phenotype or gene expression vs. gene expression. However, due to the complicated genetic backgrounds of perennial fruit, the correlation coefficients between phenotypes and genes are usually not high in fruit quality studies. In this study, a cluster-based correlation analysis method (C-CorA) is presented for fruit RNA-Seq analysis. C-CorA is composed of two main parts: the clustering analysis and the correlation analysis. The algorithm is described and then integrated into the MATLAB code and the C# WPF project. The C-CorA method was applied to RNA-Seq datasets of loquat (Eriobotrya japonica) fruit stored or ripened under different conditions. Low temperature conditioning or heat treatment of loquat fruit can alleviate the extent of lignification that occurs because of postharvest storage under low temperatures (0 °C). The C-CorA method generated correlation coefficients and identified many candidate genes correlated with lignification, including EjCAD3 and EjCAD4 and transcription factors such as MYB (UN00328). C-CorA is an effective new method for the correlation analysis of various types of data with different dimensions and can be applied to RNA-Seq data for candidate gene detection in fruit quality studies.

A robust, cost-effective and widely applicable whole-genome sequencing protocol for capripoxviruses

•

Robust method for the genomic characterization of all Capripoxviruses.

•

Pre-sequencing enrichment method based on targeted long-range PCR amplification.

•

Method applicable to low titre samples such as blood samples and vaccine batches.

•

Viral DNA enrichment compatible with various sample types and sequencing platforms.

•

Complete coding genome sequencing evaluated on three different sequencing platforms.

The diseases caused by capripoxviruses (CaPVs) are of major economic concern in sheep, goat and cattle as they are inexorably spreading into non-endemic regions. As CaPV strains are serologically indistinguishable and genetically highly homologous, typing closely related strains can only be achieved by whole genome sequencing. Unfortunately the number of publicly available genomes remains low as most sequencing methods rely on virus isolation. Therefore, we developed a robust, cost-effective and widely applicable method that allows to generate (nearly) complete CaPV genomes directly from clinical samples or commercial vaccine batches. A set of pan-CaPVs long-range PCRs spanning the entire genome was designed to generate PCR amplicons that can be sequenced on commonly used high-throughput sequencing platforms: MiSeq (Illumina), RSII (PacBio) and MinION (Oxford Nanopore Technologies). The robustness of the LR-PCR strategy was evaluated for all 3 members of CaPV directly from a variety of samples, including clinical samples (N = 7), vaccine batches (N = 6), and virus isolates (N = 2). The sequencing method described here allows to reconstruct (nearly) complete CaPV genomes in less than a week and will aid researchers studying closely-related CaPV strains worldwide.

Transcriptome analysis of Rafflesia cantleyi flower stages reveals insights into the regulation of senescence

Rafflesia is a unique plant species existing as a single flower and produces the largest flower in the world. While Rafflesia buds take up to 21 months to develop, its flowers bloom and wither within about a week. In this study, transcriptome analysis was carried out to shed light on the molecular mechanism of senescence in Rafflesia. A total of 53.3 million high quality reads were obtained from two Rafflesia cantleyi flower developmental stages and assembled to generate 64,152 unigenes. Analysis of this dataset showed that 5,166 unigenes were differentially expressed, in which 1,073 unigenes were identified as genes involved in flower senescence. Results revealed that as the flowers progress to senescence, more genes related to flower senescence were significantly over-represented compared to those related to plant growth and development. Senescence of the R. cantleyi flower activates senescence-associated genes in the transcription activity (members of the transcription factor families MYB, bHLH, NAC, and WRKY), nutrient remobilization (autophagy-related protein and transporter genes), and redox regulation (CATALASE). Most of the senescence-related genes were found to be differentially regulated, perhaps for the fine-tuning of various responses in the senescing R. cantleyi flower. Additionally, pathway analysis showed the activation of genes such as ETHYLENE RECEPTOR, ETHYLENE-INSENSITIVE 2, ETHYLENE-INSENSITIVE 3, and ETHYLENE-RESPONSIVE TRANSCRIPTION FACTOR, indicating the possible involvement of the ethylene hormone response pathway in the regulation of R. cantleyi senescence. Our results provide a model of the molecular mechanism underlying R. cantleyi flower senescence, and contribute essential information towards further understanding the biology of the Rafflesiaceae family.

Mining Unknown Porcine Protein Isoforms by Tissue-based Map of Proteome Enhances Pig Genome Annotation

A lack of the complete pig proteome has left a gap in our knowledge of the pig genome and has restricted the feasibility of using pigs as a biomedical model. In this study, we developed a tissue-based proteome map using 34 major normal pig tissues. A total of 5841 unknown protein isoforms were identified and systematically characterized, including 2225 novel protein isoforms, 669 protein isoforms from 460 genes symbolized beginning with LOC, and 2947 protein isoforms without clear NCBI annotation in the current pig reference genome. These newly identified protein isoforms were functionally annotated through profiling the pig transcriptome with high-throughput RNA sequencing of the same pig tissues, further improving the genome annotation of the corresponding protein-coding genes. Combining the well-annotated genes that have parallel expression pattern and subcellular witness, we predicted the tissue-related subcellularlocations and potential functions for these unknown proteins. Finally, we mined 3081 orthologous genes for 52.7% of unknown protein isoforms across multiple species, referring to 68 KEGG pathways as well as 23 disease signaling pathways. These findings provide valuable insights and a rich resource for enhancing studies of pig genomics and biology, as well as biomedical model application to human medicine.

Transcriptome Analysis Reveals Putative Target Genes of APETALA3-3 During Early Floral Development in Nigella damascena L

Even though petals are homoplastic structures, their identity consistently involves genes of the APETALA3 (AP3) lineage. However, the extent to which the networks downstream of AP3 are conserved in species with petals of different evolutionary origins is unknown. In Ranunculaceae, the specificity of the AP3-III lineage offers a great opportunity to identify the petal gene regulatory network in a comparative framework. Using a transcriptomic approach, we investigated putative target genes of the AP3-III ortholog NdAP3-3 in Nigella damascena at early developmental stages when petal identity is determined, and we compared our data with that from selected eudicot species. We generated a de novo reference transcriptome to carry out a differential gene expression analysis between the wild-type and mutant NdAP3-3 genotypes differing by the presence vs. absence of petals at early stages of floral development. Among the 1,620 genes that were significantly differentially expressed between the two genotypes, functional annotation suggested a large involvement of nuclear activities, including regulation of transcription, and enrichment in processes linked to cell proliferation. Comparing with Arabidopsis data, we found that highly conserved genes between the two species are enriched in homologs of direct targets of the AtAP3 protein. Integrating AP3-3 binding site data from another Ranunculaceae species, Aquilegia coerulea, allowed us to identify a set of 18 putative target genes that were conserved between the three species. Our results suggest that, despite the independent evolutionary origin of petals in core eudicots and Ranunculaceae, a small conserved set of genes determines petal identity and early development in these taxa.

Transcriptome profiling to identify tepal cell enlargement and pigmentation genes and the function of LtEXLB1 in Lilium tsingtauense

To understand the molecular mechanism underlying tepal development and pigmentation in Lilium tsingtauense Gilg, we performed whole-transcriptome profiles from closed buds at the greenish tepal stage (CBS), the full-bloom with un-horizontal tepal stage (UFS), and the completely opened bud with reflected tepal stage (RFS) of L. tsingtauense. More than 95699 transcripts were generated using a de novo assembly approach. Gene ontology and pathway analysis of the assembled transcripts revealed carbon metabolism is involved in tepal development and pigmentation. In total, 8171 differentially expression genes (DEGs) in three tepal stages were identified. Among these DEGs, ~994 genes putatively encoded transcription factors (TFs), whereas 693 putatively encoded protein kinases. Regarding hormone pathways, 51 DEGs involved in auxin biosynthesis and signalling and 10 DEGs involved in ethylene biosynthesis and signalling. We also isolated seven LtEXPANSINs, including four EXPAs, one EXPB, one EXLA and one EXLB. LtEXLB1 (GenBank: MN856627) was expressed at higher levels in UFS and RFS, compared with CBS. Silencing LtEXLB1 in leaf discs and tepals by virus-induced gene silencing significantly decreased cell expansion under rehydration conditions. Further analysis revealed that more cell numbers were existed in the abaxial and adaxial subepidermis in the silenced LtEXLB1 samples. As the first transcriptome of L. tsingtauense, the unigenes are a valuable resource for future studies on tepal development, and LtEXLB1 functions in cell expansion.

Skin transcriptional profiles in Oophaga poison frogs

Aposematic organisms advertise their defensive toxins to predators using a variety of warning signals, including bright coloration. While most Neotropical poison frogs (Dendrobatidae) rely on crypsis to avoid predators, Oophaga poison frogs from South America advertise their chemical defenses, a complex mix of diet-derived alkaloids, by using conspicuous hues. The present study aimed to characterize the skin transcriptomic profiles of South American Oophaga poison frogs. Our analyses showed very similar transcriptomic profiles for these closely related species in terms of functional annotation and relative abundance of gene ontology terms expressed. Analyses of expression profiles of Oophaga and available skin transcriptomes of cryptic anurans allowed us to propose initial hypotheses for the active sequestration of alkaloid-based chemical defenses and to highlight some genes that may be potentially involved in resistance mechanisms to avoid self-intoxication and skin coloration. In doing so, we provide an important molecular resource for the study of warning signals that will facilitate the assembly and annotation of future poison frog genomes.

Transcriptome profiling of raspberry (Rubus idaeus Var. Amira) in response to infection by tomato ringspot virus (ToRSV)

Raspberry (Rubus sp.) is a berries fruit with an ongoing agricultural and commercial interest due to its high contents of flavonoids and nutrients beneficial for human health. The growing demand for raspberries is facing great challenges associated mainly with the dispersal of diseases, which produces a decrease in productivity and fruit quality. A broad range of genomic resources is available for other Rosaceae species; however, genomic resources for species of the Rubus genus are still limited. Here, we characterize the transcriptome of the Rubus idaeus (Var. Amira) in order to 1) provide clues in the transcriptional changes of R. idaeus against tomato ringspot virus (ToRSV); and 2) generate genomic resources for this economically important species. We generate more than 200 million sequencing reads from two mRNA samples of raspberry, infected and not infected by ToRSV, using Illumina technology. After de novo assembly, we obtained 68,853 predicted protein-coding sequences of which 71.3% and 61.3% were annotated using Gene Ontology and Pfam databases, respectively. Moreover, we find 2,340 genes with differential expression between raspberries infected and not infected by ToRSV. Analysis of these genes shows functional enrichments of the oxidation-reduction process, cell wall biogenesis, terpene synthase activity, and lyase activity. These genes could be involved in the raspberry immune response through the interaction of different metabolic pathways; however, this statement needs further investigations. Up-regulation of genes encoding terpene synthases, multicopper oxidases, laccases, and beta-glucosidases might suggest that these enzymes appear to be the predominant transcriptome immune response of R. idaeus against ToRSV. Furthermore, we identify thousands of molecular markers (i.e., SSRs and SNPs), increasing considerably the genomic resources currently available for raspberries. This study is the first report on investigating the transcriptional changes of R. idaeus against ToRSV.

A zinc finger protein BBX19 interacts with ABF3 to affect drought tolerance negatively in chrysanthemum

Drought is an environmental factor that can severely influence plant development and distribution, and greatly affect the yield and economic value of crops. We characterized CmBBX19, a BBX family subgroup IV member gene, from the transcriptome database of Chrysanthemum morifolium in response to drought stress. Drought stress and ABA treatments downregulated the expression of CmBBX19. We generated CmBBX19-overexpressing (CmBBX19-OX) lines and CmBBX19-suppressing lines (CmBBX19-RNAi), and found that suppressed expression of CmBBX19 led to enhanced drought tolerance compared with the wild-type (WT) controls, while CmBBX19-OX lines exhibited reduced drought tolerance. Downstream gene analysis showed that CmBBX19 modulates drought tolerance mainly through inducing changes in the expression of ABA-dependent pathway genes, including protective protein, redox balance and cell wall biogenesis genes, such as responsive to ABA 18, peroxidase 12, and cellulose synthase-like protein G2. Moreover, CmBBX19 was shown to interact with CmABF3, a master ABA signaling component, to suppress expression of these downstream genes. We conclude that BBX19-ABF3 module functions in the regulation of drought tolerance of chrysanthemum through an ABA-dependent pathway.

The Penium margaritaceum Genome: Hallmarks of the Origins of Land Plants

The evolutionary features and molecular innovations that enabled plants to first colonize land are not well understood. Here, insights are provided through our report of the genome sequence of the unicellular alga Penium margaritaceum, a member of the Zygnematophyceae, the sister lineage to land plants. The genome has a high proportion of repeat sequences that are associated with massive segmental gene duplications, likely facilitating neofunctionalization. Compared with representatives of earlier diverging algal lineages, P. margaritaceum has expanded repertoires of gene families, signaling networks, and adaptive responses that highlight the evolutionary trajectory toward terrestrialization. These encompass a broad range of physiological processes and protective cellular features, such as flavonoid compounds and large families of modifying enzymes involved in cell wall biosynthesis, assembly, and remodeling. Transcriptome profiling further elucidated adaptations, responses, and selective pressures associated with the semi-terrestrial ecosystems of P. margaritaceum, where a simple body plan would be an advantage.

Cucurbit Genomics Database (CuGenDB): a central portal for comparative and functional genomics of cucurbit crops

The Cucurbitaceae family (cucurbit) includes several economically important crops, such as melon, cucumber, watermelon, pumpkin, squash and gourds. During the past several years, genomic and genetic data have been rapidly accumulated for cucurbits. To store, mine, analyze, integrate and disseminate these large-scale datasets and to provide a central portal for the cucurbit research and breeding community, we have developed the Cucurbit Genomics Database (CuGenDB; http://cucurbitgenomics.org) using the Tripal toolkit. The database currently contains all available genome and expressed sequence tag (EST) sequences, genetic maps, and transcriptome profiles for cucurbit species, as well as sequence annotations, biochemical pathways and comparative genomic analysis results such as synteny blocks and homologous gene pairs between different cucurbit species. A set of analysis and visualization tools and user-friendly query interfaces have been implemented in the database to facilitate the usage of these large-scale data by the community. In particular, two new tools have been developed in the database, a ‘SyntenyViewer’ to view genome synteny between different cucurbit species and an ‘RNA-Seq’ module to analyze and visualize gene expression profiles. Both tools have been packed as Tripal extension modules that can be adopted in other genomics databases developed using the Tripal system.

Transcriptome analysis of Catarina scallop (Argopecten ventricosus) juveniles treated with highly-diluted immunomodulatory compounds reveals activation of non-self-recognition system

Marine bivalve hatchery productivity is continuously challenged by apparition and propagation of new diseases, mainly those related to vibriosis. Disinfectants and antibiotics are frequently overused to prevent pathogen presence, generating a potential negative impact on the environment. Recently, the use of highly diluted compounds with immunostimulant properties in marine organisms has been trailed successfully to activate the self-protection mechanisms of marine bivalves. Despite their potential as immunostimulants, little is known about their way of action. To understand their effect, a comparative transcriptomic analysis was performed with Argopecten ventricosus juveniles. The experimental design consisted of four treatments formulated from pathogenic Vibrio lysates at two dilutions: [(T1) Vibrio parahaemolyticus and Vibrio alginolyticus 1D; (T2) V. parahaemolyticus and V. alginolyticus 7C]; minerals [(T3) PhA+SiT 7C], scorpion venom [(T4) ViT 31C]; and one control (C1) hydro-alcoholic solution (ethanol 1%). The RNA sequencing (RNAseq) analysis showed a higher modulation of differentially expressed genes (DEG) in mantle tissue compared to gill tissue. The scallops that showed a higher number of DEG related to immune response in mantle tissue corresponded to T1 (V. parahaemolyticus and V. alginolyticus lysate) and T3 (Silicea terra® - Phosphoric acid®). The transcriptome analysis allowed understanding some interactions between A. ventricosus juveniles and highly-diluted treatments.

The RhHB1/RhLOX4 module affects the dehydration tolerance of rose flowers (Rosa hybrida) by fine-tuning jasmonic acid levels

Phytohormones are key factors in plant responsiveness to abiotic and biotic stresses, and maintaining hormone homeostasis is critically important during stress responses. Cut rose (Rosa hybrida) flowers experience dehydration stress during postharvest handling, and jasmonic acid (JA) levels change as a result of this stress. However, how JA is involved in dehydration tolerance remains unclear. We investigated the functions of the JA- and dehydration-induced RhHB1 gene, which encodes a homeodomain-leucine zipper I γ-clade transcription factor, in rose flowers. Silencing RhHB1 decreased petal dehydration tolerance and resulted in a persistent increase in JA-Ile content and reduced dehydration tolerance. An elevated JA-Ile level had a detrimental effect on rose petal dehydration tolerance. RhHB1 was shown to lower the transient induction of JA-Ile accumulation in response to dehydration. In addition to transcriptomic data, we obtained evidence that RhHB1 suppresses the expression of the lipoxygenase 4 (RhLOX4) gene by directly binding to its promoter both in vivo and in vitro. We propose that increased JA-Ile levels weaken the capacity for osmotic adjustment in petal cells, resulting in reduced dehydration tolerance. In conclusion, a JA feedback loop mediated by an RhHB1/RhLOX4 regulatory module provides dehydration tolerance by fine-tuning bioactive JA levels in dehydrated flowers.

NAA and 6-BA promote accumulation of oleanolic acid by JA regulation in Achyranthes bidentata Bl

Application of plant growth regulators has become one of the most important means of improving yield and quality of medicinal plants. To understand the molecular basis of phytohormone-regulated oleanolic acid metabolism, RNA-seq was used to analyze global gene expression in Achyranthes bidentata treated with 2.0 mg/L 1-naphthaleneacetic acid (NAA) and 1.0 mg/L 6-benzyladenine (6-BA). Compared with untreated controls, the expression levels of 20,896 genes were significantly altered with phytohormone treatment. We found that 13071 (62.5%) unigenes were up-regulated, and a lot of differentially expressed genes involved in hormone or terpenoid biosynthesis, or transcription factors were significantly up-regulated. These results suggest that oleanolic acid biosynthesis induced by NAA and 6-BA occurs due to the expression of key genes involved in jasmonic acid signal transduction. This study is the first to analyze the production and hormonal regulation of medicinal A. bidentata metabolites at the molecular level. The results herein contribute to a better understanding of the regulation of oleanane-type triterpenoid saponins accumulation and define strategies to improve the yield of these useful metabolites.

Resistance to Sharka in Apricot: Comparison of Phase-Reconstructed Resistant and Susceptible Haplotypes of 'Lito' Chromosome 1 and Analysis of Candidate Genes

Sharka, a common disease among most stone fruit crops, is caused by the Plum Pox Virus (PPV). Resistant genotypes have been found in apricot (Prunus armeniaca L.), one of which-the cultivar 'Lito' heterozygous for the resistance-has been used to map a major quantitative trait locus (QTL) on linkage group 1, following a pseudo-test-cross mating design with 231 individuals. In addition, 19 SNP markers were selected from among the hundreds previously developed, which allowed the region to be limited to 236 kb on chromosome 1. A 'Lito' bacterial artificial chromosome (BAC) library was produced, screened with markers of the region, and positive BAC clones were sequenced. Resistant (R) and susceptible (S) haplotypes were assembled independently. To refine the assembly, the whole genome of 'Lito' was sequenced to high coverage (98×) using PacBio technology, enabling the development of a detailed assembly of the region that was able to predict and annotate the genes in the QTL region. The selected cultivar 'Lito' allowed not only to discriminate structural variants between the two haplotypic regions but also to distinguish specific allele expression, contributing towards mining the PPVres locus. In light of these findings, genes previously indicated (i.e., MATHd genes) to have a possible role in PPV resistance were further analyzed, and new candidates were discussed. Although the results are not conclusive, the accurate and independent assembly of R and S haplotypes of 'Lito' is a valuable resource to predict and test alternative transcription and regulation mechanisms underpinning PPV resistance.

Rosa hybrida RhERF1 and RhERF4 mediate ethylene- and auxin-regulated petal abscission by influencing pectin degradation

The timing of plant organ abscission is modulated by the balance of two hormones, ethylene and auxin, while the mechanism of organ shedding depends on the loss of middle lamella pectin in the abscission zone (AZ). However, the mechanisms involved in sensing the balance of auxin and ethylene and that affect pectin degradation during abscission are not well understood. In this study, we identified two members of the APETALA2/ethylene-responsive factor (AP2/ERF) transcription factor family in rose (Rosa hybrida), RhERF1 and RhERF4 which play a role in petal abscission. The expression of RhERF1 and RhERF4 was influenced by ethylene and auxin, respectively. Reduced expression of RhERF1 or RhERF4 was observed to accelerate petal abscission. Global expression analysis and real-time PCR assays revealed that RhERF1 and RhERF4 modulate the expression of genes encoding pectin-metabolizing enzymes. A reduction in the abundance of pectin epitopes was detected in the AZs of RhERF1 and RhERF4-silenced plants by immunofluorescence microscopy analysis. In addition, RhERF1 and RhERF4 were shown to bind to the promoter of the pectin-metabolizing gene β-GALACTOSIDASE 1 (RhBGLA1), and reduced expression of RhBGLA1 delayed petal abscission. We conclude that during petal abscission, RhERF1 and RhERF4 integrate and coordinate ethylene and auxin signals to modulate pectin metabolism, in part by regulating the expression of RhBGLA1.

Fugacium Spliced Leader Genes Identified from Stranded RNA-Seq Datasets

Trans-splicing mechanisms have been documented in many lineages that are widely distributed phylogenetically, including dinoflagellates. The spliced leader (SL) sequence itself is conserved in dinoflagellates, although its gene sequences and arrangements have diversified within or across different species. In this study, we present 18 Fugacium kawagutii SL genes identified from stranded RNA-seq reads. These genes typically have a single SL but can contain several partial SLs with lengths ranging from 103 to 292 bp. Unexpectedly, we find the SL gene transcripts contain sequences upstream of the canonical SL, suggesting that generation of mature transcripts will require additional modifications following trans-splicing. We have also identified 13 SL-like genes whose expression levels and length are comparable to Dino-SL genes. Lastly, introns in these genes were identified and a new site for Sm-protein binding was proposed. Overall, this study provides a strategy for fast identification of SL genes and identifies new sequences of F. kawagutii SL genes to supplement our understanding of trans-splicing.

Fat body-specific vitellogenin expression regulates host-seeking behaviour in the mosquito Aedes albopictus

The high vector competence of mosquitoes is intrinsically linked to their reproductive strategy because females need a vertebrate blood meal to develop large batches of eggs. However, the molecular mechanisms and pathways regulating mosquito host-seeking behaviour are largely unknown. Here, we test whether host-seeking behaviour may be linked to the female's energy reserves, with low energy levels triggering the search for a nutrient-rich blood meal. Our results demonstrate that sugar feeding delays host-seeking behaviour in the invasive tiger mosquito Aedes albopictus, but the levels of energy reserves do not correlate with changes in host-seeking behaviour. Using tissue-specific gene expression analyses, we show for the first time, to our knowledge, that sugar feeding alone induces a transient up-regulation of several vitellogenesis-related genes in the female fat body, resembling the transcriptional response after a blood meal. Specifically, high expression levels of a vitellogenin gene (Vg-2) correlated with the lowest host-seeking activity of sugar-fed females. Knocking down the Vg-2 gene via RNA interference (RNAi) restored host-seeking behaviour in these females, firmly establishing that Vg-2 gene expression has a pivotal role in regulating host-seeking behaviour in young Ae. albopictus females. The identification of a molecular mechanism regulating host-seeking behaviour in mosquitoes could pave the way for novel vector control strategies aiming to reduce the biting activity of mosquitoes. From an evolutionary perspective, this is the first demonstration of vitellogenin genes controlling feeding-related behaviours in nonsocial insects, while vitellogenins are known to regulate caste-specific foraging and brood-care behaviours in eusocial insects. Hence, this work confirms the key role of vitellogenin in controlling feeding-related behaviours in distantly related insect orders, suggesting that this function could be more ubiquitous than previously thought.

Transcriptome analysis provides insights into the regulation of metabolic processes during postharvest cold storage of loquat (Eriobotrya japonica) fruit

Loquat (Eriobotrya japonica) fruit accumulates lignin during postharvest storage under chilling conditions (0 °C), while low-temperature conditioning (LTC; 5 °C for 6 days followed by transfer to 0 °C) or heat treatment (HT; 40 °C for 4 h followed by transfer to 0 °C) can alleviate lignification. Here we compared transcriptome profiles of loquat fruit samples under LTC or HT to those stored at 0 °C at five time points from day 1 to day 8 after treatment. High-throughput transcriptome sequences were de novo assembled into 53,319 unique transcripts with an N50 length of 1306 bp. A total of 2235 differentially expressed genes were identified in LTC, and 1020 were identified in HT compared to 0 °C. Key genes in the lignin biosynthetic pathway, including EjPAL2, EjCAD1, EjCAD3, 4CL, COMT, and HCT, were responsive to LTC or HT treatment, but they showed different expression patterns during the treatments, indicating that different structural genes could regulate lignification at different treatment stages. Coexpression network analysis showed that these candidate biosynthetic genes were associated with a number of transcription factors, including those belonging to the AP2, MYB, and NAC families. Gene ontology (GO) enrichment analysis of differentially expressed genes indicated that biological processes such as stress responses, cell wall and lignin metabolism, hormone metabolism, and metal ion transport were significantly affected under LTC or HT treatment when compared to 0 °C. Our analyses provide insights into transcriptome responses to postharvest treatments in loquat fruit.

Injury affects coelomic fluid proteome of the common starfish, Asterias rubens

Echinoderms, possessing outstanding regenerative capabilities, provide a unique model system for the study of response to injury. However, little is known about the proteomic composition of coelomic fluid, an important biofluid circulating throughout the animal's body and reflecting the overall biological status of the organism. In this study, we used LC-MALDI tandem mass spectrometry to characterize the proteome of the cell-free coelomic fluid of the starfish Asterias rubens and to follow the changes occurring in response to puncture wound and blood loss. In total, 91 proteins were identified, of which 61 were extracellular soluble and 16 were bound to the plasma membrane. The most represented functional terms were 'pattern recognition receptor activity' and 'peptidase inhibitor activity'. A series of candidate proteins involved in early response to injury was revealed. Ependymin, β-microseminoprotein, serum amyloid A and avidin-like proteins, which are known to be involved in intestinal regeneration in the sea cucumber, were also identified as injury-responsive proteins. Our results expand the list of proteins potentially involved in defense and regeneration in echinoderms and demonstrate dramatic effects of injury on the coelomic fluid proteome.

GhNAC83 inhibits corm dormancy release by regulating ABA signaling and cytokinin biosynthesis in Gladiolus hybridus

Corm dormancy is an important trait for breeding in many bulb flowers, including the most cultivated Gladiolus hybridus. Gladiolus corms are modified underground stems that function as storage organs and remain dormant to survive adverse environmental conditions. Unlike seed dormancy, not much is known about corm dormancy. Here, we characterize the mechanism of corm dormancy release (CDR) in Gladiolus. We identified an important ABA (abscisic acid) signaling regulator, GhPP2C1 (protein phosphatase 2C1), by transcriptome analysis of CDR. GhPP2C1 expression increased during CDR, and silencing of GhPP2C1 expression in dormant cormels delayed CDR. Furthermore, we show that GhPP2C1 expression is directly regulated by GhNAC83, which was identified by yeast one-hybrid library screening. In planta assays show that GhNAC83 is a negative regulator of GhPP2C1, and silencing of GhNAC83 promoted CDR. As expected, silencing of GhNAC83 decreased the ABA level, but also dramatically increased cytokinin (CK; zeatin) content in cormels. Binding assays demonstrate that GhNAC83 associates with the GhIPT (ISOPENTENYLTRANSFERASE) promoter and negatively regulates zeatin biosynthesis. Taken together, our results reveal that GhNAC83 promotes ABA signaling and synthesis, and inhibits CK biosynthesis pathways, thereby inhibiting CDR. These findings demonstrate that GhNAC83 regulates the ABA and CK pathways, and therefore controls corm dormancy.

Cucurbit Genomics Database (CuGenDB): a central portal for comparative and functional genomics of cucurbit crops

The Cucurbitaceae family (cucurbit) includes several economically important crops, such as melon, cucumber, watermelon, pumpkin, squash and gourds. During the past several years, genomic and genetic data have been rapidly accumulated for cucurbits. To store, mine, analyze, integrate and disseminate these large-scale datasets and to provide a central portal for the cucurbit research and breeding community, we have developed the Cucurbit Genomics Database (CuGenDB; http://cucurbitgenomics.org) using the Tripal toolkit. The database currently contains all available genome and expressed sequence tag (EST) sequences, genetic maps, and transcriptome profiles for cucurbit species, as well as sequence annotations, biochemical pathways and comparative genomic analysis results such as synteny blocks and homologous gene pairs between different cucurbit species. A set of analysis and visualization tools and user-friendly query interfaces have been implemented in the database to facilitate the usage of these large-scale data by the community. In particular, two new tools have been developed in the database, a 'SyntenyViewer' to view genome synteny between different cucurbit species and an 'RNA-Seq' module to analyze and visualize gene expression profiles. Both tools have been packed as Tripal extension modules that can be adopted in other genomics databases developed using the Tripal system.

Long-Read Annotation: Automated Eukaryotic Genome Annotation Based on Long-Read cDNA Sequencing

Single-molecule full-length complementary DNA (cDNA) sequencing can aid genome annotation by revealing transcript structure and alternative splice forms, yet current annotation pipelines do not incorporate such information. Here we present long-read annotation (LoReAn) software, an automated annotation pipeline utilizing short- and long-read cDNA sequencing, protein evidence, and ab initio prediction to generate accurate genome annotations. Based on annotations of two fungal genomes (Verticillium dahliae and Plicaturopsis crispa) and two plant genomes (Arabidopsis [Arabidopsis thaliana] and Oryza sativa), we show that LoReAn outperforms popular annotation pipelines by integrating single-molecule cDNA-sequencing data generated from either the Pacific Biosciences or MinION sequencing platforms, correctly predicting gene structure, and capturing genes missed by other annotation pipelines.

The MADS transcription factor CmANR1 positively modulates root system development by directly regulating CmPIN2 in chrysanthemum

Plant root systems are essential for many physiological processes, including water and nutrient absorption. MADS-box transcription factor (TF) genes have been characterized as the important regulators of root development in plants; however, the underlying mechanism is largely unknown, including chrysanthemum. Here, it was found that the overexpression of CmANR1, a chrysanthemum MADS-box TF gene, promoted both adventitious root (AR) and lateral root (LR) development in chrysanthemum. Whole transcriptome sequencing analysis revealed a series of differentially expressed unigenes (DEGs) in the roots of CmANR1-transgenic chrysanthemum plants compared to wild-type plants. Functional annotation of these DEGs by alignment with Gene Ontology (GO) terms and biochemical pathway Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analysis indicated that CmANR1 TF exhibited "DNA binding" and "catalytic" activity, as well as participated in "phytohormone signal transduction". Both chromatin immunoprecipitation-polymerase chain reaction (ChIP-PCR) and gel electrophoresis mobility shift assays (EMSA) indicated the direct binding of CmPIN2 to the recognition site CArG-box motif by CmANR1. Finally, a firefly luciferase imaging assay demonstrated the transcriptional activation of CmPIN2 by CmANR1 in vivo. Overall, our results provide novel insights into the mechanisms of MADS-box TF CmANR1 modulation of both AR and LR development, which occurs by directly regulating auxin transport gene CmPIN2 in chrysanthemum.

Comparative RNA-Seq analysis reveals a critical role for brassinosteroids in rose (Rosa hybrida) petal defense against Botrytis cinerea infection

Background

One of the most popular ornamental plants worldwide, roses (Rosa sp.), are very susceptible to Botrytis gray mold disease. The necrotrophic infection of rose petals by B. cinerea causes the collapse and death of these tissues in both the growth and post-harvest stages, resulting in serious economic losses. To understand the molecular basis of rose resistance against B. cinerea, we profiled the petal transcriptome using RNA-Seq technology.

Results

We identified differentially transcribed genes (DTGs) in petals during B. cinerea infection at 30 h post inoculation (hpi) and/or 48 hpi. Gene ontology term enrichment and pathway analyses revealed that metabolic, secondary metabolite biosynthesis, plant-pathogen interaction, and plant hormone signal transduction pathways were involved. The expression of 370 cell-surface immune receptors was upregulated during infection. In addition, 188 genes encoding transcription factors were upregulated, particularly in the ERF, WRKY, bHLH, MYB, and NAC families, implying their involvement in resistance against B. cinerea. We further identified 325 upregulated DTGs in the hormone signal transduction pathways. Among them, the brassinosteroid (BR)-related genes were the most significantly enriched. To confirm the role of BR in Botrytis resistance, exogenous BR was applied to rose flowers before the inoculation of B. cinerea, which enhanced the defense response in these petals.

Conclusions

Our global transcriptome profiling provides insights into the complex gene regulatory networks mediating the rose petal response to B. cinerea. We further demonstrated the role of the phytohormone BR in the resistance of petals to necrotrophic fungal pathogens.

Electronic supplementary material

The online version of this article (10.1186/s12863-018-0668-x) contains supplementary material, which is available to authorized users.

Transcriptome sequencing and expression profiling of genes involved in the response to abiotic stress in Medicago ruthenica

Medicago ruthenica is a perennial forage legume with the remarkable ability to survive under unfavorable environmental conditions. It has been identified as an excellent species of Medicago that can adapt to various environmental stresses including low temperature, drought, and salinity. To investigate its potential as a genetic resource, we performed transcriptome sequencing and analysis in M. ruthenica under abiotic stresses. We generated >120 million reads from six cDNA libraries, resulting in 79,249 unique transcripts, most of which were highly similar to transcripts from M. truncatula (44,608, 56.3%) and alfalfa (M. sativa, 48,023, 60.6%). Based on gene expression profiles, 2,721 transcripts were identified as abiotic stress responsive genes which were predicted to be mainly involved in phytohormone signaling pathways, transcriptional regulation, and ROS-scavenging. These results suggest that they play critical roles in the response to abiotic stress. In summary, we identified genes in our transcriptome dataset involved in the regulation of the abiotic stress response in M. ruthenica which will provide a valuable resource for the future identification and functional analysis of candidate genes for adaption to unfavorable conditions. The genes identified here could be also useful for improving stress tolerance traits in alfalfa through molecular breeding in the future.

In Silico Genome-Wide Analysis Reveals the Potential Links Between Core Genome of Acidithiobacillus thiooxidans and Its Autotrophic Lifestyle

The coinage “pan-genome” was first introduced dating back to 2005, and was used to elaborate the entire gene repertoire of any given species. Core genome consists of genes shared by all bacterial strains studied and is considered to encode essential functions associated with species’ basic biology and phenotypes, yet its relatedness with bacterial lifestyle of the species remains elusive. We performed the pan-genome analysis of sulfur-oxidizing acidophile Acidithiobacillus thiooxidans as a case study to highlight species’ core genome and its relevance with autotrophic lifestyle of bacterial species. The mathematical modeling based on bacterial genomes of A. thiooxidans species, including a novel strain ZBY isolated from Zambian copper mine plus eight other recognized strains, was attempted to extrapolate the expansion of its pan-genome, suggesting that A. thiooxidans pan-genome is closed. Further investigation revealed a common set of genes, many of which were assigned to metabolic profiles, notably with respect to energy metabolism, amino acid metabolism, and carbohydrate metabolism. The predicted metabolic profiles of A. thiooxidans were characterized by the fixation of inorganic carbon, assimilation of nitrogen compounds, and aerobic oxidation of various sulfur species. Notably, several hydrogenase (H₂ase)-like genes dispersed in core genome might represent the novel classes due to the potential functional disparities, despite being closely related homologous genes that code for H₂ase. Overall, the findings shed light on the distinguishing features of A. thiooxidans genomes on a global scale, and extend the understanding of its conserved core genome pertaining to autotrophic lifestyle.

Cellular and molecular characterizations of a slow-growth variant provide insights into the fast growth of bamboo

Few studies about bamboo naturally occurring mutants have been reported so far. Using an integrated anatomy, mathematics and genomics approach, we systematically characterized the cellular and molecular basis underlying the abnormal internode development of Pseudosasa japonica var. tsutsumiana, a stable variant with dwarf and swollen internodes that are caused by the compressed spiral growth and the swollen cells in the bottom part of the internode. P. japonica var. tsutsumiana is a slow-growth variant with disorderly cell division and cell growth during the fast growth stage. Comparative transcriptome analysis identified a number of genes related to cell growth that were significantly down-regulated in the variant, including those related to auxin, vesicle transport, cytoskeleton organization and cell wall modification, consistent with the thinner cell walls and lower contents of cellulose that were found in the variant, which together with the mechanic force composed by the extrusion pressure from the neighboring fast growth cells and the weight pressure above the growing cells might finally result in the radial and irregular growth of the variant cells. This study provides key candidate genes involved in the fast growth of bamboo from a 'mutant' perspective, and supports a plausible mechanism underlying the dwarf and swollen internodes of P. japonica var. tsutsumiana.

A synthetic biology approach for consistent production of plant-made recombinant polyclonal antibodies against snake venom toxins

Antivenoms developed from the plasma of hyperimmunized animals are the only effective treatment available against snakebite envenomation but shortage of supply contributes to the high morbidity and mortality toll of this tropical disease. We describe a synthetic biology approach to affordable and cost-effective antivenom production based on plant-made recombinant polyclonal antibodies (termed pluribodies). The strategy takes advantage of virus superinfection exclusion to induce the formation of somatic expression mosaics in agroinfiltrated plants, which enables the expression of complex antibody repertoires in a highly reproducible manner. Pluribodies developed using toxin-binding genetic information captured from peripheral blood lymphocytes of hyperimmunized camels recapitulated the overall binding activity of the immune response. Furthermore, an improved plant-made antivenom (plantivenom) was formulated using an in vitro selected pluribody against Bothrops asper snake venom toxins and has been shown to neutralize a wide range of toxin activities and provide protection against lethal venom doses in mice.

Transcriptome sequence analysis and mining of SSRs in Jhar Ber (Ziziphus nummularia (Burm.f.) Wight & Arn) under drought stress

Ziziphus nummularia (Burm.f.) Wight & Arn., a perennial shrub that thrives in the arid regions, is naturally tolerant to drought. However, there are limited studies on the genomics of drought tolerance in Ziziphus sp. In this study, RNA-sequencing of one month old seedlings treated with PEG 6000 was performed using Roche GS-FLX454 Titanium pyrosequencing. A total of 367,176 raw sequence reads were generated, and upon adapter trimming and quality filtration 351,872 reads were assembled de novo into 32,739 unigenes. Further characterization of the unigenes indicated that 73.25% had significant hits in the protein database. Kyoto encyclopedia of genes and genomes database (KEGG) identified 113 metabolic pathways from the obtained unigenes. A large number of drought-responsive genes were obtained and among them differential gene expression of 16 highly induced genes was validated by qRT-PCR analysis. To develop genic-markers, 3,425 simple sequence repeats (SSRs) were identified in 2,813 unigene sequences. The data generated shall serve as an important reservoir for the identification and characterization of drought stress responsive genes for development of drought tolerant crops.

Transcriptome sequencing analysis of alfalfa reveals CBF genes potentially playing important roles in response to freezing stress

Alfalfa (Medicago sativa L.) is an important perennial forage, with high nutritional value, which is widely grown in the world. Because of low freezing tolerance, its distribution and production are threatened and limited by winter weather. To understand the complex regulation mechanisms of freezing tolerance in alfalfa, we performed transcriptome sequencing analysis under cold (4 °C) and freezing (-8 °C) stresses. More than 66 million reads were generated, and we identified 5767 transcripts differentially expressed in response to cold and/or freezing stresses. These results showed that these genes were mainly classified as response to stress, transcription regulation, hormone signaling pathway, antioxidant, nodule morphogenesis, etc., implying their important roles in response to cold and freezing stresses. Furthermore, nine CBF transcripts differentially expressed were homologous to CBF genes of Mt-FTQTL6 site, conferring freezing tolerance in M. truncatula, which indicated that a genetic mechanism controlling freezing tolerance was conservative between M. truncatula and M. sativa. In summary, this transcriptome dataset highlighted the gene regulation response to cold and/or freezing stresses in alfalfa, which provides a valuable resource for future identification and functional analysis of candidate genes in determining freezing tolerance.

Comparative transcriptomics reveal developmental turning points during embryogenesis of a hemimetabolous insect, the damselfly Ischnura elegans

Identifying transcriptional changes during embryogenesis is of crucial importance for unravelling evolutionary, molecular and cellular mechanisms that underpin patterning and morphogenesis. However, comparative studies focusing on early/embryonic stages during insect development are limited to a few taxa. Drosophila melanogaster is the paradigm for insect development, whereas comparative transcriptomic studies of embryonic stages of hemimetabolous insects are completely lacking. We reconstructed the first comparative transcriptome covering the daily embryonic developmental progression of the blue-tailed damselfly Ischnura elegans (Odonata), an ancient hemimetabolous representative. We identified a "core" set of 6,794 transcripts - shared by all embryonic stages - which are mainly involved in anatomical structure development and cellular nitrogen compound metabolic processes. We further used weighted gene co-expression network analysis to identify transcriptional changes during Odonata embryogenesis. Based on these analyses distinct clusters of transcriptional active sequences could be revealed, indicating that embryos at different development stages have their own transcriptomic profile according to the developmental events and leading to sequential reprogramming of metabolic and developmental genes. Interestingly, a major change in transcriptionally active sequences is correlated with katatrepsis (revolution) during mid-embryogenesis, a 180° rotation of the embryo within the egg and specific to hemimetabolous insects.

Transcriptomic information from Pacific white shrimp (Litopenaeus vannamei) ovary and eyestalk, and expression patterns for genes putatively involved in the reproductive process

The increased use of massive sequencing technologies has enabled the identification of several genes known to be involved in different mechanisms associated with reproduction that so far have only been studied in vertebrates and other model invertebrate species. In order to further investigate the genes involved in Litopenaeus vannamei reproduction, cDNA and SSH libraries derived from female eyestalk and gonad were produced, allowing the identification of expressed sequences tags (ESTs) that potentially have a role in the regulation of gonadal maturation. In the present study, different transcripts involved in reproduction were identified and a number of them were characterized as full-length. These transcripts were evaluated in males and females in order to establish their tissue expression profiles during developmental stages (juvenile, subadult and adult), and in the case of females, their possible association with gonad maturation was assessed through expression analysis of vitellogenin. The results indicated that the expression of vitellogenin receptor (vtgr) and minichromosome maintenance (mcm) family members in the female gonad suggest an important role during previtellogenesis. Additionally, the expression profiles of genes such as famet, igfbp and gpcr in brain tissues suggest an interaction between the insulin/insulin-like growth factor signaling pathway (IIS) and methyl farnesoate (MF) biosynthesis for control of reproduction. Furthermore, the specific expression pattern of farnesoic acid O-methyltransferase suggests that final synthesis of MF is carried out in different target tissues, where it is regulated by esterase enzymes under a tissue-specific hormonal control. Finally, the presence of a vertebrate type steroid receptor in hepatopancreas and intestine besides being highly expressed in female gonads, suggest a role of that receptor during sexual maturation.

Transcriptome Analysis of Mango (Mangifera indica L.) Fruit Epidermal Peel to Identify Putative Cuticle-Associated Genes

Mango fruit (Mangifera indica L.) are highly perishable and have a limited shelf life, due to postharvest desiccation and senescence, which limits their global distribution. Recent studies of tomato fruit suggest that these traits are influenced by the expression of genes that are associated with cuticle metabolism. However, studies of these phenomena in mango fruit are limited by the lack of genome-scale data. In order to gain insight into the mango cuticle biogenesis and identify putative cuticle-associated genes, we analyzed the transcriptomes of peels from ripe and overripe mango fruit using RNA-Seq. Approximately 400 million reads were generated and de novo assembled into 107,744 unigenes, with a mean length of 1,717 bp and with this information an online Mango RNA-Seq Database (http://bioinfo.bti.cornell.edu/cgi-bin/mango/index.cgi) which is a valuable genomic resource for molecular research into the biology of mango fruit was created. RNA-Seq analysis suggested that the pathway leading to biosynthesis of the cuticle component, cutin, is up-regulated during overripening. This data was supported by analysis of the expression of several putative cuticle-associated genes and by gravimetric and microscopic studies of cuticle deposition, revealing a complex continuous pattern of cuticle deposition during fruit development and involving substantial accumulation during ripening/overripening.

Transcriptome analysis of Polygonum minus reveals candidate genes involved in important secondary metabolic pathways of phenylpropanoids and flavonoids

BACKGROUND

Polygonum minus is an herbal plant in the Polygonaceae family which is rich in ethnomedicinal plants. The chemical composition and characteristic pungent fragrance of Polygonum minus have been extensively studied due to its culinary and medicinal properties. There are only a few transcriptome sequences available for species from this important family of medicinal plants. The limited genetic information from the public expressed sequences tag (EST) library hinders further study on molecular mechanisms underlying secondary metabolite production.

METHODS

In this study, we performed a hybrid assembly of 454 and Illumina sequencing reads from Polygonum minus root and leaf tissues, respectively, to generate a combined transcriptome library as a reference.

RESULTS

A total of 34.37 million filtered and normalized reads were assembled into 188,735 transcripts with a total length of 136.67 Mbp. We performed a similarity search against all the publicly available genome sequences and found similarity matches for 163,200 (86.5%) of Polygonum minus transcripts, largely from Arabidopsis thaliana (58.9%). Transcript abundance in the leaf and root tissues were estimated and validated through RT-qPCR of seven selected transcripts involved in the biosynthesis of phenylpropanoids and flavonoids. All the transcripts were annotated against KEGG pathways to profile transcripts related to the biosynthesis of secondary metabolites.

DISCUSSION

This comprehensive transcriptome profile will serve as a useful sequence resource for molecular genetics and evolutionary research on secondary metabolite biosynthesis in Polygonaceae family. Transcriptome assembly of Polygonum minus can be accessed at http://prims.researchfrontier.org/index.php/dataset/transcriptome.

Complete Genome Sequences of the Neethling-Like Lumpy Skin Disease Virus Strains Obtained Directly from Three Commercial Live Attenuated Vaccines

Lumpy skin disease virus (LSDV) causes an economically important disease in cattle. Here, we report the complete genome sequences of three LSDV strains obtained directly from the live attenuated vaccines: Lumpyvax (MSD Animal Health), Herbivac LS (Deltamune) and Lumpy Skin Disease Vaccine (Onderstepoort Biological Products).

VirusDetect: An automated pipeline for efficient virus discovery using deep sequencing of small RNAs

Accurate detection of viruses in plants and animals is critical for agriculture production and human health. Deep sequencing and assembly of virus-derived small interfering RNAs has proven to be a highly efficient approach for virus discovery. Here we present VirusDetect, a bioinformatics pipeline that can efficiently analyze large-scale small RNA (sRNA) datasets for both known and novel virus identification. VirusDetect performs both reference-guided assemblies through aligning sRNA sequences to a curated virus reference database and de novo assemblies of sRNA sequences with automated parameter optimization and the option of host sRNA subtraction. The assembled contigs are compared to a curated and classified reference virus database for known and novel virus identification, and evaluated for their sRNA size profiles to identify novel viruses. Extensive evaluations using plant and insect sRNA datasets suggest that VirusDetect is highly sensitive and efficient in identifying known and novel viruses. VirusDetect is freely available at http://bioinfo.bti.cornell.edu/tool/VirusDetect/.

Complete Genome Sequence of Capripoxvirus Strain KSGP 0240 from a Commercial Live Attenuated Vaccine

Capripoxviruses cause economically important diseases in domestic ruminants in regions endemic for these viruses. We report here the complete genome sequence of the KSGP 0240 vaccine strain from the live attenuated vaccine Kenyavac (JOVAC).

Complete Genome Sequence of the Goatpox Virus Strain Gorgan Obtained Directly from a Commercial Live Attenuated Vaccine

This is a report of the complete genome sequence of the goatpox virus strain Gorgan, which was obtained directly from a commercial live attenuated vaccine (Caprivac, Jordan Bio-Industries Centre).

Comparative Genomics of the Extreme Acidophile Acidithiobacillus thiooxidans Reveals Intraspecific Divergence and Niche Adaptation

Acidithiobacillus thiooxidans known for its ubiquity in diverse acidic and sulfur-bearing environments worldwide was used as the research subject in this study. To explore the genomic fluidity and intraspecific diversity of Acidithiobacillus thiooxidans (A. thiooxidans) species, comparative genomics based on nine draft genomes was performed. Phylogenomic scrutiny provided first insights into the multiple groupings of these strains, suggesting that genetic diversity might be potentially correlated with their geographic distribution as well as geochemical conditions. While these strains shared a large number of common genes, they displayed differences in gene content. Functional assignment indicated that the core genome was essential for microbial basic activities such as energy acquisition and uptake of nutrients, whereas the accessory genome was thought to be involved in niche adaptation. Comprehensive analysis of their predicted central metabolism revealed that few differences were observed among these strains. Further analyses showed evidences of relevance between environmental conditions and genomic diversification. Furthermore, a diverse pool of mobile genetic elements including insertion sequences and genomic islands in all A. thiooxidans strains probably demonstrated the frequent genetic flow (such as lateral gene transfer) in the extremely acidic environments. From another perspective, these elements might endow A. thiooxidans species with capacities to withstand the chemical constraints of their natural habitats. Taken together, our findings bring some valuable data to better understand the genomic diversity and econiche adaptation within A. thiooxidans strains.

De Novo Transcriptional Analysis of Alfalfa in Response to Saline-Alkaline Stress

Saline-alkaline stress, caused by high levels of harmful carbonate salts and high soil pH, is a major abiotic stress that affects crop productivity. Alfalfa is a widely cultivated perennial forage legume with some tolerance to biotic and abiotic stresses, especially to saline-alkaline stress. To elucidate the mechanism underlying plant saline-alkaline tolerance, we conducted transcriptome analysis of whole alfalfa seedlings treated with saline-alkaline solutions for 0 day (control), 1 day (short-term treatment), and 7 days (long-term treatment) using ion torrent sequencing technology. A transcriptome database dataset of 53,853 unigenes was generated, and 2,286 and 2,233 genes were differentially expressed in the short-term and long-term treatment, respectively. Gene ontology analysis revealed 14 highly enriched pathways and demonstrated the differential response of metabolic pathways between the short-term and long-term treatment. The expression levels of 109 and 96 transcription factors were significantly altered significantly after 1 day and 7 days of treatment, respectively. Specific responses of peroxidase, flavonoids, and the light pathway component indicated that the antioxidant capacity was one of the central mechanisms of saline-alkaline stress tolerance response in alfalfa. Among the 18 differentially expressed genes examined by real time PCR, the expression levels of eight genes, including inositol transporter, DNA binding protein, raffinose synthase, ferritin, aldo/keto reductase, glutathione S-transferase, xyloglucan endotrans glucosylase, and a NAC transcription factor, exhibited different patterns in response to saline and alkaline stress. The expression levels of the NAC transcription factor and glutathione S-transferase were altered significantly under saline stress and saline-alkaline stress; they were upregulated under saline-alkaline stress and downregulated under salt stress. Physiology assays showed an increased concentration of reactive oxygen species and malondialdehyde and a decreased content of chlorophyll, indicating that anti-oxidation and detoxification play an important role in response to saline-alkaline stress. Overall, the transcriptome analysis provided novel insights into the saline-alkaline stress tolerance response mechanisms in alfalfa.

Characterization of a New Pink-Fruited Tomato Mutant Results in the Identification of a Null Allele of the SlMYB12 Transcription Factor

The identification and characterization of new tomato (Solanum lycopersicum) mutants affected in fruit pigmentation and nutritional content can provide valuable insights into the underlying biology, as well as a source of new alleles for breeding programs. To date, all characterized pink-pigmented tomato fruit mutants appear to result from low SlMYB12 transcript levels in the fruit skin. Two new mutant lines displaying a pink fruit phenotype (pf1 and pf2) were characterized in this study. In the pf mutants, SlMYB12 transcripts accumulated to wild-type levels but exhibited the same truncation, which resulted in the absence of the essential MYB activation domain coding region. Allelism and complementation tests revealed that both pf mutants were allelic to the y locus and showed the same recessive null allele in homozygosis: Δy A set of molecular and metabolic effects, reminiscent of those observed in the Arabidopsis (Arabidopsis thaliana) myb11 myb12 myb111 triple mutant, were found in the tomato Δy mutants. To our knowledge, these have not been described previously, and our data support the idea of their being null mutants, in contrast to previously described transcriptional hypomorphic pink fruit lines. We detected a reduction in the expression of several flavonol glycosides and some associated glycosyl transferases. Transcriptome analysis further revealed that the effects of the pf mutations extended beyond the flavonoid pathway into the interface between primary and secondary metabolism. Finally, screening for Myb-binding sites in the candidate gene promoter sequences revealed that 141 of the 152 co-down-regulated genes may be direct targets of SlMYB12 regulation.

Evolutionary conservation of candidate osmoregulation genes in plant phloem sap-feeding insects

The high osmotic pressure generated by sugars in plant phloem sap is reduced in phloem-feeding aphids by sugar transformations and facilitated water flux in the gut. The genes mediating these osmoregulatory functions have been identified and validated empirically in the pea aphid Acyrthosiphon pisum: sucrase 1 (SUC1), a sucrase in glycoside hydrolase family 13 (GH13), and aquaporin 1 (AQP1), a member of the Drosophila integral protein (DRIP) family of aquaporins. Here, we describe molecular analysis of GH13 and AQP genes in phloem-feeding representatives of the four phloem-feeding groups: aphids (Myzus persicae), coccids (Planococcus citri), psyllids (Diaphorina citri, Bactericera cockerelli) and whiteflies (Bemisia tabaci MEAM1 and MED). A single candidate GH13-SUC gene and DRIP-AQP gene were identified in the genome/transcriptome of most insects tested by the criteria of sequence motif and gene expression in the gut. Exceptionally, the psyllid Ba. cockerelli transcriptome included a gut-expressed Pyrocoelia rufa integral protein (PRIP)-AQP, but has no DRIP-AQP transcripts, suggesting that PRIP-AQP is recruited for osmoregulatory function in this insect. This study indicates that phylogenetically related SUC and AQP genes may generally mediate osmoregulatory functions in these diverse phloem-feeding insects, and provides candidate genes for empirical validation and development as targets for osmotic disruption of pest species.

Identification of miRNAs and their targets through high-throughput sequencing and degradome analysis in male and female Asparagus officinalis

BACKGROUND

MicroRNAs (miRNAs), a class of non-coding small RNAs (sRNAs), regulate various biological processes. Although miRNAs have been identified and characterized in several plant species, miRNAs in Asparagus officinalis have not been reported. As a dioecious plant with homomorphic sex chromosomes, asparagus is regarded as an important model system for studying mechanisms of plant sex determination.

RESULTS

Two independent sRNA libraries from male and female asparagus plants were sequenced with Illumina sequencing, thereby generating 4.13 and 5.88 million final clean reads, respectively. Both libraries predominantly contained 24-nt sRNAs, followed by 21-nt sRNAs. Further analysis identified 154 conserved miRNAs, which belong to 26 families, and 39 novel miRNA candidates seemed to be specific to asparagus. Comparative profiling revealed that 63 miRNAs exhibited significant differential expression between male and female plants, which was confirmed by real-time quantitative PCR analysis. Among them, 37 miRNAs were significantly up-regulated in the female library, whereas the others were preferentially expressed in the male library. Furthermore, 40 target mRNAs representing 44 conserved and seven novel miRNAs were identified in asparagus through high-throughput degradome sequencing. Functional annotation showed that these target mRNAs were involved in a wide range of developmental and metabolic processes.

CONCLUSIONS

We identified a large set of conserved and specific miRNAs and compared their expression levels between male and female asparagus plants. Several asparagus miRNAs, which belong to the miR159, miR167, and miR172 families involved in reproductive organ development, were differentially expressed between male and female plants, as well as during flower development. Consistently, several predicted targets of asparagus miRNAs were associated with floral organ development. These findings suggest the potential roles of miRNAs in sex determination and reproductive developmental processes in asparagus.