R-loop 在植物中的研究方法

No abstract available

The impact of R-loops on the physiology and pathology of chromosomes has been demonstrated extensively by chromatin biology research. The progress in this field has been driven by technological advancement of R-loop mapping methods that largely relied on a single approach, DNA-RNA immunoprecipitation (DRIP). Most of the DRIP protocols use the experimental design that was developed by a few laboratories, without paying attention to the potential caveats that might affect the outcome of RNA-DNA hybrid mapping. To assess the accuracy and utility of this technology, we pursued an analytical approach to estimate inherent biases and errors in the DRIP protocol. By performing DRIP-sequencing, qPCR, and receiver operator characteristic (ROC) analysis, we tested the effect of formaldehyde fixation, cell lysis temperature, mode of genome fragmentation, and removal of free RNA on the efficacy of RNA-DNA hybrid detection and implemented workflows that were able to distinguish complex and weak DRIP signals in a noisy background with high confidence. We also show that some of the workflows perform poorly and generate random answers. Furthermore, we found that the most commonly used genome fragmentation method (restriction enzyme digestion) led to the overrepresentation of lengthy DRIP fragments over coding ORFs, and this bias was enhanced at the first exons. Biased genome sampling severely compromised mapping resolution and prevented the assignment of precise biological function to a significant fraction of R-loops. The revised workflow presented herein is established and optimized using objective ROC analyses and provides reproducible and highly specific RNA-DNA hybrid detection.

No abstract available

No abstract available

In eukaryotes, three-dimensional genome organization is critical for transcriptional regulation of gene expression. Long noncoding RNAs (lncRNAs) can modulate chromatin conformation of spatially related genomic locations within the nucleus. Here, we show that the lncRNA APOLO (AUXIN-REGULATED PROMOTER LOOP) recognizes multiple distant independent loci in the Arabidopsis thaliana genome. We found that APOLO targets are not spatially associated in the nucleus and that APOLO recognizes its targets by short sequence complementarity and the formation of DNA-RNA duplexes (R-loops). The invasion of APOLO to the target DNA decoys the plant Polycomb Repressive Complex 1 component LHP1, modulating local chromatin 3D conformation. APOLO lncRNA coordinates the expression of distal unrelated auxin-responsive genes during lateral root development in Arabidopsis. Hence, R-loop formation and chromatin protein decoy mediate trans action of lncRNAs on distant loci. VIDEO ABSTRACT.

No abstract available

The study of newly formed centromere with stable transmission ability can provide theoretical guidance for the construction of artificial chromosomes. More neocentromeres are needed to study the mechanisms of their formation. In this study, a minichromosome 7RLmini was derived from the progeny of wheat-rye 7R monosomic addition line. The minichromosome 7RLmini contained subtelomeric tandem repeats pSc119.2 and rye-specific pSc200, and it came from the distal region of the long arm of 7R chromosome. A neocentromere was formed in this minichromosome, and it did not contain centromeric repetitive sequences CCS1 and pAWRC.1. CENH3 ChIP-seq and ssDRIP-seq data confirmed that a 2.4 Mb segment from the rye 7R chromosome was involved in the neocentromere formation and enrichment of R-loops in this region. Within the 2.4 Mb segment, the GC content was higher that of AT, and a major binding position of CENH3 nucleosomes was identified on a 6 kb unknown LTR retrotransposon TE00002448. This unknown LTR retrotransposon was rye-specific and distributed through all the arms of rye chromosomes. The minichromosome exhibited stable generational transmission. A minichromosome from rye 7R with neocentromere was obtained in this study and the neocentromere was formed at the position far away from its native equivalent. This minichromosome provides additional material for the research on the mechanism of neocentromere formation. We theorize that R-loops and transposable element might be involved in the positioning of CENH3 nucleosomes in a functional neocentromere.

DNA binding proteins carry out important and diverse functions in the cell, including gene regulation, but identifying these proteins is technically challenging. In the present study, we developed a technique to capture DNA-associated proteins called reverse chromatin immunoprecipitation (R-ChIP). This technology uses a set of specific DNA probes labeled with biotin to isolate chromatin, and the DNA-associated proteins are then identified using mass spectrometry. Using R-ChIP, we identified 439 proteins that potentially bind to the promoter of the Arabidopsis thaliana gene AtCAT3 (AT1G20620). According to functional annotation, we randomly selected 5 transcription factors from these candidates, including bZIP1664, TEM1, bHLH106, BTF3, and HAT1, to verify whether they in fact bind to the AtCAT3 promoter. The binding of these 5 transcription factors was confirmed using chromatin immunoprecipitation quantitative real-time PCR and electrophoretic mobility shift assays. In addition, we improved the R-ChIP method using plants in which the DNA of interest had been transiently introduced, which does not require the T-DNA integration, and showed that this substantially improved the protein capture efficiency. These results together demonstrate that R-ChIP has a wide application to characterize chromatin composition and isolate upstream regulators of a specific gene. Xuejing Wen et al. present a new method, Reverse Chromatin Immunoprecipitation (R-ChIP), for analyzing DNA-protein interactions specifically in plant cells. They apply R-ChIP to identify proteins binding to the promoter of the Arabidopsis gene AtCAT3 and optimize the method using transient transformation of the target promoter to increase efficiency.

No abstract available

A DNA-RNA hybrid recognition sensor protein profiles native R loop structures in the genome. An R loop is a unique triple-stranded structure that participates in multiple key biological processes and is relevant to human diseases. Accurate and comprehensive R loop profiling is a prerequisite for R loops studies. However, current R loop mapping methods generate large discrepancies, therefore an independent method is in urgent need. Here, we establish an independent R loop CUT&Tag (Tn5-based cleavage under targets and tagmentation) method by combining CUT&Tag and GST-His6-2×HBD (glutathione S-transferase–hexahistidine–2× hybrid-binding domain), an artificial DNA-RNA hybrid sensor that specifically recognizes the DNA-RNA hybrids. We demonstrate that the R loop CUT&Tag is sensitive, reproducible, and convenient for native R loop mapping with high resolution, and find that the capture strategies, instead of the specificity of sensors, largely contribute to the disparities among different methods. Together, we provide an independent strategy for genomic profiling of native R loops and help resolve discrepancies among multiple R loop mapping methods.

Programmed constitutive heterochromatin silencing is essential for eukaryotic genome regulation, yet the initial step of this process is ambiguous. A large proportion of R-loops (RNA:DNA hybrids) had been unexpectedly identified within Arabidopsis pericentromeric heterochromatin with unknown functions. Through a genome-wide R-loop profiling screen, we find that DDM1 (decrease in DNA methylation 1) is the primary restrictor of pericentromeric R-loops via its RNA:DNA helicase activity. Low levels of pericentromeric R-loops resolved by DDM1 cotranscriptionally can facilitate constitutive heterochromatin silencing. Furthermore, we demonstrate that DDM1 physically excludes histone H2A variant H2A.Z and promotes H2A.W deposition for faithful heterochromatin initiation soon after R-loop clearance. The dual functions of DDM1 in R-loop resolution and H2A.Z eviction are essential for sperm nuclei structure maintenance in mature pollen. Our work unravels the cotranscriptional R-loop resolution coupled with accurate H2A variants deposition is the primary step of constitutive heterochromatin silencing in Arabidopsis, which might be conserved across eukaryotes.

Plant mitochondrial genomes undergo frequent homologous recombination (HR). Ectopic HR activity is inhibited by the HR surveillance pathway, but the underlying regulatory mechanism is unclear. Here, we show that the mitochondrial RNase H1 AtRNH1B impairs the formation of RNA:DNA hybrids (R-loops) and participates in the HR surveillance pathway in Arabidopsis thaliana. AtRNH1B suppresses ectopic HR at intermediate-sized repeats (IRs) and thus maintains mitochondrial DNA (mtDNA) replication. The RNase H1 AtRNH1C is restricted to the chloroplast; however, when cells lack AtRNH1B, transport of chloroplast AtRNH1C into the mitochondria secures HR surveillance, thus ensuring the integrity of the mitochondrial genome and allowing embryogenesis to proceed. HR surveillance is further regulated by the single-stranded DNA-binding protein ORGANELLAR SINGLE-STRANDED DNA BINDING PROTEIN1 (OSB1), which decreases the formation of R-loops. This study uncovers a facultative dual targeting mechanism between organelles and sheds light on the roles of RNase H1 in organellar genome maintenance and embryogenesis.

R-loops, three-stranded nucleic acid structures consisting of a DNA: RNA hybrid and displaced single-stranded DNA, play critical roles in gene expression and genome stability. How R-loop homeostasis is integrated into chloroplast gene expression remains largely unknown. We found an unexpected function of FtsHi1, an inner envelope membrane-bound AAA-ATPase in chloroplast R-loop homeostasis of Arabidopsis thaliana. Previously, this protein was shown to function as a component of the import motor complex for nuclear-encoded chloroplast proteins. However, this study provides evidence that FtsHi1 is an ATP-dependent helicase that efficiently unwinds both DNA-DNA and DNA-RNA duplexes, thereby preventing R-loop accumulation. Over-accumulation of R-loops could impair chloroplast transcription but not necessarily genome integrity. The dual function of FtsHi1 in both protein import and chloroplast gene expression may be important to coordinate the biogenesis of nuclear- and chloroplast-encoded subunits of multi-protein photosynthetic complexes. This study suggests a mechanical link between protein import and R-loop homeostasis in chloroplasts of higher plants.

No abstract available

No abstract available

No abstract available

Brassinosteroids, which control plant growth and development, are sensed by the membrane receptor kinase BRASSINOSTEROID INSENSITIVE 1 (BRI1). Brassinosteroid binding to the BRI1 leucine-rich repeat (LRR) domain induces heteromerisation with a SOMATIC EMBRYOGENESIS RECEPTOR KINASE (SERK)-family co-receptor. This process allows the cytoplasmic kinase domains of BRI1 and SERK to interact, trans-phosphorylate and activate each other. Here we report crystal structures of the BRI1 kinase domain in its activated form and in complex with nucleotides. BRI1 has structural features reminiscent of both serine/threonine and tyrosine kinases, providing insight into the evolution of dual-specificity kinases in plants. Phosphorylation of Thr1039, Ser1042 and Ser1044 causes formation of a catalytically competent activation loop. Mapping previously identified serine/threonine and tyrosine phosphorylation sites onto the structure, we analyse their contribution to brassinosteroid signaling. The location of known genetic missense alleles provide detailed insight into the BRI1 kinase mechanism, while our analyses are inconsistent with a previously reported guanylate cyclase activity. We identify a protein interaction surface on the C-terminal lobe of the kinase and demonstrate that the isolated BRI1, SERK2 and SERK3 cytoplasmic segments form homodimers in solution and have a weak tendency to heteromerise. We propose a model in which heterodimerisation of the BRI1 and SERK ectodomains brings their cytoplasmic kinase domains in a catalytically competent arrangement, an interaction that can be modulated by the BRI1 inhibitor protein BKI1.

No abstract available

No abstract available

No abstract available

Rc is a domestication-related gene required for red pericarp in rice (Oryza sativa). The red grain color is ubiquitous among the wild ancestors of O. sativa, in which it is closely associated with seed shattering and dormancy. Rc encodes a basic helix-loop-helix (bHLH) protein that was fine-mapped to an 18.5-kb region on rice chromosome 7 using a cross between Oryza rufipogon (red pericarp) and O. sativa cv Jefferson (white pericarp). Sequencing of the alleles from both mapping parents as well as from two independent genetic stocks of Rc revealed that the dominant red allele differed from the recessive white allele by a 14-bp deletion within exon 6 that knocked out the bHLH domain of the protein. A premature stop codon was identified in the second mutant stock that had a light red pericarp. RT-PCR experiments confirmed that the Rc gene was expressed in both red- and white-grained rice but that a shortened transcript was present in white varieties. Phylogenetic analysis, supported by comparative mapping in rice and maize (Zea mays), showed that Rc, a positive regulator of proanthocyanidin, is orthologous with INTENSIFIER1, a negative regulator of anthocyanin production in maize, and is not in the same clade as rice bHLH anthocyanin regulators.

Abstract Identifying potential molecular tags for drought tolerance is essential for achieving higher crop productivity under drought stress. We employed an integrated genomics-assisted breeding and functional genomics strategy involving association mapping, fine mapping, map-based cloning, molecular haplotyping and transcript profiling in the introgression lines (ILs)- and near isogenic lines (NILs)-based association panel and mapping population of chickpea (Cicer arietinum). This combinatorial approach delineated a bHLH (basic helix–loop–helix) transcription factor, CabHLH10 (Cicer arietinum bHLH10) underlying a major QTL, along with its derived natural alleles/haplotypes governing yield traits under drought stress in chickpea. CabHLH10 binds to a cis-regulatory G-box promoter element to modulate the expression of RD22 (responsive to desiccation 22), a drought/abscisic acid (ABA)-responsive gene (via a trans-expression QTL), and two strong yield-enhancement photosynthetic efficiency (PE) genes. This, in turn, upregulates other downstream drought-responsive and ABA signaling genes, as well as yield-enhancing PE genes, thus increasing plant adaptation to drought with reduced yield penalty. We showed that a superior allele of CabHLH10 introgressed into the NILs improved root and shoot biomass and PE, thereby enhancing yield and productivity during drought without compromising agronomic performance. Furthermore, overexpression of CabHLH10 in chickpea and Arabidopsis (Arabidopsis thaliana) conferred enhanced drought tolerance by improving root and shoot agro-morphological traits. These findings facilitate translational genomics for crop improvement and the development of genetically tailored, climate-resilient, high-yielding chickpea cultivars.

No abstract available

No abstract available

No abstract available

No abstract available

R-loops are three-stranded nucleic acid structures formed from the hybridization of RNA and DNA during nascent transcription. In 2012, Ginno et al. introduced the first R-loop mapping method, DNA:RNA immunoprecipitation (DRIP) sequencing. Since that time, dozens of studies have implemented R-loop mapping and new high-resolution techniques have been developed. The resulting datasets have tremendous potential to reveal the causes and consequences of R-loops genome-wide. However, poor quality and variability between mapping approaches pose serious barriers to the meta-analysis of these data. In our recent work, we reprocessed 693 R-loop mapping samples, devising new quality methods, defining a set of high-confidence mapping samples, and then deriving R-loop regions, consensus sites of R-loop formation. This analysis yielded the largest R-loop data resource to date along with novel computational approaches for R-loop mapping analysis. Now, we introduce RLBase, an innovative web server which builds upon those data and software by providing users with the capability to (1) explore hundreds of public R-loop mapping datasets, (2) explore consensus R-loop regions, (3) analyze user-supplied datasets to generate an HTML quality report, and (4) download all the processed data for the 693 samples we previously reprocessed and standardized. In addition to RLBase, we also describe the other software which, along with RLBase, provides a computational framework for R-loop bioinformatics. RLBase, and the rest of these software (termed “RLSuite”), are provided freely under an MIT license and made publicly available: https://gccri.bishop-lab.uthscsa.edu/rlsuite/. RLBase is directly accessible via the following URL: https://gccri.bishop-lab.uthscsa.edu/rlbase/.

R-loops are involved in transcriptional regulation, DNA and histone post-translational modifications, genome replication and genome stability. To what extent R-loop abundance and genome-wide localization is actively regulated during metazoan embryogenesis is unknown. Drosophila embryogenesis provides a powerful system to address these questions due to its well-characterized developmental program, the sudden onset of zygotic transcription and available genome-wide ChIP and transcription data sets. Here, we measure the overall abundance and genome localization of R-loops in early and late-stage embryos relative to Drosophila cultured cells. We demonstrate that absolute R-loop levels change during embryogenesis and that resolution of R-loops is critical for embryonic development. R-loop mapping by strand-specific DRIP-seq reveals that R-loop localization is plastic across development, both in the genes which form R-loops and where they localize relative to gene bodies. Importantly, these changes are not driven by changes in the transcriptional program. Negative GC skew and absolute changes in AT skew are associated with R-loop formation in Drosophila. Furthermore, we demonstrate that while some chromatin binding proteins and histone modification such as H3K27me3 are associated with R-loops throughout development, other chromatin factors associated with R-loops in a developmental specific manner. Our findings highlight the importance and developmental plasticity of R-loops during Drosophila embryogenesis.

No abstract available

No abstract available

No abstract available

No abstract available

R-loops are stable RNA-DNA hybrids that have been implicated in transcription initiation and termination, as well as in telomere homeostasis, chromatin formation, and genome replication and instability. RNA Polymerase (Pol) II transcription in the protozoan parasite Trypanosoma brucei is highly unusual: virtually all genes are co-transcribed from multigene transcription units, with mRNAs generated by linked trans-splicing and polyadenylation, and transcription initiation sites display no conserved promoter motifs. Here, we describe the genome-wide distribution of R-loops in wild type mammal-infective T. brucei and in mutants lacking RNase H1, revealing both conserved and diverged functions. Conserved localisation was found at centromeres, rRNA genes and retrotransposon-associated genes. RNA Pol II transcription initiation sites also displayed R-loops, suggesting a broadly conserved role despite the lack of promoter conservation or transcription initiation regulation. However, the most abundant sites of R-loop enrichment were within the intergenic regions of the multigene transcription units, where the hybrids coincide with sites of polyadenylation and nucleosome-depletion. Thus, instead of functioning in transcription termination, most T. brucei R-loops act in a novel role, promoting RNA Pol II movement or mRNA processing. Finally, we show there is little evidence for correlation between R-loop localisation and mapped sites of DNA replication initiation.

No abstract available

R-loops affect transcription and genome stability. Dysregulation of R-loops is related to human diseases. Genome-wide R-loop mapping typically uses the S9.6 antibody or inactive ribonuclease H, both requiring a large number of cells with varying results observed depending on the approach applied. Here, we present strand-specific kethoxal-assisted single-stranded DNA (ssDNA) sequencing (spKAS-seq) to map R-loops by taking advantage of the presence of a ssDNA in the triplex structure. We show that spKAS-seq detects R-loops and their dynamics at coding sequences, enhancers, and other intergenic regions with as few as 50,000 cells. A joint analysis of R-loops and chromatin-bound RNA binding proteins (RBPs) suggested that R-loops can be RBP binding hotspots on the chromatin.

Age‐related loss of cellular function and increased cell death are characteristic hallmarks of aging. While defects in gene expression and RNA metabolism have been linked with age‐associated human neuropathies, it is not clear how the changes that occur in aging neurons contribute to loss of gene expression homeostasis. R‐loops are RNA–DNA hybrids that typically form co‐transcriptionally via annealing of the nascent RNA to the template DNA strand, displacing the non‐template DNA strand. Dysregulation of R‐loop homeostasis has been associated with both transcriptional impairment and genome instability. Importantly, a growing body of evidence links R‐loop accumulation with cellular dysfunction, increased cell death, and chronic disease onset. Here, we characterized the R‐loop landscape in aging Drosophila melanogaster photoreceptor neurons and showed that bulk R‐loop levels increased with age. Further, genome‐wide mapping of R‐loops revealed that transcribed genes accumulated R‐loops over gene bodies during aging, which correlated with decreased expression of long and highly expressed genes. Importantly, while photoreceptor‐specific down‐regulation of Top3β, a DNA/RNA topoisomerase associated with R‐loop resolution, lead to decreased visual function, over‐expression of Top3β or nuclear‐localized RNase H1, which resolves R‐loops, enhanced positive light response during aging. Together, our studies highlight the functional link between dysregulation of R‐loop homeostasis, gene expression, and visual function during aging.

Abstract R-loops play versatile roles in many physiological and pathological processes, and are of great interest to scientists in multiple fields. However, controversy about their genomic localization and incomplete understanding of their regulatory network raise great challenges for R-loop research. Here, we present R-loopBase (https://rloopbase.nju.edu.cn) to tackle these pressing issues by systematic integration of genomics and literature data. First, based on 107 high-quality genome-wide R-loop mapping datasets generated by 11 different technologies, we present a reference set of human R-loop zones for high-confidence R-loop localization, and spot conservative genomic features associated with R-loop formation. Second, through literature mining and multi-omics analyses, we curate the most comprehensive list of R-loop regulatory proteins and their targeted R-loops in multiple species to date. These efforts help reveal a global regulatory network of R-loop dynamics and its potential links to the development of cancers and neurological diseases. Finally, we integrate billions of functional genomic annotations, and develop interactive interfaces to search, visualize, download and analyze R-loops and R-loop regulators in a well-annotated genomic context. R-loopBase allows all users, including those with little bioinformatics background to utilize these data for their own research. We anticipate R-loopBase will become a one-stop resource for the R-loop community.

R-loops are a prevalent class of non-B DNA structures that have been associated with both positive and negative cellular outcomes. DNA:RNA immunoprecipitation (DRIP) approaches based on the anti-DNA:RNA hybrid S9.6 antibody revealed that R-loops form dynamically over conserved genic hotspots. We have developed an orthogonal approach that queries R-loops via the presence of long stretches of single-stranded DNA on their looped-out strand. Non-denaturing sodium bisulfite treatment catalyzes the conversion of unpaired cytosines to uracils, creating permanent genetic tags for the position of an R-loop. Long read, single-molecule PacBio sequencing allows the identification of R-loop 'footprints' at near nucleotide resolution in a strand-specific manner on long single DNA molecules and at ultra-deep coverage. Single-molecule R-loop footprinting (SMRF-seq) revealed a strong agreement between S9.6- and bisulfite-based R-loop mapping and confirmed that R-loops form over genic hotspots, including gene bodies and terminal gene regions. Based on the largest single-molecule R-loop dataset to date, we show that individual R-loops form non-randomly, defining discrete sets of overlapping molecular clusters that pile-up through larger R-loop zones. R-loops most often map to intronic regions and their individual start and stop positions do not match with intron-exon boundaries, reinforcing the model that they form co-transcriptionally from unspliced transcripts. SMRF-seq further established that R-loop distribution patterns are not simply driven by intrinsic DNA sequence features but most likely also reflect DNA topological constraints. Overall, DRIP-based and SMRF-based approaches independently provide a complementary and congruent view of R-loop distribution, consolidating our understanding of the principles underlying R-loop formation.

Genome-wide analysis of R-loop alterations in U2OS cells deficient of DDX5, XRN2, and PRMT5 identify >50,650 DRIP-seq peaks spanning ∼4.5% of the genomic sequence. R-loops near TSS generated intergenic antisense transcription. DDX5, XRN2, and PRMT5 have been shown to resolve DNA/RNA hybrids (R-loops) at RNA polymerase II transcription termination sites at few genomic loci. Herein, we perform genome-wide R-loop mapping using classical DNA/RNA immunoprecipitation and high-throughput sequencing (DRIP-seq) of loci regulated by DDX5, XRN2, and PRMT5. We observed hundreds to thousands of R-loop gains and losses at transcribed loci in DDX5-, XRN2-, and PRMT5-deficient U2OS cells. R-loop gains were characteristic of highly transcribed genes located at gene-rich regions, whereas R-loop losses were observed in low-density gene areas. DDX5, XRN2, and PRMT5 shared many R-loop gain loci at transcription termination sites, consistent with their coordinated role in RNA polymerase II transcription termination. DDX5-depleted cells had unique R-loop gain peaks near the transcription start site that did not overlap with those of siXRN2 and siPRMT5 cells, suggesting a role for DDX5 in transcription initiation independent of XRN2 and PRMT5. Moreover, we observed that the accumulated R-loops at certain loci in siDDX5, siXRN2, and siPRMT5 cells near the transcription start site of genes led to antisense intergenic transcription. Our findings define unique and shared roles of DDX5, XRN2, and PRMT5 in DNA/RNA hybrid regulation.

RNA-Seq is a whole-transcriptome analysis method used to research biological mechanisms and functions but its use in large-scale experiments is limited by its high cost and labour requirements. In this study, we have established a high-throughput and cost-effective RNA-Seq library preparation method that does not require mRNA enrichment. The method adds unique index sequences to samples during reverse transcription (RT) that is conducted at a higher temperature (≥62 °C) to suppress RT of A-rich sequences in rRNA, and then pools all samples into a single tube. Both single-read and paired-end sequencing of libraries is enabled. We found that the pooled RT products contained large amounts of RNA, mainly rRNA, causing over-estimations of the quantity of DNA and unstable tagmentation results. Degradation of RNA before tagmentation was found to be necessary for the stable preparation of libraries. We named this protocol low-cost and easy RNA-Seq (Lasy-Seq) and used it to investigate temperature responses in Arabidopsis thaliana. We analysed how sub-ambient temperatures (10–30 °C) affected the plant transcriptomes using time-courses of RNA-Seq from plants grown in randomly fluctuating temperature conditions. Our results suggest that there are diverse mechanisms behind plant temperature responses at different time scales.

The availability of data produced from various sequencing platforms offer the possibility to answer complex questions in plant research. However, drawbacks can arise when there are gaps in the information generated, and complementary platforms are essential to obtain more comprehensive data sets relating to specific biological process, such as responses to environmental perturbations in plant systems. The investigation of transcriptional regulation raises different challenges, particularly in associating differentially expressed transcription factors with their downstream responsive genes. In this paper, we discuss the integration of transcriptional factor studies through RNA sequencing (RNA-seq) and Chromatin Immunoprecipitation sequencing (ChIP-seq). We show how the data from ChIP-seq can strengthen information generated from RNA-seq in elucidating gene regulatory mechanisms. In particular, we discuss how integration of ChIP-seq and RNA-seq data can help to unravel transcriptional regulatory networks. This review discusses recent advances in methods for studying transcriptional regulation using these two methods. It also provides guidelines for making choices in selecting specific protocols in RNA-seq pipelines for genome-wide analysis to achieve more detailed characterization of specific transcription regulatory pathways via ChIP-seq.

Background Cotton has tremendous economic value worldwide; however, its allopolyploid nature and time-consuming transformation methods have hampered the development of cotton functional genomics. The protoplast system has proven to be an important and versatile tool for functional genomics, tissue-specific marker gene identification, tracking developmental trajectories, and genome editing in plants. Nevertheless, the isolation of abundant viable protoplasts suitable for single-cell RNA sequencing (scRNA-seq) and genome editing remains a challenge in cotton. Results We established an efficient transient gene expression system using protoplasts isolated from cotton taproots. The system enables the isolation of large numbers of viable protoplasts and uses an optimized PEG-mediated transfection protocol. The highest yield (3.55 × 10^5/g) and viability (93.3%) of protoplasts were obtained from cotton roots grown in hydroponics for 72 h. The protoplasts isolated were suitable for scRNA-seq. The highest transfection efficiency (80%) was achieved when protoplasts were isolated as described above and transfected with 20 μg of plasmid for 20 min in a solution containing 200 mM Ca^2+. Our protoplast-based transient expression system is suitable for various applications, including validation the efficiency of CRISPR vectors, protein subcellular localization analysis, and protein–protein interaction studies. Conclusions The protoplast isolation and transfection protocol developed in this study is stable, versatile, and time-saving. It will accelerate functional genomics and molecular breeding in cotton.

Plant single-cell RNA-seq technology quantifies the abundance of plant transcripts at a single-cell resolution. Deciphering the transcriptomes of each plant cell, their regulation during plant cell development, and their response to environmental stresses will support the functional study of genes, the establishment of precise transcriptional programs, the prediction of more accurate gene regulatory networks, and, in the long term, the design of de novo gene pathways to enhance selected crop traits. In this review, we will discuss the opportunities, challenges, and problems, and share tentative solutions associated with the generation and analysis of plant single-cell transcriptomes. We will discuss the benefit and limitations of using plant protoplasts vs. nuclei to conduct single-cell RNA-seq experiments on various plant species and organs, the functional annotation of plant cell types based on their transcriptomic profile, the characterization of the dynamic regulation of the plant genes during cell development or in response to environmental stress, the need to characterize and integrate additional layers of -omics datasets to capture new molecular modalities at the single-cell level and reveal their causalities, the deposition and access to single-cell datasets, and the accessibility of this technology to plant scientists.

Recently, single-cell RNA sequencing (scRNA-seq) provides unprecedented power for accurately understanding gene expression regulatory mechanisms. However, scRNA-seq studies have limitations in plants, due to difficulty in protoplast isolation that requires enzymatic digestion of the cell walls from various plant tissues. Therefore, to overcome this problem, we developed a nuclei isolation approach that does not rely on Fluorescence Activated Cell Sorting (FACS). We validated the robustness of the FACS-free single-nucleus RNA sequencing (snRNA-seq) methodology in mature Arabidopsis plant tissue by comparing it to scRNA-seq results based on protoplasts extracted from the same batch of leaf materials. Sequencing results demonstrated the high quality of snRNA-seq data, as well as its utility in cell type classification and marker gene identification. This approach also showed several advantages, including the ability to use frozen samples, taking less suspension preparation time, and reducing biased cellular coverage and dissociation-induced transcriptional artifacts. Surprisingly, snRNA-seq detected two epidermal pavement cell clusters, while scRNA-seq only had one. Furthermore, we hypothesized that these two epidermal cells represent the top and lower epidermis based on differences in expression patterns of cluster-specific expressed genes. In summary, this study has advanced the application of snRNA-seq in Arabidopsis leaves and confirmed the advantages of snRNA-seq in plant research.

Nitrogen (N) is an important contributor in regulating plant growth and development as well as secondary metabolites synthesis, so as to promote the formation of tea quality and flavor. Theanine, polyphenols, and caffeine are important secondary metabolites in tea plant. In this study, the responses of Camellia sinensis roots to N deprivation and resupply were investigated by metabolome and RNA-seq analysis. N deficiency induced content increase for most amino acids (AAs) and reduction for the remaining AAs, polyphenols, and caffeine. After N recovery, the decreased AAs and polyphenols showed a varying degree of recovery in content, but caffeine did not. Meanwhile, theanine increased in content, but its related synthetic genes were down-regulated, probably due to coordination of the whole N starvation regulatory network. Flavonoids-related pathways were relatively active following N stress according to KEGG enrichment analysis. Gene co-expression analysis revealed TCS2, AMT1;1, TAT2, TS, and GOGAT as key genes, and TFs like MYB, bHLH, and NAC were also actively involved in N stress responses in C. sinensis roots. These findings facilitate the understanding of the molecular mechanism of N regulation in tea roots and provide genetic reference for improving N use efficiency in tea plant.

A high-throughput assay in transiently transformed tobacco leaves identifies enhancers, characterizes their functional elements and detects condition-specific enhancer activity. Genetic engineering of cis-regulatory elements in crop plants is a promising strategy to ensure food security. However, such engineering is currently hindered by our limited knowledge of plant cis-regulatory elements. Here, we adapted self-transcribing active regulatory region sequencing (STARR-seq)—a technology for the high-throughput identification of enhancers—for its use in transiently transformed tobacco (Nicotiana benthamiana) leaves. We demonstrate that the optimal placement in the reporter construct of enhancer sequences from a plant virus, pea (Pisum sativum) and wheat (Triticum aestivum), was just upstream of a minimal promoter and that none of these four known enhancers was active in the 3′ untranslated region of the reporter gene. The optimized assay sensitively identified small DNA regions containing each of the four enhancers, including two whose activity was stimulated by light. Furthermore, we coupled the assay to saturation mutagenesis to pinpoint functional regions within an enhancer, which we recombined to create synthetic enhancers. Our results describe an approach to define enhancer properties that can be performed in potentially any plant species or tissue transformable by Agrobacterium and that can use regulatory DNA derived from any plant genome.

High-throughput DNA sequencing (HTS) presents great opportunities for plant systematics, yet genomic complexity needs to be reduced for HTS to be effectively applied. We highlight Hyb-Seq as a promising approach, especially in light of the recent development of probes enriching 353 low-copy nuclear genes from any flowering plant taxon.

Summary Although the regulatory mechanisms of dark and light‐induced plant morphogenesis have been broadly investigated, the biological process in peanuts has not been systematically explored on single‐cell resolution. Herein, 10 cell clusters were characterized using scRNA‐seq‐identified marker genes, based on 13 409 and 11 296 single cells from 1‐week‐old peanut seedling leaves grown under dark and light conditions. 6104 genes and 50 transcription factors (TFs) displayed significant expression patterns in distinct cell clusters, which provided gene resources for profiling dark/light‐induced candidate genes. Further pseudo‐time trajectory and cell cycle evidence supported that dark repressed the cell division and perturbed normal cell cycle, especially the PORA abundances correlated with 11 TFs highly enriched in mesophyll to restrict the chlorophyllide synthesis. Additionally, light repressed the epidermis cell developmental trajectory extending by inhibiting the growth hormone pathway, and 21 TFs probably contributed to the different genes transcriptional dynamic. Eventually, peanut AHL17 was identified from the profile of differentially expressed TFs, which encoded protein located in the nucleus promoted leaf epidermal cell enlargement when ectopically overexpressed in Arabidopsis through the regulatory phytohormone pathway. Overall, our study presents the different gene atlases in peanut etiolated and green seedlings, providing novel biological insights to elucidate light‐induced leaf cell development at the single‐cell level.

Abstract Drought is a significant global environmental stress. Biostimulants offer a sustainable solution to enhance crop tolerance and mitigate productivity losses. This study assessed the impact of foliar application of ERANTHIS®, a biostimulant derived from the algae Ascophyllum nodosum and Laminaria digitata and yeast extracts, on tomato plants under mild water stress. Evaluations were conducted at 5 and 24 hours after the third treatment. Under optimal water conditions, the biostimulant showed a priming effect, with an early increase of stress markers and a timing‐specific modulation of ROS non enzymatic and enzymatic ROS scavenging activities. Under drought stress, the biostimulant later decreased stress markers, by aligning the majority of analyzed ROS scavengers closer to levels in well‐irrigated plants. Transcriptome analysis using RNA‐seq data revealed differentially expressed genes (DEGs) and multivariate data highlighted groups of co‐regulated genes (k‐means clustering). Genes involved in water channel activity, transcription regulator activity, and oxidoreductase activity were significantly modulated. Cluster analysis identified distinct gene clusters influenced by the biostimulant under optimal conditions, including early responses (cell wall modification, hormone signaling) and late responses (RNA modification, nutrient uptake process). Under water stress, early responses involved actin filament organization and MAPK signaling, while late responses were related to plasma membrane components and cell wall organization. This study, integrating biochemical and transcriptomic data, provides a comprehensive understanding of how a biostimulant primes plants under optimal conditions and mitigates water stress effects, offering valuable insights for sustainable agriculture.

Summary Branch angle (BA) is a critical morphological trait that significantly influences planting density, light interception and ultimately yield in plants. Despite its importance, the regulatory mechanism governing BA in rapeseed remains poorly understood. In this study, we generated 109 transcriptome data sets for 37 rapeseed accessions with divergent BA phenotypes. Relative to adaxial branch segments, abaxial segments accumulated higher levels of auxin and exhibited lower expression of six TCP1 homologues and one GA20ox3. A co‐expression network analysis identified two modules highly correlated with BA. The modules contained homologues to known BA control genes, such as FUL, YUCCA6, TCP1 and SGR3. Notably, a homoeologous exchange (HE), occurring at the telomeres of A09, was prevalent in large BA accessions, while an A02‐C02 HE was common in small BA accessions. In their corresponding regions, these HEs explained the formation of hub gene hotspots in the two modules. QTL‐seq analysis confirmed that the presence of a large A07‐C06 HE (~8.1 Mb) was also associated with a small BA phenotype, and BnaA07.WRKY40.b within it was predicted as candidate gene. Overexpressing BnaA07.WRKY40.b in rapeseed increased BA by up to 20°, while RNAi‐ and CRISPR‐mediated mutants (BnaA07.WRKY40.b and BnaC06.WRKY40.b) exhibited decreased BA by up to 11.4°. BnaA07.WRKY40.b was exclusively localized to the nucleus and exhibited strong expression correlations with many genes related to gravitropism and plant architecture. Taken together, our study highlights the influence of HEs on rapeseed plant architecture and confirms the role of WRKY40 homologues as novel regulators of BA.

The genome-wide expression profile of genes in different tissues/cell types and developmental stages is a vital component of many functional genomic studies. Transcriptome data obtained by RNA-sequencing (RNA-Seq) is often deposited in public databases that are made available via data portals. Data visualization is one of the first steps in assessment and hypothesis generation. However, these databases do not typically include visualization tools and establishing one is not trivial for users who are not computational experts. This, as well as the various formats in which data is commonly deposited, makes the processes of data access, sharing and utility more difficult. Our goal was to provide a simple and user-friendly repository that meets these needs for data-sets from major agricultural crops. AgriSeqDB (https://expression.latrobe.edu.au/agriseqdb) is a database for viewing, analysing and interpreting developmental and tissue/cell-specific transcriptome data from several species, including major agricultural crops such as wheat, rice, maize, barley and tomato. The disparate manner in which public transcriptome data is often warehoused and the challenge of visualizing raw data are both major hurdles to data reuse. The popular eFP browser does an excellent job of presenting transcriptome data in an easily interpretable view, but previous implementation has been mostly on a case-by-case basis. Here we present an integrated visualisation database of transcriptome data-sets from six species that did not previously have public-facing visualisations. We combine the eFP browser, for gene-by-gene investigation, with the Degust browser, which enables visualisation of all transcripts across multiple samples. The two visualisation interfaces launch from the same point, enabling users to easily switch between analysis modes. The tools allow users, even those without bioinformatics expertise, to mine into data-sets and understand the behaviour of transcripts of interest across samples and time. We have also incorporated an additional graphic download option to simplify incorporation into presentations or publications. Powered by eFP and Degust browsers, AgriSeqDB is a quick and easy-to-use platform for data analysis and visualization in five crops and Arabidopsis. Furthermore, it provides a tool that makes it easy for researchers to share their data-sets, promoting research collaborations and data-set reuse.

Summary Terminal drought is a major constraint to chickpea productivity. Two component traits responsible for reduction in yield under drought stress include reduction in seeds size and root length/root density. QTL‐seq approach, therefore, was used to identify candidate genomic regions for 100‐seed weight (100SDW) and total dry root weight to total plant dry weight ratio (RTR) under rainfed conditions. Genomewide SNP profiling of extreme phenotypic bulks from the ICC 4958 × ICC 1882 population identified two significant genomic regions, one on CaLG01 (1.08 Mb) and another on CaLG04 (2.7 Mb) linkage groups for 100SDW. Similarly, one significant genomic region on CaLG04 (1.10 Mb) was identified for RTR. Comprehensive analysis revealed four and five putative candidate genes associated with 100SDW and RTR, respectively. Subsequently, two genes (Ca_04364 and Ca_04607) for 100SDW and one gene (Ca_04586) for RTR were validated using CAPS/dCAPS markers. Identified candidate genomic regions and genes may be useful for molecular breeding for chickpea improvement.

No abstract available

The halophyte Mesembryanthemum crystallinum (common or crystalline ice plant) is a useful model for studying molecular mechanisms of salt tolerance. The morphology, physiology, metabolism, and gene expression of ice plant have been studied and large-scale analyses of gene expression profiling have drawn an outline of salt tolerance in ice plant. A rapid root growth to a sudden increase in salinity was observed in ice plant seedlings. Using a fluorescent dye to detect Na+, we found that ice plant roots respond to an increased flux of Na+ by either secreting or storing Na+ in specialized cells. High-throughput sequencing was used to identify small RNA profiles in 3-day-old seedlings treated with or without 200 mM NaCl. In total, 135 conserved miRNAs belonging to 21 families were found. The hairpin precursor of 19 conserved mcr-miRNAs and 12 novel mcr-miRNAs were identified. After 6 h of salt stress, the expression of most mcr-miRNAs showed decreased relative abundance, whereas the expression of their corresponding target genes showed increased mRNA relative abundance. The cognate target genes are involved in a broad range of biological processes: transcription factors that regulate growth and development, enzymes that catalyze miRNA biogenesis for the most conserved mcr-miRNA, and proteins that are involved in ion homeostasis and drought-stress responses for some novel mcr-miRNAs. Analyses of the functions of target genes revealed that cellular processes, including growth and development, metabolism, and ion transport activity are likely to be enhanced in roots under salt stress. The expression of eleven conserved miRNAs and two novel miRNAs were correlated reciprocally with predicted targets within hours after salt stress exposure. Several conserved miRNAs have been known to regulate root elongation, root apical meristem activity, and lateral root formation. Based upon the expression pattern of miRNA and target genes in combination with the observation of Na+ distribution, ice plant likely responds to increased salinity by using Na+ as an osmoticum for cell expansion and guard cell opening. Excessive Na+ could either be secreted through the root epidermis or stored in specialized leaf epidermal cells. These responses are regulated in part at the miRNA-mediated post-transcriptional level.

Clinically important anti-cancer drugs vinblastine and vincristine are solely synthesized by the terpenoid indole alkaloid (TIA) pathway in Catharanthus roseus. Anthranilate synthase (AS) is a rate-limiting enzyme in the TIA pathway. The transgenic C. roseus hairy root line overexpressing a feedback insensitive ASα subunit under the control of an inducible promoter and the ASβ subunit constitutively was previously created for the overproduction of TIAs. However, both increases and decreases in TIAs were detected after overexpressing ASα. Although genetic modification is targeted to one gene in the TIA pathway, it could trigger global transcriptional changes that can directly or indirectly affect TIA biosynthesis. In this study, Illumina sequencing and RT-qPCR were used to detect the transcriptional responses to overexpressing AS, which can increase understanding of the complex regulation of the TIA pathway and further inspire rational metabolic engineering for enhanced TIA production in C. roseus hairy roots. Overexpressing AS in C. roseus hairy roots altered the transcription of most known TIA pathway genes and regulators after 12, 24, and 48 h induction detected by RT-qPCR. Changes in the transcriptome of C. roseus hairy roots was further investigated 18 hours after ASα induction and compared to the control hairy roots using RNA-seq. A unigene set of 30,281 was obtained by de novo assembly of the sequencing reads. Comparison of the differentially expressed transcriptional profiles resulted in 2853 differentially expressed transcripts. Functional annotation of these transcripts revealed a complex and systematically transcriptome change in ASαβ hairy roots. Pathway analysis shows alterations in many pathways such as aromatic amino acid biosynthesis, jasmonic acid (JA) biosynthesis and other secondary metabolic pathways after perturbing AS. Moreover, many genes in overall stress response were differentially expressed after overexpressing ASα. The transcriptomic analysis illustrates overexpressing AS stimulates the overall stress response and affects the metabolic networks in C. roseus hairy roots. The up-regulation of endogenous JA biosynthesis pathway indicates the involvement of JA signal transduction to regulate TIA biosynthesis in ASαβ engineered roots and explained why many of the transcripts for TIA genes and regulators are seen to increase with AS overexpression.

Background Single-cell RNA sequencing (scRNA-seq) measurements of gene expression show great promise for studying the cellular heterogeneity of rice roots. How precisely annotating cell identity is a major unresolved problem in plant scRNA-seq analysis due to the inherent high dimensionality and sparsity. Results To address this challenge, we present NRTPredictor, an ensemble-learning system, to predict rice root cell stage and mine biomarkers through complete model interpretability. The performance of NRTPredictor was evaluated using a test dataset, with 98.01% accuracy and 95.45% recall. With the power of interpretability provided by NRTPredictor, our model recognizes 110 marker genes partially involved in phenylpropanoid biosynthesis. Expression patterns of rice root could be mapped by the above-mentioned candidate genes, showing the superiority of NRTPredictor. Integrated analysis of scRNA and bulk RNA-seq data revealed aberrant expression of Epidermis cell subpopulations in flooding, Pi, and salt stresses. Conclusion Taken together, our results demonstrate that NRTPredictor is a useful tool for automated prediction of rice root cell stage and provides a valuable resource for deciphering the rice root cellular heterogeneity and the molecular mechanisms of flooding, Pi, and salt stresses. Based on the proposed model, a free webserver has been established, which is available at https://www.cgris.net/nrtp .

Summary Cotton fibre is a unicellular seed trichome, and lint fibre initials per seed as a factor determines fibre yield. However, the mechanisms controlling fibre initiation from ovule epidermis are not understood well enough. Here, with single‐cell RNA sequencing (scRNA‐seq), a total of 14 535 cells were identified from cotton ovule outer integument of Xu142_LF line at four developmental stages (1.5, 1, 0.5 days before anthesis and the day of anthesis). Three major cell types, fibre, non‐fibre epidermis and outer pigment layer were identified and then verified by RNA in situ hybridization. A comparative analysis on scRNA‐seq data between Xu142 and its fibreless mutant Xu142 fl further confirmed fibre cluster definition. The developmental trajectory of fibre cell was reconstructed, and fibre cell was identified differentiated at 1 day before anthesis. Gene regulatory networks at four stages revealed the spatiotemporal pattern of core transcription factors, and MYB25‐like and HOX3 were demonstrated played key roles as commanders in fibre differentiation and tip‐biased diffuse growth respectively. A model for early development of a single fibre cell was proposed here, which sheds light on further deciphering mechanism of plant trichome and the improvement of cotton fibre yield.

Introduction Triptolide (TPL) is a promising plant-derived compound for clinical therapy of multiple human diseases; however, its application was limited considering its toxicity. Methods To explore the underlying molecular mechanism of TPL nephrotoxicity, a network pharmacology based approach was utilized to predict candidate targets related with TPL toxicity, followed by deep RNA-seq analysis to characterize the features of three transcriptional elements include protein coding genes (PCGs), long noncoding RNAs (lncRNAs) and circular RNAs (circRNAs) as well as their associations with nephrotoxicity in rats with TPL treatment. Results & Discussion Although the deeper mechanisms of TPL nephrotoxcity remain further exploration, our results suggested that c-Jun is a potential target of TPL and Per1 related circadian rhythm signaling is involved in TPL induced renal toxicity.

Plant-microbe symbioses require intense interaction and genetic coordination to successfully establish in specific cell types of the host and symbiont. Traditional RNA-seq methodologies lack the cellular resolution to fully capture these complexities, but single-cell and spatial transcriptomics (ST) are now allowing scientists to probe symbiotic interactions at an unprecedented level of detail. Here, we discuss the advantages that novel spatial and single-cell transcriptomic technologies provide in studying plant-microbe endosymbioses and highlight key recent studies. Finally, we consider the remaining limitations of applying these approaches to symbiosis research, which are mainly related to the simultaneous capture of both plant and microbial transcripts within the same cells.

BackgroundPrimary roots (radicles) represent the first visible developmental stages of the plant and are crucial for nutrient supply and the integration of environmental signals. Few studies have analyzed primary roots at a molecular level, and were mostly limited to Arabidopsis. Here we study the primary root transcriptomes of standard type, heterozygous columnar and homozygous columnar apple (Malus x domestica) by RNA-Seq and quantitative real-time PCR. The columnar growth habit is characterized by a stunted main axis and the development of short fruit spurs instead of long lateral branches. This compact growth possesses economic potential because it allows high density planting and mechanical harvesting of the trees. Its molecular basis has been identified as a nested Gypsy-44 retrotransposon insertion; however the link between the insertion and the phenotype as well as the timing of the phenotype emergence are as yet unclear. We extend the transcriptomic studies of columnar tissues to the radicles, which are the earliest developmental stage and investigate whether homozygous columnar seedlings are viable.ResultsRadicles mainly express genes associated with primary metabolism, growth and development. About 200 genes show differential regulation in a comparison of heterozygous columnar radicles with non-columnar radicles, whereas the comparison of homozygous columnar radicles with non-columnar radicles yields about 300 differentially regulated genes. Genes involved in cellulose and phenylpropanoid biosynthesis, cell wall modification, transcription and translation, ethylene and jasmonate biosynthesis are upregulated in columnar radicles. Genes in the vicinity of the columnar-specific Gypsy-44 insertion experience an especially strong differential regulation: the direct downstream neighbor, dmr6-like, is downregulated in heterozygous columnar radicles, but strongly upregulated in columnar shoot apical meristems.ConclusionsThe transcriptomic profile of primary roots reflects their pivotal role in growth and development. Homozygous columnar embryos are viable and form normal radicles under natural conditions, and selection towards heterozygous plants most likely occurs due to breeders’ preferences. Cell wall and phytohormone biosynthesis and metabolism experience differential regulation in columnar radicles. Presumably the first step of the differential regulation most likely happens within the region of the retrotransposon insertion and its tissue-specificity suggests involvement of one (or several) tissue-specific regulator(s).

BackgroundThe genus Cuscuta is a group of parasitic plants that are distributed world-wide. The process of parasitization starts with a Cuscuta plant coiling around the host stem. The parasite’s haustorial organs then establish a vascular connection allowing for access to the phloem content. The host and the parasite form new cellular connections, suggesting coordination of developmental and biochemical processes. Simultaneous monitoring of gene expression in the parasite’s and host’s tissues may shed light on the complex events occurring between the parasitic and host cells and may help to overcome experimental limitations (i.e. how to separate host tissue from Cuscuta tissue at the haustorial connection). A novel approach is to use bioinformatic analysis to classify sequencing reads as either belonging to the host or to the parasite and to characterize the expression patterns. Owing to the lack of a comprehensive genomic dataset from Cuscuta spp., such a classification has not been performed previously.ResultsWe first classified RNA-Seq reads from an interface region between the non-model parasitic plant Cuscuta japonica and the non-model host plant Impatiens balsamina. Without established reference sequences, we classified reads as originating from either of the plants by stepwise similarity search against de novo assembled transcript sets of C. japonica and I. balsamina, unigene sets of the same genus, and cDNA sequences of the same family. We then assembled de novo transcriptomes from the classified read sets. We assessed the quality of the classification by mapping reads to contigs of both plants, achieving a misclassification rate low enough (0.22-0.39%) to be used reliably for differential gene expression analysis. Finally, we applied our read classification method to RNA-Seq data from the interface between the non-model parasitic plant C. japonica and the model host plant Glycine max. Analysis of gene expression profiles at 5 parasitizing stages revealed differentially expressed genes from both C. japonica and G. max, and uncovered the coordination of cellular processes between the two plants.ConclusionsWe demonstrated that reliable identification of differentially expressed transcripts in undissected interface region of the parasite-host association is feasible and informative with respect to differential-expression patterns.

Premise of the study: Hyb-Seq, the combination of target enrichment and genome skimming, allows simultaneous data collection for low-copy nuclear genes and high-copy genomic targets for plant systematics and evolution studies. Methods and Results: Genome and transcriptome assemblies for milkweed (Asclepias syriaca) were used to design enrichment probes for 3385 exons from 768 genes (>1.6 Mbp) followed by Illumina sequencing of enriched libraries. Hyb-Seq of 12 individuals (10 Asclepias species and two related genera) resulted in at least partial assembly of 92.6% of exons and 99.7% of genes and an average assembly length >2 Mbp. Importantly, complete plastomes and nuclear ribosomal DNA cistrons were assembled using off-target reads. Phylogenomic analyses demonstrated signal conflict between genomes. Conclusions: The Hyb-Seq approach enables targeted sequencing of thousands of low-copy nuclear exons and flanking regions, as well as genome skimming of high-copy repeats and organellar genomes, to efficiently produce genome-scale data sets for phylogenomics.

With the increasing availability of massive transcriptome data in plants, it is possible to construct a comprehensive database with multiple transcriptome data and online imputation tools. Online databases with integrated, multifaceted functions for exploring published RNA-seq libraries or microarrays are available for several plant species (Darwish et al., 2013; Xia et al., 2017; Zhang et al., 2020), helping researchers take advantage of the vast collection of public datasets. Brassica napus (AACC, 2n = 38) is one of the most important oil crops, originating from a spontaneous hybridization between Brassica rapa (AA, 2n = 20) and Brassica oleracea (CC, 2n = 18) (Chalhoub et al., 2014). Recently, an online transcriptome database has been built in B. napus (Chao et al., 2020). This database includes transcriptomes from different labs with various genotypes and growth conditions, which affects the usability. A database with comprehensive transcriptomes from one B. napus genotype is greatly needed.

Summary Single‐cell RNA‐seq (scRNA‐seq) has been highlighted as a powerful tool for the description of human cell transcriptome, but the technology has not been broadly applied in plant cells. Herein, we describe the successful development of a robust protoplast cell isolation system in the peanut leaf. A total of 6,815 single cells were divided into eight cell clusters based on reported marker genes by applying scRNA‐seq. Further, a pseudo‐time analysis was used to describe the developmental trajectory and interaction network of transcription factors (TFs) of distinct cell types during leaf growth. The trajectory enabled re‐investigation of the primordium‐driven development processes of the mesophyll and epidermis. These results suggest that palisade cells likely differentiate into spongy cells, while the epidermal cells originated earlier than the primordium. Subsequently, the developed method integrated multiple technologies to efficiently validate the scRNA‐seq result in a homogenous cell population. The expression levels of several TFs were strongly correlated with epidermal ontogeny in accordance with obtained scRNA‐seq values. Additionally, peanut AHL23 (AT‐HOOK MOTIF NUCLEAR LOCALIZED PROTEIN 23), which is localized in nucleus, promoted leaf growth when ectopically expressed in Arabidopsis by modulating the phytohormone pathway. Together, our study displays that application of scRNA‐seq can provide new hypotheses regarding cell differentiation in the leaf blade of Arachis hypogaea. We believe that this approach will enable significant advances in the functional study of leaf blade cells in the allotetraploid peanut and other plant species.

Abstract Growing plants with modified cell wall compositions is a promising strategy to improve resistance to pathogens, increase biomass digestibility, and tune other important properties. In order to alter biomass architecture, a detailed knowledge of cell wall structure and biosynthesis is a prerequisite. We report here a glycan array‐based assay for the high‐throughput identification and characterization of plant cell wall biosynthetic glycosyltransferases (GTs). We demonstrate that different heterologously expressed galactosyl‐, fucosyl‐, and xylosyltransferases can transfer azido‐functionalized sugar nucleotide donors to selected synthetic plant cell wall oligosaccharides on the array and that the transferred monosaccharides can be visualized “on chip” by a 1,3‐dipolar cycloaddition reaction with an alkynyl‐modified dye. The opportunity to simultaneously screen thousands of combinations of putative GTs, nucleotide sugar donors, and oligosaccharide acceptors will dramatically accelerate plant cell wall biosynthesis research.

Studying development and physiology of growing roots is challenging due to limitations regarding cellular and subcellular analysis under controlled environmental conditions. We describe a microfluidic chip platform, called RootChip, that integrates live-cell imaging of growth and metabolism of Arabidopsis thaliana roots with rapid modulation of environmental conditions. The RootChip has separate chambers for individual regulation of the microenvironment of multiple roots from multiple seedlings in parallel. We demonstrate the utility of The RootChip by monitoring time-resolved growth and cytosolic sugar levels at subcellular resolution in plants by a genetically encoded fluorescence sensor for glucose and galactose. The RootChip can be modified for use with roots from other plant species by adapting the chamber geometry and facilitates the systematic analysis of root growth and metabolism from multiple seedlings, paving the way for large-scale phenotyping of root metabolism and signaling.

No abstract available

No abstract available

No abstract available

Gene activity is largely controlled by transcriptional regulation through the action of transcription factors and other regulators. QsMYB1 is a member of the R2R3-MYB transcription factor family related to secondary growth, and in particular, with the cork development process. In order to identify the putative gene targets of QsMYB1 across the cork oak genome we developed a ChIP-Seq strategy. Results provide direct evidence that QsMY1B targets genes encoding for enzymes involved in the lignin and suberin pathways as well as gene encoding for ABCG transporters and LTPs implicated in the transport of monomeric suberin units across the cellular membrane. These results highlight the role of QsMYB1 as a regulator of lignin and suberin biosynthesis, transport and assembly. To our knowledge, this work constitutes the first ChIP-Seq experiment performed in cork oak, a non-model plant species with a long-life cycle, and these results will contribute to deepen the knowledge about the molecular mechanisms of cork formation and differentiation.

The bacterium Agrobacterium tumefaciens has been the workhorse in plant genome engineering. Customized replacement of native tumor-inducing (Ti) plasmid elements enabled insertion of a sequence of interest called Transfer-DNA (T-DNA) into any plant genome. Although these transfer mechanisms are well understood, detailed understanding of structure and epigenomic status of insertion events was limited by current technologies. Here we applied two single-molecule technologies and analyzed Arabidopsis thaliana lines from three widely used T-DNA insertion collections (SALK, SAIL and WISC). Optical maps for four randomly selected T-DNA lines revealed between one and seven insertions/rearrangements, and the length of individual insertions from 27 to 236 kilobases. De novo nanopore sequencing-based assemblies for two segregating lines partially resolved T-DNA structures and revealed multiple translocations and exchange of chromosome arm ends. For the current TAIR10 reference genome, nanopore contigs corrected 83% of non-centromeric misassemblies. The unprecedented contiguous nucleotide-level resolution enabled an in-depth study of the epigenome at T-DNA insertion sites. SALK_059379 line T-DNA insertions were enriched for 24nt small interfering RNAs (siRNA) and dense cytosine DNA methylation, resulting in transgene silencing via the RNA-directed DNA methylation pathway. In contrast, SAIL_232 line T-DNA insertions are predominantly targeted by 21/22nt siRNAs, with DNA methylation and silencing limited to a reporter, but not the resistance gene. Additionally, we profiled the H3K4me3, H3K27me3 and H2A.Z chromatin environments around T-DNA insertions using ChIP-seq in SALK_059379, SAIL_232 and five additional T-DNA lines. We discovered various effect s ranging from complete loss of chromatin marks to the de novo incorporation of H2A.Z and trimethylation of H3K4 and H3K27 around the T-DNA integration sites. This study provides new insights into the structural impact of inserting foreign fragments into plant genomes and demonstrates the utility of state-of-the-art long-range sequencing technologies to rapidly identify unanticipated genomic changes.

No abstract available

Abstract Background and Aims Auxin response factors (ARFs) as transcription activators or repressors have important roles in plant growth and development, but knowledge about the functions of wheat ARF members is limited. A novel ARF member in wheat (Triticum aestivum), TaARF4, was identified, and its protein function, haplotype geographic distribution and allelic frequencies were investigated. Methods Tissue expression of TaARF4 was analysed by real-time PCR. Sub-cellular localization was performed using green fluorescent protein (GFP)-tagged TaARF4. Ectopic expression of TaARF4-A in arabidopsis was used to study its functions. Electrophoretic mobility shift assays (EMSAs), chromatin immunoprecipitation (ChIP) analyses and gene expression were performed to detect TaARF4 target genes. A dCAPS (derived cleaved amplified polymorphic sequence) marker developed from TaARF4-B was used to identify haplotypes and association analysis between haplotypes and agronomic traits. Key Results TaARF4-A was constitutively expressed and its protein was localized in the nucleus. Ectopic expression of TaARF4-A in arabidopsis caused abscisic acid (ABA) insensitivity, shorter primary root length and reduced plant height (PH). Through expression studies and ChIP assays, TaARF4-A was shown to regulate HB33 expression which negatively responded to ABA, and reduced root length and plant height by repressing expression of Gretchen Hagen 3 (GH3) genes that in turn upregulated indole-3-acetic acid content in arabidopsis. Association analysis showed that TaARF4-B was strongly associated with PH and root depth at the tillering, jointing and grain fill stages. Geographic distribution and allelic frequencies suggested that TaARF4-B haplotypes were selected in Chinese wheat breeding programmes. An amino acid change (threonine to alanine) at position 158 might be the cause of phenotype variation in accessions possessing different haplotypes. Conclusions Ectopic expression and association analysis indicate that TaARF4 may be involved in root length and plant height determination in wheat. This work is helpful for selection of wheat genotypes with optimal root and plant architecture.

Protein ubiquitylation profoundly expands proteome functionality and diversifies cellular signaling processes, with recent studies providing ample evidence for its importance to plant immunity. To gain a proteome-wide appreciation of ubiquitylome dynamics during immune recognition, we employed a two-step affinity enrichment protocol based on a 6His-tagged ubiquitin (Ub) variant coupled with high sensitivity mass spectrometry to identify Arabidopsis proteins rapidly ubiquitylated upon plant perception of the microbe-associated molecular pattern (MAMP) peptide flg22. The catalog from two-week-old seedlings treated for only 30 minutes with flg22 contained nearly 1,000 conjugates, 150 Ub footprints, and all seven types of Ub linkages, and included previously uncharacterized conjugates of immune components, such as RECEPTOR-LIKE KINASE 1 (RKL1) shown to negatively regulate plant immunity. In vivo ubiquitylation assays confirmed modification of several candidates upon immune elicitation, and revealed distinct modification patterns and dynamics for key immune components, including poly- and monoubiquitylation, as well as induced or reduced levels of ubiquitylation. Gene ontology and network analyses of the collection also uncovered rapid modification of the Ub-proteasome system itself, suggesting a critical auto-regulatory loop necessary for an effective MAMP-triggered immune response and subsequent disease resistance. Included targets were UBIQUITIN-CONJUGATING ENZYME 13 (UBC13) and proteasome component REGULATORY PARTICLE NON-ATPASE SUBUNIT 8b (RPN8b), whose subsequent biochemical and genetic analyses implied negative roles in immune elicitation. Collectively, our proteomic analyses further strengthened the connection between ubiquitylation and flg22-based immune signaling, identified novel components and pathways regulating plant immunity, and increased the database of ubiquitylated substrates in plants. One-sentence summary Proteome-wide catalogs of ubiquitylated proteins revealed a rapid engagement of the ubiquitin-proteasome system in Arabidopsis innate immunity.

No abstract available

No abstract available

No abstract available

Climate variability due to fluctuation in temperature is a worldwide concern that imperils crop production. The need to understand how the germplasm variation in major crops can be utilized to aid in discovering and developing breeding lines that can withstand and adapt to temperature fluctuations is more necessary than ever. Here, we analyzed the genetic variation associated with responses to thermal stresses in a sorghum association panel (SAP) representing major races and working groups to identify single nucleotide polymorphisms (SNPs) that are associated with resilience to temperature stress in a major cereal crop. The SAP exhibited extensive variation for seedling traits under cold and heat stress. Genome-wide analyses identified 30 SNPs that were strongly associated with traits measured at seedling stage under cold stress and tagged genes that act as regulators of anthocyanin expression and soluble carbohydrate metabolism. Meanwhile, 12 SNPs were significantly associated with seedling traits under heat stress and these SNPs tagged genes that function in sugar metabolism, and ion transport pathways. Evaluation of co-expression networks for genes near the significantly associated SNPs indicated complex gene interactions for cold and heat stresses in sorghum. We focused and validated the expression of four genes in the network of Sb06g025040, a basic-helix-loop-helix (bHLH) transcription factor that was proposed to be involved in purple color pigmentation of leaf, and observed that genes in this network were upregulated during cold stress in a moderately tolerant line as compared to the more sensitive line. This study facilitated the tagging of genome regions associated with variation in seedling traits of sorghum under cold and heat stress. These findings show the potential of genotype information for development of temperature resilient sorghum cultivars and further characterization of genes and their networks responsible for adaptation to thermal stresses. Knowledge on the gene networks from this research can be extended to the other cereal crops to better understand the genetic basis of resilience to temperature fluctuations during plant developmental stages.

No abstract available

In petunia flowers, the loci an1, an2, and an11 control the pigmentation of the flower by stimulating the transcription of anthocyanin biosynthetic genes. The an1 and an2 locus were recently cloned and encode a basic helix-loop-helix (bHLH) and MYB-domain transcriptional activator, respectively. Here, we report the isolation of the an11 locus by transposon tagging. RNA gel blot experiments show that an11 is expressed independently from an1 and an2 throughout plant development, as well as in tissues that do not express the anthocyanin pathway. It encodes a novel WD-repeat protein that is highly conserved even in species that do not produce anthocyanins such as yeast, nematodes, and mammals. The observation that the human an11 homolog partially complements the an11 petunia mutant in transient assays shows that sequence similarity reflects functional conservation. Overexpression of an2 in an11- petals restored the activity of a structural anthocyanin gene in transient assays, indicating that AN11 acts upstream of AN2. Cell fractionation experiments show that the bulk of the AN11 protein is localized in the cytoplasm. Taken together, this indicates that AN11 is a cytoplasmic component of a conserved signal transduction cascade that modulates AN2 function in petunia, thereby linking cellular signals with transcriptional activation.

This work provides the first three-dimensional structure of a member of the plant annexin family and correlates these findings with biochemical properties of this protein. Annexin 24(Ca32) from Capsicum annuum was purified as a native protein from bell pepper and was also prepared by recombinant techniques. To overcome the problem of precipitation of the recombinant wild-type protein in crystallization trials, two mutants were designed. Whereas an N-terminal truncation mutant turned out to be an unstable protein, the N-terminal His-tagged annexin 24(Ca32) was crystallized, and the three-dimensional structure was determined by x-ray diffraction at 2.8 Å resolution. The structure refined to an R-factor of 0.216 adopts the typical annexin fold; the detailed structure, however, is different from non-plant annexins, especially in domains I and III and in the membrane binding loops on the convex side. Within the unit cell there are two molecules per asymmetric unit, which differ in conformation of the IAB-loop. Both conformers show Trp-35 on the surface. The loop-out conformation is stabilized by tight interactions of this tryptophan with residue side chains of a symmetry-related molecule and enforced by a bound sulfate. Characterization of this plant annexin using biophysical methods revealed calcium-dependent binding to phospholipid vesicles with preference for phosphatidylcholine over phosphatidylserine and magnesium-dependent phosphodiesterase activity in vitro as shown with adenosine triphosphate as the substrate. A comparative unfolding study of recombinant annexin 24(Ca32) wild type and of the His-tag fusion protein indicates higher stability of the latter. The effect of this N-terminal modification is also visible from CD spectra. Both proteins were subjected to a FURA-2-based calcium influx assay, which gave high influx rates for the wild-type but greatly reduced influx rates for the fusion protein. We therefore conclude that the N-terminal domain is indeed a major regulatory element modulating different annexin properties by allosteric mechanisms.

No abstract available

DNA double-strand breaks (DSBs) are toxic DNA lesions, which, if not properly repaired, may lead to genomic instability, cell death and senescence. Damage-induced long non-coding RNAs (dilncRNAs) are transcribed from broken DNA ends and contribute to DNA damage response (DDR) signaling. Here we show that dilncRNAs play a role in DSB repair by homologous recombination (HR) by contributing to the recruitment of the HR proteins BRCA1, BRCA2, and RAD51, without affecting DNA-end resection. In S/G2-phase cells, dilncRNAs pair to the resected DNA ends and form DNA:RNA hybrids, which are recognized by BRCA1. We also show that BRCA2 directly interacts with RNase H2, mediates its localization to DSBs in the S/G2 cell-cycle phase, and controls DNA:RNA hybrid levels at DSBs. These results demonstrate that regulated DNA:RNA hybrid levels at DSBs contribute to HR-mediated repair. Long non-coding RNAs transcribed at DNA damaged sites can play part in DNA damage response. Here the authors reveal that damaged induced lncRNAs can form DNA:RNA hybrids at resected DNA-ends. These hybrids are involved in recruiting HR-mediated repair machinery which, in turn, controls their level at DSBs.

DNA:RNA hybrid formation is emerging as a significant cause of genome instability in biological systems ranging from bacteria to mammals. Here we describe the genome-wide distribution of DNA:RNA hybrid prone loci in Saccharomyces cerevisiae by DNA:RNA immunoprecipitation (DRIP) followed by hybridization on tiling microarray. These profiles show that DNA:RNA hybrids preferentially accumulated at rDNA, Ty1 and Ty2 transposons, telomeric repeat regions and a subset of open reading frames (ORFs). The latter are generally highly transcribed and have high GC content. Interestingly, significant DNA:RNA hybrid enrichment was also detected at genes associated with antisense transcripts. The expression of antisense-associated genes was also significantly altered upon overexpression of RNase H, which degrades the RNA in hybrids. Finally, we uncover mutant-specific differences in the DRIP profiles of a Sen1 helicase mutant, RNase H deletion mutant and Hpr1 THO complex mutant compared to wild type, suggesting different roles for these proteins in DNA:RNA hybrid biology. Our profiles of DNA:RNA hybrid prone loci provide a resource for understanding the properties of hybrid-forming regions in vivo, extend our knowledge of hybrid-mitigating enzymes, and contribute to models of antisense-mediated gene regulation. A summary of this paper was presented at the 26th International Conference on Yeast Genetics and Molecular Biology, August 2013.

The kinetic properties of Escherichia coli ribonuclease H (RNase H) were investigated using oligonucleotide substrates that consist of a short stretch of RNA, flanked on either side by DNA (DNA-RNA-DNA). In the presence of a complementary DNA strand, RNase H cleavage is restricted to the short ribonucleotide stretch of the DNA/RNA heteroduplex. The DNA-RNA-DNA substrate utilized for kinetic studies: (formula; see text) is cleaved at a single site (decreases) in the presence of a complementary DNA strand, to generate (dT)7-(rA)2-OH and p-(rA)2-(dT)9. Anion exchange high performance liquid chromatography was used to separate and quantitate the cleavage products. Under these conditions, RNase H-specific and nonspecific degradation products could be resolved. Kinetic parameters were measured under conditions of 100% hybrid formation (1.2-1.5 molar excess of complementary DNA, T much less than Tm). A linear double reciprocal plot was obtained, yielding a Km of 4.2 microM and a turnover number of 7.1 cleavages per s per RNase H monomer. The kinetic properties of substrate analogs containing varying lengths of RNA (n = 3-5) and 2'-O-methyl modifications were also investigated. Maximal turnover was observed with DNA-RNA-DNA substrates containing a minimum of four RNA residues. Kcat for the rA3 derivative was decreased by more than 100-fold. The Km appeared to decrease with the size of the internal RNA stretch (n = 3-5). No significant difference in turnover number of Km was observed when the flanking DNA was replaced with 2'-O-methyl RNA, suggesting that RNase H does not interact with this region of the heteroduplex.

No abstract available

No abstract available