The roles of DNA methylation on pH dependent i-motif (iM) formation in rice.

I-motifs (iMs) are four-stranded non-B DNA structures containing C-rich DNA sequences. The formation of iMs is sensitive to pH conditions and DNA methylation, although the extent of which is still unknown in both humans and plants. To investigate this, we here conducted iMab antibody-based immunoprecipitation and sequencing (iM-IP-seq) along with bisulfite sequencing using CK (original genomic DNA without methylation-related treatments) and hypermethylated or demethylated DNA at both pH 5.5 and 7.0 in rice, establishing a link between pH, DNA methylation and iM formation on a genome-wide scale. We found that iMs folded at pH 7.0 displayed higher methylation levels than those formed at pH 5.5. DNA demethylation and hypermethylation differently influenced iM formation at pH 7.0 and 5.5. Importantly, CG hypo-DMRs (differentially methylated regions) and CHH (H = A, C and T) hyper-DMRs alone or coordinated with CG/CHG hyper-DMRs may play determinant roles in the regulation of pH dependent iM formation. Thus, our study shows that the nature of DNA sequences alone or combined with their methylation status plays critical roles in determining pH-dependent formation of iMs. It therefore deepens the understanding of the pH and methylation dependent modulation of iM formation, which has important biological implications and practical applications.

Evaluation of Different SNP Analysis Software and Optimal Mining Process in Tree Species.

Single nucleotide polymorphism (SNP) is one of the most widely used molecular markers to help researchers understand the relationship between phenotypes and genotypes. SNP calling mainly consists of two steps, including read alignment and locus identification based on statistical models, and various software have been developed and applied in this issue. Meanwhile, in our study, very low agreement (<25%) was found among the prediction results generated by different software, which was much less consistent than expected. In order to obtain the optimal protocol of SNP mining in tree species, the algorithm principles of different alignment and SNP mining software were discussed in detail. And the prediction results were further validated based on in silico and experimental methods. In addition, hundreds of validated SNPs were provided along with some practical suggestions on program selection and accuracy improvement were provided, and we wish that these results could lay the foundation for the subsequent analysis of SNP mining.

Network and epigenetic characterization of subsets of genes specifically expressed in maize bundle sheath cells.

Bundle sheath (BS) cells exhibit dramatically structural differences and functional variations at physiological, biochemical and epigenetic levels as compared to mesophyll (M) cells in maize. The regulatory mechanisms controlling functional divergences between M and BS have been extensively investigated. However, BS cell-related regulatory networks are still not completely characterized. To address this, we herein conducted WGCNA-related co-expression assays using bulk M and BS cell RNA-seq data sets and identified a module containing 384 genes highly expressed in BS cells (including 20 hub TFs) instead of M cells. According to the hub TF centered regulatory network, we found that Dof22 and Dof30 might act as key regulators in the regulation of expression of BS-specific genes, and several MYB TFs exhibited a high collaboration with Dof TFs. By comparing the enrichment levels of histone modifications, we found that genes in the aforementioned module were more enriched with histone acetylation as compared to other BS-enriched DEGs with similar expression levels. Moreover, we found that a subset of genes functioning in photosynthesis, protein auto processing and enzymatic activities were significantly enriched with broad H3K4me3. Thus, our study provides evidence showing that regulatory network and histone modifications may play vital roles in the regulation of a subset of genes with important functions in BS cells.

Transcriptome-based gene regulatory network analyses of differential cold tolerance of two tobacco cultivars.

BACKGROUND: Cold is one of the main abiotic stresses that severely affect plant growth and development, and crop productivity as well. Transcriptional changes during cold stress have already been intensively studied in various plant species. However, the gene networks involved in the regulation of differential cold tolerance between tobacco varieties with contrasting cold resistance are quite limited. RESULTS: Here, we conducted multiple time-point transcriptomic analyses using Tai tobacco (TT, cold susceptibility) and Yan tobacco (YT, cold resistance) with contrasting cold responses. We identified similar DEGs in both cultivars after comparing with the corresponding control (without cold treatment), which were mainly involved in response to abiotic stimuli, metabolic processes, kinase activities. Through comparison of the two cultivars at each time point, in contrast to TT, YT had higher expression levels of the genes responsible for environmental stresses. By applying Weighted Gene Co-Expression Network Analysis (WGCNA), we identified two main modules: the pink module was similar while the brown module was distinct between the two cultivars. Moreover, we obtained 100 hub genes, including 11 important transcription factors (TFs) potentially involved in cold stress, 3 key TFs in the brown module and 8 key TFs in the pink module. More importantly, according to the genetic regulatory networks (GRNs) between TFs and other genes or TFs by using GENIE3, we identified 3 TFs (ABI3/VP1, ARR-B and WRKY) mainly functioning in differential cold responses between two cultivars, and 3 key TFs (GRAS, AP2-EREBP and C2H2) primarily involved in cold responses. CONCLUSION: Collectively, our study provides valuable resources for transcriptome- based gene network studies of cold responses in tobacco. It helps to reveal how key cold responsive TFs or other genes are regulated through network. It also helps to identify the potential key cold responsive genes for the genetic manipulation of tobacco cultivars with enhanced cold tolerance in the future.

The Harbinger transposon-derived gene PANDA epigenetically coordinates panicle number and grain size in rice.

Transposons significantly contribute to genome fractions in many plants. Although numerous transposon-related mutations have been identified, the evidence regarding transposon-derived genes regulating crop yield and other agronomic traits is very limited. In this study, we characterized a rice Harbinger transposon-derived gene called PANICLE NUMBER AND GRAIN SIZE (PANDA), which epigenetically coordinates panicle number and grain size. Mutation of PANDA caused reduced panicle number but increased grain size in rice, while transgenic plants overexpressing this gene showed the opposite phenotypic change. The PANDA-encoding protein can bind to the core polycomb repressive complex 2 (PRC2) components OsMSI1 and OsFIE2, and regulates the deposition of H3K27me3 in the target genes, thereby epigenetically repressing their expression. Among the target genes, both OsMADS55 and OsEMF1 were negative regulators of panicle number but positive regulators of grain size, partly explaining the involvement of PANDA in balancing panicle number and grain size. Moreover, moderate overexpression of PANDA driven by its own promoter in the indica rice cultivar can increase grain yield. Thus, our findings present a novel insight into the epigenetic control of rice yield traits by a Harbinger transposon-derived gene and provide its potential application for rice yield improvement.

Genome-wide characterization of i-motifs and their potential roles in the stability and evolution of transposable elements in rice.

I-motifs (iMs) are non-canonical DNA secondary structures that fold from cytosine (C)-rich genomic DNA regions termed putative i-motif forming sequences (PiMFSs). The structure of iMs is stabilized by hemiprotonated C-C base pairs, and their functions are now suspected in key cellular processes in human cells such as genome stability and regulation of gene transcription. In plants, their biological relevance is still largely unknown. Here, we characterized PiMFSs with high potential for i-motif formation in the rice genome by developing and applying a protocol hinging on an iMab antibody-based immunoprecipitation (IP) coupled with high-throughput sequencing (seq), consequently termed iM-IP-seq. We found that PiMFSs had intrinsic subgenomic distributions, cis-regulatory functions and an intricate relationship with DNA methylation. We indeed found that the coordination of PiMFSs with DNA methylation may affect dynamics of transposable elements (TEs) among different cultivated Oryza subpopulations or during evolution of wild rice species. Collectively, our study provides first and unique insights into the biology of iMs in plants, with potential applications in plant biotechnology for improving important agronomic rice traits.

Single-Cell Transcriptome and Network Analyses Unveil Key Transcription Factors Regulating Mesophyll Cell Development in Maize.

BACKGROUND: Maize mesophyll (M) cells play important roles in various biological processes such as photosynthesis II and secondary metabolism. Functional differentiation occurs during M-cell development, but the underlying mechanisms for regulating M-cell development are largely unknown. RESULTS: We conducted single-cell RNA sequencing (scRNA-seq) to profile transcripts in maize leaves. We then identified coregulated modules by analyzing the resulting pseudo-time-series data through gene regulatory network analyses. WRKY, ERF, NAC, MYB and Heat stress transcription factor (HSF) families were highly expressed in the early stage, whereas CONSTANS (CO)-like (COL) and ERF families were highly expressed in the late stage of M-cell development. Construction of regulatory networks revealed that these transcript factor (TF) families, especially HSF and COL, were the major players in the early and later stages of M-cell development, respectively. Integration of scRNA expression matrix with TF ChIP-seq and Hi-C further revealed regulatory interactions between these TFs and their targets. HSF1 and COL8 were primarily expressed in the leaf bases and tips, respectively, and their targets were validated with protoplast-based ChIP-qPCR, with the binding sites of HSF1 being experimentally confirmed. CONCLUSIONS: Our study provides evidence that several TF families, with the involvement of epigenetic regulation, play vital roles in the regulation of M-cell development in maize.

Epigenomic features of DNA G-quadruplexes and their roles in regulating rice gene transcription.

A DNA G-quadruplex (G4) is a non-canonical four-stranded nucleic acid structure involved in many biological processes in mammals. The current knowledge on plant DNA G4s, however, is limited; whether and how DNA G4s impact gene expression in plants is still largely unknown. Here, we applied a protocol referred to as BG4-DNA-IP-seq followed by a comprehensive characterization of DNA G4s in rice (Oryza sativa L.); we next integrated dG4s (experimentally detectable G4s) with existing omics data and found that dG4s exhibited differential DNA methylation between transposable element (TE) and non-TE genes. dG4 regions displayed genic-dependent enrichment of epigenomic signatures; finally, we showed that these sites displayed a positive association with expression of DNA G4-containing genes when located at promoters, and a negative association when located in the gene body, suggesting localization-dependent promotional/repressive roles of DNA G4s in regulating gene transcription. This study reveals interrelations between DNA G4s and epigenomic signatures, as well as implicates DNA G4s in modulating gene transcription in rice. Our study provides valuable resources for the functional characterization or bioengineering of some of key DNA G4s in rice.

Identification of recombination events in outbred species with next-generation sequencing data.

BACKGROUND: Meiotic recombination events include crossovers and non-crossovers or gene conversions. Although the rate of crossovers is often used for genetic mapping, the gene conversion events are not well studied especially in outbred species, which could produce distorted markers and thus affect the precision of genetic maps. RESULTS: We proposed a strategy for identifying gene conversion events in Populus with the next-generation sequencing (NGS) data from the two parents and their progeny in an F1 hybrid population. The strategy first involved phasing the heterozygous SNPs of the parents to obtain the parental haplotype blocks by NGS analytical tools, permitting to identify the parental gene conversion events with progeny genotypes. By incorporating available genetic linkage maps, longer haplotype blocks each corresponding to a chromosome can be created, not only allowing to detect crossover events but also possibly to locate a crossover in a small region. Our analysis revealed that gene conversions are more abundant than crossovers in Populus, with a higher probability to generate distorted markers in the regions involved than in the other regions on genome. The analytical procedures were implemented with Perl scripts as a freely available package, findGCO at https://github.com/tongchf/findGCO . CONCLUSIONS: The novel strategy and the new developed Perl package permit to identify gene conversion events with the next-generation sequencing technology in a hybrid population of outbred species. The new method revealed that in a genetic mapping population some distorted genetic markers are possibly due to the gene conversion events.

MVQTLCIM: composite interval mapping of multivariate traits in a hybrid F1 population of outbred species.

BACKGROUND: With the plummeting cost of the next-generation sequencing technologies, high-density genetic linkage maps could be constructed in a forest hybrid F1 population. However, based on such genetic maps, quantitative trait loci (QTL) mapping cannot be directly conducted with traditional statistical methods or tools because the linkage phase and segregation pattern of molecular markers are not always fixed as in inbred lines. RESULTS: We implemented the traditional composite interval mapping (CIM) method to multivariate trait data in forest trees and developed the corresponding software, mvqtlcim. Our method not only incorporated the various segregations and linkage phases of molecular markers, but also applied Takeuchi's information criterion (TIC) to discriminate the QTL segregation type among several possible alternatives. QTL mapping was performed in a hybrid F1 population of Populus deltoides and P. simonii, and 12 QTLs were detected for tree height over 6 time points. The software package allowed many options for parameters as well as parallel computing for permutation tests. The features of the software were demonstrated with the real data analysis and a large number of Monte Carlo simulations. CONCLUSIONS: We provided a powerful tool for QTL mapping of multiple or longitudinal traits in an outbred F1 population, in which the traditional software for QTL mapping cannot be used. This tool will facilitate studying of QTL mapping and thus will accelerate molecular breeding programs especially in forest trees. The tool package is freely available from https://github.com/tongchf /mvqtlcim.

Analysis of codon usage patterns in Ginkgo biloba reveals codon usage tendency from A/U-ending to G/C-ending.

As one of the most ancient tree species, the codon usage pattern analysis of Ginkgo biloba is a useful way to understand its evolutionary and genetic mechanisms. Several studies have been conducted on angiosperms, but seldom on gymnosperms. Based on RNA-Seq data of the G. biloba transcriptome, amount to 17,579 unigenes longer than 300 bp were selected and analyzed from 68,547 candidates. The codon usage pattern tended towards more frequently use of A/U-ending codons, which showed an obvious gradient progressing from gymnosperms to dicots to monocots. Meanwhile, analysis of high/low-expression unigenes revealed that high-expression unigenes tended to use G/C-ending codons together with more codon usage bias. Variation of unigenes with different functions suggested that unigenes involving in environment adaptation use G/C-ending codons more frequently with more usage bias, and these results were consistent with the conclusion that the formation of G. biloba codon usage bias was dominated by natural selection.

De novo SNP discovery and genetic linkage mapping in poplar using restriction site associated DNA and whole-genome sequencing technologies.

BACKGROUND: Restriction site associated DNA sequencing (RAD-seq), a next-generation sequencing technology, has greatly facilitated genetic linkage mapping studies in outbred species. RAD-seq is capable of discovering thousands of genetic markers for linkage mapping across many individuals, and can be applied in species with or without a reference genome. Although several analytical tools are available for RAD-seq data, alternative strategies are necessary for improving the marker quality and hence the genetic mapping accuracy. RESULTS: We demonstrate a strategy for constructing dense genetic linkage maps in hybrid forest trees by combining RAD-seq and whole-genome sequencing technologies. We performed RAD-seq of 150 progeny and whole-genome sequencing of the two parents in an F1 hybrid population of Populus deltoides × P. simonii. Two rough references were assembled from the whole-genome sequencing reads of the two parents separately. Based on the parental reference sequences, 3442 high-quality single nucleotide polymorphisms (SNPs) were identified that segregate in the ratio of 1:1. The maternal linkage map of P. deltoides was constructed with 2012 SNPs, containing 19 linkage groups and spanning 4067.16 cM of the genome with an average distance of 2.04 cM between adjacent markers, while the male map of P. simonii consisted of 1430 SNPs and the same number of linkage groups with a total length of 4356.04 cM and an average interval distance of 3.09 cM. Collinearity between the parental linkage maps and the reference genome of P. trichocarpa was also investigated. Compared with the result on the basis of the existing reference genome, our strategy identified more high-quality SNPs and generated parental linkage groups that nicely match the karyotype of Populus. CONCLUSIONS: The strategy of simultaneously using RAD and whole-genome sequencing technologies can be applied to constructing high-density genetic maps in forest trees regardless of whether a reference genome exists. The two parental linkage maps constructed here provide more accurate genetic resources for unraveling quantitative trait loci and accelerating molecular breeding programs, as well as for comparative genomics in Populus.