Corrigendum: Mining the Utricularia gibba genome for insulator-like elements for genetic engineering.

[This corrects the article DOI: 10.3389/fpls.2023.1279231.].

Protein phosphatase PP1 regulation of RNA polymerase II transcription termination and allelic exclusion of VSG genes in trypanosomes.

The genomes of Leishmania and trypanosomes are organized into polycistronic transcription units flanked by a modified DNA base J involved in promoting RNA polymerase II (Pol II) termination. We recently characterized a Leishmania complex containing a J-binding protein, PP1 protein phosphatase 1, and PP1 regulatory protein (PNUTS) that controls transcription termination potentially via dephosphorylation of Pol II by PP1. While T. brucei contains eight PP1 isoforms, none purified with the PNUTS complex, complicating the analysis of PP1 function in termination. We now demonstrate that the PP1-binding motif of TbPNUTS is required for function in termination in vivo and that TbPP1-1 modulates Pol II termination in T. brucei and dephosphorylation of the large subunit of Pol II. PP1-1 knock-down results in increased cellular levels of phosphorylated RPB1 accompanied by readthrough transcription and aberrant transcription of the chromosome by Pol II, including Pol I transcribed loci that are typically silent, such as telomeric VSG expression sites involved in antigenic variation. These results provide important insights into the mechanism underlying Pol II transcription termination in primitive eukaryotes that rely on polycistronic transcription and maintain allelic exclusion of VSG genes.

scifi-ATAC-seq: massive-scale single-cell chromatin accessibility sequencing using combinatorial fluidic indexing.

Single-cell ATAC-seq has emerged as a powerful approach for revealing candidate cis-regulatory elements genome-wide at cell-type resolution. However, current single-cell methods suffer from limited throughput and high costs. Here, we present a novel technique called scifi-ATAC-seq, single-cell combinatorial fluidic indexing ATAC-sequencing, which combines a barcoded Tn5 pre-indexing step with droplet-based single-cell ATAC-seq using the 10X Genomics platform. With scifi-ATAC-seq, up to 200,000 nuclei across multiple samples can be indexed in a single emulsion reaction, representing an approximately 20-fold increase in throughput compared to the standard 10X Genomics workflow.

Genome-wide DNA methylation dynamics following recent polyploidy in the allotetraploid Tragopogon miscellus (Asteraceae).

Polyploidy is an important evolutionary force, yet epigenetic mechanisms, such as DNA methylation, that regulate genome-wide expression of duplicated genes remain largely unknown. Here, we use Tragopogon (Asteraceae) as a model system to discover patterns and temporal dynamics of DNA methylation in recently formed polyploids. The naturally occurring allotetraploid Tragopogon miscellus formed in the last 95-100 yr from parental diploids Tragopogon dubius and T. pratensis. We profiled the DNA methylomes of these three species using whole-genome bisulfite sequencing. Genome-wide methylation levels in T. miscellus were intermediate between its diploid parents. However, nonadditive CG and CHG methylation occurred in transposable elements (TEs), with variation among TE types. Most differentially methylated regions (DMRs) showed parental legacy, but some novel DMRs were detected in the polyploid. Differentially methylated genes (DMGs) were also identified and characterized. This study provides the first assessment of both overall and locus-specific patterns of DNA methylation in a recent natural allopolyploid and shows that novel methylation variants can be generated rapidly after polyploid formation. Together, these results demonstrate that mechanisms to regulate duplicate gene expression may arise soon after allopolyploid formation and that these mechanisms vary among genes.

Investigating the cis-Regulatory Basis of C3 and C4 Photosynthesis in Grasses at Single-Cell Resolution.

While considerable knowledge exists about the enzymes pivotal for C4 photosynthesis, much less is known about cis-regulation important for specifying their expression in distinct cell types. Here, we use single-cell-indexed ATAC-seq to identify cell-type-specific accessible chromatin regions (ACRs) associated with C4 enzymes for five different grass species. This study spans four C4 species, covering three distinct photosynthetic subtypes: Zea mays and Sorghum bicolor (NADP-ME), Panicum miliaceum (NAD-ME), Urochloa fusca (PEPCK), along with the C3 outgroup Oryza sativa. We studied the cis-regulatory landscape of enzymes essential across all C4 species and those unique to C4 subtypes, measuring cell-type-specific biases for C4 enzymes using chromatin accessibility data. Integrating these data with phylogenetic trees revealed diverse co-option of gene family members between species, showcasing the various paths of C4 evolution. We mapped cell-type-specific ACRs surrounding each C4 gene, discovering that, on average, these genes have two to three cell-type-specific ACRs located beyond the core promoter. These results provide a detailed genomic map of potential regulatory sequences operating on key C4 genes. Examining the evolutionary history of these cell-type-specific ACRs revealed a spectrum of conserved and novel ACRs, even among closely related species, indicating ongoing evolution of cis-regulation at these C4 loci. This work illuminates the dynamic and complex nature of CRE evolution in C4 photosynthesis, particularly highlighting the intricate cis-regulatory evolution of key loci. Our findings offer a valuable resource for future investigations, potentially aiding in the optimization of C3 crop performance under changing climatic conditions.

The unusual predominance of maintenance DNA methylation in Spirodela polyrhiza.

Duckweeds are among the fastest reproducing plants, able to clonally divide at exponential rates. However, the genetic and epigenetic impact of clonality on plant genomes is poorly understood. 5-methylcytosine (5mC) is a modified base often described as necessary for the proper regulation of certain genes and transposons and for the maintenance of genome integrity in plants. However, the extent of this dogma is limited by the current phylogenetic sampling of land plant species diversity. Here we analyzed DNA methylomes, small RNAs, mRNA-seq, and H3K9me2 histone modification for Spirodela polyrhiza. S. polyrhiza has lost highly conserved genes involved in de novo methylation of DNA at sites often associated with repetitive DNA, and within genes, however, symmetrical DNA methylation and heterochromatin are maintained during cell division at certain transposons and repeats. Consequently, small RNAs that normally guide methylation to silence repetitive DNA like retrotransposons are diminished. Despite the loss of a highly conserved methylation pathway, and the reduction of small RNAs that normally target repetitive DNA, transposons have not proliferated in the genome, perhaps due in part to the rapid, clonal growth lifestyle of duckweeds.

Evolution of plant cell-type-specific cis -regulatory elements.

Cis -regulatory elements (CREs) are critical in regulating gene expression, and yet our understanding of CRE evolution remains a challenge. Here, we constructed a comprehensive single-cell atlas of chromatin accessibility in Oryza sativa , integrating data from 104,029 nuclei representing 128 discrete cell states across nine distinct organs. We used comparative genomics to compare cell-type resolved chromatin accessibility between O. sativa and 57,552 nuclei from four additional grass species ( Zea mays, Sorghum bicolor, Panicum miliaceum , and Urochloa fusca ). Accessible chromatin regions (ACRs) had different levels of conservation depending on the degree of cell-type specificity. We found a complex relationship between ACRs with conserved noncoding sequences, cell-type specificity, conservation, and tissue-specific switching. Additionally, we found that epidermal ACRs were less conserved compared to other cell types, potentially indicating that more rapid regulatory evolution has occurred in the L1 epidermal layer of these species. Finally, we identified and characterized a conserved subset of ACRs that overlapped the repressive histone modification H3K27me3, implicating them as potentially critical silencer CREs maintained by evolution. Collectively, this comparative genomics approach highlights the dynamics of cell-type-specific CRE evolution in plants.

Massive-scale single-cell chromatin accessibility sequencing using combinatorial fluidic indexing.

Single-cell ATAC-seq has emerged as a powerful approach for revealing candidate cis-regulatory elements genome-wide at cell-type resolution. However, current single-cell methods suffer from limited throughput and high costs. Here, we present a novel technique called single-cell combinatorial fluidic indexing ATAC-sequencing ("scifi-ATAC-seq"), which combines a barcoded Tn5 pre-indexing step with droplet-based single-cell ATAC-seq using a widely commercialized microfluidics platform (10X Genomics). With scifi-ATAC-seq, up to 200,000 nuclei across multiple samples in a single emulsion reaction can be indexed, representing a ~20-fold increase in throughput compared to the standard 10X Genomics workflow.

Mining the Utricularia gibba genome for insulator-like elements for genetic engineering.

Introduction: Gene expression is often controlled via cis-regulatory elements (CREs) that modulate the production of transcripts. For multi-gene genetic engineering and synthetic biology, precise control of transcription is crucial, both to insulate the transgenes from unwanted native regulation and to prevent readthrough or cross-regulation of transgenes within a multi-gene cassette. To prevent this activity, insulator-like elements, more properly referred to as transcriptional blockers, could be inserted to separate the transgenes so that they are independently regulated. However, only a few validated insulator-like elements are available for plants, and they tend to be larger than ideal. Methods: To identify additional potential insulator-like sequences, we conducted a genome-wide analysis of Utricularia gibba (humped bladderwort), one of the smallest known plant genomes, with genes that are naturally close together. The 10 best insulator-like candidates were evaluated in vivo for insulator-like activity. Results: We identified a total of 4,656 intergenic regions with expression profiles suggesting insulator-like activity. Comparisons of these regions across 45 other plant species (representing Monocots, Asterids, and Rosids) show low levels of syntenic conservation of these regions. Genome-wide analysis of unmethylated regions (UMRs) indicates ~87% of the targeted regions are unmethylated; however, interpretation of this is complicated because U. gibba has remarkably low levels of methylation across the genome, so that large UMRs frequently extend over multiple genes and intergenic spaces. We also could not identify any conserved motifs among our selected intergenic regions or shared with existing insulator-like elements for plants. Despite this lack of conservation, however, testing of 10 selected intergenic regions for insulator-like activity found two elements on par with a previously published element (EXOB) while being significantly smaller. Discussion: Given the small number of insulator-like elements currently available for plants, our results make a significant addition to available tools. The high hit rate (2 out of 10) also implies that more useful sequences are likely present in our selected intergenic regions; additional validation work will be required to identify which will be most useful for plant genetic engineering.

Engineering custom morpho- and chemotypes of Populus for sustainable production of biofuels, bioproducts, and biomaterials.

Humans have been modifying plant traits for thousands of years, first through selection (i.e., domestication) then modern breeding, and in the last 30 years, through biotechnology. These modifications have resulted in increased yield, more efficient agronomic practices, and enhanced quality traits. Precision knowledge of gene regulation and function through high-resolution single-cell omics technologies, coupled with the ability to engineer plant genomes at the DNA sequence, chromatin accessibility, and gene expression levels, can enable engineering of complex and complementary traits at the biosystem level. Populus spp., the primary genetic model system for woody perennials, are among the fastest growing trees in temperate zones and are important for both carbon sequestration and global carbon cycling. Ample genomic and transcriptomic resources for poplar are available including emerging single-cell omics datasets. To expand use of poplar outside of valorization of woody biomass, chassis with novel morphotypes in which stem branching and tree height are modified can be fabricated thereby leading to trees with altered leaf to wood ratios. These morphotypes can then be engineered into customized chemotypes that produce high value biofuels, bioproducts, and biomaterials not only in specific organs but also in a cell-type-specific manner. For example, the recent discovery of triterpene production in poplar leaf trichomes can be exploited using cell-type specific regulatory sequences to synthesize high value terpenes such as the jet fuel precursor bisabolene specifically in the trichomes. By spatially and temporally controlling expression, not only can pools of abundant precursors be exploited but engineered molecules can be sequestered in discrete cell structures in the leaf. The structural diversity of the hemicellulose xylan is a barrier to fully utilizing lignocellulose in biomaterial production and by leveraging cell-type-specific omics data, cell wall composition can be modified in a tailored and targeted specific manner to generate poplar wood with novel chemical features that are amenable for processing or advanced manufacturing. Precision engineering poplar as a multi-purpose sustainable feedstock highlights how genome engineering can be used to re-imagine a crop species.

Protein Phosphatase PP1 Regulation of Pol II Phosphorylation is Linked to Transcription Termination and Allelic Exclusion of VSG Genes and TERRA in Trypanosomes.

The genomes of Leishmania and trypanosomes are organized into polycistronic transcription units flanked by a modified DNA base J involved in promoting RNA polymerase II (Pol II) termination. We recently characterized a Leishmania complex containing a J-binding protein, PP1 protein phosphatase 1, and PP1 regulatory protein (PNUTS) that controls transcription termination potentially via dephosphorylation of Pol II by PP1. While T. brucei contains eight PP1 isoforms, none purified with the PNUTS complex, suggesting a unique PP1-independent mechanism of termination. We now demonstrate that the PP1-binding motif of TbPNUTS is required for function in termination in vivo and that TbPP1-1 modulates Pol II termination in T. brucei involving dephosphorylation of the C-terminal domain of the large subunit of Pol II. PP1-1 knock-down results in increased cellular levels of phosphorylated large subunit of Pol II accompanied by readthrough transcription and pervasive transcription of the entire genome by Pol II, including Pol I transcribed loci that are typically silent, such as telomeric VSG expression sites involved in antigenic variation and production of TERRA RNA. These results provide important insights into the mechanism underlying Pol II transcription termination in primitive eukaryotes that rely on polycistronic transcription and maintain allelic exclusion of VSG genes.

Enhancers associated with unstable RNAs are rare in plants.

Unstable transcripts have emerged as markers of active enhancers in vertebrates and shown to be involved in many cellular processes and medical disorders. However, their prevalence and role in plants is largely unexplored. Here, we comprehensively captured all actively initiating ("nascent") transcripts across diverse crops and other plants using capped small (cs)RNA-seq. We discovered that unstable transcripts are rare, unlike in vertebrates, and often originate from promoters. Additionally, many "distal" elements in plants initiate tissue-specific stable transcripts and are likely bone fide promoters of yet-unannotated genes or non-coding RNAs, cautioning against using genome annotations to infer "enhancers" or transcript stability. To investigate enhancer function, we integrated STARR-seq data. We found that annotated promoters, and other regions that initiate stable transcripts rather than unstable transcripts, function as stronger enhancers in plants. Our findings underscore the blurred line between promoters and enhancers and suggest that cis-regulatory elements encompass diverse structures and mechanisms in eukaryotes.

Haplotype-resolved genome assembly of Populus tremula × P. alba reveals aspen-specific megabase satellite DNA.

Populus species play a foundational role in diverse ecosystems and are important renewable feedstocks for bioenergy and bioproducts. Hybrid aspen Populus tremula × P. alba INRA 717-1B4 is a widely used transformation model in tree functional genomics and biotechnology research. As an outcrossing interspecific hybrid, its genome is riddled with sequence polymorphisms which present a challenge for sequence-sensitive analyses. Here we report a telomere-to-telomere genome for this hybrid aspen with two chromosome-scale, haplotype-resolved assemblies. We performed a comprehensive analysis of the repetitive landscape and identified both tandem repeat array-based and array-less centromeres. Unexpectedly, the most abundant satellite repeats in both haplotypes lie outside of the centromeres, consist of a 147 bp monomer PtaM147, frequently span >1 megabases, and form heterochromatic knobs. PtaM147 repeats are detected exclusively in aspens (section Populus) but PtaM147-like sequences occur in LTR-retrotransposons of closely related species, suggesting their origin from the retrotransposons. The genomic resource generated for this transformation model genotype has greatly improved the design and analysis of genome editing experiments that are highly sensitive to sequence polymorphisms. The work should motivate future hypothesis-driven research to probe into the function of the abundant and aspen-specific PtaM147 satellite DNA.

Leveraging Single-Cell Populations to Uncover the Genetic Basis of Complex Traits.

The ease and throughput of single-cell genomics have steadily improved, and its current trajectory suggests that surveying single-cell populations will become routine. We discuss the merger of quantitative genetics with single-cell genomics and emphasize how this synergizes with advantages intrinsic to plants. Single-cell population genomics provides increased detection resolution when mapping variants that control molecular traits, including gene expression or chromatin accessibility. Additionally, single-cell population genomics reveals the cell types in which variants act and, when combined with organism-level phenotype measurements, unveils which cellular contexts impact higher-order traits. Emerging technologies, notably multiomics, can facilitate the measurement of both genetic changes and genomic traits in single cells, enabling single-cell genetic experiments. The implementation of single-cell genetics will advance the investigation of the genetic architecture of complex molecular traits and provide new experimental paradigms to study eukaryotic genetics.

The role of Dmnt1 during spermatogenesis of the insect Oncopeltus fasciatus.

BACKGROUND: The function of DNA methyltransferase genes of insects is a puzzle, because an association between gene expression and methylation is not universal for insects. If the genes normally involved in cytosine methylation are not influencing gene expression, what might be their role? We previously demonstrated that gametogenesis of Oncopeltus fasciatus is interrupted at meiosis following knockdown of DNA methyltransferase 1 (Dnmt1) and this is unrelated to changes in levels of cytosine methylation. Here, using transcriptomics, we tested the hypothesis that Dmnt1 is a part of the meiotic gene pathway. Testes, which almost exclusively contain gametes at varying stages of development, were sampled at 7 days and 14 days following knockdown of Dmnt1 using RNAi. RESULTS: Using microscopy, we found actively dividing spermatocysts were reduced at both timepoints. However, as with other studies, we saw Dnmt1 knockdown resulted in condensed nuclei after mitosis-meiosis transition, and then cellular arrest. We found limited support for a functional role for Dnmt1 in our predicted cell cycle and meiotic pathways. An examination of a priori Gene Ontology terms showed no enrichment for meiosis. We then used the full data set to reveal further candidate pathways influenced by Dnmt1 for further hypotheses. Very few genes were differentially expressed at 7 days, but nearly half of all transcribed genes were differentially expressed at 14 days. We found no strong candidate pathways for how Dnmt1 knockdown was achieving its effect through Gene Ontology term overrepresentation analysis. CONCLUSIONS: We, therefore, suggest that Dmnt1 plays a role in chromosome dynamics based on our observations of condensed nuclei and cellular arrest with no specific molecular pathways disrupted.

Exploring plant cis-regulatory elements at single-cell resolution: overcoming biological and computational challenges to advance plant research.

Cis-regulatory elements (CREs) are important sequences for gene expression and for plant biological processes such as development, evolution, domestication, and stress response. However, studying CREs in plant genomes has been challenging. The totipotent nature of plant cells, coupled with the inability to maintain plant cell types in culture and the inherent technical challenges posed by the cell wall has limited our understanding of how plant cell types acquire and maintain their identities and respond to the environment via CRE usage. Advances in single-cell epigenomics have revolutionized the field of identifying cell-type-specific CREs. These new technologies have the potential to significantly advance our understanding of plant CRE biology, and shed light on how the regulatory genome gives rise to diverse plant phenomena. However, there are significant biological and computational challenges associated with analyzing single-cell epigenomic datasets. In this review, we discuss the historical and foundational underpinnings of plant single-cell research, challenges, and common pitfalls in the analysis of plant single-cell epigenomic data, and highlight biological challenges unique to plants. Additionally, we discuss how the application of single-cell epigenomic data in various contexts stands to transform our understanding of the importance of CREs in plant genomes.

Knockout of protein phosphatase 1 in Leishmania major reveals its role during RNA polymerase II transcription termination.

The genomes of kinetoplastids are organized into polycistronic transcription units that are flanked by a modified DNA base (base J, beta-D-glucosyl-hydroxymethyluracil). Previous work established a role of base J in promoting RNA polymerase II (Pol II) termination in Leishmania major and Trypanosoma brucei. We recently identified a PJW/PP1 complex in Leishmania containing a J-binding protein (JBP3), PP1 phosphatase 1, PP1 interactive-regulatory protein (PNUTS) and Wdr82. Analyses suggested the complex regulates transcription termination by recruitment to termination sites via JBP3-base J interactions and dephosphorylation of proteins, including Pol II, by PP1. However, we never addressed the role of PP1, the sole catalytic component, in Pol II transcription termination. We now demonstrate that deletion of the PP1 component of the PJW/PP1 complex in L. major, PP1-8e, leads to readthrough transcription at the 3'-end of polycistronic gene arrays. We show PP1-8e has in vitro phosphatase activity that is lost upon mutation of a key catalytic residue and associates with PNUTS via the conserved RVxF motif. Additionally, purified PJW complex with associated PP1-8e, but not complex lacking PP1-8e, led to dephosphorylation of Pol II, suggesting a direct role of PNUTS/PP1 holoenzymes in regulating transcription termination via dephosphorylating Pol II in the nucleus.

DNMT1 mutant ants develop normally but have disrupted oogenesis.

Although DNA methylation is an important gene regulatory mechanism in mammals, its function in arthropods remains poorly understood. Studies in eusocial insects have argued for its role in caste development by regulating gene expression and splicing. However, such findings are not always consistent across studies, and have therefore remained controversial. Here we use CRISPR/Cas9 to mutate the maintenance DNA methyltransferase DNMT1 in the clonal raider ant, Ooceraea biroi. Mutants have greatly reduced DNA methylation, but no obvious developmental phenotypes, demonstrating that, unlike mammals, ants can undergo normal development without DNMT1 or DNA methylation. Additionally, we find no evidence of DNA methylation regulating caste development. However, mutants are sterile, whereas in wild-type ants, DNMT1 is localized to the ovaries and maternally provisioned into nascent oocytes. This supports the idea that DNMT1 plays a crucial but unknown role in the insect germline.

Smar2C2: A Simple and Efficient Protocol for the Identification of Transcription Start Sites.

Promoters and the noncoding sequences that drive their function are fundamental aspects of genes that are critical to their regulation. The transcription preinitiation complex binds and assembles on promoters where it facilitates transcription. The transcription start site (TSS) is located downstream of the promoter sequence and is defined as the location in the genome where polymerase begins transcribing DNA into RNA. Knowing the location of TSSs is useful for annotation of genes, identification of non-coding sequences important to gene regulation, detection of alternative TSSs, and understanding of 5' UTR content. Several existing techniques make it possible to accurately identify TSSs, but are often difficult to perform experimentally, require large amounts of input RNA, or are unable to identify a large number of TSSs from a single sample. Many of these protocols take advantage of template switching reverse transcriptases (TSRTs), which reliably place an adaptor at the 5' end of a first strand synthesis of cDNA. Here, we introduce a protocol that exploits TSRT activity combined with rolling circle amplification to identify TSSs with several unique advantages over existing methods. Sequence adaptors are placed on the 5' and 3' end of the full-length cDNA copy of a transcript. A splint compatible with those adaptors is then used to circularize the full-length cDNA. Linear DNA containing concatemers of the cDNA are generated using rolling circle amplification, and a sequencing library is formed by fragmenting the concatemers. This protocol is straightforward to execute, requiring limited bench time with relatively stable reagents. Using extremely low amounts of RNA input, this protocol produces large numbers of accurate, deduplicated TSSs genome wide. © 2023 The Authors. Current Protocols published by Wiley Periodicals LLC. Basic Protocol 1: Splint generation Basic Protocol 2: RNA extraction Basic Protocol 3: cDNA synthesis Basic Protocol 4: cDNA circularization and amplification Basic Protocol 5: Library generation.

Sequencing Illumina libraries at high accuracy on the ONT MinION using R2C2.

High-throughput short-read sequencing has taken on a central role in research and diagnostics. Hundreds of different assays take advantage of Illumina short-read sequencers, the predominant short-read sequencing technology available today. Although other short-read sequencing technologies exist, the ubiquity of Illumina sequencers in sequencing core facilities and the high capital costs of these technologies have limited their adoption. Among a new generation of sequencing technologies, Oxford Nanopore Technologies (ONT) holds a unique position because the ONT MinION, an error-prone long-read sequencer, is associated with little to no capital cost. Here we show that we can make short-read Illumina libraries compatible with the ONT MinION by using the rolling circle to concatemeric consensus (R2C2) method to circularize and amplify the short library molecules. This results in longer DNA molecules containing tandem repeats of the original short library molecules. This longer DNA is ideally suited for the ONT MinION, and after sequencing, the tandem repeats in the resulting raw reads can be converted into high-accuracy consensus reads with similar error rates to that of the Illumina MiSeq. We highlight this capability by producing and benchmarking RNA-seq, ChIP-seq, and regular and target-enriched Tn5 libraries. We also explore the use of this approach for rapid evaluation of sequencing library metrics by implementing a real-time analysis workflow.