SEPALLATA-driven MADS transcription factor tetramerization is required for inner whorl floral organ development.

MADS transcription factors are master regulators of plant reproduction and flower development. The SEPALLATA (SEP) subfamily of MADS transcription factors is required for the development of floral organs and plays roles in inflorescence architecture and development of the floral meristem. SEPALLATAs act as organizers of MADS complexes, forming both heterodimers and heterotetramers in vitro. To date, the MADS complexes characterized in angiosperm floral organ development contain at least one SEPALLATA protein. Whether DNA-binding by SEPALLATA-containing dimeric MADS complexes is sufficient for launching floral organ identity programs, however, is not clear as only defects in floral meristem determinacy were observed in tetramerization--impaired SEPALLATA mutant proteins. Here, we used a combination of genome-wide binding studies, high resolution structural studies of the SEP3/AGAMOUS (AG) tetramerization domain, structure-based mutagenesis and complementation experiments in Arabidopsis (Arabidopsis thaliana) sep1 sep2 sep3 and sep1 sep2 sep3 ag-4 plants transformed with versions of SEP3 encoding tetramerization mutants. We demonstrate that while SEP3 heterodimers can bind DNA both in vitro and in vivo and recognize the majority of SEP3 wild-type binding sites genome-wide, tetramerization is required not only for floral meristem determinacy, but also for floral organ identity in the second, third and fourth whorls.

Plant-TFClass: a structural classification for plant transcription factors.

Transcription factors (TFs) bind DNA at specific sequences to regulate gene expression. This universal process is achieved via their DNA-binding domain (DBD). In mammals, the vast diversity of DBD structural conformations and the way in which they contact DNA has been used to organize TFs in the TFClass hierarchical classification. However, the numerous DBD types present in plants but absent from mammalian genomes were missing from this classification. We reviewed DBD 3D structures and models available for plant TFs to classify most of the 56 recognized plant TF types within the TFClass framework. This extended classification adds eight new classes and 37 new families corresponding to DBD structures absent in mammals. Plant-TFClass provides a unique resource for TF comparison across families and organisms.

JASPAR 2024: 20th anniversary of the open-access database of transcription factor binding profiles.

JASPAR (https://jaspar.elixir.no/) is a widely-used open-access database presenting manually curated high-quality and non-redundant DNA-binding profiles for transcription factors (TFs) across taxa. In this 10th release and 20th-anniversary update, the CORE collection has expanded with 329 new profiles. We updated three existing profiles and provided orthogonal support for 72 profiles from the previous release's UNVALIDATED collection. Altogether, the JASPAR 2024 update provides a 20% increase in CORE profiles from the previous release. A trimming algorithm enhanced profiles by removing low information content flanking base pairs, which were likely uninformative (within the capacity of the PFM models) for TFBS predictions and modelling TF-DNA interactions. This release includes enhanced metadata, featuring a refined classification for plant TFs' structural DNA-binding domains. The new JASPAR collections prompt updates to the genomic tracks of predicted TF binding sites (TFBSs) in 8 organisms, with human and mouse tracks available as native tracks in the UCSC Genome browser. All data are available through the JASPAR web interface and programmatically through its API and the updated Bioconductor and pyJASPAR packages. Finally, a new TFBS extraction tool enables users to retrieve predicted JASPAR TFBSs intersecting their genomic regions of interest.

Identification of Plant Transcription Factor DNA-Binding Sites Using seq-DAP-seq.

The identification of genome-wide transcription factor binding sites (TFBS) is a critical step in deciphering gene and transcriptional regulatory networks. However, determining the genome-wide binding of specific TFs or TF complexes remains a technical challenge. DNA affinity purification sequencing (DAP-seq) and modifications such as sequential DAP-seq (seq-DAP-seq) are robust in vitro methods for mapping individual TF or TF complex binding sites in a genome-wide manner. DAP-seq protocols use a genomic DNA (gDNA) library from any target organism with or without amplification, allowing the determination of TF binding on naked or endogenously modified DNA, respectively. As a first step, the gDNA is fragmented to ~200 bp, end-repaired, and sequencing adaptors are added. This gDNA library can be used directly or an amplification step may be performed to remove DNA modifications such as cytosine methylation. DNA libraries are then incubated with an affinity-tagged TF or TF- complex immobilized on magnetic beads. The TF or TF complex of interest is generally produced using recombinant protein expression and purified prior to DNA affinity purification. After incubation of the DNA library with the immobilized TF of interest, multiple wash steps are performed to reduce non-specific DNA binding and the TF-DNA complexes eluted. The eluted DNA is PCR-amplified and sequenced using next-generation sequencing. The resulting sequence reads are mapped to the corresponding reference genome, identifying direct potential bound regions and binding sites of the TF or TF complex of interest. Predictive TFBS models are generated from the bound regions using downstream bioinformatics analysis pipelines. Here, we present a detailed protocol outlining the steps required for seq-DAP-seq of a heterooligomeric TF complex (Fig. 1) and briefly describe the downstream bioinformatics pipeline used to develop a robust TFBS model from sequencing data generated from a DAP-seq experiment.

Genome analysis of Parmales, the sister group of diatoms, reveals the evolutionary specialization of diatoms from phago-mixotrophs to photoautotrophs.

The order Parmales (class Bolidophyceae) is a minor group of pico-sized eukaryotic marine phytoplankton that contains species with cells surrounded by silica plates. Previous studies revealed that Parmales is a member of ochrophytes and sister to diatoms (phylum Bacillariophyta), the most successful phytoplankton group in the modern ocean. Therefore, parmalean genomes can serve as a reference to elucidate both the evolutionary events that differentiated these two lineages and the genomic basis for the ecological success of diatoms vs. the more cryptic lifestyle of parmaleans. Here, we compare the genomes of eight parmaleans and five diatoms to explore their physiological and evolutionary differences. Parmaleans are predicted to be phago-mixotrophs. By contrast, diatoms have lost genes related to phagocytosis, indicating the ecological specialization from phago-mixotrophy to photoautotrophy in their early evolution. Furthermore, diatoms show significant enrichment in gene sets involved in nutrient uptake and metabolism, including iron and silica, in comparison with parmaleans. Overall, our results suggest a strong evolutionary link between the loss of phago-mixotrophy and specialization to a silicified photoautotrophic life stage early in diatom evolution after diverging from the Parmales lineage.

Single-cytosine methylation at W-boxes repels binding of WRKY transcription factors through steric hindrance.

DNA methylation is an epigenetic mark that fine-tunes gene expression, notably by negatively or positively regulating transcription factor (TF)-DNA binding. In plants, DNA methylation has primarily been shown to inhibit TF-DNA binding. However, little is known about the underlying mechanisms. Here, we show that DNA methylation decreases the binding of several Arabidopsis (Arabidopsis thaliana) WRKY TFs to their genomic regions and their binding sites in vitro. We also provide evidence that DNA methylation at a single cytosine located in a functional core W-box motif repels DNA binding of AtWRKY40 in vitro. Using structural modelling, we further demonstrate that this cytosine interacts through van der Waals contacts with the conserved tyrosine of WRKY-DNA binding domains. Importantly, our model predicts steric hindrance when a 5-methyl group is present on this specific cytosine, thereby likely preventing tight binding of WRKY-DNA binding domains. Finally, because the WRKY motif and the residues involved in DNA contacts are conserved across Arabidopsis and rice (Oryza sativa) WRKY TFs, we propose that this methylation-dependent WRKY-DNA binding inhibitory mechanism could be widespread across plant species.

The F-box protein UFO controls flower development by redirecting the master transcription factor LEAFY to new cis-elements.

In angiosperms, flower development requires the combined action of the transcription factor LEAFY (LFY) and the ubiquitin ligase adaptor F-box protein, UNUSUAL FLORAL ORGANS (UFO), but the molecular mechanism underlying this synergy has remained unknown. Here we show in transient assays and stable transgenic plants that the connection to ubiquitination pathways suggested by the UFO F-box domain is mostly dispensable. On the basis of biochemical and genome-wide studies, we establish that UFO instead acts by forming an active transcriptional complex with LFY at newly discovered regulatory elements. Structural characterization of the LFY-UFO-DNA complex by cryo-electron microscopy further demonstrates that UFO performs this function by directly interacting with both LFY and DNA. Finally, we propose that this complex might have a deep evolutionary origin, largely predating flowering plants. This work reveals a unique mechanism of an F-box protein directly modulating the DNA binding specificity of a master transcription factor.

Functional Differentiation among the Arabidopsis Phosphatidylinositol 4-Phosphate 5-Kinase Genes PIP5K1, PIP5K2 and PIP5K3.

Phosphatidylinositol 4-phosphate 5-kinase (PIP5K) is involved in regulating various cellular processes through the signaling function of its product, phosphatidylinositol (4,5)-bisphosphate. Higher plants encode a large number of PIP5Ks forming distinct clades in their molecular phylogenetic tree. Although biological functions of PIP5K genes have been analyzed intensively in Arabidopsis thaliana, it remains unclear how those functions differ across clades of paralogs. We performed comparative functional analysis of the Arabidopsis genes encoding PIP5K1, PIP5K2 and PIP5K3, of which the first two and the last belong to closely related but distinct clades, to clarify their conserved and/or differentiated functions. Genetic analysis with their single and multiple mutants revealed that PIP5K1 and PIP5K3 have non-overlapping functions, with the former in total plant growth and the latter in root hair elongation, whereas PIP5K2 redundantly functions in both phenomena. This pattern of functional redundancy is explainable in terms of the overlapping pattern of their promoter activities. In transformation rescue experiments, PIP5K3 promoter-directed PIP5K1-YFP completely rescued the short-root-hair phenotype of pip5k3. However, PIP5K3-YFP could substitute for PIP5K1-YFP only partially in rescuing the severe dwarfism of pip5k1pip5k2 when directed by the PIP5K1 promoter. Phylogenetic analysis of angiosperm PIP5Ks revealed that PIP5K3 orthologs have a faster rate of diversification in their amino-acid sequences compared with PIP5K1/2 orthologs after they arose through a eudicot-specific duplication event. These findings suggest that PIP5K3 specialized to promote root hair elongation and lost some of the protein-encoded functions retained by PIP5K1 and PIP5K2, whereas PIP5K1 differentiated from PIP5K2 only in its promoter-directed expression pattern.

Adaptive evolution of viruses infecting marine microalgae (haptophytes), from acute infections to stable coexistence.

Collectively known as phytoplankton, photosynthetic microbes form the base of the marine food web, and account for up to half of the primary production on Earth. Haptophytes are key components of this phytoplankton community, playing important roles both as primary producers and as mixotrophs that graze on bacteria and protists. Viruses influence the ecology and diversity of phytoplankton in the ocean, with the majority of microalgae-virus interactions described as 'boom and bust' dynamics, which are characteristic of acute virus-host systems. Most haptophytes are, however, part of highly diverse communities and occur at low densities, decreasing their chance of being infected by viruses with high host specificity. Viruses infecting these microalgae have been isolated in the laboratory, and there are several characteristics that distinguish them from acute viruses infecting bloom-forming haptophytes. Herein we synthesise what is known of viruses infecting haptophyte hosts in the ocean, discuss the adaptive evolution of haptophyte-infecting viruses -from those that cause acute infections to those that stably coexist with their host - and identify traits of importance for successful survival in the ocean.

Year-round dynamics of amplicon sequence variant communities differ among eukaryotes, Imitervirales and prokaryotes in a coastal ecosystem.

Coastal microbial communities are affected by seasonal environmental change, biotic interactions and fluctuating nutrient availability. We investigated the seasonal dynamics of communities of eukaryotes, a major group of double-stranded DNA viruses that infect eukaryotes (order Imitervirales; phylum Nucleocytoviricota), and prokaryotes in the Uranouchi Inlet, Kochi, Japan. We performed metabarcoding using ribosomal RNA genes and viral polB genes as markers in 43 seawater samples collected over 20 months. Eukaryotes, prokaryotes and Imitervirales communities characterized by the compositions of amplicon sequence variants (ASVs) showed synchronic seasonal cycles. However, the community dynamics showed intriguing differences in several aspects, such as the recovery rate after a year. We also showed that the differences in community dynamics were at least partially explained by differences in recurrence/persistence levels of individual ASVs among eukaryotes, prokaryotes and Imitervirales. Prokaryotic ASVs were the most persistent, followed by eukaryotic ASVs and Imitervirales ASVs, which were the least persistent. We argue that the differences in the specificity of interactions (virus-eukaryote vs prokaryote-eukaryote) as well as the niche breadth of community members were at the origin of the distinct community dynamics among eukaryotes, their viruses and prokaryotes.

JASPAR 2022: the 9th release of the open-access database of transcription factor binding profiles.

JASPAR (http://jaspar.genereg.net/) is an open-access database containing manually curated, non-redundant transcription factor (TF) binding profiles for TFs across six taxonomic groups. In this 9th release, we expanded the CORE collection with 341 new profiles (148 for plants, 101 for vertebrates, 85 for urochordates, and 7 for insects), which corresponds to a 19% expansion over the previous release. We added 298 new profiles to the Unvalidated collection when no orthogonal evidence was found in the literature. All the profiles were clustered to provide familial binding profiles for each taxonomic group. Moreover, we revised the structural classification of DNA binding domains to consider plant-specific TFs. This release introduces word clouds to represent the scientific knowledge associated with each TF. We updated the genome tracks of TFBSs predicted with JASPAR profiles in eight organisms; the human and mouse TFBS predictions can be visualized as native tracks in the UCSC Genome Browser. Finally, we provide a new tool to perform JASPAR TFBS enrichment analysis in user-provided genomic regions. All the data is accessible through the JASPAR website, its associated RESTful API, the R/Bioconductor data package, and a new Python package, pyJASPAR, that facilitates serverless access to the data.

The intervening domain is required for DNA-binding and functional identity of plant MADS transcription factors.

The MADS transcription factors (TF) are an ancient eukaryotic protein family. In plants, the family is divided into two main lineages. Here, we demonstrate that DNA binding in both lineages absolutely requires a short amino acid sequence C-terminal to the MADS domain (M domain) called the Intervening domain (I domain) that was previously defined only in type II lineage MADS. Structural elucidation of the MI domains from the floral regulator, SEPALLATA3 (SEP3), shows a conserved fold with the I domain acting to stabilise the M domain. Using the floral organ identity MADS TFs, SEP3, APETALA1 (AP1) and AGAMOUS (AG), domain swapping demonstrate that the I domain alters genome-wide DNA-binding specificity and dimerisation specificity. Introducing AG carrying the I domain of AP1 in the Arabidopsis ap1 mutant resulted in strong complementation and restoration of first and second whorl organs. Taken together, these data demonstrate that the I domain acts as an integral part of the DNA-binding domain and significantly contributes to the functional identity of the MADS TF.

Quantitative Assessment of Nucleocytoplasmic Large DNA Virus and Host Interactions Predicted by Co-occurrence Analyses.

Nucleocytoplasmic large DNA viruses (NCLDVs) are highly diverse and abundant in marine environments. However, the knowledge of their hosts is limited because only a few NCLDVs have been isolated so far. Taking advantage of the recent large-scale marine metagenomics census, in silico host prediction approaches are expected to fill the gap and further expand our knowledge of virus-host relationships for unknown NCLDVs. In this study, we built co-occurrence networks of NCLDVs and eukaryotic taxa to predict virus-host interactions using Tara Oceans sequencing data. Using the positive likelihood ratio to assess the performance of host prediction for NCLDVs, we benchmarked several co-occurrence approaches and demonstrated an increase in the odds ratio of predicting true positive relationships 4-fold compared to random host predictions. To further refine host predictions from high-dimensional co-occurrence networks, we developed a phylogeny-informed filtering method, Taxon Interaction Mapper, and showed it further improved the prediction performance by 12-fold. Finally, we inferred virophage-NCLDV networks to corroborate that co-occurrence approaches are effective for predicting interacting partners of NCLDVs in marine environments.IMPORTANCE NCLDVs can infect a wide range of eukaryotes, although their life cycle is less dependent on hosts compared to other viruses. However, our understanding of NCLDV-host systems is highly limited because few of these viruses have been isolated so far. Co-occurrence information has been assumed to be useful to predict virus-host interactions. In this study, we quantitatively show the effectiveness of co-occurrence inference for NCLDV host prediction. We also improve the prediction performance with a phylogeny-guided method, which leads to a concise list of candidate host lineages for three NCLDV families. Our results underpin the usage of co-occurrence approaches for the metagenomic exploration of the ecology of this diverse group of viruses.

The LEAFY floral regulator displays pioneer transcription factor properties.

Pioneer transcription factors (TFs) are a special category of TFs with the capacity to bind to closed chromatin regions in which DNA is wrapped around histones and may be highly methylated. Subsequently, pioneer TFs are able to modify the chromatin state to initiate gene expression. In plants, LEAFY (LFY) is a master floral regulator and has been suggested to act as a pioneer TF in Arabidopsis. Here, we demonstrate that LFY is able to bind both methylated and non-methylated DNA using a combination of in vitro genome-wide binding experiments and structural modeling. Comparisons between regions bound by LFY in vivo and chromatin accessibility data suggest that a subset of LFY bound regions is occupied by nucleosomes. We confirm that LFY is able to bind nucleosomal DNA in vitro using reconstituted nucleosomes. Finally, we show that constitutive LFY expression in seedling tissues is sufficient to induce chromatin accessibility in the LFY direct target genes APETALA1 and AGAMOUS. Taken together, our study suggests that LFY possesses key pioneer TF features that contribute to launching the floral gene expression program.

A persistent giant algal virus, with a unique morphology, encodes an unprecedented number of genes involved in energy metabolism.

Viruses have long been viewed as entities possessing extremely limited metabolic capacities. Over the last decade, however, this view has been challenged, as metabolic genes have been identified in viruses possessing large genomes and virions-the synthesis of which is energetically demanding. Here, we unveil peculiar phenotypic and genomic features of Prymnesium kappa virus RF01 (PkV RF01), a giant virus of the Mimiviridae family. We found that this virus encodes an unprecedented number of proteins involved in energy metabolism, such as all four succinate dehydrogenase (SDH) subunits (A-D) as well as key enzymes in the ß-oxidation pathway. The SDHA gene was transcribed upon infection, indicating that the viral SDH is actively used by the virus- potentially to modulate its host's energy metabolism. We detected orthologous SDHA and SDHB genes in numerous genome fragments from uncultivated marine Mimiviridae viruses, which suggests that the viral SDH is widespread in oceans. PkV RF01 was less virulent compared with other cultured prymnesioviruses, a phenomenon possibly linked to the metabolic capacity of this virus and suggestive of relatively long co-evolution with its hosts. It also has a unique morphology, compared to other characterized viruses in the Mimiviridae family. Finally, we found that PkV RF01 is the only alga-infecting Mimiviridae virus encoding two aminoacyl-tRNA synthetases and enzymes corresponding to an entire base-excision repair pathway, as seen in heterotroph-infecting Mimiviridae These Mimiviridae encoded-enzymes were found to be monophyletic and branching at the root of the eukaryotic tree of life. This placement suggests that the last common ancestor of Mimiviridae was endowed with a large, complex genome prior to the divergence of known extant eukaryotes.IMPORTANCE Viruses on Earth are tremendously diverse in terms of morphology, functionality, and genomic composition. Over the last decade, the conceptual gap separating viruses and cellular life has tightened because of the detection of metabolic genes in viral genomes that express complex virus phenotypes upon infection. Here, we describe Prymnesium kappa virus RF01, a large alga-infecting virus with a unique morphology, an atypical infection profile, and an unprecedented number of genes involved in energy metabolism (such as the tricarboxylic (TCA) cycle and the ß-oxidation pathway). Moreover, we show that the gene corresponding to one of these enzymes (the succinate dehydrogenase subunit A) is transcribed during infection and is widespread among marine viruses. This discovery provides evidence that a virus has the potential to actively regulate energy metabolism with its own gene.

Eukaryotic virus composition can predict the efficiency of carbon export in the global ocean.

The biological carbon pump, in which carbon fixed by photosynthesis is exported to the deep ocean through sinking, is a major process in Earth's carbon cycle. The proportion of primary production that is exported is termed the carbon export efficiency (CEE). Based on in-lab or regional scale observations, viruses were previously suggested to affect the CEE (i.e., viral "shunt" and "shuttle"). In this study, we tested associations between viral community composition and CEE measured at a global scale. A regression model based on relative abundance of viral marker genes explained 67% of the variation in CEE. Viruses with high importance in the model were predicted to infect ecologically important hosts. These results are consistent with the view that the viral shunt and shuttle functions at a large scale and further imply that viruses likely act in this process in a way dependent on their hosts and ecosystem dynamics.

Ecological Structuring of Temperate Bacteriophages in the Inflammatory Bowel Disease-Affected Gut.

The aim of this study was to elucidate the ecological structure of the human gut temperate bacteriophage community and its role in inflammatory bowel disease (IBD). Temperate bacteriophages make up a large proportion of the human gut microbiota and are likely to play a role in IBD pathogenesis. However, many of these bacteriophages await characterization in reference databases. Therefore, we conducted a large-scale reconstruction of temperate bacteriophage and bacterial genomes from the whole-metagenome sequence data generated by the IBD Multi'omics Database project. By associating phages with their hosts via genome comparisons, we found that temperate bacteriophages infect a phylogenetically wide range of bacteria. The majority of variance in bacteriophage community composition was explained by variation among individuals, but differences in the abundance of temperate bacteriophages were identified between IBD and non-IBD patients. Of note, in active ulcerative colitis patients, temperate bacteriophages infecting Bacteroides uniformis and Bacteroides thetaiotaomicron-two species experimentally proven to be beneficial to gut homeostasis-were over-represented, whereas their hosts were under-represented in comparison with non-IBD patients. Supporting the mounting evidence that gut viral community plays a vital role in IBD, our results show potential association between temperate bacteriophages and IBD pathogenesis.

Biogeography of marine giant viruses reveals their interplay with eukaryotes and ecological functions.

Nucleocytoplasmic large DNA viruses (NCLDVs) are ubiquitous in marine environments and infect diverse eukaryotes. However, little is known about their biogeography and ecology in the ocean. By leveraging the Tara Oceans pole-to-pole metagenomic data set, we investigated the distribution of NCLDVs across size fractions, depths and biomes, as well as their associations with eukaryotic communities. Our analyses reveal a heterogeneous distribution of NCLDVs across oceans, and a higher proportion of unique NCLDVs in the polar biomes. The community structures of NCLDV families correlate with specific eukaryotic lineages, including many photosynthetic groups. NCLDV communities are generally distinct between surface and mesopelagic zones, but at some locations they exhibit a high similarity between the two depths. This vertical similarity correlates to surface phytoplankton biomass but not to physical mixing processes, which suggests a potential role of vertical transport in structuring mesopelagic NCLDV communities. These results underscore the importance of the interactions between NCLDVs and eukaryotes in biogeochemical processes in the ocean.

An Optimized Metabarcoding Method for Mimiviridae.

Mimiviridae is a group of viruses with large genomes and virions. Ecological relevance of Mimiviridae in marine environments has been increasingly recognized through the discoveries of novel isolates and metagenomic studies. To facilitate ecological profiling of Mimiviridae, we previously proposed a meta-barcoding approach based on 82 degenerate primer pairs (i.e., MEGAPRIMER) targeting the DNA polymerase gene of Mimiviridae. The method detected a larger number of operational taxonomic units (OTUs) in environmental samples than previous methods. However, it required large quantities of DNA and was laborious due to the use of individual primer pairs. Here, we examined coastal seawater samples using varying PCR conditions and purification protocols to streamline the MEGAPRIMER method. Mixing primer pairs in "cocktails" reduced the required amount of environmental DNA by 90%, while reproducing the results obtained by the original protocol. We compared the results obtained by the meta-barcoding approach with quantifications using qPCR for selected OTUs. This revealed possible amplification biases among different OTUs, but the frequency profiles for individual OTUs across multiple samples were similar to those obtained by qPCR. We anticipate that the newly developed MEGAPRIMER protocols will be useful for ecological investigation of Mimiviridae in a larger set of environmental samples.

KofamKOALA: KEGG Ortholog assignment based on profile HMM and adaptive score threshold.

SUMMARY: KofamKOALA is a web server to assign KEGG Orthologs (KOs) to protein sequences by homology search against a database of profile hidden Markov models (KOfam) with pre-computed adaptive score thresholds. KofamKOALA is faster than existing KO assignment tools with its accuracy being comparable to the best performing tools. Function annotation by KofamKOALA helps linking genes to KEGG resources such as the KEGG pathway maps and facilitates molecular network reconstruction. AVAILABILITY AND IMPLEMENTATION: KofamKOALA, KofamScan and KOfam are freely available from GenomeNet (https://www.genome.jp/tools/kofamkoala/). SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.