Epigenetic and transcriptional consequences in the endosperm of chemically induced transposon mobilization in Arabidopsis.

Genomic imprinting, an epigenetic phenomenon leading to parent-of-origin-specific gene expression, has independently evolved in the endosperm of flowering plants and the placenta of mammals-tissues crucial for nurturing embryos. While transposable elements (TEs) frequently colocalize with imprinted genes and are implicated in imprinting establishment, direct investigations of the impact of de novo TE transposition on genomic imprinting remain scarce. In this study, we explored the effects of chemically induced transposition of the Copia element ONSEN on genomic imprinting in Arabidopsis thaliana. Through the combination of chemical TE mobilization and doubled haploid induction, we generated a line with 40 new ONSEN copies. Our findings reveal a preferential targeting of maternally expressed genes (MEGs) for transposition, aligning with the colocalization of H2A.Z and H3K27me3 in MEGs-both previously identified as promoters of ONSEN insertions. Additionally, we demonstrate that chemically-induced DNA hypomethylation induces global transcriptional deregulation in the endosperm, leading to the breakdown of MEG imprinting. This study provides insights into the consequences of chemically induced TE remobilization in the endosperm, revealing that chemically-induced epigenome changes can have long-term consequences on imprinted gene expression.

Reliable genotyping of recombinant genomes using a robust hidden Markov model.

Meiotic recombination is an essential mechanism during sexual reproduction and includes the exchange of chromosome segments between homologous chromosomes. New allelic combinations are transmitted to the new generation, introducing novel genetic variation in the offspring genomes. With the improvement of high-throughput whole-genome sequencing technologies, large numbers of recombinant individuals can now be sequenced with low sequencing depth at low costs, necessitating computational methods for reconstructing their haplotypes. The main challenge is the uncertainty in haplotype calling that arises from the low information content of a single genomic position. Straightforward sliding window-based approaches are difficult to tune and fail to place recombination breakpoints precisely. Hidden Markov model (HMM)-based approaches, on the other hand, tend to over-segment the genome. Here, we present RTIGER, an HMM-based model that exploits in a mathematically precise way the fact that true chromosome segments typically have a certain minimum length. We further separate the task of identifying the correct haplotype sequence from the accurate placement of haplotype borders, thereby maximizing the accuracy of border positions. By comparing segmentations based on simulated data with known underlying haplotypes, we highlight the reasons for RTIGER outperforming traditional segmentation approaches. We then analyze the meiotic recombination pattern of segregants of 2 Arabidopsis (Arabidopsis thaliana) accessions and a previously described hyper-recombining mutant. RTIGER is available as an R package with an efficient Julia implementation of the core algorithm.

Anno genominis XX: 20 years of Arabidopsis genomics.

Twenty years ago, the Arabidopsis thaliana genome sequence was published. This was an important moment as it was the first sequenced plant genome and explicitly brought plant science into the genomics era. At the time, this was not only an outstanding technological achievement, but it was characterized by a superb global collaboration. The Arabidopsis genome was the seed for plant genomic research. Here, we review the development of numerous resources based on the genome that have enabled discoveries across plant species, which has enhanced our understanding of how plants function and interact with their environments.

Transposition and duplication of MADS-domain transcription factor genes in annual and perennial Arabis species modulates flowering.

The timing of reproduction is an adaptive trait in many organisms. In plants, the timing, duration, and intensity of flowering differ between annual and perennial species. To identify interspecies variation in these traits, we studied introgression lines derived from hybridization of annual and perennial species, Arabis montbretiana and Arabis alpina, respectively. Recombination mapping identified two tandem A. montbretiana genes encoding MADS-domain transcription factors that confer extreme late flowering on A. alpina These genes are related to the MADS AFFECTING FLOWERING (MAF) cluster of floral repressors of other Brassicaceae species and were named A. montbretiana (Am) MAF-RELATED (MAR) genes. AmMAR1 but not AmMAR2 prevented floral induction at the shoot apex of A. alpina, strongly enhancing the effect of the MAF cluster, and MAR1 is absent from the genomes of all A. alpina accessions analyzed. Exposure of plants to cold (vernalization) represses AmMAR1 transcription and overcomes its inhibition of flowering. Assembly of the tandem arrays of MAR and MAF genes of six A. alpina accessions and three related species using PacBio long-sequence reads demonstrated that the MARs arose within the Arabis genus by interchromosomal transposition of a MAF1-like gene followed by tandem duplication. Time-resolved comparative RNA-sequencing (RNA-seq) suggested that AmMAR1 may be retained in A. montbretiana to enhance the effect of the AmMAF cluster and extend the duration of vernalization required for flowering. Our results demonstrate that MAF genes transposed independently in different Brassicaceae lineages and suggest that they were retained to modulate adaptive flowering responses that differ even among closely related species.

Barley Ror1 encodes a class XI myosin required for mlo-based broad-spectrum resistance to the fungal powdery mildew pathogen.

Loss-of-function alleles of plant MLO genes confer broad-spectrum resistance to powdery mildews in many eudicot and monocot species. Although barley (Hordeum vulgare) mlo mutants have been used in agriculture for more than 40 years, understanding of the molecular principles underlying this type of disease resistance remains fragmentary. Forward genetic screens in barley have revealed mutations in two Required for mlo resistance (Ror) genes that partially impair immunity conferred by mlo mutants. While Ror2 encodes a soluble N-ethylmaleimide-sensitive factor-attached protein receptor (SNARE), the identity of Ror1, located at the pericentromeric region of barley chromosome 1H, remained elusive. We report the identification of Ror1 based on combined barley genomic sequence information and transcriptomic data from ror1 mutant plants. Ror1 encodes the barley class XI myosin Myo11A (HORVU.MOREX.r3.1HG0046420). Single amino acid substitutions of this myosin, deduced from non-functional ror1 mutant alleles, map to the nucleotide-binding region and the interface between the relay-helix and the converter domain of the motor protein. Ror1 myosin accumulates transiently in the course of powdery mildew infection. Functional fluorophore-labeled Ror1 variants associate with mobile intracellular compartments that partially colocalize with peroxisomes. Single-cell expression of the Ror1 tail region causes a dominant-negative effect that phenocopies ror1 loss-of-function mutants. We define a myosin motor for the establishment of mlo-mediated resistance, suggesting that motor protein-driven intracellular transport processes are critical for extracellular immunity, possibly through the targeted transfer of antifungal and/or cell wall cargoes to pathogen contact sites.

plotsr: visualizing structural similarities and rearrangements between multiple genomes.

SUMMARY: Third-generation genome sequencing technologies have led to a sharp increase in the number of high-quality genome assemblies. This allows the comparison of multiple assembled genomes of individual species and demands new tools for visualizing their structural properties. Here, we present plotsr, an efficient tool to visualize structural similarities and rearrangements between genomes. It can be used to compare genomes on chromosome level or to zoom in on any selected region. In addition, plotsr can augment the visualization with regional identifiers (e.g. genes or genomic markers) or histogram tracks for continuous features (e.g. GC content or polymorphism density). AVAILABILITY AND IMPLEMENTATION: plotsr is implemented as a python package and uses the standard matplotlib library for plotting. It is freely available under the MIT license at GitHub (https://github.com/schneebergerlab/plotsr) and bioconda (https://anaconda.org/bioconda/plotsr). SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Mutagenesis of a Quintuple Mutant Impaired in Environmental Responses Reveals Roles for CHROMATIN REMODELING4 in the Arabidopsis Floral Transition.

Several pathways conferring environmental flowering responses in Arabidopsis (Arabidopsis thaliana) converge on developmental processes that mediate the floral transition in the shoot apical meristem. Many characterized mutations disrupt these environmental responses, but downstream developmental processes have been more refractory to mutagenesis. Here, we constructed a quintuple mutant impaired in several environmental pathways and showed that it possesses severely reduced flowering responses to changes in photoperiod and ambient temperature. RNA-sequencing (RNA-seq) analysis of the quintuple mutant showed that the expression of genes encoding gibberellin biosynthesis enzymes and transcription factors involved in the age pathway correlates with flowering. Mutagenesis of the quintuple mutant generated two late-flowering mutants, quintuple ems1 (qem1) and qem2 The mutated genes were identified by isogenic mapping and transgenic complementation. The qem1 mutant is an allele of the gibberellin 20-oxidase gene ga20ox2, confirming the importance of gibberellin for flowering in the absence of environmental responses. By contrast, qem2 is impaired in CHROMATIN REMODELING4 (CHR4), which has not been genetically implicated in floral induction. Using co-immunoprecipitation, RNA-seq, and chromatin immunoprecipitation sequencing, we show that CHR4 interacts with transcription factors involved in floral meristem identity and affects the expression of key floral regulators. Therefore, CHR4 mediates the response to endogenous flowering pathways in the inflorescence meristem to promote floral identity.

The Evolutionary Dynamics of Genetic Incompatibilities Introduced by Duplicated Genes in Arabidopsis thaliana.

Although gene duplications provide genetic backup and allow genomic changes under relaxed selection, they may potentially limit gene flow. When different copies of a duplicated gene are pseudofunctionalized in different genotypes, genetic incompatibilities can arise in their hybrid offspring. Although such cases have been reported after manual crosses, it remains unclear whether they occur in nature and how they affect natural populations. Here, we identified four duplicated-gene based incompatibilities including one previously not reported within an artificial Arabidopsis intercross population. Unexpectedly, however, for each of the genetic incompatibilities we also identified the incompatible alleles in natural populations based on the genomes of 1,135 Arabidopsis accessions published by the 1001 Genomes Project. Using the presence of incompatible allele combinations as phenotypes for GWAS, we mapped genomic regions that included additional gene copies which likely rescue the genetic incompatibility. Reconstructing the geographic origins and evolutionary trajectories of the individual alleles suggested that incompatible alleles frequently coexist, even in geographically closed regions, and that their effects can be overcome by additional gene copies collectively shaping the evolutionary dynamics of duplicated genes during population history.

FLOWERING REPRESSOR AAA+ ATPase 1 is a novel regulator of perennial flowering in Arabis alpina.

Arabis alpina is a polycarpic perennial, in which PERPETUAL FLOWERING1 (PEP1) regulates flowering and perennial traits in a vernalization-dependent manner. Mutagenesis screens of the pep1 mutant established the role of other flowering time regulators in PEP1-parallel pathways. Here we characterized three allelic enhancers of pep1 (eop002, 085 and 091) which flower early. We mapped the causal mutations and complemented mutants with the identified gene. Using quantitative reverse transcriptase PCR and reporter lines, we determined the protein spatiotemporal expression patterns and localization within the cell. We also characterized its role in Arabidopsis thaliana using CRISPR and in A. alpina by introgressing mutant alleles into a wild-type background. These mutants carried lesions in an AAA+ ATPase of unknown function, FLOWERING REPRESSOR AAA+ ATPase 1 (AaFRAT1). AaFRAT1 was detected in the vasculature of young leaf primordia and the rib zone of flowering shoot apical meristems. At the subcellular level, AaFRAT1 was localized at the interphase between the endoplasmic reticulum and peroxisomes. Introgression lines carrying Aafrat1 alleles required less vernalization to flower and reduced number of vegetative axillary branches. By contrast, A. thaliana CRISPR lines showed weak flowering phenotypes. AaFRAT1 contributes to flowering time regulation and the perennial growth habit of A. alpina.

Genetic and molecular analysis of trichome development in Arabis alpina.

The genetic and molecular analysis of trichome development in Arabidopsis thaliana has generated a detailed knowledge about the underlying regulatory genes and networks. However, how rapidly these mechanisms diverge during evolution is unknown. To address this problem, we used an unbiased forward genetic approach to identify most genes involved in trichome development in the related crucifer species Arabisalpina In general, we found most trichome mutant classes known in A. thaliana We identified orthologous genes of the relevant A. thaliana genes by sequence similarity and synteny and sequenced candidate genes in the A. alpina mutants. While in most cases we found a highly similar gene-phenotype relationship as known from Arabidopsis, there were also striking differences in the regulation of trichome patterning, differentiation, and morphogenesis. Our analysis of trichome patterning suggests that the formation of two classes of trichomes is regulated differentially by the homeodomain transcription factor AaGL2 Moreover, we show that overexpression of the GL3 basic helix-loop-helix transcription factor in A. alpina leads to the opposite phenotype as described in A. thaliana Mathematical modeling helps to explain how this nonintuitive behavior can be explained by different ratios of GL3 and GL1 in the two species.

Mutations in the miR396 binding site of the growth-regulating factor gene VvGRF4 modulate inflorescence architecture in grapevine.

Bunch rot caused by Botrytis cinerea infections is a notorious problem in grapevine cultivation. To produce high quality fruits, grapevine plants are treated with fungicides, which is cost intensive and harmful to the environment. Conversely, loose cluster bunches show a considerably enhanced physical resilience to bunch diseases. With the aim to identify genetic determinants that modulate the development of bunch architecture, we have compared loose and compact 'Pinot noir' clones. Loose cluster architecture was found to be correlated with increased berry size, elongated rachis and elongated pedicels. Using transcriptome analysis in combination with whole genome sequencing, we have identified a growth-regulating factor gene, VvGRF4, upregulated and harbours heterozygous mutations in the loose cluster clones. At late stages of inflorescence development, the mRNA pools of loose cluster clones contain predominantly mRNAs derived from the mutated alleles, which are resistant to miR396 degradation. Expression of the VvGRF4 gene and its mutated variants in Arabidopsis demonstrates that it promotes pedicel elongation. Taken together, VvGRF4 modulates bunch architecture in grapevine 'Pinot noir' clones. This trait can be introduced into other cultivars using marker-assisted breeding or CRISPR-Cas9 technology. Related growth-regulating factors or other genes of the same pathway may have similar functions.

A mis-regulated cyclic nucleotide-gated channel mediates cytosolic calcium elevation and activates immunity in Arabidopsis.

Calcium (Ca2+ ) is a second messenger for plant cell surface and intracellular receptors mediating pattern-triggered and effector-triggered immunity (respectively, PTI and ETI). Several CYCLIC NUCLEOTIDE-GATED CHANNELS (CNGCs) were shown to control transient cytosolic Ca2+ influx upon PTI activation. The contributions of specific CNGC members to PTI and ETI remain unclear. ENHANCED DISEASE SUSCEPTIBLITY1 (EDS1) regulates ETI signaling. In an Arabidopsis genetic screen for suppressors of eds1, we identify a recessive gain-of-function mutation in CNGC20, denoted cngc20-4, which partially restores disease resistance in eds1. cngc20-4 enhances PTI responses and ETI hypersensitive cell death. A cngc20-4 single mutant exhibits autoimmunity, which is dependent on genetically parallel EDS1 and salicylic acid (SA) pathways. CNGC20 self-associates, forms heteromeric complexes with CNGC19, and is phosphorylated and stabilized by BOTRYTIS INDUCED KINASE1 (BIK1). The cngc20-4 L371F exchange on a predicted transmembrane channel inward surface does not disrupt these interactions but leads to increased cytosolic Ca2+ accumulation, consistent with mis-regulation of CNGC20 Ca2+ -permeable channel activity. Our data show that ectopic Ca2+ influx caused by a mutant form of CNGC20 in cngc20-4 affects both PTI and ETI responses. We conclude that tight control of the CNGC20 Ca2+ ion channel is important for regulated immunity.

Natural variation identifies a Pxy gene controlling vascular organisation and formation of nodules and lateral roots in Lotus japonicus.

Forward and reverse genetics using the model legumes Lotus japonicus and Medicago truncatula have been instrumental in identifying the essential genes governing legume-rhizobia symbiosis. However, little information is known about the effects of intraspecific variation on symbiotic signalling. Here, we use quantitative trait locus sequencing (QTL-seq) to investigate the genetic basis of the differentiated phenotypic responses shown by the Lotus accessions Gifu and MG20 to inoculation with the Mesorhizobium loti exoU mutant that produces truncated exopolysaccharides. We identified through genetic complementation the Pxy gene as a component of this differential exoU response. Lotus Pxy encodes a leucine-rich repeat receptor-like kinase similar to Arabidopsis thaliana PXY, which regulates stem vascular development. We show that Lotus pxy insertion mutants displayed defects in root and stem vascular organisation, as well as lateral root and nodule formation. Our work links Pxy to de novo organogenesis in the root, highlights the genetic overlap between regulation of lateral root and nodule formation, and demonstrates that natural variation in Pxy affects nodulation signalling.

The U1 snRNP Subunit LUC7 Modulates Plant Development and Stress Responses via Regulation of Alternative Splicing.

Introns are removed by the spliceosome, a large macromolecular complex composed of five ribonucleoprotein subcomplexes (U snRNPs). The U1 snRNP, which binds to 5' splice sites, plays an essential role in early steps of the splicing reaction. Here, we show that Arabidopsis thaliana LETHAL UNLESS CBC7 (LUC7) proteins, which are encoded by a three-member gene family in Arabidopsis, are important for plant development and stress resistance. We show that LUC7 is a U1 snRNP accessory protein by RNA immunoprecipitation experiments and LUC7 protein complex purifications. Transcriptome analyses revealed that LUC7 proteins are not only important for constitutive splicing, but also affect hundreds of alternative splicing events. Interestingly, LUC7 proteins specifically promote splicing of a subset of terminal introns. Splicing of LUC7-dependent introns is a prerequisite for nuclear export, and some splicing events are modulated by stress in a LUC7-dependent manner. Taken together, our results highlight the importance of the U1 snRNP component LUC7 in splicing regulation and suggest a previously unrecognized role of a U1 snRNP accessory factor in terminal intron splicing.

Improving and correcting the contiguity of long-read genome assemblies of three plant species using optical mapping and chromosome conformation capture data.

Long-read sequencing can overcome the weaknesses of short reads in the assembly of eukaryotic genomes; however, at present additional scaffolding is needed to achieve chromosome-level assemblies. We generated Pacific Biosciences (PacBio) long-read data of the genomes of three relatives of the model plant Arabidopsis thaliana and assembled all three genomes into only a few hundred contigs. To improve the contiguities of these assemblies, we generated BioNano Genomics optical mapping and Dovetail Genomics chromosome conformation capture data for genome scaffolding. Despite their technical differences, optical mapping and chromosome conformation capture performed similarly and doubled N50 values. After improving both integration methods, assembly contiguity reached chromosome-arm-levels. We rigorously assessed the quality of contigs and scaffolds using Illumina mate-pair libraries and genetic map information. This showed that PacBio assemblies have high sequence accuracy but can contain several misassemblies, which join unlinked regions of the genome. Most, but not all, of these misjoints were removed during the integration of the optical mapping and chromosome conformation capture data. Even though none of the centromeres were fully assembled, the scaffolds revealed large parts of some centromeric regions, even including some of the heterochromatic regions, which are not present in gold standard reference sequences.

Using next-generation sequencing to isolate mutant genes from forward genetic screens.

The long-lasting success of forward genetic screens relies on the simple molecular basis of the characterized phenotypes, which are typically caused by mutations in single genes. Mapping the location of causal mutations using genetic crosses has traditionally been a complex, multistep procedure, but next-generation sequencing now allows the rapid identification of causal mutations at single-nucleotide resolution even in complex genetic backgrounds. Recent advances of this mapping-by-sequencing approach include methods that are independent of reference genome sequences, genetic crosses and any kind of linkage information, which make forward genetics amenable for species that have not been considered for forward genetic screens so far.

Ctf4-related protein recruits LHP1-PRC2 to maintain H3K27me3 levels in dividing cells in Arabidopsis thaliana.

Polycomb Repressive Complex (PRC) 2 catalyzes the H3K27me3 modification that warrants inheritance of a repressive chromatin structure during cell division, thereby assuring stable target gene repression in differentiated cells. It is still under investigation how H3K27me3 is passed on from maternal to filial strands during DNA replication; however, cell division can reinforce H3K27me3 coverage at target regions. To identify novel factors involved in the Polycomb pathway in plants, we performed a forward genetic screen for enhancers of the like heterochromatin protein 1 (lhp1) mutant, which shows relatively mild phenotypic alterations compared with other plant PRC mutants. We mapped enhancer of lhp1 (eol) 1 to a gene related to yeast Chromosome transmission fidelity 4 (Ctf4) based on phylogenetic analysis, structural similarities, physical interaction with the CMG helicase component SLD5, and an expression pattern confined to actively dividing cells. A combination of eol1 with the curly leaf (clf) allele, carrying a mutation in the catalytic core of PRC2, strongly enhanced the clf phenotype; furthermore, H3K27me3 coverage at target genes was strongly reduced in eol1 clf double mutants compared with clf single mutants. EOL1 physically interacted with CLF, its partially redundant paralog SWINGER (SWN), and LHP1. We propose that EOL1 interacts with LHP1-PRC2 complexes during replication and thereby participates in maintaining the H3K27me3 mark at target genes.

findGSE: estimating genome size variation within human and Arabidopsis using k-mer frequencies.

Motivation: Analyzing k-mer frequencies in whole-genome sequencing data is becoming a common method for estimating genome size (GS). However, it remains uninvestigated how accurate the method is, especially if it can capture intra-species GS variation. Results: We present findGSE, which fits skew normal distributions to k-mer frequencies to estimate GS. findGSE outperformed existing tools in an extensive simulation study. Estimating GSs of 89 Arabidopsis thaliana accessions, findGSE showed the highest capability in capturing GS variations. In an application with 71 female and 71 male human individuals, findGSE delivered an average of 3039 Mb as haploid human GS, while female genomes were on average 41 Mb larger than male genomes, in astonishing agreement with size difference of the X and Y chromosomes. Further analysis showed that human GS variations link to geographical patterns and significant differences between populations, which can be explained by variable abundances of LINE-1 retrotransposons. Availability and implementation: R package of findGSE is freely available at https://github.com/schneebergerlab/findGSE and supported on linux and Mac systems. Contact: schneeberger@mpipz.mpg.de. Supplementary information: Supplementary data are available at Bioinformatics online.

Strengths and potential pitfalls of hay transfer for ecological restoration revealed by RAD-seq analysis in floodplain Arabis species.

Achieving high intraspecific genetic diversity is a critical goal in ecological restoration as it increases the adaptive potential and long-term resilience of populations. Thus, we investigated genetic diversity within and between pristine sites in a fossil floodplain and compared it to sites restored by hay transfer between 1997 and 2014. RAD-seq genotyping revealed that the stenoecious floodplain species Arabis nemorensis is co-occurring with individuals that, based on ploidy, ITS-sequencing and morphology, probably belong to the close relative Arabis sagittata, which has a documented preference for dry calcareous grasslands but has not been reported in floodplain meadows. We show that hay transfer maintains genetic diversity for both species. Additionally, in A. sagittata, transfer from multiple genetically isolated pristine sites resulted in restored sites with increased diversity and admixed local genotypes. In A. nemorensis, transfer did not create novel admixture dynamics because genetic diversity between pristine sites was less differentiated. Thus, the effects of hay transfer on genetic diversity also depend on the genetic make-up of the donor communities of each species, especially when local material is mixed. Our results demonstrate the efficiency of hay transfer for habitat restoration and emphasize the importance of prerestoration characterization of microgeographic patterns of intraspecific diversity of the community to guarantee that restoration practices reach their goal, that is maximize the adaptive potential of the entire restored plant community. Overlooking these patterns may alter the balance between species in the community. Additionally, our comparison of summary statistics obtained from de novo- and reference-based RAD-seq pipelines shows that the genomic impact of restoration can be reliably monitored in species lacking prior genomic knowledge.

Complementary Activities of TELOMERE REPEAT BINDING Proteins and Polycomb Group Complexes in Transcriptional Regulation of Target Genes.

In multicellular organisms, Polycomb Repressive Complex 1 (PRC1) and PRC2 repress target genes through histone modification and chromatin compaction. Arabidopsis thaliana mutants strongly compromised in the pathway cannot develop differentiated organs. LIKE HETEROCHROMATIN PROTEIN1 (LHP1) is so far the only known plant PRC1 component that directly binds to H3K27me3, the histone modification set by PRC2, and also associates genome-wide with trimethylation of lysine 27 of histone H3 (H3K27me3). Surprisingly, lhp1 mutants show relatively mild phenotypic alterations. To explain this paradox, we screened for genetic enhancers of lhp1 mutants to identify novel components repressing target genes together with, or in parallel to, LHP1. Two enhancing mutations were mapped to TELOMERE REPEAT BINDING PROTEIN1 (TRB1) and its paralog TRB3. We show that TRB1 binds to thousands of genomic sites containing telobox or related cis-elements with a significant increase of sites and strength of binding in the lhp1 background. Furthermore, in combination with lhp1, but not alone, trb1 mutants show increased transcription of LHP1 targets, such as floral meristem identity genes, which are more likely to be bound by TRB1 in the lhp1 background. By contrast, expression of a subset of LHP1-independent TRB1 target genes, many involved in primary metabolism, is decreased in the absence of TRB1 alone. Thus, TRB1 is a bivalent transcriptional modulator that maintains downregulation of Polycomb Group (PcG) target genes in lhp1 mutants, while it sustains high expression of targets that are regulated independently of PcG.