Genome-wide comparative analysis of the IQD gene families in Arabidopsis thaliana and Oryza sativa

Quantitative Proteomic Analysis of Brassica Napus Reveals Intersections Between Nutrient Deficiency Responses

Macronutrients such as nitrogen (N), phosphorus (P), potassium (K) and sulphur (S) are critical for plant growth and development. Field-grown canola (Brassica napus L.) is supplemented with fertilizers to maximize plant productivity, while deficiency in these nutrients can cause significant yield loss. A holistic understanding of the interplay between these nutrient deficiency responses in a single study and canola cultivar is thus far lacking, hindering efforts to increase the nutrient use efficiency of this important oil seed crop. To address this, we performed a comparative quantitative proteomic analysis of both shoot and root tissue harvested from soil-grown canola plants experiencing either nitrogen, phosphorus, potassium or sulphur deficiency. Our data provide critically needed insights into the shared and distinct molecular responses to macronutrient deficiencies in canola. Importantly, we find more conserved responses to the four different nutrient deficiencies in canola roots, with more distinct proteome changes in aboveground tissue. Our results establish a foundation for a more comprehensive understanding of the shared and distinct nutrient deficiency response mechanisms of canola plants and pave the way for future breeding efforts.

Genome-wide association mapping of bruchid resistance loci in soybean

Soybean is a globally important industrial, food, and cash crop. Despite its importance in present and future economies, its production is severely hampered by bruchids (Callosobruchus chinensis), a destructive storage insect pest, causing considerable yield losses. Therefore, the identification of genomic regions and candidate genes associated with bruchid resistance in soybean is crucial as it helps breeders to develop new soybean varieties with improved resistance and quality. In this study, 6 multi-locus methods of the mrMLM model for genome-wide association study were used to dissect the genetic architecture of bruchid resistance on 4traits: percentage adult bruchid emergence (PBE), percentage weight loss (PWL), median development period (MDP), and Dobie susceptibility index (DSI) on 100 diverse soybean genotypes, genotyped with 14,469 single-nucleotide polymorphism (SNP) markers. Using the best linear unbiased predictors (BLUPs), 13 quantitative trait nucleotides (QTNs) were identified by the mrMLM model, of which rs16_14976250 was associated with more than 1 bruchid resistance traits. As a result, the identified QTNs linked with resistance traits can be employed in marker-assisted breeding for the accurate and rapid screening of soybean genotypes for resistance to bruchids. Moreover, a gene search on the Phytozome soybean reference genome identified 27 potential candidate genes located within a window of 478.45 kb upstream and downstream of the most reliable QTNs. These candidate genes exhibit molecular and biological functionalities associated with various soybean resistance mechanisms and, therefore, could be incorporated into the farmers' preferred soybean varieties that are susceptible to bruchids.

Comprehensive Analysis of the GiTCP Gene Family and Its Expression Under UV-B Radiation in Glycyrrhiza inflata Bat

TCP is a plant-specific transcription factor that plays an important role in plant growth and development. In this study, we used bioinformatics to identify the entire genome of the TCP gene family in Glycyrrhiza inflata Bat, and we analyzed the expression characteristics of GiTCP genes under UV-B radiation using qRT-PCR. The results were as follows: (1) 24 members of the TCP gene family were identified in G. inflata, evenly distributed on its 24 chromosomes. (2) The GiTCP genes contained 0-4 introns and 0-5 exons. (3) The GiTCP genes were phylogenetically divided into three subfamilies-PCF, CIN, and CYC/TB1, with 14, 9, and 1 GiTCP proteins, respectively. (4) A covariance analysis showed that two pairs of GiTCP genes underwent a fragmentary duplication event. (5) A cis-element analysis showed that the cis-responsive elements of the GiTCP genes' promoter regions were mainly comprised of light-responsive, stress-responsive, hormone-regulated, growth and development, and metabolic-regulated elements. (6) A protein network interaction analysis revealed a total of 14 functional molecules of TCPs and 8 potential interacting proteins directly related to GiTCP proteins. (7) GO annotation showed that the GiTCP genes were mainly enriched in BP, CC, and MF groups and had corresponding functions. (8) RNA-seq and qRT-PCR further indicated that GiTCP3, 6, 7, 8, 12, 14, 17, 23, and 24 were up- or down-regulated in G. inflata after UV-B radiation, demonstrating that these genes responded to UV-B radiation in G. inflata. (9) Subcellular localization analysis showed that the GiTCP8 protein was localized in the nucleus. The results of this study provide a basis for further exploration of the function of the GiTCP gene family in the growth and development of G. inflata.

Transcriptome-Based Gene Modules and Soluble Sugar Content Analyses Reveal the Defense Response of Cotton Leaves to Verticillium dahliae

Verticillium dahliae is a soil-borne phytopathogenic fungus causing destructive Verticillium wilt disease that greatly threats cotton production worldwide. The mechanism of cotton resistance to Verticillium wilt is very complex and requires further research. In this study, RNA-sequencing was used to investigate the defense responses of cotton leaves using varieties resistant (Zhongzhimian 2, or Z2) or susceptible (Xinluzao 7, or X7) to V. dahliae. The leaf samples were collected at 48 and 72 hpi (hours post infection) from the two varieties infected by V. dahliae (strain Vd991) or treated by water. Compared to X7, Z2 had less genes responsive to V. dahliae infection at 72 hpi and had no DEGs (differentially expressed genes) at 48 hpi. WGCNA (Weighted Gene Co-Expression Network Analysis) revealed seven key gene modules which were responsible for the resistance of Z2 and susceptibility of X7. KEGG (Kyoto Encyclopedia of Genes and Genomes) analysis of these modules found that several reported disease resistance pathways were found to be up-regulated in Z2, with some of those pathways down-regulated in X7. Unexpectedly, several photosynthesis-related pathways were significantly up-regulated in the leaves of X7 infected by V. dahliae, leading to different profiles of glucose content, which was significantly decreased at 72 hpi and 48 hpi in X7 and Z2, respectively. These results suggest that the leaves of resistant varieties have a slower and different response to V. dahliae compared to those of the susceptible variety, as well as that the translocation of sugars produced by photosynthesis in cotton leaves might vary between the two varieties. Additionally, several HUB genes regulating disease response were identified, including NDR1/HIN1-like protein 12, DELLA protein, cytochrome P450 family protein and LRR receptor-like serine/threonine-protein kinase genes, which have been reported to be related to disease resistance in other plants, which might serve as potential candidates for breeding cotton disease resistance.

Evolution of plant metabolism: the state-of-the-art

Immense chemical diversity is one of the hallmark features of plants. This chemo-diversity is mainly underpinned by a highly complex and biodiverse biochemical machinery. Plant metabolic enzymes originated and were inherited from their eukaryotic and prokaryotic ancestors and further diversified by the unprecedentedly high rates of gene duplication and functionalization experienced in land plants. Unlike prokaryotic microbes, which display frequent horizontal gene transfer events and multiple inputs of energy and organic carbon, land plants predominantly rely on organic carbon generated from CO2 and have experienced relatively few gene transfers during their recent evolutionary history. As such, plant metabolic networks have evolved in a stepwise manner using existing networks as a starting point and under various evolutionary constraints. That said, until recently, the evolution of only a handful of metabolic traits had been extensively investigated and as such, the evolution of metabolism has received a fraction of the attention of, the evolution of development, for example. Advances in metabolomics and next-generation sequencing have, however, recently led to a deeper understanding of how a wide range of plant primary and specialized (secondary) metabolic pathways have evolved both as a consequence of natural selection and of domestication and crop improvement processes. This article is part of the theme issue 'The evolution of plant metabolism'.

A Novel Function of GW5 on Controlling the Early Growth Vigor and its Haplotype Effect on Shoot Dry Weight and Grain Size in Rice (Oryza sativa L.)

Strong early growth vigor is an essential target in both direct seeded rice breeding and high-yielding rice breeding for rice varieties with relatively short growth duration in the double-cropping region. Shoot dry weight (SDW) is one of the important traits associated with early growth vigor, and breeders have been working to improve this trait. Finding stable QTLs or functional genes for SDW is crucial for improving the early growth vigor by implementing molecular breeding in rice. Here, a genome-wide association analysis revealed that the QTL for SDW, qSDW-5, was stably detected in the three cultivation methods commonly used in production practice. Through gene-based haplotype analysis of the annotated genes within the putative region of qSDW-5, and validated by gene expression and knockout transgenic experiments, LOC_Os05g09520, which is identical to the reported GW5/GSE5 controlling grain width (GW) and thousand grain weight (TGW) was identified as the causal gene for qSDW-5. Five main haplotypes of LOC_Os05g09520 were identified in the diverse international rice collection used in this study and their effects on SDW, GW and TGW were analyzed. Phenotypic comparisons of the major haplotypes of LOC_Os05g09520 in the three subpopulations (indica, japonica and aus) revealed the same patterns of wider GW and higher TGW along with higher SDW. Furtherly, the haplotype analysis of 138 rice varieties/lines widely used in southern China showed that 97.8% of the cultivars/lines carry Hap2LOC_Os05g09520. These results not only provide a promising gene source for the molecular breeding of rice varieties with strong early growth vigor, but also elucidate the effect of the LOC_Os05g09520 haplotypes on SDW, GW, and TGW in rice. Importantly, this study provides direct genetic evidence that these three traits are significantly correlated, and suggests a breeding strategy for developing high-yielding and slender grain-shaped indica cultivars with strong early growth vigor.

Arabidopsis Calmodulin-like Proteins CML13 and CML14 Interact with Calmodulin-Binding Transcriptional Activators and Function in Salinity Stress Response

Eukaryotic cells use calcium ions (Ca2+) as second messengers, particularly in response to abiotic and biotic stresses. These signals are detected by Ca2+ sensor proteins, such as calmodulin (CaM), which regulate the downstream target proteins. Plants also possess many CaM-like proteins (CMLs), most of which remain unstudied. We recently demonstrated that Arabidopsis CML13 and CML14 interact with proteins containing isoleucine/glutamine (IQ) domains, including CaM-binding transcriptional activators (CAMTAs). Here, we show that CaM, CML13 and CML14 bind all six members of the Arabidopsis CAMTA family. Using a combination of in planta and in vitro protein-interaction assays, we tested 11 members of the CaM/CML family and demonstrated that only CaM, CML13 and CML14 bind to CAMTA IQ domains. CaM, CML13 and CML14 showed Ca2+-independent binding to the IQ region of CAMTA6 and CAMTA3, and CAMTA6 in vitro exhibited some specificity toward individual IQ domains within CAMTA6 in split-luciferase in planta assays. We show that cml13 mutants exhibited enhanced salinity tolerance during germination compared to wild-type plants, a phenotype similar to camta6 mutants. In contrast, plants overexpressing CML13-GFP or CML14-GFP in the wild-type background showed increased NaCl sensitivity. Under mannitol stress, cml13 mutants were more susceptible than camta6 mutants or wild-type plants. The phenotype of cml13 mutants could be rescued with the wild-type CML13 gene. Several salinity-marker genes under CAMTA6 control were similarly misregulated in both camta6 and cml13 mutants, further supporting a role for CML13 in CAMTA6 function. Collectively, our data suggest that CML13 and CML14 participate in abiotic stress signaling as CAMTA effectors.

Genome-Wide Identification of the IQM Gene Family and Their Transcriptional Responses to Abiotic Stresses in Kiwifruit (Actinidia eriantha)

IQM is a plant-specific calcium-binding protein that plays a pivotal role in various aspects of plant growth response to stressors. We investigated the IQM gene family and its expression patterns under diverse abiotic stresses and conducted a comprehensive analysis and characterization of the AeIQMs, including protein structure, genomic location, phylogenetic relationships, gene expression profiles, salt tolerance, and expression patterns of this gene family under different abiotic stresses. Based on phylogenetic analysis, these 10 AeIQMs were classified into three distinct subfamilies (I-III). Analysis of the protein motifs revealed a considerable level of conservation among these AeIQM proteins within their respective subfamilies in kiwifruit. The genomic distribution of the 10 AeIQM genes spanned across eight chromosomes, where four pairs of IQM gene duplicates were associated with segmental duplication events. qRT-PCR analysis revealed diverse expression patterns of these AeIQM genes under different hormone treatments, and most AeIQMs showed inducibility by salt stress. Further investigations indicated that overexpression of AeIQMs in yeast significantly enhanced salt tolerance. These findings suggest that AeIQM genes might be involved in hormonal signal transduction and response to abiotic stress in Actinidia eriantha. In summary, this study provides valuable insights into the physiological functions of IQMs in kiwifruit.

The role of intrinsic disorder in binding of plant microtubule-associated proteins to the cytoskeleton

Microtubules (MTs) represent one of the main components of the eukaryotic cytoskeleton and support numerous critical cellular functions. MTs are in principle tube-like structures that can grow and shrink in a highly dynamic manner; a process largely controlled by microtubule-associated proteins (MAPs). Plant MAPs are a phylogenetically diverse group of proteins that nonetheless share many common biophysical characteristics and often contain large stretches of intrinsic protein disorder. These intrinsically disordered regions are determinants of many MAP-MT interactions, in which structural flexibility enables low-affinity protein-protein interactions that enable a fine-tuned regulation of MT cytoskeleton dynamics. Notably, intrinsic disorder is one of the major obstacles in functional and structural studies of MAPs and represents the principal present-day challenge to decipher how MAPs interact with MTs. Here, we review plant MAPs from an intrinsic protein disorder perspective, by providing a complete and up-to-date summary of all currently known members, and address the current and future challenges in functional and structural characterization of MAPs.

IQ domain-containing protein ZmIQD27 modulates water transport in maize

Plant metaxylem vessels provide physical support to promote upright growth and the transport of water and nutrients. A detailed characterization of the molecular network controlling metaxylem development is lacking. However, knowledge of the events that regulate metaxylem development could contribute to the development of germplasm with improved yield. In this paper, we screened an EMS-induced B73 mutant library, which covers 92% of maize (Zea mays) genes, to identify drought-sensitive phenotypes. Three mutants were identified, named iqd27-1, iqd27-2, and iqd27-3, and genetic crosses showed that they were allelic to each other. The causal gene in these 3 mutants encodes the IQ domain-containing protein ZmIQD27. Our study showed that defective metaxylem vessel development likely causes the drought sensitivity and abnormal water transport phenotypes in the iqd27 mutants. ZmIQD27 was expressed in the root meristematic zone where secondary cell wall deposition is initiated, and loss-of-function iqd27 mutants exhibited a microtubular arrangement disorder. We propose that association of functional ZmIQD27 with microtubules is essential for correct targeted deposition of the building blocks for secondary cell wall development in maize.

Molecular and genetic regulations of fleshy fruit shape and lessons from Arabidopsis and rice

Fleshy fruit shape is an important external quality trait influencing the usage of fruits and consumer preference. Thus, modification of fruit shape has become one of the major objectives for crop improvement. However, the underlying mechanisms of fruit shape regulation are poorly understood. In this review we summarize recent progress in the genetic basis of fleshy fruit shape regulation using tomato, cucumber, and peach as examples. Comparative analyses suggest that the OFP-TRM (OVATE Family Protein - TONNEAU1 Recruiting Motif) and IQD (IQ67 domain) pathways are probably conserved in regulating fruit shape by primarily modulating cell division patterns across fleshy fruit species. Interestingly, cucumber homologs of FRUITFULL (FUL1), CRABS CLAW (CRC) and 1-aminocyclopropane-1-carboxylate synthase 2 (ACS2) were found to regulate fruit elongation. We also outline the recent progress in fruit shape regulation mediated by OFP-TRM and IQD pathways in Arabidopsis and rice, and propose that the OFP-TRM pathway and IQD pathway coordinate regulate fruit shape through integration of phytohormones, including brassinosteroids, gibberellic acids, and auxin, and microtubule organization. In addition, functional redundancy and divergence of the members of each of the OFP, TRM, and IQD families are also shown. This review provides a general overview of current knowledge in fruit shape regulation and discusses the possible mechanisms that need to be addressed in future studies.

The cotton protein GhIQD21 interacts with GhCaM7 and modulates organ morphogenesis in Arabidopsis by influencing microtubule stability

KEY MESSAGE

GhIQD21, a cotton IQ67-domain protein, interacts with GhCaM7 and alters organ shape in Arabidopsis by modulating microtubule stability. Calcium ion (Ca²⁺) and the calcium sensor calmodulin play crucial roles in the growth and development of plants. GhCaM7, a calmodulin in upland cotton (Gossypium hirsutum L.), is highly expressed in cotton fiber cells during the rapid elongation period and plays an important role in fiber cell development. In this study, we screened for GhCaM7-interacting proteins and identified GhIQD21, which contains a typical IQ67-domain. GhIQD21 was preferentially expressed at the fiber rapid elongation stage, and the protein localized to microtubules (MTs). Ectopic expression of GhIQD21 in Arabidopsis resulted in shorter leaves, petals, siliques, and plant height, thicker inflorescences, and more trichomes when compared with wild type (WT). Further investigation indicated that the morphogenesis of leaf epidermal cells and silique cells was altered. There was less consistency in the orientation of cortical microtubules in cotyledon and hypocotyl epidermal cells. Furthermore, compared with WT, transgenic seedling hypocotyls were more sensitive to oryzalin, a MT depolymerization drug. These results indicated that GhIQD21 is a GhCaM7-interacting protein located in MTs and that it plays a role in plant growth and potentially cotton fiber development. This study provides a foundation for further studies of the function and regulatory mechanism of GhIQD21 in fiber cell development.

Genome-wide identification and analysis of the IQM gene family in soybean

IQM, a plant-specific calmodulin-binding protein, plays multiple roles in plant growth and development. Although a comprehensive analysis has been carried out on the IQM family genes in Arabidopsis and rice, the number and functions of IQM genes in other species have not been explored. In this study, we identified 15 members of the soybean (Glycine max) IQM gene family using BLASTP tools. These members were distributed on 12 soybean chromosomes and constitute six pairs caused by fragment duplication events. According to phylogeny, the 15 genes were divided into three subfamilies (I, II, and III), and members of the same subfamily had similar gene and protein structures. Yeast two-hybrid experiments revealed that the IQ motif is critical for the binding of GmIQM proteins to GmCaM, and its function is conserved in soybean, Arabidopsis, and rice. Based on real-time PCR, the soybean IQM genes were strongly induced by PEG and NaCl, suggesting their important biological functions in abiotic stress responses. Overall, this genome-wide analysis of the soybean IQM gene family lays a solid theoretical foundation for further research on the functions of GmIQM genes and could serve as a reference for the improvement and breeding of soybean stress resistance traits.

Genome-Wide Identification of the SUN Gene Family in Melon (Cucumis melo) and Functional Characterization of Two CmSUN Genes in Regulating Fruit Shape Variation

Melon (Cucumis melo) is an important economic crop cultivated worldwide. A unique SUN gene family plays a crucial role in regulating plant growth and fruit development, but many SUN family genes and their function have not been well-characterized in melon. In the present study, we performed genome-wide identification and bioinformatics analysis and identified 24 CmSUN family genes that contain integrated and conserved IQ67 domain in the melon genome. Transcriptome data analysis and qRT-PCR results showed that most CmSUNs are specifically enriched in melon reproductive organs, such as young flowers and ovaries. Through genetic transformation in melons, we found that overexpression of CmSUN23-24 and CmSUN25-26-27c led to an increased fruit shape index, suggesting that they act as essential regulators in melon fruit shape variation. Subcellular localization revealed that the CmSUN23-24 protein is located in the cytoplasmic membrane. A direct interaction between CmSUN23-24 and a Calmodulin protein CmCaM5 was found by yeast two-hybrid assay, which indicated their participation in the calcium signal transduction pathway in regulating plant growth. These findings revealed the molecular characteristics, expression profile, and functional pattern of the CmSUN genes, and may provide the theoretical basis for the genetic improvement of melon fruit breeding.

The wheat TaIQD3D-6 gene encodes a microtubule-associated protein and regulates cell morphogenesis in Arabidopsis

A plethora of microtubule-associated proteins (MAPs) modulate the dynamics of microtubules (MTs) to ensure the proper elaboration of developmental programs in plants. Among the plant-specific MAPs are the IQ67 domain (IQD) family proteins. Despite the great progress in elucidating IQD protein functions, the majority of IQD proteins, especially IQDs in crop species, remain to be functionally explored. In this study, we identified 78 putative IQD family genes in the genome of hexaploid wheat (Triticum aestivum). Phylogenetic analysis of wheat and Arabidopsis IQDs supports the previous notion that the expansion of the IQD family coincides with plant terrestrialization. Further characterization of one TaIQD, TaIQD3D-6, revealed that TaIQD3D-6 directly binds to MTs and free tubulins in vitro and is associated with cortical MTs in interphase cells in vivo. Overexpressing TaIQD3D-6 in Arabidopsis leads to a spectrum of phenotypes that are indicative of perturbed MT homeostasis, including spiral growth, hypersensitivity to MT-destabilizing drugs, defects in cell morphogenesis, and altered organization of cMT arrays. Finally, we determined that TaIQD3D-6-GFP localizes to the expanding cell plate during cytokinesis and the overexpression of TaIQD3D-6 interferes with asymmetric cell division in the stomatal lineage in Arabidopsis. In summary, our findings establish that TaIQD3D-6 is a MAP that regulates plant cell and organ morphogenesis and provide new insights into the functions of crop IQD proteins.

Genome-wide analysis of IQD proteins and ectopic expression of watermelon ClIQD24 in tomato suggests its important role in regulating fruit shape

The plant-specific IQ67 domain (IQD) is the largest class of calmodulin targets found in plants, and plays an important role in many biological processes, especially fruit development processes. However, the functional role of IQD proteins in the development of watermelon (Citrullus lanatus) shape remains unknown, as the IQD protein family in watermelon has not been systematically characterized. Herein, we elucidated the gene structures, chromosomal locations, evolutionary divergence, and functions of 35 IQD genes in the watermelon genome. The transcript profiles and quantitative real-time PCR analysis at different stages of fruit development showed that the ClIQD24 gene was highly expressed on 0 days after pollination. Furthermore, we found that the ectopic overexpression of ClIQD24 promoted tomato fruit elongation, thereby revealing the significance of ClIQD24 in the progression of watermelon shape. Our study will serve as a reference for further investigations on the molecular mechanisms underlying watermelon fruit shape formation.

VvSUN may act in the auxin pathway to regulate fruit shape in grape

Fruit shape is an essential agronomic feature in many crops. We identified and functionally characterized an auxin pathway-related gene, VvSUN. VvSUN, which belongs to the SUN/IQ67-DOMAIN (IQD) family, localizes to the plasma membrane and chloroplast and may be involved in controlling fruit shape through auxin. It is highly expressed in the ovary, and the expression level 1 week before the anthesis stage is positively correlated with the fruit shape index. Functional analyses illustrated that VvSUN gene overexpression in tomato and tobacco plants changed fruit/pod shape. The VvSUN promoter directly bound to VvARF6 in yeast and activated ß-glucuronidase (GUS) activity by indole-3-acetic acid (IAA) treatments in grapevine leaves, indicating that VvSUN functions are in coordination with auxin. Further analysis of 35S::VvSUN transgenic tomato ovaries showed that the fruit shape changes caused by VvSUN were predominantly caused by variations in cell number in longitudinal directions by regulating endogenous auxin levels via polar transport and/or auxin signal transduction process variations. Moreover, enrichment of the 35S::VvSUN transgenic tomato differentially expressed genes was found in a variety of biological processes, including primary metabolic process, transmembrane transport, calcium ion binding, cytoskeletal protein binding, tubulin binding, and microtubule-based movement. Using weighted gene co-expression network analysis (WGCNA), we confirmed that this plant hormone signal transduction may play a crucial role in controlling fruit shape. As a consequence, it is possible that VvSUN acts as a hub gene, altering cellular auxin levels and the plant hormone signal transduction pathway, which plays a role in cell division patterns, leading to anisotropic growth of the ovary and, ultimately, an elongated fruit shape.

AAQSP increases mapping resolution of stable QTLs through applying NGS-BSA in multiple genetic backgrounds

In this study, we present AAQSP as an extension of existing NGS-BSA applications for identifying stable QTLs at high resolution. GhPAP16 and GhIQD14 fine mapped on chromosome D09 of upland cotton are identified as important candidate genes for lint percentage (LP). Bulked segregant analysis combined with next generation sequencing (NGS-BSA) allows rapid identification of genome sequence differences responsible for phenotypic variation. The NGS-BSA approach applied to crops mainly depends on comparing two bulked DNA samples of individuals from an F₂ population. Since some F₂ individuals still maintain high heterozygosity, heterosis will exert complications in pursuing NGS-BSA in such populations. In addition, the genetic background influences the stability of gene expression in crops, so some QTLs mapped in one segregating population may not be widely applied in crop improvement. The AAQSP (Association Analysis of QTL-seq on Semi-homologous Populations) reported in our study combines the optimized scheme of constructing BSA bulks with NGS-BSA analysis in two (or more) different parental genetic backgrounds for isolating the stable QTLs. With application of AAQSP strategy and construction of a high-density linkage map, we have successfully identified a QTL significantly related to lint percentage (LP) in cultivated upland cotton, followed by map-based cloning to dissect two candidate genes, GhPAP16 and GhIQD14. This study demonstrated that AAQSP can efficiently identify stable QTLs for complex traits of interest, and thus accelerate the genetic improvement of upland cotton and other crop plants.

Genomic insights into genetic diploidization in the homosporous fern Adiantum nelumboides

Whole genome duplication has been recognized as a major process in speciation of land plants, especially in ferns. Whereas genome downsizing contributes greatly to the post-genome shock responses of polyploid flowering plants, diploidization of polyploid ferns diverges by maintaining most of the duplicated DNA and is thus expected to be dominated by genic processes. As a consequence, fern genomes provide excellent opportunities to study ecological speciation enforced by expansion of protein families via polyploidy. To test the key predictions of this hypothesis, we reported the de novo genome sequence of Adiantum nelumboides, a tetraploid homosporous fern. The obtained draft genome had a size of 6.27 Gb assembled into 11,767 scaffolds with the contig N50 of 1.37 Mb. Repetitive DNA sequences contributed with about 81.7%, a remarkably high proportion of the genome. With 69,568 the number of predicted protein-coding genes exceeded those reported in most other land plant genomes. Intragenomic synteny analyses recovered 443 blocks with the average block size of 1.29 Mb and the average gene content of 16 genes. The results are consistent with the hypothesis of high ancestral chromosome number, lack of substantial genome downsizing, and dominance of genic diploidization. As expected in the calciphilous plants, a notable number of detected genes were involved in calcium uptake and transport. In summary, the genome sequence of a tetraploid homosporous fern not only provides access to a genomic resource of a derived fern, but also supports the hypothesis of maintenance of high chromosome numbers and duplicated DNA in young polyploid ferns.

Microtubule-associated IQD9 orchestrates cellulose patterning in seed mucilage

Arabidopsis seeds release large capsules of mucilaginous polysaccharides, which are shaped by an intricate network of cellulosic microfibrils. Cellulose synthase complexes are guided by the microtubule cytoskeleton, but it is unclear which proteins mediate this process in the seed coat epidermis. Using reverse genetics, we identified IQ67 DOMAIN 9 (IQD9) and KINESIN LIGHT CHAIN-RELATED 1 (KLCR1) as two highly expressed genes during seed development and comprehensively characterized their roles in cell wall polysaccharide biosynthesis. Mutations in IQD9 as well as in KLCR1 lead to compact mucilage capsules with aberrant cellulose distribution, which can be rescued by transgene complementation. IQD9 physically interacts with KLCR1 and localizes to cortical microtubules (MTs) to maintain their organization in seed coat epidermal (SCE) cells. IQD9 as well as a previously identified TONNEAU1 (TON1) RECRUITING MOTIF 4 (TRM4) protein act to maintain cellulose synthase velocity. Our results demonstrate that IQD9, KLCR1 and TRM4 are MT-associated proteins that are required for seed mucilage architecture. This study provides the first direct evidence that members of the IQD, KLCR and TRM families have overlapping roles in cell wall biosynthesis. Therefore, SCE cells provide an attractive system to further decipher the complex genetic regulation of polarized cellulose deposition.

Division site determination during asymmetric cell division in plants

During development, both animals and plants exploit asymmetric cell division (ACD) to increase tissue complexity, a process that usually generates cells dissimilar in size, morphology, and fate. Plants lack the key regulators that control ACD in animals. Instead, plants have evolved two unique cytoskeletal structures to tackle this problem: the preprophase band (PPB) and phragmoplast. The assembly of the PPB and phragmoplast and their contributions to division plane orientation have been extensively studied. However, how the division plane is positioned off the cell center during asymmetric division is poorly understood. Over the past 20 years, emerging evidence points to a critical role for polarly localized membrane proteins in this process. Although many of these proteins are species- or cell type specific, and the molecular mechanism underlying division asymmetry is not fully understood, common features such as morphological changes in cells, cytoskeletal dynamics, and nuclear positioning have been observed. In this review, we provide updates on polarity establishment and nuclear positioning during ACD in plants. Together with previous findings about symmetrically dividing cells and the emerging roles of developmental cues, we aim to offer evolutionary insight into a common framework for asymmetric division-site determination and highlight directions for future work.

IQD1 Involvement in Hormonal Signaling and General Defense Responses Against Botrytis cinerea

IQ Domain 1 (IQD1) is a novel Arabidopsis thaliana calmodulin-binding protein, which was found to be a positive regulator of glucosinolate (GS) accumulation and plant defense responses against insects. We demonstrate here that the IQD1 overexpressing line (IQD1 ^OXP ) was also more resistant also to the necrotrophic fungus Botrytis cinerea, whereas an IQD1 knockout line (iqd1-1) was much more sensitive. Furthermore, we showed that IQD1 is up-regulated by jasmonic acid (JA) and downregulated by salicylic acid (SA). A comparison of whole transcriptome expression between iqd1-1 and wild type plants revealed a substantial downregulation of genes involved in plant defense and hormone regulation. Further examination revealed a marked reduction of SA and increases in the levels of ethylene, JA and abscisic acid response genes in the iqd1-1 line. Moreover, quantification of SA, JA, and abscisic acids in IQD1 ^OXP and iqd1-1 lines relative to the wild type, showed a significant reduction in endogenous JA levels in the knockout line, simultaneously with increased SA levels. Relations between IQD1 ^OXP and mutants defective in plant-hormone response indicated that IQD1 cannot rescue the absence of NPR1 or impaired SA accumulation in the NahG line. IQD1 cannot rescue ein2 or eto1 mutations connected to the ethylene pathway involved in both defense responses against B. cinerea and in regulating GS accumulation. Furthermore, IQD1cannot rescue the aos, coi1 or jar1mutations, all involved in the defense response against B. cinerea and it depends on JAR1 to control indole glucosinolate accumulation. We also found that in the B. cinerea, which infected the iqd1-1 mutant, the most abundant upregulated group of proteins is involved in the degradation of complex carbohydrates, as correlated with the sensitivity of this mutant. In summary, our results suggest that IQD1 is an important A. thaliana defensive protein against B. cinerea that is integrated into several important pathways, such as those involved in plant defense and hormone responses.

Identification of the Wheat (Triticum aestivum) IQD Gene Family and an Expression Analysis of Candidate Genes Associated with Seed Dormancy and Germination

The IQ67 Domain (IQD) gene family plays important roles in plant developmental processes and stress responses. Although IQDs have been characterized in model plants, little is known about their functions in wheat (Triticum aestivum), especially their roles in the regulation of seed dormancy and germination. Here, we identified 73 members of the IQD gene family from the wheat genome and phylogenetically separated them into six major groups. Gene structure and conserved domain analyses suggested that most members of each group had similar structures. A chromosome positional analysis showed that TaIQDs were unevenly located on 18 wheat chromosomes. A synteny analysis indicated that segmental duplications played significant roles in TaIQD expansion, and that the IQD gene family underwent strong purifying selection during evolution. Furthermore, a large number of hormone, light, and abiotic stress response elements were discovered in the promoters of TaIQDs, implying their functional diversity. Microarray data for 50 TaIQDs showed different expression levels in 13 wheat tissues. Transcriptome data and a quantitative real-time PCR analysis of wheat varieties with contrasting seed dormancy and germination phenotypes further revealed that seven genes (TaIQD4/-28/-32/-58/-64/-69/-71) likely participated in seed dormancy and germination through the abscisic acid-signaling pathway. The study results provide valuable information for cloning and a functional investigation of candidate genes controlling wheat seed dormancy and germination; consequently, they increase our understanding of the complex regulatory networks affecting these two traits.

Genome-wide identification, expression analysis and evolutionary relationships of the IQ67-domain gene family in common wheat (Triticum aestivum L.) and its progenitors

BACKGROUND

The plant-specific IQ67-domain (IQD) gene family plays an important role in plant development and stress responses. However, little is known about the IQD family in common wheat (Triticum aestivum L), an agriculturally important crop that provides more than 20% of the calories and protein consumed in the modern human diet.

RESULTS

We identified 125 IQDs in the wheat genome and divided them into four subgroups by phylogenetic analysis. The IQDs belonging to the same subgroup had similar exon-intron structure and conserved motif composition. Polyploidization contributed significantly to the expansion of IQD genes in wheat. Characterization of the expression profile of these genes revealed that a few T. aestivum (Ta)IQDs showed high tissue-specificity. The stress-induced expression pattern also revealed a potential role of TaIQDs in environmental adaptation, as TaIQD-2A-2, TaIQD-3A-9 and TaIQD-1A-7 were significantly induced by cold, drought and heat stresses, and could be candidates for future functional characterization. In addition, IQD genes in the A, B and D subgenomes displayed an asymmetric evolutionary pattern, as evidenced by their different gain or loss of member genes, expression levels and nucleotide diversity.

CONCLUSIONS

This study elucidated the potential biological functions and evolutionary relationships of the IQD gene family in wheat and revealed the divergent fates of IQD genes during polyploidization.

Genetic and Molecular Regulation Mechanisms in the Formation and Development of Vegetable Fruit Shape

Vegetable crops have a long history of cultivation worldwide and rich germplasm resources. With its continuous development and progress, molecular biology technology has been applied to various fields of vegetable crop research. Fruit is an important organ in vegetable crops, and fruit shape can affect the yield and commercialization of vegetables. In nature, fruits show differences in size and shape. Based on fruit shape diversity, the growth direction and coordination mechanism of fruits remain unclear. In this review, we discuss the latest research on fruit shape. In addition, we compare the current theories on the molecular mechanisms that regulate fruit growth, size, and shape in different vegetable families.

Multiparameter in vivo imaging in plants using genetically encoded fluorescent indicator multiplexing

Biological processes are highly dynamic, and during plant growth, development, and environmental interactions, they occur and influence each other on diverse spatiotemporal scales. Understanding plant physiology on an organismic scale requires analyzing biological processes from various perspectives, down to the cellular and molecular levels. Ideally, such analyses should be conducted on intact and living plant tissues. Fluorescent protein (FP)-based in vivo biosensing using genetically encoded fluorescent indicators (GEFIs) is a state-of-the-art methodology for directly monitoring cellular ion, redox, sugar, hormone, ATP and phosphatidic acid dynamics, and protein kinase activities in plants. The steadily growing number of diverse but technically compatible genetically encoded biosensors, the development of dual-reporting indicators, and recent achievements in plate-reader-based analyses now allow for GEFI multiplexing: the simultaneous recording of multiple GEFIs in a single experiment. This in turn enables in vivo multiparameter analyses: the simultaneous recording of various biological processes in living organisms. Here, we provide an update on currently established direct FP-based biosensors in plants, discuss their functional principles, and highlight important biological findings accomplished by employing various approaches of GEFI-based multiplexing. We also discuss challenges and provide advice for FP-based biosensor analyses in plants.

Complex cellular and molecular events determining fruit size

The understanding of plant organ-size determination represents an important challenge, especially because of the significant role of plants as food and renewable energy sources and the increasing need for plant-derived products. Most of the knowledge on the regulation of organ growth and the molecular network controlling cell division and cell expansion, the main drivers of growth, is derived from arabidopsis. The increasing use of crops such as tomato for research is now bringing essential information on the mechanisms underlying size control in agronomically important organs. This review describes our current knowledge, still very scarce, of the cellular and molecular mechanisms governing tomato fruit size and proposes future research to better understand the regulation of growth in this important crop.

Genome-Wide Analysis of the IQM Gene Family in Rice (Oryza sativa L.)

Members of the IQM (IQ-Motif Containing) gene family are involved in plant growth and developmental processes, biotic and abiotic stress response. To systematically analyze the IQM gene family and their expression profiles under diverse biotic and abiotic stresses, we identified 8 IQM genes in the rice genome. In the current study, the whole genome identification and characterization of OsIQMs, including the gene and protein structure, genome localization, phylogenetic relationship, gene expression and yeast two-hybrid were performed. Eight IQM genes were classified into three subfamilies (I-III) according to the phylogenetic analysis. Gene structure and protein motif analyses showed that these IQM genes are relatively conserved within each subfamily of rice. The 8 OsIQM genes are distributed on seven out of the twelve chromosomes, with three IQM gene pairs involved in segmental duplication events. The evolutionary patterns analysis revealed that the IQM genes underwent a large-scale event within the last 20 to 9 million years. In addition, quantitative real-time PCR analysis of eight OsIQMs genes displayed different expression patterns at different developmental stages and in different tissues as well as showed that most IQM genes were responsive to PEG, NaCl, jasmonic acid (JA), abscisic acid (ABA) treatment, suggesting their crucial roles in biotic, and abiotic stress response. Additionally, a yeast two-hybrid assay showed that OsIQMs can interact with OsCaMs, and the IQ motif of OsIQMs is required for OsIQMs to combine with OsCaMs. Our results will be valuable to further characterize the important biological functions of rice IQM genes.

Genome-Wide Identification and Characterization of the Potato IQD Family During Development and Stress

Calmodulin-binding proteins belong to the IQ67 domain (IQD) gene family and play essential roles in plant development and stress responses. However, the role of IQD gene family in potato (Solanum tuberosum L.) is yet to be known. In the present study, 23 StIQDs were identified in the potato genome and named StIQD1 to StIQD23. They were unevenly distributed on 10 of the 12 chromosomes. Phylogenetic analysis divided the IQDs into four subfamilies (IQD I–IV). StIQDs found in three of the four subfamilies. Synteny analysis confirmed that potato and tomato shared a close evolutionary relationship. Besides, RNA-Seq data analysis revealed that the expression of 19 of the 23 StIQDs was detected in at least one of the 12 tissues, and some of which showed a tissue-specific pattern. Quantitative reverse transcriptase–polymerase chain reaction results further confirmed that 14 StIQDs responded differently to various abiotic stresses, including drought, extreme temperature, and CaCl₂ treatment, suggesting their significance in stress response. This study presents a comprehensive overview of the potato IQD gene family and lays a foundation for further analysis of the StIQDs functions in plant development and stress response.

Genome wide analysis of IQD gene family in diploid and tetraploid species of cotton (Gossypium spp.)

Calmodulin (CaM) is considered as the most significant Ca²⁺ signaling messenger that mediate various biochemical and physiological reactions. IQ domain (IQD) proteins are plant specific CML/CaM calcium binding which are characterized by domains of 67 amino acids. 50, 50, 94, and 99 IQD genes were detected from G. arboreum (A2), G. raimondii (D5), G. barbadense (AD2) and G. hirsutum (AD1) respectively. Existence of more orthologous genes in cotton species than Arabidopsis, advocated that polyploidization produced new cotton specific orthologous gene clusters. Duplication of gene events depicts that IQD gene family of cotton evolution was under strong purifying selection. G. hirsutum exhibited high level synteny. GarIQD25 exhibited high expression in stem, root, flower, ovule and fiber in G. arboreum. In G. raimondii, GraIQD03 demonstrated upregulation across stem, ovule, fiber and seed. GbaIQD11 and GbaIQD62 exhibited upregulation in fiber development in G. barbadense. GhiIQD69 recognized as main candidate genes for plant parts, floral tissues, fiber and ovule development. Promotor analysis identified cis-regulatory elements were involved in plant growth and development. Overwhelmingly, present study paves the way to better understand the evolution of cotton IQD genes and lays a foundation for future investigation of IQD in cotton.

Ectopic expression of GhIQD14 (cotton IQ67 domain-containing protein 14) causes twisted organ and modulates secondary wall formation in Arabidopsis

In plants, although KNOX genes are known to regulate secondary cell wall (SCW) formation, their protein-regulating mechanisms remain largely unknown. Here, we showed that GhKNL1, which regulates SCW formation and fiber development in cotton, could interact with an IQ67 domain containing protein (GhIQD14) in yeast. Confocal observation showed that GhIQD14 was localized to the microtubules. In Arabidopsis, ectopic expression of GhIQD14 caused hypocotyls to be sensitive to microtubule depolymerization agent, organ twisting of seedlings, trichomes, rosette leaves, and capsules, as well as severely irregular xylem vessels and thicker interfascicular fiber cell walls in the inflorescence stem. Furthermore, we found that GhIQD14 interacted with AtKNAT7 in vivo, and instantaneous co-expression of GhIQD14 and AtKNAT7 in tobacco showed that GhIQD14 weakened the distribution of AtKNAT7 in the nucleus, bringing it into the microtubules, thus affecting the SCW formation related genes expression. Our results suggested that GhIQD14 might be involved in the morphological development and SCW formation in cotton.

IQ67 DOMAIN proteins facilitate preprophase band formation and division-plane orientation

Spatiotemporal control of cell division is essential for the growth and development of multicellular organisms. In plant cells, proper cell plate insertion during cytokinesis relies on the premitotic establishment of the division plane at the cell cortex. Two plant-specific cytoskeleton arrays, the preprophase band (PPB) and the phragmoplast, play important roles in division-plane orientation and cell plate formation, respectively1. Microtubule organization and dynamics and their communication with membranes at the cortex and cell plate are coordinated by multiple, mostly distinct microtubule-associated proteins2. How division-plane selection and establishment are linked, however, is still unknown. Here, we report members of the Arabidopsis IQ67 DOMAIN (IQD) family3 as microtubule-targeted proteins that localize to the PPB and phragmoplast and additionally reside at the cell plate and a polarized cortical region including the cortical division zone (CDZ). IQDs physically interact with PHRAGMOPLAST ORIENTING KINESIN (POK) proteins4,5 and PLECKSTRIN HOMOLOGY GTPase ACTIVATING (PHGAP) proteins6, which are core components of the CDZ1. The loss of IQD function impairs PPB formation and affects CDZ recruitment of POKs and PHGAPs, resulting in division-plane positioning defects. We propose that IQDs act as cellular scaffolds that facilitate PPB formation and CDZ set-up during symmetric cell division.

The domain of unknown function 4005 (DUF4005) in an Arabidopsis IQD protein functions in microtubule binding

The dynamic responses of microtubules (MTs) to internal and external signals are modulated by a plethora of microtubule-associated proteins (MAPs). In higher plants, many plant-specific MAPs have emerged during evolution as advantageous to their sessile lifestyle. Some members of the IQ67 domain (IQD) protein family have been shown to be plant-specific MAPs. However, the mechanisms of interaction between IQD proteins and MTs remain elusive. Here we demonstrate that the domain of unknown function 4005 (DUF4005) of the Arabidopsis IQD family protein ABS6/AtIQD16 is a novel MT-binding domain. Cosedimentation assays showed that the DUF4005 domain binds directly to MTs in vitro. GFP-labeled DUF4005 also decorates all types of MT arrays tested in vivo. Furthermore, we showed that a conserved stretch of 15 amino acid residues within the DUF4005 domain, which shares sequence similarity with the C-terminal MT-binding domain of human MAP Kif18A, is required for the binding to MTs. Transgenic lines overexpressing the DUF4005 domain displayed a spectrum of developmental defects, including spiral growth and stunted growth at the organismal level. At the cellular level, DUF4005 overexpression caused defects in epidermal pavement cell and trichome morphogenesis, as well as abnormal anisotropic cell elongation in the hypocotyls of dark-grown seedlings. These data establish that the DUF4005 domain of ABS6/AtIQD16 is a new MT-binding domain, overexpression of which perturbs MT homeostasis in plants. Our findings provide new insights into the MT-binding mechanisms of plant IQD proteins.

ABNORMAL SHOOT 6 interacts with KATANIN 1 and SHADE AVOIDANCE 4 to promote cortical microtubule severing and ordering in Arabidopsis

Plant interphase cortical microtubules (cMTs) mediate anisotropic cell expansion in response to environmental and developmental cues. In Arabidopsis thaliana, KATANIN 1 (KTN1), the p60 catalytic subunit of the conserved MT-severing enzyme katanin, is essential for cMT ordering and anisotropic cell expansion. However, the regulation of KTN1-mediated cMT severing and ordering remains unclear. In this work, we report that the Arabidopsis IQ67 DOMAIN (IQD) family gene ABNORMAL SHOOT 6 (ABS6) encodes a MT-associated protein. Overexpression of ABS6 leads to elongated cotyledons, directional pavement cell expansion, and highly ordered transverse cMT arrays. Genetic suppressor analysis revealed that ABS6-mediated cMT ordering is dependent on KTN1 and SHADE AVOIDANCE 4 (SAV4). Live imaging of cMT dynamics showed that both ABS6 and SAV4 function as positive regulators of cMT severing. Furthermore, ABS6 directly interacts with KTN1 and SAV4 and promotes their recruitment to the cMTs. Finally, analysis of loss-of-function mutant combinations showed that ABS6, SAV4, and KTN1 work together to ensure the robust ethylene response in the apical hook of dark-grown seedlings. Together, our findings establish ABS6 and SAV4 as positive regulators of cMT severing and ordering, and highlight the role of cMT dynamics in fine-tuning differential growth in plants.

Exploring genetic architecture for pod-related traits in soybean using image-based phenotyping

Mature pod color (PC) and pod size (PS) served as important characteristics are used in the soybean breeding programs. However, manual phenotyping of such complex traits is time-consuming, laborious, and expensive for breeders. Here, we collected pod images from two different populations, namely, a soybean association panel (SAP) consisting of 187 accessions and an inter-specific recombinant inbred line (RIL) population containing 284 RILs. An image-based phenotyping method was developed and used to extract the pod color- and size-related parameters from images. Genome-wide association study (GWAS) and linkage mapping were performed to decipher the genetic control of pod color- and size-related traits across 2 successive years. Both populations exhibited wide phenotypic variations and continuous distribution in pod color- and size-related traits, indicating quantitative polygenic inheritance of these traits. GWAS and linkage mapping approaches identified the two major quantitative trait loci (QTL) underlying the pod color parameters, i.e., qPC3 and qPC19, located to chromosomes 3 and 19, respectively, and 12 stable QTLs for pod size-related traits across nine chromosomes. Several genes residing within the genomic region of stable QTL were identified as potential candidates underlying these pod-related traits based on the gene annotation and expression profiling data. Our results provide the useful information for fine-mapping/map-based cloning of QTL and marker-assisted selection of elite varieties with desirable pod traits.

Supplementary Information

The online version contains supplementary material available at 10.1007/s11032-021-01223-2.

New Insights Into Functions of IQ67-Domain Proteins

IQ67-domain (IQD) proteins, first identified in Arabidopsis and rice, are plant-specific calmodulin-binding proteins containing highly conserved motifs. They play a critical role in plant defenses, organ development and shape, and drought tolerance. Driven by comprehensive genome identification and analysis efforts, IQDs have now been characterized in several species and have been shown to act as microtubule-associated proteins, participating in microtubule-related signaling pathways. However, the precise molecular mechanisms underpinning their biological functions remain incompletely understood. Here we review current knowledge on how IQD family members are thought to regulate plant growth and development by affecting microtubule dynamics or participating in microtubule-related signaling pathways in different plant species and propose some new insights.

Genome-Wide Identification and Characterization of Calcium Metabolism Related Gene Families in Arabidopsis thaliana and Their Regulation by Bacillus amyloliquefaciens Under High Calcium Stress

Several gene families involved in calcium signaling have been detected in plants, including calmodulin (CaM), calcium dependent protein kinases (CDPK), calcineurin B-like (CBL) and cyclic nucleotide-gated channels (CNGCs). In our previous study, we demonstrated that Bacillus amyloliquefaciens LZ04 (B. amyloliquefaciens LZ04) regulate genes involved in calcium stress in Arabidopsis thaliana (A. thaliana). Here, we aimed to explore the potential involvement of calcium-related gene families in the response of A. thaliana to calcium stress and the potential regulatory effects of B. amyloliquefaciens LZ04 on these genes. The structure, duplication, synteny, and expression profiles of 102 genes in calcium-related gene families in A. thaliana were investigated. Hidden Markov Models (HMMs) and BLASTP were used to predict candidate genes and conserved domains of the candidate genes were confirmed in SMART and NCBI CDD databases. Gene duplications and synteny were uncovered by BLASTP and phylogenetic analysis. The transcriptome expression profiles of candidate genes were investigated by strand-specific sequencing. Cluster analysis was used to find the expression profiles of calcium-related genes families under different treatment conditions. A total of 102 genes in calcium-related gene families were detected in A. thaliana genome, including 34 CDPK genes, 20 CNGC genes, 18 CIPK genes, 22 IQD genes, and 10 CBP genes. Additionally, of the 102 genes, 33 duplications (32.35%) and 26 gene pairs including 48 genes (47.06%) were detected. Treatment with B. amyloliquefaciens LZ04 enhanced the resistance of A. thaliana under high calcium stress by regulating some of the genes in the calcium-related gene families. Functional enrichment analysis revealed that the genes clustered in the 42nd expression profile which may be B. amyloliquefaciens-responsive genes under calcium stress were enriched in protein phosphorylation and protein modification process. Transcriptome data was validated by RT-PCR and the results generally corroborated the transcriptome sequencing results. These results may be useful for agricultural improvement in high calcium stress regions.

Characterization of microRNA genes from Pigeonpea (Cajanus cajan L.) and understanding their involvement in drought stress

MicroRNAs (miRNA) are non-coding 20-24 nucleotide long RNAs regulating gene expression. In this study, we have characterized and analysed expression of miRNAs in Pigeonpea by using bioinformatics and experimental tools. We identified 116 miRNAs belonging to 32 phylogenetic families. Further, transcription start sites of miRNA genes revealed abundance and unique arrangement of adenine at +1 and thymine at -1 position. Promoter analysis exhibited presence of 19 most prevalent motifs which comprises majorly of TATA box and MYC domains. In total, 252 miRNA-targets were identified and found to be involved in various developmental processes and stress responses. Moreover, genome-wide localization studies demonstrated clustering of cca-miRNA 395 and 169 genes. The tandem and segmental duplication events were observed suggesting miRNA genes have been originated parallelly with protein coding genes. The expression analysis revealed induction of cca-miR169a, 398a and 408 miRNAs under drought stress highlighting their involvement. Conversely, down-regulation of their putative targets (NFYA, SOD, and UCLA, respectively) confirmed regulatory role of miRNAs in their expression emphasising the negative relationship between these miRNAs and targets in Pigeonpea. This study reports vast repertoire of miRNA genes which further can be experimentally characterized to elucidate their functions in various biological processes and can be recommended for Pigeonpea improvement programs.

Comparative transcriptome analyses define genes and gene modules differing between two Populus genotypes with contrasting stem growth rates

BACKGROUND

Wood provides an important biomass resource for biofuel production around the world. The radial growth of tree stems is central to biomass production for forestry and biofuels, but it is challenging to dissect genetically because it is a complex trait influenced by many genes. In this study, we adopted methods of physiology, transcriptomics and genetics to investigate the regulatory mechanisms of tree radial growth and wood development.

RESULTS

Physiological comparison showed that two Populus genotypes presented different rates of radial growth of stems and accumulation of woody biomass. A comparative transcriptional network approach was used to define and characterize functional differences between two Populus genotypes. Analyses of transcript profiles from wood-forming tissue of the two genotypes showed that 1542, 2295 and 2110 genes were differentially expressed in the pre-growth, fast-growth and post-growth stages, respectively. The co-expression analyses identified modules of co-expressed genes that displayed distinct expression profiles. Modules were further characterized by correlating transcript levels with genotypes and physiological traits. The results showed enrichment of genes that participated in cell cycle and division, whose expression change was consistent with the variation of radial growth rates. Genes related to secondary vascular development were up-regulated in the faster-growing genotype in the pre-growth stage. We characterized a BEL1-like (BELL) transcription factor, PeuBELL15, which was up-regulated in the faster-growing genotype. Analyses of transgenic Populus overexpressing as well as CRISPR/Cas9-induced mutants for BELL15 showed that PeuBELL15 improved accumulation of glucan and lignin, and it promoted secondary vascular growth by regulating the expression of genes relevant for cellulose synthases and lignin biosynthesis.

CONCLUSIONS

This study illustrated that active division and expansion of vascular cambium cells and secondary cell wall deposition of xylem cells contribute to stem radial increment and biomass accumulation, and it identified relevant genes for these complex growth traits, including a BELL transcription factor gene PeuBELL15. This provides genetic resources for improving and breeding elite genotypes with fast growth and high wood biomass.

Rice microtubule-associated protein IQ67-DOMAIN14 regulates grain shape by modulating microtubule cytoskeleton dynamics

Cortical microtubule (MT) arrays play a critical role in plant cell shape determination by defining the direction of cell expansion. As plants continuously adapt to ever-changing environmental conditions, multiple environmental and developmental inputs need to be translated into changes of the MT cytoskeleton. Here, we identify and functionally characterize an auxin-inducible and MT-localized protein OsIQ67-DOMAIN14 (OsIQD14), which is highly expressed in rice seed hull cells. We show that while deficiency of OsIQD14 results in short and wide seeds and increases overall yield, overexpression leads to narrow and long seeds, caused by changed MT alignment. We further show that OsIQD14-mediated MT reordering is regulated by specifically affecting MT dynamics, and ectopic expression of OsIQD14 in Arabidopsis could change the cell shape both in pavement cells and in hypocotyl cells. Additionally, OsIQD14 activity is tightly controlled by calmodulin proteins, providing an alternative way to modify the OsIQD14 activity. Our results indicate that OsIQD14 acts as a key factor in regulating MT rearrangements in rice hull cells and hence the grain shape, and allows effective local cell shape manipulation to improve the rice yield trait.

Comprehensive Sequence Analysis of IQD Gene Family and their Expression Profiling in Grapevine (Vitis vinifera)

The plant-specific IQ67-domain (IQD) protein family members are downstream targets of calcium sensors, known to regulate plant growth and lateral organ polarity, and basal defense response against environmental cues. No systematic study of IQD gene family has been performed on grapevine. The public availability of grapevine genome enables us to perform identification, phylogeny, chromosomal orientation, and gene structure analysis of the IQD genes in grapevine. We identified 49 VvIQD genes (VvIQD1–VvIQD49) and further classified them into eight subgroups based on phylogenetic relationships. The 49 VvIQD genes were assigned to 19 different chromosomal positions. The collinear relationship between grapevine and Arabidopsis IQDs (VvIQD and AtIQD), and within grapevine VvIQDs, was highly conserved. In addition, most of duplicated gene pairs showed Ka/Ks ratio less than 1.00, indicating purifying selection within these gene pairs, implying functional discrepancy after duplication. Transcription profiling of VvIQD genes shed light on their specific role in grapevine tissue and organ development. The qRT-PCR validation of the 49 VvIQD genes in grape berry tissue from cultivars with distinct berry shape during developmental phases suggested candidate genes involved in the shape of grape berries. The subcellular prediction of VvIQD22, VvIQD23, VvIQD38, and VvIQD49 genes validated their localization in the nucleus and plasma membrane. The VvIQD49 protein interaction with VvCaM2 was also verified by bimolecular fluorescence complementation (BiFC) analysis in the plasma membrane. Our findings will be valuable for the functional genomic studies for desirable shape development of grape berries.

Genetic architecture of fruit size and shape variation in cucurbits: a comparative perspective

The Cucurbitaceae family hosts many economically important fruit vegetables (cucurbits) such as cucumber, melon, watermelon, pumpkin/squash, and various gourds. The cucurbits are probably best known for the diverse fruit sizes and shapes, but little is known about their genetic basis and molecular regulation. Here, we reviewed the literature on fruit size (FS), shape (FSI), and fruit weight (FW) QTL identified in cucumber, melon, and watermelon, from which 150 consensus QTL for these traits were inferred. Genome-wide survey of the three cucurbit genomes identified 253 homologs of eight classes of fruit or grain size/weight-related genes cloned in Arabidopsis, tomato, and rice that encode proteins containing the characteristic CNR (cell number regulator), CSR (cell size regulator), CYP78A (cytochrome P450), SUN, OVATE, TRM (TONNEAU1 Recruiting Motif), YABBY, and WOX domains. Alignment of the consensus QTL with candidate gene homologs revealed widespread structure and function conservation of fruit size/shape gene homologs in cucurbits, which was exemplified with the fruit size/shape candidate genes CsSUN25-26-27a and CsTRM5 in cucumber, CmOFP1a in melon, and ClSUN25-26-27a in watermelon. In cucurbits, the andromonoecy (for 1-aminocyclopropane-1-carboxylate synthase) and the carpel number (for CLAVATA3) loci are known to have pleiotropic effects on fruit shape, which may complicate identification of fruit size/shape candidate genes in these regions. The present work illustrates the power of comparative analysis in understanding the genetic architecture of fruit size/shape variation, which may facilitate QTL mapping and cloning for fruit size-related traits in cucurbits. The limitations and perspectives of this approach are also discussed.

Knockdown of GhIQD31 and GhIQD32 increases drought and salt stress sensitivity in Gossypium hirsutum

Drought, salinity and cold stresses have a major impact on cotton production, thus identification and utilization of plant genes vital for plant improvement Whole-genome identification and functional characterizations of the IQ67-domain (IQD) protein family was carried out in which 148, 77, and 79 IQD genes were identified in Gossypium hirsutum, G. raimondii, and G. arboreum. The entire IQD proteins had varied physiochemical properties, however; their grand hydropathy values were negative, which demonstrated that the proteins were hydrophilic, a property common among the proteins encoded by various stresses responsive genes, such as the late embryogenesis abundant (LEA) proteins. The IQD proteins were predicted to be majorly sublocalized in the nucleus; moreover, various cis-regulatory elements with higher role in enhancing abiotic stress tolerance were detected. RNA-seq and RT-qPCR analysis revealed two key genes, Gh_D06G0014 and Gh_A09G1608 with significantly higher upregulation across the various tissues under drought, salt and cold stress. Knockdown of the two genes negatively affected the ability of G. hirsutum to tolerate the effects of the three stress factors, being all the antioxidant assayed were significantly low concentrations compared to the oxidizing enzymes in VIGS plants under stress, furthermore, morphological and physiological traits were all negatively affected in VIGS plants. Expression levels of GhLEA2, GhCDK_F4, GPCR (TOM1) and Gh_A05G2067 (TH), the stress responsive genes were all downregulated in the VIGS plants, but significantly upregulated in WT and positively controlled plants. The results demonstrated that the IQD genes could be responsible for enhancing drought, salt and cold stress tolerance in cotton.

GW5-Like, a homolog of GW5, negatively regulates grain width, weight and salt resistance in rice

Grain size is an important determinant of yield potential in crops. We previously demonstrated that natural mutations in the regulatory sequences of qSW5/GW5 confer grain width diversity in rice. However, the biological function of a GW5 homolog, named GW5-Like (GW5L), remains unknown. In this study, we report on GW5L knockout mutants in Kitaake, a japonica cultivar (cv.) considered to have a weak gw5 variant allele that confers shorter and wider grains. GW5L is evenly expressed in various tissues, and its protein product is localized to the plasma membrane. Biochemical assays verified that GW5L functions in a similar fashion to GW5. It positively regulates brassinosteroid (BR) signaling through repression of the phosphorylation activity of GSK2. Genetic data show that GW5L overexpression in either Kitaake or a GW5 knockout line, Kasa^orf3 (indica cv. Kasalath background), causes more slender, longer grains relative to the wild-type. We also show that GW5L could confer salt stress resistance through an association with calmodulin protein OsCaM1-1. These findings identify GW5L as a negative regulator of both grain size and salt stress tolerance, and provide a potential target for breeders to improve grain yield and salt stress resistance in rice.

Enhancement of Glucosinolate Production in Watercress ( Nasturtium officinale) Hairy Roots by Overexpressing Cabbage Transcription Factors

Glucosinolates are secondary metabolites that play important roles in plant defense and human health, as their production in plants is enhanced by overexpressing transcription factors. Here, four cabbage transcription factors (IQD1-1, IQD1-2, MYB29-1, and MYB29-2) affecting genes in both aliphatic and indolic glucosinolates biosynthetic pathways and increasing glucosinolates accumulation were overexpressed in watercress. Five IQD1-1, six IQD1-2, five MYB29-1, six MYB29-2, and one GUS hairy root lines were created. The expression of all genes involved in glucosinolates biosynthesis was higher in transgenic lines than in the GUS hairy root line, in agreement with total glucosinolates contents, determined by high-performance liquid chromatography. In transgenic IQD1-1 (1), IQD1-2 (4), MYB29-1 (2), and MYB29-2 (1) hairy root lines, total glucosinolates were 3.39-, 3.04-, 2.58-, and 4.69-fold higher than those in the GUS hairy root lines, respectively. These results suggest a central regulatory function for IQD1-1, IQD1-2, MYB29-1, and MYB29-2 transcription factors in glucosinolates biosynthesis in watercress hairy roots.

Genome-wide analysis of the Chinese cabbage IQD gene family and the response of BrIQD5 in drought resistance

KEY MESSAGE

Thirty-five IQD genes were identified and analysed in Chinese cabbage and BrIQD5 transgenic plants enhanced the drought resistance of plants. The IQD (IQ67-domain) family plays an important role in various abiotic stress responses in plant species. However, the roles of IQD genes in the Chinese cabbage response to abiotic stress remain unclear. Here, 35 IQD genes, from BrIQD1 to BrIQD35, were identified in Chinese cabbage (Brassica rapa ssp. pekinensis). Based on the phylogenetic analysis, these genes were clustered into three subfamilies (I-III), and members within the same subfamilies shared conserved exon-intron distribution and motif composition. The 35 BrIQD genes were unevenly distributed on 9 of the 10 chromosomes with 4 segmental duplication events. Ka/Ks ratios showed that the duplicated BrIQDs had mainly experienced strong purifying selection. Quantitative real-time polymerase chain reaction of 35 BrIQDs under PEG6000 indicated that BrIQD5 was significantly induced by PEG6000. To verify BrIQD5 function, BrIQD5 was heterologously overexpressed in tobacco and was silenced in Chinese cabbage. BrIQD5-overexpressed plants showed more tolerance to drought stress than wild-type plants, while BrIQD5-silenced plants in Chinese cabbage showed decreased drought tolerance. Additionally, six BrIQD5 potential interactive proteins were isolated by the yeast two-hybrid assay, including BrCaMa, BrCaMb and four other stress-related proteins. Motif IQ1 of BrIQD5 is important for the interaction with BrCaMa and BrCaMb, and the isoleucine in motif IQ1 is an essential amino acid for calmodulin binding to BrIQD5. The identification and cloning of the new Chinese cabbage drought tolerance genes will promote the drought-resistant breeding of Chinese cabbage and help to better understand the mechanism of IQD involved in the drought tolerance of plants.

Genome-wide identification and functional analysis of the cyclic nucleotide-gated channel gene family in Chinese cabbage

Cyclic nucleotide-gated channels (CNGCs) are a class of nonselective cationic channels that are widely found in animals and plants. Plant CNGCs participate in numerous biological functions that vary from development to stress tolerance. Most CNGC genes have been identified in plant genomes, but no such comprehensive study has yet been conducted on Chinese cabbage. In this study, thirty BrCNGC genes were identified, divided into five groups, and used for evolutionary analysis. We assigned names of all individual CNGC members on the basis of phylogenetic relationship with A. thaliana CNGCs. All BrCNGC genes were randomly distributed on chromosomes, and the A08 chromosome did not carry any CNGC gene. The CNGC genes of Chinese cabbage and A. thaliana from the same group displayed similar conserved motifs and gene structures. Especially the closer the homology, the higher the similarity. Quantitative expression analysis showed that most of the CNGC genes were expressed under four stresses, indicating that they play a key role in the stress response of Chinese cabbage. Expression patterns of 12 BrCNGC in the roots, stems, leaves, flowers, and siliques showed that BrCNGC8 and BrCNGC16 were specifically expressed only in flowers but not in other parts. This study lays a theoretical foundation for future research on the function of the CNGC gene family in Chinese cabbage.

Calcium- and calmodulin-regulated microtubule-associated proteins as signal-integration hubs at the plasma membrane-cytoskeleton nexus

Plant growth and development are a genetically predetermined series of events but can change dramatically in response to environmental stimuli, involving perpetual pattern formation and reprogramming of development. The rate of growth is determined by cell division and subsequent cell expansion, which are restricted and controlled by the cell wall-plasma membrane-cytoskeleton continuum, and are coordinated by intricate networks that facilitate intra- and intercellular communication. An essential role in cellular signaling is played by calcium ions, which act as universal second messengers that transduce, integrate, and multiply incoming signals during numerous plant growth processes, in part by regulation of the microtubule cytoskeleton. In this review, we highlight recent advances in the understanding of calcium-mediated regulation of microtubule-associated proteins, their function at the microtubule cytoskeleton, and their potential role as hubs in crosstalk with other signaling pathways.

Microtubule-associated protein IQ67 DOMAIN5 regulates morphogenesis of leaf pavement cells in Arabidopsis thaliana

Plant microtubules form a highly dynamic intracellular network with important roles for regulating cell division, cell proliferation, and cell morphology. Their organization and dynamics are co-ordinated by various microtubule-associated proteins (MAPs) that integrate environmental and developmental stimuli to fine-tune and adjust cytoskeletal arrays. IQ67 DOMAIN (IQD) proteins recently emerged as a class of plant-specific MAPs with largely unknown functions. Here, using a reverse genetics approach, we characterize Arabidopsis IQD5 in terms of its expression domains, subcellular localization, and biological functions. We show that IQD5 is expressed mostly in vegetative tissues, where it localizes to cortical microtubule arrays. Our phenotypic analysis of iqd5 loss-of-function lines reveals functions of IQD5 in pavement cell (PC) shape morphogenesis. Histochemical analysis of cell wall composition further suggests reduced rates of cellulose deposition in anticlinal cell walls, which correlate with reduced anisotropic expansion. Lastly, we demonstrate IQD5-dependent recruitment of calmodulin calcium sensors to cortical microtubule arrays and provide first evidence for important roles for calcium in regulation of PC morphogenesis. Our work identifies IQD5 as a novel player in PC shape regulation and, for the first time, links calcium signaling to developmental processes that regulate anisotropic growth in PCs.

Microtubule-associated protein IQ67 DOMAIN5 regulates morphogenesis of leaf pavement cells in Arabidopsis thaliana

Plant microtubules form a highly dynamic intracellular network with important roles for regulating cell division, cell proliferation, and cell morphology. Their organization and dynamics are co-ordinated by various microtubule-associated proteins (MAPs) that integrate environmental and developmental stimuli to fine-tune and adjust cytoskeletal arrays. IQ67 DOMAIN (IQD) proteins recently emerged as a class of plant-specific MAPs with largely unknown functions. Here, using a reverse genetics approach, we characterize Arabidopsis IQD5 in terms of its expression domains, subcellular localization, and biological functions. We show that IQD5 is expressed mostly in vegetative tissues, where it localizes to cortical microtubule arrays. Our phenotypic analysis of iqd5 loss-of-function lines reveals functions of IQD5 in pavement cell (PC) shape morphogenesis. Histochemical analysis of cell wall composition further suggests reduced rates of cellulose deposition in anticlinal cell walls, which correlate with reduced anisotropic expansion. Lastly, we demonstrate IQD5-dependent recruitment of calmodulin calcium sensors to cortical microtubule arrays and provide first evidence for important roles for calcium in regulation of PC morphogenesis. Our work identifies IQD5 as a novel player in PC shape regulation and, for the first time, links calcium signaling to developmental processes that regulate anisotropic growth in PCs.