The microbial community mediated biogeochemical cycles play important role in global C-cycle and display a sensitive response to environmental changes. Limited information is available on microbial composition and functional diversity controlling biogeochemical cycles in the riverine environment. The Ganga River water and sediment samples were studied for environmental gene tags with reference to carbohydrate metabolism, photoheterotrophy and chemolithotrophy using high throughput shotgun metagenomic sequencing and functional annotation. The diversity of environmental gene tags specific microbial community was annotated against reference sequence database using Kaiju taxonomic classifier. The metagenomic analyses revealed that the river harbored a broad range of carbohydrate and energy metabolism genes. The in-depth investigation of metagenomic data revealed that the enzymes associated with reverse TCA cycle, Calvin-Benson cycle enzyme RuBisCO, starch and sucrose metabolism genes were highly abundant. The enzymes associated with sulfur metabolism such as EC:2.7.7.4 (sulfate to ammonium per sulfate), EC:1.8.1.2, EC:1.8.7.1 (sulfite to H2S) were prevalent in both the class of samples. The principal component analysis of the functional profiles revealed that the water and sediment samples were clustered distinctly suggesting that both the sites had variable abundance of functional genes and associated microbiota. The taxonomic classification showed abundance of Proteobacteria, Actinobacteria and Bacteroidetes phyla. Also, the metagenomic study showed the presence of purple sulfur bacteria viz. Thiodictyon, Nitrosococcus and purple non-sulfur bacteria viz. Bradyrhozobium and Rhodobacter. The study demonstrates that the Ganga River microbiome has prevalence of functional genes involved in carbohydrate anabolism and catabolism, and CO2 fixation with great prospects in cellulose and sulfide degrading enzyme production and characterization.

Organisms use a variety of carbohydrates and metabolic pathways in order to capitalize in their specific environments. Depending upon their habitat, organism employs different types of transporters to maintain the cellular nutritional balance via central metabolism. A major contributor in this process in bacteria is a carbohydrate ABC transporter. The focus of this study is to get an insight into the carbohydrate transport and metabolism of a hot-spring-dwelling bacterium Thermus thermophilus HB8. We applied high-throughput data-mining approaches for identification and characterization of carbohydrate ABC transporters in T. thermophilus HB8. This enabled the identification of 11 putative carbohydrate ABC transport systems. To identify the cognate ligands for these transporters, functional annotation was performed. However, scarcity of homologous-protein's function hinders the process of functional annotation. Thus, to overcome this limitation, we integrated the functional annotation of carbohydrate ABC transporters with their metabolic analysis. Our results demonstrate that out of 11 putative carbohydrate ABC transporters, six are involved in the sugar (four for monosaccharides and polysaccharides-degraded products and two for osmotic regulation), four in phospholipid precursor (namely UgpABCE) and the remaining one in purine uptake. Further, analysis suggests the existence of sharing mechanism of transmembrane domains (TMDs) and/or nucleotide-binding domains (NBDs) among the 11 carbohydrate ABC transporters.

To investigate the variations in gene expression in grass carp under high-temperature stress, two libraries were constructed from a high-temperature treatment group (T33) and a control group (T27) and sequenced using Illumina sequencing technology. The results showed that sequencing generated a total of 279,398,348 raw reads, approximately 40.7-51.8 M clean reads were obtained from each library, and the percentage of uniquely mapped transcripts ranged from 80.13 to 84.58%. A total of 260 differentially expressed genes (DEGs) were identified under high-temperature stress, among which 84 genes were upregulated and 176 genes were downregulated. Ten DEGs were randomly selected for quantitative RT-PCR (qRT-PCR) analysis, and the results confirmed that the transcriptome analysis was reliable. Furthermore, the DEGs were subjected to Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analysis, and the results showed that most of the DEGs were involved in protein, lipid and carbohydrate metabolism. Moreover, plasma urea nitrogen (Urea) and triglyceride (TG) contents were significantly lower in the high-temperature treatment group than in the control group (P < 0.01). In summary, these results indicated that high-temperature stress could inhibit protein synthesis, decrease fatty acid synthesis, and weaken carbohydrate metabolism in juvenile grass carp.

Forkhead box O1 (foxo1) is a transcription factor and plays important roles in glucose metabolism. In the present study, foxo1 in turbot Scophthalmus maximus was cloned and characterized. The siRNA of foxo1 was used to investigate the functions of foxo1 in turbot hepatocytes glucose metabolism. After that, a 10-week feeding trial with two different dietary carbohydrate levels (15% and 21%, respectively) was conducted to analyze the function of foxo1 in glucose metabolism in vivo. Results showed that the foxo1 was identified as 2176 bp (base pair) with a 2025 bp open reading frame, which encoded 675 amino acids. Sequence analysis showed that foxo1 of turbot was highly homologous to most of fishes. Tissue distribution analysis revealed that the highest expression of foxo1 was in liver. After in vitro analysis, foxo1-specific small interfering RNA (sifoxo1) treatment significantly decreased the expressions of cytosolic phosphoenolpyruvate carboxykinase (cpepck) and glucose-6-phosphatase1(g6pase1) in primary hepatocytes. Expression of mitochondrial phosphoenolpyruvate carboxykinase (mpepck) was not significantly inhibited. In contrast, the expression of glucose-6-phosphatase2 (g6pase2) increased significantly. After the in vivo study (feeding trial), with the decreased expression of foxo1 in turbot due to high dietary carbohydrate level (21%), the expression of g6pase2 was significantly upregulated. However, the expression of glucokinase (gk) was not changed significantly. These increased the level of blood glucose and hepatic glycogen. In conclusion, data from both in vitro (primary hepatocytes) and in vivo (feeding trial) showed that downregulated foxo1 in turbot could not result in significant depression of gluconeogenesis and activation of glycolysis. This could be one of the reasons why intake of high level of carbohydrate resulted in prolonged hyperglycemia in turbot.

In old animals, bone's ability to adapt its mass and architecture to functional load-bearing requirements is diminished, resulting in bone loss characteristic of osteoporosis. Here we investigate transcriptomic changes associated with this impaired adaptive response. Young adult (19-week-old) and aged (19-month-old) female mice were subjected to unilateral axial tibial loading and their cortical shells harvested for microarray analysis between 1h and 24h following loading (36 mice per age group, 6 mice per loading group at 6 time points). In non-loaded aged bones, down-regulated genes are enriched for MAPK, Wnt and cell cycle components, including E2F1. E2F1 is the transcription factor most closely associated with genes down-regulated by ageing and is down-regulated at the protein level in osteocytes. Genes up-regulated in aged bone are enriched for carbohydrate metabolism, TNFα and TGFβ superfamily components. Loading stimulates rapid and sustained transcriptional responses in both age groups. However, genes related to proliferation are predominantly up-regulated in the young and down-regulated in the aged following loading, whereas those implicated in bioenergetics are down-regulated in the young and up-regulated in the aged. Networks of inter-related transcription factors regulated by E2F1 are loading-responsive in both age groups. Loading regulates genes involved in similar signalling cascades in both age groups, but these responses are more sustained in the young than aged. From this we conclude that cells in aged bone retain the capability to sense and transduce loading-related stimuli, but their ability to translate acute responses into functionally relevant outcomes is diminished.

Rainbow trout have, as salmonid fish species, undergone sequential genome duplication events in their evolutionary history. In addition to a teleost-specific whole genome duplication approximately 320-350 million years ago, rainbow trout and salmonids in general underwent an additional salmonid lineage-specific genome duplication event approximately 80 million years ago. Through the recent sequencing of salmonid genome sequences, including the rainbow trout, the identification and study of duplicated genes has become available. A particular focus of interest has been the evolution and regulation of rainbow trout gluconeogenic genes, as recent molecular and gene expression evidence points to a possible contribution of previously uncharacterized gluconeogenic gene paralogues to the rainbow trout long-studied glucose intolerant phenotype. Since the publication of the initial rainbow trout genome draft, resequencing and annotation have further improved genome coverage. Taking advantage of these recent improvements, we here identify a salmonid-specific genome duplication of ancestral mitochondrial phosphoenolpyruvate carboxykinase 2 isoenzyme, we termed pck2a and pck2b. Cytosolic phosphoenolpyruvate carboxykinase (Pck1) and, more recently mitochondrial Pck2, are considered to be the rate-limiting enzymes in de novo gluconeogenesis. Following in silico confirmation of salmonid pck2a and pck2b evolutionary history, we simultaneously profiled cytosolic pck1 and mitochondrial pck2a and pck2b expression in rainbow trout liver under several experimental conditions known to regulate hepatic gluconeogenesis. Cytosolic pck1 abundance was increased by nutritional (diets with a high protein to carbohydrate ratio compared to diets with a low carbohydrate to protein ratio) and glucoregulatory endocrine factors (glucagon and cortisol), revealing that the well-described transcriptional regulation of pck1 in mammals is present in rainbow trout. Conversely, and in contrast to mammals, we here describe endocrine regulation of pck2a (decrease in abundance in response to glucagon infusion), and nutritional, social-status-dependent and hypoxia-dependent regulation of pck2b. Specifically, pck2b transcript abundance increased in trout fed a diet with a low protein to carbohydrate ratio compared to a diet with a high protein to carbohydrate ratio, in dominant fish compared to subordinate fish as well as hypoxia. This specific and differential expression of rainbow trout pck2 ohnologues is indicative of functional diversification, and possible functional consequences are discussed in light of the recently highlighted gluconeogenic roles of mitochondrial pck2 in mammalian models.

Enterococcus faecalis is one of the main components of symbiotic bacteria in the intestine of silkworm (Bombyx mori L.). The abundance of E. faecalis in the intestine of silkworm is affected by fluoride exposure. However, the response mechanism of E. faecalis toward fluoride remains largely unknown. In this study, a strain of E. faecalis (named TV4), which is a symbiotic bacteria of silkworm, was isolated and characterized. Inhibition assay showed that fluoride can significantly inhibit the growth of the TV4 strain (P < 0.05) after culture for 4 h. Finally, Illumina X-Ten platform was used to investigate the response mechanism of E. faecalis TV4 under fluoride exposure. We found that the TV4 strain demonstrated significant changes in its carbohydrate transport and metabolism and energy metabolism. The transcriptome sequencing results revealed that 237 genes were differentially expressed for TV4 grown after fluoride exposure, i.e., 92 genes were differentially up-regulated and 145 genes were differentially down-regulated. Many of the down-regulated genes were involved in cell carbohydrate transport and metabolism and energy production, whereas the up-regulated genes were mostly related to ethanolamine utilization and amino acid synthesis and metabolism. Our results revealed that strain TV4 reduced its carbohydrate metabolism and energy metabolism and increased ethanolamine utilization and amino acid metabolism to adapt and survive under fluoride exposure. This study enhances our understanding about the response mechanism of E. faecalis after fluoride exposure and has important implications for investigations on the three-way interaction among fluoride, symbiotic bacteria and silkworm.

Peroxisome proliferator-activated receptors (PPARs) include the nuclear receptor superfamily of ligand-activated transcription factors involved in several metabolic processes, including carbohydrate and lipid metabolism.

Ganoderma lucidum is a valuable basidiomycete with numerous pharmacological compounds, which is widely consumed throughout China. We previously found that the polysaccharide content of Ganoderma lucidum fruiting bodies could be significantly improved by 45.63% with treatment of 42 °C heat stress (HS) for 2 h. To further investigate genes involved in HS response and explore the mechanisms of HS regulating the carbohydrate metabolism in Ganoderma lucidum, high-throughput RNA-Seq was conducted to analyse the difference between control and heat-treated mycelia at transcriptome level. We sequenced six cDNA libraries with three from control group (mycelia cultivated at 28 °C) and three from heat-treated group (mycelia subjected to 42 °C for 2 h). A total of 99,899 transcripts were generated using Trinity method and 59,136 unigenes were annotated by seven public databases. Among them, 2790 genes were identified to be differential expressed genes (DEGs) under HS condition, which included 1991 up-regulated and 799 down-regulated. 176 DEGs were then manually classified into five main responsive-related categories according to their putative functions and possible metabolic pathways. These groups include stress resistance-related factors; protein assembly, transportation and degradation; signal transduction; carbohydrate metabolism and energy provision-related process; other related functions, suggesting that a series of metabolic pathways in Ganoderma lucidum are activated by HS and the response mechanism involves a complex molecular network which needs further study. Remarkably, 48 DEGs were found to regulate carbohydrate metabolism, both in carbohydrate hydrolysis for energy provision and polysaccharide synthesis. In summary, this comprehensive transcriptome analysis will provide enlarged resource for further investigation into the molecular mechanisms of basidiomycete under HS condition.

The fish insulin/insulin-like growth factor (IGF) pathway has weak control over carbohydrate metabolism. To understand the molecular basis for the metabolic diversity, we characterized the forkhead box transcription factor O1A (FoxO1A), a downstream target of the insulin/IGF pathway, in torafugu Takifugu rubripes. The cloned torafugu FoxO1A cDNA contained all conserved features critical for its transcriptional activity and a unique unspliced intron encoding a poly-glutamine stretch. Torafugu FoxO1A showed the IGF-dependent nuclear exclusion and in vitro binding to the well-conserved FoxO1 binding site, DAF-16 binding element (DBE), but failed to bind to the insulin-responsive element by which mammalian FoxO1 mediates insulin effects. The subsequent in silico genomic screening provided a list of 587 potential torafugu FoxO1A target genes containing the DBE. Some carbohydrate metabolic genes regulated by FoxO1 in mammals were not included in the list. We further identified about 250 potential fish FoxO1 target genes by integrating results of the DBE screening against fish metagenome that contained 262 species. Neuronal processes appeared to be the common major function of fish FoxO1, although further annotation of the potential target genes is required. These results provide a part of the molecular basis underlying the weak association between the insulin/IGF pathway and carbohydrate metabolism in fish.

It is known that the male hypogonadism plays an important role in regulating adipose metabolism. In the present study, fifteen pairs of full male sibs were divided into a castrated group and an intact group with a paired experiment design. The pigs were slaughtered at an age of 175days. The carcass characteristics and fat deposit of the studied animal were measured, and the hormone and serum lipid levels of the peripheral blood samples were determined, and the differentially expressed genes of the back fat between the two groups were screened with porcine genome array. Our results showed that the absence of male gonadal steroids attributed to castration significantly raised the serum lipid levels and increased fat accumulation in the pigs. A total of 225 differentially expressed genes were identified between the boars and barrows and 135 of them were upregulated. The analysis of Gene Ontology categories and KEGG pathway indicated that these differentially expressed genes were mainly involved in metabolism of lipid, carbohydrate, amino acid, xenobiotics biodegradation, and immune diseases pathways. Our results indicated that there were higher capacity of fatty acid of synthesis, enhanced uptaking capacity of fatty acids and cholesterol, inhibited lipolysis, and enhanced carbohydrate metabolism in the adipose tissue of barrows compared to boars. The findings of the present study provide new insight into the mechanisms of adipose metabolism induced by hormone deficiency.

The human gut microbiota in long-living people has been characterized, however, its metabolic potential is still largely unknown in this group. In this study, the gut microbiota was assessed in 37 Chinese long-living participants (aged 90 + years) by metagenomic sequencing of stool samples. Participants were categorized into two groups, healthy long-living (n = 28) and unhealthy long-living (n = 9). Gut microbiota composition and function were compared among these two groups. We found that the gut microbiota in the healthy long-living group was significantly separated from the unhealthy group. The healthy long-living group contained a higher abundance of Bacteroidetes and more functional pathways in energy metabolism, glycan biosynthesis and metabolism, metabolism of cofactors and vitamins, and biosynthesis of other secondary metabolites. The unhealthy group contained a higher abundance of Streptococcus and other pathogenic bacteria, and also contained more functional pathways for xenobiotics biodegradation and metabolism than the healthy group. Additionally, the unhealthy group had decreased levels of carbohydrate-active enzymes, including host-glycan and fiber degrading enzymes, and an increase in starch-degrading enzymes. In conclusion, the gut microbiota of unhealthy long-living people contains more pathogenic bacteria, and the overall gut microbiota may be in an unhealthy state, "dysbiosis", which leads to a decrease in carbohydrate digestion, glycan and thiamine (B1) metabolites, and fatty acid biosynthesis.

n-3 highly unsaturated fatty acids (n-3 HUFAs) have been shown to suppress lipid accumulation and improve protein utilization in grass carp; however, little is known about the underlying molecular mechanism. Hence, we analyzed the hepatopancreas transcriptome of grass carp (Ctenopharyngodon idellus) fed either lard oil (LO) or fish oil (FO) diets. RNA-seq data showed that 125 genes were significantly up-regulated and 107 were significantly down-regulated in the FO group. Among them, 17 lipid metabolism related genes, 12 carbohydrate metabolism related genes, and 34 protein metabolism related genes were selected. Lipid metabolism related genes, such as very long-chain acyl-CoA synthetase (ACSVL),carnitine O-palmitoyltransferase 1 (CPT1) and carnitine-acylcarnitine translocase (CACT), were up-regulated in the FO group. But the genes of diacylglycerol O-acyltransferase 2 (DGAT2) and stearoyl-CoA desaturase (SCD) were down-regulated. Down-regulation of glycolysis related genes, such as 6-phosphofructokinase (PFK), phosphoglycerate kinase (PGK) and pyruvate dehydrogenase kinase (PDK), added with up-regulation of gluconeogenesis related genes, such as phosphoenolpyruvate carboxykinase (PEPCK) and glucose-6-phosphatase (G6Pase), suggests lower utilization of carbohydrate of the FO group. Besides, dietary FO also influenced the protein metabolism related genes, such as up-regulation of genes involved in digestion of dietary protein, mRNA transcription, protein translation and amino acid utilization, down-regulation of genes involved in mRNA degradation and ubiquitination of protein. Interestingly, the up-regulation of mitochondrial uncoupling protein 2 (UCP2) and down-regulation of oxidative phosphorylation related genes (cytochrome c oxidase subunit 4 isoform 2 [COX4I2], HIG1 domain family member 1A [HIGD1A] and cytochrome-b5 reductase [CYB5R]) suggest that energy metabolism may be also influenced by dietary fatty acid composition. These findings presented here provide a comprehensive understanding of the molecular mechanisms governing the effects of fish oil in grass carp.

Taenia pisiformis is one of the most widespread gastrointestinal parasites and its larvae (cysticercosis) causes significant economic loss to rabbit industry. No efficient drug is available for this disease to date. To better understand its genomics, we assembled a 211-Mb high quality genome of T. pisiformis at chromosome level with a scaffold N50 size of 20 Mbp. Totally, 12,097 protein-coding genes was predicted from the genome. Genome-level phylogenetic analysis confirmed the taxonomic affiliations with other tapeworms and revealed that T. pisiformis diverged from its closely related relative T. hydatigena ∼ 14.6 Mya. Comparative genomic analyses revealed that the T. pisiformis genome was characterized by adaptive features of strong positive selection signals from carbohydrate/lipid metabolism and body surface integrity, and of expanded gene families related to metabolism of amino acids and lipids. The high-quality genome of T. pisiformis constitutes a resource for the comparative genomics and for further applications in general parasitology.

Nitrogen is a key nutrient for all cell forms. Most organisms respond to nitrogen scarcity by slowing down their growth rate. On the contrary, our previous studies have shown that Papiliotrema laurentii strain RY1 has a robust growth under nitrogen starvation. To understand the global regulation that leads to such an extraordinary response, we undertook a de novo approach for transcriptome analysis of the yeast. Close to 33 million sequence reads of high quality for nitrogen limited and enriched condition were generated using Illumina NextSeq500. Trinity analysis and clustered transcripts annotation of the reads produced 17,611 unigenes, out of which 14,157 could be annotated. Gene Ontology term analysis generated 44.92% cellular component terms, 39.81% molecular function terms and 15.24% biological process terms. The most over represented pathways in general were translation, carbohydrate metabolism, amino acid metabolism, general metabolism, folding, sorting, degradation followed by transport and catabolism, nucleotide metabolism, replication and repair, transcription and lipid metabolism. A total of 4256 Single Sequence Repeats were identified. Differential gene expression analysis detected 996 P-significant transcripts to reveal transmembrane transport, lipid homeostasis, fatty acid catabolism and translation as the enriched terms which could be essential for Papiliotrema laurentii strain RY1 to adapt during nitrogen deprivation. Transcriptome data was validated by quantitative real-time PCR analysis of twelve transcripts. To the best of our knowledge, this is the first report of Papiliotrema laurentii strain RY1 transcriptome which would play a pivotal role in understanding the biochemistry of the yeast under acute nitrogen stress and this study would be encouraging to initiate extensive investigations into this Papiliotrema system.

The endosperm of the cereal crop is an important nutrient source for humans. It also acts as a critical integrator of plant seed growth and development. Despite its importance, the comprehensive understanding in regulating of endosperm development in rice remains elusive. Here, we performed a genomic survey comprising the identification and functional characterization of the endosperm-specific genes (OsEnS) in rice using Affymetrix microarray data and Gene Ontology (GO) analysis. A total of 151 endosperm-specific genes were identified, and the expression patterns of 13 selected genes were confirmed by qRT-PCR analysis. Promoter regions of the endosperm-specific expression genes were analyzed by PLACE Signal Scan Search. The results indicated that some motifs were involved in endosperm-specific expression regulation, and some cis-elements were responsible for hormone regulation. The bootstrap analysis indicated that the RY repeat (CATGCA box) was over-represented in promoter regions of endosperm-specific expression genes. GO analysis indicated that these genes could be classified into 12 groups, namely, transcription factor, stress/defense, seed storage protein (SSP), carbohydrate and energy metabolism, seed maturation, protein metabolism, lipid metabolism, transport, cell wall related, hormone related, signal transduction, and one unclassified group. Taken together, our results provide informative clues for further functional characterization of the endosperm-specific genes, which facilitate the understanding of the molecular mechanism in rice endosperm development.

Topping is an important agronomic practice that significantly impacts the yield of various crop plants. Topping and the regulation of axillary shoot outgrowth are common agronomic practices in tobacco. However, the effects of topping on gene expression in tobacco remain unknown. We applied the Illumina HiSeq™ 2000 platform and analyzed differentially expressed genes (DEGs) from untopped and topped plants to study the global changes in gene expression in response to topping. We found that the number of DEGs varied from 7609 to 18,770 based on the reads per kilobase per million mapped reads (RPKM) values. The Gene Ontology (GO) enrichment analysis revealed that the cellular carbohydrate metabolic process and the disaccharide metabolic process, which may contribute to starch accumulation and stress/defense, were overrepresented terms for the DEGs. Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis revealed that many DEGs were involved in starch and sucrose metabolism, glycolysis/gluconeogenesis, pyruvate metabolism, and plant hormone signal transduction, among other processes. The knowledge gained will improve our understanding of the processes of axillary shoot formation and enlargement at the transcriptional level. This study lays a solid foundation for future studies on molecular mechanisms underlying the growth of axillary shoots.

Animals recognize the availability of nutrients and regulate the intake and storage of these nutrients accordingly. However, the molecular mechanisms underlying nutrient sensing and subsequent changes in behavior and metabolism are not fully understood. Mlx interactor (Mio), the Drosophila homolog of carbohydrate response element binding protein (ChREBP), functions as a transcription factor in the fat body of the fly to control triglyceride storage as well as feeding, suggesting that Mio may act in a nutrient-sensing pathway to coordinate food consumption and metabolism. Here, we show that Mio functions in neurons in Drosophila to regulate feeding and nutrient storage. Pan-neuronal disruption of Mio function leads to increased triglyceride and glycogen storage, and this phenotype is not due to increased food consumption. Interestingly, targeted disruption of Mio specifically in the insulin-producing cells (IPCs) has little effect on nutrient storage, but increases food consumption suggesting that Mio acts in these neurons to control feeding behavior. Since Mio is a transcription factor, one possible way Mio may act in the IPCs to control feeding is through regulating the expression of Drosophila insulin-like peptides (dilps) or drosulfakinin (dsk), neuropeptides produced in the IPCs. Consistent with this hypothesis, IPC-specific knockdown of Mio leads to an increase in dilp3 expression, while not affecting dilp2, 5 or dsk levels. Together, this study indicates a new function for Mio in the Drosophila brain and specifically in the IPCs, controlling neuropeptide gene expression, feeding and metabolism in accordance with nutrient availability.

Drought is one of the most severe threats to the growth, development and yield of plant. In order to unravel the molecular basis underlying the high tolerance of pitaya (Hylocereus undatus) to drought stress, suppression subtractive hybridization (SSH) and cDNA microarray approaches were firstly combined to identify the potential important or novel genes involved in the plant responses to drought stress. The forward (drought over drought-free) and reverse (drought-free over drought) suppression subtractive cDNA libraries were constructed using in vitro shoots of cultivar 'Zihonglong' exposed to drought stress and drought-free (control). A total of 2112 clones, among which half were from either forward or reverse SSH library, were randomly picked up to construct a pitaya cDNA microarray. Microarray analysis was carried out to verify the expression fluctuations of this set of clones upon drought treatment compared with the controls. A total of 309 expressed sequence tags (ESTs), 153 from forward library and 156 from reverse library, were obtained, and 138 unique ESTs were identified after sequencing by clustering and blast analyses, which included genes that had been previously reported as responsive to water stress as well as some functionally unknown genes. Thirty six genes were mapped to 47 KEGG pathways, including carbohydrate metabolism, lipid metabolism, energy metabolism, nucleotide metabolism, and amino acid metabolism of pitaya. Expression analysis of the selected ESTs by reverse transcriptase polymerase chain reaction (RT-PCR) corroborated the results of differential screening. Moreover, time-course expression patterns of these selected ESTs further confirmed that they were closely responsive to drought treatment. Among the differentially expressed genes (DEGs), many are related to stress tolerances including drought tolerance. Thereby, the mechanism of drought tolerance of this pitaya genotype is a very complex physiological and biochemical process, in which multiple metabolism pathways and many genes were implicated. The data gained herein provide an insight into the mechanism underlying the drought stress tolerance of pitaya, as well as may facilitate the screening of candidate genes for drought tolerance.

We present a novel genome-scale metabolic model iWV1213 of Mucor circinelloides, which is an oleaginous fungus for industrial applications. The model contains 1213 genes, 1413 metabolites and 1326 metabolic reactions across different compartments. We demonstrate that iWV1213 is able to accurately predict the growth rates of M. circinelloides on various nutrient sources and culture conditions using Flux Balance Analysis and Phenotypic Phase Plane analysis. Comparative analysis of three oleaginous genome-scale models, including M. circinelloides (iWV1213), Mortierella alpina (iCY1106) and Yarrowia lipolytica (iYL619_PCP) revealed that iWV1213 possesses a higher number of genes involved in carbohydrate, amino acid, and lipid metabolisms that might contribute to its versatility in nutrient utilization. Moreover, the identification of unique and common active reactions among the Zygomycetes oleaginous models using Flux Variability Analysis unveiled a set of gene/enzyme candidates as metabolic engineering targets for cellular improvement. Thus, iWV1213 offers a powerful metabolic engineering tool for multi-level omics analysis, enabling strain optimization as a cell factory platform of lipid-based production.