Transcriptome sequencing and differential gene expression analysis of delayed gland morphogenesis in Gossypium australe during seed germination

Identification of MYB gene family and functional analysis of GhMYB4 in cotton (Gossypium spp.)

Myeloblastosis (MYB) transcription factors (TFs) form a large gene family involved in a variety of biological processes in plants. Little is known about their roles in the development of cotton pigment glands. In this study, 646 MYB members were identified in Gossypium hirsutum genome and phylogenetic classification was analyzed. Evolution analysis revealed assymetric evolution of GhMYBs during polyploidization and sequence divergence of MYBs in G. hirustum was preferentially happend in D sub-genome. WGCNA (weighted gene co-expression network analysis) showed that four modules had potential relationship with gland development or gossypol biosynthesis in cotton. Eight differentially expressed GhMYB genes were identified by screening transcriptome data of three pairs of glanded and glandless cotton lines. Of these, four were selected as candidate genes for cotton pigment gland formation or gossypol biosynthesis by qRT-PCR assay. Silencing of GH_A11G1361 (GhMYB4) downregulated expression of multiple genes in gossypol biosynthesis pathway, indicating it could be involved in gossypol biosynthesis. The potential protein interaction network suggests that several MYBs may have indirect interaction with GhMYC2-like, a key regulator of pigment gland formation. Our study was the systematic analysis of MYB genes in cotton pigment gland development, providing candidate genes for further study on the roles of cotton MYB genes in pigment gland formation, gossypol biosynthesis and future crop plant improvement.

Comparative Transcriptome Analysis Reveals Genes Associated with the Gossypol Synthesis and Gland Morphogenesis in Gossypium hirsutum

Gossypium hirsutum is an important source of natural textile fibers. Gossypol, which is a sesquiterpenoid compound mainly existing in the cotton pigment glands, can facilitate resistance to the stress from diseases and pests. The level of gossypol in the cotton is positively correlated to the quantity of pigment glands. However, the underlying regulatory mechanisms of gossypol synthesis and gland morphogenesis are still poorly understood, especially from a transcriptional perspective. The transcripts of young leaves and ovules at 30 DPA of the glanded plants and glandless plants were studied by RNA-Seq and 865 million clean reads were obtained. A total of 34,426 differentially expressed genes (DEGs) were identified through comparative transcriptome analysis. Genes related to gossypol synthesis or gland morphogenesis displayed significant differential expression between the two cultivars. Functional annotation revealed that the candidate genes related to catalytic activity, the biosynthesis of secondary metabolites, and biomolecular decomposition processes. Our work herein unveiled several potential candidate genes related to gossypol synthesis or gland morphogenesis and may provide useful clues for a breeding program of cotton cultivars with low cottonseed gossypol contents.

Identification of superior haplotypes in a diverse natural population for breeding desirable plant height in soybean

KEY MESSAGE

Plant height of soybean is associated with a haplotype block on chromosome 19, which classified 211 soybean accessions into five distinct groups showing significant differences for the target trait. Genetic variation is pivotal for crop improvement. Natural populations are precious genetic resources. However, efficient strategies for the targeted utilization of these resources for quantitative traits, such as plant height (PH), are scarce. Being an important agronomic trait associated with soybean yield and quality, it is imperative to unravel the genetic mechanisms underlying PH in soybean. Here, a genome-wide association study (GWAS) was performed to identify single nucleotide polymorphisms (SNPs) significantly associated with PH in a natural population of 211 cultivated soybeans, which was genotyped with NJAU 355 K Soy SNP Array and evaluated across six environments. A total of 128 SNPs distributed across 17 chromosomes were found to be significantly associated with PH across six environments and a combined environment. Three significant SNPs were consistently identified in at least three environments on Chr.02 (AX-93958260), Chr.17 (AX-94154834), and Chr.19 (AX-93897200). Genomic regions of ~ 130 kb flanking these three consistent SNPs were considered as stable QTLs, which included 169 genes. Of these, 22 genes (including Dt1) were prioritized and defined as putative candidates controlling PH. The genomic region flanking 12 most significant SNPs was in strong linkage disequilibrium (LD). These SNPs formed a single haplotype block containing five haplotypes for PH, namely Hap-A, Hap-B, Hap-C, Hap-D, and Hap-E. Deployment of such superior haplotypes in breeding programs will enable development of improved soybean varieties with desirable plant height.

An Overview of Cotton Gland Development and Its Transcriptional Regulation

Cotton refers to species in the genus Gossypium that bear spinnable seed coat fibers. A total of 50 species in the genus Gossypium have been described to date. Of these, only four species, viz. Gossypium, hirsutum, G. barbadense, G. arboretum, and G. herbaceum are cultivated; the rest are wild. The black dot-like structures on the surfaces of cotton organs or tissues, such as the leaves, stem, calyx, bracts, and boll surface, are called gossypol glands or pigment glands, which store terpenoid aldehydes, including gossypol. The cotton (Gossypium hirsutum) pigment gland is a distinctive structure that stores gossypol and its derivatives. It provides an ideal system for studying cell differentiation and organogenesis. However, only a few genes involved in the process of gland formation have been identified to date, and the molecular mechanisms underlying gland initiation remain unclear. The terpenoid aldehydes in the lysigenous glands of Gossypium species are important secondary phytoalexins (with gossypol being the most important) and one of the main defenses of plants against pests and diseases. Here, we review recent research on the development of gossypol glands in Gossypium species, the regulation of the terpenoid aldehyde biosynthesis pathway, discoveries from genetic engineering studies, and future research directions.

Transcriptome Analysis Identified Coordinated Control of Key Pathways Regulating Cellular Physiology and Metabolism upon Aspergillus flavus Infection Resulting in Reduced Aflatoxin Production in Groundnut

Aflatoxin-affected groundnut or peanut presents a major global health issue to both commercial and subsistence farming. Therefore, understanding the genetic and molecular mechanisms associated with resistance to aflatoxin production during host–pathogen interactions is crucial for breeding groundnut cultivars with minimal level of aflatoxin contamination. Here, we performed gene expression profiling to better understand the mechanisms involved in reduction and prevention of aflatoxin contamination resulting from Aspergillus flavus infection in groundnut seeds. RNA sequencing (RNA-Seq) of 16 samples from different time points during infection (24 h, 48 h, 72 h and the 7th day after inoculation) in U 4-7-5 (resistant) and JL 24 (susceptible) genotypes yielded 840.5 million raw reads with an average of 52.5 million reads per sample. A total of 1779 unique differentially expressed genes (DEGs) were identified. Furthermore, comprehensive analysis revealed several pathways, such as disease resistance, hormone biosynthetic signaling, flavonoid biosynthesis, reactive oxygen species (ROS) detoxifying, cell wall metabolism and catabolizing and seed germination. We also detected several highly upregulated transcription factors, such as ARF, DBB, MYB, NAC and C2H2 in the resistant genotype in comparison to the susceptible genotype after inoculation. Moreover, RNA-Seq analysis suggested the occurrence of coordinated control of key pathways controlling cellular physiology and metabolism upon A. flavus infection, resulting in reduced aflatoxin production.

Resistance effect of flavonols and toxicology analysis of hexabromocyclododecane based on soil-microbe-plant system

The toxicity characteristics of HBCD and resistance mechanism of flavonols are investigated based on physiological and metagenomic analysis. Toxicology research of HBCD on Arabidopsis thaliana (Col and fls1-3) not only shows the toxic effect of HBCD on plants, but also indicates that flavonols could improve plant resistance to HBCD, including root length, shoot biomass and chlorophyll content. Analysis of eggNOG and GO enrichment demonstrates that HBCD has toxic effect on both gene expression and protein function, which concentrates on energy production - conversion and amino acid transport - metabolism. Differential expressed genes in flavonols-treated groups indicates that flavonols regulate the metabolism of amino acids, cofactors and vitamins, which is involved in plant defense system against oxidative damage caused by HBCD stress. HBCD is believed to affect the synthesis of proteins via genes expression of ribosome biogenesis process. Flavonols could strengthen the plant resistance and alleviate toxic effect under HBCD stress.

The conserved regulatory basis of mRNA contributions to the early Drosophila embryo differs between the maternal and zygotic genomes

The gene products that drive early development are critical for setting up developmental trajectories in all animals. The earliest stages of development are fueled by maternally provided mRNAs until the zygote can take over transcription of its own genome. In early development, both maternally deposited and zygotically transcribed gene products have been well characterized in model systems. Previously, we demonstrated that across the genus Drosophila, maternal and zygotic mRNAs are largely conserved but also showed a surprising amount of change across species, with more differences evolving at the zygotic stage than the maternal stage. In this study, we use comparative methods to elucidate the regulatory mechanisms underlying maternal deposition and zygotic transcription across species. Through motif analysis, we discovered considerable conservation of regulatory mechanisms associated with maternal transcription, as compared to zygotic transcription. We also found that the regulatory mechanisms active in the maternal and zygotic genomes are quite different. For maternally deposited genes, we uncovered many signals that are consistent with transcriptional regulation at the level of chromatin state through factors enriched in the ovary, rather than precisely controlled gene-specific factors. For genes expressed only by the zygotic genome, we found evidence for previously identified regulators such as Zelda and GAGA-factor, with multiple analyses pointing toward gene-specific regulation. The observed mechanisms of regulation are consistent with what is known about regulation in these two genomes: during oogenesis, the maternal genome is optimized to quickly produce a large volume of transcripts to provide to the oocyte; after zygotic genome activation, mechanisms are employed to activate transcription of specific genes in a spatiotemporally precise manner. Thus the genetic architecture of the maternal and zygotic genomes, and the specific requirements for the transcripts present at each stage of embryogenesis, determine the regulatory mechanisms responsible for transcripts present at these stages.

Reconstruction of the full-length transcriptome atlas using PacBio Iso-Seq provides insight into the alternative splicing in Gossypium australe

BACKGROUND

Gossypium australe F. Mueller (2n = 2x = 26, G₂ genome) possesses valuable characteristics. For example, the delayed gland morphogenesis trait causes cottonseed protein and oil to be edible while retaining resistance to biotic stress. However, the lack of gene sequences and their alternative splicing (AS) in G. australe remain unclear, hindering to explore species-specific biological morphogenesis.

RESULTS

Here, we report the first sequencing of the full-length transcriptome of the Australian wild cotton species, G. australe, using Pacific Biosciences single-molecule long-read isoform sequencing (Iso-Seq) from the pooled cDNA of ten tissues to identify transcript loci and splice isoforms. We reconstructed the G. australe full-length transcriptome and identified 25,246 genes, 86 pre-miRNAs and 1468 lncRNAs. Most genes (12,832, 50.83%) exhibited two or more isoforms, suggesting a high degree of transcriptome complexity in G. australe. A total of 31,448 AS events in five major types were found among the 9944 gene loci. Among these five major types, intron retention was the most frequent, accounting for 68.85% of AS events. 29,718 polyadenylation sites were detected from 14,536 genes, 7900 of which have alternative polyadenylation sites (APA). In addition, based on our AS events annotations, RNA-Seq short reads from germinating seeds showed that differential expression of these events occurred during seed germination. Ten AS events that were randomly selected were further confirmed by RT-PCR amplification in leaf and germinating seeds.

CONCLUSIONS

The reconstructed gene sequences and their AS in G. australe would provide information for exploring beneficial characteristics in G. australe.

Identification of cotton MOTHER OF FT AND TFL1 homologs, GhMFT1 and GhMFT2, involved in seed germination

Plant phosphatidylethanolamine-binding protein (PEBP) is comprised of three clades: FLOWERING LOCUS T (FT), TERMINAL FLOWER1 (TFL1) and MOTHER OF FT AND TFL1 (MFT). FT/TFL1-like clades regulate identities of the determinate and indeterminate meristems, and ultimately affect flowering time and plant architecture. MFT is generally considered to be the ancestor of FT/TFL1, but its function is not well understood. Here, two MFT homoeologous gene pairs in Gossypium hirsutum, GhMFT1-A/D and GhMFT2-A/D, were identified by genome-wide identification of MFT-like genes. Detailed expression analysis revealed that GhMFT1 and GhMFT2 homoeologous genes were predominately expressed in ovules, and their expression increased remarkably during ovule development but decreased quickly during seed germination. Expressions of GhMFT1 and GhMFT2 homoeologous genes in germinating seeds were upregulated in response to abscisic acid (ABA), and their expressions also responded to gibberellin (GA). In addition, ectopic overexpression of GhMFT1 and GhMFT2 in Arabidopsis inhibited seed germination at the early stage. Gene transcription analysis showed that ABA metabolism genes ABA-INSENSITIVE3 (ABI3) and ABI5, GA signal transduction pathway genes REPRESSOR OF ga1-3 (RGA) and RGA-LIKE2 (RGL2) were all upregulated in the 35S:GhMFT1 and 35S:GhMFT2 transgenic Arabidopsis seeds. GhMFT1 and GhMFT2 localize in the cytoplasm and nucleus, and both interact with a cotton bZIP transcription factor GhFD, suggesting that both of GhMFT1, 2 have similar intracellular regulation mechanisms. Taken together, the results suggest that GhMFT1 and GhMFT2 may act redundantly and differentially in the regulation of seed germination.

Comparative transcriptome analyses on terpenoids metabolism in field- and mountain-cultivated ginseng roots

BACKGROUND

There exist differences in morphological traits and phytochemical compositions between field- and mountain-cultivated Panax ginseng (FCG and MCG), which might be attributed to variations of terpenoids metabolism adapting to different growth conditions. The present work aims to uncover these variations.

RESULTS

Among 26,648 differentially expressed genes, 496 genes distributed in seven dominant terpenoids pathways were identified. Diterpenoids and triterpenoids biosynthesis genes were significantly higher-expressed in FCG root. Conversely, biosynthesis of carotenoids was significantly more active in MCG root. Additionally, terpenoids backbones, monoterpenoids, sesquiterpenoids, and terpenoid-quinones biosyntheses were neither obviously inclined. Our determination also revealed that there were more gibberellins and steroids accumulated in FCG root which might be responsible for its quick vegetative growth, and enriched abscisic acid and germacrenes as well as protopanaxatriol-type ginsenosides might be major causes of enhanced stress-resistance in MCG root.

CONCLUSIONS

The study firstly provided an overview of terpenoids metabolism in roots of FCG and MCG in elucidating the underlying mechanisms for their different morphological appearances and phytochemical compositions.

Iron and callose homeostatic regulation in rice roots under low phosphorus

BACKGROUND

Phosphorus (Pi) deficiency induces root morphological remodeling in plants. The primary root length of rice increased under Pi deficiency stress; however, the underlying mechanism is not well understood. In this study, transcriptome analysis (RNA-seq) and Real-time quantitative PCR (qRT-PCR) techniques were combined with the determination of physiological and biochemical indexes to research the regulation mechanisms of iron (Fe) accumulation and callose deposition in rice roots, to illuminate the relationship between Fe accumulation and primary root growth under Pi deficient conditions.

RESULTS

Induced expression of LPR1 genes was observed under low Pi, which also caused Fe accumulation, resulting in iron plaque formation on the root surface in rice; however, in contrast to Arabidopsis, low Pi promoted primary root lengthening in rice. This might be due to Fe accumulation and callose deposition being still appropriately regulated under low Pi. The down-regulated expression of Fe-uptake-related key genes (including IRT, NAS, NAAT, YSLs, OsNRAMP1, ZIPs, ARF, and Rabs) inhibited iron uptake pathways I, II, and III in rice roots under low Pi conditions. In contrast, due to the up-regulated expression of the VITs gene, Fe was increasingly stored in both root vacuoles and cell walls. Furthermore, due to induced expression and increased activity of β-1-3 glucanase, callose deposition was more controlled in low Pi treated rice roots. In addition, low Pi and low Fe treatment still caused primary root lengthening.

CONCLUSIONS

The obtained results indicate that Low phosphorus induces iron and callose homeostatic regulation in rice roots. Because of the Fe homeostatic regulation, Fe plays a small role in rice root morphological remodeling under low Pi.

De Novo Transcriptome Assembly and Population Genetic Analyses for an Endangered Chinese Endemic Acer miaotaiense (Aceraceae)

Acer miaotaiense (P. C. Tsoong) is a rare and highly endangered plant in China. Because of the lack of genomic information and the limited number of available molecular markers, there are insufficient tools to determine the genetic diversity of this species. Here, 93,305 unigenes were obtained by multiple assembled contigs with a transcriptome sequencing program. Furthermore, 12,819 expressed sequence tag-derived simple sequence repeat (EST-SSR) markers were generated, 300 were randomly selected and synthesized, 19 primer pairs were identified as highly polymorphic (average number of alleles (N_a) = 8, expected heterozygosity (H_e) = 0.635, polymorphism information content (PIC) = 0.604) and were further used for population genetic analysis. All 261 samples were grouped into two genetic clusters by UPGMA, a principal component analyses and a STRUCTURE analyses. A moderate level of genetic differentiation (genetic differentiation index (F_st) = 0.059–0.116, gene flow = 1.904–3.993) among the populations and the major genetic variance (81.01%) within populations were revealed by the AMOVA. Based on the results, scientific conservation strategies should be established using in situ and ex situ conservation strategies. The study provides useful genetic information for the protection of precious wild resources and for further research on the origin and evolution of this endangered plant and its related species.

Genome-wide analysis of codon usage bias in four sequenced cotton species

Codon usage bias (CUB) is an important evolutionary feature in a genome which provides important information for studying organism evolution, gene function and exogenous gene expression. The CUB and its shaping factors in the nuclear genomes of four sequenced cotton species, G. arboreum (A₂), G. raimondii (D₅), G. hirsutum (AD₁) and G. barbadense (AD₂) were analyzed in the present study. The effective number of codons (ENC) analysis showed the CUB was weak in these four species and the four subgenomes of the two tetraploids. Codon composition analysis revealed these four species preferred to use pyrimidine-rich codons more frequently than purine-rich codons. Correlation analysis indicated that the base content at the third position of codons affect the degree of codon preference. PR2-bias plot and ENC-plot analyses revealed that the CUB patterns in these genomes and subgenomes were influenced by combined effects of translational selection, directional mutation and other factors. The translational selection (P2) analysis results, together with the non-significant correlation between GC12 and GC3, further revealed that translational selection played the dominant role over mutation pressure in the codon usage bias. Through relative synonymous codon usage (RSCU) analysis, we detected 25 high frequency codons preferred to end with T or A, and 31 low frequency codons inclined to end with C or G in these four species and four subgenomes. Finally, 19 to 26 optimal codons with 19 common ones were determined for each species and subgenomes, which preferred to end with A or T. We concluded that the codon usage bias was weak and the translation selection was the main shaping factor in nuclear genes of these four cotton genomes and four subgenomes.

Genome-wide comparative transcriptome analysis of CMS-D2 and its maintainer and restorer lines in upland cotton

BACKGROUND

Cytoplasmic male sterility (CMS) conferred by the cytoplasm from Gossypium harknessii (D2) is an important system for hybrid seed production in Upland cotton (G. hirsutum). The male sterility of CMS-D2 (i.e., A line) can be restored to fertility by a restorer (i.e., R line) carrying the restorer gene Rf1 transferred from the D2 nuclear genome. However, the molecular mechanisms of CMS-D2 and its restoration are poorly understood.

RESULTS

In this study, a genome-wide comparative transcriptome analysis was performed to identify differentially expressed genes (DEGs) in flower buds among the isogenic fertile R line and sterile A line derived from a backcross population (BC8F1) and the recurrent parent, i.e., the maintainer (B line). A total of 1464 DEGs were identified among the three isogenic lines, and the Rf1-carrying Chr_D05 and its homeologous Chr_A05 had more DEGs than other chromosomes. The results of GO and KEGG enrichment analysis showed differences in circadian rhythm between the fertile and sterile lines. Eleven DEGs were selected for validation using qRT-PCR, confirming the accuracy of the RNA-seq results.

CONCLUSIONS

Through genome-wide comparative transcriptome analysis, the differential expression profiles of CMS-D2 and its maintainer and restorer lines in Upland cotton were identified. Our results provide an important foundation for further studies into the molecular mechanisms of the interactions between the restorer gene Rf1 and the CMS-D2 cytoplasm.

ccNET: Database of co-expression networks with functional modules for diploid and polyploid Gossypium

Plant genera with both diploid and polyploid species are a common evolutionary occurrence. Polyploids, especially allopolyploids such as cotton and wheat, are a great model system for heterosis research. Here, we have integrated genome sequences and transcriptome data of Gossypium species to construct co-expression networks and identified functional modules from different cotton species, including 1155 and 1884 modules in G. arboreum and G. hirsutum, respectively. We overlayed the gene expression results onto the co-expression network. We further provided network comparison analysis for orthologous genes across the diploid and allotetraploid Gossypium. We also constructed miRNA-target networks and predicted PPI networks for both cotton species. Furthermore, we integrated in-house ChIP-seq data of histone modification (H3K4me3) together with cis-element analysis and gene sets enrichment analysis tools for studying possible gene regulatory mechanism in Gossypium species. Finally, we have constructed an online ccNET database (http://structuralbiology.cau.edu.cn/gossypium) for comparative gene functional analyses at a multi-dimensional network and epigenomic level across diploid and polyploid Gossypium species. The ccNET database will be beneficial for community to yield novel insights into gene/module functions during cotton development and stress response, and might be useful for studying conservation and diversity in other polyploid plants, such as T. aestivum and Brassica napus.

Global Transcriptomic Analysis Reveals the Mechanism of Phelipanche aegyptiaca Seed Germination

Phelipanche aegyptiaca is one of the most destructive root parasitic plants of Orobanchaceae. This plant has significant impacts on crop yields worldwide. Conditioned and host root stimulants, in particular, strigolactones, are needed for unique seed germination. However, no extensive study on this phenomenon has been conducted because of insufficient genomic information. Deep RNA sequencing, including de novo assembly and functional annotation was performed on P. aegyptiaca germinating seeds. The assembled transcriptome was used to analyze transcriptional dynamics during seed germination. Key gene categories involved were identified. A total of 274,964 transcripts were determined, and 53,921 unigenes were annotated according to the NR, GO, COG, KOG, and KEGG databases. Overall, 5324 differentially expressed genes among dormant, conditioned, and GR24-treated seeds were identified. GO and KEGG enrichment analyses demonstrated numerous DEGs related to DNA, RNA, and protein repair and biosynthesis, as well as carbohydrate and energy metabolism. Moreover, ABA and ethylene were found to play important roles in this process. GR24 application resulted in dramatic changes in ABA and ethylene-associated genes. Fluridone, a carotenoid biosynthesis inhibitor, alone could induce P. aegyptiaca seed germination. In addition, conditioning was probably not the indispensable stage for P. aegyptiaca, because the transcript level variation of MAX2 and KAI2 genes (relate to strigolactone signaling) was not up-regulated by conditioning treatment.

Differential Network Analysis Reveals Evolutionary Complexity in Secondary Metabolism of Rauvolfia serpentina over Catharanthus roseus

Comparative co-expression analysis of multiple species using high-throughput data is an integrative approach to determine the uniformity as well as diversification in biological processes. Rauvolfia serpentina and Catharanthus roseus, both members of Apocyanacae family, are reported to have remedial properties against multiple diseases. Despite of sharing upstream of terpenoid indole alkaloid pathway, there is significant diversity in tissue-specific synthesis and accumulation of specialized metabolites in these plants. This led us to implement comparative co-expression network analysis to investigate the modules and genes responsible for differential tissue-specific expression as well as species-specific synthesis of metabolites. Toward these goals differential network analysis was implemented to identify candidate genes responsible for diversification of metabolites profile. Three genes were identified with significant difference in connectivity leading to differential regulatory behavior between these plants. These genes may be responsible for diversification of secondary metabolism, and thereby for species-specific metabolite synthesis. The network robustness of R. serpentina, determined based on topological properties, was also complemented by comparison of gene-metabolite networks of both plants, and may have evolved to have complex metabolic mechanisms as compared to C. roseus under the influence of various stimuli. This study reveals evolution of complexity in secondary metabolism of R. serpentina, and key genes that contribute toward diversification of specific metabolites.

De novo assembly and characterization of germinating lettuce seed transcriptome using Illumina paired-end sequencing

At supraoptimal temperature, germination of lettuce (Lactuca sativa L.) seeds exhibits a typical germination thermoinhibition, which can be alleviated by sodium nitroprusside (SNP) in a nitric oxide-dependent manner. However, the molecular mechanism of seed germination thermoinhibition and its alleviation by SNP are poorly understood. In the present study, the lettuce seeds imbibed at optimal temperature in water or at supraoptimal temperature with or without 100 μM SNP for different periods of time were used as experimental materials, the total RNA was extracted and sequenced, we gained 147,271,347 raw reads using Illumina paired-end sequencing technique and assembled the transcriptome of germinating lettuce seeds. A total of 51,792 unigenes with a mean length of 849 nucleotides were obtained. Of these unigenes, a total of 29,542 unigenes were annotated by sequence similarity searching in four databases, NCBI non-redundant protein database, SwissProt protein database, euKaryotic Ortholog Groups database, and NCBI nucleotide database. Among the annotated unigenes, 22,276 unigenes were assigned to Gene Ontology database. When all the annotated unigenes were searched against the Kyoto Encyclopedia of Genes and Genomes Pathway database, a total of 8,810 unigenes were mapped to 5 main categories including 260 pathways. We first obtained a lot of unigenes encoding proteins involved in abscisic acid (ABA) signaling in lettuce, including 11 ABA receptors, 94 protein phosphatase 2Cs and 16 sucrose non-fermenting 1-related protein kinases. These results will help us to better understand the molecular mechanism of seed germination, thermoinhibition of seed germination and its alleviation by SNP.

Novel microRNA-like viral small regulatory RNAs arising during human hepatitis A virus infection

MicroRNAs (miRNAs), including host miRNAs and viral miRNAs, play vital roles in regulating host-virus interactions. DNA viruses encode miRNAs that regulate the viral life cycle. However, it is generally believed that cytoplasmic RNA viruses do not encode miRNAs, owing to inaccessible cellular miRNA processing machinery. Here, we provide a comprehensive genome-wide analysis and identification of miRNAs that were derived from hepatitis A virus (HAV; Hu/China/H2/1982), which is a typical cytoplasmic RNA virus. Using deep-sequencing and in silico approaches, we identified 2 novel virally encoded miRNAs, named hav-miR-1-5p and hav-miR-2-5p. Both of the novel virally encoded miRNAs were clearly detected in infected cells. Analysis of Dicer enzyme silencing demonstrated that HAV-derived miRNA biogenesis is Dicer dependent. Furthermore, we confirmed that HAV mature miRNAs were generated from viral miRNA precursors (pre-miRNAs) in host cells. Notably, naturally derived HAV miRNAs were biologically and functionally active and induced post-transcriptional gene silencing (PTGS). Genomic location analysis revealed novel miRNAs located in the coding region of the viral genome. Overall, our results show that HAV naturally generates functional miRNA-like small regulatory RNAs during infection. This is the first report of miRNAs derived from the coding region of genomic RNA of a cytoplasmic RNA virus. These observations demonstrate that a cytoplasmic RNA virus can naturally generate functional miRNAs, as DNA viruses do. These findings also contribute to improved understanding of host-RNA virus interactions mediated by RNA virus-derived miRNAs.

On the optimal trimming of high-throughput mRNA sequence data

The widespread and rapid adoption of high-throughput sequencing technologies has afforded researchers the opportunity to gain a deep understanding of genome level processes that underlie evolutionary change, and perhaps more importantly, the links between genotype and phenotype. In particular, researchers interested in functional biology and adaptation have used these technologies to sequence mRNA transcriptomes of specific tissues, which in turn are often compared to other tissues, or other individuals with different phenotypes. While these techniques are extremely powerful, careful attention to data quality is required. In particular, because high-throughput sequencing is more error-prone than traditional Sanger sequencing, quality trimming of sequence reads should be an important step in all data processing pipelines. While several software packages for quality trimming exist, no general guidelines for the specifics of trimming have been developed. Here, using empirically derived sequence data, I provide general recommendations regarding the optimal strength of trimming, specifically in mRNA-Seq studies. Although very aggressive quality trimming is common, this study suggests that a more gentle trimming, specifically of those nucleotides whose Phred score <2 or <5, is optimal for most studies across a wide variety of metrics.