Analysis of the global transcriptome and miRNAome associated with seed dormancy during seed maturation in rice (Oryza sativa L. cv. Nipponbare).

BACKGROUND: Seed dormancy is a biological mechanism that prevents germination until favorable conditions for the subsequent generation of plants are encountered. Therefore, this mechanism must be effectively established during seed maturation. Studies investigating the transcriptome and miRNAome of rice embryos and endosperms at various maturation stages to evaluate seed dormancy are limited. This study aimed to compare the transcriptome and miRNAome of rice seeds during seed maturation. RESULTS: Oryza sativa L. cv. Nipponbare seeds were sampled for embryos and endosperms at three maturation stages: 30, 45, and 60 days after heading (DAH). The pre-harvest sprouting (PHS) assay was conducted to assess the level of dormancy in the seeds at each maturation stage. At 60 DAH, the PHS rate was significantly increased compared to those at 30 and 45 DAH, indicating that the dormancy is broken during the later maturation stage (45 DAH to 60 DAH). However, the largest number of differentially expressed genes (DEGs) and differentially expressed miRNAs (DEmiRs) were identified between 30 and 60 DAH in the embryo and endosperm, implying that the gradual changes in genes and miRNAs from 30 to 60 DAH may play a significant role in breaking seed dormancy. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analyses confirmed that DEGs related to plant hormones were most abundant in the embryo during 45 DAH to 60 DAH and 30 DAH to 60 DAH transitions. Alternatively, most of the DEGs in the endosperm were related to energy and abiotic stress. MapMan analysis and quantitative real-time polymerase chain reaction identified four newly profiled auxin-related genes (OsSAUR6/12/23/25) and one ethylene-related gene (OsERF087), which may be involved in seed dormancy during maturation. Additionally, miRNA target prediction (psRNATarget) and degradome dataset (TarDB) indicated a potential association between osa-miR531b and ethylene biosynthesis gene (OsACO4), along with osa-miR390-5p and the abscisic acid (ABA) exporter-related gene (OsMATE19) as factors involved in seed dormancy. CONCLUSIONS: Analysis of the transcriptome and miRNAome of rice embryos and endosperms during seed maturation provided new insights into seed dormancy, particularly its relationship with plant hormones such as ABA, auxin, and ethylene.

Distinct regulatory pathways contribute to dynamic CHH methylation patterns in transposable elements throughout Arabidopsis embryogenesis.

CHH methylation (mCHH) increases gradually during embryogenesis across dicotyledonous plants, indicating conserved mechanisms of targeting and conferral. Although it is suggested that methylation increase during embryogenesis enhances transposable element silencing, the detailed epigenetic pathways underlying this process remain unclear. In Arabidopsis, mCHH is regulated by both small RNA-dependent DNA methylation (RdDM) and RNA-independent Chromomethylase 2 (CMT2) pathways. Here, we conducted DNA methylome profiling at five stages of Arabidopsis embryogenesis, and classified mCHH regions into groups based on their dependency on different methylation pathways. Our analysis revealed that the gradual increase in mCHH in embryos coincided with the expansion of small RNA expression and regional mCHH spreading to nearby sites at numerous loci. We identified distinct methylation dynamics in different groups of mCHH targets, which vary according to transposon length, location, and cytosine frequency. Finally, we highlight the characteristics of transposable element loci that are targeted by different mCHH machinery, showing that short, heterochromatic TEs with lower mCHG levels are enriched in loci that switch from CMT2 regulation in leaves, to RdDM regulation during embryogenesis. Our findings highlight the interplay between the length, location, and cytosine frequency of transposons and the mCHH machinery in modulating mCHH dynamics during embryogenesis.

Contribution of RdDM to the ecotype-specific differential methylation on conserved as well as highly variable regions between Arabidopsis ecotypes.

BACKGROUND: Several studies showed genome-wide DNA methylation during Arabidopsis embryogenesis and germination. Although it has been known that the change of DNA methylation mainly occurs at CHH context mediated by small RNA-directed DNA methylation pathway during seed ripening and germination, the causality of the methylation difference exhibited in natural Arabidopsis ecotypes has not been thoroughly studied. RESULTS: In this study we compared DNA methylation difference using comparative pairwise multi-omics dynamics in Columbia-0 (Col) and Cape Verde Island (Cvi) ecotypes. Arabidopsis genome was divided into two regions, common regions in both ecotypes and Col-specific regions, depending on the reads mapping of whole genome bisulfite sequencing libraries from both ecotypes. Ecotype comparison was conducted within common regions and the levels of DNA methylation on common regions and Col-specific regions were also compared. we confirmed transcriptome were relatively dynamic in stage-wise whereas the DNA methylome and small RNAome were more ecotype-dependent. While the global CG methylation remains steady during maturation and germination, we found genic CG methylation differs the most between the two accessions. We also found that ecotype-specific differentially methylated regions (eDMR) are positively correlated with ecotype-specifically expressed 24-nt small RNA clusters. In addition, we discovered that Col-specific regions enriched with transposable elements (TEs) and structural variants that tend to become hypermethylated, and TEs in Col-specific regions were longer in size, more pericentromeric, and more hypermethylated than those in the common regions. Through the analysis of RdDM machinery mutants, we confirmed methylation on Col-specific region as well as on eDMRs in common region are contributed by RdDM pathway. Lastly, we demonstrated that highly variable sequences between ecotypes (HOT regions) were also affected by RdDM-mediated regulation. CONCLUSIONS: Through ecotype comparison, we revealed differences and similarities of their transcriptome, methylome and small RNAome both in global and local regions. We validated the contribution of RdDM causing differential methylation of common regions. Hypermethylated ecotype-specific regions contributed by RNA-directed DNA methylation pathway largely depend on the presence of TEs and copy-gain structural variations. These ecotype-specific regions are frequently associated with HOT regions, providing evolutionary insights into the epigenome dynamics within a species.

MiR1885 Regulates Disease Tolerance Genes in Brassica rapa during Early Infection with Plasmodiophora brassicae.

Clubroot caused by Plasmodiophora brassicae is a severe disease of cruciferous crops that decreases crop quality and productivity. Several clubroot resistance-related quantitative trait loci and candidate genes have been identified. However, the underlying regulatory mechanism, the interrelationships among genes, and how genes are regulated remain unexplored. MicroRNAs (miRNAs) are attracting attention as regulators of gene expression, including during biotic stress responses. The main objective of this study was to understand how miRNAs regulate clubroot resistance-related genes in P. brassicae-infected Brassica rapa. Two Brassica miRNAs, Bra-miR1885a and Bra-miR1885b, were revealed to target TIR-NBS genes. In non-infected plants, both miRNAs were expressed at low levels to maintain the balance between plant development and basal immunity. However, their expression levels increased in P. brassicae-infected plants. Both miRNAs down-regulated the expression of the TIR-NBS genes Bra019412 and Bra019410, which are located at a clubroot resistance-related quantitative trait locus. The Bra-miR1885-mediated down-regulation of both genes was detected for up to 15 days post-inoculation in the clubroot-resistant line CR Shinki and in the clubroot-susceptible line 94SK. A qRT-PCR analysis revealed Bra019412 expression was negatively regulated by miR1885. Both Bra019412 and Bra019410 were more highly expressed in CR Shinki than in 94SK; the same expression pattern was detected in multiple clubroot-resistant and clubroot-susceptible inbred lines. A 5' rapid amplification of cDNA ends analysis confirmed the cleavage of Bra019412 by Bra-miR1885b. Thus, miR1885s potentially regulate TIR-NBS gene expression during P. brassicae infections of B. rapa.

Identification of Genes and MicroRNAs Affecting Pre-harvest Sprouting in Rice (Oryza sativa L.) by Transcriptome and Small RNAome Analyses.

Pre-harvest sprouting (PHS) is one of the primary problems associated with seed dormancy in rice (Oryza sativa L.). It causes yield loss and reduces grain quality under unpredictable humid conditions at the ripening stage, thus affecting the economic value of the rice crop. To resolve this issue, understanding the molecular mechanism underlying seed dormancy in rice is important. Recent studies have shown that seed dormancy is affected by a large number of genes associated with plant hormone regulation. However, understanding regarding the effect of heat stress on seed dormancy and plant hormones is limited. This study compared the transcriptome and small RNAome of the seed embryo and endosperm of two contrasting japonica rice accessions, PHS susceptible (with low seed dormancy) and PHS resistant (with high seed dormancy), at three different maturation stages. We found that 9,068 genes and 35 microRNAs (miRNAs) were differentially expressed in the embryo, whereas 360 genes were differentially expressed in the endosperm. Furthermore, we identified and verified the candidate genes associated with seed dormancy and heat stress-related responses in rice using quantitative real-time PCR. We newly discovered eight hormone-related genes, four heat shock protein-related genes, and two miRNAs potentially involved in PHS. These findings provide a strong foundation for understanding the dynamics of transcriptome and small RNAome of hormone- and heat stress-related genes, which affect PHS during seed maturation.

Whole genome sequencing reveals a frameshift mutation and a large deletion in YY1AP1 in a girl with a panvascular artery disease.

BACKGROUND: Rare diseases are pathologies that affect less than 1 in 2000 people. They are difficult to diagnose due to their low frequency and their often highly heterogeneous symptoms. Rare diseases have in general a high impact on the quality of life and life expectancy of patients, which are in general children or young people. The advent of high-throughput sequencing techniques has improved diagnosis in several different areas, from pediatrics, achieving a diagnostic rate of 41% with whole genome sequencing (WGS) and 36% with whole exome sequencing, to neurology, achieving a diagnostic rate between 47 and 48.5% with WGS. This evidence has encouraged our group to pursue a molecular diagnosis using WGS for this and several other patients with rare diseases. RESULTS: We used whole genome sequencing to achieve a molecular diagnosis of a 7-year-old girl with a severe panvascular artery disease that remained for several years undiagnosed. We found a frameshift variant in one copy and a large deletion involving two exons in the other copy of a gene called YY1AP1. This gene is related to Grange syndrome, a recessive rare disease, whose symptoms include stenosis or occlusion of multiple arteries, congenital heart defects, brachydactyly, syndactyly, bone fragility, and learning disabilities. Bioinformatic analyses propose these mutations as the most likely cause of the disease, according to its frequency, in silico predictors, conservation analyses, and effect on the protein product. Additionally, we confirmed one mutation in each parent, supporting a compound heterozygous status in the child. CONCLUSIONS: In general, we think that this finding can contribute to the use of whole genome sequencing as a diagnosis tool of rare diseases, and in particular, it can enhance the set of known mutations associated with different diseases.

Transcriptomic analyses of rice (Oryza sativa) genes and non-coding RNAs under nitrogen starvation using multiple omics technologies.

BACKGROUND: Nitrogen (N) is a key macronutrient essential for plant growth, and its availability has a strong influence on crop development. The application of synthetic N fertilizers on crops has increased substantially in recent decades; however, the applied N is not fully utilized due to the low N use efficiency of crops. To overcome this limitation, it is important to understand the genome-wide responses and functions of key genes and potential regulatory factors in N metabolism. RESULTS: Here, we characterized changes in the rice (Oryza sativa) transcriptome, including genes, newly identified putative long non-coding RNAs (lncRNAs), and microRNAs (miRNAs) and their target mRNAs in response to N starvation using four different transcriptome approaches. Analysis of rice genes involved in N metabolism and/or transport using strand-specific RNA-Seq identified 2588 novel putative lncRNA encoding loci. Analysis of previously published RNA-Seq datasets revealed a group of N starvation-responsive lncRNAs showing differential expression under other abiotic stress conditions. Poly A-primed sequencing (2P-Seq) revealed alternatively polyadenylated isoforms of N starvation-responsive lncRNAs and provided precise 3' end information on the transcript models of these lncRNAs. Analysis of small RNA-Seq data identified N starvation-responsive miRNAs and down-regulation of miR169 family members, causing de-repression of NF-YA, as confirmed by strand-specific RNA-Seq and qRT-PCR. Moreover, we profiled the N starvation-responsive down-regulation of root-specific miRNA, osa-miR444a.4-3p, and Degradome sequencing confirmed MADS25 as a novel target gene. CONCLUSIONS: In this study, we used a combination of multiple RNA-Seq analyses to extensively profile the expression of genes, newly identified lncRNAs, and microRNAs in N-starved rice roots and shoots. Data generated in this study provide an in-depth understanding of the regulatory pathways modulated by N starvation-responsive miRNAs. The results of comprehensive, large-scale data analysis provide valuable information on multiple aspects of the rice transcriptome, which may be useful in understanding the responses of rice plants to changes in the N supply status of soil.

The Ancient Phosphatidylinositol 3-Kinase Signaling System Is a Master Regulator of Energy and Carbon Metabolism in Algae.

Algae undergo a complete metabolic transformation under stress by arresting cell growth, inducing autophagy and hyper-accumulating biofuel precursors such as triacylglycerols and starch. However, the regulatory mechanisms behind this stress-induced transformation are still unclear. Here, we use biochemical, mutational, and "omics" approaches to demonstrate that PI3K signaling mediates the homeostasis of energy molecules and influences carbon metabolism in algae. In Chlamydomonas reinhardtii, the inhibition and knockdown (KD) of algal class III PI3K led to significantly decreased cell growth, altered cell morphology, and higher lipid and starch contents. Lipid profiling of wild-type and PI3K KD lines showed significantly reduced membrane lipid breakdown under nitrogen starvation (-N) in the KD. RNA-seq and network analyses showed that under -N conditions, the KD line carried out lipogenesis rather than lipid hydrolysis by initiating de novo fatty acid biosynthesis, which was supported by tricarboxylic acid cycle down-regulation and via acetyl-CoA synthesis from glycolysis. Remarkably, autophagic responses did not have primacy over inositide signaling in algae, unlike in mammals and vascular plants. The mutant displayed a fundamental shift in intracellular energy flux, analogous to that in tumor cells. The high free fatty acid levels and reduced mitochondrial ATP generation led to decreased cell viability. These results indicate that the PI3K signal transduction pathway is the metabolic gatekeeper restraining biofuel yields, thus maintaining fitness and viability under stress in algae. This study demonstrates the existence of homeostasis between starch and lipid synthesis controlled by lipid signaling in algae and expands our understanding of such processes, with biotechnological and evolutionary implications.

Non-canonical targets destabilize microRNAs in human Argonautes.

Although much is known about microRNA (miRNA) biogenesis and the mechanism by which miRNAs regulate their targets, little is known about the regulation of miRNA stability. Mature miRNAs are stabilized by binding to Argonaute (Ago) proteins, the core components of the RNA-induced silencing complex (RISC). Recent studies suggest that interactions between miRNAs and their highly complementary target RNAs promote release of miRNAs from Ago proteins, and this in turn can lead to destabilization of miRNAs. However, the physiological triggers of miRNA destabilization with molecular mechanisms remain largely unknown. Here, using an in vitro system that consists of a minimal human Ago2-RISC in HEK293T cell lysates, we sought to understand how miRNAs are destabilized by their targets. Strikingly, we showed that miRNA destabilization is dramatically enhanced by an interaction with seedless, non-canonical targets. We then showed that this process entails not only unloading of miRNAs from Ago, but also 3΄ end destabilization of miRNAs occurred within Ago. Furthermore, our mutation analysis indicates that conformational changes in the hinge region of the Ago PAZ domain are likely to be the main driving force of the miRNA destabilization. Our collective results suggest that non-canonical targets may provide a stability control mechanism in the regulation of miRNAs in humans.

Genome-wide characterization of long intergenic non-coding RNAs (lincRNAs) provides new insight into viral diseases in honey bees Apis cerana and Apis mellifera.

BACKGROUND: Long non-coding RNAs (lncRNAs) are a class of RNAs that do not encode proteins. Recently, lncRNAs have gained special attention for their roles in various biological process and diseases. RESULTS: In an attempt to identify long intergenic non-coding RNAs (lincRNAs) and their possible involvement in honey bee development and diseases, we analyzed RNA-seq datasets generated from Asian honey bee (Apis cerana) and western honey bee (Apis mellifera). We identified 2470 lincRNAs with an average length of 1011 bp from A. cerana and 1514 lincRNAs with an average length of 790 bp in A. mellifera. Comparative analysis revealed that 5 % of the total lincRNAs derived from both species are unique in each species. Our comparative digital gene expression analysis revealed a high degree of tissue-specific expression among the seven major tissues of honey bee, different from mRNA expression patterns. A total of 863 (57 %) and 464 (18 %) lincRNAs showed tissue-dependent expression in A. mellifera and A. cerana, respectively, most preferentially in ovary and fat body tissues. Importantly, we identified 11 lincRNAs that are specifically regulated upon viral infection in honey bees, and 10 of them appear to play roles during infection with various viruses. CONCLUSIONS: This study provides the first comprehensive set of lincRNAs for honey bees and opens the door to discover lincRNAs associated with biological and hormone signaling pathways as well as various diseases of honey bee.

Higher biomass productivity of microalgae in an attached growth system, using wastewater.

Although most algae cultivation systems are operated in suspended culture, an attached growth system can offer several advantages over suspended systems. Algal cultivation becomes light-limited as the microalgal concentration increases in the suspended system; on the other hand, sunlight penetrates deeper and stronger in attached systems owing to the more transparent water. Such higher availability of sunlight makes it possible to operate a raceway pond deeper than usual, resulting in a higher areal productivity. The attached system achieved 2.8-times higher biomass productivity and total lipid productivity of 9.1 g m(-2) day(-1) and 1.9 g m(-2) day(-1), respectively, than the suspended system. Biomass productivity can be further increased by optimization of the culture conditions. Moreover, algal biomass harvesting and dewatering were made simpler and cheaper in attached systems, because mesh-type substrates with attached microalgae were easily removed from the culture and the remaining treated wastewater could be discharged directly. When the algal biomass was dewatered using natural sunlight, the palmitic acid (C16:0) content increased by 16% compared with the freeze-drying method. There was no great difference in other fatty acid composition. Therefore, the attached system for algal cultivation is a promising cultivation system for mass biodiesel production.

Lacibacter daechungensis sp. nov., isolated from deep freshwater of a reservoir.

A Gram-stain-negative, aerobic, motile by gliding, non-spore-forming, rod-shaped, orange-pigmented bacterium, designated strain H32-4(T), was isolated from 32 m deep water of Daechung reservoir in Daejeon, Republic of Korea. Based on the nucleotide sequence of the 16S rRNA gene, the closest neighbouring type strain was Lacibacter cauensis NJ-8(T) with which strain H32-4(T) shared 98.9% sequence similarity. The most abundant fatty acids in whole cells of strain H32-4(T) were C15 : 0 iso (40.6 %), C17 : 0 iso 3-OH (22.4%), summed feature 3 (C16:1ω7c and/or C16:1ω6c; 9.3%) and C15:0 (6.6 %). The predominant menaquinone was MK-7. The G+C content of the genomic DNA of strain H32-4(T) was 45.7 mol%. Thus, these combined genotypic and phenotypic data supported the conclusion that strain H32-4(T) represents a novel species of the genus Lacibacter, for which the name Lacibacter daechungensis sp. nov. is proposed. The type strain is H32-4(T) ( = KCTC 32395(T) = JCM 19172(T)).