Publications

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.

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.

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.

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.

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.

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.

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.

Alternative polyadenylation (APA) is an essential regulatory mechanism for gene expression. The next generation sequencing provides ample opportunity to precisely delineate APA sites genome-wide. Various methods for profiling transcriptome-wide poly(A) sites were developed. By comparing available methods, the ways for adding sequencing adaptors to fit with the Illumina sequencing platform are different. These methods have identified more than 50% genes that undergo APA in eukaryotes. However, due to the unbalanced PCR during library preparation, accurate quantification of poly(A) sites is still a challenge. Here, we describe an updated poly(A) tag sequencing method that incorporates unique molecular identifier (UMI) into the adaptor for removing quantification bias induced by PCR duplicates. Hence, quantification of poly(A) site usages can be achieved by counting UMIs. This protocol, quantifying poly(A) tag sequencing (QPAT-seq), can be finished in 1 day with reduced cost, and is particularly useful for application with a large number of samples.

The Avicennia marina is a mangrove species widely distributed throughout the tropical and subtropical intertidal wetlands. To adapt to adverse tidal waves and hypoxia environments, A. marina has evolved a sophisticated root system to better secure itself on the muddy soil with downward-grown anchor roots and upward-grown aerial roots, called pneumatophores. However, the process behind the development of a negative-gravitropic pneumatophore is not understood. Paraffin sections reveal anatomical differences among the shoots, anchor roots, and gas exchanging pneumatophores, clearly reflecting their functional diversions. The pneumatophore, in particular, contains abundant aerenchyma tissues and a thin cap structure at the tip. Transcriptomic analyses of both anchor roots and pneumatophores were performed to elucidate gene expression dynamics during the formation of pneumatophores. The results show that the plant hormone auxin regulates multiple different root initiations. The auxin related gene IAA19 plays a key role in pneumatophore development while the interaction of ethylene and abscisic acid is important for aerenchyma formation. Moreover, the molecular mechanisms behind pneumatophore anti-gravitropic growth may be regulated by the reduced strength of the statolith formation signaling pathway. These results shed light on the mechanistic understanding of pneumatophore formation in mangrove plants.

Heterochromatin is widespread in eukaryotic genomes and has diverse impacts depending on its genomic context. Previous studies have shown that a protein complex, the ASI1-AIPP1-EDM2 (AAE) complex, participates in polyadenylation regulation of several intronic heterochromatin-containing genes. However, the genome-wide functions of AAE are still unknown. Here, we show that the ASI1 and EDM2 mostly target the common genomic regions on a genome-wide level and preferentially interacts with genetic heterochromatin. Polyadenylation (poly(A) sequencing reveals that AAE complex has a substantial influence on poly(A) site usage of heterochromatin-containing genes, including not only intronic heterochromatin-containing genes but also the genes showing overlap with heterochromatin. Intriguingly, AAE is also involved in the alternative splicing regulation of a number of heterochromatin-overlapping genes, such as the disease resistance gene RPP4. We provided evidence that genic heterochromatin is indispensable for the recruitment of AAE in polyadenylation and splicing regulation. In addition to conferring RNA processing regulation at genic heterochromatin-containing genes, AAE also targets some transposable elements (TEs) outside of genes (including TEs sandwiched by genes and island TEs) for epigenetic silencing. Our results reveal new functions of AAE in RNA processing and epigenetic silencing, and thus represent important advances in epigenetic regulation.