The RNA helicase, eIF4A-1, is required for ovule development and cell size homeostasis in Arabidopsis

Deciphering the possible role of RNA-helicase genes mechanism in response to abiotic stresses in rapeseed (Brassica napus L.)

BACKGROUND

Plants mediate several defense mechanisms to withstand abiotic stresses. Several gene families respond to stress as well as multiple transcription factors to minimize abiotic stresses without minimizing their effects on performance potential. RNA helicase (RH) is one of the foremost critical gene families that can play an influential role in tolerating abiotic stresses in plants. However, little knowledge is present about this protein family in rapeseed (canola). Here, we performed a comprehensive survey analysis of the RH protein family in rapeseed (Brassica napus L.).

RESULTS

A total of 133 BnRHs genes have been discovered in this study. By phylogenetic analysis, RHs genes were divided into one main group and a subgroup. Examination of the chromosomal position of the identified genes showed that most of the genes (27%) were located on chromosome 3. All 133 identified sequences contained the main DEXDC domain, the HELICC domain, and a number of sub-domains. The results of biological process studies showed that about 17% of the proteins acted as RHs, 22% as ATP binding, and 14% as mRNA binding. Each part of the conserved motifs, communication network, and three-dimensional structure of the proteins were examined separately. The results showed that the RWC in leaf tissue decreased with higher levels of drought stress and in both root and leaf tissues sodium concentration was increased upon increased levels of salt stress treatments. The proline content were found to be increased in leaf and root with the increased level of stress treatment. Finally, the expression patterns of eight selected RHs genes that have been exposed to drought, salinity, cold, heat and cadmium stresses were investigated by qPCR. The results showed the effect of genes under stress. Examination of gene expression in the Hayola #4815 cultivar showed that all primers except primer #79 had less expression in both leaves and roots than the control level.

CONCLUSIONS

New finding from the study have been presented new insights for better understanding the function and possible mechanism of RH in response to abiotic stress in rapeseed.

Transcriptome-wide association study coupled with eQTL analysis reveals the genetic connection between gene expression and flowering time in Arabidopsis

Genome-wide association study (GWAS) has improved our understanding of complex traits, but challenges exist in distinguishing causation versus association caused by linkage disequilibrium. Instead, transcriptome-wide association studies (TWAS) detect direct associations between expression levels and phenotypic variations, providing an opportunity to better prioritize candidate genes. To assess the feasibility of TWAS, we investigated the association between transcriptomes, genomes, and various traits in Arabidopsis, including flowering time. The associated genes formerly known to regulate growth allometry or metabolite production were first identified by TWAS. Next, for flowering time, six TWAS-newly identified genes were functionally validated. Analysis of the expression quantitative trait locus (eQTL) further revealed a trans-regulatory hotspot affecting the expression of several TWAS-identified genes. The hotspot covers the FRIGIDA (FRI) gene body, which possesses multiple haplotypes differentially affecting the expression of downstream genes, such as FLOWERING LOCUS C (FLC) and SUPPRESSOR OF OVEREXPRESSION OF CO 1 (SOC1). We also revealed multiple independent paths towards the loss of function of FRI in natural accessions. Altogether, this study demonstrates the potential of combining TWAS with eQTL analysis to identify important regulatory modules of FRI-FLC-SOC1 for quantitative traits in natural populations.

Protein interactome of 3',5'-cAMP reveals its role in regulating the actin cytoskeleton

Identification of protein interactors is ideally suited for the functional characterization of small molecules. 3',5'-cAMP is an evolutionary ancient signaling metabolite largely uncharacterized in plants. To tap into the physiological roles of 3',5'-cAMP, we used a chemo-proteomics approach, thermal proteome profiling (TPP), for the unbiased identification of 3',5'-cAMP protein targets. TPP measures shifts in the protein thermal stability upon ligand binding. Comprehensive proteomics analysis yielded a list of 51 proteins significantly altered in their thermal stability upon incubation with 3',5'-cAMP. The list contained metabolic enzymes, ribosomal subunits, translation initiation factors, and proteins associated with the regulation of plant growth such as CELL DIVISION CYCLE 48. To functionally validate obtained results, we focused on the role of 3',5'-cAMP in regulating the actin cytoskeleton suggested by the presence of actin among the 51 identified proteins. 3',5'-cAMP supplementation affected actin organization by inducing actin-bundling. Consistent with these results, the increase in 3',5'-cAMP levels, obtained either by feeding or by chemical modulation of 3',5'-cAMP metabolism, was sufficient to partially rescue the short hypocotyl phenotype of the actin2 actin7 mutant, severely compromised in actin level. The observed rescue was specific to 3',5'-cAMP, as demonstrated using a positional isomer 2',3'-cAMP, and true for the nanomolar 3',5'-cAMP concentrations reported for plant cells. In vitro characterization of the 3',5'-cAMP-actin pairing argues against a direct interaction between actin and 3',5'-cAMP. Alternative mechanisms by which 3',5'-cAMP would affect actin dynamics, such as by interfering with calcium signaling, are discussed. In summary, our work provides a specific resource, 3',5'-cAMP interactome, as well as functional insight into 3',5'-cAMP-mediated regulation in plants.

Plant phase extraction: A method for enhanced discovery of the RNA-binding proteome and its dynamics in plants

RNA-binding proteins (RBPs) play critical roles in posttranscriptional gene regulation. Current methods of systematically profiling RBPs in plants have been predominantly limited to proteins interacting with polyadenylated (poly(A)) RNAs. We developed a method called plant phase extraction (PPE), which yielded a highly comprehensive RNA-binding proteome (RBPome), uncovering 2,517 RBPs from Arabidopsis (Arabidopsis thaliana) leaf and root samples with a highly diverse array of RNA-binding domains. We identified traditional RBPs that participate in various aspects of RNA metabolism and a plethora of nonclassical proteins moonlighting as RBPs. We uncovered constitutive and tissue-specific RBPs essential for normal development and, more importantly, revealed RBPs crucial for salinity stress responses from a RBP-RNA dynamics perspective. Remarkably, 40% of the RBPs are non-poly(A) RBPs that were not previously annotated as RBPs, signifying the advantage of PPE in unbiasedly retrieving RBPs. We propose that intrinsically disordered regions contribute to their nonclassical binding and provide evidence that enzymatic domains from metabolic enzymes have additional roles in RNA binding. Taken together, our findings demonstrate that PPE is an impactful approach for identifying RBPs from complex plant tissues and pave the way for investigating RBP functions under different physiological and stress conditions at the posttranscriptional level.

Functions and mechanisms of RNA helicases in plants

RNA helicases (RHs) are a family of ubiquitous enzymes that alter RNA structures and remodel ribonucleoprotein complexes typically using energy from the hydrolysis of ATP. RHs are involved in various aspects of RNA processing and metabolism, exemplified by transcriptional regulation, pre-mRNA splicing, miRNA biogenesis, liquid-liquid phase separation, and rRNA biogenesis, among other molecular processes. Through these mechanisms, RHs contribute to vegetative and reproductive growth, as well as abiotic and biotic stress responses throughout the life cycle in plants. In this review, we systematically characterize RH-featured domains and signature motifs in Arabidopsis. We also summarize the functions and mechanisms of RHs in various biological processes in plants with a focus on DEAD-box and DEAH-box RNA helicases, aiming to present the latest understanding of RHs in plant biology.

LkARF7 and LkARF19 overexpression promote adventitious root formation in a heterologous poplar model by positively regulating LkBBM1

Cuttage propagation involves adventitious root formation induced by auxin. In our previous study, Larix kaempferi BABY BOOM 1 (LkBBM1), which is known to regulate adventitious root formation, was affected by auxin. However, the relationship between LkBBM1 and auxin remains unclear. Auxin response factors (ARFs) are a class of important transcription factors in the auxin signaling pathway and modulate the expression of early auxin-responsive genes by binding to auxin response elements. In the present study, we identified 14 L. kaempferi ARFs (LkARFs), and found LkARF7 and LkARF19 bound to LkBBM1 promoter and enhanced its transcription using yeast one-hybrid, ChIP-qPCR, and dual-luciferase assays. In addition, the treatment with naphthalene acetic acid promoted the expression of LkARF7 and LkARF19. We also found that overexpression of these two genes in poplar promoted adventitious root formation. Furthermore, LkARF19 interacted with the DEAD-box ATP-dependent RNA helicase 53-like protein to form a heterodimer to regulate adventitious root formation. Altogether, our results reveal an additional regulatory mechanism underlying the control of adventitious root formation by auxin.

Integrated multi-locus genome-wide association studies and transcriptome analysis for seed yield and yield-related traits in Brassica napus

Rapeseed (Brassica napus L.), the third largest oil crop, is an important source of vegetable oil and biofuel for the world. Although the breeding and yield has been improved, rapeseed still has the lowest yield compared with other major crops. Thus, increasing rapeseed yield is essential for the high demand of vegetable oil and high-quality protein for live stocks. Silique number per plant (SN), seed per pod (SP), and 1000-seed weight (SW) are the three important factors for seed yield in rapeseed. Some yield-related traits, including plant height (PH), flowering time (FT), primary branch number (BN) and silique number per inflorescence (SI) also affect the yield per plant (YP). Using six multi-locus genome-wide association study (ML-GWAS) approaches, a total of 908 yield-related quantitative trait nucleotides (QTNs) were identified in a panel consisting of 403 rapeseed core accessions based on whole-genome sequencing. Integration of ML-GWAS with transcriptome analysis, 79 candidate genes, including BnaA09g39790D (RNA helicase), BnaA09g39950D (Lipase) and BnaC09g25980D (SWEET7), were further identified and twelve genes were validated by qRT-PCRs to affect the SW or SP in rapeseed. The distribution of superior alleles from nineteen stable QTNs in 20 elite rapeseed accessions suggested that the high-yielding accessions contained more superior alleles. These results would contribute to a further understanding of the genetic basis of yield-related traits and could be used for crop improvement in B. napus.

Genetic dissection of the natural variation of ovule number per ovary in oilseed rape germplasm (Brassica napus L.)

Oilseed rape is one of the world's largest oil and industrial crops, providing humans with various products, such as vegetable oil and biofuel. Ovules are the direct precursors of seeds, and ovule number per ovary (ONPO) largely determines seed number per fruit that affects both yield and fitness of seed crops. The ONPO shows wide variation in oilseed rape, whereas the underlying genes and mechanisms are poorly known. The present study performed the genetic, physiological and transcriptomic analyses of ovule number per ovary using an association panel and the extreme lines. The ONPO of 327 accessions planted in four environments showed a large variation from 19.2 to 43.8, indicating a great potential for the further genetic improvement of ovule number. The genome-wide association study (GWAS) identified a total of 43 significant SNP markers. Further, these SNPs were integrated into 18 association loci, which were distributed on chromosomes A01, A03, A06, A07, A09, C01, C03, C06, C07, and C09, explaining 4.3-11.5% of the phenotypic variance. The ONPO decreased as their appearance order on the inflorescence and was associated with the level of several types of endogenous phytohormones but not related to leaf area and photosynthetic rate. Comparative transcriptomic analysis identified a total of 4,449 DEGs enriched in 30 classes, including DNA, RNA, protein, signaling, transport, development, cell wall, lipid metabolism, and secondary metabolism. Nearly half of DEGs were involved in the known pathways in regulating ovule number, of which 12 were homologous to know ovule number regulating genes, indicating a strong link between the identified DEGs and ovule number. A total of 73 DEGs were located within the genomic regions of association loci, of which six were identified as candidates based on functional annotation. These results provide useful information for the further genetic improvement of ovule and seed number in oilseed rape.

Complex developmental and transcriptional dynamics underlie pollinator-driven evolutionary transitions in nectar spur morphology in Aquilegia (columbine)

PREMISE

Determining the developmental programs underlying morphological variation is key to elucidating the evolutionary processes that generated the stunning biodiversity of the angiosperms. Here, we characterized the developmental and transcriptional dynamics of the elaborate petal nectar spur of Aquilegia (columbine) in species with contrasting pollination syndromes and spur morphologies.

METHODS

We collected petal epidermal cell number and length data across four Aquilegia species, two with short, curved nectar spurs of the bee-pollination syndrome and two with long, straight spurs of the hummingbird-pollination syndrome. We also performed RNA-seq on A. brevistyla (bee) and A. canadensis (hummingbird) distal and proximal spur compartments at multiple developmental stages. Finally, we intersected these data sets with a previous QTL mapping study on spur length and shape to identify new candidate loci.

RESULTS

The differential growth between the proximal and distal surfaces of curved spurs is primarily driven by differential cell division. However, independent transitions to straight spurs in the hummingbird syndrome have evolved by increasing differential cell elongation between spur surfaces. The RNA-seq data reveal these tissues to be transcriptionally distinct and point to auxin signaling as being involved with the differential cell elongation responsible for the evolution of straight spurs. We identify several promising candidate genes for future study.

CONCLUSIONS

Our study, taken together with previous work in Aquilegia, reveals the complexity of the developmental mechanisms underlying trait variation in this system. The framework we established here will lead to exciting future work examining candidate genes and processes involved in the rapid radiation of the genus.

Recurrent Plant-Specific Duplications of KNL2 and Its Conserved Function as a Kinetochore Assembly Factor

KINETOCHORE NULL2 (KNL2) plays key role in the recognition of centromeres and new CENH3 deposition. To gain insight into the origin and diversification of the KNL2 gene, we reconstructed its evolutionary history in the plant kingdom. Our results indicate that the KNL2 gene in plants underwent three independent ancient duplications in ferns, grasses and eudicots. Additionally, we demonstrated that previously unclassified KNL2 genes could be divided into two clades αKNL2 and βKNL2 in eudicots and γKNL2 and δKNL2 in grasses, respectively. KNL2s of all clades encode the conserved SANTA domain, but only the αKNL2 and γKNL2 groups additionally encode the CENPC-k motif. In the more numerous eudicot sequences, signatures of positive selection were found in both αKNL2 and βKNL2 clades, suggesting recent or ongoing adaptation. The confirmed centromeric localization of βKNL2 and mutant analysis suggests that it participates in loading of new CENH3, similarly to αKNL2. A high rate of seed abortion was found in heterozygous βKNL2 plants and the germinated homozygous mutants did not develop beyond the seedling stage. Taken together, our study provides a new understanding of the evolutionary diversification of the plant kinetochore assembly gene KNL2, and suggests that the plant-specific duplicated KNL2 genes are involved in centromere and/or kinetochore assembly for preserving genome stability.

Genome-Wide Analysis of DEAD-box RNA Helicase Family in Wheat (Triticum aestivum) and Functional Identification of TaDEAD-box57 in Abiotic Stress Responses

DEAD-box RNA helicases constitute the largest subfamily of RNA helicase superfamily 2 (SF2), and play crucial roles in plant growth, development, and abiotic stress responses. Wheat is one of the most important cereal crops in worldwide, and abiotic stresses greatly restrict its production. So far, the DEAD-box RNA helicase family has yet to be characterized in wheat. Here, we performed a comprehensive genome-wide analysis of the DEAD-box RNA helicase family in wheat, including phylogenetic relationships, chromosomal distribution, duplication events, and protein motifs. A total of 141 TaDEAD-box genes were identified and found to be unevenly distributed across all 21 chromosomes. Whole genome/segmental duplication was identified as the likely main driving factor for expansion of the TaDEAD-box family. Expression patterns of the 141 TaDEAD-box genes were compared across different tissues and under abiotic stresses to identify genes to be important in growth or stress responses. TaDEAD-box57-3B was significantly up-regulated under multiple abiotic stresses, and was therefore selected for further analysis. TaDEAD-box57-3B was localized to the cytoplasm and plasma membrane. Ectopic expression of TaDEAD-box57-3B in Arabidopsis improved tolerance to drought and salt stress as measured by germination rates, root lengths, fresh weights, and survival rates. Transgenic lines also showed higher levels of proline and chlorophyll and lower levels of malonaldehyde (MDA) than WT plants in response to drought or salt stress. In response to cold stress, the transgenic lines showed significantly better growth and higher survival rates than WT plants. These results indicate that TaDEAD-box57-3B may increase tolerance to drought, salt, and cold stress in transgenic plants through regulating the degree of membrane lipid peroxidation. This study provides new insights for understanding evolution and function in the TaDEAD-box gene family.

Molecular Network for Regulation of Ovule Number in Plants

In seed-bearing plants, the ovule ("small egg") is the organ within the gynoecium that develops into a seed after fertilization. The gynoecium located in the inner compartment of the flower turns into a fruit. The number of ovules in the ovary determines the upper limit or the potential of seed number per fruit in plants, greatly affecting the final seed yield. Ovule number is an important adaptive characteristic for plant evolution and an agronomic trait for crop improvement. Therefore, understanding the mechanism and pathways of ovule number regulation becomes a significant research aspect in plant science. This review summarizes the ovule number regulators and their regulatory mechanisms and pathways. Specially, an integrated molecular network for ovule number regulation is constructed, in which phytohormones played a central role, followed by transcription factors, enzymes, other protein and micro-RNA. Of them, AUX, BR and CK are positive regulator of ovule number, whereas GA acts negatively on it. Interestingly, many ovule number regulators have conserved functions across several plant taxa, which should be the targets of genetic improvement via breeding or gene editing. Many ovule number regulators identified to date are involved in the diverse biological process, such as ovule primordia formation, ovule initiation, patterning, and morphogenesis. The relations between ovule number and related characteristics/traits especially of gynoecium/fruit size, ovule fertility, and final seed number, as well as upcoming research questions, are also discussed. In summary, this review provides a general overview of the present finding in ovule number regulation, which represents a more comprehensive and in-depth cognition on it.

The nuclear GUCT domain-containing DEAD-box RNA helicases govern gametophytic and sporophytic development in Physcomitrium patens

KEY MESSAGE

In Physcomitrium patens, PpRH1/PpRH2 are GUCT-domain-containing DEAD-BOX RNA helicases localize to the nucleus. They are implicated in cell and tissue development in all stages of the moss life cycle.

ABSTRACT

The DEAD-box-containing RNA helicase family encompasses a large and functionally important group of enzymes involved in cellular processes committed to the metabolism of RNA, including its transcription, processing, transport, translation and decay. Studies indicate this protein family has implied roles in plant vegetative and reproductive developmental processes as well as response to environmental stresses such has cold and high salinity. We focus here on a small conserved sub-group of GUCT domain-containing RNA helicase in the moss Physcomitrium patens. Phylogenetic analysis shows that RNA helicases containing the GUCT domain form a distinct conserved clade across the green lineage. In this clade, the P. patens genome possesses two closely related paralogues RNA helicases predicted to be nuclear, PpRH1 and PpRH2. Using in-locus gene fluorescent tagging we show that PpRH1 is localized to the nucleus in protonema. Analysis of PpRH1 and PpRH2 deletions, individually and together, indicates their potential roles in protonema, gametophore and sporophyte cellular and tissue development in P. patens. Additionally, the ultrastructural analysis of phyllid chloroplasts in Δrh2 and Δrh1/2 shows distinct starch granule accumulation under standard growth conditions associated with changes in photosynthetic activity parameters. We could not detect effects of either temperature or stress on protonema growth or PpRH1 and PpRH2 expression. Together, these results suggest that nuclear GUCT-containing RNA helicases play a role primarily in developmental processes directly or indirectly linked to photosynthesis activity in the moss P. patens.

SUPPLEMENTARY INFORMATION

The online version contains supplementary material available at 10.1007/s11103-021-01152-w.

eIF4A, a target of siRNA derived from rice stripe virus, negatively regulates antiviral autophagy by interacting with ATG5 in Nicotiana benthamiana

Autophagy is induced by viral infection and has antiviral functions in plants, but the underlying mechanism is poorly understood. We previously identified a viral small interfering RNA (vsiRNA) derived from rice stripe virus (RSV) RNA4 that contributes to the leaf-twisting and stunting symptoms caused by this virus by targeting the host eukaryotic translation initiation factor 4A (eIF4A) mRNA for silencing. In addition, autophagy plays antiviral roles by degrading RSV p3 protein, a suppressor of RNA silencing. Here, we demonstrate that eIF4A acts as a negative regulator of autophagy in Nicotiana benthamiana. Silencing of NbeIF4A activated autophagy and inhibited RSV infection by facilitating autophagic degradation of p3. Further analysis showed that NbeIF4A interacts with NbATG5 and interferes with its interaction with ATG12. Overexpression of NbeIF4A suppressed NbATG5-activated autophagy. Moreover, expression of vsiRNA-4A, which targets NbeIF4A mRNA for cleavage, induced autophagy by silencing NbeIF4A. Finally, we demonstrate that eIF4A from rice, the natural host of RSV, also interacts with OsATG5 and suppresses OsATG5-activated autophagy, pointing to the conserved function of eIF4A as a negative regulator of antiviral autophagy. Taken together, these results reveal that eIF4A negatively regulates antiviral autophagy by interacting with ATG5 and that its mRNA is recognized by a virus-derived siRNA, resulting in its silencing, which induces autophagy against viral infection.

DeepPod: a convolutional neural network based quantification of fruit number in Arabidopsis

Background

High-throughput phenotyping based on non-destructive imaging has great potential in plant biology and breeding programs. However, efficient feature extraction and quantification from image data remains a bottleneck that needs to be addressed. Advances in sensor technology have led to the increasing use of imaging to monitor and measure a range of plants including the model Arabidopsis thaliana. These extensive datasets contain diverse trait information, but feature extraction is often still implemented using approaches requiring substantial manual input.

Results

The computational detection and segmentation of individual fruits from images is a challenging task, for which we have developed DeepPod, a patch-based 2-phase deep learning framework. The associated manual annotation task is simple and cost-effective without the need for detailed segmentation or bounding boxes. Convolutional neural networks (CNNs) are used for classifying different parts of the plant inflorescence, including the tip, base, and body of the siliques and the stem inflorescence. In a post-processing step, different parts of the same silique are joined together for silique detection and localization, whilst taking into account possible overlapping among the siliques. The proposed framework is further validated on a separate test dataset of 2,408 images. Comparisons of the CNN-based prediction with manual counting (R² = 0.90) showed the desired capability of methods for estimating silique number.

Conclusions

The DeepPod framework provides a rapid and accurate estimate of fruit number in a model system widely used by biologists to investigate many fundemental processes underlying growth and reproduction

The vascular targeted citrus FLOWERING LOCUS T3 gene promotes non-inductive early flowering in transgenic Carrizo rootstocks and grafted juvenile scions

Shortening the juvenile stage in citrus and inducing early flowering has been the focus of several citrus genetic improvement programs. FLOWERING LOCUS T (FT) is a small phloem-translocated protein that regulates precocious flowering. In this study, two populations of transgenic Carrizo citrange rootstocks expressing either Citrus clementina FT1 or FT3 genes under the control of the Arabidopsis thaliana phloem specific SUCROSE SYNTHASE 2 (AtSUC2) promoter were developed. The transgenic plants were morphologically similar to the non-transgenic controls (non-transgenic Carrizo citrange), however, only AtSUC2-CcFT3 was capable of inducing precocious flowers. The transgenic lines produced flowers 16 months after transformation and flower buds appeared 30-40 days on juvenile immature scions grafted onto transgenic rootstock. Gene expression analysis revealed that the expression of SUPPRESSOR OF OVEREXPRESSION OF CONSTANS 1 (SOC1) and APETALA1 (AP1) were enhanced in the transgenics. Transcriptome profiling of a selected transgenic line showed the induction of genes in different groups including: genes from the flowering induction pathway, APETALA2/ETHYLENE RESPONSE FACTOR (AP2/ERF) family genes, and jasmonic acid (JA) pathway genes. Altogether, our results suggested that ectopic expression of CcFT3 in phloem tissues of Carrizo citrange triggered the expression of several genes to mediate early flowering.

Diverse Allyl Glucosinolate Catabolites Independently Influence Root Growth and Development

Glucosinolates (GSLs) are sulfur-containing defense metabolites produced in the Brassicales, including the model plant Arabidopsis (Arabidopsis thaliana). Previous work suggests that specific GSLs may function as signals to provide direct feedback regulation within the plant to calibrate defense and growth. These GSLs include allyl-GSL, a defense metabolite that is one of the most widespread GSLs in Brassicaceae and has also been associated with growth inhibition. Here we show that at least three separate potential catabolic products of allyl-GSL or closely related compounds affect growth and development by altering different mechanisms influencing plant development. Two of the catabolites, raphanusamic acid and 3-butenoic acid, differentially affect processes downstream of the auxin signaling cascade. Another catabolite, acrylic acid, affects meristem development by influencing the progression of the cell cycle. These independent signaling events propagated by the different catabolites enable the plant to execute a specific response that is optimal to any given environment.

Morphological and proteomic analysis of young spikes reveals new insights into the formation of multiple-pistil wheat

A new multiple-pistil wheat mutant germplasm with more than one pistil in a floret was obtained from natural mutagenesis. This mutant can develop 2-3 grains in a glume after pollination and has a significant grain number advantage compared with normal wheat. However, the basis of the formation of multiple-pistil wheat has thus far not been well established. In this study, we first performed a continuous phenotypic observation of the floral meristem (FM) in multiple-pistil wheat. The results indicated that the secondary pistils are derived from extra stem cells that fail to terminate normally between the carpel primordium and the lodicule primordium. To further probe the potential molecular basis for the formation of secondary pistils, comparative proteomic analyses were conducted. A total of 334 differentially abundant proteins (DAPs) were identified using isobaric tags for relative and absolute quantification (iTRAQ), among which 131 proteins were highly abundant and 203 proteins were less abundant in the young spikes of multiple-pistil wheat. The DAPs, located primarily in the cell, were involved in the translation and the metabolisms of carbohydrate, nucleotide, and amino acid. Differential expression analysis showed that TaHUA2, TaRF2a, TaCHR12 and TaHEN2 may play vital roles in the regulation of wheat flower organ number. In general, the DAPs support the phenotypic analysis results at the molecular level. In combination, these results reveal new insights into the formation of multiple-pistil wheat and provide possible targets for further research on the regulation of floral organ number in wheat.

Gene Function Rather than Reproductive Mode Drives the Evolution of RNA Helicases in Sexual and Apomictic Boechera

In higher plants, sexual and asexual reproductions through seeds (apomixis) have evolved as alternative strategies. Evolutionary advantages leading to coexistence of both reproductive modes are currently not well understood. It is expected that accumulation of deleterious mutations leads to a rapid elimination of apomictic lineages from populations. In this line, apomixis originated repeatedly, likely from deregulation of the sexual pathway, leading to alterations in the development of reproductive lineages (germlines) in apomicts as compared with sexual plants. This potentially involves mutations in genes controlling reproduction.

Increasing evidence suggests that RNA helicases are crucial regulators of germline development. To gain insights into the evolution of 58 members of this diverse gene family in sexual and apomictic plants, we applied target enrichment combined with next-generation sequencing to identify allelic variants from 24 accessions of the genus Boechera, comprising sexual, facultative, and obligate apomicts. Interestingly, allelic variants from apomicts did not show consistently increased mutation frequency. Either sequences were highly conserved in any accession, or allelic variants preferentially harbored mutations in evolutionary less conserved C- and N-terminal domains, or presented high mutation load independent of the reproductive mode. Only for a few genes allelic variants harboring deleterious mutations were only identified in apomicts. To test if high sequence conservation correlates with roles in fundamental cellular or developmental processes, we analyzed Arabidopsis thaliana mutant lines in VASA-LIKE (VASL), and identified pleiotropic defects during ovule and reproductive development. This indicates that also in apomicts mechanisms of selection are in place based on gene function.

The DEAD-box RNA helicase eIF4A regulates plant development and interacts with the hnRNP LIF2L1 in Physcomitrella patens

eIF4A is a RNA-stimulated ATPase and helicase. Besides its key role in regulating cap-dependent translation initiation in eukaryotes, it also performs specific functions in regulating cell cycle progression, plant growth and abiotic stress tolerance. Flowering plants encode three eIF4A paralogues, eIF4A1, eIF4A2 and eIF4A3 that share conserved sequence motifs but differ in functions. To date, however, no information is available on eIF4A in basal land plants. In this study we report that genome of the moss Physcomitrella patens encodes multiple eIF4A genes. The encoded proteins possess the highly conserved motifs characteristic of the DEAD box helicases. Spatial expression analysis shows these genes to be ubiquitously expressed in all tissue types with Pp3c6_1080V3.1 showing high expression in filamentous protonemata. Targeted deletion of conserved core motifs in Pp3c6_1080V3.1 slowed protonemata growth and resulted in dwarfing of leafy gametophores suggesting a role for Pp3c6_1080V3.1 in regulating cell division/elongation. Rapid and strong induction of Pp3c6_1080V3.1 under salt stress and slow recovery of knockout plants upon exposure to high salt further suggest Pp3c6_1080V3.1 to be involved in stress management in P. patens. Protein-protein interaction studies that show Pp3c6_1080V3.1 to interact with the Physcomitrella heterogenous ribonucleoprotein, LIF2L1, a transcriptional regulator of stress-responsive genes in Arabidopsis. The results presented in this study provide insight into evolutionary conserved functions of eIF4A and shed light on the novel link between eIF4A activities and stress mitigation pathways/RNA metabolic processes in P. patens.

NHX1 and eIF4A1-stacked transgenic sweetpotato shows enhanced tolerance to drought stress

Co-expression of Na⁺/H⁺ antiporter NHX1 and DEAD-box RNA helicase eIF4A1 from Arabidopsis positively regulates drought stress tolerance by improving ROS scavenging capacity and maintaining membrane integrity in sweetpotato. Plants evolve multiple strategies for stress adaptation in nature. To improve sweetpotato resistance to drought stress, transgenic sweetpotato plants overexpressing the Arabidopsis Na⁺/H⁺ antiporter, NHX1, and the translation initiation factor elF4A1 were characterized for phenotypic traits and physiological performance. Without drought treatment, the NHX1-elF4A1 stacked lines (NE lines) showed normal, vigorous growth comparable to the WT plants. The NE plants showed dense green foliage with delayed leaf senescence and developed more roots than WT plants under drought treatment for 18 days. Compared to WT plants, higher level of reactive oxygen scavenging capacity was detected in NE lines as indicated by reduced H₂O₂ accumulation as well as increased superoxide dismutase activity and proline content. The relative ion leakage and malondialdehyde content were reduced in NE plants, indicating improved maintenance of intact membranes system. Both NE plants and NHX1-overexpressing plants (N lines) showed larger aerial parts and well-developed root system compared to WT plants under the drought stress conditions, likely due to the improved antioxidant capacity. The NE plants showed better ROS scavenging than N-line plants. All N- and NE-line plants produced normal storage roots with similar yields as WT in the field under normal growth conditions. These results demonstrated the potential to enhance sweetpotato productivity through stacking genes that are involved in ion compartmentalization and translation initiation.

DEAD-box RNA helicase Belle posttranscriptionally promotes gene expression in an ATPase activity-dependent manner

Drosophila Belle (human ortholog DDX3) is a conserved DEAD-box RNA helicase implicated in regulating gene expression. However, the molecular mechanisms by which Belle/DDX3 regulates gene expression are poorly understood. Here we performed systematic mutational analysis to determine the contributions of conserved motifs within Belle to its in vivo function. We found that Belle RNA-binding and RNA-unwinding activities and intrinsically disordered regions (IDRs) are required for Belle in vivo function. Expression of Belle ATPase mutants that cannot bind, hydrolyze, or release ATP resulted in dominant toxic phenotypes. Mechanistically, we discovered that Belle up-regulates reporter protein level when tethered to reporter mRNA, without corresponding changes at the mRNA level, indicating that Belle promotes translation of mRNA that it binds. Belle ATPase activity and amino-terminal IDR were required for this translational promotion activity. We also found that ectopic ovary expression of dominant Belle ATPase mutants decreases levels of cyclin proteins, including Cyclin B, without corresponding changes in their mRNA levels. Finally, we found that Belle binds endogenous cyclin B mRNA. We propose that Belle promotes translation of specific target mRNAs, including cyclin B mRNA, in an ATPase activity-dependent manner.

Comparative analysis of microRNA expression profiles of adult Schistosoma japonicum isolated from water buffalo and yellow cattle

BACKGROUND

Yellow cattle and water buffalo are important natural reservoir hosts and the main transmission sources of Schistosoma japonicum in endemic areas of China. The worms from the two hosts have marked differences in general worm morphology and ultrastructure, gene transcription and protein expression profiles.

RESULTS

To investigate microRNAs (miRNAs) involved in the regulation of schistosome development and survival, we compared miRNA expression profiles of adult schistosomes derived from yellow cattle and water buffalo by using high-throughput sequencing with Illumina Hiseq Xten. Schistosoma japonicum from water buffalo and yellow cattle yielded 63.78 million and 63.21 million reads, respectively, of which nearly 50% and 49% could be mapped to selected miRNAs in miRbase. A total of 206 miRNAs were identified, namely 79 previously annotated miRNAs of S. japonicum and 127 miRNAs that matched with the S. japonicum genome and were highly similar to the annotated miRNAs from other organisms. Among the 79 miRNAs, five (sja-miR-124-3p, sja-miR-219-5p, sja-miR-2e-3p, sja-miR-7-3p and sja-miR-3490) were significantly upregulated in the schistosomes from water buffalo compared with those from yellow cattle. A total of 268 potential target genes were predicted for these five differentially expressed miRNAs. Eleven differentially expressed targets were confirmed by qRT-PCR among 15 tested targets, one of which was further validated through dual-luciferase reporter assay. Among the 127 'possible' S. japonicum miRNAs, ten were significantly differentially expressed in the schistosomes from these two hosts.

CONCLUSIONS

These results highlight the important roles of miRNAs in regulating the development and survival of schistosomes in water buffalo and yellow cattle and facilitate understanding of the miRNA regulatory mechanisms in schistosomes derived from different susceptible hosts.

Auxin Modulated Initiation of Lateral Roots Is Linked to Pericycle Cell Length in Maize

Auxin is essential for the regulation of root system architecture by controlling primary root elongation and lateral root (LR) formation. Exogenous auxin has been reported to inhibit primary root elongation and promote the formation of LRs. In this study, LR formation in the Zea mays primary root was quantitatively evaluated after exogenous auxin treatment by comparing the effects of auxin on two selected zones elongated either before or after auxin application. We determined two main variables in both zones: the LR density per unit of root length (LRD), and the mean phloem pericycle cell length. The total number of phloem pericycle cells (PPCs) per unit of root length was then calculated. Considering that each LR primordium is initiated from four founder cells (FCs), the percentage of PPCs (%PPC) that behave as FCs in a specific root zone was estimated by dividing the number of pericycle cells by four times the LRD. This index was utilized to describe LR initiation. Root zones elongated in the presence of a synthetic auxin (1-naphthalene acetic acid, NAA) at low concentrations (0.01 μM) showed reduced cell length and increased LRD. However, a high concentration of NAA (0.1 μM) strongly reduced both cell length and LRD. In contrast, both low and high levels of NAA stimulated LRD in zones elongated before auxin application. Analysis of the percentage of FCs in the phloem pericycle in zones elongated in the presence or absence of NAA showed that low concentrations of NAA increased the %PFC, indicating that LR initiation is promoted at new sites; however, high concentrations of NAA elicited a considerable reduction in this variable in zones developed in the presence of auxin. As these zones are composed of short pericycle cells, we propose that short pericycle cells are incapable to participate in LR primordium initiation and that auxin modulated initiation of LRs is linked to pericycle cell length.

The MATH-BTB Protein TaMAB2 Accumulates in Ubiquitin-Containing Foci and Interacts With the Translation Initiation Machinery in Arabidopsis

MATH-BTB proteins are known to act as substrate-specific adaptors of CUL3-based E3 ligases in the ubiquitin proteasome pathway. Their BTB domain binds to CUL3 scaffold proteins and the less conserved MATH domain targets a highly diverse collection of substrate proteins to promote their ubiquitination and subsequent degradation. In plants, a significant expansion of the MATH-BTB family occurred in the grasses. Here, we report analysis of TaMAB2, a MATH-BTB protein transiently expressed at the onset of embryogenesis in wheat. Due to difficulties in studying its role in zygotes and early embryos, we have overexpressed TaMAB2 in Arabidopsis to generate gain-of-function mutants and to elucidate interaction partners and substrates. Overexpression plants showed severe growth defects as well as disorganization of microtubule bundles indicating that TaMAB2 interacts with substrates in Arabidopsis. In tobacco BY-2 cells, TaMAB2 showed a microtubule and ubiquitin-associated cytoplasmic localization pattern in form of foci. Its direct interaction with CUL3 suggests functions in targeting specific substrates for ubiquitin-dependent degradation. Although direct interactions with tubulin could not be confimed, tandem affinity purification of TaMAB2 interactors point towards cytoskeletal proteins including tubulin and actin as well as the translation initiation machinery. The idenification of various subunits of eucaryotic translation initiation factors eIF3 and eIF4 as TaMAB2 interactors indicate regulation of translation initiation as a major function during onset of embryogenesis in plants.

Emerging Roles of the Nuclear Cap-Binding Complex in Abiotic Stress Responses

Plant nuclear CBC consisted of two subunits (CBP20 and CBP80) is involved in both conserved processes related to RNA metabolism and simultaneously in extremely dynamic plant stress response.

Comparative expression profiling of AtRAD5B and AtNDL1: Hints towards a role in G protein mediated signaling

Arabidopsis AtRAD5B encodes for a putative helicase of the class SWItch/Sucrose Non-Fermentable (SWI/SNF) ATPases. We identified AtRAD5B as an interactor of N-MYC DOWNREGULATED-LIKE1 (AtNDL1) in a yeast two-hybrid screen. AtNDL1 is a G protein signaling component which regulates auxin transport and gradients together with GTP binding protein beta 1 (AGB1). Auxin gradients are known to recruit SWI/SNF remodeling complexes to the chromatin and regulate expression of genes involved in flower and leaf formation. In current study, a comparative spatial and temporal co-expression/localization analysis of AtNDL1, AGB1 with AtRAD5B was carried out in order to explore the possibility of their coexistence in a common signaling network. Translational fusion (GUS) of AtNDL1 and AtRAD5B in seedlings and reproductive organs revealed that both shared similar expression patterns with the highest expression observed in male reproductive organs. Moreover, they shared similar domains of localization in roots, suggesting their potential functioning together in reproductive and root development processes. This study predicts the existence of a signaling network involving AtNDL1, AGB1 with AtRAD5B.

In Planta Determination of the mRNA-Binding Proteome of Arabidopsis Etiolated Seedlings

RNA binding proteins (RBPs) control the fate and expression of a transcriptome. Despite this fundamental importance, our understanding of plant RBPs is rudimentary, being mainly derived via bioinformatic extrapolation from other kingdoms. Here, we adapted the mRNA-protein interactome capture method to investigate the RNA binding proteome in planta. From Arabidopsis thaliana etiolated seedlings, we captured more than 700 proteins, including 300 with high confidence that we have defined as the At-RBP set. Approximately 75% of these At-RBPs are bioinformatically linked with RNA biology, containing a diversity of canonical RNA binding domains (RBDs). As no prior experimental RNA binding evidence exists for the majority of these proteins, their capture now authenticates them as RBPs. Moreover, we identified protein families harboring emerging and potentially novel RBDs, including WHIRLY, LIM, ALBA, DUF1296, and YTH domain-containing proteins, the latter being homologous to animal RNA methylation readers. Other At-RBP set proteins include major signaling proteins, cytoskeleton-associated proteins, membrane transporters, and enzymes, suggesting the scope and function of RNA-protein interactions within a plant cell is much broader than previously appreciated. Therefore, our foundation data set has provided an unbiased insight into the RNA binding proteome of plants, on which future investigations into plant RBPs can be based.

eIF4A RNA Helicase Associates with Cyclin-Dependent Protein Kinase A in Proliferating Cells and Is Modulated by Phosphorylation

Eukaryotic initiation factor 4A (eIF4A) is a highly conserved RNA-stimulated ATPase and helicase involved in the initiation of messenger RNA translation. Previously, we found that eIF4A interacts with cyclin-dependent kinase A (CDKA), the plant ortholog of mammalian CDK1. Here, we show that this interaction occurs only in proliferating cells where the two proteins coassociate with 5'-cap-binding protein complexes, eIF4F or the plant-specific eIFiso4F. CDKA phosphorylates eIF4A on a conserved threonine residue (threonine-164) within the RNA-binding motif 1b TPGR. In vivo, a phospho-null (APGR) variant of the Arabidopsis (Arabidopsis thaliana) eIF4A1 protein retains the ability to functionally complement a mutant (eif4a1) plant line lacking eIF4A1, whereas a phosphomimetic (EPGR) variant fails to complement. The phospho-null variant (APGR) rescues the slow growth rate of roots and rosettes, together with the ovule-abortion and late-flowering phenotypes. In vitro, wild-type recombinant eIF4A1 and its phospho-null variant both support translation in cell-free wheat germ extracts dependent upon eIF4A, but the phosphomimetic variant does not support translation and also was deficient in ATP hydrolysis and helicase activity. These observations suggest a mechanism whereby CDK phosphorylation has the potential to down-regulate eIF4A activity and thereby affect translation.