Efficient plant genome engineering using a probiotic sourced CRISPR-Cas9 system.

Among CRISPR-Cas genome editing systems, Streptococcus pyogenes Cas9 (SpCas9), sourced from a human pathogen, is the most widely used. Here, through in silico data mining, we have established an efficient plant genome engineering system using CRISPR-Cas9 from probiotic Lactobacillus rhamnosus. We have confirmed the predicted 5'-NGAAA-3' PAM via a bacterial PAM depletion assay and showcased its exceptional editing efficiency in rice, wheat, tomato, and Larix cells, surpassing LbCas12a, SpCas9-NG, and SpRY when targeting the identical sequences. In stable rice lines, LrCas9 facilitates multiplexed gene knockout through coding sequence editing and achieves gene knockdown via targeted promoter deletion, demonstrating high specificity. We have also developed LrCas9-derived cytosine and adenine base editors, expanding base editing capabilities. Finally, by harnessing LrCas9's A/T-rich PAM targeting preference, we have created efficient CRISPR interference and activation systems in plants. Together, our work establishes CRISPR-LrCas9 as an efficient and user-friendly genome engineering tool for diverse applications in crops and beyond.

Genome-wide chromatin accessibility landscape and dynamics of transcription factor networks during ovule and fiber development in cotton.

BACKGROUND: The development of cotton fiber is regulated by the orchestrated binding of regulatory proteins to cis-regulatory elements associated with developmental genes. The cis-trans regulatory dynamics occurred throughout the course of cotton fiber development are elusive. Here we generated genome-wide high-resolution DNase I hypersensitive sites (DHSs) maps to understand the regulatory mechanisms of cotton ovule and fiber development. RESULTS: We generated DNase I hypersensitive site (DHS) profiles from cotton ovules at 0 and 3 days post anthesis (DPA) and fibers at 8, 12, 15, and 18 DPA. We obtained a total of 1185 million reads and identified a total of 199,351 DHSs through ~ 30% unique mapping reads. It should be noted that more than half of DNase-seq reads mapped multiple genome locations and were not analyzed in order to achieve a high specificity of peak profile and to avoid bias from repetitive genomic regions. Distinct chromatin accessibilities were observed in the ovules (0 and 3 DPA) compared to the fiber elongation stages (8, 12, 15, and 18 DPA). Besides, the chromatin accessibility during ovules was particularly elevated in genomic regions enriched with transposable elements (TEs) and genes in TE-enriched regions were involved in ovule cell division. We analyzed cis-regulatory modules and revealed the influence of hormones on fiber development from the regulatory divergence of transcription factor (TF) motifs. Finally, we constructed a reliable regulatory network of TFs related to ovule and fiber development based on chromatin accessibility and gene co-expression network. From this network, we discovered a novel TF, WRKY46, which may shape fiber development by regulating the lignin content. CONCLUSIONS: Our results not only reveal the contribution of TEs in fiber development, but also predict and validate the TFs related to fiber development, which will benefit the research of cotton fiber molecular breeding.

The Landscapes of Gluten Regulatory Network in Elite Wheat Cultivars Contrasting in Gluten Strength.

Yangmai-13 (YM13) is a wheat cultivar with weak gluten fractions. In contrast, Zhenmai-168 (ZM168) is an elite wheat cultivar known for its strong gluten fractions and has been widely used in a number of breeding programs. However, the genetic mechanisms underlying the gluten signatures of ZM168 remain largely unclear. To address this, we combined RNA-seq and PacBio full-length sequencing technology to unveil the potential mechanisms of ZM168 grain quality. A total of 44,709 transcripts were identified in Y13N (YM13 treated with nitrogen) and 51,942 transcripts in Z168N (ZM168 treated with nitrogen), including 28,016 and 28,626 novel isoforms in Y13N and Z168N, respectively. Five hundred and eighty-four differential alternative splicing (AS) events and 491 long noncoding RNAs (lncRNAs) were discovered. Incorporating the sodium-dodecyl-sulfate (SDS) sedimentation volume (SSV) trait, both weighted gene coexpression network analysis (WGCNA) and multiscale embedded gene coexpression network analysis (MEGENA) were employed for network construction and prediction of key drivers. Fifteen new candidates have emerged in association with SSV, including 4 transcription factors (TFs) and 11 transcripts that partake in the post-translational modification pathway. The transcriptome atlas provides new perspectives on wheat grain quality and would be beneficial for developing promising strategies for breeding programs.

Bioinformatic Prediction of Bulked Oligonucleotide Probes for FISH Using Chorus2.

Fluorescence in situ hybridization (FISH) provides great conveniences for detection and visualization of specific genomic segments. Oligonucleotide (Oligo)-based FISH further broadened the applications in plant cytogenetics researches. High-specific single-copy oligo probes are essential for successful oligo-FISH experiments. Here, we introduce the bioinformatic pipeline to design genome-scaled single-copy oligos and filter repeat-related probes with Chorus2 software. Robust probes are accessible for both well-assembled genome and species without a reference genome based on this pipeline.

An efficient CRISPR-Cas12a promoter editing system for crop improvement.

Promoter editing represents an innovative approach to introduce quantitative trait variation (QTV) in crops. However, an efficient promoter editing system for QTV needs to be established. Here we develop a CRISPR-Cas12a promoter editing (CAPE) system that combines a promoter key-region estimating model and an efficient CRISPR-Cas12a-based multiplexed or singular editing system. CAPE is benchmarked in rice to produce QTV continuums for grain starch content and size by targeting OsGBSS1 and OsGS3, respectively. We then apply CAPE for promoter editing of OsD18, a gene encoding GA3ox in the gibberellin biosynthesis pathway. The resulting lines carry a QTV continuum of semidwarfism without significantly compromising grain measures. Field trials demonstrated that the OsD18 promoter editing lines have the same yield performance and antilodging phenotype as the Green Revolution OsSD1 mutants in different genetic backgrounds. Hence, promoter editing of OsD18 generates a quantitative Green Revolution trait. Together, we demonstrate a CAPE-based promoter editing and tuning pipeline for efficient production of useful QTV continuum in crops.

Chromosome painting reveals inter-chromosomal rearrangements and evolution of subgenome D of wheat.

Aegilops species represent the most important gene pool for breeding bread wheat (Triticum aestivum). Thus, understanding the genome evolution, including chromosomal structural rearrangements and syntenic relationships among Aegilops species or between Aegilops and wheat, is important for both basic genome research and practical breeding applications. In the present study, we attempted to develop subgenome D-specific fluorescence in situ hybridization (FISH) probes by selecting D-specific oligonucleotides based on the reference genome of Chinese Spring. The oligo-based chromosome painting probes consisted of approximately 26 000 oligos per chromosome and their specificity was confirmed in both diploid and polyploid species containing the D subgenome. Two previously reported translocations involving two D chromosomes have been confirmed in wheat varieties and their derived lines. We demonstrate that the oligo painting probes can be used not only to identify the translocations involving D subgenome chromosomes, but also to determine the precise positions of chromosomal breakpoints. Chromosome painting of 56 accessions of Ae. tauschii from different origins led us to identify two novel translocations: a reciprocal 3D-7D translocation in two accessions and a complex 4D-5D-7D translocation in one accession. Painting probes were also used to analyze chromosomes from more diverse Aegilops species. These probes produced FISH signals in four different genomes. Chromosome rearrangements were identified in Aegilops umbellulata, Aegilops markgrafii, and Aegilops uniaristata, thus providing syntenic information that will be valuable for the application of these wild species in wheat breeding.

The Methylation Inhibitor 5-Aza-2'-Deoxycytidine Induces Genome-Wide Hypomethylation in Rice.

DNA methylation is a conserved epigenetic modification which is vital for regulating gene expression and maintaining genome stability in both mammals and plants. Homozygous mutation of rice methyltransferase 1 (met1) gene can cause host death in rice, making it difficult to obtain plant material needed for hypomethylation research. To circumvent this challenge, the methylation inhibitor, 5-Aza-2'-deoxycytidine (AzaD), is used as a cytosine nucleoside analogue to reduce genome wide hypomethylation and is widely used in hypomethylation research. However, how AzaD affects plant methylation profiles at the genome scale is largely unknown. Here, we treated rice seedlings with AzaD and compared the AzaD treatment with osmet1-2 mutants, illustrating that there are similar CG hypomethylation and distribution throughout the whole genome. Along with global methylation loss class I transposable elements (TEs) which are farther from genes compared with class II TEs, were more significantly activated, and the RNA-directed DNA Methylation (RdDM) pathway was activated in specific genomic regions to compensate for severe CG loss. Overall, our results suggest that AzaD is an effective DNA methylation inhibitor that can influence genome wide methylation and cause a series of epigenetic variations.

De novo centromere formation in pericentromeric region of rice chromosome 8.

Neocentromeres develop when kinetochores assemble de novo at DNA loci that are not previously associated with CenH3 nucleosomes, and can rescue rearranged chromosomes that have lost a functional centromere. The molecular mechanisms associated with neocentromere formation in plants have been elusive. Here, we developed a Xian (indica) rice line with poor growth performance in the field due to approximately 272 kb deletion that spans centromeric DNA sequences, including the centromeric satellite repeat CentO, in the centromere of chromosome 8 (Cen8). The CENH3-binding domains were expanded downstream of the original CentO position in Cen8, which revealed a de novo centromere formation in rice. The neocentromere formation avoids chromosomal regions containing functional genes. Meanwhile, canonical histone H3 was replaced by CENH3 in the regions with low CENH3 levels, and the CenH3 nucleosomes in these regions became more periodic. In addition, we identified active genes in the deleted centromeric region, which are essential for chloroplast growth and development. In summary, our results provide valuable insights into neocentromere formation and show that functional genes exist in the centromeric regions of plant chromosomes.

Genome-wide analyses of PAM-relaxed Cas9 genome editors reveal substantial off-target effects by ABE8e in rice.

PAM-relaxed Cas9 nucleases, cytosine base editors and adenine base editors are promising tools for precise genome editing in plants. However, their genome-wide off-target effects are largely unexplored. Here, we conduct whole-genome sequencing (WGS) analyses of transgenic plants edited by xCas9, Cas9-NGv1, Cas9-NG, SpRY, nCas9-NG-PmCDA1, nSpRY-PmCDA1 and nSpRY-ABE8e in rice. Our results reveal that Cas9 nuclease and base editors, when coupled with the same guide RNA (gRNA), prefer distinct gRNA-dependent off-target sites. De novo generated gRNAs by SpRY editors lead to additional, but insubstantial, off-target mutations. Strikingly, ABE8e results in ~500 genome-wide A-to-G off-target mutations at TA motif sites per transgenic plant. ABE8e's preference for the TA motif is also observed at the target sites. Finally, we investigate the timeline and mechanism of somaclonal variation due to tissue culture, which chiefly contributes to the background mutations. This study provides a comprehensive understanding on the scale and mechanisms of off-target and background mutations occurring during PAM-relaxed genome editing in plants.

CRISPR-BETS: a base-editing design tool for generating stop codons.

Cytosine base editors (CBEs) can install a predefined stop codon at the target site, representing a more predictable and neater method for creating genetic knockouts without altering the genome size. Due to the enhanced predictability of the editing outcomes, it is also more efficient to obtain homozygous mutants in the first generation. With the recent advancement of CBEs on improved editing activity, purify and specificity in plants and animals, base editing has become a more appealing technology for generating knockouts. However, there is a lack of design tools that can aid the adoption of CBEs for achieving such a purpose, especially in plants. Here, we developed a user-friendly design tool named CRISPR-BETS (base editing to stop), which helps with guide RNA (gRNA) design for introducing stop codons in the protein-coding genes of interest. We demonstrated in rice and tomato that CRISPR-BETS is easy-to-use, and its generated gRNAs are highly specific and efficient for generating stop codons and obtaining homozygous knockout lines. While we tailored the tool for the plant research community, CRISPR-BETS can also serve non-plant species.

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.

Single Copy Oligonucleotide Fluorescence In Situ Hybridization Probe Design Platforms: Development, Application and Evaluation.

Oligonucleotides fluorescence in situ hybridization (Oligo-FISH) is an emerging technology and is an important tool in research areas such as detection of chromosome variation, identification of allopolyploid, and deciphering of three-dimensional (3D) genome structures. Based on the demand for highly efficient oligo probes for oligo-FISH experiments, increasing numbers of tools have been developed for probe design in recent years. Obsolete oligonucleotide design tools have been adapted for oligo-FISH probe design because of their similar considerations. With the development of DNA sequencing and large-scale synthesis, novel tools have been designed to increase the specificity of designed oligo probes and enable genome-scale oligo probe design, which has greatly improved the application of single copy oligo-FISH. Despite this, few studies have introduced the development of the oligo-FISH probe design tools and their application in FISH experiments systematically. Besides, a comprehensive comparison and evaluation is lacking for the available tools. In this review, we provide an overview of the oligo-FISH probe design process, summarize the development and application of the available tools, evaluate several state-of-art tools, and eventually provide guidance for single copy oligo-FISH probe design.

Improved plant cytosine base editors with high editing activity, purity, and specificity.

Cytosine base editors (CBEs) are great additions to the expanding genome editing toolbox. To improve C-to-T base editing in plants, we first compared seven cytidine deaminases in the BE3-like configuration in rice. We found A3A/Y130F-CBE_V01 resulted in the highest C-to-T base editing efficiency in both rice and Arabidopsis. Furthermore, we demonstrated this A3A/Y130F cytidine deaminase could be used to improve iSpyMacCas9-mediated C-to-T base editing at A-rich PAMs. To showcase its applications, we first applied A3A/Y130F-CBE_V01 for multiplexed editing to generate microRNA-resistant mRNA transcripts as well as pre-mature stop codons in multiple seed trait genes. In addition, we harnessed A3A/Y130F-CBE_V01 for efficient artificial evolution of novel ALS and EPSPS alleles which conferred herbicide resistance in rice. To further improve C-to-T base editing, multiple CBE_V02, CBE_V03 and CBE_V04 systems were developed and tested in rice protoplasts. The CBE_V04 systems were found to have improved editing activity and purity with focal recruitment of more uracil DNA glycosylase inhibitors (UGIs) by the engineered single guide RNA 2.0 scaffold. Finally, we used whole-genome sequencing (WGS) to compare six CBE_V01 systems and four CBE_V04 systems for genome-wide off-target effects in rice. Different levels of cytidine deaminase-dependent and sgRNA-independent off-target effects were indeed revealed by WGS among edited lines by these CBE systems. We also investigated genome-wide sgRNA-dependent off-target effects by different CBEs in rice. This comprehensive study compared 21 different CBE systems, and benchmarked PmCDA1-CBE_V04 and A3A/Y130F-CBE_V04 as next-generation plant CBEs with high editing efficiency, purity, and specificity.

Chorus2: design of genome-scale oligonucleotide-based probes for fluorescence in situ hybridization.

Oligonucleotide (oligo)-fluorescence in situ hybridization (FISH) has rapidly becoming the new generation of FISH technique in plant molecular cytogenetics research. Genome-scale identification of single-copy oligos is the foundation of successful oligo-FISH experiments. Here, we introduce Chorus2, a software that is developed specifically for oligo selection. We demonstrate that Chorus2 is highly effective to remove all repetitive elements in selection of single-copy oligos, which is critical for the development of successful FISH probes. Chorus2 is more effective than Chorus, the original version of the pipeline, and OligoMiner for repeat removal. Chorus2 allows to select oligos that are conserved among related species, which extends the usage of oligo-FISH probes among phylogenetically related plant species. We also implemented a new function in Chorus2 that allows development of FISH probes from plant species without an assembled genome. We anticipate that Chorus2 can be used in plants as well as in mammalian and other non-plant species. Chorus2 will broadly facilitate the design of FISH probes for various types of application in molecular cytogenetics research.

Evaluation and application of tools for the identification of known microRNAs in plants.

MicroRNAs (miRNAs), endogenous non-coding RNA regulators, post-transcriptionally inhibit the expression of their target genes. Several tools have been developed for predicting annotated known miRNAs, but there is no consensus about how to select the most suitable method for any given species. In this study, eight miRNA prediction tools (mirnovo, miRPlant, miRDeep-P2, miRExpress, miRkwood, miRDeep2, miR-PREFeR, and sRNAbench) were selected for evaluation. High-throughput small RNA sequencing data from four plant species (including C3 and C4 species, and both monocots and dicots, i.e., Arabidopsis thaliana, Oryza sativa, Triticum aestivum, and Zea mays) were used for the analysis. The sensitivity, accuracy, area under the curve, consistency, duration, and RAM usage of the known miRNA predictions were evaluated for each tool. The miRNA annotations were obtained using miRBase and sRNAanno. Algorithms, such as random forest, BLAST, and receiver operating characteristic curves, were used to evaluate accuracy. Of the tools evaluated, sRNAbench was found to be the most accurate, miRDeep-P2 was the most sensitive, miRDeep-P2 was the fastest, and miRkwood had the highest memory usage. Due to its large genome size, only three tools were able to successfully predict known miRNAs in wheat (Triticum aestivum). Our results enable us to recommend the tool best suited to a variety of researcher needs, which we hope will reduce confusion and enhance future work.

PAM-less plant genome editing using a CRISPR-SpRY toolbox.

The rapid development of the CRISPR-Cas9, -Cas12a and -Cas12b genome editing systems has greatly fuelled basic and translational plant research1-6. DNA targeting by these Cas nucleases is restricted by their preferred protospacer adjacent motifs (PAMs). The PAM requirement for the most popular Streptococcus pyogenes Cas9 (SpCas9) is NGG (N = A, T, C, G)7, limiting its targeting scope to GC-rich regions. Here, we demonstrate genome editing at relaxed PAM sites in rice (a monocot) and the Dahurian larch (a coniferous tree), using an engineered SpRY Cas9 variant8. Highly efficient targeted mutagenesis can be readily achieved by SpRY at relaxed PAM sites in the Dahurian larch protoplasts and in rice transgenic lines through non-homologous end joining (NHEJ). Furthermore, an SpRY-based cytosine base editor was developed and demonstrated by directed evolution of new herbicide resistant OsALS alleles in rice. Similarly, a highly active SpRY adenine base editor was developed based on ABE8e (ref. 9) and SpRY-ABE8e was able to target relaxed PAM sites in rice plants, achieving up to 79% editing efficiency with high product purity. Thus, the SpRY toolbox breaks a PAM restriction barrier in plant genome engineering by enabling DNA editing in a PAM-less fashion. Evidence was also provided for secondary off-target effects by de novo generated single guide RNAs (sgRNAs) due to SpRY-mediated transfer DNA self-editing, which calls for more sophisticated programmes for designing highly specific sgRNAs when implementing the SpRY genome editing toolbox.

Analysis of Off-Target Mutations in CRISPR-Edited Rice Plants Using Whole-Genome Sequencing.

The CRISPR/Cas systems have become the most widely used tool for genome editing in plants and beyond. However, CRISPR/Cas systems may cause unexpected off-target mutations due to sgRNA recognizing highly homologous DNA sequence elsewhere in the genome. Whole-genome sequencing (WGS) can be used to identify on- and off-target mutation. Here, we describe a pipeline of analyzing WGS data using a series of open source software for analysis of off-target mutations in CRISPR-edited rice plants. In this pipeline, the adapter is trimmed using SKEWER. Then, the cleaned reads are mapped to reference genome by applying BWA. To avoid mapping bias, the GATK is used to realign reads near indels (insertions and deletions) and recalibrate base quality controls. Whole-genome single nucleotide variations (SNVs) and indels are detected by LoFreq*, Mutect2, VarScan2, and Pindel. Last, SNVs and indels are compared with in silico off-target sites using Cas-OFFinder.

Computational approaches for effective CRISPR guide RNA design and evaluation.

The Clustered Regularly Interspaced Short Palindromic Repeat (CRISPR)/ CRISPR-associated (Cas) system has emerged as the main technology for gene editing. Successful editing by CRISPR requires an appropriate Cas protein and guide RNA. However, low cleavage efficiency and off-target effects hamper the development and application of CRISPR/Cas systems. To predict cleavage efficiency and specificity, numerous computational approaches have been developed for scoring guide RNAs. Most scores are empirical or trained by experimental datasets, and scores are implemented using various computational methods. Herein, we discuss these approaches, focusing mainly on the features or computational methods they utilise. Furthermore, we summarise these tools and give some suggestions for their usage. We also recommend three versatile web-based tools with user-friendly interfaces and preferable functions. The review provides a comprehensive and up-to-date overview of computational approaches for guide RNA design that could help users to select the optimal tools for their research.

Genome-wide Profiling of Histone Lysine Butyrylation Reveals its Role in the Positive Regulation of Gene Transcription in Rice.

BACKGROUND: Histone modifications play important roles in growth and development of rice (Oryza sativa L.). Lysine butyrylation (Kbu) with a four-carbon chain is a newly-discovered histone acylation modification in rice. MAIN BODY: In this study, we performed chromatin immunoprecipitation sequencing (ChIP-seq) analyses, the result showed that major enrichment of histone Kbu located in genebody regions of rice genome, especially in exons. The enrichment level of Kbu histone modification is positively correlated with gene expression. Furthermore, we compared Kbu with DNase-seq and other histone modifications in rice. We found that 60.06% Kub enriched region co-located with DHSs in intergenic regions. The similar profiles were detected among Kbu and several acetylation modifications such as H3K4ac, H3K9ac, and H3K23ac, indicating that Kbu modification is an active signal of transcription. Genes with both histone Kbu and one other acetylation also had significantly increased expression compared with genes with only one acetylation. Gene Ontology (GO) enrichment analysis revealed that these genes with histone Kbu can regulate multiple metabolic process in different rice varieties. CONCLUSION: Our study showed that the lysine butyrylation modificaiton may promote gene expression as histone acetylation and will provide resources for futher studies on histone Kbu and other epigenetic modifications in plants.