Transcription factor evolution in eukaryotes and the assembly of the regulatory toolkit in multicellular lineages

Recent advances in improving yield and immunity through transcription factor engineering

Transcription factors (TFs) function as master regulators in multiple signaling pathways and govern diverse developmental and adaptive processes in plants. Some TFs identified in crop plants play critical roles in regulating yield through changes in plant architecture, including roots, stems, leaves, flowers, fruits, and grains. Although altering crop architecture can increase yields, the extent of yield enhancement is frequently hampered by diseases. Developing new crop varieties with improved yields and enhanced disease resistance remains challenging because immune system activation often impairs plant growth. Recently, approaches using TF engineering have made substantial progress in elevating both growth performance and disease resistance. However, most of these techniques rely on traditional transgenic methods. This review highlights discoveries in the last decade demonstrating improvements in growth performance, yield and immunity through TF engineering. We focus mainly on changes in plant architecture related to improved yield and disease resistance. We conclude with perspectives on the potential application of these discoveries for generating desirable crop traits by merging the most noteworthy biotechnology approaches, such as clustered regularly interspaced small palindromic repeats (CRISPR)/CRISPR-associated protein 9-mediated genome editing, with canonical molecular biology.

The expanded Bostrychia moritziana genome unveils evolution in the most diverse and complex order of red algae

Red algae are an ancient eukaryotic lineage that were among the first to evolve multicellularity. Although they share a common origin with modern-day plants and display complex multicellular development, comprehensive genome data from the most highly evolved red algal groups remain scarce. Here, we present a chromosome-level genome assembly of Bostrychia moritziana, a complex red seaweed in the Rhodomelaceae family of the Ceramiales-the largest and most diverse order of red algae. Contrary to the view that red algal genomes are typically small, we report significant genome size expansion in Bostrychia and other Ceramiales, which represents one of at least three independent expansion events in red algal evolution. Our analyses suggest that these expansions do not involve polyploidy or ancient whole-genome duplications, but in Bostrychia rather stem from the proliferation of a single lineage of giant Plavaka DNA transposons. Consistent with its enlarged genome, Bostrychia has an increased gene content shaped by de novo gene emergence and amplified gene families in common with other Ceramiales, providing insight into the genetic adaptations underpinning this successful and species-rich order. Finally, our sex-specific assemblies resolve the UV sex chromosomes in Bostrychia, which feature expanded gene-rich sex-linked regions. Notably, each sex chromosome harbors a three amino acid loop extension homeodomain (TALE-HD) transcription factor orthologous to ancient regulators of haploid-diploid transitions in other multicellular lineages. Together, our findings offer a unique perspective of the genomic adaptations driving red algal diversity and demonstrate how this red seaweed lineage can provide insight into the evolutionary origins and universal principles underpinning complex multicellularity.

Evolution of sexual systems and regressive evolution in Riccia

The genus Riccia is the most speciose of the Marchantiopsida, the complex thalloid liverworts. In contrast to the predicted ancestral liverwort and the majority of extant liverworts that are dioicous, most Riccia species are monoicous. Both the gametophyte and the sporophyte generations have undergone regressive evolution in Riccia, with the sporophyte lacking many adaptations for spore dispersal. Phylogenetic reconstructions suggest multiple transitions between dioicy and monoicy, and vice versa, within the Riccia. We analysed the genomes of two Riccia species and the transcriptomes of two additional Riccia species to explore the evolution of sexual systems and investigate genomic signatures of regressive evolution. Genomic analyses suggest a transition from dioicy to monoicy at the base of the Riccia clade and a subsequent reversion to dioicy with the concomitant evolution of a new sex chromosome, followed by further transitions back to monoicy. Additionally, Riccia species exhibit significant losses of transcription factors known to control aspects of sporophyte and gametophyte development. Comparisons with the monotypic sister genus, Ricciocarpos, suggest independent evolutions to both monoicy and reductive morphologies in the two genera. The reductive morphology has facilitated their colonization of extreme habitats but has perhaps canalized them to remain there.

Development necessitates evolutionarily conserved factors

Early-stage generalised transcription factors in biological development are often evolutionarily conserved across species. Here, we find for the first time that similar factors functionally emerge in an alternative medium of development. Through comprehensively analysing a Neural Cellular Automata (NCA) model of morphogenesis, we find multiple properties of the hidden units that are functionally analogous to early factors in biological development. We test the generalisation abilities of our model through transfer learning of other morphologies and find that developmental strategies learnt by the model are reused to grow new body forms by conserving its early generalised factors. Our paper therefore provides evidence that nature did not become locked into one arbitrary method of developing multicellular organisms: the use of early generalised factors as fundamental control mechanisms and the resulting necessity for evolutionary conservation of those factors may be fundamental to development, regardless of the details of how development is implemented.

The OsZHD1 and OsZHD2, Two Zinc Finger Homeobox Transcription Factor, Redundantly Control Grain Size by Influencing Cell Proliferation in Rice

Grain size is vital determinant for grain yield and quality, which specified by its three-dimensional structure of seeds (length, width and thickness). The ZINC FINGER-HOMEODOMAIN (ZHD) proteins play critical roles in plant growth and development. However, the information regarding the function in reproductive development of ZHD proteins is scarce. Here, we deeply characterized the phenotype of oszhd1, oszhd2, and oszhd1oszhd2. The single mutants of OsZHD1/2 were similar with wild type. Nevertheless, the double mutant displayed dwarfism and smaller reproductive organs, and shorter, narrower, and thinner grain size. oszhd1oszhd2 revealed a significant decrease in total cell length and number, and single cell width in outer parenchyma; reducing the average width of longitudinal epidermal cells, but the length were increased in outer and inner glumes of oszhd1oszhd2 compared with wild-type, oszhd1-1, oszhd2-1, respectively. OsZHD1 and OsZHD2 encoded the nucleus protein and were distributed predominately in stem and the developing spikelets, asserting their roles in grain size. Meanwhile, yeast two-hybrid, bimolecular fluorescence complementation, and Co-immunoprecipitation assay clarified that OsZHD1 could directly interacted with OsZHD2. The differential expression analysis showed that 839 DEGs, which were down-regulated in oszhd1oszhd2 than wild type and single mutants, were mainly enriched in secondary metabolite biosynthetic, integral component of membrane, and transporter activity pathway. Moreover, it is reliable that the altered expression of cell cycle and expansion-related and grain size-related genes were observed in RNA-seq data, highly consistent with the qRT-PCR results. Altogether, our results suggest that OsZHD1/2 are functional redundancy and involved in regulating grain size by influencing cell proliferation in rice.

Evolution of Transcription Factor-containing Superfamilies in Eukaryotes

Regulation of gene expression helps determine various phenotypes in most cellular life forms. It is orchestrated at different levels and at the point of transcription initiation by transcription factors (TFs). TFs bind to DNA through domains that are evolutionarily related, by shared membership of the same superfamilies (TF-SFs), to those found in other nucleic acid binding and protein-binding functions (nTFs for non-TFs). Here we ask how TF DNA binding sequence families in eukaryotes have evolved in relation to their nTF relatives. TF numbers scale by power law with the total number of protein-coding genes differently in different clades, with fungi usually showing sub-linear powers whereas chordates show super-linear scaling. The LECA probably encoded a complex regulatory machinery with both TFs and nTFs, but with an excess of nTFs when compared to the relative distribution of TFs and nTFs in extant organisms. Losses drive the evolution of TFs and nTFs, with the possible exception of TFs in animals for some tree topologies. TFs are highly dynamic in evolution, showing higher gain and loss rates than nTFs in some TF-SFs though both are conserved to similar extents. Gains of TFs and nTFs are driven by the appearance of a large number of new sequence clusters in a small number of nodes, which determine the presence of as many as a third of extant TFs and nTFs as well as the relative presence of TFs and nTFs. Whereas nodes showing explosion of TF numbers belong to multicellular clades, those for nTFs lie among the fungi and the protists.

Molecular and Functional Convergences Associated with Complex Multicellularity in Eukarya

A key trait of Eukarya is the independent evolution of complex multicellularity in animals, land plants, fungi, brown algae, and red algae. This phenotype is characterized by the initial exaptation of cell-cell adhesion genes followed by the emergence of mechanisms for cell-cell communication, together with the expansion of transcription factor gene families responsible for cell and tissue identity. The number of cell types is commonly used as a quantitative proxy for biological complexity in comparative genomics studies. While expansions of individual gene families have been associated with variations in the number of cell types within individual complex multicellular lineages, the molecular and functional roles responsible for the independent evolution of complex multicellular across Eukarya remain poorly understood. We employed a phylogeny-aware strategy to conduct a genomic-scale search for associations between the number of cell types and the abundance of genomic components across a phylogenetically diverse set of 81 eukaryotic species, including species from all complex multicellular lineages. Our annotation schemas represent 2 complimentary aspects of genomic information: homology, represented by conserved sequences, and function, represented by Gene Ontology terms. We found many gene families sharing common biological themes that define complex multicellular to be independently expanded in 2 or more complex multicellular lineages, such as components of the extracellular matrix, cell-cell communication mechanisms, and developmental pathways. Additionally, we describe many previously unknown associations of biological themes and biological complexity, such as expansions of genes playing roles in wound response, immunity, cell migration, regulatory processes, and response to natural rhythms. Together, our findings unveil a set of functional and molecular convergences independently expanded in complex multicellular lineages likely due to the common selective pressures in their lifestyles.

Evolutionary trajectory of transcription factors and selection of targets for metabolic engineering

Transcription factors (TFs) provide potentially powerful tools for plant metabolic engineering as they often control multiple genes in a metabolic pathway. However, selecting the best TF for a particular pathway has been challenging, and the selection often relies significantly on phylogenetic relationships. Here, we offer examples where evolutionary relationships have facilitated the selection of the suitable TFs, alongside situations where such relationships are misleading from the perspective of metabolic engineering. We argue that the evolutionary trajectory of a particular TF might be a better indicator than protein sequence homology alone in helping decide the best targets for plant metabolic engineering efforts. This article is part of the theme issue 'The evolution of plant metabolism'.

The emergence of Sox and POU transcription factors predates the origins of animal stem cells

Stem cells are a hallmark of animal multicellularity. Sox and POU transcription factors are associated with stemness and were believed to be animal innovations, reported absent in their unicellular relatives. Here we describe unicellular Sox and POU factors. Choanoflagellate and filasterean Sox proteins have DNA-binding specificity similar to mammalian Sox2. Choanoflagellate-but not filasterean-Sox can replace Sox2 to reprogram mouse somatic cells into induced pluripotent stem cells (iPSCs) through interacting with the mouse POU member Oct4. In contrast, choanoflagellate POU has a distinct DNA-binding profile and cannot generate iPSCs. Ancestrally reconstructed Sox proteins indicate that iPSC formation capacity is pervasive among resurrected sequences, thus loss of Sox2-like properties fostered Sox family subfunctionalization. Our findings imply that the evolution of animal stem cells might have involved the exaptation of a pre-existing set of transcription factors, where pre-animal Sox was biochemically similar to extant Sox, whilst POU factors required evolutionary innovations.

START domains generate paralog-specific regulons from a single network architecture

Functional divergence of transcription factors (TFs) has driven cellular and organismal complexity throughout evolution, but its mechanistic drivers remain poorly understood. Here we test for new mechanisms using CORONA (CNA) and PHABULOSA (PHB), two functionally diverged paralogs in the CLASS III HOMEODOMAIN LEUCINE ZIPPER (HD-ZIPIII) family of TFs. We show that virtually all genes bound by PHB ( ~ 99%) are also bound by CNA, ruling out occupation of distinct sets of genes as a mechanism of functional divergence. Further, genes bound and regulated by both paralogs are almost always regulated in the same direction, ruling out opposite regulation of shared targets as a mechanistic driver. Functional divergence of CNA and PHB instead results from differential usage of shared binding sites, with hundreds of uniquely regulated genes emerging from a commonly bound genetic network. Regulation of a given gene by CNA or PHB is thus a function of whether a bound site is considered 'responsive' versus 'non-responsive' by each paralog. Discrimination between responsive and non-responsive sites is controlled, at least in part, by their lipid binding START domain. This suggests a model in which HD-ZIPIII TFs use information integrated by their START domain to generate paralog-specific transcriptional outcomes from a shared network architecture. Taken together, our study identifies a mechanism of HD-ZIPIII TF paralog divergence and proposes the ubiquitously distributed START evolutionary module as a driver of functional divergence.

Reconstructing the last common ancestor of all eukaryotes

Understanding the origin of eukaryotic cells is one of the most difficult problems in all of biology. A key challenge relevant to the question of eukaryogenesis is reconstructing the gene repertoire of the last eukaryotic common ancestor (LECA). As data sets grow, sketching an accurate genomics-informed picture of early eukaryotic cellular complexity requires provision of analytical resources and a commitment to data sharing. Here, we summarise progress towards understanding the biology of LECA and outline a community approach to inferring its wider gene repertoire. Once assembled, a robust LECA gene set will be a useful tool for evaluating alternative hypotheses about the origin of eukaryotes and understanding the evolution of traits in all descendant lineages, with relevance in diverse fields such as cell biology, microbial ecology, biotechnology, agriculture, and medicine. In this Consensus View, we put forth the status quo and an agreed path forward to reconstruct LECA's gene content.

A multicellular developmental program in a close animal relative

All animals develop from a single-celled zygote into a complex multicellular organism through a series of precisely orchestrated processes1,2. Despite the remarkable conservation of early embryogenesis across animals, the evolutionary origins of how and when this process first emerged remain elusive. Here, by combining time-resolved imaging and transcriptomic profiling, we show that single cells of the ichthyosporean Chromosphaera perkinsii-a close relative that diverged from animals about 1 billion years ago3,4-undergo symmetry breaking and develop through cleavage divisions to produce a prolonged multicellular colony with distinct co-existing cell types. Our findings about the autonomous and palintomic developmental program of C. perkinsii hint that such multicellular development either is much older than previously thought or evolved convergently in ichthyosporeans.

Comprehensive analysis of computational approaches in plant transcription factors binding regions discovery

Transcription factors (TFs) are regulatory proteins which bind to a specific DNA region known as the transcription factor binding regions (TFBRs) to regulate the rate of transcription process. The identification of TFBRs has been made possible by a number of experimental and computational techniques established during the past few years. The process of TFBR identification involves peak identification in the binding data, followed by the identification of motif characteristics. Using the same binding data attempts have been made to raise computational models to identify such binding regions which could save time and resources spent for binding experiments. These computational approaches depend a lot on what way they learn and how. These existing computational approaches are skewed heavily around human TFBRs discovery, while plants have drastically different genomic setup for regulation which these approaches have grossly ignored. Here, we provide a comprehensive study of the current state of the matters in plant specific TF discovery algorithms. While doing so, we encountered several software tools' issues rendering the tools not useable to researches. We fixed them and have also provided the corrected scripts for such tools. We expect this study to serve as a guide for better understanding of software tools' approaches for plant specific TFBRs discovery and the care to be taken while applying them, especially during cross-species applications. The corrected scripts of these software tools are made available at https://github.com/SCBB-LAB/Comparative-analysis-of-plant-TFBS-software.

The homeodomain regulates stable DNA binding of prostate cancer target ONECUT2

The CUT and homeodomain are ubiquitous DNA binding elements often tandemly arranged in multiple transcription factor families. However, how the CUT and homeodomain work concertedly to bind DNA remains unknown. Using ONECUT2, a driver and therapeutic target of advanced prostate cancer, we show that while the CUT initiates DNA binding, the homeodomain thermodynamically stabilizes the ONECUT2-DNA complex through allosteric modulation of CUT. We identify an arginine pair in the ONECUT family homeodomain that can adapt to DNA sequence variations. Base interactions by this ONECUT family-specific arginine pair as well as the evolutionarily conserved residues are critical for optimal DNA binding and ONECUT2 transcriptional activity in a prostate cancer model. The evolutionarily conserved base interactions additionally determine the ONECUT2-DNA binding energetics. These findings provide insights into the cooperative DNA binding by CUT-homeodomain proteins.

CONSTANS alters the circadian clock in Arabidopsis thaliana

Plants are sessile organisms that have acquired highly plastic developmental strategies to adapt to the environment. Among these processes, the floral transition is essential to ensure reproductive success and is finely regulated by several internal and external genetic networks. The photoperiodic pathway, which controls plant response to day length, is one of the most important pathways controlling flowering. In Arabidopsis photoperiodic flowering, CONSTANS (CO) is the central gene activating the expression of the florigen FLOWERING LOCUS T (FT) in the leaves at the end of a long day. The circadian clock strongly regulates CO expression. However, to date, no evidence has been reported regarding a feedback loop from the photoperiod pathway back to the circadian clock. Using transcriptional networks, we have identified relevant network motifs regulating the interplay between the circadian clock and the photoperiod pathway. Gene expression, chromatin immunoprecipitation experiments, and phenotypic analysis allowed us to elucidate the role of CO over the circadian clock. Plants with altered CO expression showed a different internal clock period, measured by daily leaf rhythmic movements. We showed that CO upregulates the expression of key genes related to the circadian clock, such as CCA1, LHY, PRR5, and GI, at the end of a long day by binding to specific sites on their promoters. Moreover, a high number of PRR5-repressed target genes are upregulated by CO, and this could explain the phase transition promoted by CO. The CO-PRR5 complex interacts with the bZIP transcription factor HY5 and helps to localize the complex in the promoters of clock genes. Taken together, our results indicate that there may be a feedback loop in which CO communicates back to the circadian clock, providing seasonal information to the circadian system.

Transcription factor binding specificities of the oomycete Phytophthora infestans reflect conserved and divergent evolutionary patterns and predict function

BACKGROUND

Identifying the DNA-binding specificities of transcription factors (TF) is central to understanding gene networks that regulate growth and development. Such knowledge is lacking in oomycetes, a microbial eukaryotic lineage within the stramenopile group. Oomycetes include many important plant and animal pathogens such as the potato and tomato blight agent Phytophthora infestans, which is a tractable model for studying life-stage differentiation within the group.

RESULTS

Mining of the P. infestans genome identified 197 genes encoding proteins belonging to 22 TF families. Their chromosomal distribution was consistent with family expansions through unequal crossing-over, which were likely ancient since each family had similar sizes in most oomycetes. Most TFs exhibited dynamic changes in RNA levels through the P. infestans life cycle. The DNA-binding preferences of 123 proteins were assayed using protein-binding oligonucleotide microarrays, which succeeded with 73 proteins from 14 families. Binding sites predicted for representatives of the families were validated by electrophoretic mobility shift or chromatin immunoprecipitation assays. Consistent with the substantial evolutionary distance of oomycetes from traditional model organisms, only a subset of the DNA-binding preferences resembled those of human or plant orthologs. Phylogenetic analyses of the TF families within P. infestans often discriminated clades with canonical and novel DNA targets. Paralogs with similar binding preferences frequently had distinct patterns of expression suggestive of functional divergence. TFs were predicted to either drive life stage-specific expression or serve as general activators based on the representation of their binding sites within total or developmentally-regulated promoters. This projection was confirmed for one TF using synthetic and mutated promoters fused to reporter genes in vivo.

CONCLUSIONS

We established a large dataset of binding specificities for P. infestans TFs, representing the first in the stramenopile group. This resource provides a basis for understanding transcriptional regulation by linking TFs with their targets, which should help delineate the molecular components of processes such as sporulation and host infection. Our work also yielded insight into TF evolution during the eukaryotic radiation, revealing both functional conservation as well as diversification across kingdoms.

Gene expression networks regulated by human personality

Genome-wide association studies of human personality have been carried out, but transcription of the whole genome has not been studied in relation to personality in humans. We collected genome-wide expression profiles of adults to characterize the regulation of expression and function in genes related to human personality. We devised an innovative multi-omic approach to network analysis to identify the key control elements and interactions in multi-modular networks. We identified sets of transcribed genes that were co-expressed in specific brain regions with genes known to be associated with personality. Then we identified the minimum networks for the co-localized genes using bioinformatic resources. Subjects were 459 adults from the Young Finns Study who completed the Temperament and Character Inventory and provided peripheral blood for genomic and transcriptomic analysis. We identified an extrinsic network of 45 regulatory genes from seed genes in brain regions involved in self-regulation of emotional reactivity to extracellular stimuli (e.g., self-regulation of anxiety) and an intrinsic network of 43 regulatory genes from seed genes in brain regions involved in self-regulation of interpretations of meaning (e.g., production of concepts and language). We discovered that interactions between the two networks were coordinated by a control hub of 3 miRNAs and 3 protein-coding genes shared by both. Interactions of the control hub with proteins and ncRNAs identified more than 100 genes that overlap directly with known personality-related genes and more than another 4000 genes that interact indirectly. We conclude that the six-gene hub is the crux of an integrative network that orchestrates information-transfer throughout a multi-modular system of over 4000 genes enriched in liquid-liquid-phase-separation (LLPS)-related RNAs, diverse transcription factors, and hominid-specific miRNAs and lncRNAs. Gene expression networks associated with human personality regulate neuronal plasticity, epigenesis, and adaptive functioning by the interactions of salience and meaning in self-awareness.

Evolution of Whirly1 in the angiosperms: sequence, splicing, and expression in a clade of early transitional mycoheterotrophic orchids

The plastid-targeted transcription factor Whirly1 (WHY1) has been implicated in chloroplast biogenesis, plastid genome stability, and fungal defense response, which together represent characteristics of interest for the study of autotrophic losses across the angiosperms. While gene loss in the plastid and nuclear genomes has been well studied in mycoheterotrophic plants, the evolution of the molecular mechanisms impacting genome stability is completely unknown. Here, we characterize the evolution of WHY1 in four early transitional mycoheterotrophic orchid species in the genus Corallorhiza by synthesizing the results of phylogenetic, transcriptomic, and comparative genomic analyses with WHY1 genomic sequences sampled from 21 orders of angiosperms. We found an increased number of non-canonical WHY1 isoforms assembled from all but the greenest Corallorhiza species, including intron retention in some isoforms. Within Corallorhiza, phylotranscriptomic analyses revealed the presence of tissue-specific differential expression of WHY1 in only the most photosynthetically capable species and a coincident increase in the number of non-canonical WHY1 isoforms assembled from fully mycoheterotrophic species. Gene- and codon-level tests of WHY1 selective regimes did not infer significant signal of either relaxed selection or episodic diversifying selection in Corallorhiza but did so for relaxed selection in the late-stage full mycoheterotrophic orchids Epipogium aphyllum and Gastrodia elata. Additionally, nucleotide substitutions that most likely impact the function of WHY1, such as nonsense mutations, were only observed in late-stage mycoheterotrophs. We propose that our findings suggest that splicing and expression changes may precede the selective shifts we inferred for late-stage mycoheterotrophic species, which therefore does not support a primary role for WHY1 in the transition to mycoheterotrophy in the Orchidaceae. Taken together, this study provides the most comprehensive view of WHY1 evolution across the angiosperms to date.

Stability and dynamics of extradenticle modulates its function

Extradenticle (EXD) is a partner protein of the HOX transcription factors and plays an important role in the development of Drosophila. It confers increased affinity and specificity of DNA-binding to the HOX proteins. However, the DNA-binding homeodomain of EXD has a significantly weaker affinity to DNA compared to the HOX homeodomains. Here, we show that a glycine residue (G290) in the middle of the EXD DNA-binding helix primarily results in this weaker binding. Glycine destabilizes helices. To probe its role in the stability and function of the protein, G290 was mutated to alanine. The intrinsic stability of the DNA-binding helix increased in the G290A mutant as observed by NMR studies and molecular dynamics (MD) simulation. Also, NMR dynamics and MD simulation show that dynamic motions present in the wild-type protein are quenched in the mutant. This in turn resulted in increased stability of the entire homeodomain (ΔΔGG→A of -2.6 kcal/mol). Increased protein stability resulted in three-fold better DNA-binding affinity of the mutant as compared to the wild-type protein. Molecular mechanics with generalized Born and surface area solvation (MMGBSA) analysis of our MD simulation on DNA-bound models of both wild-type and mutant proteins shows that the contribution to binding is enhanced for most of the interface residues in the mutant compared to the wild-type. Interestingly, the flexible N-terminal arm makes more stable contact with the DNA minor groove in the mutant. We found that the two interaction sites i.e. the DNA-binding helix and the unstructured N-terminal arm influence each other via the bound DNA. These results provide an interesting conundrum: alanine at position 290 enhances both the stability and the DNA-binding affinity of the protein, however, evolution prefers glycine at this position. We have provided several plausible explanations for this apparent conundrum. The function of the EXD as a HOX co-factor requires its ability to discriminate similar DNA sequences, which is most likely comprom.

Evolutionary divergence of subgenomes in common carp provides insights into speciation and allopolyploid success

Hybridization and polyploidization have made great contributions to speciation, heterosis, and agricultural production within plants, but there is still limited understanding and utilization in animals. Subgenome structure and expression reorganization and cooperation post hybridization and polyploidization are essential for speciation and allopolyploid success. However, the mechanisms have not yet been comprehensively assessed in animals. Here, we produced a high-fidelity reference genome sequence for common carp, a typical allotetraploid fish species cultured worldwide. This genome enabled in-depth analysis of the evolution of subgenome architecture and expression responses. Most genes were expressed with subgenome biases, with a trend of transition from the expression of subgenome A during the early stages to that of subgenome B during the late stages of embryonic development. While subgenome A evolved more rapidly, subgenome B contributed to a greater level of expression during development and under stressful conditions. Stable dominant patterns for homoeologous gene pairs both during development and under thermal stress suggest a potential fixed heterosis in the allotetraploid genome. Preferentially expressing either copy of a homoeologous gene at higher levels to confer development and response to stress indicates the dominant effect of heterosis. The plasticity of subgenomes and their shifting of dominant expression during early development, and in response to stressful conditions, provide novel insights into the molecular basis of the successful speciation, evolution, and heterosis of the allotetraploid common carp.

Species-Specific Transcription Factors Associated with Long Terminal Repeat Promoters of Endogenous Retroviruses: A Comprehensive Review

Endogenous retroviruses (ERVs) became a part of the eukaryotic genome through endogenization millions of years ago. Moreover, they have lost their innate capability of virulence or replication. Nevertheless, in eukaryotic cells, they actively engage in various activities that may be advantageous or disadvantageous to the cells. The mechanisms by which transcription is triggered and implicated in cellular processes are complex. Owing to the diversity in the expression of transcription factors (TFs) in cells and the TF-binding motifs of viruses, the comprehensibility of ERV initiation and its impact on cellular functions are unclear. Currently, several factors are known to be related to their initiation. TFs that bind to the viral long-terminal repeat (LTR) are critical initiators. This review discusses the TFs shown to actively associate with ERV stimulation across species such as humans, mice, pigs, monkeys, zebrafish, Drosophila, and yeast. A comprehensive summary of the expression of previously reported TFs may aid in identifying similarities between animal species and endogenous viruses. Moreover, an in-depth understanding of ERV expression will assist in elucidating their physiological roles in eukaryotic cell development and in clarifying their relationship with endogenous retrovirus-associated diseases.

Participation of the TBP-associated factors (TAFs) in cell differentiation

The understanding of the mechanisms that regulate gene expression to establish differentiation programs and determine cell lineages, is one of the major challenges in Developmental Biology. Besides the participation of tissue-specific transcription factors and epigenetic processes, the role of general transcription factors has been ignored. Only in recent years, there have been scarce studies that address this issue. Here, we review the studies on the biological activity of some TATA-box binding protein (TBP)-associated factors (TAFs) during the proliferation of stem/progenitor cells and their involvement in cell differentiation. Particularly, the accumulated evidence suggests that TAF4, TAF4b, TAF7L, TAF8, TAF9, and TAF10, among others, participate in nervous system development, adipogenesis, myogenesis, and epidermal differentiation; while TAF1, TAF7, TAF15 may be involved in the regulation of stem cell proliferative abilities and cell cycle progression. On the other hand, evidence suggests that TBP variants such as TBPL1 and TBPL2 might be regulating some developmental processes such as germ cell maturation and differentiation, myogenesis, or ventral specification during development. Our analysis shows that it is necessary to study in greater depth the biological function of these factors and its participation in the assembly of specific transcription complexes that contribute to the differential gene expression that gives rise to the great diversity of cell types existing in an organism. The understanding of TAFs' regulation might lead to the development of new therapies for patients which suffer from mutations, alterations, and dysregulation of these essential elements of the transcriptional machinery.

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.

Divergence time estimates for the hypoxia-inducible factor-1 alpha (HIF1α) reveal an ancient emergence of animals in low-oxygen environments

Unveiling the tempo and mode of animal evolution is necessary to understand the links between environmental changes and biological innovation. Although the earliest unambiguous metazoan fossils date to the late Ediacaran period, molecular clock estimates agree that the last common ancestor (LCA) of all extant animals emerged ~850 Ma, in the Tonian period, before the oldest evidence for widespread ocean oxygenation at ~635-560 Ma in the Ediacaran period. Metazoans are aerobic organisms, that is, they are dependent on oxygen to survive. In low-oxygen conditions, most animals have an evolutionarily conserved pathway for maintaining oxygen homeostasis that triggers physiological changes in gene expression via the hypoxia-inducible factor (HIFa). However, here we confirm the absence of the characteristic HIFa protein domain responsible for the oxygen sensing of HIFa in sponges and ctenophores, indicating the LCA of metazoans lacked the functional protein domain as well, and so could have maintained their transcription levels unaltered under the very low-oxygen concentrations of their environments. Using Bayesian relaxed molecular clock dating, we inferred that the ancestral gene lineage responsible for HIFa arose in the Mesoproterozoic Era, ~1273 Ma (Credibility Interval 957-1621 Ma), consistent with the idea that important genetic machinery associated with animals evolved much earlier than the LCA of animals. Our data suggest at least two duplication events in the evolutionary history of HIFa, which generated three vertebrate paralogs, products of the two successive whole-genome duplications that occurred in the vertebrate LCA. Overall, our results support the hypothesis of a pre-Tonian emergence of metazoans under low-oxygen conditions, and an increase in oxygen response elements during animal evolution.

Topological structures and syntenic conservation in sea anemone genomes

There is currently little information about the evolution of gene clusters, genome architectures and karyotypes in early branching animals. Slowly evolving anthozoan cnidarians can be particularly informative about the evolution of these genome features. Here we report chromosome-level genome assemblies of two related anthozoans, the sea anemones Nematostella vectensis and Scolanthus callimorphus. We find a robust set of 15 chromosomes with a clear one-to-one correspondence between the two species. Both genomes show chromosomal conservation, allowing us to reconstruct ancestral cnidarian and metazoan chromosomal blocks, consisting of at least 19 and 16 ancestral linkage groups, respectively. We show that, in contrast to Bilateria, the Hox and NK clusters of investigated cnidarians are largely disintegrated, despite the presence of staggered hox/gbx expression in Nematostella. This loss of microsynteny conservation may be facilitated by shorter distances between cis-regulatory sequences and their cognate transcriptional start sites. We find no clear evidence for topologically associated domains, suggesting fundamental differences in long-range gene regulation compared to vertebrates. These data suggest that large sets of ancestral metazoan genes have been retained in ancestral linkage groups of some extant lineages; yet, higher order gene regulation with associated 3D architecture may have evolved only after the cnidarian-bilaterian split.

Genomes of fungi and relatives reveal delayed loss of ancestral gene families and evolution of key fungal traits

Fungi are ecologically important heterotrophs that have radiated into most niches on Earth and fulfil key ecological services. Despite intense interest in their origins, major genomic trends of their evolutionary route from a unicellular opisthokont ancestor to derived multicellular fungi remain poorly known. Here we provide a highly resolved genome-wide catalogue of gene family changes across fungal evolution inferred from the genomes of 123 fungi and relatives. We show that a dominant trend in early fungal evolution has been the gradual shedding of protist genes and the punctuated emergence of innovation by two main gene duplication events. We find that the gene content of non-Dikarya fungi resembles that of unicellular opisthokonts in many respects, owing to the conservation of protist genes in their genomes. The most rapidly duplicating gene groups included extracellular proteins and transcription factors, as well as ones linked to the coordination of nutrient uptake with growth, highlighting the transition to a sessile osmotrophic feeding strategy and subsequent lifestyle evolution as important elements of early fungal history. These results suggest that the genomes of pre-fungal ancestors evolved into the typical filamentous fungal genome by a combination of gradual gene loss, turnover and several large duplication events rather than by abrupt changes. Consequently, the taxonomically defined Fungi represents a genomically non-uniform assemblage of species.

The compact genome of the sponge Oopsacas minuta (Hexactinellida) is lacking key metazoan core genes

BACKGROUND

Explaining the emergence of the hallmarks of bilaterians is a central focus of evolutionary developmental biology-evodevo-and evolutionary genomics. For this purpose, we must both expand and also refine our knowledge of non-bilaterian genomes, especially by studying early branching animals, in particular those in the metazoan phylum Porifera.

RESULTS

We present a comprehensive analysis of the first whole genome of a glass sponge, Oopsacas minuta, a member of the Hexactinellida. Studying this class of sponge is evolutionary relevant because it differs from the three other Porifera classes in terms of development, tissue organization, ecology, and physiology. Although O. minuta does not exhibit drastic body simplifications, its genome is among the smallest of animal genomes sequenced so far, and surprisingly lacks several metazoan core genes (including Wnt and several key transcription factors). Our study also provides the complete genome of a symbiotic Archaea dominating the associated microbial community: a new Thaumarchaeota species.

CONCLUSIONS

The genome of the glass sponge O. minuta differs from all other available sponge genomes by its compactness and smaller number of encoded proteins. The unexpected loss of numerous genes previously considered ancestral and pivotal for metazoan morphogenetic processes most likely reflects the peculiar syncytial tissue organization in this group. Our work further documents the importance of convergence during animal evolution, with multiple convergent evolution of septate-like junctions, electrical-signaling and multiciliated cells in metazoans.

The homeodomain drives favorable DNA binding energetics of prostate cancer target ONECUT2

The ONECUT transcription factors feature a CUT and a homeodomain, evolutionarily conserved elements that bind DNA cooperatively, but the process remains mechanistically enigmatic. Using an integrative DNA binding analysis of ONECUT2, a driver of aggressive prostate cancer, we show that the homeodomain energetically stabilizes the ONECUT2-DNA complex through allosteric modulation of CUT. Further, evolutionarily conserved base-interactions in both the CUT and homeodomain are necessary for the favorable thermodynamics. We have identified a novel arginine pair unique to the ONECUT family homeodomain that can adapt to DNA sequence variations. Base interactions in general, including by this arginine pair, are critical for optimal DNA binding and transcription in a prostate cancer model. These findings provide fundamental insights into DNA binding by CUT-homeodomain proteins with potential therapeutic implications.

One-Sentence Summary

Base-specific interactions regulate homeodomain-mediated stabilization of DNA binding by the ONECUT2 transcription factor.

The START domain potentiates HD-ZIPIII transcriptional activity

The CLASS III HOMEODOMAIN-LEUCINE ZIPPER (HD-ZIPIII) transcription factors (TFs) were repeatedly deployed over 725 million years of evolution to regulate central developmental innovations. The START domain of this pivotal class of developmental regulators was recognized over 20 years ago, but its putative ligands and functional contributions remain unknown. Here, we demonstrate that the START domain promotes HD-ZIPIII TF homodimerization and increases transcriptional potency. Effects on transcriptional output can be ported onto heterologous TFs, consistent with principles of evolution via domain capture. We also show the START domain binds several species of phospholipids, and that mutations in conserved residues perturbing ligand binding and/or its downstream conformational readout abolish HD-ZIPIII DNA-binding competence. Our data present a model in which the START domain potentiates transcriptional activity and uses ligand-induced conformational change to render HD-ZIPIII dimers competent to bind DNA. These findings resolve a long-standing mystery in plant development and highlight the flexible and diverse regulatory potential coded within this widely distributed evolutionary module.

Uncovering gene-family founder events during major evolutionary transitions in animals, plants and fungi using GenEra

We present GenEra ( https://github.com/josuebarrera/GenEra ), a DIAMOND-fueled gene-family founder inference framework that addresses previously raised limitations and biases in genomic phylostratigraphy, such as homology detection failure. GenEra also reduces computational time from several months to a few days for any genome of interest. We analyze the emergence of taxonomically restricted gene families during major evolutionary transitions in plants, animals, and fungi. Our results indicate that the impact of homology detection failure on inferred patterns of gene emergence is lineage-dependent, suggesting that plants are more prone to evolve novelty through the emergence of new genes compared to animals and fungi.

Integrating transcription and splicing into cell fate: Transcription factors on the block

Transcription factors (TFs) are present in all life forms and conserved across great evolutionary distances in eukaryotes. From yeast to complex multicellular organisms, they are pivotal players of cell fate decision by orchestrating gene expression at diverse molecular layers. Notably, TFs fine-tune gene expression by coordinating RNA fate at both the expression and splicing levels. They regulate alternative splicing, an essential mechanism for cell plasticity, allowing the production of many mRNA and protein isoforms in precise cell and tissue contexts. Despite this apparent role in splicing, how TFs integrate transcription and splicing to ultimately orchestrate diverse cell functions and cell fate decisions remains puzzling. We depict substantial studies in various model organisms underlining the key role of TFs in alternative splicing for promoting tissue-specific functions and cell fate. Furthermore, we emphasize recent advances describing the molecular link between the transcriptional and splicing activities of TFs. As TFs can bind both DNA and/or RNA to regulate transcription and splicing, we further discuss their flexibility and compatibility for DNA and RNA substrates. Finally, we propose several models integrating transcription and splicing activities of TFs in the coordination and diversification of cell and tissue identities. This article is categorized under: RNA Processing > Splicing Regulation/Alternative Splicing RNA Interactions with Proteins and Other Molecules > Protein-RNA Interactions: Functional Implications RNA Processing > Splicing Mechanisms.

The complex network of transcription factors, immune checkpoint inhibitors and stemness features in colorectal cancer: A recent update

Cancer immunity is regulated by several mechanisms that include co-stimulatory and/or co-inhibitory molecules known as immune checkpoints expressed by the immune cells. In colorectal cancer (CRC), CTLA-4, LAG3, TIM-3 and PD-1 are the major co-inhibitory checkpoints involved in tumor development and progression. On the other hand, the deregulation of transcription factors and cancer stem cells activity plays a major role in the development of drug resistance and in the spread of metastatic disease in CRC. In this review, we describe how the modulation of such transcription factors affects the response of CRC to therapies. We also focus on the role of cancer stem cells in tumor metastasis and chemoresistance and discuss both preclinical and clinical approaches for targeting stem cells to prevent their tumorigenic effect. Finally, we provide an update on the clinical applications of immune checkpoint inhibitors in CRC and discuss the regulatory effects of transcription factors on the expression of the immune inhibitory checkpoints with specific focus on the PD-1 and PD-L1 molecules.

The premetazoan ancestry of the synaptic toolkit and appearance of first neurons

Neurons, especially when coupled with muscles, allow animals to interact with and navigate through their environment in ways unique to life on earth. Found in all major animal lineages except sponges and placozoans, nervous systems range widely in organization and complexity, with neurons possibly representing the most diverse cell-type. This diversity has led to much debate over the evolutionary origin of neurons as well as synapses, which allow for the directed transmission of information. The broad phylogenetic distribution of neurons and presence of many of the defining components outside of animals suggests an early origin of this cell type, potentially in the time between the first animal and the last common ancestor of extant animals. Here, we highlight the occurrence and function of key aspects of neurons outside of animals as well as recent findings from non-bilaterian animals in order to make predictions about when and how the first neuron(s) arose during animal evolution and their relationship to those found in extant lineages. With advancing technologies in single cell transcriptomics and proteomics as well as expanding functional techniques in non-bilaterian animals and the close relatives of animals, it is an exciting time to begin unraveling the complex evolutionary history of this fascinating animal cell type.

Predicted landscape of RETINOBLASTOMA-RELATED LxCxE-mediated interactions across the Chloroplastida

The colonization of land by a single streptophyte algae lineage some 450 million years ago has been linked to multiple key innovations such as three-dimensional growth, alternation of generations, the presence of stomata, as well as innovations inherent to the birth of major plant lineages, such as the origins of vascular tissues, roots, seeds and flowers. Multicellularity, which evolved multiple times in the Chloroplastida coupled with precise spatiotemporal control of proliferation and differentiation were instrumental for the evolution of these traits. RETINOBLASTOMA-RELATED (RBR), the plant homolog of the metazoan Retinoblastoma protein (pRB), is a highly conserved and multifunctional core cell cycle regulator that has been implicated in the evolution of multicellularity in the green lineage as well as in plant multicellularity-related processes such as proliferation, differentiation, stem cell regulation and asymmetric cell division. RBR fulfills these roles through context-specific protein-protein interactions with proteins containing the Leu-x-Cys-x-Glu (LxCxE) short-linear motif (SLiM); however, how RBR-LxCxE interactions have changed throughout major innovations in the Viridiplantae kingdom is a question that remains unexplored. Here, we employ an in silico evo-devo approach to predict and analyze potential RBR-LxCxE interactions in different representative species of key Chloroplastida lineages, providing a valuable resource for deciphering RBR-LxCxE multiple functions. Furthermore, our analyses suggest that RBR-LxCxE interactions are an important component of RBR functions and that interactions with chromatin modifiers/remodelers, DNA replication and repair machinery are highly conserved throughout the Viridiplantae, while LxCxE interactions with transcriptional regulators likely diversified throughout the water-to-land transition.

S1PR1 induces metabolic reprogramming of ceramide in vascular endothelial cells, affecting hepatocellular carcinoma angiogenesis and progression

Angiogenesis is a fundamental process underlying the occurrence, growth and metastasis of hepatocellular carcinoma (HCC), a prevalent tumour type with an extremely poor prognosis due to abundant vasculature. However, the underlying mechanism of angiogenesis in HCC remains largely unknown. Herein, we found that sphingosine-1-phosphate receptor 1 (S1PR1) plays an important role in HCC angiogenesis. S1PR1 was found to be selectively and highly expressed in the blood vessels of HCC tissues compared with those of paratumour tissues. Functionally, high expression of S1PR1 in endothelial cells (ECs) promoted angiogenesis and progression of HCC in vitro and in vivo. Mechanistically, proangiogenic factors (S1P, IL-6, VEGFA) in conditioned medium from HCC cells induced the upregulation of S1PR1 in ECs via the phosphorylation of STAT3 at Y705. Further study also revealed that S1PR1 promotes angiogenesis by decreasing ceramide levels via CerS3 downregulation. Interestingly, we demonstrated that S1PR1 downregulates CerS3 by inducing CerS6 translocation into the nucleus to inhibit CerS3 at the transcriptional level in ECs. In addition, we found that a high concentration of Lenvatinib significantly downregulated the expression of S1PR1 and obviously enhanced S1PR1 knockdown-mediated angiogenesis inhibition, indicating that S1PR1 may be a target by which Lenvatinib combats angiogenesis in HCC. Thus, S1PR1 may be an important target for suppressing angiogenesis in HCC, and inhibiting S1PR1 is a promising approach to antitumor therapy in HCC.

Zinc Finger-Homeodomain and Mini Zinc Finger proteins are key players in plant growth and responses to environmental stresses

The ZINC FINGER-HOMEODOMAIN (ZHD) protein family is a plant-specific family of transcription factors containing two conserved motifs: a non-canonical C5H3 zinc finger domain (ZF) and a DNA-binding homeodomain (HD). The MINI ZINC FINGER (MIF) proteins belong to this family, but were possibly derived from the ZHDs by losing the HD. Information regarding the function of ZHD and MIF proteins is scarce. However, different studies have shown that ZHD/MIF proteins play important roles not only in plant growth and development, but also in response to environmental stresses, including drought and pathogen attack. Here we review recent advances relative to ZHD/MIF functions in multiple species, to provide new insights into the diverse roles of these transcription factors in plants. Their mechanism of action in relation to their ability to interact with other proteins and DNA is also discussed. We then propose directions for future studies to understand better their important roles and pinpoint strategies for potential applications in crop improvement.

Structural and functional analysis of the human cone‐rod homeobox transcription factor

The cone‐rod homeobox (CRX) protein is a critical K50 homeodomain transcription factor responsible for the differentiation and maintenance of photoreceptor neurons in the vertebrate retina. Mutant alleles in the human gene encoding CRX result in a variety of distinct blinding retinopathies, including retinitis pigmentosa, cone‐rod dystrophy, and Leber congenital amaurosis. Despite the success of using in vitro biochemistry, animal models, and genomics approaches to study this clinically relevant transcription factor over the past 25 years since its initial characterization, there are no high‐resolution structures in the published literature for the CRX protein. In this study, we use bioinformatic approaches and small‐angle X‐ray scattering (SAXS) structural analysis to further understand the biochemical complexity of the human CRX homeodomain (CRX‐HD). We find that the CRX‐HD is a compact, globular monomer in solution that can specifically bind functional cis‐regulatory elements encoded upstream of retina‐specific genes. This study presents the first structural analysis of CRX, paving the way for a new approach to studying the biochemistry of this protein and its disease‐causing mutant protein variants.

Moss transcription factors regulating development and defense responses to stress

Transcription factors control gene expression, leading to regulation of biological processes that determine plant development and adaptation to the environment. Land colonization by plants occurred 450-470 million years ago and was accompanied by an increase in the complexity of transcriptional regulation associated to transcription factor gene expansions. AP2/ERF, bHLH, MYB, NAC, GRAS, and WRKY transcription factor families increased in land plants compared with algae. In angiosperms, they play crucial roles in regulating plant growth and responses to environmental stressors. However, less information is available in bryophytes and only in a few cases is the functional role of moss transcription factors in stress mechanisms known. In this review, we discuss current knowledge of the transcription factor families involved in development and defense responses to stress in mosses and other bryophytes. By exploring and analysing the Physcomitrium patens public database and published transcriptional profiles, we show that a high number of AP2/ERF, bHLH, MYB, NAC, GRAS, and WRKY genes are differentially expressed in response to abiotic stresses and during biotic interactions. Expression profiles together with a comprehensive analysis provide insights into relevant transcription factors involved in moss defenses, and hint at distinct and conserved biological roles between bryophytes and angiosperms.

Origin, Expansion, and Divergence of ETHYLENE-INSENSITIVE 3 (EIN3)/EIN3-LIKE Transcription Factors During Streptophytes Evolution

The transition of plants to land required several regulatory adaptive mechanisms. Little is known about these mechanisms, but they no doubt involved the evolution of transcription factor (TF) families. ETHYLENE-INSENSITIVE 3 (EIN3)/EIN3-LIKE (EIL) transcription factors (TFs) are core components of the ethylene signaling pathway that play important roles in almost every aspect of plant development and environmental responses by regulating the transcription of numerous genes. However, the evolutionary history of EIN3/EIL TFs, which are present in a wide range of streptophytes, is still not clear. Here, to explore the evolution and functions of EIN3/EIL TFs, we performed phylogenetic analysis of these TFs and investigated their gene and protein structures as well as sequence features. Our results suggest that the EIN3/EIL TF family was already was already present in the ancestor of streptophyte algae. Phylogenetic analysis divided the EIN3/EIL TFs into three groups (Group A-C). Analysis of gene and protein structure revealed that most of the structural features of these TFs had already formed in ancient lineages. Further investigation suggested that all groups have undergone several duplication events related to whole-genome duplications in plants, generating multiple, functionally diverse gene copies. Therefore, as plants colonized terrestrial habitats and evolved key traits, the EIN3/EIL TF family expanded broadly via multiple duplication events, which could have promoted their fundamental neo- and sub-functionalization to help plants adapt to terrestrial life. Our findings shed light on the functional evolution of the EIN3/EIL TF family in the streptophytes.

Rel/NF-κB Transcription Factors Emerged at the Onset of Opisthokonts

The Rel/NF-κB transcription factor family has myriad roles in immunity, development, and differentiation in animals, and was considered a key innovation for animal multicellularity. Rel homology domain-containing proteins were previously hypothesized to have originated in a last common ancestor of animals and some of their closest unicellular relatives. However, key taxa were missing from previous analyses, necessitating a systematic investigation into the distribution and evolution of these proteins. Here, we address this knowledge gap by surveying taxonomically broad data from eukaryotes, with a special emphasis on lineages closely related to animals. We report an earlier origin for Rel/NF-κB proteins than previously described, in the last common ancestor of animals and fungi, and show that even in the sister group to fungi, these proteins contain elements that in animals are necessary for the subcellular regulation of Rel/NF-κB.

MYB transcription factor family in sweet cherry (Prunus avium L.): genome-wide investigation, evolution, structure, characterization and expression patterns

BACK GROUND

MYB Transcription factors (TFs) are most imperative and largest gene family in plants, which participate in development, metabolism, defense, differentiation and stress response. The MYB TFs has been studied in various plant species. However, comprehensive studies of MYB gene family in the sweet cherry (Prunus avium L.) are still unknown.

RESULTS

In the current study, a total of 69 MYB genes were investigated from sweet cherry genome and classified into 28 subfamilies (C1-C28 based on phylogenetic and structural analysis). Microcollinearity analysis revealed that dispersed duplication (DSD) events might play an important role in the MYB genes family expansion. Chromosomal localization, the synonymous (Ks) and nonsynonymous (Ka) analysis, molecular characteristics (pI, weight and length of amino acids) and subcellular localization were accomplished using several bioinformatics tools. Furthermore, the members of distinct subfamilies have diverse cis-acting regions, conserved motifs, and intron-exon architectures, indicating functional heterogeneity in the MYB family. Moreover, the transcriptomic data exposed that MYB genes might play vital role in bud dormancy. The quantitative real-time qRT-PCR was carried out and the expression pattern indicated that MYB genes significantly expressed in floral bud as compared to flower and fruit.

CONCLUSION

Our comprehensive findings provide supportive insights into the evolutions, expansion complexity and functionality of PavMYB genes. These PavMYB genes should be further investigated as they seem to be brilliant candidates for dormancy manipulation in sweet cherry.

The Non-continuum Nature of Eukaryotic Transcriptional Regulation

Eukaryotic transcription factors are versatile mediators of specificity in gene regulation. This versatility is achieved through mutual specification by context-specific DNA binding on the one hand, and identity-specific protein-protein partnerships on the other. This interactivity, known as combinatorial control, enables a repertoire of complex transcriptional outputs that are qualitatively disjoint, or non-continuum, with respect to binding affinity. This feature contrasts starkly with prokaryotic gene regulators, whose activities in general vary quantitatively in step with binding affinity. Biophysical studies on prokaryotic model systems and more recent investigations on transcription factors highlight an important role for folded state dynamics and molecular hydration in protein/DNA recognition. Analysis of molecular models of combinatorial control and recent literature in low-affinity gene regulation suggest that transcription factors harbor unique conformational dynamics that are inaccessible or unused by prokaryotic DNA-binding proteins. Thus, understanding the intrinsic dynamics involved in DNA binding and co-regulator recruitment appears to be a key to understanding how transcription factors mediate non-continuum outcomes in eukaryotic gene expression, and how such capability might have evolved from ancient, structurally conserved counterparts.

Transcription Factors in the Fungus Aspergillus nidulans: Markers of Genetic Innovation, Network Rewiring and Conflict between Genomics and Transcriptomics

Gene regulatory networks (GRNs) are shaped by the democratic/hierarchical relationships among transcription factors (TFs) and associated proteins, together with the cis-regulatory sequences (CRSs) bound by these TFs at target promoters. GRNs control all cellular processes, including metabolism, stress response, growth and development. Due to the ability to modify morphogenetic and developmental patterns, there is the consensus view that the reorganization of GRNs is a driving force of species evolution and differentiation. GRNs are rewired through events including the duplication of TF-coding genes, their divergent sequence evolution and the gain/loss/modification of CRSs. Fungi (mainly Saccharomycotina) have served as a reference kingdom for the study of GRN evolution. Here, I studied the genes predicted to encode TFs in the fungus Aspergillus nidulans (Pezizomycotina). The analysis of the expansion of different families of TFs suggests that the duplication of TFs impacts the species level, and that the expansion in Zn2Cys6 TFs is mainly due to dispersed duplication events. Comparison of genomic annotation and transcriptomic data suggest that a significant percentage of genes should be re-annotated, while many others remain silent. Finally, a new regulator of growth and development is identified and characterized. Overall, this study establishes a novel theoretical framework in synthetic biology, as the overexpression of silent TF forms would provide additional tools to assess how GRNs are rewired.

Comparative Transcriptomics Reveals Distinct Patterns of Gene Expression Conservation through Vertebrate Embryogenesis

Despite life's diversity, studies of variation often remind us of our shared evolutionary past. Abundant genome sequencing and analyses of gene regulatory networks illustrate that genes and entire pathways are conserved, reused, and elaborated in the evolution of diversity. Predating these discoveries, 19th-century embryologists observed that though morphology at birth varies tremendously, certain stages of vertebrate embryogenesis appear remarkably similar across vertebrates. In the mid to late 20th century, anatomical variability of early and late-stage embryos and conservation of mid-stages embryos (the "phylotypic" stage) was named the hourglass model of diversification. This model has found mixed support in recent analyses comparing gene expression across species possibly owing to differences in species, embryonic stages, and gene sets compared. We compare 186 microarray and RNA-seq data sets covering embryogenesis in six vertebrate species. We use an unbiased clustering approach to group stages of embryogenesis by transcriptomic similarity and ask whether gene expression similarity of clustered embryonic stages deviates from a null expectation. We characterize expression conservation patterns of each gene at each evolutionary node after correcting for phylogenetic nonindependence. We find significant enrichment of genes exhibiting early conservation, hourglass, late conservation patterns in both microarray and RNA-seq data sets. Enrichment of genes showing patterned conservation through embryogenesis indicates diversification of embryogenesis may be temporally constrained. However, the circumstances under which each pattern emerges remain unknown and require both broad evolutionary sampling and systematic examination of embryogenesis across species.

Structural variation affecting DNA backbone interactions underlies adaptation of B3 DNA binding domains to constraints imposed by protein architecture

Functional and architectural diversification of transcription factor families has played a central role in the independent evolution of complex development in plants and animals. Here, we investigate the role of architectural constraints on evolution of B3 DNA binding domains that regulate plant embryogenesis. B3 domains of ABI3, FUS3, LEC2 and VAL1 proteins recognize the same cis-element. Complex architectures of ABI3 and VAL1 integrate cis-element recognition with other signals, whereas LEC2 and FUS3 have reduced architectures conducive to roles as pioneer activators. In yeast and plant in vivo assays, B3 domain functions correlate with architectural complexity of the parent transcription factor rather than phylogenetic relatedness. In a complex architecture, attenuated ABI3-B3 and VAL1-B3 activities enable integration of cis-element recognition with hormone signaling, whereas hyper-active LEC2-B3 and FUS3-B3 over-ride hormonal control. Three clade-specific amino acid substitutions (β4-triad) implicated in interactions with the DNA backbone account for divergence of LEC2-B3 and ABI3-B3. We find a striking correlation between differences in in vitro DNA binding affinity and in vivo activities of B3 domains in plants and yeast. Our results highlight the role of DNA backbone interactions that preserve DNA sequence specificity in adaptation of B3 domains to functional constraints associated with domain architecture.

Direct molecular evidence for an ancient, conserved developmental toolkit controlling post-transcriptional gene regulation in land plants

In plants, miRNA production is orchestrated by a suite of proteins that control transcription of the pri-miRNA gene, post-transcriptional processing and nuclear export of the mature miRNA. Post-transcriptional processing of miRNAs is controlled by a pair of physically interacting proteins, hyponastic leaves 1 (HYL1) and Dicer-like 1 (DCL1). However, the evolutionary history and structural basis of the HYL1–DCL1 interaction is unknown. Here we use ancestral sequence reconstruction and functional characterization of ancestral HYL1 in vitro and in Arabidopsis thaliana to better understand the origin and evolution of the HYL1–DCL1 interaction and its impact on miRNA production and plant development. We found the ancestral plant HYL1 evolved high affinity for both double-stranded RNA (dsRNA) and its DCL1 partner before the divergence of mosses from seed plants (∼500 Ma), and these high-affinity interactions remained largely conserved throughout plant evolutionary history. Structural modeling and molecular binding experiments suggest that the second of two dsRNA-binding motifs (DSRMs) in HYL1 may interact tightly with the first of two C-terminal DCL1 DSRMs to mediate the HYL1–DCL1 physical interaction necessary for efficient miRNA production. Transgenic expression of the nearly 200 Ma-old ancestral flowering-plant HYL1 in A. thaliana was sufficient to rescue many key aspects of plant development disrupted by HYL1⁻ knockout and restored near-native miRNA production, suggesting that the functional partnership of HYL1–DCL1 originated very early in and was strongly conserved throughout the evolutionary history of terrestrial plants. Overall, our results are consistent with a model in which miRNA-based gene regulation evolved as part of a conserved plant “developmental toolkit.”

Cnidarian-bilaterian comparison reveals the ancestral regulatory logic of the β-catenin dependent axial patterning

In animals, body axis patterning is based on the concentration-dependent interpretation of graded morphogen signals, which enables correct positioning of the anatomical structures. The most ancient axis patterning system acting across animal phyla relies on β-catenin signaling, which directs gastrulation, and patterns the main body axis. However, within Bilateria, the patterning logic varies significantly between protostomes and deuterostomes. To deduce the ancestral principles of β-catenin-dependent axial patterning, we investigate the oral–aboral axis patterning in the sea anemone Nematostella—a member of the bilaterian sister group Cnidaria. Here we elucidate the regulatory logic by which more orally expressed β-catenin targets repress more aborally expressed β-catenin targets, and progressively restrict the initially global, maternally provided aboral identity. Similar regulatory logic of β-catenin-dependent patterning in Nematostella and deuterostomes suggests a common evolutionary origin of these processes and the equivalence of the cnidarian oral–aboral and the bilaterian posterior–anterior body axes.

The authors show in Nematostella that the more orally expressed β-catenin targets repress the more aborally expressed β-catenin targets, thus patterning the oral-aboral axis. This likely represents the common mechanism of β-catenin-dependent axial patterning shared by Cnidaria and Bilateria.

Mechanisms of enhancer action: the known and the unknown

Differential gene expression mechanisms ensure cellular differentiation and plasticity to shape ontogenetic and phylogenetic diversity of cell types. A key regulator of differential gene expression programs are the enhancers, the gene-distal cis-regulatory sequences that govern spatiotemporal and quantitative expression dynamics of target genes. Enhancers are widely believed to physically contact the target promoters to effect transcriptional activation. However, our understanding of the full complement of regulatory proteins and the definitive mechanics of enhancer action is incomplete. Here, we review recent findings to present some emerging concepts on enhancer action and also outline a set of outstanding questions.

Liver regeneration observed across the different classes of vertebrates from an evolutionary perspective

The liver is a key organ that performs diverse functions such as metabolic processing of nutrients or disposal of dangerous substances (xenobiotics). Accordingly, it seems to be protected by several mechanisms throughout the life of organisms, one of which is compensatory hyperplasia, also known as liver regeneration. This review is a recapitulation of the scientific reports describing the different ways in which the various classes of vertebrates deal with liver injuries, where since mammals have an improved molecular toolkit, exhibit optimized regeneration of the liver compared to lower vertebrates. The main molecules involved in the compensatory process, such as proinflammatory and inhibitory cytokines, are analyzed across vertebrates with an evolutionary perspective. In addition, the possible significance of this mechanism is discussed in the context of the long life span of vertebrates, especially in the case of mammals.