Genome-wide analysis of homeobox gene family in legumes: identification, gene duplication and expression profiling

Genome wide identification and characterization of the amino acid transporter (AAT) genes regulating seed protein content in chickpea (Cicer arietinum L.)

Amino acid transporters (AATs), besides, being a crucial component for nutrient partitioning system are also vital for growth and development of the plants and stress resilience. In order to understand the role of AAT genes in seed quality proteins, a comprehensive analysis of AAT gene family was carried out in chickpea leading to identification of 109 AAT genes, representing 10 subfamilies with random distribution across the chickpea genome. Several important stress responsive cis-regulatory elements like Myb, ABRE, ERE were detected in the promoter region of these CaAAT genes. Most of the genes belonging to the same sub-families shared the intron-exon distribution pattern owing to their conserved nature. Random distribution of these CaAAT genes was observed on plasma membrane, vacuolar membrane, Endoplasmic reticulum and Golgi membranes, which may be associated to distinct biochemical pathways. In total 92 out 109 CaAAT genes arise as result of duplication, among which segmental duplication was more prominent over tandem duplication. As expected, the phylogenetic tree was divided into 2 major clades, and further sub-divided into different sub-families. Among the 109 CaAAT genes, 25 were found to be interacting with 25 miRNAs, many miRNAs like miR156, miR159 and miR164 were interacting only with single AAT genes. Tissues specific expression pattern of many CaAAT genes was observed like CaAAP7 and CaAVT18 in nodules, CaAAP17, CaAVT5 and CaCAT9 in vegetative tissues while CaCAT10 and CaAAP23 in seed related tissues as per the expression analysis. Mature seed transcriptome data revealed that genotypes having high protein content (ICC 8397, ICC 13461) showed low CaAATs expression as compared to the genotypes having low protein content (FG 212, BG 3054). Amino acid profiling of these genotypes revealed a significant difference in amount of essential and non-essential amino acids, probably due to differential expression of CaAATs. Thus, the present study provides insights into the biological role of AAT genes in chickpea, which will facilitate their functional characterization and role in various developmental stages, stress responses and involvement in nutritional quality enhancement.

Overexpression of LaGRAS enhances phosphorus acquisition via increased root growth of phosphorus-deficient white lupin

The GRAS transcription factors play an indispensable role in plant growth and responses to environmental stresses. The GRAS gene family has extensively been explored in various plant species; however, the comprehensive investigation of GRAS genes in white lupin remains insufficient. In this study, bioinformatics analysis of white lupin genome revealed 51 LaGRAS genes distributed into 10 distinct phylogenetic clades. Gene structure analyses revealed that LaGRAS proteins were considerably conserved among the same subfamilies. Notably, 25 segmental duplications and a single tandem duplication showed that segmental duplication was the major driving force for the expansion of GRAS genes in white lupin. Moreover, LaGRAS genes exhibited preferential expression in young cluster root and mature cluster roots and may play key roles in nutrient acquisition, particularly phosphorus (P). To validate this, RT-qPCR analysis of white lupin plants grown under +P (normal P) and -P (P deficiency) conditions elucidated significant differences in the transcript level of GRAS genes. Among them, LaGRAS38 and LaGRAS39 were identified as potential candidates with induced expression in MCR under -P. Additionally, white lupin transgenic hairy root overexpressing OE-LaGRAS38 and OE-LaGRAS39 showed increased root growth, and P concentration in root and leaf compared to those with empty vector control, suggesting their role in P acquisition. We believe this comprehensive analysis of GRAS members in white lupin is a first step in exploring their role in the regulation of root growth, tissue development, and ultimately improving P use efficiency in legume crops under natural environments.

Genome-Wide Characterization and Expression Profiling of HD-Zip Genes in ABA-Mediated Processes in Fragaria vesca

Members of homeodomain-leucine zipper (HD-Zip) transcription factors can play their roles by modulating abscisic acid (ABA) signaling in Arabidopsis. So far, our knowledge of the functions of HD-Zips in woodland strawberries (Fragaria vesca), a model plant for studying ABA-mediated fruit ripening, is limited. Here, we identified a total of 31 HD-Zip genes (FveHDZ1-31) in F. vesca, and classified them into four subfamilies (I to IV). Promoter analyses show that the ABA-responsive element, ABRE, is prevalent in the promoters of subfamily I and II FveHDZs. RT-qPCR results demonstrate that 10 of the 14 investigated FveHDZs were consistently >1.5-fold up-regulated or down-regulated in expression in response to exogenous ABA, dehydration, and ABA-induced senescence in leaves. Five of the six consistently up-regulated genes are from subfamily I and II. Thereinto, FveHDZ4, and 20 also exhibited significantly enhanced expression along with increased ABA content during fruit ripening. In yeast one-hybrid assays, FveHDZ4 proteins could bind the promoter of an ABA signaling gene FvePP2C6. Collectively, our results strongly support that the FveHDZs, particularly those from subfamilies I and II, are involved in the ABA-mediated processes in F. vesca, providing a basis for further functional characterization of the HD-Zips in strawberry and other plants.

Identification of Brachypodium distachyon B3 genes reveals that BdB3-54 regulates primary root growth

B3 is a class of plant-specific transcription factors with important roles in plant development and growth. Here, we identified 69 B3 transcription factors in Brachypodium distachyon that were unevenly distributed across all five chromosomes. The ARF, REM, LAV, and RAV subfamilies were grouped based on sequence characteristics and phylogenetic relationships. The phylogenetically related members in the B3 family shared conserved domains and gene structures. Expression profiles showed that B3 genes were widely expressed in different tissues and varied in response to different abiotic stresses. BdB3-54 protein from the REM subfamily was located in the nucleus by subcellular localization and processed transcriptional activation activity. Overexpression of BdB3-54 in Arabidopsis increased primary root length. Our study provides a basis for further research on the functions of BdB3 genes.

Microevolution, speciation and macroevolution in rhizobia: Genomic mechanisms and selective patterns

Nodule bacteria (rhizobia), N₂-fixing symbionts of leguminous plants, represent an excellent model to study the fundamental issues of evolutionary biology, including the tradeoff between microevolution, speciation, and macroevolution, which remains poorly understood for free-living organisms. Taxonomically, rhizobia are extremely diverse: they are represented by nearly a dozen families of α-proteobacteria (Rhizobiales) and by some β-proteobacteria. Their genomes are composed of core parts, including house-keeping genes (hkg), and of accessory parts, including symbiotically specialized (sym) genes. In multipartite genomes of evolutionary advanced fast-growing species (Rhizobiaceae), sym genes are clustered on extra-chromosomal replicons (megaplasmids, chromids), facilitating gene transfer in plant-associated microbial communities. In this review, we demonstrate that in rhizobia, microevolution and speciation involve different genomic and ecological mechanisms: the first one is based on the diversification of sym genes occurring under the impacts of host-induced natural selection (including its disruptive, frequency-dependent and group forms); the second one-on the diversification of hkgs under the impacts of unknown factors. By contrast, macroevolution represents the polyphyletic origin of super-species taxa, which are dependent on the transfer of sym genes from rhizobia to various soil-borne bacteria. Since the expression of newly acquired sym genes on foreign genomic backgrounds is usually restricted, conversion of resulted recombinants into the novel rhizobia species involves post-transfer genetic changes. They are presumably supported by host-induced selective processes resulting in the sequential derepression of nod genes responsible for nodulation and of nif/fix genes responsible for symbiotic N₂ fixation.

An integrated transcriptome mapping the regulatory network of coding and long non-coding RNAs provides a genomics resource in chickpea

Large-scale transcriptome analysis can provide a systems-level understanding of biological processes. To accelerate functional genomic studies in chickpea, we perform a comprehensive transcriptome analysis to generate full-length transcriptome and expression atlas of protein-coding genes (PCGs) and long non-coding RNAs (lncRNAs) from 32 different tissues/organs via deep sequencing. The high-depth RNA-seq dataset reveal expression dynamics and tissue-specificity along with associated biological functions of PCGs and lncRNAs during development. The coexpression network analysis reveal modules associated with a particular tissue or a set of related tissues. The components of transcriptional regulatory networks (TRNs), including transcription factors, their cognate cis-regulatory motifs, and target PCGs/lncRNAs that determine developmental programs of different tissues/organs, are identified. Several candidate tissue-specific and abiotic stress-responsive transcripts associated with quantitative trait loci that determine important agronomic traits are also identified. These results provide an important resource to advance functional/translational genomic and genetic studies during chickpea development and environmental conditions.

A full-length transcriptome and expression atlas of protein-coding genes and long non-coding RNAs is generated in chickpea. Components of transcriptional regulatory networks and candidate tissue-specific transcripts associated with quantitative trait loci are identified.

Genome-wide identification, in silico characterization and expression analysis of the RNA helicase gene family in chickpea (C. arietinum L.)

The RNA helicases are an important class of enzymes which are known to influence almost every aspect of RNA metabolism. The majority of RNA helicases belong to the SF2 (superfamily 2) superfamily, members of which are further categorized into three separate subfamilies i.e., the DEAD, DEAH and DExD/H-box subfamilies. In chickpea, these RNA helicases have not been characterized until now. A genome-wide analysis across the chickpea genome led to the identification of a total of 150 RNA helicase genes which included 50 DEAD, 33 DEAH and 67 DExD/H-box genes. These were distributed across all the eight chromosomes, with highest number on chromosome 4 (26) and least on chromosome 8 (8). Gene duplication analysis resulted in identification of 15 paralogous gene pairs with Ka/Ks values < 1, indicating towards the genes being under purifying selection during the course of evolution. The promoter regions of the RNA helicase genes were enriched in cis-acting elements like the light and ABA-responsive elements. The drought responsiveness of the genes was analysed by studying the expression profiles of few of these genes, in two different genotypes, the cultivated variety ICC 8261 (kabuli, C. arietinum) and the wild accession ILWC 292 (C. reticulatum), through qRT-PCR. These genotypes were selected based on their drought responsiveness in a field experiment, where it was observed that the percentage (%) reduction in relative water content (RWC) and membrane stability index (MSI) for the drought stressed plants after withholding water for 24 days, over the control or well-watered plants, was least for both the genotypes. The genes CaDEAD50 and CaDExD/H66 were identified as drought-responsive RNA helicase genes in chickpea. The protein encoded by the CaDExD/H66 gene shares a high degree of homology with one of the CLSY (CLASSY) proteins of A. thaliana. We hypothesize that this gene could possibly be involved in regulation of DNA methylation levels in chickpea by regulating siRNA production, in conjunction with other proteins like the Argonaute, RNA dependent RNA polymerases and Dicer-like proteins.

The Roles of BLH Transcription Factors in Plant Development and Environmental Response

Despite recent advancements in plant molecular biology and biotechnology, providing enough, and safe, food for an increasing world population remains a challenge. The research into plant development and environmental adaptability has attracted more and more attention from various countries. The transcription of some genes, regulated by transcript factors (TFs), and their response to biological and abiotic stresses, are activated or inhibited during plant development; examples include, rooting, flowering, fruit ripening, drought, flooding, high temperature, pathogen infection, etc. Therefore, the screening and characterization of transcription factors have increasingly become a hot topic in the field of plant research. BLH/BELL (BEL1-like homeodomain) transcription factors belong to a subfamily of the TALE (three-amino-acid-loop-extension) superfamily and its members are involved in the regulation of many vital biological processes, during plant development and environmental response. This review focuses on the advances in our understanding of the function of BLH/BELL TFs in different plants and their involvement in the development of meristems, flower, fruit, plant morphogenesis, plant cell wall structure, the response to the environment, including light and plant resistance to stress, biosynthesis and signaling of ABA (Abscisic acid), IAA (Indoleacetic acid), GA (Gibberellic Acid) and JA (Jasmonic Acid). We discuss the theoretical basis and potential regulatory models for BLH/BELL TFs' action and provide a comprehensive view of their multiple roles in modulating different aspects of plant development and response to environmental stress and phytohormones. We also present the value of BLHs in the molecular breeding of improved crop varieties and the future research direction of the BLH gene family.

Genome-wide identification and expression analysis of the GRAS gene family in response to drought stress in chickpea (Cicer arietinum L.)

The GRAS (gibberellic acid insensitive, repressor of GAI and scarecrow) transcription factors (TFs) regulate diverse biological processes involved in plant growth and development. These TFs are also known to regulate gene expression in response to various abiotic stress factors like cold, drought, etc. In chickpea one of the most devastating abiotic stress factors is terminal drought. The GRAS TF family has not been characterized in chickpea (Cicer arietinum L.) until now. In this study, we report 46 GRAS TF genes (CaGRAS genes) in the chickpea genome. The CaGRAS proteins were categorized into nine subfamilies based on their phylogenetic relationship with known GRAS members of Arabidopsis and soybean. The PAT subfamily was the largest consisting of ten CaGRAS members whereas the LAS subfamily was the smallest with only one member. Gene duplication analysis revealed that segmental duplication was the primary reason for the expansion of this gene family within the chickpea genome. The gene expression levels of CaGRAS genes were analysed using two different chickpea varieties contrasting for drought tolerance trait, i.e., ICC 4958 (drought tolerant) and ICC 1882 (drought sensitive). On exposure to drought stress, the two chickpea genotypes, exhibited differential drought response, which was quantified and estimated in terms of differences in leaf relative water content (RWC). The well-watered or control plants of the drought tolerant variety were able to maintain a higher leaf RWC by the end of the drought stress period, whereas the control plants of the drought sensitive variety continued to show a decline in leaf RWC. The two genotypes also differed in their root morphologies, under well-watered and drought stress conditions. The gene expression analysis revealed a potential role of PAT, SCR, SCL3 and SHR GRAS members in the regulation of differential response to drought, in the root tissues, for both the genotypes. CaGRAS 12 (SCR) was identified as a drought-responsive GRAS TF gene, which could serve as a potential candidate gene for utilization in developing chickpea varieties with improved drought tolerance. This study demonstrates the drought-responsive expression of CaGRAS genes in chickpea and also describes the morpho-physiological response of chickpea plants to drought stress conditions.

SUPPLEMENTARY INFORMATION

The online version contains supplementary material available at 10.1007/s13205-021-03104-z.

Genome-Wide Identification, Expression, and Interaction Analysis of BEL-Like Homeodomain Gene Family in Peach

BEL1-like homeodomain (BLH) family genes as homeodomain transcription factors are found ubiquitously in plants to play important regulatory roles in reproductive development, morphological development, and stress response. Although BLH proteins have been reported in some species, there is little information about BLH genes in peach. In this study, we identified 11 peach PpBLH genes based on the conserved domain. Phylogenetic analysis suggested that the PpBLH proteins could be divided into five groups, which might be involved in different aspects of morphogenesis. Genomics structure analysis revealed that there were four exons in the PpBLH gene, and the length of the third exon was 61 bp. Chromosomal location analysis showed that the PpBLH genes were not distributed uniformly on six chromosomes. Promoter analysis showed that the promoter sequences of six PpBLH genes contained multiple cis-acting elements for hormones and stress. Six PpBLH genes were cloned by RT-PCR, and PpBLH1, PpBLH4, and PpBLH7 showed different expression patterns in the tested fruits under common temperature and high temperature. Y2H results indicated that PpBLH7 andPpBLH10 interacted with the PpOFP6 protein, and PpBLH1 interacted with the PpOFP1, PpOFP2, PpOFP4, and PpOFP13 proteins. These results provide new insight for further study of PpBLH genes, and construction of regulatory networks of PpBLH proteins in the growth, development, and stress response of peach.

Genome-wide investigation of the ZF-HD gene family in two varieties of alfalfa (Medicago sativa L.) and its expression pattern under alkaline stress

Background

Zinc finger homeodomain (ZHD) protein is a plant-specific transcription factor and a potential regulator of phosphoenolpyruvate carboxylase (PEPCase)-coding genes, and it also participates in plant growth regulation and abiotic stress responses. To study the function of MsZF-HD genes in the alkaline stress response, this paper assessed biological information and performed transcriptome analysis of the MsZF-HD gene family by using the genomes of two different varieties of alfalfa (XinJiangDa Ye and Zhongmu No. 1).

Results

In total, 49 and 11 MsZF-HD genes were identified in the two different varieties respectively, including the alleles of XinJiangDa Ye. According to their phylogenetic relationships, the 60 MsZF-HD genes were divided into 5 ZHD subfamilies and 1 MIF subfamily. A total of 88.3% of MsZF-HD genes do not contain introns and are unevenly distributed among the 6 chromosomes of alfalfa. A collinearity analysis indicated that 26 genes of XinJiangDa Ye have no orthologous genes in Zhongmu No. 1, although these genes (such as ZHD-X1–2, ZHD-X3–2 and ZHD-X4–2) have homologous genes in Arabidopsis thaliana, Medicago truncatula and Glycine max. Through RNA-seq and qRT–PCR verification, it was found that MsZF-HD genes are downregulated to participate in the alkaline stress response.

Conclusion

The results of this study may lay the foundation for the cloning and functional study of MsZF-HD genes and provide a theoretical basis for revealing the difference between XinJiangDa Ye and Zhongmu No. 1 at the genome level.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-022-08309-x.

Gene cluster from plant to microbes: Their role in genome architecture, organism's development, specialized metabolism and drug discovery

Plants and microbes fulfil our daily requirements through different high-value chemicals, e.g., nutraceuticals, pharmaceuticals, cosmetics, and through varieties of fruits, crops, vegetables, and many more. Utmost care would therefore be taken for growth, development and sustainability of these important crops and medicinal plants and microbes. Homeobox genes and HOX clusters and their recently characterized expanded family members, including newly discovered homeobox, WOX gene from medicinal herb, Panax ginseng, significantly contributes in the growth and development of these organisms. On the other hand, secondary metabolites produced through secondary metabolism of plants and microbes are used as organisms defense as well as drugs/drug-like molecules for humans. Both the developmental HOX cluster and the biosynthetic gene-cluster (BGC) for secondary metabolites are organised in organisms genome. Genome mining and genomewide analysis of these clusters will definitely identify and characterize many more important molecules from unexplored plants and microbes and underexplored human microbiota and the evolution studies of these clusters will indicate their source of origin. Although genomics revolution now continues at a pace, till date only few hundred plant genome sequences are available. However, next-generation sequencing (NGS) technology now in market and may be applied even for plants with recalcitrant genomes, eventually may discover genomic potential towards production of secondary metabolites of diverse plants and micro-organisms present in the environment and microbiota. Additionally, the development of tools for genome mining e.g., antiSMASH, plantiSMASH, and more and more computational approaches that predicts hundreds of secondary metabolite BGCs will be discussed.

Genome-Wide Analysis of the Homeobox Gene Family and Identification of Drought-Responsive Members in Populus trichocarpa

Homeobox (HB) genes play critical roles in the regulation of plant morphogenesis, growth and development. Here, we identified a total of 156 PtrHB genes from the Populus trichocarpa genome. According to the topologies and taxonomy of the phylogenetic tree constructed by Arabidopsis thaliana HB members, all PtrHB proteins were divided into six subgroups, namely HD-ZIP, ZF-HD, HB-PHD, TALE, WOX and HB-OTHERS. Multiple alignments of conserved homeodomains (HDs) revealed the conserved loci of each subgroup, while gene structure analysis showed similar exon–intron gene structures, and motif analysis indicated the similarity of motif number and pattern in the same subgroup. Promoter analysis indicated that the promoters of PtrHB genes contain a series of cis-acting regulatory elements involved in responding to various abiotic stresses, indicating that PtrHBs had potential functions in these processes. Collinearity analysis revealed that there are 96 pairs of 127 PtrHB genes mainly distributing on Chromosomes 1, 2, and 5. We analyzed the spatio-temporal expression patterns of PtrHB genes, and the virus-induced gene silencing (VIGS) of PtrHB3 gene resulted in the compromised tolerance of poplar seedlings to mannitol treatment. The bioinformatics on PtrHB family and preliminary exploration of drought-responsive genes can provide support for further study of the family in woody plants, especially in drought-related biological processes. It also provides a direction for developing new varieties of poplar with drought resistance. Overall, our results provided significant information for further functional analysis of PtrHB genes in poplar and demonstrated that PtrHB3 is a dominant gene regulating tolerance to water stress treatment in poplar seedlings.

Genome-wide identification of ZF-HD gene family in Triticum aestivum: Molecular evolution mechanism and function analysis

ZF-HD family genes play important roles in plant growth and development. Studies about the whole genome analysis of ZF-HD gene family have been reported in some plant species. In this study, the whole genome identification and expression profile of the ZF-HD gene family were analyzed for the first time in wheat. A total of 37 TaZF-HD genes were identified and divided into TaMIF and TaZHD subfamilies according to the conserved domain. The phylogeny tree of the TaZF-HD proteins was further divided into six groups based on the phylogenetic relationship. The 37 TaZF-HDs were distributed on 18 of 21 chromosomes, and almost all the genes had no introns. Gene duplication and Ka/Ks analysis showed that the gene family may have experienced powerful purification selection pressure during wheat evolution. The qRT-PCR analysis showed that TaZF-HD genes had significant expression patterns in different biotic stress and abiotic stress. Through subcellular localization experiments, we found that TaZHD6-3B was located in the nucleus, while TaMIF4-5D was located in the cell membrane and nucleus. Our research contributes to a comprehensive understanding of the TaZF-HD family, provides a new perspective for further research on the biological functions of TaZF-HD genes in wheat.

Identification of a Potential Gene for Elevating ω-3 Concentration and Its Efficiency for Improving the ω-6/ω-3 Ratio in Soybean

This present study was to identify a novel candidate gene that contributes to the elevated α-linolenic acid (ALA, ω-3) concentration in PE2166 from mutagenesis of Pungsannamul. Major loci qALA5_1 and qALA5_2 were detected on chromosome 5 of soybean through quantitative trait loci mapping analyses of recombinant inbred lines. With next-generation sequencing of parental lines and Pungsannamul and recombinant analyses, a potential gene, Glyma.05g221500 (HD), controlling elevated ALA concentration was identified. HD is a homeodomain-like transcriptional regulator that may regulate the expression level of microsomal ω-3 fatty acid desaturase (FAD3) genes responsible for the conversion of linoleic acid into ALA in the fatty acid biosynthetic pathway. In addition, we hypothesized that a combination of mutant alleles, HD, and either of microsomal delta-12 fatty acid desaturase 2-1 (FAD2-1) could reduce the ω-6/ω-3 ratio. In populations where HD, FAD2-1A, and FAD2-1B genes were segregated, a combination of a hd allele from PE2166 and either of the variant FAD2-1 alleles was sufficient to reduce the ω-6/ω-3 ratio in seeds.

Development and validation of a HOXB8 gene-based prognostic model and nomogram for colorectal cancer patients

BACKGROUND

At present, colorectal cancer (CRC) is still associated with a high rate of recurrence and distant metastasis with a poor prognosis. HOXB8 gene is related to the tumorigenesis and development in CRC.

AIM

To explore the prognostic value of HOXB8 gene in CRC patients, and provide a novel insight into the monitoring of disease progression and cancer recurrence in patients with high-risk CRC.

METHODS

The mRNA sequencing data of HOXB8 in CRC patients was downloaded from The Cancer Genome Atlas database. Then, we analyzed the relationship between HOXB8 expression and clinicopathologic features in CRC, and performed survival analysis based on HOXB8 expression. Univariate and multivariate Cox regression analyese were performed for identifying prognostic factors for CRC, and then a nomogram was established and evaluated by concordance index, calibration curve, and decision curve analysis (DCA).

RESULTS

HOXB8 mRNA expression was significantly correlated with CRC tumor tissue (P < 0.001), right-side CRC (P < 0.001), T stage (P = 0.024), and M stage (P = 0.0074). Survival analysis showed that overexpression of HOXB8 was associated with a poor progression-free survival (PFS) in CRC patients (P = 0.0019). Univariate and multivariate COX analyses suggested that the expression level of HOXB8 [HR: 1.539 (1.066-2.221), P = 0.021] and TNM stage were independent prognostic factors for PFS of CRC patients. A nomogram was established to predict 3- and 5-year PFS of CRC patients based on four factors including HOXB8 expression and TNM stage. The concordance index was 0.735, suggesting good discrimination; the calibration curve and DCA showed that the nomogram had good predictive power and clinical practicability.

CONCLUSION

The expression of HOXB8 is significantly related to the prognosis of CRC patients, and it has appreciated predictive ability for disease progression and cancer recurrence in CRC patients. HOXB8 could act as a potential biomarker to identify high-risk CRC patients and become a novel therapeutic target and prognostic indicator for CRC.

Genome-wide survey of the amino acid transporter gene family in wheat (Triticum aestivum L.): Identification, expression analysis and response to abiotic stress

Amino acid transporters (AATs), which transport amino acids across cell membranes, play important roles in alleviating plant damage under stresses and in plant growth. To data, little is known about the AAT genes in wheat because of its complex genome. In this study, a total of 296 AAT genes were identified from the latest wheat genome sequence (IWGSC v1.1) and classified into 12 distinct subfamilies based upon their sequence composition and phylogenetic relationship. The expansion of the wheat AAT family was mainly the results of whole-genome duplication (WGD) and tandem events. The unequal expansion of different subfamilies brought new features to TaAATs. TaAATs were highly expressed and exhibited distinct expression patterns in different tissues. On the basis of homology and expression pattern analysis, we identified several wheat AAT family members that may affect grain quality. In addition, TaAAP3, TaATLa2 and TaATLb13 exhibited sustained expression in response to drought and high-temperature stress. These genes are involved in the response of wheat to abiotic stress by regulating the transport and distribution of amino acids. Overall, our results help to understand the complexity of TaAATs and provide a theoretical basis for further identification and utilization of AATs in wheat and other crop species.

Genome-wide identification and expression profiling of HD-ZIP gene family in Medicago truncatula

The homeodomain-leucine zipper (HD-ZIP) transcription factors are important regulators in various developmental processes and responses to environmental stimuli. Currently, little information is available for HD-ZIP gene family in Medicago truncatula. Here we perform a genome-wide analysis of HD-ZIP gene family in M. truncatula. Totally 52 M. truncatula HD-ZIPs (MtHDZs) were identified and classified into four distinctive subfamilies (I to IV). Members clustered in the same subfamily shared similar gene structure and protein motifs. Fifty-one MtHDZs were non-evenly distributed on eight chromosomes. Segmental duplication and purifying selection mainly contributed to the expansion and retention of M. truncatula HD-ZIP gene family. Expression profiling using the publicly available microarray data revealed that MtHDZ genes exhibited distinctive tissue-specific patterns and divergent responses to drought and salt stresses. In addition, the expression profile between each paralogous pair diverged differentially. Our results identified potential targets for the genetic improvement of abiotic stress tolerance in Medicago.

Developing Climate-Resilient Chickpea Involving Physiological and Molecular Approaches With a Focus on Temperature and Drought Stresses

Chickpea is one of the most economically important food legumes, and a significant source of proteins. It is cultivated in more than 50 countries across Asia, Africa, Europe, Australia, North America, and South America. Chickpea production is limited by various abiotic stresses (cold, heat, drought, salt, etc.). Being a winter-season crop in northern south Asia and some parts of the Australia, chickpea faces low-temperature stress (0-15°C) during the reproductive stage that causes substantial loss of flowers, and thus pods, to inhibit its yield potential by 30-40%. The winter-sown chickpea in the Mediterranean, however, faces cold stress at vegetative stage. In late-sown environments, chickpea faces high-temperature stress during reproductive and pod filling stages, causing considerable yield losses. Both the low and the high temperatures reduce pollen viability, pollen germination on the stigma, and pollen tube growth resulting in poor pod set. Chickpea also experiences drought stress at various growth stages; terminal drought, along with heat stress at flowering and seed filling can reduce yields by 40-45%. In southern Australia and northern regions of south Asia, lack of chilling tolerance in cultivars delays flowering and pod set, and the crop is usually exposed to terminal drought. The incidences of temperature extremes (cold and heat) as well as inconsistent rainfall patterns are expected to increase in near future owing to climate change thereby necessitating the development of stress-tolerant and climate-resilient chickpea cultivars having region specific traits, which perform well under drought, heat, and/or low-temperature stress. Different approaches, such as genetic variability, genomic selection, molecular markers involving quantitative trait loci (QTLs), whole genome sequencing, and transcriptomics analysis have been exploited to improve chickpea production in extreme environments. Biotechnological tools have broadened our understanding of genetic basis as well as plants' responses to abiotic stresses in chickpea, and have opened opportunities to develop stress tolerant chickpea.

Genome-Wide DNA Methylation Analysis of Soybean Curled-Cotyledons Mutant and Functional Evaluation of a Homeodomain-Leucine Zipper (HD-Zip) I Gene GmHDZ20

DNA methylation is a major, conserved epigenetic modification that influences many biological processes. Cotyledons are specialized tissues that provide nutrition for seedlings at the early developmental stage. To investigate the patterns of genomic DNA methylation of germinated cotyledons in soybean (Glycine max) and its effect on cotyledon development, we performed a genome-wide comparative analysis of DNA methylation between the soybean curled-cotyledons (cco) mutant, which has abnormal cotyledons, and its corresponding wild type (WT) by whole-genome bisulfite sequencing. The cco mutant was methylated at more sites but at a slightly lower level overall than the WT on the whole-genome level. A total of 46 CG-, 92 CHG-, and 9723 CHH- (H = A, C, or T) differentially methylated genes (DMGs) were identified in cotyledons. Notably, hypomethylated CHH-DMGs were enriched in the gene ontology term "sequence-specific DNA binding transcription factor activity." We selected a DMG encoding a homeodomain-leucine zipper (HD-Zip) I subgroup transcription factor (GmHDZ20) for further functional characterization. GmHDZ20 localized to the nucleus and was highly expressed in leaf and cotyledon tissues. Constitutive expression of GmHDZ20 in Arabidopsis thaliana led to serrated rosette leaves, shorter siliques, and reduced seed number per silique. A yeast two-hybrid assay revealed that GmHDZ20 physically interacted with three proteins associated with multiple aspects of plant growth. Collectively, our results provide a comprehensive study of soybean DNA methylation in normal and aberrant cotyledons, which will be useful for the identification of specific DMGs that participate in cotyledon development, and also provide a foundation for future in-depth functional study of GmHDZ20 in soybean.

Identification of Homeobox Genes Associated with Lignification and Their Expression Patterns in Bamboo Shoots

: Homeobox (HB) genes play critical roles in regulating various aspects of plant growth and development. However, little is known about HB genes in bamboo. In this study, a total of 115 HB genes (PeHB001‒PeHB115) were identified from moso bamboo (Phyllostachys edulis) and grouped into 13 distinct classes (BEL, DDT, HD-ZIP I‒IV, KNOX, NDX, PHD, PINTOX, PLINC, SAWADEE, and WOX) based on the conserved domains and phylogenetic analysis. The number of members in the different classes ranged from 2 to 24, and they usually varied in terms of exon‒intron distribution pattern and length. There were 20 conserved motifs found in 115 PeHBs, with motif 1 being the most common. Gene ontology (GO) analysis showed that PeHBs had diverse molecular functions, with 19 PeHBs being annotated as having xylem development, xylem, and phloem pattern formation functions. Co-expression network analysis showed that 10 of the 19 PeHBs had co-expression correlations, and three members of the KNOX class were hub proteins that interacted with other transcription factors (TFs) such as MYB, bHLH, and OVATE, which were associated with lignin synthesis. Yeast two-hybridization results further proved that PeHB037 (BEL class) interacted with PeHB057 (KNOX class). Transcriptome expression profiling indicated that all PeHBs except PeHB017 were expressed in at least one of the seven tissues of moso bamboo, and 90 PeHBs were expressed in all the tissues. The qRT-PCR results of the 19 PeHBs showed that most of them were upregulated in shoots as the height increased. Moreover, a KNOX binding site was found in the promoters of the key genes involved in lignin synthesis such as Pe4CL, PeC3H, PeCCR, and PeCOMT, which had positive expression correlations with five KNOX genes. Similar results were found in winter bamboo shoots with prolonged storage time, which was consistent with the degree of lignification. These results provide basic data on PeHBs in moso bamboo, which will be helpful for future functional research on PeHBs with positive regulatory roles in the process of lignification.

The grapevine homeobox gene VvHB58 influences seed and fruit development through multiple hormonal signaling pathways

BACKGROUND

The homeobox transcription factor has a diversity of functions during plant growth and development process. Previous transcriptome analyses of seed development in grape hybrids suggested that specific homeodomain transcription factors are involved in seed development in seedless cultivars. However, the molecular mechanism of homeobox gene regulating seed development in grape is rarely reported.

RESULTS

Here, we report that the grapevine VvHB58 gene, encoding a homeodomain-leucine zipper (HD-Zip I) transcription factor, participates in regulating fruit size and seed number. The VvHB58 gene was differentially expressed during seed development between seedless and seeded cultivars. Subcellular localization assays revealed that the VvHB58 protein was located in the nucleus. Transgenic expression of VvHB58 in tomato led to loss of apical dominance, a reduction in fruit pericarp expansion, reduced fruit size and seed number, and larger endosperm cells. Analysis of the cytosine methylation levels within the VvHB58 promoter indicated that the differential expression during seed development between seedless and seeded grapes may be caused by different transcriptional regulatory mechanisms rather than promoter DNA methylation. Measurements of five classic endogenous hormones and expression analysis of hormone-related genes between VvHB58 transgenic and nontransgenic control plants showed that expression of VvHB58 resulted in significant changes in auxin, gibberellin and ethylene signaling pathways. Additionally, several DNA methylation-related genes were expressed differentially during seed development stages in seedless and seeded grapes, suggesting changes in methylation levels during seed development may be associated with seed abortion.

CONCLUSION

VvHB58 has a potential function in regulating fruit and seed development by impacting multiple hormonal pathways. These results expand understanding of homeodomain transcription factors and potential regulatory mechanism of seed development in grapevine, and provided insights into molecular breeding for grapes.

Small Grain and Dwarf 2, encoding an HD-Zip II family transcription factor, regulates plant development by modulating gibberellin biosynthesis in rice

Homeodomain leucine zipper (HD-Zip) proteins are transcription factors that regulate plant development. Bioactive gibberellin (GA) is a key endogenous hormone that participates in plant growth. However, the relationship between HD-Zip genes and modulation of GA biosynthesis in rice remains elusive. Here, we identified a rice mutant, designated as small grain and dwarf 2 (sgd2), which had reduced height and grain size compared with the wild type. Cytological observations indicated that the defective phenotype was mainly due to decreased cell length. Map-based cloning and complementation tests demonstrated that a 9 bp deletion in a homeodomain leucine zipper (HD-Zip) II family transcription factor was responsible for the sgd2 mutant phenotype. Expression of SGD2 was pronounced in developing panicles, and its protein was localized in nucleus. Luciferase reporter system and transactivation assays in yeast suggested that SGD2 functioned as a transcriptional repressor. High performance liquid chromatography assays showed that the endogenous GA₁ level in the sgd2 mutant was dramatically decreased, and exogenous GA₃ recovered the second leaf sheath to normal length. Results of qRT-PCR showed that the expression levels of genes positively regulating GA-biosynthesis were mostly down-regulated in the mutant. Our data identified the role of an HD-Zip transcription factor that affects rice plant development by modulating gibberellin biosynthesis.

Genomic Organization of the B3-Domain Transcription Factor Family in Grapevine (Vitis vinifera L.) and Expression during Seed Development in Seedless and Seeded Cultivars

Members of the plant-specific B3-domain transcription factor family have important and varied functions, especially with respect to vegetative and reproductive growth. Although B3 genes have been studied in many other plants, there is limited information on the genomic organization and expression of B3 genes in grapevine (Vitis vinifera L.). In this study, we identified 50 B3 genes in the grapevine genome and analyzed these genes in terms of chromosomal location and syntenic relationships, intron–exon organization, and promoter cis-element content. We also analyzed the presumed proteins in terms of domain structure and phylogenetic relationships. Based on the results, we classified these genes into five subfamilies. The syntenic relationships suggest that approximately half of the genes resulted from genome duplication, contributing to the expansion of the B3 family in grapevine. The analysis of cis-element composition suggested that most of these genes may function in response to hormones, light, and stress. We also analyzed expression of members of the B3 family in various structures of grapevine plants, including the seed during seed development. Many B3 genes were expressed preferentially in one or more structures of the developed plant, suggesting specific roles in growth and development. Furthermore, several of the genes were expressed differentially in early developing seeds from representative seeded and seedless cultivars, suggesting a role in seed development or abortion. The results of this study provide a foundation for functional analysis of B3 genes and new resources for future molecular breeding of grapevine.

A comprehensive analysis of the B3 superfamily identifies tissue-specific and stress-responsive genes in chickpea (Cicer arietinum L.)

The aim of this study was to provide a comprehensive analysis of the plant-specific B3 domain-containing transcription factors (TFs) in chickpea. Scanning of the chickpea genome resulted in the identification of 51 B3 domain-containing TFs that were located on seven out of eight chickpea chromosomes. Based on the presence of additional domains other than the B3 domain, the candidates were classified into four subfamilies, i.e., ARF (24), REM (19), LAV (6) and RAV (2). Phylogenetic analysis classified them into four groups in which members of the same group had similar intron-exon organization and motif composition. Genome duplication analysis of the candidate B3 genes revealed an event of segmental duplication that was instrumental in the expansion of the B3 gene family. Ka/Ks analysis showed that the B3 gene family was under purifying selection. Further, chickpea B3 genes showed maximum orthology with Medicago followed by soybean and Arabidopsis. Promoter analyses of the B3 genes led to the identification of several tissue-specific and stress-responsive cis-regulatory elements. Expression profiling of the candidate B3 genes using publicly available RNA-seq data of several chickpea tissues indicated their putative role in plant development and abiotic stress response. These findings were further validated by real-time expression analysis. Overall, this study provides a comprehensive analysis of the B3 domain-containing proteins in chickpea that would aid in devising strategies for crop manipulation in chickpea.

A siphonous macroalgal genome suggests convergent functions of homeobox genes in algae and land plants

Genome evolution and development of unicellular, multinucleate macroalgae (siphonous algae) are poorly known, although various multicellular organisms have been studied extensively. To understand macroalgal developmental evolution, we assembled the ∼26 Mb genome of a siphonous green alga, Caulerpa lentillifera, with high contiguity, containing 9,311 protein-coding genes. Molecular phylogeny using 107 nuclear genes indicates that the diversification of the class Ulvophyceae, including C. lentillifera, occurred before the split of the Chlorophyceae and Trebouxiophyceae. Compared with other green algae, the TALE superclass of homeobox genes, which expanded in land plants, shows a series of lineage-specific duplications in this siphonous macroalga. Plant hormone signalling components were also expanded in a lineage-specific manner. Expanded transport regulators, which show spatially different expression, suggest that the structural patterning strategy of a multinucleate cell depends on diversification of nuclear pore proteins. These results not only imply functional convergence of duplicated genes among green plants, but also provide insight into evolutionary roots of green plants. Based on the present results, we propose cellular and molecular mechanisms involved in the structural differentiation in the siphonous alga.

Genome-wide identification of GRF transcription factors in soybean and expression analysis of GmGRF family under shade stress

BACKGROUND

The Growth-regulating factor (GRF) family encodes plant-specific transcription factors which contain two conserved domains, QLQ and WRC. Members of this family play vital roles in plant development and stress response processes. Although GRFs have been identified in various plant species, we still know little about the GRF family in soybean (Glycine max).

RESULTS

In the present study, 22 GmGRFs distributed on 14 chromosomes and one scaffold were identified by searching soybean genome database and were clustered into five subgroups according to their phylogenetic relationships. GmGRFs belonging to the same subgroup shared a similar motif composition and gene structure. Synteny analysis revealed that large-scale duplications played key roles in the expansion of the GmGRF family. Tissue-specific expression data showed that GmGRFs were strongly expressed in growing tissues, including the shoot apical meristems, developing seeds and flowers, indicating that GmGRFs play critical roles in plant growth and development. On the basis of expression analysis of GmGRFs under shade conditions, we found that all GmGRFs responded to shade stress. Most GmGRFs were down-regulated in soybean leaves after shade treatment.

CONCLUSIONS

Taken together, this research systematically analyzed the characterization of the GmGRF family and its primary roles in soybean development and shade stress response. Further studies of the function of the GmGRFs in the growth, development and stress tolerance of soybean, especially under shade stress, will be valuable.

Phaseolus vulgaris genome possesses CAMTA genes, and phavuCAMTA1 contributes to the drought tolerance

Calmodulin-binding transcription activators (CAMTAs) are a family of transcription factors that play an important role in plants' response to the various biotic and abiotic stresses. The common bean (Phaseolus vulgaris L.) is one of the most important crops in the world and plays a pivotal role in sustainable agriculture. To date, the composition of CAMTA genes in genomes of Phaseolus species and their role in resistance to drought stress are not known. In this study, five PhavuCAMTA genes were characterized in common bean genome through bioinformatics analysis, the morphological and biochemical response of 23 Ph.vulgaris genotypes to different levels of drought stress were evaluated and the expression patterns of PhCAMTA1 in the leaf tissues of sensitive and tolerant genotypes were analysed. Gene structure, protein domain organization and phylogenetic analyses showed that the CAMTAs of Phaseolus were structurally similar and clustered into three groups as other plant CAMTAs. Genotypes showeda differential response to drought stress. Thus, the plant height, number of nodes and flower, total chlorophyll and total protein content, and activity of antioxidant enzymes (ascorbate peroxidase and catalase) in plants significantly influenced by genotype and drought stress interaction. Moreover, the resistant and susceptible genotypes were identified according to three-dimensional plots and the expression patterns of PhavuCAMTA1 gene were studied using real-time quantitative polymerase chain reaction. The results of the present study serve as the basis for future functional studies on the Phaseolus CAMTA family.

A HD-ZIP II HOMEBOX transcription factor, PpHB.G7, mediates ethylene biosynthesis during fruit ripening in peach

Homeobox transcription factors belong to a superfamily that has been widely studied in plant growth and development, but little is known regarding their role in fruit development and ripening. Using a genome-wide expression analysis of homeobox (HB) genes and quantitative real-time PCR, a HD-ZIP II member, PpHB.G7, which presented higher levels of expression in ripening fruits than in developing fruits in all of the tested cultivars, was isolated from peach. Transient transformations showed that PpHB.G7 affects ethylene production and the expression of ethylene biosynthesis genes (PpACS1 and PpACO1). Both dual-luciferase and yeast one-hybrid assays confirmed that PpHB.G7 interacts with the promoters of PpACS1 and PpACO1. Thus, PpHB.G7 mediates ethylene biosynthesis by stimulating PpACS1 and PpACO1 activities. Furthermore, we also found that the other eight HB genes were differentially expressed in the developing fruits, with seven of these genes belonging to the HD-ZIP family. These results suggest that the HB genes in the HD-ZIP family play important roles in fruit development and ripening.

Zinc Finger-Homeodomain Transcriptional Factors (ZHDs) in Upland Cotton (Gossypium hirsutum): Genome-Wide Identification and Expression Analysis in Fiber Development

Zinc finger-homeodomain (ZHD) genes encode a family of plant-specific transcription factors that not only participate in the regulation of plant growth and development but also play an important role in the response to abiotic stress. The ZHD gene family has been studied in several model plants, including Solanum lycopersicum, Zea mays, Oryza sativa, and Arabidopsis thaliana. However, a comprehensive study of the genes of the ZHD family and their roles in fiber development and pigmentation in upland cotton has not been completed. To address this gap, we selected a brown fiber cultivar for our study; brown color in cotton is one of the most desired colors in the textile industry. The natural colored fibers require less processing and little dying, thereby eliminating dye costs and chemical residues. Using bioinformatics approaches, we identified 37 GhZHD genes from Gossypium hirsutum and then divided these genes into seven groups based on their phylogeny. The GhZHD genes were mostly conserved in each subfamily with minor variations in motif distribution and gene structure. These genes were largely distributed on 19 of the 26 upland cotton chromosomes. Among the Gossypium genomes, the paralogs and orthologs of the GhZHD genes were identified and further characterized. Furthermore, among the paralogs, we observed that the ZHD family duplications in Gossypium genomes (G. hirsutum, G. arboreum, and G. raimondii) were probably derived from segmental duplication or genome-wide duplication (GWD) events. Through a combination of qRT-PCR and proanthocyanidins (PA) accumulation analyses in brown cotton fibers, we concluded that the candidate genes involved in early fiber development and fiber pigment synthesis include the following: GhZHD29, GhZHD35, GhZHD30, GhZHD31, GhZHD11, GhZHD27, GhZHD18, GhZHD15, GhZHD16, GhZHD22, GhZHD6, GhZHD33, GhZHD13, GhZHD5, and GhZHD23. This study delivers insights into the evolution of the GhZHD genes in brown cotton, serves as a valuable resource for further studies, and identifies the conditions necessary for improving the quality of brown cotton fiber.

Legume genetic resources and transcriptome dynamics under abiotic stress conditions

Grain legumes are an important source of nutrition and income for billions of consumers and farmers around the world. However, the low productivity of new legume varieties, due to the limited genetic diversity available for legume breeding programmes and poor policymaker support, combined with an increasingly unpredictable global climate is resulting in a large gap between current yields and the increasing demand for legumes as food. Hence, there is a need for novel approaches to develop new high-yielding legume cultivars that are able to cope with a range of environmental stressors. Next-generation technologies are providing the tools that could enable the more rapid and cost-effective genomic and transcriptomic studies for most major crops, allowing the identification of key functional and regulatory genes involved in abiotic stress resistance. In this review, we provide an overview of the recent achievements regarding abiotic stress resistance in a wide range of legume crops and highlight the transcriptomic and miRNA approaches that have been used. In addition, we critically evaluate the availability and importance of legume genetic resources with desirable abiotic stress resistance traits.

Current advances in chickpea genomics: applications and future perspectives

Chickpea genomics promises to illuminate our understanding of genome organization, structural variations, evolutionary and domestication-related insights and fundamental biology of legume crops. Unprecedented advancements of next generation sequencing (NGS) technologies have enabled in decoding of multiple chickpea genome sequences and generating huge genomic resources in chickpea both at functional and structural level. This review is aimed to update the current progress of chickpea genomics ranging from high density linkage map development, genome-wide association studies (GWAS), functional genomics resources for various traits, emerging role of abiotic stress responsive coding and non-coding RNAs after the completion of draft chickpea genome sequences. Additionally, the current efforts of whole genome re-sequencing (WGRS) approach of global chickpea germplasm to capture the global genetic diversity existing in the historically released varieties across the world and increasing the resolution of the previously identified candidate gene(s) of breeding importance have been discussed. Thus, the outcomes of these genomics resources will assist in genomics-assisted selection and facilitate breeding of climate-resilient chickpea cultivars for sustainable agriculture.

Genome-wide identification, expansion, and evolution analysis of homeobox genes and their expression profiles during root development in carrot

The homeobox gene family, a large family represented by transcription factors, has been implicated in secondary growth, early embryo patterning, and hormone response pathways in plants. However, reports about the information and evolutionary history of the homeobox gene family in carrot are limited. In the present study, a total of 130 homeobox family genes were identified in the carrot genome. Specific codomain and phylogenetic analyses revealed that the genes were classified into 14 subgroups. Whole genome and proximal duplication participated in the homeobox gene family expansion in carrot. Purifying selection also contributed to the evolution of carrot homeobox genes. In Gene Ontology (GO) analysis, most members of the HD-ZIP III and IV subfamilies were found to have a lipid binding (GO:0008289) term. Most HD-ZIP III and IV genes also harbored a steroidogenic acute regulatory protein-related lipid transfer (START) domain. These results suggested that the HD-ZIP III and IV subfamilies might be related to lipid transfer. Transcriptome and quantitative real-time PCR (RT-qPCR) data indicated that members of the WOX and KNOX subfamilies were likely implicated in carrot root development. Our study provided a useful basis for further studies on the complexity and function of the homeobox gene family in carrot.

Genome-wide identification and expression analyses of the homeobox transcription factor family during ovule development in seedless and seeded grapes

Seedless grapes are of considerable importance for the raisin and table grape industries. Previous transcriptome analyses of seed development in grape revealed that genes encoding homeobox transcription factors were differentially regulated in seedless compared with seeded grape during seed development. In the present study, we identified a total of 73 homeobox-like genes in the grapevine genome and analyzed the genomic content and expression profiles of these genes. Based on domain architecture and phylogenetic analyses grape homeobox genes can be classified into eleven subfamilies. An analysis of the exon-intron structures and conserved motifs provided further insight into the evolutionary relationships between these genes. Evaluation of synteny indicated that segmental and tandem duplications have contributed greatly to the expansion of the grape homeobox gene superfamily. Synteny analysis between the grape and Arabidopsis genomes provided a potential functional relevance for these genes. The tissue-specific expression patterns of homeobox genes suggested roles in both vegetative and reproductive tissues. Expression profiling of these genes during the course of ovule development in seeded and seedless cultivars suggested a potential role in ovule abortion associated with seedlessness. This study will facilitate the functional analysis of these genes and provide new resources for molecular breeding of seedless grapes.

Over-Expression of OsHOX24 Confers Enhanced Susceptibility to Abiotic Stresses in Transgenic Rice via Modulating Stress-Responsive Gene Expression

Homeobox transcription factors play critical roles in plant development and abiotic stress responses. In the present study, we raised rice transgenics over-expressing stress-responsive OsHOX24 gene (rice homeodomain-leucine zipper I sub-family member) and analyzed their response to various abiotic stresses at different stages of development. At the seed germination stage, rice transgenics over-expressing OsHOX24 exhibited enhanced sensitivity to abiotic stress conditions and abscisic acid as compared to wild-type (WT). OsHOX24 over-expression rice seedlings showed reduced root and shoot growth under salinity and desiccation stress (DS) conditions. Various physiological and phenotypic assays confirmed higher susceptibility of rice transgenics toward abiotic stresses as compared to WT at mature and reproductive stages of rice development too. Global gene expression profiling revealed differential regulation of several genes in the transgenic plants under control and DS conditions. Many of these differentially expressed genes were found to be involved in transcriptional regulatory activities, besides carbohydrate, nucleic acid and lipid metabolic processes and response to abiotic stress and hormones. Taken together, our findings highlighted the role of OsHOX24 in regulation of abiotic stress responses via modulating the expression of stress-responsive genes in rice.

Genome-wide development and deployment of informative intron-spanning and intron-length polymorphism markers for genomics-assisted breeding applications in chickpea

The discovery and large-scale genotyping of informative gene-based markers is essential for rapid delineation of genes/QTLs governing stress tolerance and yield component traits in order to drive genetic enhancement in chickpea. A genome-wide 119169 and 110491 ISM (intron-spanning markers) from 23129 desi and 20386 kabuli protein-coding genes and 7454 in silico InDel (insertion-deletion) (1-45-bp)-based ILP (intron-length polymorphism) markers from 3283 genes were developed that were structurally and functionally annotated on eight chromosomes and unanchored scaffolds of chickpea. A much higher amplification efficiency (83%) and intra-specific polymorphic potential (86%) detected by these markers than that of other sequence-based genetic markers among desi and kabuli chickpea accessions was apparent even by a cost-effective agarose gel-based assay. The genome-wide physically mapped 1718 ILP markers assayed a wider level of functional genetic diversity (19-81%) and well-defined phylogenetics among domesticated chickpea accessions. The gene-derived 1424 ILP markers were anchored on a high-density (inter-marker distance: 0.65cM) desi intra-specific genetic linkage map/functional transcript map (ICC 4958×ICC 2263) of chickpea. This reference genetic map identified six major genomic regions harbouring six robust QTLs mapped on five chromosomes, which explained 11-23% seed weight trait variation (7.6-10.5 LOD) in chickpea. The integration of high-resolution QTL mapping with differential expression profiling detected six including one potential serine carboxypeptidase gene with ILP markers (linked tightly to the major seed weight QTLs) exhibiting seed-specific expression as well as pronounced up-regulation especially in seeds of high (ICC 4958) as compared to low (ICC 2263) seed weight mapping parental accessions. The marker information generated in the present study was made publicly accessible through a user-friendly web-resource, "Chickpea ISM-ILP Marker Database". The designing of multiple ISM and ILP markers (2-5 markers/gene) from an individual gene (transcription factor) with numerous aforementioned desirable genetic attributes can widen the user-preference to select suitable primer combination for simultaneous large-scale assaying of functional allelic variation, natural allelic diversity, molecular mapping and expression profiling of genes among chickpea accessions. This will essentially accelerate the identification of functionally relevant molecular tags regulating vital agronomic traits for genomics-assisted crop improvement by optimal resource expenses in chickpea.

Transcriptome analyses reveal genotype- and developmental stage-specific molecular responses to drought and salinity stresses in chickpea

Drought and salinity are the major factors that limit chickpea production worldwide. We performed whole transcriptome analyses of chickpea genotypes to investigate the molecular basis of drought and salinity stress response/adaptation. Phenotypic analyses confirmed the contrasting responses of the chickpea genotypes to drought or salinity stress. RNA-seq of the roots of drought and salinity related genotypes was carried out under control and stress conditions at vegetative and/or reproductive stages. Comparative analysis of the transcriptomes revealed divergent gene expression in the chickpea genotypes at different developmental stages. We identified a total of 4954 and 5545 genes exclusively regulated in drought-tolerant and salinity-tolerant genotypes, respectively. A significant fraction (~47%) of the transcription factor encoding genes showed differential expression under stress. The key enzymes involved in metabolic pathways, such as carbohydrate metabolism, photosynthesis, lipid metabolism, generation of precursor metabolites/energy, protein modification, redox homeostasis and cell wall component biogenesis, were affected by drought and/or salinity stresses. Interestingly, transcript isoforms showed expression specificity across the chickpea genotypes and/or developmental stages as illustrated by the AP2-EREBP family members. Our findings provide insights into the transcriptome dynamics and components of regulatory network associated with drought and salinity stress responses in chickpea.

Characterization of Rice Homeobox Genes, OsHOX22 and OsHOX24, and Over-expression of OsHOX24 in Transgenic Arabidopsis Suggest Their Role in Abiotic Stress Response

Homeobox transcription factors are well known regulators of plant growth and development. In this study, we carried out functional analysis of two candidate stress-responsive HD-ZIP I class homeobox genes from rice, OsHOX22, and OsHOX24. These genes were highly up-regulated under various abiotic stress conditions at different stages of rice development, including seedling, mature and reproductive stages. The transcript levels of these genes were enhanced significantly in the presence of plant hormones, including abscisic acid (ABA), auxin, salicylic acid, and gibberellic acid. The recombinant full-length and truncated homeobox proteins were found to be localized in the nucleus. Electrophoretic mobility shift assay established the binding of these homeobox proteins with specific DNA sequences, AH1 (CAAT(A/T)ATTG) and AH2 (CAAT(C/G)ATTG). Transactivation assays in yeast revealed the transcriptional activation potential of full-length OsHOX22 and OsHOX24 proteins. Homo- and hetero-dimerization capabilities of these proteins have also been demonstrated. Further, we identified putative novel interacting proteins of OsHOX22 and OsHOX24 via yeast-two hybrid analysis. Over-expression of OsHOX24 imparted higher sensitivity to stress hormone, ABA, and abiotic stresses in the transgenic Arabidopsis plants as revealed by various physiological and phenotypic assays. Microarray analysis revealed differential expression of several stress-responsive genes in transgenic lines as compared to wild-type. Many of these genes were found to be involved in transcriptional regulation and various metabolic pathways. Altogether, our results suggest the possible role of OsHOX22/OsHOX24 homeobox proteins as negative regulators in abiotic stress responses.

Genome-Wide Identification, Classification, and Expression Analysis of Amino Acid Transporter Gene Family in Glycine Max

Amino acid transporters (AATs) play important roles in transporting amino acid across cellular membranes and are essential for plant growth and development. To date, the AAT gene family in soybean (Glycine max L.) has not been characterized. In this study, we identified 189 AAT genes from the entire soybean genomic sequence, and classified them into 12 distinct subfamilies based upon their sequence composition and phylogenetic positions. To further investigate the functions of these genes, we analyzed the chromosome distributions, gene structures, duplication patterns, phylogenetic tree, tissue expression patterns of the 189 AAT genes in soybean. We found that a large number of AAT genes in soybean were expanded via gene duplication, 46 and 36 GmAAT genes were WGD/segmental and tandemly duplicated, respectively. Further comprehensive analyses of the expression profiles of GmAAT genes in various stages of vegetative and reproductive development showed that soybean AAT genes exhibited preferential or distinct expression patterns among different tissues. Overall, our study provides a framework for further analysis of the biological functions of AAT genes in either soybean or other crops.