Environment-induced heritable variations are common in Arabidopsis thaliana.

Parental or ancestral environments can induce heritable phenotypic changes, but whether such environment-induced heritable changes are a common phenomenon remains unexplored. Here, we subject 14 genotypes of Arabidopsis thaliana to 10 different environmental treatments and observe phenotypic and genome-wide gene expression changes over four successive generations. We find that all treatments caused heritable phenotypic and gene expression changes, with a substantial proportion stably transmitted over all observed generations. Intriguingly, the susceptibility of a genotype to environmental inductions could be predicted based on the transposon abundance in the genome. Our study thus challenges the classic view that the environment only participates in the selection of heritable variation and suggests that the environment can play a significant role in generating of heritable variations.

Anchoring and catalytic insights into bilayer C4N3 material for lithium-selenium batteries: a first-principles study.

In the present work, a theoretical design for the viability of bilayer C4N3 (bi-C4N3) as a promising host material for Li-Se battery was conducted utilizing first-principles calculations. The AA- and AB-stacking configurations of bilayer C4N3 can effectively inhibit the shuttling of high-order polyselenides through the synergistic effect of physical confinement and strong Li-N bonds. Compared to conventional electrolytes, the AA- and AB-stacking bilayer C4N3 demonstrate enhanced adsorption capabilities for the polyselenides. The anchored structures of Se8 or Li2Sen (n = 1, 2, 4, 6, 8) molecules within the bilayer C4N3 exhibit high electrical conductivities, which are beneficial for enhancing the electrochemical performance. The catalytic effects of AA- and AB-stacking bilayer C4N3 were investigated by the reduction of Se8 and the energy barrier associated with the decomposition of Li2Se. The AA- and AB-stacking bilayer C4N3 can significantly decrease the activation barrier and promote the decomposition of Li2Se. The mean square displacement (MSD) curves reveal the pronounceably sluggish Li-ions diffusions in polyselenides within the AA- and AB-stacking bilayer C4N3, which in turn demonstrates the notable prospects in mitigating the shuttle effect.

Stress responses of plants through transcriptome plasticity by mRNA alternative polyadenylation.

The sessile nature of plants confines their responsiveness to changing environmental conditions. Gene expression regulation becomes a paramount mechanism for plants to adjust their physiological and morphological behaviors. Alternative polyadenylation (APA) is known for its capacity to augment transcriptome diversity and plasticity, thereby furnishing an additional set of tools for modulating gene expression. APA has also been demonstrated to exhibit intimate associations with plant stress responses. In this study, we review APA dynamic features and consequences in plants subjected to both biotic and abiotic stresses. These stresses include adverse environmental stresses, and pathogenic attacks, such as cadmium toxicity, high salt, hypoxia, oxidative stress, cold, heat shock, along with bacterial, fungal, and viral infections. We analyzed the overarching research framework employed to elucidate plant APA response and the alignment of polyadenylation site transitions with the modulation of gene expression levels within the ambit of each stress condition. We also proposed a general APA model where transacting factors, including poly(A) factors, epigenetic regulators, RNA m6A modification factors, and phase separation proteins, assume pivotal roles in APA related transcriptome plasticity during stress response in plants.

Inhibitor AN3661 reveals biological functions of Arabidopsis CLEAVAGE and POLYADENYLATION SPECIFICITY FACTOR 73.

Cleavage and polyadenylation specificity factor (CPSF) is a protein complex that plays an essential biochemical role in mRNA 3'-end formation, including poly(A) signal recognition and cleavage at the poly(A) site. However, its biological functions at the organismal level are mostly unknown in multicellular eukaryotes. The study of plant CPSF73 has been hampered by the lethality of Arabidopsis (Arabidopsis thaliana) homozygous mutants of AtCPSF73-I and AtCPSF73-II. Here, we used poly(A) tag sequencing to investigate the roles of AtCPSF73-I and AtCPSF73-II in Arabidopsis treated with AN3661, an antimalarial drug with specificity for parasite CPSF73 that is homologous to plant CPSF73. Direct seed germination on an AN3661-containing medium was lethal; however, 7-d-old seedlings treated with AN3661 survived. AN3661 targeted AtCPSF73-I and AtCPSF73-II, inhibiting growth through coordinating gene expression and poly(A) site choice. Functional enrichment analysis revealed that the accumulation of ethylene and auxin jointly inhibited primary root growth. AN3661 affected poly(A) signal recognition, resulted in lower U-rich signal usage, caused transcriptional readthrough, and increased the distal poly(A) site usage. Many microRNA targets were found in the 3' untranslated region lengthened transcripts; these miRNAs may indirectly regulate the expression of these targets. Overall, this work demonstrates that AtCPSF73 plays important part in co-transcriptional regulation, affecting growth, and development in Arabidopsis.

Effects of comprehensive landscape patterns on water quality and identification of key metrics thresholds causing its abrupt changes.

Comprehensive landscape patterns influence water quality with multiple factors, complex processes, and scale dependence. However, studies identifying landscape thresholds causing abrupt water quality changes and characterizing the contribution of topography to water quality are still limited. Exploring the impact mechanisms of natural geographical and landscape characteristics on spatial and seasonal water quality variations is conducive to watershed water resource protection and ecosystem restoration. Based on water quality monitoring data of Minjiahe River in the typical headwater area of the upstream Dan River in China from 2019 to 2021, we employed redundancy analysis, partial redundancy analysis, and nonparametric change-point analysis to analyze the relationship between stream water quality and multi-spatial scale comprehensive landscape patterns, to obtain the interactive and independent contributions of different landscape categories at multi-spatial scales on water quality, and to find the key landscape threshold leading to abrupt changes in water quality. Results showed that landscape configuration, landscape composition, and topographic factors collectively explain over 89.1% of water quality variation. Most seasonal variations in water quality were primarily caused by landscape configuration. The landscape composition was mainly responsible for the differences in water quality variations among spatial scales. The topographic factors made the least independent contribution and had a potential impact on overall water quality variation. In order to protect the water quality of streams, it is more reasonable to regulate the landscape at different scales. At the sub-catchment scale, interspersion and juxtaposition index (IJI) and landscape shape index (LSI) should be controlled below 82% and 22. At the 100 m riparian scale, farmland, urban land, IJI, and LSI should be controlled below 29%, 6.5%, 92%, and 26, respectively. Our results provide important guidance for optimizing landscape patterns and water conservation in the watershed.

SPR9 encodes a 60 S ribosomal protein that modulates panicle spreading and affects resistance to false smut in rice (Oryza sativa. L).

BACKGROUND: The architecture of inflorescence in crops is a key agronomic feature determining grain yield and thus has been a major target trait of cereal domestication. RESULTS: In this study, we show that a simple spreading panicle change in rice panicle shape, controlled by the Spreading Panicle 9 (SPR9) locus, also has a significant impact on the resistance to rice false smut (RFS). Meanwhile, we mapped a novel spr9 mutant gene between markers Indel5-18 and Indel5-22 encompassing a genomic region of 43-kb with six candidate genes. Through gene prediction and cDNA sequencing, we confirmed that LOC_Os05g38520 is the target gene in the spr9 mutant, which encodes 60 S ribosomal protein L36-2. Further analysis showed that the spr9 mutant is caused by a 1 bp deletion in the first exon that resulted in premature termination. Knockout experiments showed that the SPR9 gene is responsible for the spreading panicle phenotype of the spr9 mutant. Interestingly, the spr9 mutant was found to improve resistance to RFS without affecting major agronomic traits. Taken together, our results revealed that the spr9 allele has good application prospects in rice breeding for disease resistance and panicle improvement. CONCLUSIONS: We report the map-based cloning and functional characterization of SPR9, which encodes a 60 S ribosomal protein that regulates spreading panicles and affects the resistance to false smut in rice.

Microglia-containing cerebral organoids derived from induced pluripotent stem cells for the study of neurological diseases.

Microglia play an important role in neuroinflammation and neurodegeneration. Here, we report an approach for generating microglia-containing cerebral organoids derived from human pluripotent stem cells involving the supplementation of growth factors (FGF, EGF, heparin) and 10% CO2 culture conditions. Using this platform, Western Pacific Amyotrophic Lateral Sclerosis and Parkinsonism-Dementia Complex (ALS-PDC) cerebral organoids were generated from patient-derived induced pluripotent stem cells (iPSCs). These ALS-PDC-affected organoids had more reactive astrocytes and M1 microglia, and had fewer M2 microglia than their unaffected counterparts, leading to impaired microglia-mediated phagocytosis. RNA-seq analysis of ALS-PDC and control organoids indicated that the most significant changes were microglia- and astrocyte-related genes (IFITM1/2, TGF-ß, and GFAP). The most significantly downregulated pathway was type I interferon signaling. Interferon-gamma supplementation increased IFITM expression, enhanced microglia-mediated phagocytosis, and reduced beta-amyloid accumulation in ALS-PDC-affected network. The results demonstrated the feasibility of using microglia-containing organoids for the study of neurodegenerative diseases.

Differential alternative polyadenylation of homoeologous genes of allohexaploid wheat ABD subgenomes during drought stress response.

Because allohexaploid wheat genome contains ABD subgenomes, how the expression of homoeologous genes is coordinated remains largely unknown, particularly at the co-transcriptional level. Alternative polyadenylation (APA) is an important part of co-transcriptional regulation, which is crucial in developmental processes and stress responses. Drought stress is a major threat to the stable yield of wheat. Focusing on APA, we used poly(A) tag sequencing to track poly(A) site dynamics in wheat under drought stress. The results showed that drought stress led to extensive APA involving 37-47% of differentially expressed genes in wheat. Significant poly(A) site switching was found in stress-responsive genes. Interestingly, homoeologous genes exhibit unequal numbers of poly(A) sites, divergent APA patterns with tissue specificity and time-course dynamics, and distinct 3'-UTR length changes. Moreover, differentially expressed transcripts in leaves and roots used different poly(A) signals, the up- and downregulated isoforms had distinct preferences for non-canonical poly(A) sites. Genes that encode key polyadenylation factors showed differential expression patterns under drought stress. In summary, poly(A) signals and the changes in core poly(A) factors may widely affect the selection of poly(A) sites and gene expression levels during the response to drought stress, and divergent APA patterns among homoeologous genes add extensive plasticity to this responsive network. These results not only reveal the significant role of APA in drought stress response, but also provide a fresh perspective on how homoeologous genes contribute to adaptability through transcriptome diversity. In addition, this work provides information about the ends of transcripts for a better annotation of the wheat genome.

Coupling epigenetics and RNA polyadenylation: missing links.

Precise regulation of gene expression is crucial for plant survival. As a cotranscriptional regulatory mechanism, pre-mRNA polyadenylation is essential for fine-tuning gene expression. Polyadenylation can be alternatively projected at various sites of a transcript, which contributes to transcriptome diversity. Epigenetic modification is another mechanism of transcriptional control. Recent studies have uncovered crosstalk between cotranscriptional polyadenylation processes and both epigenomic and epitranscriptomic markers. Genetic analyses have demonstrated that DNA methylation, histone modifications, and epitranscriptomic modification are involved in regulating polyadenylation in plants. Here we summarize current understanding of the links between epigenetics and polyadenylation and their novel biological efficacy for plant development and environmental responses. Unresolved issues and future directions are discussed to shed light on the field.

Analysis of a rice blast resistance gene Pita-Fuhui2663 and development of selection marker.

Rice blast is a detrimental rice disease caused by the fungus Magnaporthe oryzae. Here, we identified a resistance gene from the rice cultivar Fuhui 2663 which is resistant to the rice blast isolate KJ201. Through isolated population analyses and sequencing approaches, the candidate gene was traced to chromosome 12. With the use of a map-based cloning strategy, the resistance gene was ultimately mapped to an 80-kb resistance locus region containing the Pita gene. Candidate gene prediction and cDNA sequencing indicated that the target resistance gene in Fuhui 2663 was allelic to Pita, thus being referred to as Pita-Fuhui2663 hereafter. Further analysis showed that the Fuhui 2663 protein had one amino acid change: Ala (A) residue 918 in Pita-Fuhui2663 was replaced by Ser (S) in Pita-S, leading to a significant change in the 3D structure of the Pita-S protein. CRISPR/Cas9 knockout experiments confirmed that Pita-Fuhui2663 is responsible for the resistance phenotype of Fuhui 2663. Importantly, Pita-Fuhui2663 did not affect the main agronomic traits of the variety compared to the Pita gene as verified by knockout experiments, indicative of potential applications of Pita-Fuhui2663 in broader breeding programs. Furthermore, a Pita-Fuhui2663-dCAPS molecular marker with good specificity and high efficiency was developed to facilitate rice breeding for resistance to this devastating disease.

A high-resolution single-molecule sequencing-based Arabidopsis transcriptome using novel methods of Iso-seq analysis.

BACKGROUND: Accurate and comprehensive annotation of transcript sequences is essential for transcript quantification and differential gene and transcript expression analysis. Single-molecule long-read sequencing technologies provide improved integrity of transcript structures including alternative splicing, and transcription start and polyadenylation sites. However, accuracy is significantly affected by sequencing errors, mRNA degradation, or incomplete cDNA synthesis. RESULTS: We present a new and comprehensive Arabidopsis thaliana Reference Transcript Dataset 3 (AtRTD3). AtRTD3 contains over 169,000 transcripts-twice that of the best current Arabidopsis transcriptome and including over 1500 novel genes. Seventy-eight percent of transcripts are from Iso-seq with accurately defined splice junctions and transcription start and end sites. We develop novel methods to determine splice junctions and transcription start and end sites accurately. Mismatch profiles around splice junctions provide a powerful feature to distinguish correct splice junctions and remove false splice junctions. Stratified approaches identify high-confidence transcription start and end sites and remove fragmentary transcripts due to degradation. AtRTD3 is a major improvement over existing transcriptomes as demonstrated by analysis of an Arabidopsis cold response RNA-seq time-series. AtRTD3 provides higher resolution of transcript expression profiling and identifies cold-induced differential transcription start and polyadenylation site usage. CONCLUSIONS: AtRTD3 is the most comprehensive Arabidopsis transcriptome currently. It improves the precision of differential gene and transcript expression, differential alternative splicing, and transcription start/end site usage analysis from RNA-seq data. The novel methods for identifying accurate splice junctions and transcription start/end sites are widely applicable and will improve single-molecule sequencing analysis from any species.

Plastid development of albino viviparous propagules in the woody mangrove species of Kandelia obovata.

The process of plastids developing into chloroplasts is critical for plants to survive. However, this process in woody plants is less understood. Kandelia obovata Sheue, Liu & Yong is a viviparous mangrove species; the seeds germinate on the maternal tree, and the hypocotyls continue to develop into mature propagules. We identified rare albino propagules through field observation among normal green and brown ones. Toward unveiling the propagule plastid development mechanism, albino propagule leaves only have etioplasts, low photosynthesis rates, and drastically reduced chlorophyll a/b and carotenoid contents, but with increased superoxide dismutase activities. To identify candidate genes controlling propagule plastid development, a genome-wide association study (GWAS) was performed between the albino and green propagules. Twenty-five significant single nucleotide polymorphisms (SNPs) were associated with albino propagule plastid development, the most significant SNPs being located on chromosomes 1 and 5. Significant differentially expressed genes were identified in porphyrin and chlorophyll metabolisms, carotenoid and flavonoid biosynthesis by combining transcriptome and GWAS data. In particular, KoDELLAs, encoding a transcription factor and KoCHS, encoding chalcone synthase, may be essential to regulate the albino propagules plastid development through weakened chlorophyll and flavonoid biosynthesis pathways while promoting chlorophyll degradation. Our results provide insights into genetic mechanisms regulating propagule plastid development in woody plants.

Signatures of mRNA Alternative Polyadenylation in Arabidopsis Leaf Development.

Alternative polyadenylation (APA) of pre-mRNA is an important co-transcriptional mechanism that modulates gene expression, leading to transcriptomic and functional diversities. The role of APA in Arabidopsis leaf development, however, remains elusive. We applied a poly(A)-tag sequencing (PAT-seq) technique to characterize APA-mediated regulation events in cotyledon and in five stages of true leaf development. Over 60% APA was identified in genes expressed in leaves, consistent with the results in previous publications. However, a reduced APA level was detected in younger leaves, reaching 44% in the 18th true leaf. Importantly, we also found that >70% of the poly(A) site usages were altered in the second true leaf relative to the cotyledon. Compared with the cotyledon, more genes in the second true leaf tended to use the distal site of 3'UTR, but this was not found in pairwise comparison among other true leaves. In addition, a significant APA gene was found to be decreased in a pairwise comparison among true leaves, including differentially expressed genes. The APA genes identified herein were associated with specific biological processes, including metabolic and cellular processes and response to stimuli and hormones. These results provide a new insight into the regulation of Arabidopsis leaf development through APA.

Divergence in the Regulation of the Salt Tolerant Response Between Arabidopsis thaliana and Its Halophytic Relative Eutrema salsugineum by mRNA Alternative Polyadenylation.

Salt tolerance is an important mechanism by which plants can adapt to a saline environment. To understand the process of salt tolerance, we performed global analyses of mRNA alternative polyadenylation (APA), an important regulatory mechanism during eukaryotic gene expression, in Arabidopsis thaliana and its halophytic relative Eutrema salsugineum with regard to their responses to salt stress. Analyses showed that while APA occurs commonly in both Arabidopsis and Eutrema, Eutrema possesses fewer APA genes than Arabidopsis (47% vs. 54%). However, the proportion of APA genes was significantly increased in Arabidopsis under salt stress but not in Eutrema. This indicated that Arabidopsis is more sensitive to salt stress and that Eutrema exhibits an innate response to such conditions. Both species utilized distal poly(A) sites under salt stress; however, only eight genes were found to overlap when their 3' untranslated region (UTR) lengthen genes were compared, thus revealing their distinct responses to salt stress. In Arabidopsis, genes that use distal poly(A) sites were enriched in response to salt stress. However, in Eutrema, the use of poly(A) sites was less affected and fewer genes were enriched. The transcripts with upregulated poly(A) sites in Arabidopsis showed enriched pathways in plant hormone signal transduction, starch and sucrose metabolism, and fatty acid elongation; in Eutrema, biosynthetic pathways (stilbenoid, diarylheptanoid, and gingerol) and metabolic pathways (arginine and proline) showed enrichment. APA was associated with 42% and 29% of the differentially expressed genes (DE genes) in Arabidopsis and Eutrema experiencing salt stress, respectively. Salt specific poly(A) sites and salt-inducible APA events were identified in both species; notably, some salt tolerance-related genes and transcription factor genes exhibited differential APA patterns, such as CIPK21 and LEA4-5. Our results suggest that adapted species exhibit more orderly response at the RNA maturation step under salt stress, while more salt-specific poly(A) sites were activated in Arabidopsis to cope with salinity conditions. Collectively, our findings not only highlight the importance of APA in the regulation of gene expression in response to salt stress, but also provide a new perspective on how salt-sensitive and salt-tolerant species perform differently under stress conditions through transcriptome diversity.

Phenotypic and Methylome Responses to Salt Stress in Arabidopsis thaliana Natural Accessions.

Salt stress threatens plant growth, development and crop yields, and has become a critical global environmental issue. Increasing evidence has suggested that the epigenetic mechanism such as DNA methylation can mediate plant response to salt stress through transcriptional regulation and transposable element (TE) silencing. However, studies exploring genome-wide methylation dynamics under salt stress remain limited, in particular, for studies on multiple genotypes. Here, we adopted four natural accessions of the model species Arabidopsis thaliana and investigated the phenotypic and genome-wide methylation responses to salt stress through whole-genome bisulfite sequencing (WGBS). We found that salt stress significantly changed plant phenotypes, including plant height, rosette diameter, fruit number, and aboveground biomass, and the change in biomass tended to depend on accessions. Methylation analysis revealed that genome-wide methylation patterns depended primarily on accessions, and salt stress caused significant methylation changes in ∼ 0.1% cytosines over the genomes. About 33.5% of these salt-induced differential methylated cytosines (DMCs) were located to transposable elements (TEs). These salt-induced DMCs were mainly hypermethylated and accession-specific. TEs annotated to have DMCs (DMC-TEs) across accessions were found mostly belonged to the superfamily of Gypsy, a type II transposon, indicating a convergent DMC dynamic on TEs across different genetic backgrounds. Moreover, 8.0% of salt-induced DMCs were located in gene bodies and their proximal regulatory regions. These DMCs were also accession-specific, and genes annotated to have DMCs (DMC-genes) appeared to be more accession-specific than DMC-TEs. Intriguingly, both accession-specific DMC-genes and DMC-genes shared by multiple accessions were enriched in similar functions, including methylation, gene silencing, chemical homeostasis, polysaccharide catabolic process, and pathways relating to shifts between vegetative growth and reproduction. These results indicate that, across different genetic backgrounds, methylation changes may have convergent functions in post-transcriptional, physiological, and phenotypic modulation under salt stress. These convergent methylation dynamics across accession may be autonomous from genetic variation or due to convergent genetic changes, which requires further exploration. Our study provides a more comprehensive picture of genome-wide methylation dynamics under salt stress, and highlights the importance of exploring stress response mechanisms from diverse genetic backgrounds.

Genome-wide identification and characterization of aquaporins in mangrove plant Kandelia obovata and its role in response to the intertidal environment.

Aquaporins (AQPs) play important roles in plant growth, development and tolerance to environmental stresses. To understand the role of AQPs in the mangrove plant Kandelia obovata, which has the ability to acquire water from seawater, we identified 34 AQPs in the K. obovata genome and analysed their structural features. Phylogenetic analysis revealed that KoAQPs are homologous to AQPs of Populus and Arabidopsis, which are evolutionarily conserved. The key amino acid residues were used to assess water-transport ability. Analysis of cis-acting elements in the promoters indicated that KoAQPs may be stress- and hormone-responsive. Subcellular localization of KoAQPs in yeast showed most KoAQPs function in the membrane system. That transgenic yeast with increased cell volume showed that some KoAQPs have significant water-transport activity, and the substrate sensitivity assay indicates that some KoAQPs can transport H2 O2 . The transcriptome data were used to analyze the expression patterns of KoAQPs in different tissues and developing fruits of K. obovata. In addition, real-time quantitative PCR analyses combined transcriptome data showed that KoAQPs have complex responses to environmental factors, including salinity, flooding and cold. Collectively, the transport of water and solutes by KoAQPs contributed to the adaptation of K. obovata to the coastal intertidal environment.

Uninterrupted embryonic growth leading to viviparous propagule formation in woody mangrove.

Vivipary is a rare sexual reproduction phenomenon where embryos germinate directly on the maternal plants. However, it is a common genetic event of woody mangroves in the Rhizophoraceae family. The ecological benefits of vivipary in mangroves include the nurturing of seedlings in harsh coastal and saline environments, but the genetic and molecular mechanisms of vivipary remain unclear. Here we investigate the viviparous embryo development and germination processes in mangrove Kandelia obovata by a transcriptomic approach. Many key biological pathways and functional genes were enriched in different tissues and stages, contributing to vivipary. Reduced production of abscisic acid set a non-dormant condition for the embryo to germinate directly. Genes involved in the metabolism of and response to other phytohormones (gibberellic acid, brassinosteroids, cytokinin, and auxin) are expressed precociously in the axis of non-vivipary stages, thus promoting the embryo to grow through the seed coat. Network analysis of these genes identified the central regulatory roles of LEC1 and FUS3, which maintain embryo identity in Arabidopsis. Moreover, photosynthesis related pathways were significantly up-regulated in viviparous embryos, and substance transporter genes were highly expressed in the seed coat, suggesting a partial self-provision and maternal nursing. We conclude that the viviparous phenomenon is a combinatorial result of precocious loss of dormancy and enhanced germination potential during viviparous seed development. These results shed light on the relationship between seed development and germination, where the continual growth of the embryo replaces a biphasic phenomenon until a mature propagule is established.

Global gene expression signatures in response to citrate-coated silver nanoparticles exposure.

Silver nanoparticles (AgNPs) are widely used in medical and commercial products for their unique antibacterial functions. However, the impact of AgNPs on human neural development is not well understood. To investigate the effect of AgNPs on human neural development, various doses of 20 nm citrate-coated AgNP (AgSC) were administered to human embryonic stem cell derived neural progenitors during the neuronal differentiation. Immunofluorescence staining with neuronal progenitor markers SOX2 (sex determining region Y-box 2) and Nestin (VI intermediate filament protein) showed that AgSC inhibited rosette formation, neuronal progenitor proliferation, and neurite outgrowth. Furthermore, AgSC promoted astrocyte activation and neuronal apoptosis. These adverse effects can be partially recovered with ascorbic acid. A genome-wide transcriptome analysis of both AgSC treated and untreated samples indicated that the most up-graduated genes were a group of Metallothionein (1F, 1E, 2A) proteins, a metal-binding protein that plays an essential role in metal homeostasis, heavy metal detoxification, and cellular anti-oxidative defence. The most significantly down-regulated genes were neuronal differentiation 6 (NEUROD6) and fork head box G1 (FOXG1). GO analyse indicated that the regulation of cholesterol biosynthetic process, neuron differentiation, synapse organization and pattern specification, oliogenesis, and neuronal apoptosis were the most impacted biological processes. KEGG pathway analyse showed that the most significantly impacted pathways were C5 isoprenoid, axon guidance, Notch, WNT, RAS-MAPK signalling pathways, lysosome, and apoptosis. Our data suggests that AgSCs interfered with metal homeostasis and cholesterol biosynthesis which induced oxidative stress, inhibited neurogenesis, axon guidance, and promoted apoptosis. Supplementation with ascorbic acid could act as an antioxidant to prevent AgSC-mediated neurotoxicity.

Intragenic heterochromatin-mediated alternative polyadenylation modulates miRNA and pollen development in rice.

Despite a much higher proportion of intragenic heterochromatin-containing genes in crop genomes, the importance of intragenic heterochromatin in crop development remains unclear. Intragenic heterochromatin can be recognised by a protein complex, ASI1-AIPP1-EDM2 (AAE) complex, to regulate alternative polyadenylation. Here, we investigated the impact of rice ASI1 on global poly(A) site usage through poly(A) sequencing and ASI1-dependent regulation on rice development. We found that OsASI1 is essential for rice pollen development and flowering. OsASI1 dysfunction has an important impact on global poly(A) site usage, which is closely related to heterochromatin marks. Intriguingly, OsASI1 interacts with the intronic heterochromatin of OsXRNL, a nuclear XRN family exonuclease gene involved in the processing of an miRNA precursor, to promote the processing of full-length OsXRNL and regulate miRNA abundance. We found that OsASI1-mediated regulation of pollen development partially depends on OsXRNL. Finally, we characterised the rice AAE complex and its involvement in alternative polyadenylation and pollen development. Our findings help to elucidate an epigenetic mechanism governing miRNA abundance and rice development, and provide a valuable resource for studying the epigenetic mechanisms of many important processes in crops.

QuantifyPoly(A): reshaping alternative polyadenylation landscapes of eukaryotes with weighted density peak clustering.

The dynamic choice of different polyadenylation sites in a gene is referred to as alternative polyadenylation, which functions in many important biological processes. Large-scale messenger RNA 3' end sequencing has revealed that cleavage sites for polyadenylation are presented with microheterogeneity. To date, the conventional determination of polyadenylation site clusters is subjective and arbitrary, leading to inaccurate annotations. Here, we present a weighted density peak clustering method, QuantifyPoly(A), to accurately quantify genome-wide polyadenylation choices. Applying QuantifyPoly(A) on published 3' end sequencing datasets from both animals and plants, their polyadenylation profiles are reshaped into myriads of novel polyadenylation site clusters. Most of these novel polyadenylation site clusters show significantly dynamic usage across different biological samples or associate with binding sites of trans-acting factors. Upstream sequences of these clusters are enriched with polyadenylation signals UGUA, UAAA and/or AAUAAA in a species-dependent manner. Polyadenylation site clusters also exhibit species specificity, while plants ones generally show higher microheterogeneity than that of animals. QuantifyPoly(A) is broadly applicable to any types of 3' end sequencing data and species for accurate quantification and construction of the complex and dynamic polyadenylation landscape and enables us to decode alternative polyadenylation events invisible to conventional methods at a much higher resolution.