Selection of reference genes for quantitative gene expression normalization in flax (Linum usitatissimum L.)

Dual Transcriptomic and Metabolomic Analysis of Elicited Flax Sheds Light on the Kinetics of Immune Defense Activation Against the Biotrophic Pathogen Oidium lini

Flax (Linum usitatissimum) grown under controlled conditions displayed genotype-dependent resistance to powdery mildew (Oidium lini) following COS-OGA (comprising chitosan- and pectin-derived oligomers) elicitor application. The present study reveals a two-step immune response in plants preventively challenged with the elicitor: an initial, rapid response characterized by the transcription of defense genes whose protein products act in contact with or within the cell wall, where biotrophic pathogens initially thrive, followed by a prolonged activation of cell wall peroxidases and accumulation of secondary metabolites. Thus, dozens of genes encoding membrane receptors, pathogenesis-related proteins, and wall peroxidases were initially overexpressed. Repeated COS-OGA treatments had a transient effect on the transcriptome response while cumulatively remodeling the metabolome over time, with a minimum of two applications required for maximal metabolomic shifts. Secondary metabolites, in particular terpenoids and phenylpropanoids, emerged as major components of this secondary defense response alongside pathogenesis-related proteins and wall peroxidases. The sustained accumulation of secondary metabolites, even after cessation of elicitation, contrasted with the short-lived transcriptomic response. Wall peroxidase enzyme activity also exhibited cumulative effects, increasing strongly for weeks after a third elicitor treatment. This underscores the plasticity of the plant immune response in the face of a potential infection, and the need for repeated preventive applications to achieve the full protective potential of the elicitor.

Molecular Mechanism of Resistance to Alternaria alternata Apple Pathotype in Apple by Alternative Splicing of Transcription Factor MdMYB6-like

As a fruit tree with great economic value, apple is widely cultivated in China. However, apple leaf spot disease causes significant damage to apple quality and economic value. In our study, we found that MdMYB6-like is a transcription factor without auto-activation activity and with three alternative spliced variants. Among them, MdMYB6-like-β responded positively to the pathogen infection. Overexpression of MdMYB6-like-β increased the lignin content of leaves and improved the pathogenic resistance of apple flesh callus. In addition, all three alternative spliced variants of MdMYB6-like could bind to the promoter of MdBGLU H. Therefore, we believe that MdMYB6-like plays an important role in the infection process of the pathogen and lays a solid foundation for breeding disease-resistant cultivars of apple in the future.

Reliable reference gene selection for quantitative real-time PCR (qRT-PCR) in floral developmental phases of dioecious species Coccinia grandis

The flowering process is intricate and regulated by a combination of external and internal factors. Delving into gene expression research has the potential to enhance our comprehension of the molecular foundations underlying floral development. Because of its accuracy, specificity, reproducibility, and efficiency, qRT-PCR is now a biological research tool for studying expression pattern of desired genes. The gene expression investigations using qRT-PCR required a reference gene with relatively uniform expression levels in multiple biological samples, including different developmental stages, tissues, and experimental conditions. In this study, experimental sets offloral and floral organ development in the male and female plants of C. grandis, a dioecious Cucurbitaceae species, qRT-PCR profiling was performed using six reference genes as internal control with B-class floral identity gene, PISTILLATA (PI). To analyse the data, algorithms such as geNorm, NormFinder, RefFinder, and BestKeeper were used to pick out the best internal controls from a group of candidates. The optimal reference gene for qRT-PCR studies with floral samples has been recommended as β-actin combined with β-tubulin. This is the first report on the validation of candidate reference genes across flower developmental stages in the dioecious species C. grandis, which will provide basic data for research on the molecular mechanism underlying flower development in this species and lay the groundwork for similar studies in other related species.

Genome-Wide Analysis of Flax (Linum usitatissimum L.) Growth-Regulating Factor (GRF) Transcription Factors

Flax is an important cash crop globally with a variety of commercial uses. It has been widely used for fiber, oil, nutrition, feed and in composite materials. Growth regulatory factor (GRF) is a transcription factor family unique to plants, and is involved in regulating many processes of growth and development. Bioinformatics analysis of the GRF family in flax predicted 17 LuGRF genes, which all contained the characteristic QLQ and WRC domains. Equally, 15 of 17 LuGRFs (88%) are predicted to be regulated by lus-miR396 miRNA. Phylogenetic analysis of GRFs from flax and several other well-characterized species defined five clades; LuGRF genes were found in four clades. Most LuGRF gene promoters contained cis-regulatory elements known to be responsive to hormones and stress. The chromosomal locations and collinearity of LuGRF genes were also analyzed. The three-dimensional structure of LuGRF proteins was predicted using homology modeling. The transcript expression data indicated that most LuGRF family members were highly expressed in flax fruit and embryos, whereas LuGRF3, LuGRF12 and LuGRF16 were enriched in response to salt stress. Real-time quantitative fluorescent PCR (qRT-PCR) showed that both LuGRF1 and LuGRF11 were up-regulated under ABA and MeJA stimuli, indicating that these genes were involved in defense. LuGRF1 was demonstrated to be localized to the nucleus as expected for a transcription factor. These results provide a basis for further exploration of the molecular mechanism of LuGRF gene function and obtaining improved flax breeding lines.

Phosphorus and naphthalene acetic acid increased the seed yield by regulating carbon and nitrogen assimilation of flax

To evaluate the impact of phosphorus (P) combined with exogenous NAA on flax yield, enhance flax P utilization efficiency and productivity, minimize resource inputs and mitigate negative environmental and human effects. Therefore, it is crucial to comprehend the physiological and biochemical responses of flax to P and naphthylacetic acid (NAA) in order to guide future agronomic management strategies for increasing seed yield. A randomized complete block design trial was conducted under semi-arid conditions in Northwest China, using a factorial split-plot to investigate the effects of three P (0, 67.5, and 135.0 kg P2O5 ha-1) and three exogenous spray NAA levels (0, 20, and 40 mg NAA L-1) on sucrose phosphate synthase (SPS) and diphosphoribulose carboxylase (Rubisco) activities as well as nitrogen (N) and P accumulation and translocation in flax. Results indicated that the SPS and Rubisco activities, N and P accumulation at flowering and maturity along with assimilation and translocation post-flowering, fruiting branches per plant, tillers per plant, capsules per plant, and seed yield were 95, 105, 14, 27, 55, 15, 13, 110, 103, 82, 16, 61, 8, and 13% greater in the P treatments compared to those in the zero P treatment, respectively. Moreover, those characteristics were observed to be greater with exogenous spray NAA treatments compared to that no spray NAA treatment. Additionally, the maximum SPS and Rubisco activities, N and P accumulation, assimilation post-flowering and translocation, capsules per plant, and seed yield were achieved with the application of 67.5 kg P2O5 ha-1 with 20 mg NAA L-1. Therefore, these findings demonstrate that the appropriate combination of P fertilizer and spray NAA is an effective agronomic management strategy for regulating carbon and nitrogen assimilation by maintaining photosynthetic efficiency in plants to increase flax productivity.

Genome-Wide Identification and Expression Analysis of Auxin Response Factor Gene Family in Linum usitatissimum

Auxin response factors (ARFs) are critical components of the auxin signaling pathway, and are involved in diverse plant biological processes. However, ARF genes have not been investigated in flax (Linum usitatissimum L.), an important oilseed and fiber crop. In this study, we comprehensively analyzed the ARF gene family and identified 33 LuARF genes unevenly distributed on the 13 chromosomes of Longya-10, an oil-use flax variety. Detailed analysis revealed wide variation among the ARF family members and predicted nuclear localization for all proteins. Nineteen LuARFs contained a complete ARF structure, including DBD, MR, and CTD, whereas the other fourteen lacked the CTD. Phylogenetic analysis grouped the LuARFs into four (I-V) clades. Combined with sequence analysis, the LuARFs from the same clade showed structural conservation, implying functional redundancy. Duplication analysis identified twenty-seven whole-genome-duplicated LuARF genes and four tandem-duplicated LuARF genes. These duplicated gene pairs' K_a/K_s ratios suggested a strong purifying selection pressure on the LuARF genes. Collinearity analysis revealed that about half of the LuARF genes had homologs in other species, indicating a relatively conserved nature of the ARFs. The promoter analysis identified numerous hormone- and stress-related elements, and the qRT-PCR experiment revealed that all LuARF genes were responsive to phytohormone (IAA, GA3, and NAA) and stress (PEG, NaCl, cold, and heat) treatments. Finally, expression profiling of LuARF genes in different tissues by qRT-PCR indicated their specific functions in stem or capsule growth. Thus, our findings suggest the potential functions of LuARFs in flax growth and response to an exogenous stimulus, providing a basis for further functional studies on these genes.

Strigolactone signaling gene from soybean GmMAX2a enhances the drought and salt-alkaline resistance in Arabidopsis via regulating transcriptional profiles of stress-related genes

Strigolactone (SL) is a new plant hormone, which not only plays an important role in stimulating seed germination, plant branching, and regulating root development, but also plays an important role in the response of plants to abiotic stresses. In this study, the full-length cDNA of a soybean SL signal transduction gene (GmMAX2a) was isolated, cloned and revealed an important role in abiotic stress responses. Tissue-specific expression analysis by qRT-PCR indicated that GmMAX2a was expressed in all tissues of soybean, but highest expression was detected in seedling stems. Moreover, upregulation of GmMAX2a transcript expression under salt, alkali, and drought conditions were noted at different time points in soybean leaves compared to roots. Additionally, histochemical GUS staining studies revealed the deep staining in P_GmMAX2a: GUS transgenic lines compared to WT indicating active involvement of GmMAX2a promoter region to stress responses. To further investigate the function of GmMAX2a gene in transgenic Arabidopsis, Petri-plate experiments were performed and GmMAX2a OX lines appeared with longer roots and improved fresh biomass compared to WT plants to NaCl, NaHCO₃, and mannitol supplementation. Furthermore, the expression of several stress-related genes such as RD29B, SOS1, NXH1, AtRD22, KIN1, COR15A, RD29A, COR47, H+-APase, NADP-ME, NCED3, and P5CS were significantly high in GmMAX2a OX plants after stress treatment compared to WT plants. In conclusion, GmMAX2a improves soybean tolerance towards abiotic stresses (salt, alkali, and drought). Hence, GmMAX2a can be considered a candidate gene for transgenic breeding against various abiotic stresses in plants.

Drought and Oxidative Stress in Flax (Linum usitatissimum L.) Entails Harnessing Non-Canonical Reference Gene for Precise Quantification of qRT-PCR Gene Expression

Flax (Linum usitatissimum L.) is a self-pollinating, annual, diploid crop grown for multi-utility purposes for its quality oil, shining bast fiber, and industrial solvent. Being a cool (Rabi) season crop, it is affected by unprecedented climatic changes such as high temperature, drought, and associated oxidative stress that, globally, impede its growth, production, and productivity. To precisely assess the imperative changes that are inflicted by drought and associated oxidative stress, gene expression profiling of predominant drought-responsive genes (AREB, DREB/CBF, and ARR) was carried out by qRT-PCR. Nevertheless, for normalization/quantification of data obtained from qRT-PCR results, a stable reference gene is mandatory. Here, we evaluated a panel of four reference genes (Actin, EF1a, ETIF5A, and UBQ) and assessed their suitability as stable reference genes for the normalization of gene expression data obtained during drought-induced oxidative stress in flax. Taking together, from the canonical expression of the proposed reference genes in three different genotypes, we report that EF1a as a stand-alone and EF1a and ETIF5A in tandem are suitable reference genes to be used for the real-time visualization of cellular impact of drought and oxidative stress on flax.

Anthocyanin accumulation in grape berry flesh is associated with an alternative splicing variant of VvMYBA1

Anthocyanins are flavonoids that contribute to the color of grape berries and are an essential component of grape berry and wine quality. Anthocyanin accumulation in grape berries is dependent on the coordinated expression of genes encoding enzymes in the anthocyanin pathway that are principally regulated at the transcriptional level, with VvMYBA1 as the main transcriptional regulator in grapes. Alternative splicing (AS) events in VvMYBA1, however, have not been examined. In the present study, VvMYBA1-L, an AS variant of VvMYBA1, was identified in 'ZhongShan-Hong' (ZS-H) and its offspring. The AS variant is characterized by a deletion in the third exon of the open reading frame (ORF) of VvMYBA1-L, resulting in the early termination of the encoded protein. Overexpression of VvMYBA1-L in grape berries resulted in delayed flesh coloration and ectopic overexpression of VvMYBA1-L in tobacco inhibited the coloration of petals. Yeast two-hybrid and bimolecular fluorescence complementation (BiFC) assays revealed that VvMYBA1-L interacts with VvMYBA1. Dual luciferase assays indicated that co-infiltration of VvMYC1 and VvMYBA1 significantly activated the promoter regulated expression of VvCHS3, VvDFR, VvUFGT, and VvF3H. In the presence of VvMYBA1-L, however, the induction effect of VvMYBA1 on the indicated promoters was significantly inhibited. Our findings provide insight into the essential role of VvMYBA1 and its variant, VvMYBA1-L, in regulating anthocyanin accumulation in grape berry flesh.

Phenotypic analysis of Longya-10 × pale flax hybrid progeny and identification of candidate genes regulating prostrate/erect growth in flax plants

Flax is a dual-purpose crop that is important for oil and fiber production. The growth habit is one of the crucial targets of selection during flax domestication. Wild hybridization between cultivated flax and wild flax can produce superior germplasms for flax breeding and facilitate the study of the genetic mechanism underlying agronomically important traits. In this study, we used pale flax, Linum grandiflorum, and L. perenne to pollinate Longya-10. Only pale flax interspecific hybrids were obtained, and the trait analysis of the F₁ and F₂ generations showed that the traits analyzed in this study exhibited disparate genetic characteristics. In the F₁ generation, only one trait, i.e., the number of capsules per plant (140) showed significant heterosis, while the characteristics of other traits were closely associated with those of the parents or a decline in hybrid phenotypes. The traits of the F₂ generation were widely separated, and the variation coefficient ranged from 9.96% to 146.15%. The quantitative trait locus underlying growth habit was preliminarily found to be situated on chromosome 2 through Bulked-segregant analysis sequencing. Then linkage mapping analysis was performed to fine-map GH2.1 to a 23.5-kb interval containing 4 genes. Among them, L.us.o.m.scaffold22.109 and L.us.o.m.scaffold22.112 contained nonsynonymous SNPs with Δindex=1. Combined with the qRT-PCR results, the two genes might be possible candidate genes for GH2.1. This study will contribute to the development of important germplasms for flax breeding, which would facilitate the elucidation of the genetic mechanisms regulating the growth habit and development of an ideal architecture for the flax plant.

Genome-wide identification and characterization of NHL gene family in response to alkaline stress, ABA and MEJA treatments in wild soybean (Glycine soja)

Background

NDR1/HIN1-like (NHL) family genes are known to be involved in pathogen induced plant responses to biotic stress. Even though the NHL family genes have been identified and characterized in plant defense responses in some plants, the roles of these genes associated with the plant abiotic stress tolerance in wild soybean is not fully established yet, especially in response to alkaline stress.

Methods

We identified the potential NHL family genes by using the Hidden Markov model and wild soybean genome. The maximum-likelihood phylogenetic tree and conserved motifs were generated by using the MEME online server and MEGA 7.0 software, respectively. Furthermore, the syntenic analysis was generated with Circos-0.69. Then we used the PlantCARE online software to predict and analyze the regulatory cis-acting elements in promoter regions. Hierarchical clustering trees was generated using TM4: MeV4.9 software. Additionally, the expression levels of NHL family genes under alkaline stress, ABA and MEJA treatment were identified by qRT-PCR.

Results

In this study, we identified 59 potential NHL family genes in wild soybean. We identified that wild soybean NHL family genes could be mainly classified into five groups as well as exist with conserved motifs. Syntenic analysis of NHL family genes revealed genes location on 18 chromosomes and presence of 65 pairs of duplication genes. Moreover, NHL family genes consisted of a variety of putative hormone-related and abiotic stress responsive elements, where numbers of methyl jasmonate (MeJA) and abscisic acid (ABA) responsive elements were significantly larger than other elements. We confirmed the regulatory roles of NHL family genes in response to alkaline stress, ABA and MEJA treatment. In conclusion, we identified and provided valuable information on the wild soybean NHL family genes, and established a foundation to further explore the potential roles of NHL family genes in crosstalk with MeJA or ABA signal transduction mechanisms under alkaline stress.

RNA-Seq combined with population-level analysis reveals important candidate genes related to seed size in flax (Linum usitatissimum L.)

Seed size is a key determinant of crop yields. Understanding the regulatory mechanisms of seed size is beneficial for improving flax seed yield. In this study, the development of large flax seeds lagged behind that of small seeds, and 1,751 protein-coding genes were differentially expressed in early seeds, torpedo-stage embryos, and endosperms of CIli2719 and Z11637 using RNA sequencing. Homologous alignment revealed that 129 differentially expressed genes (DEGs) in flax were homologous with 71 known seed size-related genes in Arabidopsis thaliana and rice (Oryza sativa L.). These DEGs controlled seed size through multiple processes and factors, among which phytohormone pathways and transcription factors were the most important. Moreover, 54 DEGs were found to be associated with seed size and weight in a DEG-based association study. Nucleotide diversity (π) analysis of seed size-related candidate DEGs by homologous alignment and association analysis showed that the π values decreased significantly during flax acclimation from oil to fiber flax, suggesting that some seed size-related candidate genes were selected in this acclimation process. These results provide important resources and genetic foundation for further research on seed size regulation and seed improvement in flax.

Overexpression of WRINKLED1 improves the weight and oil content in seeds of flax (Linum usitatissimum L.)

Seeds of flax (Linum usitatissimum L.) are highly rich in both oil and linolenic acid (LIN). It is crucial for flax agricultural production to identify positive regulators of fatty acid biosynthesis. In this study, we find that WRINKLED1 transcription factors play important positive roles during flax seed oil accumulation. Two WRINKLED1 genes, LuWRI1a and LuWRI1b, were cloned from flax, and LuWRI1a was found be expressed predominantly in developing seeds during maturation. Overexpression of LuWRI1a increased seed size, weight, and oil content in Arabidopsis and increased seed storage oil content in transgenic flax without affecting seed production or seed oil quality. The rise in oil content in transgenic flax seeds was primarily attributable to the increase in seed weight, according to a correlational analysis. Furthermore, overexpression or interference of LuWRI1a upregulated the expression of genes in the fatty acid biosynthesis pathway and LAFL genes, and the expression level of WRI1 was highly significantly positively associated between L1L, LEC1, and BCCP2. Our findings give a theoretical scientific foundation for the future application of genetic engineering to enhance the oil content of plant seeds.

Selection and validation of reference genes for quantitative real-time PCR normalization in Psoralea corylifolia (Babchi) under various abiotic stress

Psoralea corylifolia L. is a popular herb and has long been used in traditional Ayurvedic and Chinese medicine owing to its extensive pharmacological activities, especially in the treatment of various shin diseases. To date, the systematic evaluation and selection of the optimum reference genes for gene expression analysis of P. corylifolia were not reported. Reverse transcription-quantitative polymerase chain reaction (RT-qPCR) is a method for gene expression quantification. Selecting appropriate reference genes is the prerequisite for accurate normalization of RT-qPCR results. This study assessed the expression stability of 10 candidate reference genes under different abiotic stresses. First, amplification primers for RT-qPCR were designed and received testing and optimization. Then, expression data from each candidate gene were evaluated based on three statistical algorithms, and their results were further integrated into a comprehensive ranking based on the geometric mean. Additionally, one target gene, i.e., 1-aminocyclopropane-1-carboxylate oxidase (ACO), was used to validate the effectiveness of the selected reference. Our analysis suggested that thioredoxin-like protein YLS8 (YLS8), TIP41-like family protein (TIP41), and cyclophilin 2 (CYP2) genes provided superior expression normalization under different abiotic stresses. Overall, this work constitutes the first effort to select optimal endogenous controls for RT-qPCR studies of gene expression in P. corylifolia. It also provides a reasonable normalization standard and basis for further analysis of the gene expression of bioactive components in P. corylifolia.

Identification and Validation of Reliable Reference Genes for Gene Expression Studies in Koelreuteria paniculata

RT-qPCR is considered a rapid and reliable technique for analyzing gene expression. This technique is commonly used to analyze the expression of various genes at diverse transcriptional levels in different samples. However, few studies have characterized ornamental Koelreuteria species for reliable reference genes. In this study, eight reference genes were evaluated as controls in RT-qPCR with SYBR green to quantify gene expression in different Koelreuteria paniculata samples. All selected reference genes showed a broad range of Ct values in all samples, which was supportive of their variable expression. Our results showed significant variation in the stable expression of K. paniculata genes. Sample data, analyzed using geNorm, NormFinder, and BestKeeper, showed that phospholipase (PLA2) and β-actin (ACT) were the most suitable and statistically reliable reference genes, whereas ribosomal protein L13 (RPL13) and elongation factor 1-α (EF1α) were less stable and unsuitable for use as internal controls. To compare gene expression levels, two or more reference genes should be used for data normalization. Thus, the stability and expression of both PLA2 and ACT were believed to provide better normalization and quantification of the transcript levels for gene expression studies in K. paniculata.

Selection and Evaluation of Candidate Reference Genes for Quantitative Real-Time PCR in Aboveground Tissues and Drought Conditions in Rhododendron Delavayi

There has been no systematic identification and screening of candidate reference genes for normalization of quantitative real-time PCR (qRT-PCR) results in Rhododendron delavayi to date. Therefore, the present study used GAPDH, Act, EF1, Tub-, Tub-5, UEC1, TATA, TATA-2, UEP, TIP41, and Ubiquitin to predict their stabilities on different aboveground tissues (matured leaves (ML), stem tips (STM), and flower buds (FB)) at different developmental stages (young and adult plants) using five statistical algorithms: Delta Ct method, BestKeeper, geNorm, Normfinder, and RefFinder. The findings were confirmed using ML obtained from plants that had been stressed by drought. By using RefFinder with ML samples collected under drought conditions, it was determined that the top five most stable reference genes were GAPDH > UEC1 > Actin > Tubulin- > Tubulin—5, whereas the least stable reference gene was Ubiquitin. In addition, under control conditions, UEC1, UEC2, Actin, and GAPDH were selected as the highest stable potential reference genes at the juvenile stage of R. delavayi with ML and STM. When ML and STM were combined with drought-stressed samples, TIP41, GAPDH, or their combination proved to be the most effective qRT-PCR primers. The findings will aid in the improvement of the precision and reliability of qRT-PCR data and laying the groundwork for future gene functional studies in R. delavayi.

Selection and Validation of Reference Genes for Quantitative Real-Time PCR Normalization in Athetis dissimilis (Lepidoptera: Noctuidae) Under Different Conditions

Reference genes are the key to study gene expression patterns using quantitative real-time PCR (qRT-PCR). No studies on the reference genes of Athetis dissimilis, an important agricultural pest, have been reported. In order to determine the reference genes for qRT-PCR normalization in A. dissimilis under different conditions, 10 candidate genes [18S ribosomal protein (18S), 28S ribosomal protein (28S), arginine kinase (AK), elongation factor 1 alpha (EF1-α), glyceraldehyde-3-phosphate dehydrogenase (GAPDH), ribosomal protein L32 (RPL32), ribosomal protein L40 (RPL40), alpha-tubulin (α-TUB), beta-actin (β-ACT), and beta-tubulin (β-TUB)] of A. dissimilis were selected to evaluate their stability as reference genes under different biotic and abiotic conditions by using five tools, geNorm, NormFinder, BestKeeper, ΔCt, and RefFinder. Furthermore, CSP1 and superoxide dismutase (SOD) were used as target genes to validate the candidate reference genes. The results showed that different reference genes were needed under different experimental conditions, among which, EF-1α, RPL40, and 18S are most suitable reference genes for studying genes related development stages of A. dissimilis, RPL40 and α-TUB for larval tissues, α-TUB and 28S for adult tissues, EF-1α and β-ACT for insecticidal treatments, β-ACT and 28S for temperature treatments, EF-1α and β-ACT for starvation treatments, RPL40 and 18S for dietary treatments, and 18S, 28S, and α-TUB for all the samples. These results provide suitable reference genes for studying gene expression in A. dissimilis under different experimental conditions, and also lay the foundation for further research into the function of related genes in A. dissimilis.

Selection and Validation of Reliable Reference Genes for qRT-PCR Normalization of Bursaphelenchus xylophilus from Different Temperature Conditions and Developmental Stages

Quantitative reverse transcription polymerase chain reaction (qRT-PCR) is a powerful technique for studying gene expression. The key to quantitative accuracy depends on the stability of the reference genes used for data normalization under different experimental conditions. Pine wood nematode (PWN; Bursaphelenchus xylophilus) is the causal agent of the devastating pine wilt disease (PWD). Extensive and prompt research is needed to understand the molecular mechanism of PWD, but identification of the reference PWN genes for standardized qRT-PCR has not been reported yet. We have analyzed eight candidate reference genes of PWN across different temperature conditions and developmental stages. Delta Ct method, GeNorm, NormFinder, BestKeeper, and RefFinder algorithms were used to evaluate the stability of expression of these genes. Finally, we use arginine kinase gene (AK) in different temperatures and heat shock protein 90 (HSP90) in different developmental stages to confirm the stability of these genes. UBCE and EF1γ were most stable across different temperature treatments, whereas EF1γ and Actin were most stable across different developmental stages. In general, these results indicate that EF1γ is the most stable gene for qRT-PCR under different conditions. The systematic analysis of qRT-PCR reference gene selection will be helpful for future functional analysis and exploration of B. xylophilus genetic resources.

Selection and Validation of Reference Genes for Quantitative Real-Time PCR Analysis of Development and Tissue-Dependent Flower Color Formation in Cymbidium lowianum

The development and tissue-dependent color formation of the horticultural plant results in various color pattern flowers. Anthocyanins and carotenoids contribute to the red and yellow colors, respectively. In this study, quantitative real-time polymerase chain reaction (qRT-PCR) is used to analyze the expression profiles of anthocyanin and carotenoids biosynthesis genes in Cymbidium lowianum (Rchb.f.) Rchb.f. Appropriate reference gene selection and validation are required before normalization of gene expression in qRT-PCR analysis. Thus, we firstly selected 12 candidate reference genes from transcriptome data, and used geNorm and Normfinder to evaluate their expression stability in lip (divided into abaxial and adaxial), petal, and sepal of the bud and flower of C. lowianum. Our results show that the two most stable reference genes in different tissues of C. lowianum bud and flower are EF1δ and 60S, the most unstable reference gene is 26S. The expression profiles of the CHS and BCH genes were similar to FPKM value profiles after normalization to the two most stable reference genes, EF1δ and 60S, with the upregulated CHS and BCH expression in flower stage, indicating that the ABP and CBP were activated across the stages of flower development. However, when the most unstable reference gene, 26S, was used to normalize the qRT-PCR data, the expression profiles of CHS and BCH differed from FPKM value profiles, indicating the necessity of selecting stable reference genes. Moreover, CHS and BCH expression was highest in the abaxial lip and adaxial lip, respectively, indicating that the ABP and CBP were activated in abaxial and adaxial lip, respectively, resulting in a presence of red or yellow segments in abaxial and adaxial lip. This study is the first to provide reference genes in C. lowianum, and also provide useful information for studies that aim to understand the molecular mechanisms of flower color formation in C. lowianum.

Genome-wide analysis of genomic imprinting in the endosperm and allelic variation in flax

Genomic imprinting is an epigenetic phenomenon that causes biased expression of maternally and paternally inherited alleles. In flowering plants, genomic imprinting predominantly occurs in the triploid endosperm and plays a vital role in seed development. In this study, we identified 248 candidate imprinted genes including 114 maternally expressed imprinted genes (MEGs) and 134 paternally expressed imprinted genes (PEGs) in flax (Linum usitatissimum L.) endosperm using deep RNA sequencing. These imprinted genes were neither clustered in specific chromosomal regions nor well conserved among flax and other plant species. MEGs tended to be expressed specifically in the endosperm, whereas the expression of PEGs was not tissue-specific. Imprinted single nucleotide polymorphisms differentiated 200 flax cultivars into the oil flax, oil-fiber dual purpose flax and fiber flax subgroups, suggesting that genomic imprinting contributed to intraspecific variation in flax. The nucleotide diversity of imprinted genes in the oil flax subgroup was significantly higher than that in the fiber flax subgroup, indicating that some imprinted genes underwent positive selection during flax domestication from oil flax to fiber flax. Moreover, imprinted genes that underwent positive selection were related to flax functions. Thirteen imprinted genes related to flax seed size and weight were identified using a candidate gene-based association study. Therefore, our study provides information for further exploration of the function and genomic variation of imprinted genes in the flax population.

Genome-wide identification and expression analysis of the WRKY transcription factor family in flax (Linum usitatissimum L.)

BACKGROUND

Members of the WRKY protein family, one of the largest transcription factor families in plants, are involved in plant growth and development, signal transduction, senescence, and stress resistance. However, little information is available about WRKY transcription factors in flax (Linum usitatissimum L.).

RESULTS

In this study, comprehensive genome-wide characterization of the flax WRKY gene family was conducted that led to prediction of 102 LuWRKY genes. Based on bioinformatics-based predictions of structural and phylogenetic features of encoded LuWRKY proteins, 95 LuWRKYs were classified into three main groups (Group I, II, and III); Group II LuWRKYs were further assigned to five subgroups (IIa-e), while seven unique LuWRKYs (LuWRKYs 96-102) could not be assigned to any group. Most LuWRKY proteins within a given subgroup shared similar motif compositions, while a high degree of motif composition variability was apparent between subgroups. Using RNA-seq data, expression patterns of the 102 predicted LuWRKY genes were also investigated. Expression profiling data demonstrated that most genes associated with cellulose, hemicellulose, or lignin content were predominantly expressed in stems, roots, and less in leaves. However, most genes associated with stress responses were predominantly expressed in leaves and exhibited distinctly higher expression levels in developmental stages 1 and 8 than during other stages.

CONCLUSIONS

Ultimately, the present study provides a comprehensive analysis of predicted flax WRKY family genes to guide future investigations to reveal functions of LuWRKY proteins during plant growth, development, and stress responses.

Identification and Evaluation of Suitable Reference Genes for RT-qPCR Analysis in Hippodamia variegata (Coleoptera: Coccinellidae) Under Different Biotic and Abiotic Conditions

Reverse transcriptase-quantitative polymerase chain reaction (RT-qPCR) is an accurate and convenient technique for quantifying expression levels of the target genes. Selection of the appropriate reference gene is of the vital importance for RT-qPCR analysis. Hippodamia variegata is one of the most important predatory natural enemies of aphids. Recently, transcriptome and genome sequencings of H. variegata facilitate the gene functional studies. However, there has been rare investigation on the detection of stably expressed reference genes in H. variegata. In the current study, by using five analytical tools (Delta Ct, geNorm, NormFinder, BestKeeper, and RefFinder), eight candidate reference genes, namely, Actin, EF1α, RPL7, RPL18, RPS23, Tubulin-α, Tubulin-β, and TufA, were evaluated under four experimental conditions including developmental stages, tissues, temperatures, and diets. As a result, a specific set of reference genes were recommended for each experimental condition. These findings will help to improve the accuracy and reliability of RT-qPCR data, and lay a foundation for further exploration on the gene function of H. variegata.

Genome-wide identification and expression analysis of detoxification efflux carriers (DTX) genes family under abiotic stresses in flax

The detoxification efflux carriers (DTX)/multidrug and toxic compound extrusion (MATE) transporters encompass an ancient gene family of secondary transporters involved in the process of plant detoxification. A genome-wide analysis of these transporters was carried out in order to better understand the transport of secondary metabolites in flaxseed genome (Linum usitassimum). A total of 73 genes coding for DTX/MATE transporters were identified. Gene structure, protein domain and motif organization were found to be notably conserved over the distinct phylogenetic groups, showing the evolutionary significant role of each class. Gene ontology (GO) annotation revealed a link to transporter activities, response to stimulus and localizations. The presence of various hormone and stress-responsive cis-regulatory elements in promoter regions could be directly correlated with the alteration of their transcripts. Tertiary structure showed conservation for pore size and constrains in the pore, which indicate their involvement in the exclusion of toxic substances from the cell. MicroRNA target analysis revealed that LuDTXs genes were targeted by different classes of miRNA families. Twelve LuDTX genes were chosen for further quantitative real-time polymerase chain reaction analysis in response to cold, salinity and cadmium stress at 0, 6, 12 and 24 hours after treatment. Altogether, the identified members of the DTX gene family, their expression profile, phylogenetic and miRNAs analysis might provide opportunities for future functional validation of this important gene family in flax.

Appropriate Reference Genes for RT-qPCR Normalization in Various Organs of Anemone flaccida Fr. Schmidt at Different Growing Stages

Anemone flaccida Fr. Schmidt is a traditional medicinal herb in southwestern China and has multiple pharmacological effects on bruise injuries and rheumatoid arthritis (RA). A new drug with a good curative effect on RA has recently been developed from the extract of A. flaccida rhizomes, of which the main medicinal ingredients are triterpenoid saponins. Due to excessive exploitation, the wild population has been scarce and endangered in a few of its natural habitats and research on the cultivation of the plant commenced. Studies on the gene expressions related to the biosynthesis of triterpenoid saponins are not only helpful for understanding the effects of environmental factors on the medicinal ingredient accumulations but also necessary for monitoring the herb quality of the cultivated plants. Reverse transcription quantitative polymerase chain reaction (RT-qPCR) as a sensitive and powerful technique has been widely used to detect gene expression across tissues in plants at different stages; however, its accuracy and reliability depend largely on the reference gene selection. In this study, the expressions of 10 candidate reference genes were evaluated in various organs of the wild and cultivated plants at different stages, using the algorithms of geNorm, NormFinder and BestKeeper, respectively. The purpose of this study was to identify the suitable reference genes for RT-qPCR detection in A. flaccida. The results showed that two reference genes were sufficient for RT-qPCR data normalization in A. flaccida. PUBQ and ETIF1a can be used as suitable reference genes in most organs at various stages because of their expression stabilitywhereas the PUBQ and EF1Α genes were desirable in the rhizomes of the plant at the vegetative stage.

Selection and Validation of Reference Genes for Quantitative RT-PCR Analysis in Corylus heterophylla Fisch. × Corylus avellana L

(1) Background: the species of Corylus have sporophytic type of self-incompatibility. Several genes related to recognition reaction between pollen and stigma have been identified in hazelnuts. To better understand the self-incompatibility (SI) response, we screened the suitable reference genes by using quantitative real-time reverse transcription PCR (qRT-PCR) analysis in hazelnut for the first time. (2) Methods: the major cultivar "Dawei" was used as material. A total of 12 candidate genes were identified and their expression profiles were compared among different tissues and in response to various treatments (different times after self- and cross-pollination) by RT-qPCR. The expression stability of these 12 candidate reference genes was evaluated using geNorm, NormFinder, BestKeeper, Delta Ct, and RefFinder programs. (3) Results: the comprehensive ranking of RefFinder indicated that ChaActin, VvActin, ChaUBQ14, and ChaEF1-α were the most suitable reference genes. According to the stability analysis of 12 candidate reference genes for each sample group based on four software packages, ChaActin and ChaEF1-α were most stable in different times after self-pollination and 4 h after self- and cross-pollination, respectively. To further validate the suitability of the reference genes identified in this study, CavPrx, which the expression profiles in Corylus have been reported, was quantified by using ChaActin and ChaEF1-α as reference genes. (4) Conclusions: our study of reference genes selection in hazelnut shows that the two reference genes, ChaActin and ChaEF1-α, are suitable for the evaluation of gene expression, and can be used for the analysis of pollen-pistil interaction in Corylus. The results supply a reliable foundation for accurate gene quantifications in Corylus species, which will facilitate the studies related to the reproductive biology in Corylus.

Identification of Tomato Ve1 Homologous Proteins in Flax and Assessment for Race-Specific Resistance in Two Fiber FlaxCultivars against Verticillium dahliae Race 1

In the last decade, the soil borne fungal pathogen Verticillium dahliae has had an increasingly strong effect on fiber flax (Linum usitatissimum L.), thus causing important yield losses in Normandy, France. Race-specific resistance against V. dahliae race 1 is determined by tomato Ve1, a leucine-rich repeat (LRR) receptor-like protein (RLP). Furthermore, homologous proteins have been found in various plant families. Herein, four homologs of tomato Ve1 were identified in the flax proteome database. The selected proteins were named LuVe11, LuVe12, LuVe13 and LuVe14 and were compared to other Ve1. Sequence alignments and phylogenic analysis were conducted and detected a high similarity in the content of amino acids and that of the Verticillium spp. race 1 resistance protein cluster. Annotations on the primary structure of these homologs reveal several features of tomato Ve1, including numerous copies of a 28 amino acids consensus motif [XXIXNLXXLXXLXLSXNXLSGXIP] in the LRR domain. An in vivo assay was performed using V. dahliae race 1 on susceptible and tolerant fiber flax cultivars. Despite the presence of homologous genes and the stronger expression of LuVe11 compared to controls, both cultivars exhibited symptoms and the pathogen was observed within the stem. Amino acid substitutions within the segments of the LRR domain could likely affect the ligand binding and thus the race-specific resistance. The results of this study indicate that complex approaches including pathogenicity tests, microscopic observations and gene expression should be implemented for assessing race-specific resistance mediated by Ve1 within the large collection of flax genotypes.

Selection and validation of reference genes for quantitative gene expression normalization in Taxus spp

Quantitative real-time PCR (qRT-PCR) is commonly used to measure gene expression to further explore gene function, while suitable reference genes must be stably expressed under different experimental conditions to obtain accurate and reproducible data for relative quantification. Taxol or paclitaxel is an important anticancer compound mainly identified in Taxus spp. The molecular mechanism of the regulation of taxol biosynthesis is current research goal. However, in the case of Taxus spp., few reports were published on screening suitable reference genes as internal controls for qRT-PCR. Here, eight reference genes were selected as candidate reference genes for further study. Common statistical algorithms geNorm, NormFinder, BestKeeper, ΔCt, and RefFinder were used to analyze the data from samples collected from a cell line of Taxus × media under various experimental conditions and from tissues of Taxus chinensis var. mairei. The expression patterns of TcMYC under salicylic acid treatment differed significantly, with the best and worst reference genes in the cell line. This study screened out suitable reference genes (GAPDH1 and SAND) under different treatments and tissues for the accurate and reliable normalization of the qRT-PCR expression data of Taxus spp. At the same time, this study will aid future research on taxol biosynthesis-related genes expression in Taxus spp., and can also be directly used to other related species.

Identification of two compounds able to improve flax resistance towards Fusarium oxysporum infection

Plants are surrounded by a diverse range of microorganisms that causes serious crop losses and requires the use of pesticides. Flax is a major crop in Normandy used for its fibres and is regularly challenged by the pathogenic fungus Fusarium oxysporum (Fo) f. sp. lini. To protect themselves, plants use "innate immunity" as a first line of defense level against pathogens. Activation of plant defense with elicitors could be an alternative for crop plant protection. A previous work was conducted by screening a chemical library and led to the identification of compounds able to activate defense responses in Arabidopsis thaliana. Four compounds were tested for their abilities to improve resistance of two flax varieties against Fo. Two of them, one natural (holaphyllamine or HPA) and one synthetic (M4), neither affected flax nor Fo growth. HPA and M4 induced oxidative burst and callose deposition. Furthermore, HPA and M4 caused changes in the expression patterns of defense-related genes coding a glucanase and a chitinase-like. Finally, plants pre-treated with HPA or M4 exhibited a significant decrease in the disease symptoms. Together, these findings demonstrate that HPA and M4 are able to activate defense responses in flax and improve its resistance against Fo infection.

The Toolbox for Fiber Flax Breeding: A Pipeline From Gene Expression to Fiber Quality

The goal of any plant breeding program is to improve quality of a target crop. Crop quality is a comprehensive feature largely determined by biological background. To improve the quality parameters of crops grown for the production of fiber, a functional approach was used to search for genes suitable for the effective manipulation of technical fiber quality. A key step was to identify genes with tissue and stage-specific pattern of expression in the developing fibers. In the current study, we investigated the relationship between gene expression evaluated in bast fibers of developing flax plants and the quality parameters of technical fibers measured after plant harvesting. Based on previously published transcriptomic data, two sets of genes that are upregulated in fibers during intrusive growth and tertiary cell wall deposition were selected. The expression level of the selected genes and fiber quality parameters were measured in fiber flax, linseed (oil flax) cultivars, and wild species that differ in type of yield and fiber quality parameters. Based on gene expression data, linear regression models for technical stem length, fiber tensile strength, and fiber flexibility were constructed, resulting in the identification of genes that have high potential for manipulating fiber quality. Chromosomal localization and single nucleotide polymorphism distribution in the selected genes were characterized for the efficacy of their use in conventional breeding and genome editing programs. Transcriptome-based selection is a highly targeted functional approach that could be used during the development of new cultivars of various crops.

Identification and Expression Analysis of Strigolactone Biosynthetic and Signaling Genes in Response to Salt and Alkaline Stresses in Soybean (Glycine max)

Strigolactones (SLs) are the major plant hormones that play important roles in regulating organ development and environmental stress tolerance in plants. Even though the SL-related genes have been identified and well characterized in some plants, the information of SL-related genes in soybean is not fully established yet, especially in response to salt and alkaline stresses. In this study, we identified nine SL biosynthetic genes that include two D27, two CCD7, two CCD8, and three MAX1, as well as seven SL signaling genes that comprised two D14, two MAX2, and three D53 in the soybean genome. We found that SL biosynthetic and signaling genes are evolutionary conserved among different species. Syntenic analysis of these genes revealed their location on nine chromosomes as well as the presence of 10 pairs of duplication genes. Moreover, plant hormone and stress-responsive elements were identified in the promoter regions of SL biosynthetic and signaling genes. By using reverse transcription quantitative real-time PCR, we confirmed that SL genes have different tissue expressions in roots, stems, and leaves. The expression profile of SL biosynthetic and signaling genes under salt and alkaline stresses further confirmed the regulatory roles of SL biosynthetic and signaling genes under stress. In conclusion, we identified and provided valuable information on the soybean SL biosynthetic and signaling genes, and established a foundation for further functional analysis of soybean SL-related genes in response to salt and alkaline stresses.

Selection and evaluation of appropriate reference genes for RT-qPCR based expression analysis in Candida tropicalis following azole treatment

Candida tropicalis arises as one of the predominant non-Candida albicans Candida (NCAC) species causing invasive candidiasis in Asian countries. A rise in reports of C. tropicalis with a parallel increase in fluconazole resistance has also been observed. The genes and underlying pathways associated with azole antifungal resistance in C. tropicalis is still not properly understood. The RT-qPCR is the most promising approach for expression analysis of target genes to understand the mechanisms of resistance. The reliability and reproducibility of this technique depend on the selection of suitable reference genes for the normalization in expression study. The present study investigated the expression stability levels of ten genes including ACT1, EF1, GAPDH, PGK1, RDN5.8, RDN18, RDN28, SDHA, TUB1, and UBC13 for their suitability in fluconazole treated/untreated C. tropicalis. The stability levels of these genes were examined by the ∆∆CT, ΔCT, Pfaffl methods and five independent software including hkgFinder, geNorm, NormFinder, BestKeeper, and RefFinder software. We report, the EF1 and ACT1 were the most stable reference genes for normalization and can be used for the gene expression analysis in C. tropicalis. To the best of our knowledge, our study is the first to select and validate the reference genes in C. tropicalis for RT-qPCR based expression analysis.

Genome-wide identification and characterization of the soybean SOD family during alkaline stress

Background

Superoxide dismutase (SOD) proteins, as one kind of the antioxidant enzymes, play critical roles in plant response to various environment stresses. Even though its functions in the oxidative stress were very well characterized, the roles of SOD family genes in regulating alkaline stress response are not fully reported.

Methods

We identified the potential family members by using Hidden Markov model and soybean genome database. The neighbor-joining phylogenetic tree and exon-intron structures were generated by using software MEGA 5.0 and GSDS online server, respectively. Furthermore, the conserved motifs were analyzed by MEME online server. The syntenic analysis was conducted using Circos-0.69. Additionally, the expression levels of soybean SOD genes under alkaline stress were identified by qRT-PCR.

Results

In this study, we identified 13 potential SOD genes in soybean genome. Phylogenetic analysis suggested that SOD genes could be classified into three subfamilies, including MnSODs (GmMSD1-2), FeSODs (GmFSD1-5) and Cu/ZnSODs (GmCSD1-6). We further investigated the gene structure, chromosomal locations and gene-duplication, conserved domains and promoter cis-elements of the soybean SOD genes. We also explored the expression profiles of soybean SOD genes in different tissues and alkaline, salt and cold stresses, based on the transcriptome data. In addition, we detected their expression patterns in roots and leaves by qRT-PCR under alkaline stress, and found that different SOD subfamily genes may play different roles in response to alkaline stress. These results also confirmed the hypothesis that the great evolutionary divergence may contribute to the potential functional diversity in soybean SOD genes. Taken together, we established a foundation for further functional characterization of soybean SOD genes in response to alkaline stress in the future.

Selection and Validation of Appropriate Reference Genes for Quantitative RT-PCR Analysis in Rubia yunnanensis Diels Based on Transcriptome Data

Real-time quantitative polymerase chain reaction (RT-qPCR) has been widely applied in gene expression and transcription abundance analysis because of its high sensitivity, good repeatability, and strong specificity. Selection of relatively stable reference genes is a precondition in order to obtain the reliable analysis results. However, little is known about evaluation of a set of reference genes through scientific experiments in Rubia plants. Here, 15 candidate reference genes were selected from R. yunnanensis transcriptome database and analyzed under abiotic stresses, hormone treatments, and different tissues. Among these 15 candidate reference genes, heterogeneous nuclear ribonucleoprotein (hnRNP), TATA binding protein (TBP), ribosomal protein L5 (RPL5), malate dehydrogenase (MDH), and elongation factor 1-alpha (EF-1α) were indicated as the five most stable reference genes by four statistical programs (geNorm, NormFinder, BestKeeper, and RefFinder). Ultimately, the validity of reference genes was confirmed by normalizing the expression of o-succinylbenzoate-CoA ligase (OSBL) and isochorismate synthase (ICS) involved in the anthraquinone biosynthesis pathway in different tissues and hormone treatments. Meanwhile, four other putative genes involved in the anthraquinone biosynthesis pathway were also normalized with the selected reference genes, which showed similar expression levels with those given by transcriptome data. This work is the first research that aims at a systematic validation on the stability of reference genes selected from R. yunnanensis transcriptome data and will be conducive to analyze gene expression in R. yunnanensis or other Rubia species.

Evaluation of duplicated reference genes for quantitative real-time PCR analysis in genome unknown hexaploid oat (Avena sativa L.)

BACKGROUND

Oat (Avena sativa L.), a hexaploid crop with unknown genome, has valuable nutritional, medicinal and pharmaceutical uses. However, no suitable RGs (reference genes) for qPCR (quantitative real-time PCR) has been documented for oat yet. Single-copy gene is often selected as RG, which is challengeable or impactable in unexplored polyploids.

RESULTS

In this study, eleven candidate RGs, including four duplicated genes, were selected from oat transcriptome. The stability and the optimal combination of these candidate RGs were assessed in 18 oat samples by using four statistical algorithms including the ΔCt method, geNorm, NormFinder and BestKeeper. The most stable RGs for "all samples", "shoots and roots of seedlings", "developing seeds" and "developing endosperms" were EIF4A (Eukaryotic initiation factor 4A-3), UBC21 (Ubiquitin-Conjugating Enzyme 21), EP (Expressed protein) and EIF4A respectively. Among these RGs, UBC21 was a four-copy duplicated gene. The reliability was validated by the expression patterns of four various genes normalized to the most and the least stable RGs in different sample sets.

CONCLUSIONS

Results provide a proof of concept that the duplicated RG is feasible for qPCR in polyploids. To our knowledge, this study is the first systematic research on the optimal RGs for accurate qPCR normalization of gene expression in different organs and tissues of oat.

Selection of valid reference genes for quantitative real-time PCR in Cotesia chilonis (Hymenoptera: Braconidae) exposed to different temperatures

In quantitative real-time PCR (qRT-PCR), data are normalized using reference genes, which helps to control for internal differences and reduce error among samples. In this study, the expression profiles of eight candidate housekeeping genes, 18S ribosomal (18S rRNA), elongation factor (EF1), glyceraldehyde-3-phosphate dehydrogenase (GAPDH), ribosomal protein L10 (RPL10), ribosomal protein L17 (RPL17), histone 3 (H3), arginine kinase (AK), amd β-Actin (ACTB), were evaluated in the parasitic wasp Cotesia chilonis in response to different temperatures. Specifically, the performance and stabilities of these genes were compared in adult wasps maintained in a growth condition at 27°C (normal storage conditions) and in adults obtained from pupae refrigerated at 4°C for five days (cold storage conditions). Data were analyzed using the ΔCt method, BestKeeper, NormFinder, and geNorm. The optimal numbers and stabilities of reference genes varied between the two temperature treatments (27°C and 4°C). In samples stored at normal developmental temperature (27°C), the requirement for normalization in response to low temperature exposures was three genes (18S, H3, AK), whereas normalization in response to high temperature exposures required only two reference genes (GAPDH, ACTB). In samples stored at cold temperature (4°C), for low temperature exposures two reference genes (RPL17, RPL10) were required for standardization, while following high temperature exposures three reference genes (18S, H3, ACTB) were needed. This study strengthens understanding of the selection of reference genes before qRT-PCR analysis in C. chilonis. The reference genes identified here will facilitate further investigations of the biological characteristics of this important parasitoid.

Screening and Evaluation of Stable Reference Genes for Quantitative Real-Time Polymerase Chain Reaction (qRT-PCR) Analysis in Chinese Fir Roots under Water, Phosphorus, and Nitrogen Stresses

Chinese fir (Cunninghamia lanceolata) is an economical important timber species widely planted in southeastern Asia. Decline in yield and productivity during successive rotation is believed to be linked with abiotic stress, such as drought stress and nitrogen (N) and phosphorus (P) starvation. Molecular breeding could be an option to develop tolerant genotypes. For gene expression studies using quantitative real-time reverse transcription-polymerase chain reaction (qRT-PCR), stable reference genes are needed for normalization of gene expression under different experimental conditions. However, there is no internal reference genes identified for Chinese fir under abiotic stresses. Thus, nine internal reference genes based on transcriptome data were selected and analyzed in the root of Chinese fir under drought stress and N and P starvation. Data were analyzed using geNorm, NormFinder, and BestKeeper, to screen and identify the best reference genes. The results showed that the UBQ and GAPDH genes were the two most stable genes under drought stress and the Actin1 and GAPDH were the two most stable genes under P starvation. Further, it was discovered that the Actin1 and UBC were the two most stable genes under N starvation among nine candidate reference genes. The gene expression of drought stress induced expression protein 14-3-3-4, the P transporter gene ClPht1;3, and the nitrate transporter gene NRT1.1 were used to verify the stability of the selected reference genes under drought stress and P and N starvation, respectively, and the results revealed that the screened reference genes were sufficient to normalize expression of the target genes. In conclusion, the results demonstrate that the stability of reference genes was closely related to the external conditions and reference genes applied to the roots of Chinese fir under different abiotic stress treatments were different. Our data will facilitate further studies on stress ecology and gene function analysis in Chinese fir.

Flax rhamnogalacturonan lyases: phylogeny, differential expression and modeling of protein structure

Rhamnogalacturonan lyases (RGLs; EC 4.2.2.23) degrade the rhamnogalacturonan I (RG-I) backbone of pectins present in the plant cell wall. These enzymes belong to polysaccharide lyase family 4, members of which are mainly from plants and plant pathogens. RGLs are investigated, as a rule, as pathogen 'weapons' for plant cell wall degradation and subsequent infection. Despite the presence of genes annotated as RGLs in plant genomes and the presence of substrates for enzyme activity in plant cells, evidence supporting the involvement of this enzyme in certain processes is limited. The differential expression of some RGL genes in flax (Linum usitatissimum L.) tissues, revealed in our previous work, prompted us to carry out a total revision (phylogenetic analysis, analysis of expression and protein structure modeling) of all the sequences of flax predicted as coding for RGLs. Comparison of the expressions of LusRGL in various tissues of flax stem revealed that LusRGLs belong to distinct phylogenetic clades, which correspond to two co-expression groups. One of these groups comprised LusRGL6-A and LusRGL6-B genes and was specifically upregulated in flax fibers during deposition of the tertiary cell wall, which has complex RG-I as a key noncellulosic component. The results of homology modeling and docking demonstrated that the topology of the LusRGL6-A catalytic site allowed binding to the RG-I ligand. These findings lead us to suggest the presence of RGL activity in planta and the involvement of special isoforms of RGLs in the modification of RG-I of the tertiary cell wall in plant fibers.

Combined genome-wide association analysis and transcriptome sequencing to identify candidate genes for flax seed fatty acid metabolism

Flax seeds have a high oil content and are rich in unsaturated fatty acids, which have advantageous effects in preventing chronic diseases, such as cardiovascular diseases. At present, flax seeds are mainly developed for oil. Therefore, it is of practical significance to identify the candidate genes of fatty acid metabolism in flax seeds for breeding flax seeds with high oil content. In the present study, a natural population of flax containing 224 samples planted in 3 different environments was studied. The genome-wide association analysis (GWAS) of seed fatty acid content was conducted based on specific length amplified fragment sequencing (SLAF-seq) data. Transcriptome sequencing (RNA-seq) of samples from 3 different periods (14 d, 21 d and 28 d after anthesis) during seed development of the low oil variety Shuangya 4 and the high oil variety NEW was performed. The candidate genes for seed fatty acid metabolism were identified by combined analysis of these 2 methods. GWAS detected 16 SNP loci significantly associated with seed fatty acid content, and RNA-seq analysis identified 11,802 differentially expressed genes between high and low oil samples. Pathway enrichment analysis revealed that some differentially expressed genes were classified into fatty acid-related pathways. After comparison of these differentially expressed genes with the Kyoto Encyclopedia of Genes and Genomes (KEGG) database, 20 genes homologous to other species were obtained. After analysis, 10 candidate genes were screened by GWAS and RNA-seq screening. Of these 10 genes, qRT-PCR assays using flax seeds in 5 different developmental stages showed that the expression levels of 6 candidate genes were significantly correlated with 5 fatty acid contents in seeds of the high oil variety NEW. Through metabolic pathway analysis found that 6 genes were involved in important fatty acid metabolic pathways, and some of them also have upstream and downstream regulation relations. The present study combined GWAS and RNA-seq methods to identify candidate genes for fatty acid metabolism in flax seeds, which provided reference for screening of candidate genes with complex traits.

Selection of reference genes for flowering pathway analysis in the masting plants, Celmisia lyallii and Chionochloa pallens, under variable environmental conditions

Mast flowering is characterised by mass synchronised flowering at irregular intervals over a wide geographical area. An understanding of the molecular drivers of mast flowering requires expression analysis of key developmentally regulated gene(s). Reverse transcription-quantitative PCR is the gold standard technique used to assess expression of target gene(s) and to validate high-throughput sequencing data. Selection and validation of appropriate reference gene(s), used as normalisation factors in transcript abundance analysis, is an essential step to avoid ambiguous expression results. Eight candidate reference genes were assessed to select the best internal normalisation factors in naturally growing masting plants Chionochloa pallens and Celmisia lyallii. Statistical packages geNorm, Normfinder, BestKeeper, ΔCt and RefFinder were used to determine the expression stability in plants translocated to different altitudes and sampled across the season. GAPDH and PP2a in Celmisia and ExP and THP in Chionochloa were found to be the best pairs of reference genes for normalisation of the gene expression data. Our study revealed environmentally-induced changes in reference gene expression, information that will be utilised as we investigate flowering phenology of masting plants under global climatic change.

Characterization of LuWRKY36, a flax transcription factor promoting secoisolariciresinol biosynthesis in response to Fusarium oxysporum elicitors in Linum usitatissimum L. hairy roots

The involvement of a WRKY transcription factor in the regulation of lignan biosynthesis in flax using a hairy root system is described. Secoisolariciresinol is the main flax lignan synthesized by action of LuPLR1 (pinoresinol-lariciresinol reductase 1). LuPLR1 gene promoter deletion experiments have revealed a promoter region containing W boxes potentially responsible for the response to Fusarium oxysporum. W boxes are bound by WRKY transcription factors that play a role in the response to stress. A candidate WRKY transcription factor, LuWRKY36, was isolated from both abscisic acid and Fusarium elicitor-treated flax cell cDNA libraries. This transcription factors contains two WRKY DNA-binding domains and is a homolog of AtWRKY33. Different approaches confirmed LuWRKY36 binding to a W box located in the LuPLR1 promoter occurring through a unique direct interaction mediated by its N-terminal WRKY domain. Our results propose that the positive regulator action of LuWRKY36 on the LuPLR1 gene regulation and lignan biosynthesis in response to biotic stress is positively mediated by abscisic acid and inhibited by ethylene. Additionally, we demonstrate a differential Fusarium elicitor response in susceptible and resistant flax cultivars, seen as a faster and stronger LuPLR1 gene expression response accompanied with higher secoisolariciresinol accumulation in HR of the resistant cultivar.

The control exerted by ABA on lignan biosynthesis in flax (Linum usitatissimum L.) is modulated by a Ca2+ signal transduction involving the calmodulin-like LuCML15b

The LuPLR1 gene encodes a pinoresinol lariciresinol reductase responsible for the biosynthesis of (+)-secoisolariciresinol, a cancer chemopreventive lignan, highly accumulated in the seedcoat of flax (Linum usitatissimum L.). Abscisic acid (ABA) plays a key role in the regulation of LuPLR1 gene expression and lignan accumulation in both seeds and cell suspensions, which require two cis-acting elements (ABRE and MYB2) for this regulation. Ca²⁺ is a universal secondary messenger involved in a wide range of physiological processes including ABA signaling. Therefore, Ca²⁺ may be involved as a mediator of LuPLR1 gene expression and lignan biosynthesis regulation exerted by ABA. To test the potential implication of Ca²⁺ signaling, a pharmacological approach was conducted using both flax cell suspensions and maturing seed systems coupled with a ß-glucuronidase reporter gene experiment, RT-qPCR analysis, lignan quantification as well as Ca²⁺ fluorescence imaging. Exogenous ABA application results in an increase in the intracellular Ca²⁺ cytosolic concentration, originating mainly from the extracellular medium. Promoter-reporter deletion experiments suggest that the ABRE and MYB2 cis-acting elements of the LuPLR1 gene promoter functioned as Ca²⁺-sensitive sequences involved in the ABA-mediated regulation. The use of specific inhibitors pointed the crucial roles of the Ca²⁺ sensors calmodulin-like proteins and Ca²⁺-dependent protein kinases in this regulation. This regulation appeared conserved in the two different studied systems, i.e. cell suspensions and maturing seeds. A calmodulin-like, LuCML15b, identified from gene network analysis is proposed as a key player involved in this signal transduction since RNAi experiments provided direct evidences of this role. Taken together, these results provide new information on the regulation of plant defense and human health-promoting compounds, which could be used to optimize their production.

Identification and validation of reference genes for real-time quantitative RT-PCR analysis in jute

BACKGROUND

With the availability of genome sequences, gene expression analysis of jute has drawn considerable attention for understanding the regulatory mechanisms of fiber development and improving fiber quality. Gene expression profiles of a target gene can provide valuable clues towards the understanding of its biological function. Reverse transcription quantitative real-time PCR (qRT-PCR) is the best method for targeted gene expression analysis due to its sensitivity and reproducibility. However, calculating relative expression requires reference genes, which must be stable across various biological conditions. For this purposes, 11 prospective genes namely, 28S RNA, ACT7, CYP, EF1A, EF2, ETIF3E, GAPDH, PP2Ac, PTB, UBC2 and UBI1 were evaluated for their potential use as reference genes in jute.

RESULTS

The expression stabilities of eleven prospective genes were analyzed in various jute plant tissues, such as the root, stick, bark, leaf, flower, seed and fiber, as well as under abiotic (waterlogged, drought and salinity) and biotic stress (infestation with Macrophomina phaseolina) conditions with different time points. All 11 genes were variably expressed in different tissues and stress conditions. To find suitable reference genes in different sample sets, a comprehensive approach based on four statistical algorithms such as GeNorm, BestKeeper, NormFinder the ΔCt was used. The PP2Ac and EF2 genes were the most stably expressed across the different tissues. ACT7 and UBC2 were suitable reference genes under drought stress, and CYP and PP2Ac were the most appropriate after inoculation with Macrophomina phaseolina. Under salinity stress, PP2Ac and UBC2 were the best genes, and ACT7 and PP2Ac were the most suitable under waterlogged conditions.

CONCLUSION

Expression stability of reference genes from jute varied in different tissues and selected experimental conditions. Our results provide a valuable resource for the accurate normalization of gene expression experiments in fiber research for important bast fiber crops.

Reference genes for accurate evaluation of expression levels in Trichophyton interdigitale grown under different carbon sources, pH levels and phosphate levels

Tinea pedis is a type of dermatophytosis caused by anthropophilic keratinolytic fungi such as Trichophyton interdigitale. Quantitative reverse transcription PCR (RT-qPCR) is a reliable and reproducible technique for measuring changes in target gene expression across various biological conditions. A crucial aspect of accurate normalization is the choice of appropriate internal controls. To identify reference genes for accurate evaluation of expression levels in T. interdigitale, the transcription levels of eight candidate reference genes (adp-rf, β-act, ef1-α, gapdh, psm1, sdha, rpl2 and ubc) and one target gene (Tri m4) were analysed by RT-qPCR after growing the dermatophyte under different environmental conditions. The results obtained from expression stability evaluations with NormFinder, geNorm, BestKeeper, and RefFinder software demonstrated that adp-rf and psm1 were the most stable internal control genes across all experimental conditions. The present study constitutes the first report of the identification and validation of reference genes for RT-qPCR normalization for T. interdigitale grown under different environmental conditions resembling the conditions encountered by fungi during invasion of skin.

Genomic insights into HSFs as candidate genes for high-temperature stress adaptation and gene editing with minimal off-target effects in flax

Flax (Linum usitatissimum) is a cool season crop commercially cultivated for seed oil and stem fibre production. A comprehensive characterization of the heat shock factor (HSF) candidate genes in flax can accelerate genetic improvement and adaptive breeding for high temperature stress tolerance. We report the genome-wide identification of 34 putative HSF genes from the flax genome, which we mapped on 14 of the 15 chromosomes. Through comparative homology analysis, we classified these genes into three broad groups, and sub-groups. The arrangement of HSF-specific protein motifs, DNA-binding domain (DBD) and hydrophobic heptad repeat (HR-A/B), and exon-intron boundaries substantiated the phylogenetic separation of these genes. Orthologous relationships and evolutionary analysis revealed that the co-evolution of the LusHSF genes was due to recent genome duplication events. Digital and RT-qPCR analyses provided significant evidence of the differential expression of the LusHSF genes in various tissues, at various developmental stages, and in response to high-temperature stress. The co-localization of diverse cis-acting elements in the promoters of the LusHSF genes further emphasized their regulatory roles in the abiotic stress response. We further confirmed DNA-binding sites on the LusHSF proteins and designed guide RNA sequences for gene editing with minimal off-target effects. These results will hasten functional investigations of LusHSFs or assist in devising genome engineering strategies to develop high-temperature stress tolerant flax cultivars.

Selection and Validation of Reference Genes for Gene Expression Studies by RT-PCR in Dalbergia odorifera

Perennial tree Dalbergia odorifera T. Chen could form the precious heartwood used to produce chinese traditional medicine, rosewood furniture and fragrances. However the formation of heartwood is time-consuming and low efficient, leading to the severe destruction of its wild resources. Thus, it is urgent to study the molecular mechanism of heartwood formation in D. odorifera. But till now, there is no report about the reference gene selection in this species. In this study, the expression stability of nine candidate reference genes were evaluated across different tissues and stems treated by wound and chemical stimulators. Four algorithms were applied to obtain the robust genes. The results support HIS2, GAPDH, and CYP to be the most stable reference genes in samples under different wound treatments while DNAj was the least stable. In different tissues, HIS2, UBQ, and RPL were the most stable reference genes while DNAj was the least stable. The selected reference genes were validated through the normalization of the qRT-PCR data of six heartwood related genes in terpene biosynthesis pathway and ethylene signal pathway. The results showed that their expression levels were accurate when they were normalized by the most stable reference gene HIS2, or by the combination of the two or three most stable reference genes. These results demonstrated that these selected reference genes are reliable.