High-throughput Tag Sequencing Research Articles

Gene expression profiling reveals the effects of light on adventitious root formation in lotus seedlings (Nelumbo nucifera Gaertn.)

BackgroundLotus is an aquatic horticultural crop that is widely cultivated in most regions of China and is used as an important off-season vegetable. The principal root of lotus is degenerated, and adventitious roots (ARs) are irreplaceable for plant growth. We found that no ARs formed under darkness and that exposure to high-intensity light significantly promoted the development of root primordia. Four differential expression libraries based on three light intensities were constructed to monitor metabolic changes, especially in indole-3-acetic acid (IAA) and sugar metabolism.ResultsAR formation was significantly affected by light, and high light intensity accelerated AR development. Metabolic changes during AR formation under different light intensities were evaluated using gene expression profiling by high-throughput tag-sequencing. More than 2.2 × 104 genes were obtained in each library; the expression level of most genes was between 0.01 and 100 (FPKF value). Libraries constructed from plants grown under darkness (D/CK), under 5000 lx (E/CK), and under 20,000 lx (F/CK) contained 1739, 1683, and 1462 upregulated genes and 1533, 995, and 834 downregulated genes, respectively, when compared to those in the initial state (CK). Additionally, we found that 1454 and 478 genes had altered expression in a comparison of libraries D/CK and F/CK. Gene transcription between libraries D/F ranged from a 5-fold decrease to a 5-fold increase. Twenty differentially expressed genes (DEGs) were involved in the signal transduction pathway, 28 DEGs were related to the IAA response, and 35 DEGs were involved in sugar metabolism. We observed that the IAA content was enhanced after seed germination, even in darkness; this was responsible for AR formation. We also observed that sucrose could eliminate the negative effect of 150 μMol IAA during AR development.ConclusionsAR formation was regulated by IAA, even in the dark, where induction and developmental processes could also be completed. In addition, 36 genes displayed altered expression in carbohydrate metabolism and ucrose metabolism was involved in AR development (expressed stage) according to gene expression and content change characteristics.

Comparative transcriptome analysis revealed the cooperative regulation of sucrose and IAA on adventitious root formation in lotus (Nelumbo nucifera Gaertn)

BackgroundIn China, lotus is an important cultivated crop with multiple applications in ornaments, food, and environmental purification. Adventitious roots (ARs), a secondary root is necessary for the uptake of nutrition and water as the lotus principle root is underdeveloped. Therefore, AR formation in seedlings is very important for lotus breeding due to its effect on plant early growth. As lotus ARs formation was significantly affected by sucrose treatment, we analyzed the expression of genes and miRNAs upon treatment with differential concentrations of sucrose, and a crosstalk between sucrose and IAA was also identified.ResultsNotably, 20 mg/L sucrose promoted the ARs development, whereas 60 mg/L sucrose inhibited the formation of ARs. To investigate the regulatory pathway during ARs formation, the expression of genes and miRNAs was evaluated by high-throughput tag-sequencing. We observed that the expression of 5438, 5184, and 5345 genes was enhanced in the GL20/CK0, GL60/CK0, and CK1/CK0 libraries, respectively. Further, the expression of 73, 78, and 71 miRNAs was upregulated in the ZT20/MCK0, ZT60/MCK0, and MCK1/MCK0 libraries, respectively. Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis revealed that most of the differentially expressed genes and miRNAs in the GL20/GL60 and ZT20/ZT60 libraries were involved in signal transduction. A large number of these genes (29) and miRNAs (53) were associated with plant hormone metabolism. We observed an association between five miRNAs (miR160, miR156a-5p, miR397-5p_1, miR396a and miR167d) and nine genes (auxin response factor, protein brassinosteroid insensitive 1, laccase, and peroxidase 27) in the ZT20/ ZT60 libraries during ARs formation. Quantitative polymerase chain reaction (qRT-PCR) was used to validate the high-throughput tag-sequencing data.ConclusionsWe found that the expression of many critical genes involved in IAA synthesis and IAA transport was changed after treatment with various concentration of sucrose. Based on the change of these genes expression, IAA and sucrose content, we concluded that sucrose and IAA cooperatively regulated ARs formation. Sucrose affected ARs formation by improving IAA content at induction stage, and increased sucrose content might be also required for ARs development according to the changes tendency after application of exogenous IAA.

Expression and functional analysis of the propamocarb-related gene CsMAPEG in cucumber

BackgroundPropamocarb (PM) is one of the main pesticides used for controlling cucumber downy mildew. However, due to its volatility and internal absorption, PM can easily form pesticide residues on cucumber fruits that seriously endanger human health and pollute the environment. The breeding of new cucumber varieties with a low abundance of PM residues via genetic methods constitutes an effective strategy for reducing pesticide residues and improving cucumber safety and quality. To help elucidate the molecular mechanism resulting in a low PM residue abundance in cucumber, we used the cucumber cultivar ‘D0351’ (which has the lowest PM residue content) as the test material and identified genes related to low PM residue abundance through high-throughput tag-sequencing (Tag-Seq).ResultsCsMAPEG was constitutively expressed and showed both varietal and organizational differences. This gene was strongly expressed in ‘D0351’. The expression levels of CsMAPEG in different cucumber tissues under PM stress were as follows: fruit>leaf>stem>root. CsMAPEG can respond to salicylic acid (SA), gibberellin (GA) and Corynespora cassiicola Wei (Cor) stress and thus plays an important regulatory role in plant responses to abiotic and biological stresses. The PM residue abundance in the fruits of CsMAPEG-overexpressing plants was lower than those found in antisense CsMAPEG plants and wild-type plants at all tested time points. The results revealed that CsMAPEG played a positive role in reducing the PM residue abundance. A CsMAPEG sense construct increased the contents of SOD, POD and GST in cucumber fruits, enhanced the degradation and metabolism of PM in cucumber, and thus effectively reduced the pesticide residue abundance in cucumber fruits.ConclusionsThe expression patterns of CsMAPEG in cucumber cultivars with high and low pesticide residue abundances and a transgenic verification analysis showed that CsMAPEG can actively respond to PM stress and effectively reduce the PM residue abundance in cucumber fruits. The results of this study will help researchers further elucidate the mechanism responsible for a low PM residue abundance in cucumber and lay a foundation for the breeding of new agricultural cucumber varieties with low pesticide residue abundances.

Expression and Functional Analysis of the Propamocarb-Related Gene CsMCF in Cucumber.

Propamocarb (PM), a carbamate fungicide, can effectively control downy mildew on cucumber. However, due to the large-scale and high-dose use of this fungicide, PM residues have become a major problem in cucumber production. In this report, the cucumber cultivar “D0351” (with the lowest residual PM content) and the cucumber cultivar “D9320” (with the highest residual PM content) were used as experimental materials. The candidate gene CsMCF, which is related to a low residual PM content in cucumber, was screened by high-throughput tag-sequencing (Tag-Seq) and PM analysis, and its role in reducing PM residue in cucumber was explored. CsMCF was cloned and obtained. This gene contains an open reading frame of 1026 bp, encodes 341 amino acids and contains 3 Mito-carr domains. The encoded protein is a hydrophobic protein with 4 distinct transmembrane structures but no signal peptide cleavage sites. The subcellular localization of the protein is the cytoplasm. Evolutionary tree analysis showed that CsMCF had the highest homology to a gene from the melon Cucumis melo L. (XM_008464998.2). The core elements of the promoter include cis-acting elements, such as those related to salicylic acid (SA), jasmonic acid (JA), gibberellin (GA), and abscisic acid (ABA). Following PM treatment, CsMCF was significantly upregulated at most time points in different parts of the fruit, leaf, stem and root of “D0351,” while expression was downregulated at most time points in the fruit, leaf and stem of “D9320.” The order of the expression levels in different cucumber organs was as follows: fruit>leaf > stem > root. CsMCF was specifically expressed in the stems and leaves of “D0351.” The PM residues in CsMCF (+)-overexpressing T0 and T1 cucumber fruits were significantly lower than those in the wild type, while the PM residues in CsMCF (-)-overexpressing T0 and T1 cucumber fruits were significantly higher. The qRT-PCR results showed that CsMCF can respond to biotic and abiotic stresses, actively respond to PM treatment and play a role in reducing PM residues in cucumber fruits.

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

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

Transcriptome analysis of miRNAs expression reveals novel insights into adventitious root formation in lotus (Nelumbo nucifera Gaertn.).

MicroRNA (miRNA)-regulated gene expression plays an important role in various plant metabolic processes. Although adventitious roots are critical to plant growth in lotus, the role of miRNA in AR formation remains unclear. Expression profiling of miRNAs was carried out during three different developmental stages of ARs in lotus: no induction of AR stage, initial stage of ARs, and maximum number of ARs. These data are referenced with the whole lotus genome as already identified through high-throughput tag-sequencing. 1.3 × 107 tags were achieved, of which 11,035,798, 11,436,062, and 12,542,392 clean tags were obtained from each stage, respectively. miRNA analysis revealed that miRNAs were less than 10% among all small RNAs. In total, 310 miRNAs (240 up-regulated and 70 down-regulated miRNAs) exhibited expression changes from the no induction stage to the initial stage. Moreover, expression of 140 miRNAs was increased and that of 123 miRNAs was decreased between the initial and maximum AR stages, mostly by ~ - 4-4-fold. miRNAs involved in metabolic pathways differed between the initial stage/no induction stage and the maximum number stage/initial stage. Several miRNAs in the initial stage/no induction stage were related to plant hormone metabolism and pyruvate and MAPK pathways, while major miRNAs in the maximum number stage/initial stage were involved in carbohydrate metabolism. All differentially expressed miRNAs associated with AR formation from the initial stage to maximum stage were also analyzed. The expression of 16 miRNAs was determined using qRT-PCR. This work provides a general insight into miRNA regulation during AR formation in lotus.

Transcriptome analysis of nitric oxide-responsive genes in upland cotton (Gossypium hirsutum).

Nitric oxide (NO) is an important signaling molecule with diverse physiological functions in plants. It is therefore important to characterize the downstream genes and signal transduction networks modulated by NO. Here, we identified 1,932 differentially expressed genes (DEGs) responding to NO in upland cotton using high throughput tag sequencing. The results of quantitative real-time polymerase chain reaction (qRT-PCR) analysis of 25 DEGs showed good consistency. Gene Ontology (GO) and KEGG pathway were analyzed to gain a better understanding of these DEGs. We identified 157 DEGs belonging to 36 transcription factor (TF) families and 72 DEGs related to eight plant hormones, among which several TF families and hormones were involved in stress responses. Hydrogen peroxide and malondialdehyde (MDA) contents were increased, as well related genes after treatment with sodium nitroprusside (SNP) (an NO donor), suggesting a role for NO in the plant stress response. Finally, we compared of the current and previous data indicating a massive number of NO-responsive genes at the large-scale transcriptome level. This study evaluated the landscape of NO-responsive genes in cotton and identified the involvement of NO in the stress response. Some of the identified DEGs represent good candidates for further functional analysis in cotton.

Expression and functional analysis of the Propamocarb-related gene CsDIR16 in cucumbers

BackgroundCucumber downy mildew is among the most important diseases that can disrupt cucumber production. Propamocarb, also known as propyl-[3-(dimethylamino)propyl]carbamate (PM), is a systemic carbamate fungicide pesticide that is widely applied in agricultural production because of its high efficiency of pathogens control, especially cucumber downy mildew. However, residual PM can remain in cucumbers after the disease has been controlled. To explore the molecular mechanisms of PM retention, cucumber cultivars ‘D9320’ (with the highest residual PM content) and ‘D0351’ (lowest residual PM content) were studied. High-throughput tag-sequencing (Tag-Seq) results showed that the CsDIR16 gene was related to PM residue, which was verified using transgenic technology.ResultsWe investigated the activity of a dirigent cucumber protein encoded by the CsDIR16 in gene response to stress induced by PM treatment. Gene-expression levels of CsDIR16 were up-regulated in the fruits, leaves, and stems of ‘D0351’ plants in response to PM treatment. However, in cultivar ‘D9320’, CsDIR16 levels were down-regulated in the leaves and stems after PM treatment, with no statistically significant differences observed in the fruits. Induction by jasmonic acid, abscisic acid, polyethylene glycol 4000, NaCl, and Corynespora cassiicola Wei (Cor) resulted in CsDIR16 up-regulation in ‘D0351’ and ‘D9320’. Expression after salicylic acid treatment was up-regulated in ‘D0351’, but was down-regulated in ‘D9320’. CsDIR16 overexpression lowered PM residues, and these were more rapidly reduced in CsDIR16(+) transgenic ‘D9320’ plants than in wild-type ‘D9320’ and CsDIR16(−) transgenic plants.ConclusionsAnalyses of the CsDIR16-expression patterns in the cucumber cultivars with the highest and lowest levels of PM residue, and transgenic validation indicated that CsDIR16 plays a positive role in reducing PM residues. The findings of this study help understand the regulatory mechanisms occurring in response to PM stress in cucumbers and in establishing the genetic basis for developing low-pesticide residue cucumber cultivars.

Photosynthetic Picoeukaryotes in the Land-Fast Ice of the White Sea, Russia.

The White Sea is a unique marine environment combining features of temperate and Arctic seas. The composition and abundance of photosynthetic picoeukaryotes (PPEs) were investigated in the land-fast ice of the White Sea, Russia, in March 2013 and 2014. High-throughput tag sequencing (Illumina MiSeq system) of the V4 region of the 18S rRNA gene was used to reveal the diversity of PPE ice community. The integrated PPE abundance varied from 11×106 cells/m2 to 364×106 cells/m2; the integrated biomass ranged from 0.02 to 0.26mg С/m2. The composition of sea-ice PPEs was represented by 16 algae genera belonging to eight classes and three super-groups. Chlorophyta, especially Mamiellophyceae, dominated among ice PPEs. The detailed analysis revealed the latent diversity of Micromonas and Mantоniella. Micromonas clade E2 revealed in the subarctic White Sea ice indicates that the area of distribution of this species is wider than previously thought. We suppose there exists a new Micromonas clade F. Micromonas clade C and Minutocellulus polymorphus were first discovered in the ice and extend the modern concept of sympagic communities' diversity generally and highlights the importance of further targeting subarctic sea ice for microbial study.

Combined use of network inference tools identifies ecologically meaningful bacterial associations in a paddy soil

High-throughput sequencing technologies have recently made it possible to interrogate the phylogenetic diversity of soils at considerable depth. This ability has led to the development of many computational tools to infer interaction networks from environmental samples. Although such tools have widely been used, they have more often served as a visual means to compare microbial communities across environmental gradients than as a means to appreciate microbial interactions associated with certain ecological processes. Previous studies have often regarded a subnetwork (module) as a functional unit but its functionality in ecological context has never been evidenced. To make better use of these tools in appreciating microbial interactions, we propose the combinational use of different inference tools. This ensemble approach permits the use of more independent predictors and the removal of tool-specific predictions in order to increase prediction accuracy. The purpose of the present study is to identify ecologically meaningful bacterial associations using multi-tool approach. Soil samples were collected in time series from experimental paddy rice plots. Bacterial communities were characterized by high-throughput tag sequencing of 16S rRNA gene fragments. We used three tools, Co-occurrence Network inference (CoNet), Molecular Ecological Network Analysis (MENA) and extended Local Similarity Analysis (eLSA), to infer networks from abundance profiles, partitioned the networks into modules, screened for the modules with ≥50% of genus-/species-level nodes, captured the modules that were derived from different tools and shared ≥ 50% of order-level nodes (tool-agreed modules) and tested their robustness against the changes in the tool parameters. By these procedures, two three-tool-agreed modules were found. One represented a guild that is phenotypically associated with aerobic respiration and fermentation and the other represented a guild phenotypically associated with metal/sulphur cycles, all of which are essential processes of water-submerged paddy soils that are mediated by bacteria. These data suggested that the linked members in a module were functionally associated taxa that work together to achieve a distinct function or an ecological process, and thus were ecologically meaningful to the environment. We selected three linked species from a three-tool-agreed module and validated their interactions using co-culture methods. Results showed that the interaction type between Janthinobacterium lividum and Leuconostoc lactis in the two-species mixture was validated to be ambivalent, positive for one partner and negative for the other. However, this type of interaction was not retained when a third party Lactococcus piscium was introduced, signifying the complexity of multi-species interactions. Validation results suggested that the selected species were interacting partners in laboratory but the validated interaction types were different from those inferred. By multi-tool approach, we also found that highly linked nodes, which are often referred to as “keystone species” and are frequently interpreted as the species playing important roles in soils, are tool dependent. Among top ten highly linked nodes, only four are conserved across three tools. These results suggest more research is required on the ecological significance of degree-based identification of keystone species. Overall, the present study highlights the potential utility of combined use of inference tools to identify ecologically meaningful bacterial associations in soils and other environmental samples. It is interesting to see what type of ecologically meaningful bacterial associations can be found in other soils.

Transcriptome Analysis of Gene Expression during Chinese Water Chestnut Storage Organ Formation.

The product organ (storage organ; corm) of the Chinese water chestnut has become a very popular food in Asian countries because of its unique nutritional value. Corm formation is a complex biological process, and extensive whole genome analysis of transcripts during corm development has not been carried out. In this study, four corm libraries at different developmental stages were constructed, and gene expression was identified using a high-throughput tag sequencing technique. Approximately 4.9 million tags were sequenced, and 4,371,386, 4,372,602, 4,782,494, and 5,276,540 clean tags, including 119,676, 110,701, 100,089, and 101,239 distinct tags, respectively, were obtained after removal of low-quality tags from each library. More than 39% of the distinct tags were unambiguous and could be mapped to reference genes, while 40% were unambiguous tag-mapped genes. After mapping their functions in existing databases, a total of 11,592, 10,949, 10,585, and 7,111 genes were annotated from the B1, B2, B3, and B4 libraries, respectively. Analysis of the differentially expressed genes (DEGs) in B1/B2, B2/B3, and B3/B4 libraries showed that most of the DEGs at the B1/B2 stages were involved in carbohydrate and hormone metabolism, while the majority of DEGs were involved in energy metabolism and carbohydrate metabolism at the B2/B3 and B3/B4 stages. All of the upregulated transcription factors and 9 important genes related to product organ formation in the above four stages were also identified. The expression changes of nine of the identified DEGs were validated using a quantitative PCR approach. This study provides a comprehensive understanding of gene expression during corm formation in the Chinese water chestnut.

Dynamic changes in the composition of photosynthetic picoeukaryotes in the northwestern Pacific Ocean revealed by high-throughput tag sequencing of plastid 16S rRNA genes

Photosynthetic picoeukaryotes (PPEs) are major oceanic primary producers. However, the diversity of such communities remains poorly understood, especially in the northwestern (NW) Pacific. We investigated the abundance and diversity of PPEs, and recorded environmental variables, along a transect from the coast to the open Pacific Ocean. High-throughput tag sequencing (using the MiSeq system) revealed the diversity of plastid 16S rRNA genes. The dominant PPEs changed at the class level along the transect. Prymnesiophyceae were the only dominant PPEs in the warm pool of the NW Pacific, but Mamiellophyceae dominated in coastal waters of the East China Sea. Phylogenetically, most Prymnesiophyceae sequences could not be resolved at lower taxonomic levels because no close relatives have been cultured. Within the Mamiellophyceae, the genera Micromonas and Ostreococcus dominated in marginal coastal areas affected by open water, whereas Bathycoccus dominated in the lower euphotic depths of oligotrophic open waters. Cryptophyceae and Phaeocystis (of the Prymnesiophyceae) dominated in areas affected principally by coastal water. We also defined the biogeographical distributions of Chrysophyceae, prasinophytes, Bacillariophyceaea and Pelagophyceae. These distributions were influenced by temperature, salinity and chlorophyll a and nutrient concentrations.

High-Throughput Tag-Sequencing Analysis of Early Events Induced by Ochratoxin A in HepG-2 Cells.

Ochratoxin A (OTA) is produced by fungi of the species Aspergillus and Penicillium. OTA has displayed hepatotoxicity in mammals. Although recent studies have indicated that OTA influences liver function, little is known regarding its impact on differential early liver toxicity. In this study, we report high-throughput tag-sequencing (Tag-seq) analysis of the transcriptome using Solexa Analyzer platform after 4 h of OTA treatment on HepG-2 cells. The analyses of differentially expressed genes revealed the substantial changes. A total of 21,449 genes were identified and quantified, with 2726 displaying significantly altered expression levels. Expression level data were then integrated with a network of gene-gene interactions, and biological pathways to obtain a systems-level view of changes in the transcriptome that occur with OTA resistance. Our data suggest that OTA exposure leads to an imbalance in zinc finger expression and shed light on splicing factor and mitochondrial-based mechanisms.

Identification of differentially expressed genes using digital gene expression profiles in Pyrus pyrifolia Nakai cv. Hosui bud release following early defoliation

“Returning bloom” is a major problem in the warm pear-producing regions during autumn–winter after the pear harvest and can cause tremendous pear production losses in the following year. The molecular mechanism causing this phenomenon is a mystery. We used high-throughput tag-sequencing (tag-seq) technology to analyze the gene expression profiling of pear (Pyrus pyrifolia Nakai) floral buds at three time points for completely defoliated plants after harvest. A total of 13.2–15.8 million reads were obtained for each of the six libraries, and approximately 61.81–67.11 % were mapped to Pyrus bretschneideri Rehd. genome sequences. A total of 2713 differentially expressed genes, categorized into 35 functional groups, were associated with 251 Kyoto Encyclopedia of Genes and Genomes pathways. Gene Ontology and pathway enrichment analysis showed differentially expressed genes linked to photosynthesis, hormone/signaling, and carbohydrate biosynthesis and metabolism. Six differentially expressed genes were selected for quantitative real-time PCR analysis, and the expression patterns were consistent with the tag-seq results. We concluded that bud release after artificial defoliation required compensatory mechanisms, which induce genes specifically expressed in association with the process of photosynthesis, endogenous hormones, carbohydrate, and reactive oxygen species. The differential expressions of genes found in pear tree defoliated by hand should play important roles in the reversion phase from paradormancy back to growth.

Transcriptome analysis of blood orange (Citrus sinensis) following fruit bagging treatment by digital gene expression profiling

SummaryBagging has been widely used in fruit crops to improve fruit quality. In this study, the internal qualities of blood orange (Citrus sinensis cv. ‘Tarocco’ ) fruit were significantly affected by bagging treatment. Total flavonoid (TFL) and total anthocyanin (TAN) concentrations increased, while total soluble solids (TSS), total phenolics (TPH), and ascorbic acid (AsA) concentrations and titratable acidity (TA) decreased in response to bagging treatment. To elucidate the molecular mode of action of bagging on fruit, high-throughput tag-sequencing (Tag-seq) analysis was used to analyse gene expression profiles in blood orange fruit in response to a bagging treatment. Over 21 × 106 total clean reads were obtained per library. Approx. 53.7 – 71.7% and 3.1 – 6.4% of the clean reads were mapped onto citrus genes or intergenic regions, respectively. However, approx. 25.2 – 39.9% of the clean reads failed to align with the citrus genome. Overall, bagging treatment resulted in an increase in transcripts involved in a range of metabolic pathways such as hormone synthesis and signalling, transcriptional regulation processes, the tricarboxylic acid (TCA) cycle, sucrose and starch metabolism, ascorbate metabolism, and the phenylpropanoid pathway, many of which are related to fruit quality. Analysis of the differentially expressed genes revealed that decreased TSS concentrations, TA, and AsA concentrations could mainly be attributed to the increased expression of genes involved in catabolic pathways, rather than in anabolic pathways. Competition for limited amounts of substrates for the flavonoid, phenolics, and anthocyanin pathways may have led to an increase in TAN and TFL concentrations, and a decrease in TPH concentration. In addition, our results showed that bagging treatment imposed a variety of stresses on fruit, thereby leading to an array of biochemical, physiological, and molecular changes. One gene (bHLH) that may be involved in anthocyanin biosynthesis in blood orange fruit has been identified and is discussed. These results reveal the extensive biochemical and genetic mechanisms operating in citrus fruit in response to a bagging treatment.

Transcriptome profilings of female Schistosoma japonicum reveal significant differential expression of genes after pairing

Pairing of Schistosoma japonicum initiates female development, leads to female sexual maturation, and maintains this mature state. To understand the mechanism involved in these processes, we studied parasites isolated from single- and double-sex cercariae-infected mice using deep-sequencing analysis, Solexa, to uncover pair-regulated transcriptional profiles. In this study, we report the results of high-throughput tag-sequencing (Tag-seq) analysis of the transcriptome of female worms 18 and 23 days postsingle- and double-sex infections. We sequenced over 3 million tags, obtained a total of 14,034, 27,251, 22,755, and 22,555 distinct tags corresponding to 5,773, 9,794, 8,885, and 8,870 tag-mapped genes for 23-day-old female schistosomula from double-sex infections (23DSI), 23-day-old female schistosomula from single-sex infections (23SSI), 18-day-old female schistosomula from double-sex infections (18DSI), and 18-day-old female schistosomula from single-sex infections (18SSI), respectively. Analyses of differentially expressed genes revealed similarities in the gene expression profiles between 18SSI and 18DSI as well as rational differential gene expression between 18SSI and 23SSI. However, fewer upregulated genes were found in 23DSI compared with 18DSI. Of the 3,446 differentially expressed genes between 23DSI and 23SSI, 2,913 genes were upregulated in 23SSI, whereas only 533 genes were upregulated in 23DSI. In these upregulated genes in 23DSI, phosphoglycerate mutase, superoxide dismutase, egg antigen, ribosomal proteins, ferritin-1 heavy chain, and eukaryotic translation initiation factor 2 were detected. Detection of these genes suggests that gene expression in 23DSI is specialized for functions such as promotion and maintenance of female sexual maturation and egg production. Quantitative real-time (RT)-PCR analysis confirmed the Solexa results, thereby supporting the reliability of the system. Our results offer new insights into the biological significance of pairing, which directs the expression of genes specific for sexual maturation and egg production.

Transcriptome analysis reveals differentially expressed genes associated with propamocarb response in cucumber (Cucumis sativus L.) fruit

High-throughput tag-sequencing (Tag-seq) from Illumina analysis, which is based on the Solexa Genome Analyzer platform, was applied to analyze the gene expression profiling of propamocarb (PM) treatment and control in cucumber fruit. Approximately 3.6 million complete clean sequence tags at PM treatment or control library were obtained with approximately 0.1 million distinct clean tag sequences. Approximately 41.79–43.15 % of the distinct clean tags were mapped unambiguously to the unigene database, and 32.54–33.46 % of the distinct clean tags were mapped to the cucumber genome database. The profiling analysis of the differentially expressed genes revealed the up-regulation of 546 genes and the down-regulation of 185 genes with PM response. Furthermore, the differentially expressed genes mainly linked to pesticide detoxication, response to stress/stimulus, transporter/signaling, and some important transcription factors. Finally, quantitative real-time polymerase chain reaction (qRT-PCR) and reverse transcription polymerase chain reaction (RT-PCR) using 16 genes independently verified the tag-mapped results. The present study reveals the comprehensive mechanisms of PM response in cucumber fruit.

Integrated detection of both 5-mC and 5-hmC by high-throughput tag sequencing technology highlights methylation reprogramming of bivalent genes during cellular differentiation

5-methylcytosine (5-mC) can be oxidized to 5-hydroxymethylcytosine (5-hmC). Genome-wide profiling of 5-hmC thus far indicates 5-hmC may not only be an intermediate form of DNA demethylation but could also constitute an epigenetic mark per se. Here we describe a cost-effective and selective method to detect both the hydroxymethylation and methylation status of cytosines in a subset of cytosines in the human genome. This method involves the selective glucosylation of 5-hmC residues, short-Sequence tag generation and high-throughput sequencing. We tested this method by screening H9 human embryonic stem cells and their differentiated embroid body cells, and found that differential hydroxymethylation preferentially occurs in bivalent genes during cellular differentiation. Especially, our results support hydroxymethylation can regulate key transcription regulators with bivalent marks through demethylation and affect cellular decision on choosing active or inactive state of these genes upon cellular differentiation. Future application of this technology would enable us to uncover the status of methylation and hydroxymethylation in dynamic biological processes and disease development in multiple biological samples.

Identification of differentially expressed genes relevant to corm formation in Sagittaria trifolia.

Sagittaria trifolia is a good model of wetland plants to elucidate the formation of corm. However, few studies have been conducted to uncover the complexity of gene expression involved in corm formation. In this study, high-throughput tag-sequencing based on Solexa Genome Analyzer Platform was applied to monitor the changes in gene expression with three libraries of differentially expressed genes (DEGs) (C1 library: stolon stage, C2 library: initial swelling stage and C3 library: swelling stage) during corm formation in Sagittaria trifolia. Approximately 6.0 million tags were sequenced, and 5854021, 5983454, and 5761079 clean tags including 138319, 116804, and 101739 distinct tags were obtained after removal of low quality tags from each library, respectively. About 46% distinct tags were unambiguous tags mapping to the reference genes, and 33% were unambiguous tag-mapped genes. Totally, 20575, 19807, and 18438 were annotated in C1, C2, and C3 libraries, respectively, after mapping their functions in existing databases. In addition, we found that profiling of gene expression in C1/C2 and C2/C3 libraries were different among most of the selected 20 DEGs. Most DEGs in C1/C2 libraries were relevant to hormone synthesis and response; energy metabolism and stress response, while most of the genes in C2/C3 libraries were involved in carbohydrate metabolism. All up-regulated transcriptional factors and 16 important genes relevant to corm formation in three libraries were also identified. To further analyze the expression of 9 genes, from the results of tag-sequencing, qRT-PCR was applied. In summary, this study provides a comprehensive understanding of gene expression, during the formation of corm in Sagittaria trifolia.

Protists in Arctic drift and land-fast sea ice.

Global climate change is having profound impacts on polar ice with changes in the duration and extent of both land-fast ice and drift ice, which is part of the polar ice pack. Sea ice is a distinct habitat and the morphologically identifiable sympagic community living within sea ice can be readily distinguished from pelagic species. Sympagic metazoa and diatoms have been studied extensively since they can be identified using microscopy techniques. However, non-diatom eukaryotic cells living in ice have received much less attention despite taxa such as the dinoflagellate Polarella and the cercozoan Cryothecomonas being isolated from sea ice. Other small flagellates have also been reported, suggesting complex microbial food webs. Since smaller flagellates are fragile, often poorly preserved, and are difficult for non-experts to identify, we applied high throughput tag sequencing of the V4 region of the 18S rRNA gene to investigate the eukaryotic microbiome within the ice. The sea ice communities were diverse (190 taxa) and included many heterotrophic and mixotrophic species. Dinoflagellates (43 taxa), diatoms (29 taxa) and cercozoans (12 taxa) accounted for ~80% of the sequences. The sympagic communities living within drift ice and land-fast ice harbored taxonomically distinct communities and we highlight specific taxa of dinoflagellates and diatoms that may be indicators of land-fast and drift ice.