piRNA generation is associated with the pioneer round of translation in stem cells.

Much insight has been gained on how stem cells maintain genomic integrity, but less attention has been paid to how they maintain their transcriptome. Here, we report that the PIWI protein SMEDWI-1 plays a role in the filtering of dysfunctional transcripts from the transcriptome of planarian stem cells. SMEDWI-1 accomplishes this through association with the ribosomes during the pioneer round of translation, and processing of poorly translated transcripts into piRNAs. This results in the removal of such transcripts from the cytoplasmic pool and at the same time creates a dynamic pool of small RNAs for post-transcriptional surveillance through the piRNA pathway. Loss of SMEDWI-1 results in elevated levels of several non-coding transcripts, including rRNAs, snRNAs and pseudogene mRNAs, while reducing levels of several coding transcripts. In the absence of SMEDWI-1, stem cell colonies are delayed in their expansion and a higher fraction of descendants exit the stem cell state, indicating that this transcriptomic sanitation mediated by SMEDWI-1 is essential to maintain stem cell health. This study presents a new model for the function of PIWI proteins in stem cell maintenance, that complements their role in transposon repression, and proposes a new biogenesis pathway for piRNAs in stem cells.

H3.1K27me1 maintains transcriptional silencing and genome stability by preventing GCN5-mediated histone acetylation.

Epigenetic mechanisms play diverse roles in the regulation of genome stability in eukaryotes. In Arabidopsis thaliana, genome stability is maintained during DNA replication by the H3.1K27 methyltransferases ARABIDOPSIS TRITHORAX-RELATED PROTEIN 5 (ATXR5) and ATXR6, which catalyze the deposition of K27me1 on replication-dependent H3.1 variants. The loss of H3.1K27me1 in atxr5 atxr6 double mutants leads to heterochromatin defects, including transcriptional de-repression and genomic instability, but the molecular mechanisms involved remain largely unknown. In this study, we identified the transcriptional co-activator and conserved histone acetyltransferase GCN5 as a mediator of transcriptional de-repression and genomic instability in the absence of H3.1K27me1. GCN5 is part of a SAGA-like complex in plants that requires the GCN5-interacting protein ADA2b and the chromatin remodeler CHR6 to mediate the heterochromatic defects in atxr5 atxr6 mutants. Our results also indicate that Arabidopsis GCN5 acetylates multiple lysine residues on H3.1 variants, but H3.1K27 and H3.1K36 play essential functions in inducing genomic instability in the absence of H3.1K27me1. Finally, we show that H3.1K36 acetylation by GCN5 is negatively regulated by H3.1K27me1 in vitro. Overall, this work reveals a key molecular role for H3.1K27me1 in maintaining transcriptional silencing and genome stability in heterochromatin by restricting GCN5-mediated histone acetylation in plants.

MtING2 encodes an ING domain PHD finger protein which affects Medicago growth, flowering, global patterns of H3K4me3, and gene expression.

Flowering of the reference legume Medicago truncatula is promoted by winter cold (vernalization) followed by long-day photoperiods (VLD) similar to winter annual Arabidopsis. However, Medicago lacks FLC and CO, key regulators of Arabidopsis VLD flowering. Most plants have two INHIBITOR OF GROWTH (ING) genes (ING1 and ING2), encoding proteins with an ING domain with two anti-parallel alpha-helices and a plant homeodomain (PHD) finger, but their genetic role has not been previously described. In Medicago, Mting1 gene-edited mutants developed and flowered normally, but an Mting2-1 Tnt1 insertion mutant and gene-edited Mting2 mutants had developmental abnormalities including delayed flowering particularly in VLD, compact architecture, abnormal leaves with extra leaflets but no trichomes, and smaller seeds and barrels. Mting2 mutants had reduced expression of activators of flowering, including the FT-like gene MtFTa1, and increased expression of the candidate repressor MtTFL1c, consistent with the delayed flowering of the mutant. MtING2 overexpression complemented Mting2-1, but did not accelerate flowering in wild type. The MtING2 PHD finger bound H3K4me2/3 peptides weakly in vitro, but analysis of gene-edited mutants indicated that it was dispensable to MtING2 function in wild-type plants. RNA sequencing experiments indicated that >7000 genes are mis-expressed in the Mting2-1 mutant, consistent with its strong mutant phenotypes. Interestingly, ChIP-seq analysis identified >5000 novel H3K4me3 locations in the genome of Mting2-1 mutants compared to wild type R108. Overall, our mutant study has uncovered an important physiological role of a plant ING2 gene in development, flowering, and gene expression, which likely involves an epigenetic mechanism.

Analysis of Hi-C data using SIP effectively identifies loops in organisms from C. elegans to mammals.

Chromatin loops are a major component of 3D nuclear organization, visually apparent as intense point-to-point interactions in Hi-C maps. Identification of these loops is a critical part of most Hi-C analyses. However, current methods often miss visually evident CTCF loops in Hi-C data sets from mammals, and they completely fail to identify high intensity loops in other organisms. We present SIP, Significant Interaction Peak caller, and SIPMeta, which are platform independent programs to identify and characterize these loops in a time- and memory-efficient manner. We show that SIP is resistant to noise and sequencing depth, and can be used to detect loops that were previously missed in human cells as well as loops in other organisms. SIPMeta corrects for a common visualization artifact by accounting for Manhattan distance to create average plots of Hi-C and HiChIP data. We then demonstrate that the use of SIP and SIPMeta can lead to biological insights by characterizing the contribution of several transcription factors to CTCF loop stability in human cells. We also annotate loops associated with the SMC component of the dosage compensation complex (DCC) in Caenorhabditis elegans and demonstrate that loop anchors represent bidirectional blocks for symmetrical loop extrusion. This is in contrast to the asymmetrical extrusion until unidirectional blockage by CTCF that is presumed to occur in mammals. Using HiChIP and multiway ligation events, we then show that DCC loops form a network of strong interactions that may contribute to X Chromosome-wide condensation in C. elegans hermaphrodites.

The Histone Chaperone Network Is Highly Conserved in Physarum polycephalum.

The nucleosome is composed of histones and DNA. Prior to their deposition on chromatin, histones are shielded by specialized and diverse proteins known as histone chaperones. They escort histones during their entire cellular life and ensure their proper incorporation in chromatin. Physarum polycephalum is a Mycetozoan, a clade located at the crown of the eukaryotic tree. We previously found that histones, which are highly conserved between plants and animals, are also highly conserved in Physarum. However, histone chaperones differ significantly between animal and plant kingdoms, and this thus probed us to further study the conservation of histone chaperones in Physarum and their evolution relative to animal and plants. Most of the known histone chaperones and their functional domains are conserved as well as key residues required for histone and chaperone interactions. Physarum is divergent from yeast, plants and animals, but PpHIRA, PpCABIN1 and PpSPT6 are similar in structure to plant orthologues. PpFACT is closely related to the yeast complex, and the Physarum genome encodes the animal-specific APFL chaperone. Furthermore, we performed RNA sequencing to monitor chaperone expression during the cell cycle and uncovered two distinct patterns during S-phase. In summary, our study demonstrates the conserved role of histone chaperones in handling histones in an early-branching eukaryote.

NODeJ: an ImageJ plugin for 3D segmentation of nuclear objects.

BACKGROUND: The three-dimensional nuclear arrangement of chromatin impacts many cellular processes operating at the DNA level in animal and plant systems. Chromatin organization is a dynamic process that can be affected by biotic and abiotic stresses. Three-dimensional imaging technology allows to follow these dynamic changes, but only a few semi-automated processing methods currently exist for quantitative analysis of the 3D chromatin organization. RESULTS: We present an automated method, Nuclear Object DetectionJ (NODeJ), developed as an imageJ plugin. This program segments and analyzes high intensity domains in nuclei from 3D images. NODeJ performs a Laplacian convolution on the mask of a nucleus to enhance the contrast of intra-nuclear objects and allow their detection. We reanalyzed public datasets and determined that NODeJ is able to accurately identify heterochromatin domains from a diverse set of Arabidopsis thaliana nuclei stained with DAPI or Hoechst. NODeJ is also able to detect signals in nuclei from DNA FISH experiments, allowing for the analysis of specific targets of interest. CONCLUSION AND AVAILABILITY: NODeJ allows for efficient automated analysis of subnuclear structures by avoiding the semi-automated steps, resulting in reduced processing time and analytical bias. NODeJ is written in Java and provided as an ImageJ plugin with a command line option to perform more high-throughput analyses. NODeJ can be downloaded from https://gitlab.com/axpoulet/image2danalysis/-/releases with source code, documentation and further information avaliable at https://gitlab.com/axpoulet/image2danalysis . The images used in this study are publicly available at https://www.brookes.ac.uk/indepth/images/ and https://doi.org/10.15454/1HSOIE .

RT States: systematic annotation of the human genome using cell type-specific replication timing programs.

MOTIVATION: The replication timing (RT) program has been linked to many key biological processes including cell fate commitment, 3D chromatin organization and transcription regulation. Significant technology progress now allows to characterize the RT program in the entire human genome in a high-throughput and high-resolution fashion. These experiments suggest that RT changes dynamically during development in coordination with gene activity. Since RT is such a fundamental biological process, we believe that an effective quantitative profile of the local RT program from a diverse set of cell types in various developmental stages and lineages can provide crucial biological insights for a genomic locus. RESULTS: In this study, we explored recurrent and spatially coherent combinatorial profiles from 42 RT programs collected from multiple lineages at diverse differentiation states. We found that a Hidden Markov Model with 15 hidden states provide a good model to describe these genome-wide RT profiling data. Each of the hidden state represents a unique combination of RT profiles across different cell types which we refer to as 'RT states'. To understand the biological properties of these RT states, we inspected their relationship with chromatin states, gene expression, functional annotation and 3D chromosomal organization. We found that the newly defined RT states possess interesting genome-wide functional properties that add complementary information to the existing annotation of the human genome. AVAILABILITY AND IMPLEMENTATION: R scripts for inferring HMM models and Perl scripts for further analysis are available https://github.com/PouletAxel/script_HMM_Replication_timing. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Arabidopsis ATRX Modulates H3.3 Occupancy and Fine-Tunes Gene Expression.

Histones are essential components of the nucleosome, the major chromatin subunit that structures linear DNA molecules and regulates access of other proteins to DNA. Specific histone chaperone complexes control the correct deposition of canonical histones and their variants to modulate nucleosome structure and stability. In this study, we characterize the Arabidopsis thaliana Alpha Thalassemia-mental Retardation X-linked (ATRX) ortholog and show that ATRX is involved in histone H3 deposition. Arabidopsis ATRX mutant alleles are viable, but show developmental defects and reduced fertility. Their combination with mutants of the histone H3.3 chaperone HIRA (Histone Regulator A) results in impaired plant survival, suggesting that HIRA and ATRX function in complementary histone deposition pathways. Indeed, ATRX loss of function alters cellular histone H3.3 pools and in consequence modulates the H3.1/H3.3 balance in the cell. H3.3 levels are affected especially at genes characterized by elevated H3.3 occupancy, including the 45S ribosomal DNA (45S rDNA) loci, where loss of ATRX results in altered expression of specific 45S rDNA sequence variants. At the genome-wide scale, our data indicate that ATRX modifies gene expression concomitantly to H3.3 deposition at a set of genes characterized both by elevated H3.3 occupancy and high expression. Together, our results show that ATRX is involved in H3.3 deposition and emphasize the role of histone chaperones in adjusting genome expression.

Genetic and epigenetic variation in 5S ribosomal RNA genes reveals genome dynamics in Arabidopsis thaliana.

Organized in tandem repeat arrays in most eukaryotes and transcribed by RNA polymerase III, expression of 5S rRNA genes is under epigenetic control. To unveil mechanisms of transcriptional regulation, we obtained here in depth sequence information on 5S rRNA genes from the Arabidopsis thaliana genome and identified differential enrichment in epigenetic marks between the three 5S rDNA loci situated on chromosomes 3, 4 and 5. We reveal the chromosome 5 locus as the major source of an atypical, long 5S rRNA transcript characteristic of an open chromatin structure. 5S rRNA genes from this locus translocated in the Landsberg erecta ecotype as shown by linkage mapping and chromosome-specific FISH analysis. These variations in 5S rDNA locus organization cause changes in the spatial arrangement of chromosomes in the nucleus. Furthermore, 5S rRNA gene arrangements are highly dynamic with alterations in chromosomal positions through translocations in certain mutants of the RNA-directed DNA methylation pathway and important copy number variations among ecotypes. Finally, variations in 5S rRNA gene sequence, chromatin organization and transcripts indicate differential usage of 5S rDNA loci in distinct ecotypes. We suggest that both the usage of existing and new 5S rDNA loci resulting from translocations may impact neighboring chromatin organization.

The LINC complex contributes to heterochromatin organisation and transcriptional gene silencing in plants.

The linker of nucleoskeleton and cytoskeleton (LINC) complex is an evolutionarily well-conserved protein bridge connecting the cytoplasmic and nuclear compartments across the nuclear membrane. While recent data support its function in nuclear morphology and meiosis, its involvement in chromatin organisation has not been studied in plants. Here, 3D imaging methods have been used to investigate nuclear morphology and chromatin organisation in interphase nuclei of the model plant Arabidopsis thaliana in which heterochromatin clusters in conspicuous chromatin domains called chromocentres. Chromocentres form a repressive chromatin environment contributing to transcriptional silencing of repeated sequences, a general mechanism needed for genome stability. Quantitative measurements of the 3D position of chromocentres indicate their close proximity to the nuclear periphery but that their position varies with nuclear volume and can be altered in specific mutants affecting the LINC complex. Finally, we propose that the plant LINC complex contributes to proper heterochromatin organisation and positioning at the nuclear periphery, since its alteration is associated with the release of transcriptional silencing as well as decompaction of heterochromatic sequences.

Replication-coupled histone H3.1 deposition determines nucleosome composition and heterochromatin dynamics during Arabidopsis seedling development.

Developmental phase transitions are often characterized by changes in the chromatin landscape and heterochromatin reorganization. In Arabidopsis, clustering of repetitive heterochromatic loci into so-called chromocenters is an important determinant of chromosome organization in nuclear space. Here, we investigated the molecular mechanisms involved in chromocenter formation during the switch from a heterotrophic to a photosynthetically competent state during early seedling development. We characterized the spatial organization and chromatin features at centromeric and pericentromeric repeats and identified mutant contexts with impaired chromocenter formation. We find that clustering of repetitive DNA loci into chromocenters takes place in a precise temporal window and results in reinforced transcriptional repression. Although repetitive sequences are enriched in H3K9me2 and linker histone H1 before repeat clustering, chromocenter formation involves increasing enrichment in H3.1 as well as H2A.W histone variants, hallmarks of heterochromatin. These processes are severely affected in mutants impaired in replication-coupled histone assembly mediated by CHROMATIN ASSEMBLY FACTOR 1 (CAF-1). We further reveal that histone deposition by CAF-1 is required for efficient H3K9me2 enrichment at repetitive sequences during chromocenter formation. Taken together, we show that chromocenter assembly during post-germination development requires dynamic changes in nucleosome composition and histone post-translational modifications orchestrated by the replication-coupled H3.1 deposition machinery.

The cell wall of Arabidopsis thaliana influences actin network dynamics.

In plant cells, molecular connections link the cell wall-plasma membrane-actin cytoskeleton to form a continuum. It is hypothesized that the cell wall provides stable anchor points around which the actin cytoskeleton remodels. Here we use live cell imaging of fluorescently labelled marker proteins to quantify the organization and dynamics of the actin cytoskeleton and to determine the impact of disrupting connections within the continuum. Labelling of the actin cytoskeleton with green fluorescent protein (GFP)-fimbrin actin-binding domain 2 (FABD2) resulted in a network composed of fine filaments and thicker bundles that appeared as a highly dynamic remodelling meshwork. This differed substantially from the GFP-Lifeact-labelled network that appeared much more sparse with thick bundles that underwent 'simple movement', in which the bundles slightly change position, but in such a manner that the structure of the network was not substantially altered during the time of observation. Label-dependent differences in actin network morphology and remodelling necessitated development of two new image analysis techniques. The first of these, 'pairwise image subtraction', was applied to measurement of the more rapidly remodelling actin network labelled with GFP-FABD2, while the second, 'cumulative fluorescence intensity', was used to measure bulk remodelling of the actin cytoskeleton when labelled with GFP-Lifeact. In each case, these analysis techniques show that the actin cytoskeleton has a decreased rate of bulk remodelling when the cell wall-plasma membrane-actin continuum is disrupted either by plasmolysis or with isoxaben, a drug that specifically inhibits cellulose deposition. Changes in the rate of actin remodelling also affect its functionality, as observed by alteration in Golgi body motility.

Bioinformatics Analysis of Phylogeny and Transcription of TAA/YUC Auxin Biosynthetic Genes.

Auxin is a main plant growth hormone crucial in a multitude of developmental processes in plants. Auxin biosynthesis via the tryptophan aminotransferase of arabidopsis (TAA)/YUCCA (YUC) route involving tryptophan aminotransferases and YUC flavin-dependent monooxygenases that produce the auxin indole-3-acetic acid (IAA) from tryptophan is currently the most researched auxin biosynthetic pathway. Previous data showed that, in maize and arabidopsis, TAA/YUC-dependent auxin biosynthesis can be detected in endoplasmic reticulum (ER) microsomal fractions, and a subset of auxin biosynthetic proteins are localized to the ER, mainly due to transmembrane domains (TMD). The phylogeny presented here for TAA/TAR (tryptophan aminotransferase related) and YUC proteins analyses phylogenetic groups as well as transmembrane domains for ER-membrane localisation. In addition, RNAseq datasets are analysed for transcript abundance of YUC and TAA/TAR proteins in Arabidopsis thaliana. We show that ER membrane localisation for TAA/YUC proteins involved in auxin biosynthesis is already present early on in the evolution of mosses and club mosses. ER membrane anchored YUC proteins can mainly be found in roots, while cytosolic proteins are more abundant in the shoot. The distribution between the different phylogenetic classes in root and shoot may well originate from gene duplications, and the phylogenetic groups detected also overlap with the biological function.

The histone chaperone complex HIR maintains nucleosome occupancy and counterbalances impaired histone deposition in CAF-1 complex mutants.

Chromatin organization is essential for coordinated gene expression, genome stability, and inheritance of epigenetic information. The main components involved in chromatin assembly are specific complexes such as Chromatin Assembly Factor 1 (CAF-1) and Histone Regulator (HIR), which deposit histones in a DNA synthesis-dependent or -independent manner, respectively. Here, we characterize the role of the plant orthologs Histone Regulator A (HIRA), Ubinuclein (UBN) and Calcineurin Binding protein 1 (CABIN1), which constitute the HIR complex. Arabidopsis loss-of-function mutants for the various subunits of the complex are viable, but hira mutants show reduced fertility. We show that loss of HIRA reduces extractable histone H3 protein levels and decreases nucleosome occupancy at both actively transcribed genes and heterochromatic regions. Concomitantly, HIRA contributes to maintenance of silencing of pericentromeric repeats and certain transposons. A genetic analysis based on crosses between mutants deficient in subunits of the CAF-1 and HIR complexes showed that simultaneous loss of both the CAF-1 and HIR histone H3 chaperone complexes severely affects plant survival, growth and reproductive development. Our results suggest that HIRA partially rescues impaired histone deposition in fas mutants to preserve nucleosome occupancy, implying plasticity in histone variant interaction and deposition.

NucleusJ: an ImageJ plugin for quantifying 3D images of interphase nuclei.

UNLABELLED: NucleusJ is a simple and user-friendly ImageJ plugin dedicated to the characterization of nuclear morphology and chromatin organization in 3D. Starting from image stacks, the nuclear boundary is delimited by combining the Otsu segmentation method with optimization of nuclear sphericity. Chromatin domains are segmented by partitioning the nucleus using a 3D watershed algorithm and by thresholding a contrast measure over the resulting regions. As output, NucleusJ quantifies 15 parameters including shape and size of nuclei as well as intra-nuclear objects and their position within the nucleus. A step-by-step documentation is available for self-training, together with data sets of nuclei with different nuclear organization. AVAILABILITY AND IMPLEMENTATION: Dataset of nuclei is available at https://www.gred-clermont.fr/media/WorkDirectory.zip. NucleusJ is available at http://imagejdocu.tudor.lu/doku.php?id=plugin:stacks:nuclear_analysis_plugin:start. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

A novel family of plant nuclear envelope-associated proteins.

This paper describes the characterisation of a new family of higher plant nuclear envelope-associated proteins (NEAPs) that interact with other proteins of the nuclear envelope. In the model plant Arabidopsis thaliana, the family consists of three genes expressed ubiquitously (AtNEAP1-3) and a pseudogene (AtNEAP4). NEAPs consist of extensive coiled-coil domains, followed by a nuclear localisation signal and a C-terminal predicted transmembrane domain. Domain deletion mutants confirm the presence of a functional nuclear localisation signal and transmembrane domain. AtNEAP proteins localise to the nuclear periphery as part of stable protein complexes, are able to form homo- and heteromers, and interact with the SUN domain proteins AtSUN1 and AtSUN2, involved in the linker of nucleoskeleton and cytoskeleton (LINC) complex. An A. thaliana cDNA library screen identified a putative transcription factor called AtbZIP18 as a novel interactor of AtNEAP1, which suggest a connection between NEAP and chromatin. An Atneap1 Atneap3 double-knockout mutant showed reduced root growth, and altered nuclear morphology and chromatin structure. Thus AtNEAPs are suggested as inner nuclear membrane-anchored coiled-coil proteins with roles in maintaining nuclear morphology and chromatin structure.

Evolutionary history of Methyltransferase 1 genes in hexaploid wheat.

BACKGROUND: Plant and animal methyltransferases are key enzymes involved in DNA methylation at cytosine residues, required for gene expression control and genome stability. Taking advantage of the new sequence surveys of the wheat genome recently released by the International Wheat Genome Sequencing Consortium, we identified and characterized MET1 genes in the hexaploid wheat Triticum aestivum (TaMET1). RESULTS: Nine TaMET1 genes were identified and mapped on homoeologous chromosome groups 2A/2B/2D, 5A/5B/5D and 7A/7B/7D. Synteny analysis and evolution rates suggest that the genome organization of TaMET1 genes results from a whole genome duplication shared within the grass family, and a second gene duplication, which occurred specifically in the Triticeae tribe prior to the speciation of diploid wheat. Higher expression levels were observed for TaMET1 homoeologous group 2 genes compared to group 5 and 7, indicating that group 2 homoeologous genes are predominant at the transcriptional level, while group 5 evolved into pseudogenes. We show the connection between low expression levels, elevated evolution rates and unexpected enrichment in CG-dinucleotides (CG-rich isochores) at putative promoter regions of homoeologous group 5 and 7, but not of group 2 TaMET1 genes. Bisulfite sequencing reveals that these CG-rich isochores are highly methylated in a CG context, which is the expected target of TaMET1. CONCLUSIONS: We retraced the evolutionary history of MET1 genes in wheat, explaining the predominance of group 2 homoeologous genes and suggest CG-DNA methylation as one of the mechanisms involved in wheat genome dynamics.

Gene-edited Mtsoc1 triple mutant Medicago plants do not flower.

Optimized flowering time is an important trait that ensures successful plant adaptation and crop productivity. SOC1-like genes encode MADS transcription factors, which are known to play important roles in flowering control in many plants. This includes the best-characterized eudicot model Arabidopsis thaliana (Arabidopsis), where SOC1 promotes flowering and functions as a floral integrator gene integrating signals from different flowering-time regulatory pathways. Medicago truncatula (Medicago) is a temperate reference legume with strong genomic and genetic resources used to study flowering pathways in legumes. Interestingly, despite responding to similar floral-inductive cues of extended cold (vernalization) followed by warm long days (VLD), such as in winter annual Arabidopsis, Medicago lacks FLC and CO which are key regulators of flowering in Arabidopsis. Unlike Arabidopsis with one SOC1 gene, multiple gene duplication events have given rise to three MtSOC1 paralogs within the Medicago genus in legumes: one Fabaceae group A SOC1 gene, MtSOC1a, and two tandemly repeated Fabaceae group B SOC1 genes, MtSOC1b and MtSOC1c. Previously, we showed that MtSOC1a has unique functions in floral promotion in Medicago. The Mtsoc1a Tnt1 retroelement insertion single mutant showed moderately delayed flowering in long- and short-day photoperiods, with and without prior vernalization, compared to the wild-type. In contrast, Mtsoc1b Tnt1 single mutants did not have altered flowering time or flower development, indicating that it was redundant in an otherwise wild-type background. Here, we describe the generation of Mtsoc1a Mtsoc1b Mtsoc1c triple mutant lines using CRISPR-Cas9 gene editing. We studied two independent triple mutant lines that segregated plants that did not flower and were bushy under floral inductive VLD. Genotyping indicated that these non-flowering plants were homozygous for the predicted strong mutant alleles of the three MtSOC1 genes. Gene expression analyses using RNA-seq and RT-qPCR indicated that these plants remained vegetative. Overall, the non-flowering triple mutants were dramatically different from the single Mtsoc1a mutant and the Arabidopsis soc1 mutant; implicating multiple MtSOC1 genes in critical overlapping roles in the transition to flowering in Medicago.

A comparative roadmap of PIWI-interacting RNAs across seven species reveals insights into de novo piRNA-precursor formation in mammals.

PIWI-interacting RNAs (piRNAs) play a crucial role in safeguarding genome integrity by silencing mobile genetic elements. From flies to humans, piRNAs originate from long single-stranded precursors encoded by genomic piRNA clusters. How piRNA clusters form to adapt to genomic invaders and evolve to maintain protection remain key outstanding questions. Here, we generate a roadmap of piRNA clusters across seven species that highlights both similarities and variations. In mammals, we identify transcriptional readthrough as a mechanism to generate piRNAs from transposon insertions (piCs) downstream of genes (DoG). Together with the well-known stress-dependent DoG transcripts, our findings suggest a molecular mechanism for the formation of piRNA clusters in response to retroviral invasion. Finally, we identify a class of dynamic piRNA clusters in humans, underscoring unique features of human germ cell biology. Our results advance the understanding of conserved principles and species-specific variations in piRNA biology and provide tools for future studies.

Chromatin analysis of adult pluripotent stem cells reveals a unique stemness maintenance strategy.

Many highly regenerative organisms maintain adult pluripotent stem cells throughout their life, but how the long-term maintenance of pluripotency is accomplished is unclear. To decipher the regulatory logic of adult pluripotent stem cells, we analyzed the chromatin organization of stem cell genes in the planarian Schmidtea mediterranea. We identify a special chromatin state of stem cell genes, which is distinct from that of tissue-specific genes and resembles constitutive genes. Where tissue-specific promoters have detectable transcription factor binding sites, the promoters of stem cell-specific genes instead have sequence features that broadly decrease nucleosome binding affinity. This genic organization makes pluripotency-related gene expression the default state in these cells, which is maintained by the activity of chromatin remodelers ISWI and SNF2 in the stem cells.