Publications

2019

Yin J, Zhou M, Lin Z, Li QQ, Zhang Y-Y. Transgenerational effects benefit offspring across diverse environments: a meta-analysis in plants and animals.. Ecology letters. 2019;22(11):1976–1986. doi:10.1111/ele.13373

The adaptive value of transgenerational effects (the ancestor environmental effects on offspring) in changing environments has received much attention in recent years, but the related empirical evidence remains equivocal. Here, we conducted a meta-analysis summarising 139 experimental studies in plants and animals with 1170 effect sizes to investigate the generality of transgenerational effects across taxa, traits, and environmental contexts. It was found that transgenerational effects generally enhanced offspring performance in response to both stressful and benign conditions. The strongest effects are in annual plants and invertebrates, whereas vertebrates appear to benefit mostly under benign conditions, and perennial plants show hardly any transgenerational responses at all. These differences among taxonomic/life-history groups possibly reflect that vertebrates can avoid stressful conditions through their mobility, and longer-lived plants have alternative strategies. In addition to environmental contexts and taxonomic/life-history groups, transgenerational effects also varied among traits and developmental stages of ancestors and offspring, but the effects were similarly strong across three generations of offspring. By way of a more comprehensive data set and a different effect size, our results differ from those of a recent meta-analysis, suggesting that transgenerational effects are widespread, strong and persistent and can substantially impact the responses of plants and animals to changing environments.

2018

Hong L, Ye C, Lin J, Fu H, Wu X, Li QQ. Alternative polyadenylation is involved in auxin-based plant growth and development.. The Plant journal : for cell and molecular biology. 2018;93(2):246–258. doi:10.1111/tpj.13771

Auxin is widely involved in plant growth and development. However, the molecular mechanism on how auxin carries out this work is unclear. In particular, the effect of auxin on pre-mRNA post-transcriptional regulation is mostly unknown. By using a poly(A) tag (PAT) sequencing approach, mRNA alternative polyadenylation (APA) profiles after auxin treatment were revealed. We showed that hundreds of poly(A) site clusters (PACs) are affected by auxin at the transcriptome level, where auxin reduces PAC distribution in 5'-untranslated region (UTR), but increases in the 3'UTR. APA site usage frequencies of 42 genes were switched by auxin, suggesting that auxin affects the choice of poly(A) sites. Furthermore, poly(A) signal selection was altered after auxin treatment. For example, a mutant of poly(A) signal binding protein CPSF30 showed altered sensitivity to auxin treatment, indicating interactions between auxin and the poly(A) signal recognition machinery. We also found that auxin activity on lateral root development is likely mediated by altered expression of ARF7, ARF19 and IAA14 through poly(A) site switches. Our results shed light on the molecular mechanisms of auxin responses relative to its interactions with mRNA polyadenylation.

Hong L, Su W, Zhang Y, Ye C, Shen Y, Li QQ. Transcriptome profiling during mangrove viviparity in response to abscisic acid.. Scientific reports. 2018;8(1):770. doi:10.1038/s41598-018-19236-x

Mangrove plants adapt to coastal tidal mudflats with specially evolved viviparity seed development. However, very little is known about the genetic and molecular mechanisms of mangrove viviparity. Here, we tested a hypothesis that plant hormone abscisic acid (ABA) plays a significant role in precocious germination of viviparous Kandelia obovata seeds by exogenous applications. Through transcriptome analysis of ABA treated seeds, it was found that ABA repressed mangrove fruit growth and development, and there were thousands of genes differentially expressed. As a result, dynamics of the pathways were dramatically altered. In particular, "Plant hormone signal transduction" and "MAPK signaling pathway" were represented significantly. Among differentially expressed genes, some key genes of ABA signal transduction were induced, while ABA biosynthesis genes were repressed. Take ABI1 and ABI2, key negative regulators in ABA signal pathway, as examples, homologous alignment and a phylogenetic tree in various species showed that ABI1 and ABI2 are highly conserved among various species. The functional similarity of these genes was confirmed by transgenic work in Arabidopsis. Taken together, ABA inhibited mangrove viviparity, but mangroves developed a mechanism to prevent accidently increase of ABA in the harsh environment for maintaining viviparous reproductive strategy.

Ye C, Long Y, Ji G, Li QQ, Wu X. APAtrap: identification and quantification of alternative polyadenylation sites from RNA-seq data.. Bioinformatics (Oxford, England). 2018;34(11):1841–1849. doi:10.1093/bioinformatics/bty029

MOTIVATION: Alternative polyadenylation (APA) has been increasingly recognized as a crucial mechanism that contributes to transcriptome diversity and gene expression regulation. As RNA-seq has become a routine protocol for transcriptome analysis, it is of great interest to leverage such unprecedented collection of RNA-seq data by new computational methods to extract and quantify APA dynamics in these transcriptomes. However, research progress in this area has been relatively limited. Conventional methods rely on either transcript assembly to determine transcript 3' ends or annotated poly(A) sites. Moreover, they can neither identify more than two poly(A) sites in a gene nor detect dynamic APA site usage considering more than two poly(A) sites.

RESULTS: We developed an approach called APAtrap based on the mean squared error model to identify and quantify APA sites from RNA-seq data. APAtrap is capable of identifying novel 3' UTRs and 3' UTR extensions, which contributes to locating potential poly(A) sites in previously overlooked regions and improving genome annotations. APAtrap also aims to tally all potential poly(A) sites and detect genes with differential APA site usages between conditions. Extensive comparisons of APAtrap with two other latest methods, ChangePoint and DaPars, using various RNA-seq datasets from simulation studies, human and Arabidopsis demonstrate the efficacy and flexibility of APAtrap for any organisms with an annotated genome.

AVAILABILITY AND IMPLEMENTATION: Freely available for download at https://apatrap.sourceforge.io.

CONTACT: liqq@xmu.edu.cn or xhuister@xmu.edu.cn.

SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Repar N, Li H, Aguilar JS, Li QQ, Drobne D, Hong Y. Silver nanoparticles induce neurotoxicity in a human embryonic stem cell-derived neuron and astrocyte network.. Nanotoxicology. 2018;12(2):104–116. doi:10.1080/17435390.2018.1425497

Silver nanoparticles (AgNPs) are among the most extensively used nanoparticles and are found in a variety of products. This ubiquity leads to inevitable exposure to these particles in everyday life. However, the effects of AgNPs on neuron and astrocyte networks are still largely unknown. In this study, we used neurons and astrocytes derived from human embryonic stem cells as a cellular model to study the neurotoxicity that is induced by citrate-coated AgNPs (AgSCs). Immunostaining with the astrocyte and neuron markers, glial fibrillary acidic protein and microtubule-associated protein-2 (MAP2), respectively, showed that exposure to AgSCs at the concentration of 0.1 µg/mL increased the astrocyte/neuron ratio. In contrast, a higher concentration of AgSCs (5.0 µg/ml) significantly changed the morphology of astrocytes. These results suggest that astrocytes are sensitive to AgSC exposure and that low concentrations of AgSCs promote astrogenesis. Furthermore, our results showed that AgSCs reduced neurite outgrowth, decreased the expression of postsynaptic density protein 95 and synaptophysin, and induced neurodegeneration in a concentration-dependent manner. Our findings additionally suggest that the expression and phosphorylation status of MAP2 isoforms, as modulated by the activation of the Akt/glycogen synthase kinase-3/caspase-3 signaling pathway, may play an important role in AgSC-mediated neurotoxicity. We also found that AgNO3 exposure only slightly reduced neurite outgrowth and had little effect on MAP2 expression, suggesting that AgSCs and AgNO3 have different neuronal toxicity mechanisms. In addition, most of these effects were reduced when the cell culture was co-treated with AgSCs and the antioxidant ascorbic acid, which implies that oxidative stress is the major cause of AgSC-mediated astrocytic/neuronal toxicity and that antioxidants may have a neuroprotective effect.

2017

Irizarry KJL, Chan A, Kettle D, Kezian S, Ma D, Palacios L, Li QQ, Keeler CL, Drechsler Y. Bioinformatics Analysis of Chicken miRNAs Associated with Monocyte to Macrophage Differentiation and Subsequent IFNγ Stimulated Activation.. MicroRNA (Shariqah, United Arab Emirates). 2017;6(1):53–70. doi:10.2174/2211536605666161129122803

BACKGROUND: The goal of this project was to characterize the molecular and cellular roles of various gene targets regulated by miRNAs identified in differentiating and stimulating avian macrophages. Once a monocyte arrives to a site of infection, local signals induce a redistribution of resources into a macrophage phenotype. This may involve upregulating pathogen pattern recognizing receptors and increasing the efficiency of lysosomal biogenesis, while simultaneously recycling components involved in circulatory migration and leukocyte extravasation. a monocyte tooled with chemokine surface receptors and an internal cytoskeletal structure geared towards mobility may efficiently sense, react, and migrate toward a site of infection.

METHODS: Peripheral blood derived monocytes were purified and cultured from young chickens. RNA sequencing was performed on both peripheral blood monocytes during differentiation into macrophages and on mature macrophages following stimulation with interferon gamma. A set of microRNAs were identified and investigated using bioinformatics methods to ascertain their potential role in avian macrophage biology.

RESULTS: Among a number of miRNAs that are found to be expressed in avian macrophages, we focused on eight specific miRNAs (miR-1618, miR-1586, miR-1633, miR-1627, miR-1646, miR-1649, miR-1610, miR-1647) associated with macrophage differentiation and activation. Expression profiles of microRNAs were characterized during differentiation and activation. Candidate miRNA targets were implicated in processes including Wnt signaling, ubiquitination, PPAR mediated macrophage function, vesicle mediated cytokine trafficking, and WD40 domain protein functions.

CONCLUSION: A global theme for macrophage function that may be modulated by microRNAs is the comprehensive redistribution of the cell's protein repertoire. This redistribution involves two processes: 1) the degradation and recycling of unneeded cytoplasmic and membrane components and 2) the mobilization of newly synthesized cellular components via vesicular trafficking. Generally, it appears that macrophages need to closely regulate gene expression for differentiation to be able to activate successfully in response to a pathogen. This is a process in which miRNAs participate by affecting several pathways critical for both, differentiation and activation.

Pre-mRNA alternative splicing and alternative polyadenylation have been implicated to play important roles during eukaryotic gene expression. However, much remains unknown regarding the regulatory mechanisms and the interactions of these two processes in plants. Here we focus on an Arabidopsis gene OXT6 (Oxidative Tolerant-6) that has been demonstrated to encode two proteins through alternative splicing and alternative polyadenylation. Specifically, alternative polyadenylation at Intron-2 of OXT6 produces a transcript coding for AtCPSF30, an Arabidopsis ortholog of 30 kDa subunit of the Cleavage and Polyadenylation Specificity Factor. On the other hand, alternative splicing of Intron-2 generates a longer transcript encoding a protein named AtC30Y, a polypeptide including most part of AtCPSF30 and a YT521B domain. To investigate the expression outcome of OXT6 in plants, a set of mutations were constructed to alter the splicing and polyadenylation patterns of OXT6. Analysis of transgenic plants bearing these mutations by quantitative RT-PCR revealed a competition relationship between these two processes. Moreover, when both splice sites and poly(A) signals were mutated, polyadenylation became the preferred mode of OXT6 processing. These results demonstrate the interplay between alternative splicing and alternative polyadenylation, and it is their concerted actions that define a gene's expression outcome.

Lin J, Xu R, Wu X, Shen Y, Li QQ. Role of cleavage and polyadenylation specificity factor 100: anchoring poly(A) sites and modulating transcription termination.. The Plant journal : for cell and molecular biology. 2017;91(5):829–839. doi:10.1111/tpj.13611

CPSF100 is a core component of the cleavage and polyadenylation specificity factor (CPSF) complex for 3'-end formation of mRNA, but it still has no clear functional assignment. CPSF100 was reported to play a role in RNA silencing and promote flowering in Arabidopsis. However, the molecular mechanisms underlying these phenomena are not fully understood. Our genetics analyses indicate that plants with a hypomorphic mutant of CPSF100 (esp5) show defects in embryogenesis, reduced seed production or altered root morphology. To unravel this puzzle, we employed a poly(A) tag sequencing protocol and uncovered a different poly(A) profile in esp5. This transcriptome-wide analysis revealed alternative polyadenylation of thousands of genes, most of which result in transcriptional read-through in protein-coding genes. AtCPSF100 also affects poly(A) signal recognition on the far-upstream elements; in particular it prefers less U-rich sequences. Importantly, AtCPSF100 was found to exert its functions through the change of poly(A) sites on genes encoding binding proteins, such as nucleotide-binding, RNA-binding and poly(U)-binding proteins. In addition, through its interaction with RNA Polymerase II C-terminal domain (CTD) and affecting the expression level of CTD phosphatase-like 3 (CPL3), AtCPSF100 is shown to potentially ensure transcriptional termination by dephosphorylation of Ser2 on the CTD. These data suggest a key role for CPSF100 in locating poly(A) sites and affecting transcription termination.

Guo C, Spinelli M, Ye C, Li QQ, Liang C. Genome-Wide Comparative Analysis of Miniature Inverted Repeat Transposable Elements in 19 Arabidopsis thaliana Ecotype Accessions.. Scientific reports. 2017;7(1):2634. doi:10.1038/s41598-017-02855-1

Miniature inverted repeat transposable elements (MITEs) are prevalent in eukaryotic genomes. They are known to critically influence the process of genome evolution and play a role in gene regulation. As the first study concentrated in the transposition activities of MITEs among different ecotype accessions within a species, we conducted a genome-wide comparative analysis by characterizing and comparing MITEs in 19 Arabidopsis thaliana accessions. A total of 343485 MITE putative sequences, including canonical, diverse and partial ones, were delineated from all 19 accessions. Within the entire population of MITEs sequences, 80.7% of them were previously unclassified MITEs, demonstrating a different genomic distribution and functionality compared to the classified MITEs. The interactions between MITEs and homologous genes across 19 accessions provided a fine source for analyzing MITE transposition activities and their impacts on genome evolution. Moreover, a significant proportion of MITEs were found located in the last exon of genes besides the ordinary intron locality, thus potentially modifying the end of genes. Finally, analysis of the impact of MITEs on gene expression suggests that migrations of MITEs have no detectable effect on the expression level for host genes across accessions.

2016

Guo C, Spinelli M, Liu M, Li QQ, Liang C. A Genome-wide Study of "Non-3UTR" Polyadenylation Sites in Arabidopsis thaliana.. Scientific reports. 2016;6:28060. doi:10.1038/srep28060

Alternative polyadenylation has been recognized as a key contributor of gene expression regulation by generating different transcript isoforms with altered 3' ends. Although polyadenylation is well known for marking the end of a 3' UTR, an increasing number of studies have reported previously less-addressed polyadenylation events located in other parts of genes in many eukaryotic organisms. These other locations include 5' UTRs, introns and coding sequences (termed herein as non-3UTR), as well as antisense and intergenic polyadenlation. Focusing on the non-3UTR polyadenylation sites (n3PASs), we detected and characterized more than 11000 n3PAS clusters in the Arabidopsis genome using poly(A)-tag sequencing data (PAT-Seq). Further analyses suggested that the occurrence of these n3PASs were positively correlated with certain characteristics of their respective host genes, including the presence of spliced, diminutive or diverse beginning of 5' UTRs, number of introns and whether introns have extreme lengths. The interaction of the host genes with surrounding genetic elements, like a convergently overlapped gene and associated transposable element, may contribute to the generation of a n3PAS as well. Collectively, these results provide a better understanding of n3PASs, and offer some new insights of the underlying mechanisms for non-3UTR polyadenylation and its regulation in plants.