Genome-Wide Identification of the PHD-Finger Family Genes and Their Responses to Environmental Stresses in Oryza sativa L

Genome-Wide Identification and Expression Analysis of the PHD Finger Gene Family in Pea (Pisum sativum)

The plant homeodomain finger (PHD finger) protein, a type of zinc finger protein extensively distributed in eukaryotes, plays diverse roles in regulating plant growth and development. While PHD finger proteins have been identified in various species, their functions remain largely unexplored in pea (Pisum sativum). In this study, we identified 84 members of the PHD finger gene family in pea, which displayed an uneven distribution across seven chromosomes. Through a comprehensive analysis using data from Arabidopsis thaliana and Medicago truncatula, we categorized the PHD finger proteins into 20 subfamilies via phylogenetic tree analysis. Each subfamily exhibited distinct variations in terms of quantity, genetic structure, conserved domains, and physical and chemical properties. Collinearity analysis revealed conserved evolutionary relationships among the PHD finger genes across the three different species. Furthermore, we identified the conserved and important roles of the subfamily M members in anther development. RT-qPCR and in situ hybridization revealed high expression of the pea subfamily M members PsPHD11 and PsPHD16 in microspores and the tapetum layer. In conclusion, this analysis of the PHD finger family in pea provides valuable guidance for future research on the biological roles of PHD finger proteins in pea and other leguminous plants.

Genome-wide identification of the plant homeodomain-finger family in rye and ScPHD5 functions in cold tolerance and flowering time

KEY MESSAGE

111 PHD genes were newly identified in rye genome and ScPHD5's role in regulating cold tolerance and flowering time was suggested. Plant homeodomain (PHD)-finger proteins regulate the physical properties of chromatin and control plant development and stress tolerance. Although rye (Secale cereale L.) is a major winter crop, PHD-finger proteins in rye have not been studied. Here, we identified 111 PHD genes in the rye genome that exhibited diverse gene and protein sequence structures. Phylogenetic tree analysis revealed that PHDs were genetically close in monocots and diverged from those in dicots. Duplication and synteny analyses demonstrated that ScPHDs have undergone several duplications during evolution and that high synteny is conserved among the Triticeae species. Tissue-specific and abiotic stress-responsive gene expression analyses indicated that ScPHDs were highly expressed in spikelets and developing seeds and were responsive to cold and drought stress. One of these genes, ScPHD5, was selected for further functional characterization. ScPHD5 was highly expressed in the spike tissues and was localized in the nuclei of rye protoplasts and tobacco leaves. ScPHD5-overexpressing Brachypodium was more tolerant to freezing stress than wild-type (WT), with increased CBF and COR gene expression. Additionally, these transgenic plants displayed an extremely early flowering phenotype that flowered more than two weeks earlier than the WT, and vernalization genes, rather than photoperiod genes, were increased in the WT. RNA-seq analysis revealed that diverse stress response genes, including HSPs, HSFs, LEAs, and MADS-box genes, were also upregulated in transgenic plants. Our study will help elucidate the roles of PHD genes in plant development and abiotic stress tolerance in rye.

Transcriptome wide identification, characterization and expression analysis of PHD gene family in Crocus sativus

Crocus sativus L., of the Iridaceae family, yields world's most prized spice, saffron. Saffron is well known for its distinctive aroma, odour and colour, which are imputed to the presence of some specific glycosylated apocarotenoids. Even though the main biosynthetic pathway and most of the enzymes leading to apocarotenoid production have been identified, the regulatory mechanisms that govern the developmental stage and tissue specific production of apocarotenoids in Crocus remain comparatively unravelled. Towards this, we report identification, and characterization of plant homeodomain (PHD) finger transcription factor family in Crocus sativus. We also report cloning and characterisation of CstPHD27 from C. sativus. CstPHD27 recorded highest expression in stigma throughout flower development. CstPHD27 exhibited expression pattern which corresponded to the apocarotenoid accumulation in Crocus stigmas. CstPHD27 is nuclear localized and transcriptionally active in yeast Y187 strain. Over-expression of CstPHD27 in Crocus stigmas enhanced apocarotenoid content by upregulating the biosynthetic pathway genes. This report on PHD finger transcription factor family from C. sativus may offer a basis for elucidating role of this gene family in this traditionally and industrially prized crop.

Supplementary Information

The online version contains supplementary material available at 10.1007/s12298-024-01410-3.

PHD finger proteins function in plant development and abiotic stress responses: an overview

The plant homeodomain (PHD) finger with a conserved Cys4-His-Cys3 motif is a common zinc-binding domain, which is widely present in all eukaryotic genomes. The PHD finger is the "reader" domain of methylation marks in histone H3 and plays a role in the regulation of gene expression patterns. Numerous proteins containing the PHD finger have been found in plants. In this review, we summarize the functional studies on PHD finger proteins in plant growth and development and responses to abiotic stresses in recent years. Some PHD finger proteins, such as VIN3, VILs, and Ehd3, are involved in the regulation of flowering time, while some PHD finger proteins participate in the pollen development, for example, MS, TIP3, and MMD1. Furthermore, other PHD finger proteins regulate the plant tolerance to abiotic stresses, including Alfin1, ALs, and AtSIZ1. Research suggests that PHD finger proteins, as an essential transcription regulator family, play critical roles in various plant biological processes, which is helpful in understanding the molecular mechanisms of novel PHD finger proteins to perform specific function.

Characterization of Plant Homeodomain Transcription Factor Genes Involved in Flower Development and Multiple Abiotic Stress Response in Pepper

Plant homeodomain (PHD) transcription factor genes are involved in plant development and in a plant's response to stress. However, there are few reports about this gene family in peppers (Capsicum annuum L.). In this study, the pepper inbred line "Zunla-1" was used as the reference genome, and a total of 43 PHD genes were identified, and systematic analysis was performed to study the chromosomal location, evolutionary relationship, gene structure, domains, and upstream cis-regulatory elements of the CaPHD genes. The fewest CaPHD genes were located on chromosome 4, while the most were on chromosome 3. Genes with similar gene structures and domains were clustered together. Expression analysis showed that the expression of CaPHD genes was quite different in different tissues and in response to various stress treatments. The expression of CaPHD17 was different in the early stage of flower bud development in the near-isogenic cytoplasmic male-sterile inbred and the maintainer inbred lines. It is speculated that this gene is involved in the development of male sterility in pepper. CaPHD37 was significantly upregulated in leaves and roots after heat stress, and it is speculated that CaPHD37 plays an important role in tolerating heat stress in pepper; in addition, CaPHD9, CaPHD10, CaPHD11, CaPHD17, CaPHD19, CaPHD20, and CaPHD43 were not sensitive to abiotic stress or hormonal factors. This study will provide the basis for further research into the function of CaPHD genes in plant development and responses to abiotic stresses and hormones.

PHD-finger family genes in wheat (Triticum aestivum L.): Evolutionary conservatism, functional diversification, and active expression in abiotic stress

Plant homeodomain (PHD) transcription factors (TFs) are a class of proteins with conserved Cys4-His-Cys3 domains that play important roles in plant growth and development and in response to abiotic stresses. Although characterization of PHDs has been performed in plants, little is known about their function in wheat (Triticum aestivum L.), especially under stress conditions. In the present study, 244 TaPHDs were identified in wheat using comparative genomics. We renamed them TaPHD1-244 based on their chromosomal distribution, and almost all PHD proteins were predicted to be located in the nucleus. According to the unrooted neighbor-joining phylogenetic tree, gene structure, and motif analyses, PHD genes were divided into four clades. A total of 149 TaPHD genes were assigned to arise from duplication events. Furthermore, 230 gene pairs came from wheat itself, and 119, 186, 168, 7, 2, and 6 gene pairs came from six other species (Hordeum vulgareto, Zea mays, Oryza sativa, Arabidopsis thaliana, Brassica rapa, and Gossypium raimondii, respectively). A total of 548 interacting protein branches were identified to be involved in the protein interaction network. Tissue-specific expression pattern analysis showed that TaPHDs were highly expressed in the stigma and ovary during flowering, suggesting that the TaPHD gene plays an active role in the reproductive growth of wheat. In addition, the qRT-PCR results further confirmed that these TaPHD genes are involved in the abiotic stress response of wheat. In conclusion, our study provides a theoretical basis for deciphering the molecular functions of TaPHDs, particularly in response to abiotic stress.

Identification of novel PHD-finger genes in pepper by genomic re-annotation and comparative analyses

BACKGROUND

The plant homeodomain (PHD)-finger gene family that belongs to zinc-finger genes, plays an important role in epigenetics by regulating gene expression in eukaryotes. However, inaccurate annotation of PHD-finger genes hinders further downstream comparative, evolutionary, and functional studies.

RESULTS

We performed genome-wide re-annotation in Arabidopsis thaliana (Arabidopsis), Oryza sativa (rice), Capsicum annuum (pepper), Solanum tuberosum (potato), and Solanum lycopersicum (tomato) to better understand the role of PHD-finger genes in these species. Our investigation identified 875 PHD-finger genes, of which 225 (26% of total) were newly identified, including 57 (54%) novel PHD-finger genes in pepper. The PHD-finger genes of the five plant species have various integrated domains that may be responsible for the diversification of structures and functions of these genes. Evolutionary analyses suggest that PHD-finger genes were expanded recently by lineage-specific duplication, especially in pepper and potato, resulting in diverse repertoires of PHD-finger genes among the species. We validated the expression of six newly identified PHD-finger genes in pepper with qRT-PCR. Transcriptome analyses suggest potential functions of PHD-finger genes in response to various abiotic stresses in pepper.

CONCLUSIONS

Our data, including the updated annotation of PHD-finger genes, provide useful information for further evolutionary and functional analyses to better understand the roles of the PHD-finger gene family in pepper.

Rootstock effects on scion gene expression in maritime pine

Pines are the dominant conifers in Mediterranean forests. As long-lived sessile organisms that seasonally have to cope with drought periods, they have developed a variety of adaptive responses. However, during last decades, highly intense and long-lasting drought events could have contributed to decay and mortality of the most susceptible trees. Among conifer species, Pinus pinaster Ait. shows remarkable ability to adapt to different environments. Previous molecular analysis of a full-sib family designed to study drought response led us to find active transcriptional activity of stress-responding genes even without water deprivation in tolerant genotypes. To improve our knowledge about communication between above- and below-ground organs of maritime pine, we have analyzed four graft-type constructions using two siblings as rootstocks and their progenitors, Gal 1056 and Oria 6, as scions. Transcriptomic profiles of needles from both scions were modified by the rootstock they were grafted on. However, the most significant differential gene expression was observed in drought-sensitive Gal 1056, while in drought-tolerant Oria 6, differential gene expression was very much lower. Furthermore, both scions grafted onto drought-tolerant rootstocks showed activation of genes involved in tolerance to abiotic stress, and is most remarkable in Oria 6 grafts where higher accumulation of transcripts involved in phytohormone action, transcriptional regulation, photosynthesis and signaling has been found. Additionally, processes, such as those related to secondary metabolism, were mainly associated with the scion genotype. This study provides pioneering information about rootstock effects on scion gene expression in conifers.

Genome-Wide Association Mapping of Salinity Tolerance at the Seedling Stage in a Panel of Vietnamese Landraces Reveals New Valuable QTLs for Salinity Stress Tolerance Breeding in Rice

Rice tolerance to salinity stress involves diverse and complementary mechanisms, such as the regulation of genome expression, activation of specific ion-transport systems to manage excess sodium at the cell or plant level, and anatomical changes that avoid sodium penetration into the inner tissues of the plant. These complementary mechanisms can act synergistically to improve salinity tolerance in the plant, which is then interesting in breeding programs to pyramidize complementary QTLs (quantitative trait loci), to improve salinity stress tolerance of the plant at different developmental stages and in different environments. This approach presupposes the identification of salinity tolerance QTLs associated with different mechanisms involved in salinity tolerance, which requires the greatest possible genetic diversity to be explored. To contribute to this goal, we screened an original panel of 179 Vietnamese rice landraces genotyped with 21,623 SNP markers for salinity stress tolerance under 100 mM NaCl treatment, at the seedling stage, with the aim of identifying new QTLs involved in the salinity stress tolerance via a genome-wide association study (GWAS). Nine salinity tolerance-related traits, including the salt injury score, chlorophyll and water content, and K⁺ and Na⁺ contents were measured in leaves. GWAS analysis allowed the identification of 26 QTLs. Interestingly, ten of them were associated with several different traits, which indicates that these QTLs act pleiotropically to control the different levels of plant responses to salinity stress. Twenty-one identified QTLs colocalized with known QTLs. Several genes within these QTLs have functions related to salinity stress tolerance and are mainly involved in gene regulation, signal transduction or hormone signaling. Our study provides promising QTLs for breeding programs to enhance salinity tolerance and identifies candidate genes that should be further functionally studied to better understand salinity tolerance mechanisms in rice.

Identification of conserved genes linked to responses to abiotic stresses in leaves among different plant species

As a consequence of global climate change, certain stress factors that have a negative impact on crop productivity such as heat, cold, drought and salinity are becoming increasingly prevalent. We conducted a meta-analysis to identify genes conserved across plant species involved in (1) general abiotic stress conditions, and (2) specific and unique abiotic stress factors (drought, salinity, extreme temperature) in leaf tissues. We collected raw data and re-analysed eight RNA-Seq studies using our previously published bioinformatic pipeline. A total of 68 samples were analysed. Gene set enrichment analysis was performed using MapMan and PageMan whereas DAVID (Database for Annotation, Visualisation and Integrated Discovery) was used for metabolic process enrichment analysis. We identified of a total of 5122 differentially expressed genes when considering all abiotic stresses (3895 were upregulated and 1227 were downregulated). Jasmonate-related genes were more commonly upregulated by drought, whereas gibberellin downregulation was a key signal for drought and heat. In contrast, cold stress clearly upregulated genes involved in ABA (abscisic acid), cytokinin and gibberellins. A gene (non-phototrophic hypocotyl) involved in IAA (indoleacetic acid) response was induced by heat. Regarding secondary metabolism, as expected, MVA pathway (mevalonate pathway), terpenoids and alkaloids were generally upregulated by all different stresses. However, flavonoids, lignin and lignans were more repressed by heat (cinnamoyl coA reductase 1 and isopentenyl pyrophosphatase). Cold stress drastically modulated genes involved in terpenoid and alkaloids. Relating to transcription factors, AP2-EREBP, MADS-box, WRKY22, MYB, homoebox genes members were significantly modulated by drought stress whereas cold stress enhanced AP2-EREBPs, bZIP members, MYB7, BELL 1 and one bHLH member. C2C2-CO-LIKE, MADS-box and a homeobox (HOMEOBOX3) were mostly repressed in response to heat. Gene set enrichment analysis showed that ubiquitin-mediated protein degradation was enhanced by heat, which unexpectedly repressed glutaredoxin genes. Cold stress mostly upregulated MAP kinases (mitogen-activated protein kinase). Findings of this work will allow the identification of new molecular markers conserved across crops linked to major genes involved in quantitative agronomic traits affected by different abiotic stress.

Genome-wide analysis of PHD finger gene family and identification of potential miRNA and their PHD finger gene specific targets in Oryza sativa indica

Rice (Oryza sativa L.) is one of the most important cereal crops for one third of the world population. However, the grain quality as well as yield of rice is severely affected by various abiotic stresses. Environmental stresses affect the expression of various microRNAs (miRNAs) which in turn negatively regulate gene expression at the post-transcriptional level either by degrading the target mRNA genes or suppressing translation in plants. Plant homeo-domain (PHD) finger proteins are known to be involved in the plant response to salinity stress. In the present study, we identified 44 putative OsPHD finger genes in Oryza sativa Indica, using Ensembl Plants Database. Using computational approach, potential miRNAs that target OsPHD finger genes were identified. Out of the 44 OsPHD finger genes only three OsPHD finger genes i.e., OsPHD2, OsPHD35 and OsPHD11, were found to be targeted by five newly identified putative miRNAs i.e., ath-miRf10010-akr, ath-miRf10110-akr, osa-miR1857–3p, osa-miRf10863-akr, and osa-miRf11806-akr. This is the first report of these five identified miRNAs on targeting PHD finger in Oryza sativa Indica. Further, expression analysis of 44 PHD finger genes under salinity was also performed using quantitative Real-Time PCR. The expression profile of 8 genes were found to be differentially regulated, among them two genes were significantly up regulated i.e., OsPHD6 and OsPHD12. In silico protein-protein interaction analysis using STRING database showed interaction of the OsPHD finger proteins with other protein partners that are directly or indirectly involved in development and abiotic stress tolerance.

•

Total of 44 Plant homeo-domain (PHD) finger proteins were identified & classified into 10 groups in Oryza sativa Indica.

•

This is the first report showing 5 newly identified putative miRNAs targeting three OsPHD genes i.e., OsPHD2, 11 and 35.

•

Expression analysis of PHD finger genes showed up-regulation of the 2 genes OsPHD 6 & 12 under salinity stress treatment.

•

Protein-protein network of OsPHDs showed protein partners that are involved in plant growth and abiotic stress tolerance.

Genome-Wide Identification and Characterization of Drought Stress Responsive microRNAs in Tibetan Wild Barley

Drought stress is a major obstacle to agricultural production. Tibetan wild barley with rich genetic diversity is useful for drought-tolerant improvement of cereals. MicroRNAs (miRNAs) play critical roles in controlling gene expression in response to various environment perturbations in plants. However, the genome-wide expression profiles of miRNAs and their targets in response to drought stress are largely unknown in wild barley. In this study, a polyethylene glycol (PEG) induced drought stress hydroponic experiment was performed, and the expression profiles of miRNAs from the roots of two contrasting Tibetan wild barley genotypes XZ5 (drought-tolerant) and XZ54 (drought-sensitive), and one cultivated barley Tadmor (drought-tolerant) generated by high-throughput sequencing were compared. There were 69 conserved miRNAs and 1574 novel miRNAs in the dataset of three genotypes under control and drought conditions. Among them, seven conserved miRNAs and 36 novel miRNAs showed significantly genotype-specific expression patterns in response to drought stress. And 12 miRNAs were further regarded as drought tolerant associated miRNAs in XZ5, which mostly participate in gene expression, metabolism, signaling and transportation, suggesting that they and their target genes play important roles in plant drought tolerance. This is the first comparation study on the miRNA transcriptome in the roots of two Tibetan wild barley genotypes differing in drought tolerance and one drought tolerant cultivar in response to PEG treatment. Further results revealed the candidate drought tolerant miRNAs and target genes in the miRNA regulation mechanism in wild barley under drought stress. Our findings provide valuable understandings for the functional characterization of miRNAs in drought tolerance.

Ectopic expression of Medicago truncatula homeodomain finger protein, MtPHD6, enhances drought tolerance in Arabidopsis

BACKGROUND

The plant homeodomain (PHD) finger is a Cys4HisCys3-type zinc finger which promotes protein-protein interactions and binds to the cis-acting elements in the promoter regions of target genes. In Medicago truncatula, five PHD homologues with full-length sequence were identified. However, the detailed function of PHD genes was not fully addressed.

RESULTS

In this study, we characterized the function of MtPHD6 during plant responses to drought stress. MtPHD6 was highly induced by drought stress. Ectopic expression of MtPHD6 in Arabidopsis enhanced tolerance to osmotic and drought stresses. MtPHD6 transgenic plants exhibited decreased water loss rate, MDA and ROS contents, and increased leaf water content and antioxidant enzyme activities under drought condition. Global transcriptomic analysis revealed that MtPHD6 reprogramed transcriptional networks in transgenic plants. Expression levels of ABA receptor PYR/PYLs, ZINC FINGER, AP2/EREBP and WRKY transcription factors were mainly up-regulated after transformation of MtPHD6. Interaction network analysis showed that ZINC FINGER, AP2/EREBP and WRKY interacted with each other and downstream stress induced proteins.

CONCLUSIONS

We proposed that ZINC FINGER, AP2/EREBP and WRKY transcription factors were activated through ABA dependent and independent pathways to increase drought tolerance of MtPHD6 transgenic plants.

Genome wide association study of 5 agronomic traits in olive (Olea europaea L.)

Olive (Olea europaea L.) is one of the most economically and historically important fruit crops worldwide. Genetic progress for valuable agronomic traits has been slow in olive despite its importance and benefits. Advances in next generation sequencing technologies provide inexpensive and highly reproducible genotyping approaches such as Genotyping by Sequencing, enabling genome wide association study (GWAS). Here we present the first comprehensive GWAS study on olive using GBS. A total of 183 accessions (FULL panel) were genotyped using GBS, 94 from the Turkish Olive GenBank Resource (TOGR panel) and 89 from the USDA-ARS National Clonal Germplasm Repository (NCGR panel) in the USA. After filtering low quality and redundant markers, GWAS was conducted using 24,977 SNPs in FULL, TOGR and NCGR panels. In total, 52 significant associations were detected for leaf length, fruit weight, stone weight and fruit flesh to pit ratio using the MLM_K. Significant GWAS hits were mapped to their positions and 19 candidate genes were identified within a 10-kb distance of the most significant SNP. Our findings provide a framework for the development of markers and identification of candidate genes that could be used in olive breeding programs.

Genome wide survey, evolution and expression analysis of PHD finger genes reveal their diverse roles during the development and abiotic stress responses in Brassica rapa L

BACKGROUND

Plant homeodomain (PHD) finger proteins are widely present in all eukaryotes and play important roles in chromatin remodeling and transcriptional regulation. The PHD finger can specifically bind a number of histone modifications as an "epigenome reader", and mediate the activation or repression of underlying genes. Many PHD finger genes have been characterized in animals, but only few studies were conducted on plant PHD finger genes to this day. Brassica rapa (AA, 2n = 20) is an economically important vegetal, oilseed and fodder crop, and also a good model crop for functional and evolutionary studies of important gene families among Brassica species due to its close relationship to Arabidopsis thaliana.

RESULTS

We identified a total of 145 putative PHD finger proteins containing 233 PHD domains from the current version of B. rapa genome database. Gene ontology analysis showed that 67.7% of them were predicted to be located in nucleus, and 91.3% were predicted to be involved in protein binding activity. Phylogenetic, gene structure, and additional domain analyses clustered them into different groups and subgroups, reflecting their diverse functional roles during plant growth and development. Chromosomal location analysis showed that they were unevenly distributed on the 10 B. rapa chromosomes. Expression analysis from RNA-Seq data showed that 55.7% of them were constitutively expressed in all the tested tissues or organs with relatively higher expression levels reflecting their important housekeeping roles in plant growth and development, while several other members were identified as preferentially expressed in specific tissues or organs. Expression analysis of a subset of 18 B. rapa PHD finger genes under drought and salt stresses showed that all these tested members were responsive to the two abiotic stress treatments.

CONCLUSIONS

Our results reveal that the PHD finger genes play diverse roles in plant growth and development, and can serve as a source of candidate genes for genetic engineering and improvement of Brassica crops against abiotic stresses. This study provides valuable information and lays the foundation for further functional determination of PHD finger genes across the Brassica species.

Loci and candidate genes in soybean that confer resistance to Fusarium graminearum

KEY MESSAGE

Association analysis techniques were used to identify and verify twelve single nucleotide polymorphisms (SNPs) associated with Fusarium graminearum resistance. Two novel candidate genes were obtained. Fusarium graminearum causes seed and root rot and seedling damping-off of soybean, leading to severe yield loss. Presently, the genetic basis of resistance to F. graminearum is elucidated in only four soybean accessions, which is not sufficient for resistance improvement. The objective of the present study was to identify the genome-wide genetic architecture of resistance to F. graminearum in landraces and cultivated soybeans based on a growth room evaluation. The resistance levels of 314 diverse accessions were tested, and 22,888 single nucleotide polymorphisms (SNPs) with a minor allele frequency of > 0.05 were developed using the specific-locus amplified fragment sequencing (SLAF-seq) approach. Twelve SNPs were identified as associated with F. graminearum resistance, and these SNPs were located at 12 genomic regions on eight chromosomes (Chr.) and could explain 5.53-14.71% of the observed phenotypic variation. One SNP, rs9479021, located on Chr.6, overlapped with qRfg_Gm06, the known QTL for resistance to F. graminearum. The other SNPs were novel and associated with resistance to F. graminearum. Nine novel candidate genes were predicted to contribute to resistance to F. graminearum according to the haplotype and transcript abundance analysis of the candidate genes. The identified markers and resistant cultivars are valuable for the improvement of resistance to F. graminearum.