Organ-specific alternative transcripts of KNOX family class 2 homeobox genes of rice

KNOX II transcription factor HOS59 functions in regulating rice grain size

Plant Knotted1-like homeobox (KNOX) genes encode homeodomain-containing transcription factors. In rice (Oryza sativa L.), little is known about the downstream target genes of KNOX Class II subfamily proteins. Here we generated chromatin immunoprecipitation (ChIP)-sequencing datasets for HOS59, a member of the rice KNOX Class II subfamily, and characterized the genome-wide binding sites of HOS59. We conducted trait ontology (TO) analysis of 9705 identified downstream target genes, and found that multiple TO terms are related to plant structure morphology and stress traits. ChIP-quantitative PCR (qPCR) was conducted to validate some key target genes. Meanwhile, our IP-MS datasets showed that HOS59 was closely associated with BELL family proteins, some grain size regulators (OsSPL13, OsSPL16, OsSPL18, SLG, etc.), and some epigenetic modification factors such as OsAGO4α and OsAGO4β, proteins involved in small interfering RNA-mediated gene silencing. Furthermore, we employed CRISPR/Cas9 editing and transgenic approaches to generate hos59 mutants and overexpression lines, respectively. Compared with wild-type plants, the hos59 mutants have longer grains and increased glume cell length, a loose plant architecture, and drooping leaves, while the overexpression lines showed smaller grain size, erect leaves, and lower plant height. The qRT-PCR results showed that mutation of the HOS59 gene led to upregulation of some grain size-related genes such as OsSPL13, OsSPL18, and PGL2. In summary, our results indicate that HOS59 may be a repressor of the downstream target genes, negatively regulating glume cell length, rice grain size, plant architecture, etc. The identified downstream target genes and possible interaction proteins of HOS59 improve our understanding of the KNOX regulatory networks.

Exploring genetic variation among Jordanian Solanum lycopersicon L. landraces and their performance under salt stress using SSR markers

Background

Tomatoes (Solanum lycopersicon L.) are one of the main daily consumed vegetables in the human diet. Tomato has been classified as moderately sensitive to salinity at most stages of plant development, including seed germination, seedling (vegetative), and reproduction phases. In this study, we evaluated the performance and response of 39 tomato landraces from Jordan under salt stress conditions. Furthermore, the landraces were also genetically characterized using simple sequence repeat (SSR) markers.

Results

The studied morphological-related traits at the seedling stage were highly varied among landraces of which the landrace number 24 (Jo970) showed the best performance with the highest salt tolerance. The total number of amplification products produced by five primers (LEaat002, LEaat006, LEaat008, LEga003, LEta019) was 346 alleles. Primer LEta 019 produced the highest number of alleles (134) and generated the highest degree of polymorphism (100%) among landraces in addition to primers (LEaat002, LEaat006, LEaat008). The lowest dissimilarity among landraces ranged from 0.04 between accessions 25 (Jo969) and 26 (Jo981) and the highest dissimilarity (1.45) was found between accessions 39 (Jo980) and both 3 (Jo960) and 23 (Jo978). The dendrogram showed two main clusters and separated 30 landraces from the rest 9 landraces. High genetic diversity was detected (0.998) based on the average polymorphism information. Therefore, the used SSRs in the current study provide new insights to reveal the genetic variation among thirty-nine Jordanian tomato landraces. According to functional annotations of the gene-associated SSRs in tomatoes, a few of SSR markers gene-associated markers, for example, LEaat002 and LEaat008 markers are related to MEIS1 Transcription factors genes (Solyc07g007120 and Solyc07g007120.2). The LEaat006 is related to trypsin and protease inhibitor (Kunitz_legume) gene (Solyc03g020010). Also, the SSR LEga003 marker was related to the Carbonic anhydrase gene (Solyc09g010970).

Conclusions

The genetic variation of tomato landraces could be used for considering salt tolerance improvement in tomato breeding programs. 

Supplementary Information

The online version contains supplementary material available at 10.1186/s43141-022-00327-2.

Genome-Wide Identification and Characterization of TALE Superfamily Genes in Soybean (Glycine max L.)

The three-amino-acid-loop-extension (TALE) superfamily genes broadly existed in plants, which played important roles in plant growth, development and abiotic stress responses. In this study, we identified 68 Glycine max TALE (GmTALE) superfamily members. Phylogenetic analysis divided the GmTALE superfamily into the BEL1-like (BLH/BELL homeodomain) and the KNOX (KNOTTED-like homeodomain) subfamilies. Moreover, the KNOX subfamily could be further categorized into three clades (KNOX Class I, KNOX Class II and KNOX Class III). The GmTALE genes showed similarities in the gene structures in the same subfamily or clade, whose coding proteins exhibited analogous motif and conserved domain compositions. Besides, synteny analyses and evolutionary constraint evaluations of the TALE members among soybean and different species provided more clues for GmTALE superfamily evolution. The cis-element analyses in gene promoter regions and relevant gene expression profiling revealed different regulating roles of GmTALE genes during soybean plant development, saline and dehydration stresses. Genome-wide characterization, evolution, and expression profile analyses of GmTALE genes can pave the way for future gene functional research and facilitate their roles for applications in genetic improvement on soybean in saline and dehydration stresses.

Rice Ovate Family Protein 2 (OFP2) alters hormonal homeostasis and vasculature development

OFP (Ovate Family Protein) is a transcription factor family found only in plants. In dicots, OFPs control fruit shape and secondary cell wall biosynthesis. OFPs are also thought to function through interactions with KNOX and BELL transcription factors. Here, we have functionally characterized OsOFP2, a member of the OFP subgroup associated with regulating fruit shape. OsOFP2 was found to localize to the nucleus and to the cytosol. A putative nuclear export signal was identified within the OVATE domain and was required for the localization of OsOFP2 to distinct cytosolic spots. Rice plants overexpressing OsOFP2 were reduced in height and exhibited altered leaf morphology, seed shape, and positioning of vascular bundles in stems. Transcriptome analysis indicated disruptions of genes associated with vasculature development, lignin biosynthesis, and hormone homeostasis. Reduced expression of the gibberellin biosynthesis gene GA 20-oxidase 7 coincided with lower gibberellin content in OsOFP2 overexpression lines. Also, we found that OsOFP2 was expressed in plant vasculature and determined that putative vascular development KNOX and BELL proteins interact with OsOFP2. KNOX and BELL genes are known to suppress gibberellin biosynthesis through GA20ox gene regulation and can restrict lignin biosynthesis. We propose that OsOFP2 could modulate KNOX-BELL function to control diverse aspects of development including vasculature development.

Intrinsically disordered proteins and multicellular organisms

Intrinsically disordered proteins (IDPs) and IDP regions lack stable tertiary structure yet carry out numerous biological functions, especially those associated with signaling, transcription regulation, DNA condensation, cell division, and cellular differentiation. Both post-translational modifications (PTMs) and alternative splicing (AS) expand the functional repertoire of IDPs. Here we propose that an "IDP-based developmental toolkit," which is comprised of IDP regions, PTMs, especially multiple PTMs, within these IDP regions, and AS events within segments of pre-mRNA that code for these same IDP regions, allows functional diversification and environmental responsiveness for molecules that direct the development of complex metazoans.

Deep transcriptome sequencing of rhizome and aerial-shoot in Sorghum propinquum

Transcriptomic data for Sorghum propinquum, the wild-type sorghum, are limited in public databases. S. propinquum has a subterranean rhizome and transcriptome data will help in understanding the molecular mechanisms underlying rhizome formation. We sequenced the transcriptome of S. propinquum aerial-shoot and rhizome using an Illumina platform. More than 70 % of the genes in the S. propinquum genome were expressed in aerial-shoot and rhizome. The expression patterns of 1963 and 599 genes, including transcription factors, were specific or enriched in aerial-shoot and rhizome respectively, indicating their possible roles in physiological processes in these tissues. Comparative analysis revealed several cis-elements, ACGT box, GCCAC, GATC and TGACG box, which showed significantly higher abundance in aerial-shoot-specific genes. In rhizome-specific genes MYB and ROOTMOTIFTAPOX1 motifs, and 10 promoter and cytokinin-responsive elements were highly enriched. Of the S. propinquum genes, 27.9 % were identified as alternatively spliced and about 60 % of the alternative splicing (AS) events were tissue-specific, suggesting that AS played a crucial role in determining tissue-specific cellular function. The transcriptome data, especially the co-localized rhizome-enriched expressed transcripts that mapped to the publicly available rhizome-related quantitative trait loci, will contribute to gene discovery in S. propinquum and to functional studies of the sorghum genome. Deep transcriptome sequencing revealed a clear difference in the expression patterns of genes between aerial-shoot and rhizome in S. propinquum. This data set provides essential information for future studies into the molecular genetic mechanisms involved in rhizome formation.

Characterization of the variable 3' UTR and expression of the two intron-containing KIN transcripts from Capsella bursa-pastoris

KIN genes are crucial members of the cold-regulated (COR) gene family, and are exclusively involved in normal developmental processes in many organs and respond to a variety of abiotic stresses in plants. Here, we cloned and sequenced not only two completely-spliced KIN transcripts (CbKIN1-S and CbKIN2-S), but also two intron-containing KIN transcripts (CbKIN1-U and CbKIN2-U), from Capsella bursa-pastoris, a widespread plant of the Brassicaceae family. The CbKIN1-U and CbKIN2-U transcripts each contained one additional intron in the coding region compared to the corresponding CbKIN1-S and CbKIN2-S transcripts. In addition, the two intron-containing KIN transcripts were found by rapid amplification of cDNA 3' ends (3' RACE) analysis with specific primers to have variable 3' untranslated regions (3' UTRs). We also analyzed CbKIN1-U and CbKIN2-U levels in different organs and embryonic stages by quantitative polymerase chain reaction (qPCR). They were found to be expressed in middle-stage embryos and flowers. After abscisic acid (ABA) treatment, CbKIN1-U and CbKIN2-U showed strong responses in young leaves and weak responses in flowers. Levels of the two intron-containing KIN transcripts were markedly increased in young leaves when plants were exposed to cold and heat stress. Both of them showed stronger responses to ABA treatment and cold stress than that to heat stress. CbKIN1-U and CbKIN2-U share similar gene expression profiles in development and in response to exposure to different stresses, suggesting that they probably play similar biological roles in C. bursa-pastoris.

Alternative splicing studies of the reactive oxygen species gene network in Populus reveal two isoforms of high-isoelectric-point superoxide dismutase

Recent evidence has shown that alternative splicing (AS) is widely involved in the regulation of gene expression, substantially extending the diversity of numerous proteins. In this study, a subset of expressed sequence tags representing members of the reactive oxygen species gene network was selected from the PopulusDB database to investigate AS mechanisms in Populus. Examples of all known types of AS were detected, but intron retention was the most common. Interestingly, the closest Arabidopsis (Arabidopsis thaliana) homologs of half of the AS genes identified in Populus are not reportedly alternatively spliced. Two genes encoding the protein of most interest in our study (high-isoelectric-point superoxide dismutase [hipI-SOD]) have been found in black cottonwood (Populus trichocarpa), designated PthipI-SODC1 and PthipI-SODC2. Analysis of the expressed sequence tag libraries has indicated the presence of two transcripts of PthipI-SODC1 (hipI-SODC1b and hipI-SODC1s). Alignment of these sequences with the PthipI-SODC1 gene showed that hipI-SODC1b was 69 bp longer than hipI-SODC1s due to an AS event involving the use of an alternative donor splice site in the sixth intron. Transcript analysis showed that the splice variant hipI-SODC1b was differentially expressed, being clearly expressed in cambial and xylem, but not phloem, regions. In addition, immunolocalization and mass spectrometric data confirmed the presence of hipI-SOD proteins in vascular tissue. The functionalities of the spliced gene products were assessed by expressing recombinant hipI-SOD proteins and in vitro SOD activity assays.

Genome-wide identification, classification, evolutionary expansion and expression analyses of homeobox genes in rice

Homeobox genes play a critical role in regulating various aspects of plant growth and development. In the present study, we identified a total of 107 homeobox genes in the rice genome and grouped them into ten distinct subfamilies based upon their domain composition and phylogenetic analysis. A significantly large number of homeobox genes are located in the duplicated segments of the rice genome, which suggests that the expansion of homeobox gene family, in large part, might have occurred due to segmental duplications in rice. Furthermore, microarray analysis was performed to elucidate the expression profiles of these genes in different tissues and during various stages of vegetative and reproductive development. Several genes with predominant expression during various stages of panicle and seed development were identified. At least 37 homeobox genes were found to be differentially expressed significantly (more than two-fold; P < 0.05) under various abiotic stress conditions. The results of the study suggest a critical role of homeobox genes in reproductive development and abiotic stress signaling in rice, and will facilitate the selection of candidate genes of agronomic importance for functional validation.

Downstream promoter sequence of an Indian isolate of Rice tungro bacilliform virus alters tissue-specific expression in host rice and acts differentially in heterologous system

An Indian isolate of Rice tungro bacilliform virus from West Bengal (RTBV-WB) showed significant nucleotide differences in its putative promoter region when compared with a previously characterized isolate from Philippines. The transcription start site of RTBV-WB was mapped followed by assessing the activity and tissue-specificity of the full-length (FL) promoter (-231 to +645) and several of its upstream and downstream deletions by studying the expression of beta-Glucuronidase (GUS) reporter gene in transgenic rice (Oryza sativa L. subsp. indica) plants at various stages of development. In addition to the expected vascular-specific expression pattern, studied by histochemical staining, GUS enzymatic assay and northern and RT-PCR analysis, two novel patterns were revealed in some of the downstream deleted versions; a non-expressing type, representing no expression at any stage in any tissue and constitutive type, representing constitutive expression at all stages in most tissues. This indicated the presence of previously unreported positive and negative cis-regulatory elements in the downstream region. The negative element and a putative enhancer region in the upstream region specifically bound to rice nuclear proteins in vitro. The FL and its deletion derivatives were also active in heterologous systems like tobacco (Nicotiana tabacum) and wheat (Triticum durum). Expression patterns in tobacco were different from those observed in rice suggesting the importance of upstream elements in those systems and host-specific regulation of the promoter in diverse organisms. Thus, the RTBV-WB FL promoter and its derivatives contain an array of cis-elements, which control constitutive or tissue- and development-specific gene expression in a combinatorial fashion.

Differential transcription initiation and alternative RNA splicing of Knox7, a class 2 homeobox gene of maize

Knox7, a class 2 homeobox gene has been characterized in maize. A combination of experimental (3'- and 5'-RACE) and bioinformatics approaches supported the idea that Knox7 would be transcribed into two alternative transcripts by differential initiation of transcription. Sequence differences between alternative transcripts, Knox7L the larger and Knox7S the smaller, were confined to their 5' end regions and exon 1 was only found in Knox7L transcripts. Deduced proteins shared the same homeodomain, while an Ala and Ala/Gly rich domain was found only in KNOX7L protein. We hypothesize that KNOX7L and KNOX7S might regulate (differentially) the expression of the same gene(s) by binding competitively to the same cis-acting element(s). Further expression analysis using RT-PCR to amplify cDNA portions corresponding to ORFs of both Knox7 alternative transcripts showed that seven cDNA clones were probably generated by alternative splicing of Knox7L. Alignment of these sequences showed that they are in frame suggesting the existence of the corresponding proteins. Quantitative RT-PCR experiments indicated that Knox7S and Knox7L were expressed in maize embryos during germination. In the same tissue, expression of Knox7S was stimulated by light and ABA and inhibited by GA, two hormones that control germination process.

EU-OSTID: a collection of transposon insertional mutants for functional genomics in rice

A collection of 1373 unique flanking sequence tags (FSTs), generated from Ac/Ds and Ac transposon lines for reverse genetics studies, were produced in japonica and indica rice, respectively. The Ds and Ac FSTs together with the original T-DNAs were assigned a position in the rice genome sequence represented as assembled pseudomolecules, and found to be distributed evenly over the entire rice genome with a distinct bias for predicted gene-rich regions. The bias of the Ds and Ac transposon inserts for genes was exemplified by the presence of 59% of the inserts in genes annotated on the rice chromosomes and 41% present in genes transcribed as disclosed by their homology to cDNA clones. In a screen for inserts in a set of 75 well annotated transcription factors, including homeobox-containing genes, we found six Ac/Ds inserts. This high frequency of Ds and Ac inserts in genes suggests that saturated knockout mutagenesis in rice using this strategy will be efficient and possible with a lower number of inserts than expected. These FSTs and the corresponding plant lines are publicly available through OrygenesDB database and from the EU consortium members.

Radial axis differentiation in a globular embryo is marked by HAZ1, a PHD-finger homeobox gene of rice

Homeobox genes that encode transcription factors play important roles in development and differentiation of both plant and animal systems. From a cDNA library of 3-day after-pollination (DAP) rice embryos we cloned a HAZ1 cDNA that encodes a protein with a PHD-finger domain and a homeodomain. A database search showed that HAZ1 was most similar in its entire amino acid sequence to Zmhox1a (52% identity) and Zmhox1b (50%), PHD-finger family homeodomain proteins of maize. Differing from Zmhox1, overexpression of HAZ1 brought no morphological change either in tobacco or in rice. In situ hybridization showed that HAZ1 was expressed at a higher level in the outer layers of a developing embryo than in the inner parts of the embryo at 3 DAP. At 4 and 5 DAPs, the expression of HAZ1 was concentrated at the ventral part of an embryo. These results indicate that HAZ1 marks outer layer cells of a globular embryo before any morphological differentiation is discerned in it. Radial axis differentiation marked by HAZ1 is then collapsed dynamically along with embryo morphogenesis, and HAZ1 later marks the ventral surface of the embryo.

Conservation and molecular dissection of ROUGH SHEATH2 and ASYMMETRIC LEAVES1 function in leaf development

Maize ROUGH SHEATH2 (RS2) and Arabidopsis ASYMMETRIC LEAVES1 (AS1) are orthologous Myb-related genes required for leaf development and act as negative regulators of class 1 KNOTTED1-like homeobox (KNOX) genes in leaf primordia. Expression of RS2 in Arabidopsis fully complements as1 leaf phenotypes and represses the expression of the KNOX gene KNAT1 in leaves. Whereas loss of AS1 function in Arabidopsis results in rounded, lobed leaves with shorter and wider petioles, overexpression of either RS2 or AS1 results in longer and narrower leaves with longer petioles than wild type. A conserved C-terminal domain (CTD) mediates homodimerization of both RS2 and AS1 and modulates leaf shape when expressed independently of the Myb domain in Arabidopsis. Homodimerization is not absolutely required for KNAT1 repression. RS2:GFP fusion protein is biologically active, localized in discrete dynamic subnuclear foci and associates with DNA during cell division.