The involvement of wheat F-box protein gene TaFBA1 in the oxidative stress tolerance of plants

Identifying the Quantitative Trait Locus and Candidate Genes of Traits Related to Milling Quality in Rice via a Genome-Wide Association Study

Milling quality directly affects production efficiency in rice, which is closely related to the brown rice recovery (BRR), the milled rice recovery (MRR) and the head milled rice recovery (HMRR). The present study investigated these three traits in 173 germplasms in two environments, finding abundant phenotypic variation. Three QTLs for BRR, two for MRR, and three for HMRR were identified in a genome-wide association study, five of these were identified in previously reported QTLs and three were newly identified. By combining the linkage disequilibrium (LD) analyses, the candidate gene LOC_Os05g08350 was identified. It had two haplotypes with significant differences and Hap 2 increased the BRR by 4.40%. The results of the qRT-PCR showed that the expression of LOC_Os05g08350 in small-BRR accessions was significantly higher than that in large-BRR accessions at Stages 4-5 of young panicle development, reaching the maximum value at Stage 5. The increase in thickness of the spikelet hulls of the accession carrying LOC_Os05g08350TT occurred due to an increase in the cell width and the cell numbers in cross-sections of spikelet hulls. These results help to further clarify the molecular genetic mechanism of milling-quality-related traits and provide genetic germplasm materials for high-quality breeding in rice.

Physiological function and regulation of ascorbate peroxidase isoforms

Ascorbate peroxidase (APX) reduces H2O2 to H2O by utilizing ascorbate as a specific electron donor and constitutes the ascorbate-glutathione cycle in organelles of plants including chloroplasts, cytosol, mitochondria, and peroxisomes. It has been almost 40 years since APX was discovered as an important plant-specific H2O2-scavenging enzyme, during which time many research groups have conducted molecular physiological analyses. It is now clear that APX isoforms function not only just as antioxidant enzymes but also as important factors in intracellular redox regulation through the metabolism of reactive oxygen species. The function of APX isoforms is regulated at multiple steps, from the transcriptional level to post-translational modifications of enzymes, thereby allowing them to respond flexibly to ever-changing environmental factors and physiological phenomena such as cell growth and signal transduction. In this review, we summarize the physiological functions and regulation mechanisms of expression of each APX isoform.

WGCNA Reveals Hub Genes and Key Gene Regulatory Pathways of the Response of Soybean to Infection by Soybean mosaic virus

Soybean mosaic virus (SMV) is one of the main pathogens that can negatively affect soybean production and quality. To study the gene regulatory network of soybeans in response to SMV SC15, the resistant line X149 and susceptible line X97 were subjected to transcriptome analysis at 0, 2, 8, 12, 24, and 48 h post-inoculation (hpi). Differential expression analysis revealed that 10,190 differentially expressed genes (DEGs) responded to SC15 infection. Weighted gene co-expression network analysis (WGCNA) was performed to identify highly related resistance gene modules; in total, eight modules, including 2256 DEGs, were identified. Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analysis of 2256 DEGs revealed that the genes significantly clustered into resistance-related pathways, such as the plant-pathogen interaction pathway, mitogen-activated protein kinases (MAPK) signaling pathway, and plant hormone signal transduction pathway. Among these pathways, we found that the flg22, Ca2+, hydrogen peroxide (H2O2), and abscisic acid (ABA) regulatory pathways were fully covered by 36 DEGs. Among the 36 DEGs, the gene Glyma.01G225100 (protein phosphatase 2C, PP2C) in the ABA regulatory pathway, the gene Glyma.16G031900 (WRKY transcription factor 22, WRKY22) in Ca2+ and H2O2 regulatory pathways, and the gene Glyma.04G175300 (calcium-dependent protein kinase, CDPK) in Ca2+ regulatory pathways were highly connected hub genes. These results indicate that the resistance of X149 to SC15 may depend on the positive regulation of flg22, Ca2+, H2O2, and ABA regulatory pathways. Our study further showed that superoxide dismutase (SOD) activity, H2O2 content, and catalase (CAT) and peroxidase (POD) activities were significantly up-regulated in the resistant line X149 compared with those in 0 hpi. This finding indicates that the H2O2 regulatory pathway might be dependent on flg22- and Ca2+-pathway-induced ROS generation. In addition, two hub genes, Glyma.07G190100 (encoding F-box protein) and Glyma.12G185400 (encoding calmodulin-like proteins, CMLs), were also identified and they could positively regulate X149 resistance. This study provides pathways for further investigation of SMV resistance mechanisms in soybean.

Identification of Gene Responsible for Conferring Resistance against Race KN2 of Podosphaera xanthii in Melon

Powdery mildew caused by Podosphaera xanthii is a serious fungal disease which causes severe damage to melon production. Unlike with chemical fungicides, managing this disease with resistance varieties is cost effective and ecofriendly. But, the occurrence of new races and a breakdown of the existing resistance genes poses a great threat. Therefore, this study aimed to identify the resistance locus responsible for conferring resistance against P. xanthii race KN2 in melon line IML107. A bi-parental F2 population was used in this study to uncover the resistance against race KN2. Genetic analysis revealed the resistance to be monogenic and controlled by a single dominant gene in IML107. Initial marker analysis revealed the position of the gene to be located on chromosome 2 where many of the resistance gene against P. xanthii have been previously reported. Availability of the whole genome of melon and its R gene analysis facilitated the identification of a F-box type Leucine Rich Repeats (LRR) to be accountable for the resistance against race KN2 in IML107. The molecular marker developed in this study can be used for marker assisted breeding programs.

Identification and molecular evolution of the GLX genes in 21 plant species: a focus on the Gossypium hirsutum

BACKGROUND

The glyoxalase system includes glyoxalase I (GLXI), glyoxalase II (GLXII) and glyoxalase III (GLXIII), which are responsible for methylglyoxal (MG) detoxification and involved in abiotic stress responses such as drought, salinity and heavy metal.

RESULTS

In this study, a total of 620 GLX family genes were identified from 21 different plant species. The results of evolutionary analysis showed that GLX genes exist in all species from lower plants to higher plants, inferring that GLX genes might be important for plants, and GLXI and GLXII account for the majority. In addition, motif showed an expanding trend in the process of evolution. The analysis of cis-acting elements in 21 different plant species showed that the promoter region of the GLX genes were rich in phytohormones and biotic and abiotic stress-related elements, indicating that GLX genes can participate in a variety of life processes. In cotton, GLXs could be divided into two groups and most GLXIs distributed in group I, GLXIIs and GLXIIIs mainly belonged to group II, indicating that there are more similarities between GLXII and GLXIII in cotton evolution. The transcriptome data analysis and quantitative real-time PCR analysis (qRT-PCR) show that some members of GLX family would respond to high temperature treatment in G.hirsutum. The protein interaction network of GLXs in G.hirsutum implied that most members can participate in various life processes through protein interactions.

CONCLUSIONS

The results elucidated the evolutionary history of GLX family genes in plants and lay the foundation for their functions analysis in cotton.

Role of F-box E3-ubiquitin ligases in plant development and stress responses

KEY MESSAGE

F-box E3-ubiquitin ligases regulate critical biological processes in plant development and stress responses. Future research could elucidate why and how plants have acquired a large number of F-box genes. The ubiquitin-proteasome system (UPS) is a predominant regulatory mechanism employed by plants to maintain the protein turnover in the cells and involves the interplay of three classes of enzymes, E1 (ubiquitin-activating), E2 (ubiquitin-conjugating), and E3 ligases. The diverse and most prominent protein family among eukaryotes, F-box proteins, are a vital component of the multi-subunit SCF (Skp1-Cullin 1-F-box) complex among E3 ligases. Several F-box proteins with multifarious functions in different plant systems have evolved rapidly over time within closely related species, but only a small part has been characterized. We need to advance our understanding of substrate-recognition regulation and the involvement of F-box proteins in biological processes and environmental adaptation. This review presents a background of E3 ligases with particular emphasis on the F-box proteins, their structural assembly, and their mechanism of action during substrate recognition. We discuss how the F-box proteins regulate and participate in the signaling mechanisms of plant development and environmental responses. We highlight an urgent need for research on the molecular basis of the F-box E3-ubiquitin ligases in plant physiology, systems biology, and biotechnology. Further, the developments and outlooks of the potential technologies targeting the E3-ubiquitin ligases for developing crop improvement strategies have been discussed.

Insights of auxin signaling F-box genes in wheat (Triticum aestivum L.) and their dynamic expression during the leaf rust infection

The TRANSPORT INHIBITOR RESPONSE 1/AUXIN SIGNALING F-BOX (TIR1/AFB) protein serves as auxin receptor and links with Aux/IAA repressor protein leading to its degradation via SKP-Cullin-F box (SCF^TIR1/AFB) complex in the auxin signaling pathway. Present study revealed 11 TIR1/AFB genes in wheat by genome-wide search using AFB HMM profile. Phylogenetic analysis clustered these genes in two classes. Several phytohormone, abiotic, and biotic stress responsive cis-elements were detected in promoter regions of TIR1/AFB genes. These genes were localized on homoeologous chromosome groups 2, 3, and 5 showing orthologous relation with other monocot plants. Most genes were interrupted by introns and the gene products were localized in cytoplasm, nucleus, and cell organelles. TaAFB3, TaAFB5, and TaAFB8 had nuclear localization signals. The evolutionary constraint suggested paralogous sister pairs and orthologous genes went through strong purifying selection process and are slowly evolving at protein level. Functional annotation revealed all TaAFB genes participated in auxin activated signaling pathway and SCF-mediated ubiquitination process. Furthermore, in silico expression study revealed their diverse expression profiles during various developmental stages in different tissues and organs as well as during biotic and abiotic stress. QRT-PCR based studies suggested distinct expression pattern of TIR1-1, TIR1-3, TaAFB1, TaAFB2, TaAFB3, TaAFB4, TaAFB5, TaAFB7, and TaAFB8 displaying maximum expression at 24 and 48 h post inoculation in both susceptible and resistant near isogenic wheat lines infected with leaf rust pathogen. Importantly, this also reflects coordinated responses in expression patterns of wheat TIR1/AFB genes during progression stages of leaf rust infection.

Dynamic Regulation of Grapevine's microRNAs in Response to Mycorrhizal Symbiosis and High Temperature

MicroRNAs (miRNAs) are non-coding small RNAs that play crucial roles in plant development and stress responses and can regulate plant interactions with beneficial soil microorganisms such as arbuscular mycorrhizal fungi (AMF). To determine if root inoculation with distinct AMF species affected miRNA expression in grapevines subjected to high temperatures, RNA-seq was conducted in leaves of grapevines inoculated with either Rhizoglomus irregulare or Funneliformis mosseae and exposed to a high-temperature treatment (HTT) of 40 °C for 4 h per day for one week. Our results showed that mycorrhizal inoculation resulted in a better plant physiological response to HTT. Amongst the 195 identified miRNAs, 83 were considered isomiRs, suggesting that isomiRs can be biologically functional in plants. The number of differentially expressed miRNAs between temperatures was higher in mycorrhizal (28) than in non-inoculated plants (17). Several miR396 family members, which target homeobox-leucine zipper proteins, were only upregulated by HTT in mycorrhizal plants. Predicted targets of HTT-induced miRNAs in mycorrhizal plants queried to STRING DB formed networks for Cox complex, and growth and stress-related transcription factors such as SQUAMOSA promoter-binding-like-proteins, homeobox-leucine zipper proteins and auxin receptors. A further cluster related to DNA polymerase was found in R. irregulare inoculated plants. The results presented herein provide new insights into miRNA regulation in mycorrhizal grapevines under heat stress and can be the basis for functional studies of plant-AMF-stress interactions.

Genome-Wide Analysis and Characterization of Eggplant F-Box Gene Superfamily: Gene Evolution and Expression Analysis under Stress

F-box genes play an important role in plant growth and resistance to abiotic and biotic stresses. To date, systematic analysis of F-box genes and functional annotation in eggplant (Solanum melongena) is still limited. Here, we identified 389 F-box candidate genes in eggplant. The domain study of F-box candidate genes showed that the F-box domain is conserved, whereas the C-terminal domain is diverse. There are 376 SmFBX candidate genes distributed on 12 chromosomes. A collinearity analysis within the eggplant genome suggested that tandem duplication is the dominant form of F-box gene replication in eggplant. The collinearity analysis between eggplant and the three other species (Arabidopsis thaliana, rice and tomato) provides insight into the evolutionary characteristics of F-box candidate genes. In addition, we analyzed the expression of SmFBX candidate genes in different tissues under high temperature and bacterial wilt stress. The results identified several F-box candidate genes that potentially participate in eggplant heat tolerance and bacterial wilt resistance. Moreover, the yeast two-hybrid assay showed that several representative F-box candidate proteins interacted with representative Skp1 proteins. Overexpression of SmFBX131 and SmFBX230 in tobacco increased resistance to bacterial wilt. Overall, these results provide critical insights into the functional analysis of the F-box gene superfamily in eggplant and provide potentially valuable targets for heat and bacterial resistance.

Identification, Classification and Characterization Analysis of FBXL Gene in Cotton

F-box/LR (FBXL), Leucine-rich repeats in F-box proteins, belongs to the Skp1-Cullin1-F-box protein (SCF) E3 ligase family. FBXL genes play important roles in plant growth, such as plant hormones, responses to environmental stress, and floral organ development. Here, a total of 518 FBXL genes were identified and analyzed in six plant species. Phylogenetic analysis showed that AtFBXLs, VvFBXLs, and GrFBXLs were clustered into three subfamilies (Ⅰ-Ⅲ). Based on the composition of the F-box domain and carboxyl-terminal amino acid sequence, FBXL proteins were classified into three types (Type-A/-B/-C). Whole-genome duplication (WGD) along with tandem duplications and segmental contributed to the expansion of this gene family. The result indicates that four cotton species are also divided into three subfamilies. FBXLs in cotton were classified into three clades by phylogenetic and structural analyses. Furthermore, expression analyses indicated that the expression patterns of GhFBXLs in different cotton tissues were different. The highly expressed of GH_A07G2363 in 5-8 mm anthers, indicates that this gene might play a role in the reproductive process, providing candidate genes for future studies on cotton fertility materials. This study provides an original functional opinion and a useful interpretation of the FBXL protein family in cotton.

Elucidating the genetic architecture controlling antioxidant status and ionic balance in barley under salt stress

Association genetic analysis empowered us to identify candidate genes underlying natural variation of morpho-physiological, antioxidants, and grain yield-related traits in barley. Novel intriguing genomic regions were identified and dissected. Salinity stress is one of the abiotic stresses that influence the morpho-physiological, antioxidants, and yield-related traits in crop plants. The plants of a core set of 138 diverse barley accessions were analyzed after exposure to salt stress under field conditions during the reproductive phase. A genome-wide association scan (GWAS) was then conducted using 19,276 single nucleotide polymorphisms (SNPs) to uncover the genetic basis of morpho-physiological and grain-related traits. A wide range of responses to salt stress by the accessions was explored in the current study. GWAS detected 263 significantly associated SNPs with the antioxidants, K⁺/Na⁺ content ratio, and agronomic traits. Five genomic regions harbored interesting putative candidate genes within LD ± 1.2 Mbp. Choromosome 2H harbored many candidate genes associated with the antioxidants ascorbic acid (AsA) and glutathione (GSH), such as superoxide dismutase (SOD), ascorbate peroxidase (APX), and glutathione reductase (GR), under salt stress. Markedly, an A:C SNP at 153,773,211 bp on chromosome 7H is located inside the gene HORVU.MOREX.r3.7HG0676830 (153,772,300-153,774,057 bp) that was annotated as L-gulonolactone oxidase, regulating the natural variation of SOD_S and APX_S. The allelic variation at this SNP reveals a negative selection of accessions carrying the C allele, predominantly found in six-rowed spring landraces originating from Far-, Near-East, and central Asia carrying photoperiod sensitive alleles having lower activity of enzymatic antioxidants. The SNP-trait associations detected in the current study constitute a benchmark for developing molecular selection tools for antioxidant compound selection in barley.

Genome-wide identification and characterization of F-box family proteins in sweet potato and its expression analysis under abiotic stress

In this study, genome-wide characterization of F-box proteins in sweet potato yielded 243 IbFBX genes, unevenly distributed on the 15 chromosomes of sweet potato. Gene duplication analysis suggested segmental duplication as the principal factor influencing the expansive evolution of IbFBX genes in sweet potato. Phylogenetic analysis clustered F-box proteins in sweet potato, Arabidopsis, and rice into six clades (I-VI). Gene structure analysis of the IbFBX genes revealed that most of the genes within the same clade were highly conserved. Expression profiles of IbFBX family genes in 9 different tissues and under stress conditions revealed that the IbFBXs were highly upregulated or downregulated in response to salt and drought stress, suggesting their significant roles in abiotic stress response and adaptation. Knowledge of the diverse functions and expression patterns of IbFBXs presents a solid theoretical basis for annotating the functions of IbFBXs and further facilitate the molecular breeding of sweet potato.

Plant E3 Ligases and Their Role in Abiotic Stress Response

Plants, as sessile organisms, have limited means to cope with environmental changes. Consequently, they have developed complex regulatory systems to ameliorate abiotic stresses im-posed by environmental changes. One such system is the ubiquitin proteasome pathway, which utilizes E3 ligases to target proteins for proteolytic degradation via the 26S proteasome. Plants ex-press a plethora of E3 ligases that are categorized into four major groups depending on their structure. They are involved in many biological and developmental processes in plants, such as DNA repair, photomorphogenesis, phytohormones signaling, and biotic stress. Moreover, many E3 ligase targets are proteins involved in abiotic stress responses, such as salt, drought, heat, and cold. In this review, we will provide a comprehensive overview of E3 ligases and their substrates that have been connected with abiotic stress in order to illustrate the diversity and complexity of how this pathway enables plant survival under stress conditions.

Genome-wide association identifies candidate genes for drought tolerance in coast redwood and giant sequoia

Drought is a major limitation for survival and growth in plants. With more frequent and severe drought episodes occurring due to climate change, it is imperative to understand the genomic and physiological basis of drought tolerance to be able to predict how species will respond in the future. In this study, univariate and multitrait multivariate genome-wide association study methods were used to identify candidate genes in two iconic and ecosystem-dominating species of the western USA, coast redwood and giant sequoia, using 10 drought-related physiological and anatomical traits and genome-wide sequence-capture single nucleotide polymorphisms. Population-level phenotypic variation was found in carbon isotope discrimination, osmotic pressure at full turgor, xylem hydraulic diameter, and total area of transporting fibers in both species. Our study identified new 78 new marker × trait associations in coast redwood and six in giant sequoia, with genes involved in a range of metabolic, stress, and signaling pathways, among other functions. This study contributes to a better understanding of the genomic basis of drought tolerance in long-generation conifers and helps guide current and future conservation efforts in the species.

Distinct GmASMTs are involved in regulating transcription factors and signalling cross-talk across embryo development, biotic, and abiotic stress in soybean

N-Acetylserotonin O-methyltransferase (ASMT) is the final enzyme involved in melatonin biosynthesis. Identifying the expression of ASMT will reveal the regulatory role in the development and stress conditions in soybean. To identify and characterize ASMT in soybean (GmASMT), we employed genome-wide analysis, gene structure, cis-acting elements, gene expression, co-expression network analysis, and enzyme assay. We found seven pairs of segmental and tandem duplication pairs among the 44 identified GmASMTs by genome-wide analysis. Notably, co-expression network analysis reported that distinct GmASMTs are involved in various stress response. For example, GmASMT3, GmASMT44, GmASMT17, and GmASMT7 are involved in embryo development, heat, drought, aphid, and soybean cyst nematode infections, respectively. These distinct networks of GmASMTs were associated with transcription factors (NAC, MYB, WRKY, and ERF), stress signalling, isoflavone and secondary metabolites, calcium, and calmodulin proteins involved in stress regulation. Further, GmASMTs demonstrated auxin-like activities by regulating the genes involved in auxin transporter (WAT1 and NRT1/PTR) and auxin-responsive protein during developmental and biotic stress. The current study identified the key regulatory role of GmASMTs during development and stress. Hence GmASMT could be the primary target in genetic engineering for crop improvement under changing environmental conditions.

Small RNA, Transcriptome and Degradome Analysis of the Transgenerational Heat Stress Response Network in Durum Wheat

Heat stress is a major limiting factor of grain yield and quality in crops. Abiotic stresses have a transgenerational impact and the mechanistic basis is associated with epigenetic regulation. The current study presents the first systematic analysis of the transgenerational effects of post-anthesis heat stress in tetraploid wheat. Leaf physiological traits, harvest components and grain quality traits were characterized under the impact of parental and progeny heat stress. The parental heat stress treatment had a positive influence on the offspring for traits including chlorophyll content, grain weight, grain number and grain total starch content. Integrated sequencing analysis of the small RNAome, mRNA transcriptome and degradome provided the first description of the molecular networks mediating heat stress adaptation under transgenerational influence. The expression profile of 1771 microRNAs (733 being novel) and 66,559 genes was provided, with differentially expressed microRNAs and genes characterized subject to the progeny treatment, parental treatment and tissue-type factors. Gene Ontology and KEGG pathway analysis of stress responsive microRNAs-mRNA modules provided further information on their functional roles in biological processes such as hormone homeostasis, signal transduction and protein stabilization. Our results provide new insights on the molecular basis of transgenerational heat stress adaptation, which can be used for improving thermo-tolerance in breeding.

Identification of Candidate Genes for Root Traits Using Genotype-Phenotype Association Analysis of Near-Isogenic Lines in Hexaploid Wheat (Triticum aestivum L.)

Global wheat (Triticum aestivum L.) production is constrained by different biotic and abiotic stresses, which are increasing with climate change. An improved root system is essential for adaptability and sustainable wheat production. In this study, 10 pairs of near-isogenic lines (NILs)-targeting four genomic regions (GRs) on chromosome arms 4BS, 4BL, 4AS, and 7AL of hexaploid wheat-were used to phenotype root traits in a semi-hydroponic system. Seven of the 10 NIL pairs significantly differed between their isolines for 11 root traits. The NIL pairs targeting qDSI.4B.1 GR varied the most, followed by the NIL pair targeting qDT.4A.1 and QHtscc.ksu-7A GRs. For pairs 5-7 targeting qDT.4A.1 GR, pair 6 significantly differed in the most root traits. Of the 4 NIL pairs targeting qDSI.4B.1 GR, pairs 2 and 4 significantly differed in 3 and 4 root traits, respectively. Pairs 9 and 10 targeting QHtscc.ksu-7A GR significantly differed in 1 and 4 root traits, respectively. Using the wheat 90K Illumina iSelect array, we identified 15 putative candidate genes associated with different root traits in the contrasting isolines, in which two UDP-glycosyltransferase (UGT)-encoding genes, TraesCS4A02G185300 and TraesCS4A02G442700, and a leucine-rich repeat receptor-like protein kinase (LRR-RLK)-encoding gene, TraesCS4A02G330900, also showed important functions for root trait control in other crops. This study characterized, for the first time, that these GRs control root traits in wheat, and identified candidate genes, although the candidate genes will need further confirmation and validation for marker-assisted wheat breeding.

The Landscape of the Genomic Distribution and the Expression of the F-Box Genes Unveil Genome Plasticity in Hexaploid Wheat during Grain Development and in Response to Heat and Drought Stress

FBX proteins are subunits of the SCF complex (Skp1-cullin-FBX) belonging to the E3 ligase family, which is involved in the ubiquitin-proteasome 26S (UPS) pathway responsible for the post-translational protein turnover. By targeting, in a selective manner, key regulatory proteins for ubiquitination and 26S proteasome degradation, FBX proteins play a major role in plant responses to diverse developmental and stress conditions. Although studies on the genomic organization of the FBX gene family in various species have been reported, knowledge related to bread wheat (Triticum aestivum) is scarce and needs to be broadened. Using the latest assembly of the wheat genome, we identified 3670 TaFBX genes distributed non-homogeneously within the three subgenomes (A, B and D) and between the 21 chromosomes, establishing it as one of the richest gene families among plant species. Based on the presence of the five different chromosomal regions previously identified, the present study focused on the genomic distribution of the TaFBX family and the identification of differentially expressed genes during the embryogenesis stages and in response to heat and drought stress. Most of the time, when comparing the expected number of genes (taking into account the formal gene distribution on the entire wheat genome), the TaFBX family harbors a different pattern at the various stratum of observation (subgenome, chromosome, chromosomal regions). We report here that the local gene expansion of the TaFBX family must be the consequence of multiple and complex events, including tandem and small-scale duplications. Regarding the differentially expressed TaFBX genes, while the majority of the genes are localized in the distal chromosomal regions (R1 and R3), differentially expressed genes are more present in the interstitial regions (R2a and R2b) than expected, which could be an indication of the preservation of major genes in those specific chromosomal regions.

Classification and Computational Analysis of Arabidopsis thaliana Sperm Cell-Specific F-Box Protein Gene 3p.AtFBP113

In plants, F-box proteins (FBPs) constitute one of the largest superfamilies of regulatory proteins. Most F-box proteins are shown to be an integral part of SCF complexes, which carry out the degradation of proteins and regulate diverse important biological processes. Anthers and pollen development have a huge importance in crop breeding. Despite the vast diversity of FBPs in Arabidopsis male reproductive organs, their role in anther and pollen development is not much explored. Moreover, a standard nomenclature for naming FBPs is also lacking. Here, we propose a standard nomenclature for naming the FBPs of Arabidopsis thaliana uniformly and carry out a systematic analysis of sperm cell-specific FBP gene, i.e., 3p.AtFBP113 due to its reported high and preferential expression, for detailed functional annotation. The results revealed that 3p.AtFBP113 is located on the small arm of chromosome and encodes 397 amino acid long soluble, stable, and hydrophilic protein with the possibility of localization in various cellular compartments. The presence of the C-terminal F-box associated domain (FBA) with immunoglobulin-like fold anticipated its role in protein binding. Gene ontology based functional annotation and tissue-specific gene co-expression analysis further strengthened its role in protein binding and ubiquitination. Moreover, various potential post/co-translational modifications were anticipated and the predicted tertiary structure also showed the presence of characteristic domains and fold. Thus, the outcomes of the study will be useful in developing a better understating of the function of 3p.AtFBP113 during the process of pollen development, which will be helpful for targeting the gene for manipulation of male fertility that has immense importance in hybrid breeding.

F-Box Family Genes, LTSF1 and LTSF2, Regulate Low-Temperature Stress Tolerance in Pepper (Capsicum chinense)

The F-box proteins belong to a family of regulatory proteins that play key roles in the proteasomal degradation of other proteins. Plant F-box proteins are functionally diverse, and the precise roles of many such proteins in growth and development are not known. Previously, two low-temperature-sensitive F-box protein family genes (LTSF1 and LTSF2) were identified as candidates responsible for the sensitivity to low temperatures in the pepper (Capsicum chinense) cultivar ‘sy-2’. In the present study, we showed that the virus-induced gene silencing of these genes stunted plant growth and caused abnormal leaf development under low-temperature conditions, similar to what was observed in the low-temperature-sensitive ‘sy-2’ line. Protein–protein interaction analyses revealed that the LTSF1 and LTSF2 proteins interacted with S-phase kinase-associated protein 1 (SKP1), part of the Skp, Cullin, F-box-containing (SCF) complex that catalyzes the ubiquitination of proteins for degradation, suggesting a role for LTSF1 and LTSF2 in protein degradation. Furthermore, transgenic Nicotiana benthamiana plants overexpressing the pepper LTSF1 gene showed an increased tolerance to low-temperature stress and a higher expression of the genes encoding antioxidant enzymes. Taken together, these results suggest that the LTSF1 and LTSF2 F-box proteins are a functional component of the SCF complex and may positively regulate low-temperature stress tolerance by activating antioxidant-enzyme activities.

Integrated Analysis of Small RNA, Transcriptome, and Degradome Sequencing Reveals the Water-Deficit and Heat Stress Response Network in Durum Wheat

Water-deficit and heat stress negatively impact crop production. Mechanisms underlying the response of durum wheat to such stresses are not well understood. With the new durum wheat genome assembly, we conducted the first multi-omics analysis with next-generation sequencing, providing a comprehensive description of the durum wheat small RNAome (sRNAome), mRNA transcriptome, and degradome. Single and combined water-deficit and heat stress were applied to stress-tolerant and -sensitive Australian genotypes to study their response at multiple time-points during reproduction. Analysis of 120 sRNA libraries identified 523 microRNAs (miRNAs), of which 55 were novel. Differentially expressed miRNAs (DEMs) were identified that had significantly altered expression subject to stress type, genotype, and time-point. Transcriptome sequencing identified 49,436 genes, with differentially expressed genes (DEGs) linked to processes associated with hormone homeostasis, photosynthesis, and signaling. With the first durum wheat degradome report, over 100,000 transcript target sites were characterized, and new miRNA-mRNA regulatory pairs were discovered. Integrated omics analysis identified key miRNA-mRNA modules (particularly, novel pairs of miRNAs and transcription factors) with antagonistic regulatory patterns subject to different stresses. GO (Gene Ontology) and KEGG (Kyoto Encyclopedia of Genes and Genomes) enrichment analysis revealed significant roles in plant growth and stress adaptation. Our research provides novel and fundamental knowledge, at the whole-genome level, for transcriptional and post-transcriptional stress regulation in durum wheat.

Transcriptome Analysis Reveals Important Candidate Genes Related to Nutrient Reservoir, Carbohydrate Metabolism, and Defence Proteins during Grain Development of Hexaploid Bread Wheat and Its Diploid Progenitors

Wheat grain development after anthesis is an important biological process, in which major components of seeds are synthesised, and these components are further required for germination and seed vigour. We have made a comparative RNA-Seq analysis between hexaploid wheat and its individual diploid progenitors to know the major differentially expressed genes (DEGs) involved during grain development. Two libraries from each species were generated with an average of 55.63, 55.23, 68.13, and 103.81 million reads, resulting in 79.3K, 113.7K, 90.6K, and 121.3K numbers of transcripts in AA, BB, DD, and AABBDD genome species respectively. Number of expressed genes in hexaploid wheat was not proportional to its genome size, but marginally higher than that of its diploid progenitors. However, to capture all the transcripts in hexaploid wheat, sufficiently higher number of reads was required. Functional analysis of DEGs, in all the three comparisons, showed their predominance in three major classes of genes during grain development, i.e., nutrient reservoirs, carbohydrate metabolism, and defence proteins; some of them were subsequently validated through real time quantitative Reverse Transcription Polymerase Chain Reaction (qRT-PCR). Further, developmental stage-specific gene expression showed most of the defence protein genes expressed during initial developmental stages in hexaploid contrary to the diploids at later stages. Genes related to carbohydrates anabolism expressed during early stages, whereas catabolism genes expressed at later stages in all the species. However, no trend was observed in case of different nutrient reservoirs gene expression. This data could be used to study the comparative gene expression among the three diploid species and homeologue-specific expression in hexaploid.

Diverse and dynamic roles of F-box proteins in plant biology

The SCF complex is a widely studied multi-subunit ring E3 ubiquitin ligase that tags targeted proteins with ubiquitin for protein degradation by the ubiquitin 26S-proteasome system (UPS). The UPS is an important system that generally keeps cellular events tightly regulated by purging misfolded or damaged proteins and selectively degrading important regulatory proteins. The specificity of this post-translational regulation is controlled by F-box proteins (FBPs) via selective recognition of a protein-protein interaction motif at the C-terminal domain. Hence, FBPs are pivotal proteins in determining the plant response in multiple scenarios. It is not surprising that the FBP family is one of the largest protein families in the plant kingdom. In this review, the roles of FBPs, specifically in plants, are compiled to provide insights into their involvement in secondary metabolites, plant stresses, phytohormone signalling, plant developmental processes and miRNA biogenesis.

Genome-wide analysis and characterization of F-box gene family in Gossypium hirsutum L

Background

F-box proteins are substrate-recognition components of the Skp1-Rbx1-Cul1-F-box protein (SCF) ubiquitin ligases. By selectively targeting the key regulatory proteins or enzymes for ubiquitination and 26S proteasome mediated degradation, F-box proteins play diverse roles in plant growth/development and in the responses of plants to both environmental and endogenous signals. Studies of F-box proteins from the model plant Arabidopsis and from many additional plant species have demonstrated that they belong to a super gene family, and function across almost all aspects of the plant life cycle. However, systematic exploration of F-box family genes in the important fiber crop cotton (Gossypium hirsutum) has not been previously performed. The genome-wide analysis of the cotton F-box gene family is now possible thanks to the completion of several cotton genome sequencing projects.

Results

In current study, we first conducted a genome-wide investigation of cotton F-box family genes by reference to the published F-box protein sequences from other plant species. 592 F-box protein encoding genes were identified in the Gossypium hirsutume acc.TM-1 genome and, subsequently, we were able to present their gene structures, chromosomal locations, syntenic relationships with their parent species. In addition, duplication modes analysis showed that cotton F-box genes were distributed to 26 chromosomes, with the maximum number of genes being detected on chromosome 5. Although the WGD (whole-genome duplication) mode seems play a dominant role during cotton F-box gene expansion process, other duplication modes including TD (tandem duplication), PD (proximal duplication), and TRD (transposed duplication) also contribute significantly to the evolutionary expansion of cotton F-box genes. Collectively, these bioinformatic analysis suggest possible evolutionary forces underlying F-box gene diversification. Additionally, we also conducted analyses of gene ontology, and expression profiles in silico, allowing identification of F-box gene members potentially involved in hormone signal transduction.

Conclusion

The results of this study provide first insights into the Gossypium hirsutum F-box gene family, which lays the foundation for future studies of functionality, particularly those involving F-box protein family members that play a role in hormone signal transduction.

Wheat F-box Protein TaFBA1 Positively Regulates Plant Drought Tolerance but Negatively Regulates Stomatal Closure

The phytohormone abscisic acid (ABA) regulates plant growth and development, as well as responses to various stresses, such as salt and drought. The wheat TaFBA1 gene, which encodes an F-box protein, was previously identified in our laboratory by homologous cloning. We previously found that TaFBA1 expression was induced by ABA and drought stress. In this study, wild-type (WT), TaFBA1 over-expressing (OEs), TaFBA1 homologous gene mutants, and TaFBA1 recovery (Rs) Arabidopsis plants were used. We found that the germination rate, the cotyledon greening rate, the root length, and the photosynthetic performance of TaFBA1 OE plants were better than those of WT under drought and ABA conditions, but mutant plants showed the opposite trend, and overexpression of TaFBA1 in mutants can recover their phenotype. In addition, TaFBA1 was found to be a negative regulator of ABA-induced stoma movement; mRNA transcription of certain ABA signaling-related genes was lower in TaFBA1 OE plants than in WT plants following ABA treatment. Further, we found that TaFBA1 can interact with RCAR1 (an ABA receptor) and ABI5. BiFC assay showed that TaFBA1 may interact with RCAR1 in the plasma membrane. In addition, accumulation of ROS and MDA in TaFBA1 OE plants was lower than that in the WT plants after ABA and drought treatments. Based on these results, we suggest that TaFBA1-regulated ABA insensitivity may be dependent on regulating ABA-mediated gene expression through interacting with RCAR1 and ABI5. Increased antioxidant competence and decreased ROS accumulation may be an important mechanism that underlies improved drought tolerance in TaFBA1 OE plants.

Enhanced resistance to Verticillium dahliae mediated by an F-box protein GhACIF1 from Gossypium hirsutum

Avr9/Cf-9-INDUCED F-BOX1 (ACIF1) was first identified during screening of Avr9/Cf-9-elicited genes in tobacco. Further analysis revealed that ACIF1 was required for hypersensitive responses triggered by various elicitors in tobacco and tomato, indicating that it may be involved in various disease resistance. Here, we cloned its cotton (Gossypium hirsutum) homolog GhACIF1, which encodes an F-box protein. We show that GhACIF1 interacts with the putative SKP1-like protein, named GhSKP1. Disease resistance assays show that GhACIF1 enhances resistance to Verticillium dahliae in Arabidopsis plants, while silencing of GhACIF1 confers sensitivity to V. dahliae in cotton. Further analysis show that PevD1 elicitor activates hypersensitive and acquired immune response mediated by GhACIF1. Collectively, these results indicate that GhACIF1 contributes to protection against V. dahliae infection.

MicroRNA-guided regulation of heat stress response in wheat

Background

With rising global temperature, understanding plants’ adaptation to heat stress has implications in plant breeding. MicroRNAs (miRNAs) are small, non-coding, regulatory RNAs guiding gene expression at the post-transcriptional level. In this study, small RNAs and the degradome (parallel analysis of RNA ends) of leaf tissues collected from control and heat-stressed wheat plants immediately at the end of the stress period and 1 and 4 days later were analysed.

Results

Sequencing of 24 small RNA libraries produced 55.2 M reads while 404 M reads were obtained from the corresponding 24 PARE libraries. From these, 202 miRNAs were ascertained, of which mature miRNA evidence was obtained for 104 and 36 were found to be differentially expressed after heat stress. The PARE analysis identified 589 transcripts targeted by 84 of the ascertained miRNAs. PARE sequencing validated the targets of the conserved members of miRNA156, miR166 and miR393 families as squamosa promoter-binding-like, homeobox leucine-zipper and transport inhibitor responsive proteins, respectively. Heat stress responsive miRNA targeted superoxide dismutases and an array of homeobox leucine-zipper proteins, F-box proteins and protein kinases. Query of miRNA targets to interactome databases revealed a predominant association of stress responses such as signalling, antioxidant activity and ubiquitination to superoxide dismutases, F-box proteins, pentatricopeptide repeat-containing proteins and mitochondrial transcription termination factor-like proteins.

Conclusion

The interlaced data set generated in this study identified and validated heat stress regulated miRNAs and their target genes associated with thermotolerance. Such accurate identification and validation of miRNAs and their target genes are essential to develop novel regulatory gene-based breeding strategies.

Electronic supplementary material

The online version of this article (10.1186/s12864-019-5799-6) contains supplementary material, which is available to authorized users.

Arabidopsis SKP1-like protein13 (ASK13) positively regulates seed germination and seedling growth under abiotic stress

SKP1 (S-phase kinase-associated protein1) proteins are key members of the SCF (SKP-cullin-F-box protein) E3 ligase complexes that ubiquitinate target proteins and play diverse roles in plant biology. However, in comparison with other members of the SCF complex, knowledge of SKP1-like proteins is very limited in plants. In the present work, we report that Arabidopsis SKP1-like protein13 (ASK13) is differentially regulated in different organs during seed development and germination and is up-regulated in response to abiotic stress. Yeast two-hybrid library screening and subsequent assessment of in vivo interactions through bimolecular fluorescence complementation analysis revealed that ASK13 not only interacts with F-box proteins but also with other proteins that are not components of SCF complexes. Biochemical analysis demonstrated that ASK13 not only exists as a monomer but also as a homo-oligomer or heteromer with other ASK proteins. Functional analysis using ASK13 overexpression and knockdown lines showed that ASK13 positively influences seed germination and seedling growth, particularly under abiotic stress. Taken together, our data strongly suggest that apart from participation to form SCF complexes, ASK13 interacts with several other proteins and is implicated in different cellular processes distinct from protein degradation.

Wheat F-Box Protein Gene TaFBA1 Is Involved in Plant Tolerance to Heat Stress

Adverse environmental conditions, including high temperature, often affect the growth and production of crops worldwide. F-box protein, a core component of the Skp1-Cullin-F-box (SCF) E3 ligase complex, plays an important role in abiotic stress responses. A previously cloned gene from wheat, TaFBA1, encodes a homologous F-box protein. A Yeast two-Hybrid (Y2H) assay showed that TaFBA1 interacted with other SCF proteins. We found that the expression of TaFBA1 could be induced by heat stress (45°C). Overexpression of TaFBA1 enhanced heat stress tolerance in transgenic tobacco, because growth inhibition was reduced and photosynthesis increased as compared with those in the wild type (WT) plants. Furthermore, the accumulation of H₂O₂, O₂^-, and carbonyl protein decreased and cell damage was alleviated in transgenic plants under heat stress, which resulted in less oxidative damage. However, the transgenic plants contained more enzymatic antioxidants after heat stress, which might be related to the regulation of some antioxidant gene expressions. The qRT-PCR analysis showed that the overexpression of TaFBA1 upregulated the expression of genes involved in reactive oxygen species (ROS) scavenging, proline biosynthesis, and abiotic stress responses. We identified the interaction of TaFBA1 with Triticum aestivum stress responsive protein 1 (TaASRP1) by Y2H assay and bimolecular fluorescence complementation (BiFC) assay. The results suggested that TaFBA1 may improve enzymatic antioxidant levels and regulate gene expression by interacting with other proteins, such as TaASRP1, which leads to the enhanced heat stress tolerance seen in the transgenic plants.

A ζ-carotene desaturase gene, IbZDS, increases β-carotene and lutein contents and enhances salt tolerance in transgenic sweetpotato

ζ-Carotene desaturase (ZDS) is one of the key enzymes in carotenoid biosynthesis pathway. However, the ZDS gene has not been applied to carotenoid improvement of plants. Its roles in tolerance to abiotic stresses have not been reported. In this study, the IbZDS gene was isolated from storage roots of sweetpotato (Ipomoea batatas (L.) Lam.) cv. Nongdafu 14. Its overexpression significantly increased β-carotene and lutein contents and enhanced salt tolerance in transgenic sweetpotato (cv. Kokei No. 14) plants. Significant up-regulation of lycopene β-cyclase (β-LCY) and β-carotene hydroxylase (β-CHY) genes and significant down-regulation of lycopene ε-cyclase (ε-LCY) and ε-carotene hydroxylase (ε-CHY) genes were found in the transgenic plants. Abscisic acid (ABA) and proline contents and superoxide dismutase (SOD), catalase (CAT) and peroxidase (POD) activities were significantly increased, whereas malonaldehyde (MDA) content was significantly decreased in the transgenic plants under salt stress. The salt stress-responsive genes encoding pyrroline-5-carboxylate reductase (P5CR), SOD, CAT, ascorbate peroxidase (APX) and POD were found to be significantly up-regulated in the transgenic plants under salt stress. This study indicates that the IbZDS gene has the potential to be applied for improving β-carotene and lutein contents and salt tolerance in sweetpotato and other plants.

MicroRNA-mediated responses to long-term magnesium-deficiency in Citrus sinensis roots revealed by Illumina sequencing

Background

Magnesium (Mg)-deficiency occurs most frequently in strongly acidic, sandy soils. Citrus are grown mainly on acidic and strong acidic soils. Mg-deficiency causes poor fruit quality and low fruit yield in some Citrus orchards. For the first time, we investigated Mg-deficiency-responsive miRNAs in ‘Xuegan’ (Citrus sinensis) roots using Illumina sequencing in order to obtain some miRNAs presumably responsible for Citrus Mg-deficiency tolerance.

Results

We obtained 101 (69) miRNAs with increased (decreased) expression from Mg-starved roots. Our results suggested that the adaptation of Citrus roots to Mg-deficiency was related to the several aspects: (a) inhibiting root respiration and related gene expression via inducing miR158 and miR2919; (b) enhancing antioxidant system by down-regulating related miRNAs (miR780, miR6190, miR1044, miR5261 and miR1151) and the adaptation to low-phosphorus (miR6190); (c) activating transport-related genes by altering the expression of miR6190, miR6485, miR1044, miR5029 and miR3437; (d) elevating protein ubiquitination due to decreased expression levels of miR1044, miR5261, miR1151 and miR5029; (e) maintaining root growth by regulating miR5261, miR6485 and miR158 expression; and (f) triggering DNA repair (transcription regulation) by regulating miR5176 and miR6485 (miR6028, miR6190, miR6485, miR5621, miR160 and miR7708) expression. Mg-deficiency-responsive miRNAs involved in root signal transduction also had functions in Citrus Mg-deficiency tolerance.

Conclusions

We obtained several novel Mg-deficiency-responsive miRNAs (i.e., miR5261, miR158, miR6190, miR6485, miR1151 and miR1044) possibly contributing to Mg-deficiency tolerance. These results revealed some novel clues on the miRNA-mediated adaptation to nutrient deficiencies in higher plants.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-017-3999-5) contains supplementary material, which is available to authorized users.

Concentration-dependent alterations in gene expression induced by cadmium in Solanum lycopersicum

Cadmium (Cd) toxicity in agricultural soil has received significant attention because of its higher transformation in the food chain and toxicity to humans. The aim of the present study was to develop sensitive and specific biomarkers for Cd stress. Therefore, transcriptional analyses were performed to investigate concentration-response characteristics of Cd responsive genes identified from a Solanum lycopersicum microarray. The results showed that the lowest observable adverse effect concentrations (LOAECs) of Cd to S. lycopersicum were 1 mg/kg for seed germination, 8 mg/kg for root dry weight, 8 mg/kg for root elongation, and 8 mg/kg for root morphology. Furthermore, the genes were differentially expressed even at the lowest Cd concentrations (0.5 mg/kg), indicating that the detection of Cd in soil at the molecular level is a highly sensitive method. Cd in soil was positively correlated with the expression of the F-box protein PP2-B15 (r = 0.809, p < 0.01) and zinc transporter 4 (r = 0.643, p < 0.01), indicating that these two genes could be selected as indicators of soil Cd contamination.

Water deficit stress tolerance in chickpea is mediated by the contribution of integrative defence systems in different tissues of the plant

Drought induces heavy yield losses in chickpea (Cicer arietinum L.). Besides understanding the physiological and biochemical parameters contributing to drought tolerance, we need to understand the importance of one tissue in combatting drought stress-induced oxidative stress and influencing the antioxidative defence system in other tissues. The study was conducted to examine the influence of drought stress conditions on the antioxidative defence system and physiology in different tissues such as roots, leaves, nodules, pod walls and seeds at various vegetative and reproductive growth stages in two chickpea cultivars differing in rooting behaviour: ICC4958 (deep rooted) and ILC3279 (shallow rooted). The traits contributing to drought tolerance in ICC4958 were increased root area, decreased leaf area index or increase in root area, decreased leaf area; ILC3279 displayed a decrease in root area and an increase in LAI. The adaptation of ICC4958 was also accompanied by biochemical adjustments, like increases in antioxidative enzymes (superoxide dismutase, peroxidase, ascorbate peroxidase, glutathione reductase and catalase) and nonenzymatic antioxidants (ascorbic acid, proline and stress-induced proteins). However, increases in antioxidant enzymes, nonenzymatic antioxidants and proteins in ILC3279 were lower than in ICC4958. The lower malondialdehyde content and membrane permeability index in ICC4958 might be responsible for reduced damage under drought stress. Increased H2O2 content in ICC4958 was related to enhanced antioxidative defence, emphasising its role as a signalling molecule under stress. This is the first study conducted on drought stress-induced enzymatic and nonenzymatic antioxidative defence systems in underground, aboveground vegetative and reproductive tissues in chickpea cultivars differing in rooting behaviour.

Deep-sequencing transcriptome analysis of field-grown Medicago sativa L. crown buds acclimated to freezing stress

Medicago sativa L. (alfalfa) 'Zhaodong' is an important forage legume that can safely survive in northern China where winter temperatures reach as low as -30 °C. Survival of alfalfa following freezing stress depends on the amount and revival ability of crown buds. In order to investigate the molecular mechanisms of frost tolerance in alfalfa, we used transcriptome sequencing technology and bioinformatics strategies to analyze crown buds of field-grown alfalfa during winter. We statistically identified a total of 5605 differentially expressed genes (DEGs) involved in freezing stress including 1900 upregulated and 3705 downregulated DEGs. We validated 36 candidate DEGs using qPCR to confirm the accuracy of the RNA-seq data. Unlike other recent studies, this study employed alfalfa plants grown in the natural environment. Our results indicate that not only the CBF orthologs but also membrane proteins, hormone signal transduction pathways, and ubiquitin-mediated proteolysis pathways indicate the presence of a special freezing adaptation mechanism in alfalfa. The antioxidant defense system may rapidly confer freezing tolerance to alfalfa. Importantly, biosynthesis of secondary metabolites and phenylalanine metabolism, which is of potential importance in coordinating freezing tolerance with growth and development, were downregulated in subzero temperatures. The adaptive mechanism for frost tolerance is a complex multigenic process that is not well understood. This systematic analysis provided an in-depth view of stress tolerance mechanisms in alfalfa.

The Small-RNA Profiles of Almond (Prunus dulcis Mill.) Reproductive Tissues in Response to Cold Stress

Spring frost is an important environmental stress that threatens the production of Prunus trees. However, little information is available regarding molecular response of these plants to the frost stress. Using high throughput sequencing, this study was conducted to identify differentially expressed miRNAs, both the conserved and the non-conserved ones, in the reproductive tissues of almond tolerant H genotype under cold stress. Analysis of 50 to 58 million raw reads led to identification of 174 unique conserved and 59 novel microRNAs (miRNAs). Differential expression pattern analysis showed that 50 miRNA families were expressed differentially in one or both of almond reproductive tissues (anther and ovary). Out of these 50 miRNA families, 12 and 15 displayed up-regulation and down-regulation, respectively. The distribution of conserved miRNA families indicated that miR482f harbor the highest number of members. Confirmation of miRNAs expression patterns by quantitative real- time PCR (qPCR) was performed in cold tolerant (H genotype) alongside a sensitive variety (Sh12 genotype). Our analysis revealed differential expression for 9 miRNAs in anther and 3 miRNAs in ovary between these two varieties. Target prediction of miRNAs followed by differential expression analysis resulted in identification of 83 target genes, mostly transcription factors. This study comprehensively catalogued expressed miRNAs under different temperatures in two reproductive tissues (anther and ovary). Results of current study and the previous RNA-seq study, which was conducted in the same tissues by our group, provide a unique opportunity to understand the molecular basis of responses of almond to cold stress. The results can also enhance the possibility for gene manipulation to develop cold tolerant plants.