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.

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.

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.

Deficiency of the Heterogeneous Nuclear Ribonucleoprotein U locus leads to delayed hindbrain neurogenesis.

Genetic variants affecting Heterogeneous Nuclear Ribonucleoprotein U (HNRNPU) have been identified in several neurodevelopmental disorders (NDDs). HNRNPU is widely expressed in the human brain and shows the highest postnatal expression in the cerebellum. Recent studies have investigated the role of HNRNPU in cerebral cortical development, but the effects of HNRNPU deficiency on cerebellar development remain unknown. Here, we describe the molecular and cellular outcomes of HNRNPU locus deficiency during in vitro neural differentiation of patient-derived and isogenic neuroepithelial stem cells with a hindbrain profile. We demonstrate that HNRNPU deficiency leads to chromatin remodeling of A/B compartments, and transcriptional rewiring, partly by impacting exon inclusion during mRNA processing. Genomic regions affected by the chromatin restructuring and host genes of exon usage differences show a strong enrichment for genes implicated in epilepsies, intellectual disability, and autism. Lastly, we show that at the cellular level HNRNPU downregulation leads to an increased fraction of neural progenitors in the maturing neuronal population. We conclude that the HNRNPU locus is involved in delayed commitment of neural progenitors to differentiate in cell types with hindbrain profile.

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.

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.

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 .

The histone H3.1 variant regulates TONSOKU-mediated DNA repair during replication.

The tail of replication-dependent histone H3.1 varies from that of replication-independent H3.3 at the amino acid located at position 31 in plants and animals, but no function has been assigned to this residue to demonstrate a unique and conserved role for H3.1 during replication. We found that TONSOKU (TSK/TONSL), which rescues broken replication forks, specifically interacts with H3.1 via recognition of alanine 31 by its tetratricopeptide repeat domain. Our results indicate that genomic instability in the absence of ATXR5/ATXR6-catalyzed histone H3 lysine 27 monomethylation in plants depends on H3.1, TSK, and DNA polymerase theta (Pol θ). This work reveals an H3.1-specific function during replication and a common strategy used in multicellular eukaryotes for regulating post-replicative chromatin maturation and TSK, which relies on histone monomethyltransferases and reading of the H3.1 variant.

Identification and characterization of histones in Physarum polycephalum evidence a phylogenetic vicinity of Mycetozoans to the animal kingdom.

Physarum polycephalum belongs to Mycetozoans, a phylogenetic clade apart from the animal, plant and fungus kingdoms. Histones are nuclear proteins involved in genome organization and regulation and are among the most evolutionary conserved proteins within eukaryotes. Therefore, this raises the question of their conservation in Physarum and the position of this organism within the eukaryotic phylogenic tree based on histone sequences. We carried out a comprehensive study of histones in Physarum polycephalum using genomic, transcriptomic and molecular data. Our results allowed to identify the different isoforms of the core histones H2A, H2B, H3 and H4 which exhibit strong conservation of amino acid residues previously identified as subject to post-translational modifications. Furthermore, we also identified the linker histone H1, the most divergent histone, and characterized a large number of its PTMs by mass spectrometry. We also performed an in-depth investigation of histone genes and transcript structures. Histone proteins are highly conserved in Physarum and their characterization will contribute to a better understanding of the polyphyletic Mycetozoan group. Our data reinforce that P. polycephalum is evolutionary closer to animals than plants and located at the crown of the eukaryotic tree. Our study provides new insights in the evolutionary history of Physarum and eukaryote lineages.

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.

Automated 3D bio-imaging analysis of nuclear organization by NucleusJ 2.0.

NucleusJ 1.0, an ImageJ plugin, is a useful tool to analyze nuclear morphology and chromatin organization in plant and animal cells. NucleusJ 2.0 is a new release of NucleusJ, in which image processing is achieved more quickly using a command-lineuser interface. Starting with large collection of 3D nuclei, segmentation can be performed by the previously developed Otsu-modified method or by a new 3D gift-wrapping method, taking better account of nuclear indentations and unstained nucleoli. These two complementary methods are compared for their accuracy by using three types of datasets available to the community at https://www.brookes.ac.uk/indepth/images/ . Finally, NucleusJ 2.0 was evaluated using original plant genetic material by assessing its efficiency on nuclei stained with DNA dyes or after 3D-DNA Fluorescence in situ hybridization. With these improvements, NucleusJ 2.0 permits the generation of large user-curated datasets that will be useful for software benchmarking or to train convolution neural networks.

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.

Rapid Irreversible Transcriptional Reprogramming in Human Stem Cells Accompanied by Discordance between Replication Timing and Chromatin Compartment.

The temporal order of DNA replication is regulated during development and is highly correlated with gene expression, histone modifications and 3D genome architecture. We tracked changes in replication timing, gene expression, and chromatin conformation capture (Hi-C) A/B compartments over the first two cell cycles during differentiation of human embryonic stem cells to definitive endoderm. Remarkably, transcriptional programs were irreversibly reprogrammed within the first cell cycle and were largely but not universally coordinated with replication timing changes. Moreover, changes in A/B compartment and several histone modifications that normally correlate strongly with replication timing showed weak correlation during the early cell cycles of differentiation but showed increased alignment in later differentiation stages and in terminally differentiated cell lines. Thus, epigenetic cell fate transitions during early differentiation can occur despite dynamic and discordant changes in otherwise highly correlated genomic properties.

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.

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.

Computational Methods for Studying the Plant Nucleus.

The analysis of nuclear envelope components and their function has recently been progressed by the use of computational methods of analysis. The methods in this chapter provided by members of the International Plant Nucleus Consortium address the identification of novel nuclear envelope proteins and the study of structure and mobility of the nucleus. DORY2 is an upgrade of the KASH-finder DORY, and NucleusJ is used to characterize the three-dimensional structure of the nucleus in light microscope images. Finally, a method is provided for analysis of the migration of the nucleus, a key technique for exploring the function of plant nuclear proteins.

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.

Quantitative 3D Analysis of Nuclear Morphology and Heterochromatin Organization from Whole-Mount Plant Tissue Using NucleusJ.

Image analysis is a classical way to study nuclear organization. While nuclear organization used to be investigated by colorimetric or fluorescent labeling of DNA or specific nuclear compartments, new methods in microscopy imaging now enable qualitative and quantitative analyses of chromatin pattern, and nuclear size and shape. Several procedures have been developed to prepare samples in order to collect 3D images for the analysis of spatial chromatin organization, but only few preserve the positional information of the cell within its tissue context. Here, we describe a whole mount tissue preparation procedure coupled to DNA staining using the PicoGreen® intercalating agent suitable for image analysis of the nucleus in living and fixed tissues. 3D Image analysis is then performed using NucleusJ, an open source ImageJ plugin, which allows for quantifying variations in nuclear morphology such as nuclear volume, sphericity, elongation, and flatness as well as in heterochromatin content and position in respect to the nuclear periphery.

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.