Co-regulatory network analysis of the main secondary metabolite (SM) biosynthesis in Crocus sativus L.

Saffron (Crocus sativus L.) is being embraced as the most important medicinal plant and the commercial source of saffron spice. Despite the beneficial economic and medicinal properties of saffron, the regulatory mechanism of the correlation of TFs and genes related to the biosynthesis of the apocarotenoids pathway is less obvious. Realizing these regulatory hierarchies of gene expression networks related to secondary metabolites production events is the main challenge owing to the complex and extensive interactions between the genetic behaviors. Recently, high throughput expression data have been highly feasible for constructing co-regulation networks to reveal the regulated processes and identifying novel candidate hub genes in response to complex processes of the biosynthesis of secondary metabolites. Herein, we performed Weighted Gene Co-expression Network Analysis (WGCNA), a systems biology method, to identify 11 regulated modules and hub TFs related to secondary metabolites. Three specialized modules were found in the apocarotenoids pathway. Several hub TFs were identified in notable modules, including MADS, C2H2, ERF, bZIP, HD-ZIP, and zinc finger protein MYB and HB, which were potentially associated with apocarotenoid biosynthesis. Furthermore, the expression levels of six hub TFs and six co-regulated genes of apocarotenoids were validated with RT-qPCR. The results confirmed that hub TFs specially MADS, C2H2, and ERF had a high correlation (P < 0.05) and a positive effect on genes under their control in apocarotenoid biosynthesis (CCD2, GLT2, and ADH) among different C. sativus ecotypes in which the metabolite contents were assayed. Promoter analysis of the co-expressed genes of the modules involved in apocarotenoids biosynthesis pathway suggested that not only are the genes co-expressed, but also share common regulatory motifs specially related to hub TFs of each module and that they may describe their common regulation. The result can be used to engineer valuable secondary metabolites of C. sativus by manipulating the hub regulatory TFs.

Characterizing reference genes for high-fidelity gene expression analysis under different abiotic stresses and elicitor treatments in fenugreek leaves.

BACKGROUND: Quantifying gene expression is a critical aspect of applied genomics research across all organisms, and real-time PCR has emerged as a powerful tool for this purpose. However, selecting appropriate internal control genes for data normalization presents specific challenges. This study aimed to identify suitable reference genes for gene expression analysis under various conditions, encompassing salinity, low and high-temperature stresses, and different elicitor treatments. These treatments included titanium dioxide, cold plasma, 24-epibrassinolide, and melatonin, resulting in a total of 13 unique treatments and 148 treatment combinations applied to fenugreek plants. RESULTS: As per the analysis performed with the BestKeeper tool, EEF-1α, and GAPDH were recognized as the most stable reference genes under the majority of conditions. Furthermore, the GeNorm and NormFinder tools identified ß-tubulin and EEF-1α as the most stable reference genes. The findings of this research demonstrated that, although the stability of three reference genes expression was acceptable in almost all evaluated treatments, fluctuations in their expression were observed under the treatments of cold stress with TiO2 NPs application, cold plasma application with salinity stress, and cold plasma application with high-temperature stress compared to others. Simultaneously, the GeNorm analysis results demonstrated that in the mentioned treatments, relying on only one reference gene is inadequate. To corroborate the results, we examined the expression profile of the SSR gene, a pivotal gene in diosgenin biosynthesis, under all investigated treatments and treatment combinations. The outcomes suggested that employing stable reference genes yielded highly consistent results. CONCLUSIONS: The varying expression patterns of the target genes emphasize the crucial need for precise optimization of experimental conditions and selecting stable reference genes to achieve accurate results in gene expression studies utilizing real-time PCR. These findings offer valuable insights into the selection of appropriate reference genes for gene expression analysis under diverse conditions using real-time PCR.

Exogenous 24-epibrassinolide ameliorates tolerance to high-temperature by adjusting the biosynthesis of pigments, enzymatic, non-enzymatic antioxidants, and diosgenin content in fenugreek.

High-temperature stress is widely considered a main plant-growth-limiting factor. The positive effects of 24-epibrassinolide (EBR) as analogs of brassinosteroids (BRs) in modulating abiotic stresses have led this hormone to be referred to as a growth regulator in plants. The current study highlights the influence of EBR on enhancing tolerance to high-temperature and altering the diosgenin content in fenugreek. Different amounts of EBR (4, 8, and 16 µM), harvesting times (6, and 24 h), as well as temperature regimes (23 °C, and 42 °C) were, used as treatments. EBR application under normal temperature and high-temperature stress resulted in decreased malondialdehyde content and electrolyte leakage percentage, while the activity of antioxidant enzymes improved significantly. Exogenous EBR application possibly contributes to activating the nitric oxide, H2O2, and ABA-dependent pathways, enhancing the biosynthesis of abscisic acid and auxin, and regulating the signal transduction pathways, which raises fenugreek tolerance to high-temperature. The SQS (eightfold), SEP (2.8-fold), CAS (11-fold), SMT (17-fold), and SQS (sixfold) expression, considerably increased following EBR application (8 µM) compared to the control. Compared to the control, when the short-term (6 h) high-temperature stress was accompanied by EBR (8 µM), a sixfold increase in diosgenin content was achieved. Our findings highlight the potential role of exogenous 24-epibrassinolide in mitigating the high-temperature stress in fenugreek by stimulating the biosynthesis processes of enzymatic and non-enzymatic antioxidants, chlorophylls, and diosgenin. In conclusion, the current results could be of utmost importance in breeding or biotechnology-based programs of fenugreek and also in the researches related to the engineering of the biosynthesis pathway of diosgenin in this valuable plant.

Effect of redroot pigweed interference on antioxidant enzyme and light response of common bean (Phaseolus vulgaris L.) depends on cultivars and growth stages.

Redroot Pigweed (Amaranthus retroflexus L.) is an important weed that is highly competitive with common bean. Photosynthetic pigments, the activity of antioxidant enzymes, the relative expression of a number of antioxidant enzyme and light response genes, were studied in three of common bean cultivars and in V4 and R7 stages under Redroot Pigweed free and infested. The presence of weeds reduced the content of chlorophyll, relative chlorophyll and anthocyanin of common bean leaves. With the increase of weed competition, the expression of antioxidant genes and enzymes increased, which indicates the increase of their activity in order to reduce the amount of reactive oxygen species. Among the studied antioxidant enzymes, the activity of catalase and ascorbate peroxidase produced in the leaves was higher than that of superoxide dismutase. With the increase of weed interference, the expression of phytochrome interacting factor 3 (PIF3) gene as a positive regulator of light signals is increased and the expression of phytochrome rapidly regulated1 (PAR1) gene as a negative regulator is decreased. Chlorophyll a/b-binding protein (CAB1) and auxin-responsive protein IAA8 (IAA8) genes also down-regulated with increasing competition. Along with the decrease of CAB expression in the conditions of competition with weeds, the chlorophyll a, b content also decreased. Correlation between gene expression and physiological traits related to them highlights the prominent role of CWCP in maintaining yield potential.

Insights Into The Effects of Amino Acid Substitutions on The Stability of Reteplase Structure: A Molecular Dynamics Simulation Study.

Background: Reteplase (recombinant plasminogen activator, r-PA) is a recombinant protein designed to imitate the endogenous tissue plasminogen activator and catalyze the plasmin production. It is known that the application of reteplase is limited by the complex production processes and protein's stability challenges. Computational redesign of proteins has gained momentum in recent years, particularly as a powerful tool for improving protein stability and consequently its production efficiency. Hence, in the current study, we implemented computational approaches to improve r-PA conformational stability, which fairly correlates with protein's resistance to proteolysis. Objectives: The current study was developed in order to evaluate the effect of amino acid substitutions on the stability of reteplase structure using molecular dynamic simulations and computational predictions. Materials and Methods: Several web servers designed for mutation analysis were utilized to select appropriate mutations. Additionally, the experimentally reported mutation, R103S, converting wild type r-PA into non-cleavable form, was also employed. Firstly, mutant collection, consisting of 15 structures, was constructed based on the combinations of four designated mutations. Then, 3D structures were generated using MODELLER. Finally, 17 independent 20-ns molecular dynamics (MD) simulations were conducted and different analysis were performed like root-mean-square deviation (RMSD), root-mean-square fluctuations (RMSF), secondary structure analysis, number of hydrogen bonds, principal components analysis (PCA), eigenvector projection, and density analysis. Results: Predicted mutations successfully compensated the more flexible conformation caused by R103S substitution, so, improved conformational stability was analyzed from MD simulations. In particular, R103S/A286I/G322I indicated the best results and remarkably enhanced the protein stability. Conclusion: The conformational stability conferred by these mutations will probably lead to more protection of r-PA in protease-rich environments in various recombinant systems and potentially enhance its production and expression level.

Metabolic engineering of the diosgenin biosynthesis pathway in Trigonella foenum-graceum hairy root cultures.

Diosgenin as a triterpene with numbers of pharmaceutical applications has been identified in Trigonella foenum-graceum. In this survey, in order to scale up the amount of diosgenin in Fenugreek as a promising alternative of yam, ∆24-reductase as a rate limiting enzyme in diosgenin biosynthesis pathway has been overexpressed by utilizing pBI121 expression plasmid in hairy roots culture platform. The recombinant binary vector pBI121-∆24-reductase was transformed into R. rhizogenes strain ATCC 15834 to induce transgenic hairy roots in "Hamedan" as a low-diosgenin production genotype. In the transgenic hairy roots, the ∆24-reductase expression level was significantly 8.15 times overexpressed comparing to the non-transgenic hairy roots, Nonetheless the Sterol-methyltransferase, as a competitive enzyme, was 6 times downregulated. Furthermore, the expression rate of Squalene synthase, Cycloartenol synthase, C26-Hydroxylase were also increased 1.5, 1.7, 2.9 times higher than those of the non-transgenic hairy roots, respectively. The diosgenin content in the transgenic hairy root was raised 3 times up comparing to the non-transgenic hairy roots, besides it was scaled up 25-fold comparing to the diosgenin amount in "Hamedan" Leaf. As a result, the first metabolic engineering on this pathway was clearly revealed the impact of ∆24 -reductase gene in diosgenin content enhancement.

Improving water deficit tolerance of Salvia officinalis L. using putrescine.

To study the effects of foliar application of putrescine (distilled water (0), 0.75, 1.5, and 2.25 mM) and water deficit stress (20%, 40%, 60%, and 80% available soil water depletion (ASWD)) on the physiological, biochemical, and molecular attributes of Salvia officinalis L., a factorial experiment was performed in a completely randomized design with three replications in the growth chamber. The results of Real-Time quantitative polymerase chain reaction (qRT-PCR) analysis showed that putrescine concentration, irrigation regime, and the two-way interaction between irrigation regime and putrescine concentration significantly influenced cineole synthase (CS), sabinene synthase (SS), and bornyl diphosphate synthase (BPPS) relative expression. The highest concentration of 1,8-cineole, camphor, α-thujone, ß-thujone, CS, SS, and BPPS were obtained in the irrigation regime of 80% ASWD with the application of 0.75 mM putrescine. There was high correlation between expression levels of the main monoterpenes synthase and the concentration of main monoterpenes. The observed correlation between the two enzyme activities of ascorbate peroxidase (APX) and catalase (CAT) strongly suggests they have coordinated action. On the other hand, the highest peroxidase (PO) and superoxide dismutase (SOD) concentrations were obtained with the application of 0.75 mM putrescine under the irrigation regime of 40% ASWD. Putrescine showed a significant increase in LAI and RWC under water deficit stress. There was an increasing trend in endogenous putrescine when putrescine concentration was increased in all irrigation regimes. Overall, the results suggest that putrescine may act directly as a stress-protecting compound and reduced H2O2 to moderate the capacity of the antioxidative system, maintain the membrane stability, and increase secondary metabolites under water deficit stress.

Reteplase Fc-fusions produced in N. benthamiana are able to dissolve blood clots ex vivo.

Thrombolytic and fibrinolytic therapies are effective treatments to dissolve blood clots in stroke therapy. Thrombolytic drugs activate plasminogen to its cleaved form plasmin, a proteolytic enzyme that breaks the crosslinks between fibrin molecules. The FDA-approved human tissue plasminogen activator Reteplase (rPA) is a non-glycosylated protein produced in E. coli. rPA is a deletion mutant of the wild-type Alteplase that benefits from an extended plasma half-life, reduced fibrin specificity and the ability to better penetrate into blood clots. Different methods have been proposed to improve the production of rPA. Here we show for the first time the transient expression in Nicotiana benthamiana of rPA fused to the immunoglobulin fragment crystallizable (Fc) domain on an IgG1, a strategy commonly used to improve the stability of therapeutic proteins. Despite our success on the expression and purification of dimeric rPA-Fc fusions, protein instability results in high amounts of Fc-derived degradation products. We hypothesize that the "Y"- shape of dimeric Fc fusions cause steric hindrance between protein domains and leads to physical instability. Indeed, mutations of critical residues in the Fc dimerization interface allowed the expression of fully stable rPA monomeric Fc-fusions. The ability of rPA-Fc to convert plasminogen into plasmin was demonstrated by plasminogen zymography and clot lysis assay shows that rPA-Fc is able to dissolve blood clots ex vivo. Finally, we addressed concerns with the plant-specific glycosylation by modulating rPA-Fc glycosylation towards serum-like structures including α2,6-sialylated and α1,6-core fucosylated N-glycans completely devoid of plant core fucose and xylose residues.

Toward understanding of the methoxylated flavonoid biosynthesis pathway in Dracocephalum kotschyi Boiss.

Nowadays, with the development and advancement of next-generation sequencing technologies, a new path has been provided for transcriptomic studies. In this study, the transcriptome of Dracocephalum kotschyi Boiss., as an endemic and endangered plant which is contained a large amount of valuable secondary metabolites with antioxidant and anticancer properties, was sequenced. Then functional annotation and gene ontology analysis for 165,597 assembled transcripts were performed, most were associated with the metabolic pathways. This might be because there are various active biochemical pathways in this plant. Furthermore, after comprehensive transcript annotation, the putative genes involved in the main metabolic pathways of D. kotschyi were identified. Then, the biosynthetic pathway of its valuable methoxylated flavones was proposed. Finally, the accumulations of important methoxylated-flavone metabolites in three different tissues were quantified by HPLC. The relative expression of the genes involved in the proposed pathway was investigated by qRT-PCR, which indicated high expression levels in the bud tissue. The present results may lead to the design strategies to preserve the genetic diversity of endangered D. kotschyi plants and apply the new methods for engineering its valuable methoxylated-flavones pathway.

Identification and in-silico characterization of taxadien-5α-ol-O-acetyltransferase (TDAT) gene in Corylus avellana L.

Paclitaxel® (PC) is one of the most effective and profitable anti-cancer drugs. The most promising sources of this compound are natural materials such as tissue cultures of Taxus species and, more recently, hazelnut (Corylus avellana L.). A large part of the PC biosynthetic pathway in the yew tree and a few steps in the hazelnut have been identified. Since understanding the biosynthetic pathway of plant-based medicinal metabolites is an effective step toward their development and engineering, this paper aimed to identify taxadiene-5α-ol-O-acetyltransferase (TDAT) in hazelnut. TDAT is one of the key genes involved in the third step of the PC biosynthetic pathway. In this study, the TDAT gene was isolated using the nested-PCR method and then characterized. The cotyledon-derived cell mass induced with 150 µM of methyl jasmonate (MeJA) was utilized to isolate RNA and synthesize the first-strand cDNA. The full-length cDNA of TDAT is 1423 bp long and contains a 1302 bp ORF encoding 433 amino acids. The phylogenetic analysis of this gene revealed high homology with its ortholog genes in Quercus suber and Juglans regia. Bioinformatics analyses were used to predict the secondary and tertiary structures of the protein. Due to the lack of signal peptide, protein structure prediction suggested that this protein may operate at the cytoplasm. The homologous superfamily of the T5AT protein, encoded by TDAT, has two domains. The highest and lowest hydrophobicity of amino acids were found in proline 142 and lysine 56, respectively. T5AT protein fragment had 24 hydrophobic regions. The tertiary structure of this protein was designed using Modeler software (V.9.20), and its structure was verified based on the results of the Verify3D (89.46%) and ERRAT (90.3061) programs. The T5AT enzyme belongs to the superfamily of the transferase, and the amino acids histidine 164, cysteine 165, leucine 166, histidine 167, and Aspartic acid 168 resided at its active site. More characteristics of TDAT, which would aid PC engineering programs and maximize its production in hazelnut, were discussed.

Identification and expression analysis of S-alk(en)yl-L-cysteine sulfoxide lyase isoform genes and determination of allicin contents in Allium species.

Alliinase is the key enzyme in allicin biosynthesis pathway. In the current study, the identification and sequencing of alliinase genes along with determination of allicin contents were reported for Allium species with a novel report for Iranian endemic species. The presence of different isoforms in the Allium being discovered for the first time. In bulbs tissue, the highest allicin concentration was in Allium sativum, A. umbilicatum, and A. fistolosum (1.185%, 0.367%, and 0.34%, respectively), followed by A. spititatum (0.072%), A. lenkoranicum (0.055%), A. atroviolaseum (0.36%), A. rubellum (0.041%), and A. stamineum (0.007%). The highest allicin content in the leaves and roots were in A. sativum (0.13%), and A. stamineum (0.195%), respectively. The ORFs length ranged from 1416 in A. sativum (iso-alliinase2; ISA2) to 1523 bp in A. sativum (alliinase); the identity with A. sativum (alliinase) varies from 95% to 68% for A. ampeloprasum, and A. sativum (iso-alliinase1, ISA1) respectively. These data suggested that both ISA1 and ISA2 had a high expression in the roots and bulbs compared to A. sativum as the control in all species. Note that ISA1 and ISA2 were not expressed in the leaves. The results showed that isoforms expression patterns among different tissues in Allium species were variable. The presence of various isoforms is a possible explanation for the difference between the species in terms of obtained results, especially the amount of allicin.

New insights into diosgenin biosynthesis pathway and its regulation in Trigonella foenum-graecum L.

INTRODUCTION: Throughout history, thousands of medicinal and aromatic plants have been widely utilised by people worldwide. Owing to them possessing of valuable compounds with little side effects in comparison with chemical drugs, herbs have been of interest to humans for a number of purposes. Diosgenin, driven from fenugreek, Trigonella foenum-graecum L., has extensively drawn scientist's attention owing to having curable properties and being a precursor of steroid hormones synthesis. Nonetheless, complete knowledge about the biosynthesis pathway of this metabolite is still elusive. OBJECTIVE: In the present research, we isolated the full-length CDS of 14 genes involving in diosgenin formation and measured their expression rate in various genotypes, which had illustrated different amount of diosgenin. METHODOLOGY: The genes were successfully isolated, and functional motifs were also assessed using in silico approaches. RESULTS: Moreover, combining transcript and metabolite analysis revealed that there are many genes playing the role in diosgenin formation, some of which are highly influential. Among them, ∆24 -reductase, which converts cycloartenol to cycloartanol, is the first-committed and rate-limiting enzyme in this pathway. Additionally, no transcripts indicating to the presence or expression of lanosterol synthase were detected, contradicting the previous hypothesis about the biosynthetic pathway of diosgenin in fenugreek. CONCLUSION: Considering all these, therefore, we propose the most possible pathway of diosgenin. This knowledge will then pave the way toward cloning the genes as well as engineering the diosgenin biosynthesis pathway.

Comprehensive functional analysis and mapping of SSR markers in the chickpea genome (Cicer arietinum L.).

Plant molecular breeding largely depends on the relationship between molecular markers and major traits. Herein, a total of 32,962 genomic simple sequence repeats (SSRs) were detected in the whole genome of chickpea with an average density of 94.93 SSRs/Mb. Chickpea chromosomes uniformity test indicated that the genomic SSRs (gSSRs) were steadily distributed across the genome. Moreover, 48,667 transcriptome sequences were analyzed and 1949 SSR-containing transcript assembly contigs (TACs) were identified. The analysis showed that di- and trinucleotide SSRs were the most frequent SSR motifs within the transcriptome sequences. Among them, AT and TTA and AG and TTC motifs within the transcriptome showed the highest frequencies among di- and trinucleotide repeat motifs, respectively. The SSRs-containing TACs were compared to the GenBank non-redundant database using BLASTX, and subsequently, gene ontology (GO) analysis was performed using QuickGO browser to reduce complexity and highlight biological processes associated with the SSRs-containing TACs. The identified SSRs-containing TACs were categorized into 35 enriched functional-related gene group. The mapping of characterized SSRs-containing TACs onto chickpea chromosomes was performed using BLASTN. The mapping result showed that, a total of 1798 SSRs-containing TACs were mapped onto the chickpea genome. Based on the functional analysis result, 249 and 242 of the mapped SSRs-containing TACs were found in the genes encoding for putative stress-related proteins and transcription factors, respectively. The results presented here can be applied to improve and speed up the chickpea breeding programs.

Author Correction: Expression of key genes affecting artemisinin content in five Artemisia species.

A correction to this article has been published and is linked from the HTML and PDF versions of this paper. The error has been fixed in the paper.

An investigation of the HMGR gene and IPI gene expression in black caraway (Bunium persicum).

Black caraway is of great importance for its terpene compounds. Many genes are involved in the biosynthesis of secondary metabolites in medicinal plants. For this study, black caraway seeds were collected from five different regions, i.e. [Isfahan; Kerman (Khabr); Semnan; Kerman (Sirch); and Hormozgan]. The black caraway seed oil was extracted and analyzed by means of the gas chromatography method. There was a negatively significant correlation (p ≤ 0.05) observed between cuminaldehyde and gammaterpinene compounds. 3-hydroxy-3-methylglutaryl coenzyme A reductase (HMGR) and isopentenyl pyrophosphate isomerase (IPI) play an important role in the biosynthesis of secondary metabolites. Appropriate primers were designed for these genes based on the conserved regions in other plants. Amplified fragments were then sequenced. Blastn results indicated the similarity of the high RNA sequences between new sequences and other HMGR and IPI gene sequences in GenBank, and it also identified the HMGR and IPI gene sequences of B. persicum. A fragment of the HMGR gene with KJ143741 number was recorded in the gene bank. Quantitative PCR showed that the relative expression of two genes in different growth stages of B. persicum was significantly different between the germination stage and the multi-leaf stage, and also between the germination stage and the flowering stage (p < 0.05); however, there was no significant difference observed between the flowering stage and the multi-leaf stage. The results indicated that the expression of HMGR increased from the germination stage to the adult plant, and then it got stable until the flowering stage; in the same vein, the expression of IPI increased continuously from the germination stage to the flowering stage. The expression of HMGR and IPI genes occurred differently at the germination stage of five ecotypes. The Hormozgan ecotype showed the least expression rate.

Identification and expression analysis of a microRNA cluster derived from pre-ribosomal RNA in Papaver somniferum L. and Papaver bracteatum L.

Opium poppy (Papaver somniferum L.) is one of the ancient medical crops, which produces several important alkaloids such as morphine, noscapine, sanguinarine and codeine. MicroRNAs are endogenous non-coding RNAs that play important regulatory roles in plant diverse biological processes. Many plant miRNAs are encoded as single transcriptional units, in contrast to animal miRNAs, which are often clustered. Herein, using computational approaches, a total of 22 miRNA precursors were identified, which five of them were located as a clustered in pre-ribosomal RNA. Afterward, the transcript level of the precursor and the mature of clustered miRNAs in two species of the Papaveraceae family, i.e. P. somniferum L. and P. bracteatum L, were quantified by RT-PCR. With respect to obtained results, these clustered miRNAs were expressed differentially in different tissues of these species. Moreover, using target prediction and Gene Ontology (GO)-based on functional classification indicated that these miRNAs might play crucial roles in various biological processes as well as metabolic pathways. In this study, we discovered the clustered miRNA derived from pre-rRNA, which may shed some light on the importance of miRNAs in the plant kingdom.

Expression of key genes affecting artemisinin content in five Artemisia species.

Artemisinin, an effective anti-malarial drug is synthesized in the specialized 10-celled biseriate glandular trichomes of some Artemisia species. In order to have an insight into artemisinin biosynthesis in species other than A. annua, five species with different artemisinin contents were investigated for the expression of key genes that influence artemisinin content. The least relative expression of the examined terpene synthase genes accompanied with very low glandular trichome density (4 No. mm-2) and absence of artemisinin content in A. khorassanica (S2) underscored the vast metabolic capacity of glandular trichomes. A. deserti (S4) with artemisinin content of 5.13 mg g-1 DW had a very high expression of Aa-ALDH1 and Aa-CYP71AV1 and low expression of Aa-DBR2. It is possible to develop plants with high artemisinin synthesis ability by downregulating Aa-ORA in S4, which may result in the reduction of Aa-ALDH1 and Aa-CYP71AV1 genes expression and effectively change the metabolic flux to favor more of artemisinin production than artemisinic acid. Based on the results, the Aa-ABCG6 transporter may be involved in trichome development. S4 had high transcript levels and larger glandular trichomes (3.46 fold) than A. annua found in Iran (S1), which may be due to the presence of more 2C-DNA (3.48 fold) in S4 than S1.

Co-regulation analysis of co-expressed modules under cold and pathogen stress conditions in tomato.

A primary mechanism for controlling the development of multicellular organisms is transcriptional regulation, which carried out by transcription factors (TFs) that recognize and bind to their binding sites on promoter region. The distance from translation start site, order, orientation, and spacing between cis elements are key factors in the concentration of active nuclear TFs and transcriptional regulation of target genes. In this study, overrepresented motifs in cold and pathogenesis responsive genes were scanned via Gibbs sampling method, this method is based on detection of overrepresented motifs by means of a stochastic optimization strategy that searches for all possible sets of short DNA segments. Then, identified motifs were checked by TRANSFAC, PLACE and Soft Berry databases in order to identify putative TFs which, interact to the motifs. Several cis/trans regulatory elements were found using these databases. Moreover, cross-talk between cold and pathogenesis responsive genes were confirmed. Statistical analysis was used to determine distribution of identified motifs on promoter region. In addition, co-regulation analysis results, illustrated genes in pathogenesis responsive module are divided into two main groups. Also, promoter region was crunched to six subareas in order to draw the pattern of distribution of motifs in promoter subareas. The result showed the majority of motifs are concentrated on 700 nucleotides upstream of the translational start site (ATG). In contrast, this result isn't true in another group. In other words, there was no difference between total and compartmentalized regions in cold responsive genes.

Population structure, morphological and genetic diversity within and among melon (Cucumis melo L.) landraces in Iran.

ISSR markers were applied to evaluate the genetic diversity and differentiation of 270 individuals of 27 Iranian C. melo landraces of various varietal groups include vars. inodorous, cantalupensis, reticulatus, ameri, dudaim. Genetic diversity among the studied genotypes obtained by GeneAlex analysis (H = 0.08, I = 0.12, Na = 0.77, PPL = 22.6%). Cluster analysis divided Iranian melon landraces into two main cluster. Non-sweet genotype (dudaim group) was well separated from sweet genotypes (inodorous, ameri, reticulatus, cantalupensis). The most similar genotypes were BANI and TONI (0.95) and the most dissimilar ones were GER and TS (0.58). AMOVA result showed that the percentage of genetic variation among and within Iranian melon is 69% and 31%, respectively. All landraces evaluated based on 10 morphological traits which revealed the diversity of melon varietal groups. Bayesian analysis assigned ten landraces to Pop 1, eight landraces to Pop 2 and nine melon landraces to Pop 3. Bayesian and UPGMA cluster analyses demonstrated the almost related results. Our results indicated that ISSR markers technique alongside polyacrylamide gel analysis could be helpful to discriminate varieties of melon.

Assessment of Iranian endemic Artemisia khorassanica: karyological, genome size, and gene expressions involved in artemisinin production.

The species of Artemisia, one of the largest genera of the family Asteraceae, are frequently utilized for the treatment of diseases such as malaria, hepatitis, cancer, inflammation, and infections by fungi, bacteria, and viruses. Karyological studies were performed on 18 Artemisia khorassanica populations: eleven were diploid (2n = 18) and seven were tetraploid (2n = 36). The mean chromosome lengths were 3.61 and 3.84 µm for diploids and tetraploids, respectively. Two chromosome types ("m", "sm") formed karyotype formulas "18m" for diploids and "36m" and "34m + 2sm" for tetraploids. The mean 2C DNA contents were 5.91 and 11.53 pg in diploids and tetraploids, respectively. The transcription levels of key genes involved in artemisinin production were compared in diploid (B, D, H) and tetraploid (O, P, R) A. khorassanica relative to A. annua as a standard species. No artemisinin content was detected in diploid and tetraploid A. khorassanica populations. No significant diefrences were detected between diploids and tetraploids in terms of DXR , HMGR, FDS, and ADS gene expression. This implies that most of the genomic amplification likely occurs in the amount of repetitive DNA and not in unique sequences. The DBR2 gene was expressed in the diploid A. khorassanica in a low amount but silenced in the autotetraploid A. khorassanica.