Insights into cryptochrome modulation of ABA signaling to mediate dormancy regulation in Marchantia polymorpha.

The acquisition of dormancy capabilities has enabled plants to survive in adverse terrestrial environmental conditions. Dormancy accumulation and release is coupled with light signaling, which is well studied in Arabidopsis, but it is unclear in the distant nonvascular relative. We study the characteristics and function on dormancy regulation of a blue light receptor cryptochrome in Marchantia polymorpha (MpCRY). Here, we identified MpCRY via bioinformatics and mutant complement analysis. The biochemical characteristics were assessed by multiple protein-binding assays. The function of MpCRY in gemma dormancy was clarified by overexpression and mutation of MpCRY, and its mechanism was analyzed via RNA sequencing and quantitative PCR analyses associated with hormone treatment. We found that the unique MpCRY protein in M. polymorpha undergoes both blue light-promoted interaction with itself (self-interaction) and blue light-dependent phosphorylation. MpCRY has the specific characteristics of blue light-induced nuclear localization and degradation. We further demonstrated that MpCRY transcriptionally represses abscisic acid (ABA) signaling-related gene expression to suppress gemma dormancy, which is dependent on blue light signaling. Our findings indicate that MpCRY possesses specific biochemical and molecular characteristics, and modulates ABA signaling under blue light conditions to regulate gemma dormancy in M. polymorpha.

Whole-genome characterization of chronological age-associated changes in methylome and circular RNAs in moso bamboo (Phyllostachys edulis) from vegetative to floral growth.

In mammals, DNA methylation is associated with aging. However, age-related DNA methylation changes during phase transitions largely remain unstudied in plants. Moso bamboo (Phyllostachys edulis) requires a very long time to transition from the vegetative to the floral phase. To comprehensively investigate the association of DNA methylation with aging, we present here single-base-resolution DNA methylation profiles using both high-throughput bisulfite sequencing and single-molecule nanopore-based DNA sequencing, covering the long period of vegetative growth and transition to flowering in moso bamboo. We discovered that CHH methylation gradually accumulates from vegetative to reproductive growth in a time-dependent fashion. Differentially methylated regions, correlating with chronological aging, occurred preferentially at both transcription start sites and transcription termination sites. Genes with CG methylation changes showed an enrichment of Gene Ontology (GO) categories in 'vegetative to reproductive phase transition of meristem'. Combining methylation data with mRNA sequencing revealed that DNA methylation in promoters, introns and exons may have different roles in regulating gene expression. Finally, circular RNA (circRNA) sequencing revealed that the flanking introns of circRNAs are hypermethylated and enriched in long terminal repeat (LTR) retrotransposons. Together, the observations in this study provide insights into the dynamic DNA methylation and circRNA landscapes, correlating with chronological age, which paves the way to study further the impact of epigenetic factors on flowering in moso bamboo.

Profiling of circular RNA N6 -methyladenosine in moso bamboo (Phyllostachys edulis) using nanopore-based direct RNA sequencing.

N6 -methyladenosine (m6 A) is a prevalent modification in messenger RNAs and circular RNAs that play important roles in regulating various aspects of RNA metabolism. However, the occurrence of the m6 A modification in plant circular RNAs has not been reported. A widely used method to identify m6 A modifications relies on m6 A-specific antibodies followed by next-generation sequencing of precipitated RNAs (MeRIP-Seq). However, one limitation of MeRIP-Seq is that it does not provide the precise location of m6 A at single-nucleotide resolution. Although more recent sequencing techniques such as Nanopore-based direct RNA sequencing (DRS) can overcome such limitations, the technology does not allow sequencing of circular RNAs, as these molecules lack a poly(A) tail. Here, we developed a novel method to detect the precise location of m6 A modifications in circular RNAs using Nanopore DRS. We first enriched our samples for circular RNAs, which we then fragmented and sequenced on the Nanopore platform with a customized protocol. Using this method, we identified 470 unique circular RNAs from DRS reads based on the back-spliced junction region. Among exonic circular RNAs, about 10% contained m6 A sites, which mainly occurred around acceptor and donor splice sites. This study demonstrates the utility of our antibody-independent method in identifying total and methylated circular RNAs using Nanopore DRS. This method has the additional advantage of providing the exact location of m6 A sites at single-base resolution in circular RNAs or linear transcripts from non-coding RNA without poly(A) tails.

Genome-Wide Profiling of Circular RNAs in the Rapidly Growing Shoots of Moso Bamboo (Phyllostachys edulis).

Circular RNAs, including circular exonic RNAs (circRNA), circular intronic RNAs (ciRNA) and exon-intron circRNAs (EIciRNAs), are a new type of noncoding RNAs. Growing shoots of moso bamboo (Phyllostachys edulis) represent an excellent model of fast growth and their circular RNAs have not been studied yet. To understand the potential regulation of circular RNAs, we systematically characterized circular RNAs from eight different developmental stages of rapidly growing shoots. Here, we identified 895 circular RNAs including a subset of mutually inclusive circRNA. These circular RNAs were generated from 759 corresponding parental coding genes involved in cellulose, hemicellulose and lignin biosynthetic process. Gene co-expression analysis revealed that hub genes, such as DEFECTIVE IN RNA-DIRECTED DNA METHYLATION 1 (DRD1), MAINTENANCE OF METHYLATION (MOM), dicer-like 3 (DCL3) and ARGONAUTE 1 (AGO1), were significantly enriched giving rise to circular RNAs. The expression level of these circular RNAs presented correlation with its linear counterpart according to transcriptome sequencing. Further protoplast transformation experiments indicated that overexpressing circ-bHLH93 generating from transcription factor decreased its linear transcript. Finally, the expression profiles suggested that circular RNAs may have interplay with miRNAs to regulate their cognate linear mRNAs, which was further supported by overexpressing miRNA156 decreasing the transcript of circ-TRF-1 and linear transcripts of TRF-1. Taken together, the overall profile of circular RNAs provided new insight into an unexplored category of long noncoding RNA regulation in moso bamboo.

Neighbor Detection Induces Organ-Specific Transcriptomes, Revealing Patterns Underlying Hypocotyl-Specific Growth.

In response to neighbor proximity, plants increase the growth of specific organs (e.g., hypocotyls) to enhance access to sunlight. Shade enhances the activity of Phytochrome Interacting Factors (PIFs) by releasing these bHLH transcription factors from phytochrome B-mediated inhibition. PIFs promote elongation by inducing auxin production in cotyledons. In order to elucidate spatiotemporal aspects of the neighbor proximity response, we separately analyzed gene expression patterns in the major light-sensing organ (cotyledons) and in rapidly elongating hypocotyls of Arabidopsis thaliana PIFs initiate transcriptional reprogramming in both organs within 15 min, comprising regulated expression of several early auxin response genes. This suggests that hypocotyl growth is elicited by both local and distal auxin signals. We show that cotyledon-derived auxin is both necessary and sufficient to initiate hypocotyl growth, but we also provide evidence for the functional importance of the local PIF-induced response. With time, the transcriptional response diverges increasingly between organs. We identify genes whose differential expression may underlie organ-specific elongation. Finally, we uncover a growth promotion gene expression signature shared between different developmentally regulated growth processes and responses to the environment in different organs.

PtWOX11 acts as master regulator conducting the expression of key transcription factors to induce de novo shoot organogenesis in poplar.

KEY MESSAGE: WUSCHEL-RELATED HOMEOBOX 11 establishes the acquisition of pluripotency during callus formation and accomplishes de novo shoot formation by regulating key transcription factors in poplar. De novo shoot regeneration is a prerequisite for propagation and genetic engineering of elite cultivars in forestry. However, the regulatory mechanism of de novo organogenesis is poorly understood in tree species. We previously showed that WUSCHEL (WUS)-RELATED HOMEOBOX 11 (PtWOX11) of the hybrid poplar clone 84K (Populus alba × P. glandulosa) promotes de novo root formation. In this study, we found that PtWOX11 also regulates de novo shoot regeneration in poplar. The overexpression of PtWOX11 enhanced de novo shoot formation, whereas overexpression of PtWOX11 fused with the transcriptional repressor domain (PtWOX11-SRDX) or reduced expression of PtWOX11 inhibited this process, indicating that PtWOX11 promotes de novo shoot organogenesis. Although PtWOX11 promotes callus formation, overexpression of PtWOX11 and PtWOX11-SRDX also produced increased and decreased numbers of de novo shoots per unit weight, respectively, implying that PtWOX11 promotes de novo shoot organogenesis partially by regulating the intrinsic mechanism of shoot development. RNA-seq and qPCR analysis further revealed that PtWOX11 activates the expression of PLETHORA1 (PtPLT1) and PtPLT2, whose Arabidopsis paralogs establish the acquisition of pluripotency, during incubation on callus-inducing medium. Moreover, PtWOX11 activates the expression of shoot-promoting factors and meristem regulators such as CUP-SHAPED COTYLEDON2 (PtCUC2), PtCUC3, WUS and SHOOT MERISTEMLESS to fulfill shoot regeneration during incubation on shoot-inducing medium. These results suggest that PtWOX11 acts as a central regulator of the expression of key genes to cause de novo shoot formation. Our studies further provide a possible means to genetically engineer economically important tree species for their micropropagation.

Phytochrome interacting factors 4 and 5 control seedling growth in changing light conditions by directly controlling auxin signaling.

Plant growth is strongly influenced by the presence of neighbors that compete for light resources. In response to vegetational shading shade-intolerant plants such as Arabidopsis display a suite of developmental responses known as the shade-avoidance syndrome (SAS). The phytochrome B (phyB) photoreceptor is the major light sensor to mediate this adaptive response. Control of the SAS occurs in part with phyB, which controls protein abundance of phytochrome-interacting factors 4 and 5 (PIF4 and PIF5) directly. The shade-avoidance response also requires rapid biosynthesis of auxin and its transport to promote elongation growth. The identification of genome-wide PIF5-binding sites during shade avoidance revealed that this bHLH transcription factor regulates the expression of a subset of previously identified SAS genes. Moreover our study suggests that PIF4 and PIF5 regulate elongation growth by controlling directly the expression of genes that code for auxin biosynthesis and auxin signaling components.

Multi-omics of Circular RNAs and Their Responses to Hormones in Moso Bamboo (Phyllostachys edulis).

Circular RNAs (circRNAs) are endogenous non-coding RNAs with covalently closed structures, which have important functions in plants. However, their biogenesis, degradation, and function upon treatment with gibberellins (GAs) and auxins (1-naphthaleneacetic acid, NAA) remain unknown. Here, we systematically identified and characterized the expression patterns, evolutionary conservation, genomic features, and internal structures of circRNAs using RNase R-treated libraries from moso bamboo (Phyllostachys edulis) seedlings. Moreover, we investigated the biogenesis of circRNAs dependent on both cis- and trans-regulation. We explored the function of circRNAs, including their roles in regulating microRNA (miRNA)-related genes and modulating the alternative splicing of their linear counterparts. Importantly, we developed a customized degradome sequencing approach to detect miRNA-mediated cleavage of circRNAs. Finally, we presented a comprehensive view of the participation of circRNAs in the regulation of hormone metabolism upon treatment of bamboo seedlings with GA and NAA. Collectively, our study provides insights into the biogenesis, function, and miRNA-mediated degradation of circRNAs in moso bamboo.

The Function and Photoregulatory Mechanisms of Cryptochromes From Moso Bamboo (Phyllostachys edulis).

Light is one of the most important environmental factors affecting growth and geographic distribution of forestry plants. Moso bamboo is the largest temperate bamboo on earth and an important non-woody forestry species that serves not only important functions in the economy of rural areas but also carbon sequestration in the world. Due to its decades-long reproductive timing, the germplasm of moso bamboo cannot be readily improved by conventional breeding methods, arguing for a greater need to study the gene function and regulatory mechanisms of this species. We systematically studied the photoregulatory mechanisms of the moso bamboo (Phyllostachys edulis) cryptochrome 1, PheCRY1. Our results show that, similar to its Arabidopsis counterpart, the bamboo PheCRY1s are functionally restricted to the blue light inhibition of cell elongation without an apparent activity in promoting floral initiation. We demonstrate that PheCRY1s undergo light-dependent oligomerization that is inhibited by PheBIC1s, and light-dependent phosphorylation that is catalyzed by PhePPKs. We hypothesize that light-induced phosphorylation of PheCRY1s facilitate their degradation, which control availability of the PheCRY1 proteins and photosensitivity of bamboo plants. Our results demonstrate the evolutionary conservation of not only the function but also photoregulatory mechanism of PheCRY1 in this monocot forestry species.

Quantitative profiling of N6-methyladenosine at single-base resolution in stem-differentiating xylem of Populus trichocarpa using Nanopore direct RNA sequencing.

There are no comprehensive methods to identify N6-methyladenosine (m6A) at single-base resolution for every single transcript, which is necessary for the estimation of m6A abundance. We develop a new pipeline called Nanom6A for the identification and quantification of m6A modification at single-base resolution using Nanopore direct RNA sequencing based on an XGBoost model. We validate our method using methylated RNA immunoprecipitation sequencing (MeRIP-Seq) and m6A-sensitive RNA-endoribonuclease-facilitated sequencing (m6A-REF-seq), confirming high accuracy. Using this method, we provide a transcriptome-wide quantification of m6A modification in stem-differentiating xylem and reveal that different alternative polyadenylation (APA) usage shows a different ratio of m6A.

Classification of Rice Yield Using UAV-Based Hyperspectral Imagery and Lodging Feature.

High-yield rice cultivation is an effective way to address the increasing food demand worldwide. Correct classification of high-yield rice is a key step of breeding. However, manual measurements within breeding programs are time consuming and have high cost and low throughput, which limit the application in large-scale field phenotyping. In this study, we developed an accurate large-scale approach and presented the potential usage of hyperspectral data for rice yield measurement using the XGBoost algorithm to speed up the rice breeding process for many breeders. In total, 13 japonica rice lines in regional trials in northern China were divided into different categories according to the manual measurement of yield. Using an Unmanned Aerial Vehicle (UAV) platform equipped with a hyperspectral camera to capture images over multiple time series, a rice yield classification model based on the XGBoost algorithm was proposed. Four comparison experiments were carried out through the intraline test and the interline test considering lodging characteristics at the midmature stage or not. The result revealed that the degree of lodging in the midmature stage was an important feature affecting the classification accuracy of rice. Thus, we developed a low-cost, high-throughput phenotyping and nondestructive method by combining UAV-based hyperspectral measurements and machine learning for estimation of rice yield to improve rice breeding efficiency.

Associations between carabid beetles and fungi in the light of 200 years of published literature.

Describing and conserving ecological interactions woven into ecosystems is one of the great challenges of the 21st century. Here, we present a unique dataset compiling the biotic interactions between two ecologically and economically important taxa: ground beetles (Coleoptera: Carabidae) and fungi. The resulting dataset contains the carabid-fungus associations collected from 392 scientific publications, 129 countries, mostly from the Palearctic region, published over a period of 200 years. With an updated taxonomy to match the currently accepted nomenclature, 3,378 unique associations among 5,564 records were identified between 1,776 carabid and 676 fungal taxa. Ectoparasitic Laboulbeniales were the most frequent fungal group associated with carabids, especially with Trechinae. The proportion of entomopathogens was low. Three different formats of the data have been provided along with an interactive data digest platform for analytical purposes. Our database summarizes the current knowledge on biotic interactions between insects and fungi, while offering a valuable resource to test large-scale hypotheses on those interactions.

Large Scale Profiling of Protein Isoforms Using Label-Free Quantitative Proteomics Revealed the Regulation of Nonsense-Mediated Decay in Moso Bamboo (Phyllostachys edulis).

Moso bamboo is an important forest species with a variety of ecological, economic, and cultural values. However, the gene annotation information of moso bamboo is only based on the transcriptome sequencing, lacking the evidence of proteome. The lignification and fiber in moso bamboo leads to a difficulty in the extraction of protein using conventional methods, which seriously hinders research on the proteomics of moso bamboo. The purpose of this study is to establish efficient methods for extracting the total proteins from moso bamboo for following mass spectrometry-based quantitative proteome identification. Here, we have successfully established a set of efficient methods for extracting total proteins of moso bamboo followed by mass spectrometry-based label-free quantitative proteome identification, which further improved the protein annotation of moso bamboo genes. In this study, 10,376 predicted coding genes were confirmed by quantitative proteomics, accounting for 35.8% of all annotated protein-coding genes. Proteome analysis also revealed the protein-coding potential of 1015 predicted long noncoding RNA (lncRNA), accounting for 51.03% of annotated lncRNAs. Thus, mass spectrometry-based proteomics provides a reliable method for gene annotation. Especially, quantitative proteomics revealed the translation patterns of proteins in moso bamboo. In addition, the 3284 transcript isoforms from 2663 genes identified by Pacific BioSciences (PacBio) single-molecule real-time long-read isoform sequencing (Iso-Seq) was confirmed on the protein level by mass spectrometry. Furthermore, domain analysis of mass spectrometry-identified proteins encoded in the same genomic locus revealed variations in domain composition pointing towards a functional diversification of protein isoform. Finally, we found that part transcripts targeted by nonsense-mediated mRNA decay (NMD) could also be translated into proteins. In summary, proteomic analysis in this study improves the proteomics-assisted genome annotation of moso bamboo and is valuable to the large-scale research of functional genomics in moso bamboo. In summary, this study provided a theoretical basis and technical support for directional gene function analysis at the proteomics level in moso bamboo.

Analysis of Transcriptome and Epitranscriptome in Plants Using PacBio Iso-Seq and Nanopore-Based Direct RNA Sequencing.

Nanopore sequencing from Oxford Nanopore Technologies (ONT) and Pacific BioSciences (PacBio) single-molecule real-time (SMRT) long-read isoform sequencing (Iso-Seq) are revolutionizing the way transcriptomes are analyzed. These methods offer many advantages over most widely used high-throughput short-read RNA sequencing (RNA-Seq) approaches and allow a comprehensive analysis of transcriptomes in identifying full-length splice isoforms and several other post-transcriptional events. In addition, direct RNA-Seq provides valuable information about RNA modifications, which are lost during the PCR amplification step in other methods. Here, we present a comprehensive summary of important applications of these technologies in plants, including identification of complex alternative splicing (AS), full-length splice variants, fusion transcripts, and alternative polyadenylation (APA) events. Furthermore, we discuss the impact of the newly developed nanopore direct RNA-Seq in advancing epitranscriptome research in plants. Additionally, we summarize computational tools for identifying and quantifying full-length isoforms and other co/post-transcriptional events and discussed some of the limitations with these methods. Sequencing of transcriptomes using these new single-molecule long-read methods will unravel many aspects of transcriptome complexity in unprecedented ways as compared to previous short-read sequencing approaches. Analysis of plant transcriptomes with these new powerful methods that require minimum sample processing is likely to become the norm and is expected to uncover novel co/post-transcriptional gene regulatory mechanisms that control biological outcomes during plant development and in response to various stresses.

Identification and Characterization of the PEBP Family Genes in Moso Bamboo (Phyllostachys heterocycla).

Moso bamboo is one of the economically most important plants in China. Moso bamboo is a monocarpic perennial that exhibits poor and slow germination. Thus, the flowering often causes destruction of moso bamboo forestry. However, how control of flowering and seed germination are regulated in moso bamboo is largely unclear. In this study, we identified 5 members (PhFT1-5) of the phosphatidyl ethanolamine-binding proteins (PEBP) family from moso bamboo genome that regulate flowering, flower architecture and germination, and characterized the function of these PEBP family genes further in Arabidopsis. Phylogenetic analysis revealed that 3 (PhFT1, PhFT2 and PhFT3), 1 (PhFT4) and 1 (PhFT5) members belong to the TFL1-like clade, FT-like clade, and MFT-like clade, respectively. These PEBP family genes possess all structure necessary for PEBP gene function. The ectopic overexpression of PhFT4 and PhFT5 promotes flowering time in Arabidopsis, and that of PhFT1, PhFT2 and PhFT3 suppresses it. In addition, the overexpression of PhFT5 promotes seed germination rate. Interestingly, the overexpression of PhFT1 suppressed seed germination rate in Arabidopsis. The expression of PhFT1 and PhFT5 is significantly higher in seed than in tissues including leaf and shoot apical meristem, implying their function in seed germination. Taken together, our results suggested that the PEBP family genes play important roles as regulators of flowering and seed germination in moso bamboo and thereby are necessary for the sustainability of moso bamboo forest.