Search

Background: The development of coronary artery disease (CAD) is influenced by sex and genetic factors. Genome-wide association studies (GWAS) have linked genetic loci to CAD, mostly in European populations. The study aims to find sex-related genetic differences in the Iranian population., Research Design and Methods: We conducted a sex-stratified GWAS with 4519 subjects (1832 males and 2687 females) in the discovery group and 922 subjects (495 males and 427 females) in the confirmation group of an Iranian cohort. We analyzed 9,141,124 variants using a genome-wide complex trait analysis (GCTA) tool., Results: We detected distinct genetic variants associated with CAD in males: rs34952209 [OR = 1.79; p  = 5.216E-8], rs1432687863 [OR = 1.95; p  = 8.477E-8], and in females, rs7314741 [OR = 1.67; p  = 7.142-8E] positively influenced CAD risk. The CAD-associated SNPs that were obtained have been confirmed using independent samples. Rs3495229 May impact histone mark and Pou2f2 motifs, while rs7314741 in the LEM Domain Containing 3 (LEMD3) promoter may affect a regulatory motif for the STAT transcription factor. According to Roadmap and ENCODE data, Rs1432687863 is a new variant affecting CAD in males, potentially through H3K9me3 in the heart., Conclusions: Our findings highlight the role of sex-specific genetic differences in CAD development, providing novel insights into disease pathways which is not appropriate using a sex-combined strategy. [Figure: see text].

This study aims to assess the effect of familial structures on the still-missing heritability estimate and prediction accuracy of Type 2 Diabetes (T2D) using pedigree estimated risk values (ERV) and genomic ERV. We used 11,818 individuals (T2D cases: 2,210) with genotype (649,932 SNPs) and pedigree information from the ongoing periodic cohort study of the Iranian population project. We considered three different familial structure scenarios, including (i) all families, (ii) all families with ≥ 1 generation, and (iii) families with ≥ 1 generation in which both case and control individuals are presented. Comprehensive simulation strategies were implemented to quantify the difference between estimates of [Formula: see text] and [Formula: see text]. A proportion of still-missing heritability in T2D could be explained by overestimation of pedigree-based heritability due to the presence of families with individuals having only one of the two disease statuses. Our research findings underscore the significance of including families with only case/control individuals in cohort studies. The presence of such family structures (as observed in scenarios i and ii) contributes to a more accurate estimation of disease heritability, addressing the underestimation that was previously overlooked in prior research. However, when predicting disease risk, the absence of these families (as seen in scenario iii) can yield the highest prediction accuracy and the strongest correlation with Polygenic Risk Scores. Our findings represent the first evidence of the important contribution of familial structure for heritability estimations and genomic prediction studies in T2D., (© 2024. The Author(s).)

Maturity-onset diabetes of the young (MODY) is an uncommon monogenic type of diabetes mellitus. Detecting genetic variants for MODY is a necessity for precise diagnosis and treatment. The majority of MODY genetic predisposition has been documented in European populations and a lack of information is present in Iranians which leads to misdiagnosis as a consequence of defects in unknown variants. In this study, using genetic variant information of 20,002 participants from the family-based TCGS (Tehran Cardiometabolic Genetic Study) cohort, we evaluated the genetic spectrum of MODY in Iran. We concentrated on previously discovered MODY-causing genes. Genetic variants were evaluated for their pathogenicity. We discovered 6 variants that were previously reported in the ClinVar as pathogenic/likely pathogenic (P/LP) for MODY in 45 participants from 24 families (INS in 21 cases, GCK in 13, HNF1B in 8, HNF4A, HNF1A, and CEL in 1 case). One potential MODY variant with Uncertain Risk Allele in ClinVar classification was also identified, which showed complete disease penetrance (100%) in four subjects from one family. This is the first family-based study to define the genetic spectrum and estimate the prevalence of MODY in Iran. The discovered variants need to be investigated by additional studies., (© 2024. The Author(s).)

FABP2 is one of the key genes involved in obesity development across different populations. However, there is no comprehensive report about the FABP2 contribution to obesity incidence among Iranians. Hence, the present study was designed to assess the probable role of FABP2 polymorphisms in obesity incidence in the Tehran Cardio- metabolic Genetic Study (TCGS) representative Iran population. Unrelated adults who had BMI information for at least 3 consecutive phases of the TCGS cohort were included. The control and case groups were defined as individuals who always had long-term persistent normal weight (20 < BMI < 25; n = 1526) and individuals who were long-term persistent obese (30 < BMI < 35; n = 1313), respectively. The logistic regression test was used to assess the possible association between SNPs located in and around the FABP2 gene with obesity. Also, we used Haploview and SHEsis to perform haplotype analysis to detect whether or not this chromosomal region is correlated with obesity. We found a gender-dependent association between the rs10857064 FABP2 and the risk of obesity. The presence of the rs10857064-G allele could significantly increase the risk of obesity only in women, not men (OR = 1.26; 95 % CI: 1.02-1.57; p = 0.03). Through haplotype analysis, we also detected that the TG haplotype containing rs7670862 and rs10857064 could significantly enhance the risk of obesity in women, further supporting the central role of rs10857064 in women's long-term obesity risk. In the current study, we revealed that rs10857064-G FABP2 can significantly predispose women to develop obesity. It highlights the importance of different genetic variants in both genders, which could help us to distinguish more efficient obesity screening tests and treatments based on gender in the future., Competing Interests: Declaration of Competing Interest The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper., (Copyright © 2023 Elsevier B.V. All rights reserved.)

The Tehran cardiometabolic genetic study (TCGS) is a large population-based cohort study that conducts periodic follow-ups. TCGS has created a comprehensive database comprising 20,367 participants born between 1911 and 2015 selected from four main ongoing studies in a family-based longitudinal framework. The study's primary goal is to identify the potential targets for prevention and intervention for non-communicable diseases that may develop in mid-life and late life. TCGS cohort focuses on cardiovascular, endocrine, metabolic abnormalities, cancers, and some inherited diseases. Since 2017, the TCGS cohort has augmented by encoding all health-related complications, including hospitalization outcomes and self-reports according to ICD11 coding, and verifying consanguineous marriage using genetic markers. This research provides an update on the rationale and design of the study, summarizes its findings, and outlines the objectives for precision medicine., (© 2023. Springer Nature B.V.)

Due to the central role in insulin secretion, the potassium inwardly-rectifying channel subfamily J member 11 (KCNJ11) gene is one of the essential genes for type 2 diabetes (T2D) predisposition. However, the relevance of this gene to T2D development is not consistent among diverse populations. In the current study, we aim to capture the possible association of common KCNJ11 variants across Iranian adults, followed by a meta-analysis. We found that the tested variants of KCNJ11 have not contributed to T2D incidence in Iranian adults, consistent with similar insulin secretion levels among individuals with different genotypes. The integration of our results with 72 eligible published case-control studies (41,372 cases and 47,570 controls) as a meta-analysis demonstrated rs5219 and rs5215 are significantly associated with the increased T2D susceptibility under different genetic models. Nevertheless, the stratified analysis according to ethnicity showed rs5219 is involved in the T2D risk among disparate populations, including American, East Asian, European, and Greater Middle Eastern, but not South Asian. Additionally, the meta-regression analysis demonstrated that the sample size of both case and control groups was significantly associated with the magnitude of pooled genetic effect size. The present study can expand our knowledge about the KCNJ11 common variant's contributions to T2D incidence, which is valuable for designing SNP-based panels for potential clinical applications in precision medicine. It also highlights the importance of similar sample sizes for avoiding high heterogeneity and conducting a more precise meta-analysis., (© 2022. The Author(s).)

Purpose: Fat mass and obesity-associated (FTO) is considered the first locus associated with adiposity, a concerning health problem worldwide. Many studies have evaluated the relationship between the FTO variants and obesity susceptibility. While the strong association of FTO rs1421085 with the risk of obesity across populations was reported in different studies, some researchers found a lack of association of this variant with adiposity. This systematic review and meta-analysis aimed to assess the association between obesity and rs1421085 polymorphism., Methods: We systematically searched PubMed, Scopus, and Google Scholar up to June 2022 to find pertinent studies. To further assess this issue, we surveyed the probable association of rs1421085 with obesity development among Iranian adults using the logistic regression analysis, and the obtained results were used for doing meta-analysis. After selection, nine eligible studies were included in the meta-analysis through the random- and fixed-effect models to determine the combined odds ratios (OR) and 95% confidence intervals (CI)., Results: According to our meta-analysis conducted on 5169 obese and 7772 non-obese individuals using different genetic models, including recessive, dominant, over-dominant, and additive, rs1421085 could positively increase the risk of obesity under all tested genetic models. Also, we detected a high to moderate level of heterogeneity among different studies under various genetic models., Conclusion: This meta-analysis further verified the positive association of FTO rs1421085 with the risk of developing obesity., Study Registration: This study is registered as PROSPERO CRD42021220092., Level of Evidence: Level I, systematic reviews and meta-analyses., (© 2022. The Author(s), under exclusive licence to Springer Nature Switzerland AG.)

Hedgehog (Hh) is a conserved signaling pathway that is involved in embryo development as well as adult tissue maintenance and repair in invertebrates and vertebrates. Abnormal activation of this pathway in various types of malignant drug- and apoptosis-resistant tumors has made it a therapeutic target against tumorigenesis. Thus, understanding the molecular mechanisms that promote the activation or inhibition of this pathway is critical. Long non-coding RNAs (lncRNAs), a subclass of non-coding RNAs with a length of > 200 nt, affect the expression of Hh signaling components via a variety of transcriptional and post-transcriptional processes. This review focuses on the crosstalk between lncRNAs and the Hh pathway in carcinogenesis, outlines the broad role of Hh-related lncRNAs in tumor progression, and illustrates their clinical diagnostic, prognostic, and therapeutic potential in tumors., (© 2022. The Author(s), under exclusive licence to Springer Science+Business Media, LLC, part of Springer Nature.)

Background: Molecular-based studies have revealed heterogeneity in Breast cancer BC while also improving classification and treatment. However, efforts are underway to distinguish between distinct subtypes of breast cancer. In this study, the results of several microarray studies were combined to identify genes and pathways specific to each BC subtype., Methods: Meta-analysis of multiple gene expression profile datasets was screened to find differentially expressed genes (DEGs) across subtypes of BC and normal breast tissue samples. Protein-protein interaction network and gene set enrichment analysis were used to identify critical genes and pathways associated with BC subtypes. The differentially expressed genes from meta-analysis was validated using an independent comprehensive breast cancer RNA-sequencing dataset obtained from the Cancer Genome Atlas (TCGA)., Results: We identified 110 DEGs (13 DEGs in all and 97 DEGs in each subtype) across subtypes of BC. All subtypes had a small set of shared DEGs enriched in the Chemokine receptor bind chemokine pathway. Luminal A specific were enriched in the translational elongation process in mitochondria, and the enhanced process in luminal B subtypes was interferon-alpha/beta signaling. Cell cycle and mitotic DEGs were enriched in the basal-like group. All subtype-specific DEG genes (100%) were successfully validated for Luminal A, Luminal B, ERBB2, and Normal-like. However, the validation percentage for Basal-like group was 77.8%., Conclusion: Integrating researches such as a meta-analysis of gene expression might be more effective in uncovering subtype-specific DEGs and pathways than a single-study analysis. It would be more beneficial to increase the number of studies that use matched BC subtypes along with GEO profiling approaches to reach a better result regarding DEGs and reduce probable biases. However, achieving 77.8% overlap in basal-specific genes and complete concordance in specific genes related to other subtypes can implicate the strength of our analysis for discovering the subtype-specific genes., (Copyright © 2022 Elsevier B.V. All rights reserved.)

Biological processes involving environmental and genetic factors drive the interplay between age- and sex-regulating lipid profile. The relation between variations in the LPA gene with increasing the risk of coronary heart disease is dependent on population differences, sex, and age. The present study tried to do a gene candidate association analysis in people with myocardial infarction (MI) in a 22 year cohort family-based longitudinal cohort study, Tehran Cardiometabolic Genetic Study (TCGS). After adjusting p value by the FDR method, only the association of rs6415084 with the MI probability and the age-of-CHD-onset was significant in males in their middle age (p < 0.005). Surprisingly, a lack of association was observed for the rest of the markers (16 SNPs). These results revealed the moderator effects of age and sex on the association between the genetic variants (SNPs) of LPA and heart disease risk. Our observations may provide new insights into the biology that underlies lipid profile with age or the sexual dimorphism of Lp(a) metabolism. Finally, Lp(a) appears to be an independent risk factor; however, the role of sex and ethnicity is important., (© 2022. The Author(s).)

This population-based longitudinal study is the first investigation that assesses the association of common MC4R SNPs with the obesity-related parameters over time and determines the effect of risk alleles during the three adulthood life periods (early, middle, and late) in a large Iranian cohort, a population with a unique genetic make-up that has been understudied and relatively unexplored. We obtained the genotype of 5370 unrelated adults who participated in the ongoing Tehran Cardiometabolic Genetic Study (TCGS) cohort project for the common MC4R SNPs. Linear regression and linear mixed model analyses were performed to examine the effect of MC4R polymorphisms on maximum BMI and other obesity-related factors over time. We recognized that several SNPs associated with the maximum BMI and the increased BMI, waist circumference, and waist-hip ratio across Iranian adults over a lifetime. Interestingly, we found that rs9954571-A has a yet unreported protective role against obesity-related factors, including BMI, waist circumference, waist-hip ratio, and triglyceride level. Additionally, a survey of the impact of the MC4R risk score throughout the adulthood life periods indicated that the MC4R risk score is influenced both the elevated BMI and waist circumference only during the early adulthood period. Our findings can expand our knowledge about the MC4R genetic variant's contributions to adulthood obesity and highlight the importance of evaluating the genetic components affecting obesity over a lifetime, which could be considered for obesity clinical screening and treatment., (Copyright © 2021 Elsevier B.V. All rights reserved.)

RNA-protein interactions of a virus play a major role in the replication of RNA viruses. The replication and transcription of these viruses take place in the cytoplasm of the host cell; hence, there is a probability for the host RNA-viral protein and viral RNA-host protein interactions. The current study applies a high-throughput computational approach, including feature extraction and machine learning methods, to predict the affinity of protein sequences of ten viruses to three categories of RNA sequences. These categories include RNAs involved in the protein-RNA complexes stored in the RCSB database, the human miRNAs deposited at the mirBase database, and the lncRNA deposited in the LNCipedia database. The results show that evolution not only tries to conserve key viral proteins involved in the replication and transcription but also prunes their interaction capability. These proteins with specific interactions do not perturb the host cell through undesired interactions. On the other hand, the hypermutation rate of NSP3 is related to its affinity to host cell RNAs. The Gene Ontology (GO) analysis of the miRNA with affiliation to NSP3 suggests that these miRNAs show strongly significantly enriched GO terms related to the known symptoms of COVID-19. Docking and MD simulation study of the obtained miRNA through high-throughput analysis suggest a non-coding RNA (an RNA antitoxin, ToxI) as a natural aptamer drug candidate for NSP5 inhibition. Finally, a significant interplay of the host RNA-viral protein in the host cell can disrupt the host cell's system by influencing the RNA-dependent processes of the host cells, such as a differential expression in RNA. Furthermore, our results are useful to identify the side effects of mRNA-based vaccines, many of which are caused by the off-label interactions with the human lncRNAs., (Copyright © 2021 Elsevier Ltd. All rights reserved.)

The long non-coding RNAs (lncRNAs) play a critical regulatory role in the host response to the viral infection. However, little is understood about the transcriptome architecture, especially lncRNAs pattern during the SARS-CoV-2 infection. In the present study, using publicly available RNA sequencing data of bronchoalveolar lavage fluid (BALF) and peripheral blood mononuclear cells (PBMC) samples from COVID-19 patients and healthy individuals, three interesting findings highlighted: (a) More than half of the interactions between lncRNAs-PCGs of BALF samples established by three trans-acting lncRNAs (HOTAIRM1, PVT1 and AL392172.1), which also exhibited the high affinity for binding to the SARS-CoV-2 genome, suggesting the major regulatory role of these lncRNAs during the SARS-CoV-2 infection. (b) lncRNAs of MALAT1 and NEAT1 are possibly contributed to the inflammation development in the SARS-CoV-2 infected cells. (c) In contrast to the 3' part of the SARS-CoV-2 genome, the 5' part can interact with many human lncRNAs. Therefore, the mRNA-based vaccines will not show any side effects because of the off-label interactions with the human lncRNAs. Overall, the putative functionalities of lncRNAs can be promising to design the non-coding RNA-based drugs and to inspect the efficiency of vaccines to overcome the current pandemic., (© 2021 The Authors. Journal of Cellular and Molecular Medicine published by Foundation for Cellular and Molecular Medicine and John Wiley & Sons Ltd.)

Recent genome-wide association studies (GWAS) highlighted the importance of genetic variations on SLC22A3 and MIA3 genes in developing coronary heart disease (CHD) among different ethnicities. However, the influence of these variations is not recognized within the Iranian population. Hence, in the present study, we aim to investigate two key single nucleotide polymorphisms (SNPs) on CHD incidence in this population. For this purpose, from Tehran Cardiometabolic Genetic Study (TCGS), 453 individuals with CHD were selected as a case and 453 individuals as a control that matched their age and gender. After quality control of two selected SNPs, rs2048327 (SLC22A3) and rs17465637 (MIA3), we used genotyps resulted from chip-typing technology and conducted the logistic regression analysis adjusted for non-genetic risk factors to detect the possible association of these SNPs with the CHD development. Our findings demonstrated the rs2048327-G and rs17465637-C can significantly increase the risk of CHD development about two times in only males and females, respectively. Interestingly, in the male carriers of the risk allele (G) of rs2048327, the low high-density lipoprotein (HDL) level can significantly predispose them to develop coronary heart disease in the future. According to our results, paying more attention to gender and genetic markers can help more efficient coronary heart disease screening and diagnosis., (Copyright © 2021 Elsevier B.V. All rights reserved.)

The genetic variations among individuals are one of the notable factors determining disease severity and drug response. Nowadays, COVID-19 pandemic has been adversely affecting many aspects of human life. We used the Tehran Cardio-Metabolic Genetic Study (TCGS) data that is an ongoing genetic study including the whole-genome sequencing of 1200 individuals and chip genotyping of more than 15,000 participants. Here, the effect of ACE2 variations by focusing on the receptor-binding site of SARS-CoV-2 and ACE2 cleavage by TMPRSS2 protease were investigated through simulations study. After analyzing TCGS data, 570 genetic variations on the ACE2 gene, including single nucleotide polymorphisms (SNP) and insertion/deletion (INDEL) were detected. Interestingly, two observed missense variants, K26R and S331F, which only the first one was previously reported, can reduce the receptor affinity for the viral Spike protein. Moreover, our bioinformatics simulation of 3D structures and docking of proteins explains important details of ACE2-Spike and ACE2-TMPRSS2 interactions, especially the critical role of Arg652 of ACE2 for protease function of TMPRSS2 was uncovered. As our results show that the genetic variation of ACE2 can at least influence the affinity of this receptor to its partners, we need to consider the genetic variations on ACE2 as well as other genes in the pathways that contribute to the pathogenesis of COVID-19 for designing efficient drugs and vaccines.

Type 2 diabetes (T2D) is emerging as one of the serious public health issues in both developed and developing counties. Here, we surveyed the worldwide population differentiation in T2D-associated variants and assessed the genetic burden of the disease in an ongoing Tehran Cardio-Metabolic Genetic Study (TCGS) cohort represented the Iranian population. We found multiple SNPs that were significantly depleted or enriched in at least one of the five populations of 1,000 Genome Project (African, American, East Asian, European, and South Asian) as well as the Iranian population. Interestingly, TCF7L2, a well-known associated gene with T2D, harbors the highest number of enriched risk alleles almost in all populations except for East Asian, where this gene embraces the largest number of significantly depleted risk alleles. The polygenic risk score (PRS) of the enriched risk alleles was calculated for 1,867 diabetic and 2,855 non-diabetic participants in the TCGS cohort, interestingly demonstrating that the risk of developing T2D was almost two times higher in top PRS quintile compared with the lowest quintile after adjusting for other known risk factors.

Long non-coding RNAs (lncRNAs) play crucial roles in regulating gene expression in response to plant stresses. Given the importance regulatory roles of lncRNAs, providing methods for predicting the function of these molecules, especially in non-model plants, is strongly demanded by researchers. Here, we constructed a reference sequence for lncRNAs in P. vera (Pistacia vera L.) with 53220 transcripts. In total, we identified 1909 and 2802 salt responsive lncRNAs in Ghazvini, a salt tolerant cultivar, after 6 and 24 h salt treatment, respectively and 1820 lncRNAs in Sarakhs, a salt sensitive cultivar, after 6 h salt treatment. Functional analysis of these lncRNAs by several hybrid methods, revealed that salt responsive NAT-related lncRNAs associated with transcription factors, CERK1, LEA, Laccase genes and several genes involved in the hormone signaling pathways. Moreover, gene ontology (GO) enrichment analysis of salt responsive target genes related to top five selected lncRNAs showed their involvement in the regulation of ATPase, cation transporter, kinase and UDP-glycosyltransferases genes. Quantitative real-time PCR (qRT-PCR) experiment results of lncRNAs, pre-miRNAs and mature miRNAs were in accordance with our RNA-seq analysis. In the present study, a comparative analysis of differentially expressed lncRNAs and microRNA precursors between salt tolerant and sensitive pistachio cultivars provides valuable knowledge on gene expression regulation under salt stress condition.

Constitutively photomorphogenic 1 (COP1) is an E3 ubiquitin ligase that has been studied extensively in the photomorphogenesis- and light-related processes in Arabidopsis. However, the possible role of COP1 in plant drought stress response remains unknown. Hence, in the present study, the stomatal behavior as one of the key elements in plant dehydration response was investigated in Arabidopsis cop1-4 and pea light-independent photomorphogenesis (lip1) mutants. We observed that water loss rate in the cop1-4 and lip1 detached leaves was significantly much faster than wild-type, resulting from failing to reduce the stomatal aperture by the mutants. But, interestingly, the cop1-4 and lip1 isolated leaves treated with abscisic acid (ABA) as well as cop1-4 and lip1 soil-grown under drought stress could close their stomata as wild-type. Hence, COP1 plays a fundamental role in the regulation of stomatal movements in response to dehydration and its function was conserved during evolution in both Arabidopsis and pea. Further evaluations showed the cop1-4 mutant was not significantly damaged from the oxidative stress derived from soil water limiting conditions when compared to wild-type. Similarly, the up-regulation level of several key stress-responsive genes was relatively lower in cop1-4 than wild-type. Therefore, COP1 might be considered as a potential key regulator of both short-and long-term dehydration response. Multiple stress-related cis-elements were also detected in the COP1 promoter region, which supported its up-regulation in response to drought, salt, and cold stresses. Besides, we figured out the constitutively open stomata of cop1-4 in darkness can be as a result of the reduced AtMYB61 expression., (Copyright © 2018. Published by Elsevier Masson SAS.)

The tree species, Pistacia vera (P. vera) is an important commercial product that is salt-tolerant and long-lived, with a possible lifespan of over one thousand years. Gene expression analysis is an efficient method to explore the possible regulatory mechanisms underlying these characteristics. Therefore, having the most suitable set of reference genes is required for transcript level normalization under different conditions in P. vera. In the present study, we selected eight widely used reference genes, ACT, EF1α, α-TUB, β-TUB, GAPDH, CYP2, UBQ10, and 18S rRNA. Using qRT-PCR their expression was assessed in 54 different samples of three cultivars of P. vera. The samples were collected from different organs under various abiotic treatments (cold, drought, and salt) across three time points. Several statistical programs (geNorm, NormFinder, and BestKeeper) were applied to estimate the expression stability of candidate reference genes. Results obtained from the statistical analysis were then exposed to Rank aggregation package to generate a consensus gene rank. Based on our results, EF1α was found to be the superior reference gene in all samples under all abiotic treatments. In addition to EF1α, ACT and β-TUB were the second best reference genes for gene expression analysis in leaf and root. We recommended β-TUB as the second most stable gene for samples under the cold and drought treatments, while ACT holds the same position in samples analyzed under salt treatment. This report will benefit future research on the expression profiling of P. vera and other members of the Anacardiaceae family.

The quality and quantity of RNA are critical for successful downstream transcriptome-based studies such as microarrays and RNA sequencing (RNA-Seq). RNA isolation from woody plants, such as Pistacia vera, with very high amounts of polyphenols and polysaccharides is an enormous challenge. Here, we describe a highly efficient protocol that overcomes the limitations posed by poor quality and low yield of isolated RNA from pistachio and various recalcitrant woody plants. The key factors that resulted in a yield of 150 μg of high quality RNA per 200 mg of plant tissue include the elimination of phenol from the extraction buffer, raising the concentration of β-mercaptoethanol, long time incubation at 65 °C, and nucleic acid precipitation with optimized volume of NaCl and isopropyl alcohol. Also, the A260/A280 and A260/A230 of extracted RNA were about 1.9-2.1and 2.2-2.3, respectively, revealing the high purity. Since the isolated RNA passed highly stringent quality control standards for sensitive reactions, including RNA sequencing and real-time PCR, it can be considered as a reliable and cost-effective method for RNA extraction from woody plants.