Fine-tuning of p-coumaric acid synthesis to increase (2S)-naringenin production in yeast.

(2S)-Naringenin is a key precursor for biosynthesis of various high-value flavonoids and possesses a variety of nutritional and pharmaceutical properties on human health. Systematic optimization approaches have been employed to improve (2S)-naringenin production in different microbial hosts. However, very few studies have focused on the spatiotemporal distribution of (2S)-naringenin and the related pathway intermediate p-coumaric acid, which is an important factor for efficient production. Here, we first optimized the (2S)-naringenin biosynthetic pathway by alleviating the bottleneck downstream of p-coumaric acid and increasing malonyl-CoA supply, which improved (2S)-naringenin production but significant accumulation of p-coumaric acid still existed extracellularly. We thus established a dual dynamic control system through combining a malonyl-CoA biosensor regulator and an RNAi strategy, to autonomously control the synthesis of p-coumaric acid with the supply of malonyl-CoA. Furthermore, screening potential transporters led to identification of Pdr12 for improved (2S)-naringenin production and reduced accumulation of p-coumaric acid. Finally, a titer of 2.05 g/L (2S)-naringenin with negligible accumulation of p-coumaric acid was achieved in a fed batch fermentation. Our work highlights the importance of systematic control of pathway intermediates for efficient microbial production of plant natural products.

Riboswitch-guided chalcone synthase engineering and metabolic flux optimization for enhanced production of flavonoids.

Flavonoids are a group of secondary metabolites from plants that have received attention as high value-added pharmacological substances. Recently, a robust and efficient bioprocess using recombinant microbes has emerged as a promising approach to supply flavonoids. In the flavonoid biosynthetic pathway, the rate of chalcone synthesis, the first committed step, is a major bottleneck. However, chalcone synthase (CHS) engineering was difficult because of high-level conservation and the absence of effective screening tools, which are limited to overexpression or homolog-based combinatorial strategies. Furthermore, it is necessary to precisely regulate the metabolic flux for the optimum availability of malonyl-CoA, a substrate of chalcone synthesis. In this study, we engineered CHS and optimized malonyl-CoA availability to establish a platform strain for naringenin production, a key molecular scaffold for various flavonoids. First, we engineered CHS through synthetic riboswitch-based high-throughput screening of rationally designed mutant libraries. Consequently, the catalytic efficiency (kcat/Km) of the optimized CHS enzyme was 62% higher than that of the wild-type enzyme. In addition to CHS engineering, we designed genetic circuits using transcriptional repressors to fine-tune the malonyl-CoA availability. The best mutant with synergistic effects of the engineered CHS and the optimized genetic circuit produced 98.71 mg/L naringenin (12.57 mg naringenin/g glycerol), which is the highest naringenin concentration and yield from glycerol in similar culture conditions reported to date, a 2.5-fold increase compared to the parental strain. Overall, this study provides an effective strategy for efficient production of flavonoids.

Metabolic engineering of the malonyl-CoA pathway to efficiently produce malonate in Saccharomyces cerevisiae.

Malonate is a platform chemical that has been utilized to synthesize many valuable chemical compounds. Here, Saccharomyces cerevisiae was metabolically engineered to produce malonate through the malonyl-CoA pathway. To construct the key step of converting malonyl-CoA to malonate, a native mitochondrial 3-hydroxyisobutyryl-CoA hydrolase gene EHD3 was mutated to target the cytoplasm and obtain malonyl-CoA hydrolase activity. The malonyl-CoA hydrolase activity of Ehd3 was achieved by mutating the malonyl-CoA binding site F121 to I121 and the active site E124 to seven amino acids (S/T/H/K/R/N/Q). We identified that the strain with E124S mutation had the highest malonate titer with 13.6 mg/L. Genomic integration of the mutant EHD3 and ACC1** to delta sequence sites was further explored to increase their reliable expression. Accordingly, a screening method with the work flow of fluorescence detection, shake-tube fermentation, and shake-flask fermentation was constructed to screen high copy delta sequences efficiently. The malonate titer was improved to 73.55 mg/L after screening the ∼1500 integrative strains, which was increased 4.4-folds than that of the episomal strain. We further engineered the strain by regulating the expression of key enzyme in the malonyl-CoA pathway to improve the precursor supply and inhibiting its competing pathways, and the final engineered strain LMA-16 produced 187.25 mg/L in the flask, 14-fold compared with the initial episomal expression strain. Finally, the combined efforts increased the malonate titer to 1.62 g/L in fed-batch fermentation.

Heterogenous Clinical Landscape in a Consanguineous Malonic Aciduria Family.

Malonic aciduria is an extremely rare inborn error of metabolism due to malonyl-CoA decarboxylase deficiency. This enzyme is encoded by the MLYCD (Malonyl-CoA Decarboxylase) gene, and the disease has an autosomal recessive inheritance. Malonic aciduria is characterized by systemic clinical involvement, including neurologic and digestive symptoms, metabolic acidosis, hypoglycemia, failure to thrive, seizures, developmental delay, and cardiomyopathy. We describe here two index cases belonging to the same family that, despite an identical genotype, present very different clinical pictures. The first case is a boy with neonatal metabolic symptoms, abnormal brain MRI, and dilated cardiomyopathy. The second case, the cousin of the first patient in a consanguineous family, showed later symptoms, mainly with developmental delay. Both patients showed high levels of malonylcarnitine on acylcarnitine profiles and malonic acid on urinary organic acid chromatographies. The same homozygous pathogenic variant was identified, c.346C > T; p. (Gln116*). We also provide a comprehensive literature review of reported cases. A review of the literature yielded 52 cases described since 1984. The most common signs were developmental delay and cardiomyopathy. Increased levels of malonic acid and malonylcarnitine were constant. Presentations ranged from neonatal death to patients surviving past adolescence. These two cases and reported patients in the literature highlight the inter- and intrafamilial variability of malonic aciduria.

Considerations of expanded carrier screening: Lessons learned from combined malonic and methylmalonic aciduria.

BACKGROUND: Expanded carrier screening (ECS) utilizes high-throughput next-generation sequencing to evaluate an individual's carrier status for multiple conditions. Combined malonic and methylmalonic aciduria (CMAMMA) due to ACSF3 deficiency is a rare inherited disease included in such screening panels. Some cases have been reported with metabolic symptoms in childhood yet other cases describe a benign clinical course, suggesting the clinical phenotype is not well defined. METHODS/CASE REPORT: Clinical and laboratory findings during the prenatal period were obtained retrospectively from medical records. RESULTS: A 37-year-old nulliparous woman and her partner were each identified as carriers of ACSF3 variants and presented at 9 weeks gestation for prenatal genetic consultation. The couple received extensive genetic counseling and proceeded with chorionic villus sampling at 11 weeks gestation. Subsequent analysis confirmed that the fetus inherited both parental ACSF variants. The couple was devastated by the results and after reviewing options of pregnancy continuation and termination, they decided to terminate the pregnancy. Following this decision, the patient was diagnosed with acute stress disorder. CONCLUSION: This case highlights how expanded carrier screening adds complexity to reproductive decision-making. Stronger guidelines and additional research are needed to direct and evaluate the timing, composition, and implementation of ECS panels.

Isoflavone malonyl-CoA acyltransferase GmMaT2 is involved in nodulation of soybean by modifying synthesis and secretion of isoflavones.

Malonyl-CoA:flavonoid acyltransferases (MaTs) modify isoflavones, but only a few have been characterized for activity and assigned to specific physiological processes. Legume roots exude isoflavone malonates into the rhizosphere, where they are hydrolyzed into isoflavone aglycones. Soybean GmMaT2 was highly expressed in seeds, root hairs, and nodules. GmMaT2 and GmMaT4 recombinant enzymes used isoflavone 7-O-glucosides as acceptors and malonyl-CoA as an acyl donor to generate isoflavone glucoside malonates. GmMaT2 had higher activity towards isoflavone glucosides than GmMaT4. Overexpression in hairy roots of GmMaT2 and GmMaT4 produced more malonyldaidzin, malonylgenistin, and malonylglycitin, and resulted in more nodules than control. However, only GmMaT2 knockdown (KD) hairy roots showed reduced levels of malonyldaidzin, malonylgenistin, and malonylglycitin, and, likewise, reduced nodule numbers. These were consistent with the up-regulation of only GmMaT2 by rhizobial infection, and higher expression levels of early nodulation genes in GmMaT2- and GmMaT4-overexpressing roots, but lower only in GmMaT2-KD roots compared with control roots. Higher malonyl isoflavonoid levels in transgenic hairy roots were associated with higher levels of isoflavones in root exudates and more nodules, and vice versa. We suggest that GmMaT2 participates in soybean nodulation by catalyzing isoflavone malonylation and affecting malonyl isoflavone secretion for activation of Nod factor and nodulation.

Sensing of nutrients by CPT1C controls SAC1 activity to regulate AMPA receptor trafficking.

Carnitine palmitoyltransferase 1C (CPT1C) is a sensor of malonyl-CoA and is located in the ER of neurons. AMPA receptors (AMPARs) mediate fast excitatory neurotransmission in the brain and play a key role in synaptic plasticity. In the present study, we demonstrate across different metabolic stress conditions that modulate malonyl-CoA levels in cortical neurons that CPT1C regulates the trafficking of the major AMPAR subunit, GluA1, through the phosphatidyl-inositol-4-phosphate (PI(4)P) phosphatase SAC1. In normal conditions, CPT1C down-regulates SAC1 catalytic activity, allowing efficient GluA1 trafficking to the plasma membrane. However, under low malonyl-CoA levels, such as during glucose depletion, CPT1C-dependent inhibition of SAC1 is released, facilitating SAC1's translocation to ER-TGN contact sites to decrease TGN PI(4)P pools and trigger GluA1 retention at the TGN. Results reveal that GluA1 trafficking is regulated by CPT1C sensing of malonyl-CoA and provide the first report of a SAC1 inhibitor. Moreover, they shed light on how nutrients can affect synaptic function and cognition.

Gap junctions deliver malonyl-CoA from soma to germline to support embryogenesis in Caenorhabditis elegans.

Gap junctions are ubiquitous in metazoans and play critical roles in important biological processes, including electrical conduction and development. Yet, only a few defined molecules passing through gap junction channels have been linked to specific functions. We isolated gap junction channel mutants that reduce coupling between the soma and germ cells in the Caenorhabditis elegans gonad. We provide evidence that malonyl-CoA, the rate-limiting substrate for fatty acid synthesis (FAS), is produced in the soma and delivered through gap junctions to the germline; there it is used in fatty acid synthesis to critically support embryonic development. Separation of malonyl-CoA production from its site of utilization facilitates somatic control of germline development. Additionally, we demonstrate that loss of malonyl-CoA production in the intestine negatively impacts germline development independently of FAS. Our results suggest that metabolic outsourcing of malonyl-CoA may be a strategy by which the soma communicates nutritional status to the germline.

A Korean child diagnosed with malonic aciduria harboring a novel start codon mutation following presentation with dilated cardiomyopathy.

BACKGROUND: Malonic aciduria (MA, OMIM#248360) is an extremely rare inherited metabolic disorder caused by the deficiency of malonyl-CoA decarboxylase. The phenotype exhibited by patients with MA is variable, but may include symptoms, such as developmental delay in early childhood, seizures, vomiting, metabolic acidosis, hypoglycemia, ketosis, and cardiomyopathy. We describe the first case of a Korean child with MA who presented with dilated cardiomyopathy (DCMP) at the age of 3 months. METHODS AND RESULTS: A 3-month-old Korean boy visited our hospital for diagnosis and management of cardiomegaly. Newborn screening for inherited metabolic diseases showed a normal result; therefore, DCMP management was initiated. Biochemical and the MLYCD gene analyses subsequently confirmed diagnosis of MA. Elevated plasma C3DC level and excessive excretion of urinary malonate were observed, and two pathogenic variants, including a novel start codon mutation (c.1A>G), were identified in MLYCD. A low long-chain fat diet with middle-chain triglyceride formula and L-carnitine supplementation was initiated. The patient is now 5 years old and exhibits considerably improved cardiac function. CONCLUSIONS: MA can be diagnosed using newborn screening; however, negative results do not exclude the possibility of disease. Metabolic screening for differential diagnosis of infantile DCMP is recommended to rule out rare, but manageable, metabolic cardiomyopathies.

Coupling metabolic addiction with negative autoregulation to improve strain stability and pathway yield.

Metabolic addiction, an organism that is metabolically addicted with a compound to maintain its growth fitness, is an underexplored area in metabolic engineering. Microbes with heavily engineered pathways or genetic circuits tend to experience metabolic burden leading to degenerated or abortive production phenotype during long-term cultivation or scale-up. A promising solution to combat metabolic instability is to tie up the end-product with an intermediary metabolite that is essential to the growth of the producing host. Here we present a simple strategy to improve both metabolic stability and pathway yield by coupling chemical addiction with negative autoregulatory genetic circuits. Naringenin and lipids compete for the same precursor malonyl-CoA with inversed pathway yield in oleaginous yeast. Negative autoregulation of the lipogenic pathways, enabled by CRISPRi and fatty acid-inducible promoters, repartitions malonyl-CoA to favor flavonoid synthesis and increased naringenin production by 74.8%. With flavonoid-sensing transcriptional activator FdeR and yeast hybrid promoters to control leucine synthesis and cell grwoth fitness, this amino acid feedforward metabolic circuit confers a flavonoid addiction phenotype that selectively enrich the naringenin-producing pupulation in the leucine auxotrophic yeast. The engineered yeast persisted 90.9% of naringenin titer up to 324 generations. Cells without flavonoid addiction regained growth fitness but lost 94.5% of the naringenin titer after cell passage beyond 300 generations. Metabolic addiction and negative autoregulation may be generalized as basic tools to eliminate metabolic heterogeneity, improve strain stability and pathway yield in long-term and large-scale bioproduction.

A Synthetic Malonyl-CoA Metabolic Oscillator in Komagataella phaffii.

Malonyl-CoA is a key metabolic molecule that participates in a diverse range of physiological responses and can act as a building block for a variety of value-added pharmaceuticals and chemicals. The cytosolic malonyl-CoA concentration is usually very low, and thus dynamic metabolic control of malonyl-CoA variation will aid its stable formation and efficient consumption. We developed a synthetic malonyl-CoA metabolic oscillator in yeast. A synthetic regulatory protein, Prm1-FapR, was constructed by fusing a yeast transcriptional activator, Prm1, with a bacterial malonyl-CoA-sensitive transcription repressor, FapR. Two oppositely regulated biosensors were then engineered. A total of 18 hybrid promoter variants were designed, each carrying the operator sequence (fapO) of FapR and the core promoter of PAOX1 (cPAOX1), which is naturally regulated by Prm1. The promoter activities of these variants, regulated by Prm1-FapR, were tested. Through this process, a sensor for Prm1-FapR/(-52)fapO-PAOX1 carrying an activation/deactivation regulation module was built. Meanwhile, 24 promoter variants of PGAP with fapO inserted were designed and tested using the fusion regulator, giving a sensor for Prm1-FapR/PGAP-(+22) fapO that contained a repression/derepression regulation module. Both sensors were subsequently integrated into a single cell, which exhibited correct metabolic switching of eGFP and mCherry reporters following manipulation of cytosolic malonyl-CoA levels. The Prm1-FapR/(-52)fapO-PAOX1 and the Prm1-FapR/PGAP-(+22)fapO were also used to control the malonyl-CoA source and sink pathways, respectively, for the synthesis of 6-methylsalicylic acid. This finally led to an oscillatory metabolic mode of cytosolic malonyl-CoA. Such a metabolator is useful in exploring potential industrial and biomedical applications not limited by natural cellular behavior.

Coenzyme A and Its Thioester Pools in Obese Zucker and Zucker Diabetic Fatty Rats.

Feeding behavior is closely related to hypothalamic malonyl-CoA level in the brain and diet-induced obesity affects total CoA pools in liver. Herein, we performed a comprehensive analysis of the CoA pools formed in thirteen tissues of Zucker and Zucker diabetic fatty (ZDF) rats. Hypothalamic malonyl-CoA levels in obese rats remained low and were almost the same as those of lean rats, despite obese rats having much higher content of leptin, insulin, and glucose in their sera. Regardless of the fa-genotypes, larger total CoA pools were formed in the livers of ZDF rats and the size of hepatic total CoA pools in Zucker rats showed almost one tenth of the size of ZDF rats. The decreased total CoA pool sizes in Zucker rats was observed in the brown adipose tissues, while ZDF-fatty rats possessed 6% of total CoA pool in the lean rats in response to fa deficiency. This substantially lower CoA content in the obese rats would be disadvantageous to non-shivering thermogenesis. Thus, comparing the intracellular CoA behaviors between Zucker and ZDF rats, as well as the lean and fatty rats of each strain would help to elucidate features of obesity and type 2 diabetes in combination with result (s) of differential gene expression analysis and/or comparative genomics.

Optimizing Oleaginous Yeast Cell Factories for Flavonoids and Hydroxylated Flavonoids Biosynthesis.

Plants possess myriads of secondary metabolites with a broad spectrum of health-promoting benefits. To date, plant extraction is still the primary route to produce high-value natural products which inherently suffers from economics and scalability issues. Heterologous expression of plant biosynthetic gene clusters in microbial host is considered as a feasible approach to overcoming these limitations. Oleaginous yeast produces a large amount of lipid bodies, the abundant membrane structure and the lipophilic environment provide the ideal environment for the regioselectivity and stereoselectivity of many plant-derived P450 enzymes. In this work, we used modular method to construct, characterize, and optimize the flavonoid pathways in Yarrowia lipolytica. We also evaluated various precursor biosynthetic routes and unleashed the metabolic potential of Y. lipolytica to produce flavonoids and hydroxylated flavonoids. Specifically, we have identified that chalcone synthase (CHS) and cytochrome P450 reductases (CPR) were the bottlenecks of hydroxylated flavonoid production. We determined the optimal gene copy number of CHS and CPR to be 5 and 2, respectively. We further removed precursor pathway limitations by expressing genes associated with chorismate and malonyl-CoA supply. With pH and carbon-nitrogen ratio (C/N) optimization, our engineered strain produced 252.4 mg/L naringenin, 134.2 mg/L eriodictyol, and 110.5 mg/L taxifolin from glucose in shake flasks. Flavonoid and its hydroxylated derivatives are most prominently known as antioxidant and antiaging agents. These findings demonstrate our ability to harness the oleaginous yeast as the microbial workhorse to expand nature's biosynthetic potential, enabling us to bridge the gap between drug discovery and natural product manufacturing.

Expanding the Dynamic Range of a Transcription Factor-Based Biosensor in Saccharomyces cerevisiae.

Metabolite biosensors are useful tools for high-throughput screening approaches and pathway regulation approaches. An important feature of biosensors is the dynamic range. To expand the maximum dynamic range of a transcription factor-based biosensor in Saccharomyces cerevisiae, using the fapO/FapR system from Bacillus subtilis as an example case, five native promoters, including constitutive and glucose-regulated ones, were modified. By evaluating different binding site (BS) positions in the core promoters, we identified locations that resulted in a high maximum dynamic range with low expression under repressed conditions. We further identified BS positions in the upstream element region of the TEF1 promoter that did not influence the native promoter strength but resulted in repression in the presence of a chimeric repressor consisting of FapR and the yeast repressor Mig1. These modified promoters with broad dynamic ranges will provide useful information for the engineering of future biosensors and their use in complex genetic circuits.

Normal FGF-21-Serum Levels in Patients with Carnitine Palmitoyltransferase II (CPT II) Deficiency.

Fibroblast growth factor 21 (FGF-21) is known to be a biomarker for mitochondrial disorders. An upregulation of FGF-21 in serum and muscle of carnitine palmitoyltransferase I (CPT I) and carnitine palmitoyltransferase II (CPT II) knock-out mice has been reported. In human CPT II deficiency, enzyme activity and protein content are normal, but the enzyme is abnormally regulated by malonyl-CoA and is abnormally thermolabile. Citrate synthase (CS) activity is increased in patients with CPT II deficiency. This may indicate a compensatory response to an impaired function of CPT II. In this study, FGF-21 serum levels in patients with CPT II deficiency during attack free intervals and in healthy controls were measured by enzyme linked immunosorbent assay (ELISA). The data showed no significant difference between FGF-21 concentration in the serum of patients with CPT II deficiency and that in the healthy controls. The results of the present work support the hypothesis that in muscle CPT II deficiency, in contrast to the mouse knockout model, mitochondrial fatty acid utilization is not persistently reduced. Thus, FGF-21 does not seem to be a useful biomarker in the diagnosis of CPT II deficiency.

Modulation of the central carbon metabolism of Corynebacterium glutamicum improves malonyl-CoA availability and increases plant polyphenol synthesis.

In recent years microorganisms have been engineered towards synthesizing interesting plant polyphenols such as flavonoids and stilbenes from glucose. Currently, the low endogenous supply of malonyl-CoA, indispensable for plant polyphenol synthesis, impedes high product titers. Usually, limited malonyl-CoA availability during plant polyphenol production is avoided by supplementing fatty acid synthesis-inhibiting antibiotics such as cerulenin, which are known to increase the intracellular malonyl-CoA pool as a side effect. Motivated by the goal of microbial polyphenol synthesis being independent of such expensive additives, we used rational metabolic engineering approaches to modulate regulation of fatty acid synthesis and flux into the tricarboxylic acid cycle (TCA cycle) in Corynebacterium glutamicum strains capable of flavonoid and stilbene synthesis. Initial experiments showed that sole overexpression of genes coding for the native malonyl-CoA-forming acetyl-CoA carboxylase is not sufficient for increasing polyphenol production in C. glutamicum. Hence, the intracellular acetyl-CoA availability was also increased by reducing the flux into the TCA cycle through reduction of citrate synthase activity. In defined cultivation medium, the constructed C. glutamicum strains accumulated 24 mg·L -1 (0.088 mM) naringenin or 112 mg·L -1 (0.49 mM) resveratrol from glucose without supplementation of phenylpropanoid precursor molecules or any inhibitors of fatty acid synthesis.

Role of the malonyl-CoA synthetase ACSF3 in mitochondrial metabolism.

Malonyl-CoA is a central metabolite in fatty acid biochemistry. It is the rate-determining intermediate in fatty acid synthesis but is also an allosteric inhibitor of the rate-setting step in mitochondrial long-chain fatty acid oxidation. While these canonical cytoplasmic roles of malonyl-CoA have been well described, malonyl-CoA can also be generated within the mitochondrial matrix by an alternative pathway: the ATP-dependent ligation of malonate to Coenzyme A by the malonyl-CoA synthetase ACSF3. Malonate, a competitive inhibitor of succinate dehydrogenase of the TCA cycle, is a potent inhibitor of mitochondrial respiration. A major role for ACSF3 is to provide a metabolic pathway for the clearance of malonate by the generation of malonyl-CoA, which can then be decarboxylated to acetyl-CoA by malonyl-CoA decarboxylase. Additionally, ACSF3-derived malonyl-CoA can be used to malonylate lysine residues on proteins within the matrix of mitochondria, possibly adding another regulatory layer to post-translational control of mitochondrial metabolism. The discovery of ACSF3-mediated generation of malonyl-CoA defines a new mitochondrial metabolic pathway and raises new questions about how the metabolic fates of this multifunctional metabolite intersect with mitochondrial metabolism.

Combinatorial Optimization of Resveratrol Production in Engineered E. coli.

Resveratrol, a plant-derived polyphenolic compound with various health activities, is widely used in nutraceutical and food additives. Herein, combinatorial optimization of resveratrol biosynthetic pathway and intracellular environment of E. coli was carried out. By screening pathway genes from various species and exploring their expression pattern, we initially constructed resveratrol-producing strains. Further targeting at availability of malonyl-CoA through expressing ACC of Corynebacterium glutamicum and antisense inhibiting native fabD significantly increased resveratrol biosynthesis. Transport engineering for resveratrol secretion and molecular chaperones helping for folding heterologous enzymes were employed to improve the intracellular environments in remarkable degrees. By introducing PcTAL of Phanerochaete chrysosporium and tuning expression model of PcTAL, At4CL, and VvSTS, an engineered E. coli produced 57.77 mg/L of resveratrol from l-tyrosine. After integrating the above strategies, resveratrol titer reached to 238.71 mg/L from l-tyrosine. The combinatorial optimization in this study provides a promising strategy to produce valuable natural products in heterologous expression systems.

Full-Length Transcriptome Sequencing and Modular Organization Analysis of the Naringin/Neoeriocitrin-Related Gene Expression Pattern in Drynaria roosii.

Drynaria roosii (Nakaike) is a traditional Chinese medicinal fern, known as 'GuSuiBu'. The effective components, naringin and neoeriocitrin, share a highly similar chemical structure and medicinal function. Our HPLC-tandem mass spectrometry (MS/MS) results showed that the accumulation of naringin/neoeriocitrin depended on specific tissues or ages. However, little was known about the expression patterns of naringin/neoeriocitrin-related genes involved in their regulatory pathways. Due to a lack of basic genetic information, we applied a combination of single molecule real-time (SMRT) sequencing and second-generation sequencing (SGS) to generate the complete and full-length transcriptome of D. roosii. According to the SGS data, the differentially expressed gene (DEG)-based heat map analysis revealed that naringin/neoeriocitrin-related gene expression exhibited obvious tissue- and time-specific transcriptomic differences. Using the systems biology method of modular organization analysis, we clustered 16,472 DEGs into 17 gene modules and studied the relationships between modules and tissue/time point samples, as well as modules and naringin/neoeriocitrin contents. We found that naringin/neoeriocitrin-related DEGs distributed in nine distinct modules, and DEGs in these modules showed significantly different patterns of transcript abundance to be linked to specific tissues or ages. Moreover, weighted gene co-expression network analysis (WGCNA) results further identified that PAL, 4CL and C4H, and C3H and HCT acted as the major hub genes involved in naringin and neoeriocitrin synthesis, respectively, and exhibited high co-expression with MYB- and basic helix-leucine-helix (bHLH)-regulated genes. In this work, modular organization and co-expression networks elucidated the tissue and time specificity of the gene expression pattern, as well as hub genes associated with naringin/neoeriocitrin synthesis in D. roosii. Simultaneously, the comprehensive transcriptome data set provided important genetic information for further research on D. roosii.