Apprehending multicellularity: regulatory networks, genomics, and evolution

Exploring the expressiveness of abstract metabolic networks

Metabolism is characterised by chemical reactions linked to each other, creating a complex network structure. The whole metabolic network is divided into pathways of chemical reactions, such that every pathway is a metabolic function. A simplified representation of metabolism, which we call an abstract metabolic network, is a graph in which metabolic pathways are nodes and there is an edge between two nodes if their corresponding pathways share one or more compounds. The abstract metabolic network of a given organism results in a small network that requires low computational power to be analysed and makes it a suitable model to perform a large-scale comparison of organisms' metabolism. To explore the potentials and limits of such a basic representation, we considered a comprehensive set of KEGG organisms, represented through their abstract metabolic network. We performed pairwise comparisons using graph kernel methods and analyse the results through exploratory data analysis and machine learning techniques. The results show that abstract metabolic networks discriminate macro evolutionary events, indicating that they are expressive enough to capture key steps in metabolism evolution.

Genomic Variation, Evolvability, and the Paradox of Mental Illness

Twentieth-century genetics was hard put to explain the irregular behavior of neuropsychiatric disorders. Autism and schizophrenia defy a principle of natural selection; they are highly heritable but associated with low reproductive success. Nevertheless, they persist. The genetic origins of such conditions are confounded by the problem of variable expression, that is, when a given genetic aberration can lead to any one of several distinct disorders. Also, autism and schizophrenia occur on a spectrum of severity, from mild and subclinical cases to the overt and disabling. Such irregularities reflect the problem of missing heritability; although hundreds of genes may be associated with autism or schizophrenia, together they account for only a small proportion of cases. Techniques for higher resolution, genomewide analysis have begun to illuminate the irregular and unpredictable behavior of the human genome. Thus, the origins of neuropsychiatric disorders in particular and complex disease in general have been illuminated. The human genome is characterized by a high degree of structural and behavioral variability: DNA content variation, epistasis, stochasticity in gene expression, and epigenetic changes. These elements have grown more complex as evolution scaled the phylogenetic tree. They are especially pertinent to brain development and function. Genomic variability is a window on the origins of complex disease, neuropsychiatric disorders, and neurodevelopmental disorders in particular. Genomic variability, as it happens, is also the fuel of evolvability. The genomic events that presided over the evolution of the primate and hominid lineages are over-represented in patients with autism and schizophrenia, as well as intellectual disability and epilepsy. That the special qualities of the human genome that drove evolution might, in some way, contribute to neuropsychiatric disorders is a matter of no little interest.

Multifunctional Involvement of a C2H2 Zinc Finger Protein (PbZfp) in Malaria Transmission, Histone Modification, and Susceptibility to DNA Damage Response

In sexually reproducing organisms, meiosis is an essential step responsible for generation of haploid gametes from diploid somatic cells. The quest for understanding regulatory mechanisms of meiotic recombination in Plasmodium led to identification of a gene encoding a protein that contains 11 copies of C₂H₂ zinc fingers (ZnF). Reverse genetic approaches were used to create Plasmodium berghei parasites either lacking expression of full-length Plasmodium berghei zinc finger protein (PbZfp) (knockout [KO]) or expressing PbZfp lacking C-terminal zinc finger region (truncated [Trunc]). Mice infected with KO parasites survived two times longer (P < 0.0001) than mice infected with wild-type (WT) parasites. In mosquito transmission experiments, the infectivity of KO and Trunc parasites was severely compromised (>95% oocyst reduction). KO parasites revealed a total lack of trimethylation of histone 3 at several lysine residues (K4, K27, and K36) without any effect on acetylation patterns (H3K9, H3K14, and H4K16). Reduced DNA damage and reduced expression of topoisomerase-like Spo11 in the KO parasites with normal Rad51 expression further suggest a functional role for PbZfp during genetic recombination that involves DNA double-strand break (DSB) formation followed by DNA repair. These finding raise the possibility of some convergent similarities of PbZfp functions to functions of mammalian PRDM9, also a C₂H₂ ZnF protein with histone 3 lysine 4 (H3K4) methyltransferase activity. These functions include the major role played by the latter in binding recombination hotspots in the genome during meiosis and trimethylation of the associated histones and subsequent chromatin recruitment of topoisomerase-like Spo11 to catalyze DNA DSB formation and DMC1/Rad51-mediated DNA repair and homologous recombination.

Malaria parasites are haploid throughout their life cycle except for a brief time period when zygotes are produced as a result of fertilization between male and female gametes during transmission through the mosquito vector. The reciprocal recombination events that follow zygote formation ensure orderly segregation of homologous chromosomes during meiosis, creating genetic diversity among offspring. Studies presented in the current manuscript identify a novel C2H2 ZnF-containing protein exhibiting multifunctional roles in parasite virulence, mosquito transmission, and homologous recombination during meiosis. Understanding the transmission biology of malaria will result in the identification of novel targets for transmission-blocking intervention approaches.

Evolution and multiple roles of the Pancrustacea specific transcription factor zelda in insects

Gene regulatory networks (GRNs) evolve as a result of the coevolutionary processes acting on transcription factors (TFs) and the cis-regulatory modules they bind. The zinc-finger TF zelda (zld) is essential for the maternal-to-zygotic transition (MZT) in Drosophila melanogaster, where it directly binds over thousand cis-regulatory modules to regulate chromatin accessibility. D. melanogaster displays a long germ type of embryonic development, where all segments are simultaneously generated along the whole egg. However, it remains unclear if zld is also involved in the MZT of short-germ insects (including those from basal lineages) or in other biological processes. Here we show that zld is an innovation of the Pancrustacea lineage, being absent in more distant arthropods (e.g. chelicerates) and other organisms. To better understand zld´s ancestral function, we thoroughly investigated its roles in a short-germ beetle, Tribolium castaneum, using molecular biology and computational approaches. Our results demonstrate roles for zld not only during the MZT, but also in posterior segmentation and patterning of imaginal disc derived structures. Further, we also demonstrate that zld is critical for posterior segmentation in the hemipteran Rhodnius prolixus, indicating this function predates the origin of holometabolous insects and was subsequently lost in long-germ insects. Our results unveil new roles of zld in different biological contexts and suggest that changes in expression of zld (and probably other major TFs) are critical in the evolution of insect GRNs.

Pioneer transcription factors (TFs) are considered the first regulators of chromatin accessibility in fruit flies and vertebrates, modulating the expression of a large number of target genes. In fruit flies, zelda resembles a pioneer TF, being essential during early embryogenesis. However, the evolutionary origins and ancestral functions of zelda remain largely unknown. Through a number of gene silencing, microscopy and evolutionary analysis, the present work shows that zelda is an innovation of the Pancrustacea lineage, governing not only the MZT in the short-germ insect Tribolium castaneum, but also posterior segmentation and post-embryonic patterning of imaginal disc derived structures such as wings, legs and antennae. Further, zelda regulation of posterior segmentation predates the origin of insects with complete metamorphosis (holometabolous), as supported by gene silencing experiments in the kissing bug Rhodnius prolixus. We hypothesize that the emergence of zelda contributed to the evolution of gene regulatory networks and new morphological structures of insects.

Myxobacteria: Moving, Killing, Feeding, and Surviving Together

Myxococcus xanthus, like other myxobacteria, is a social bacterium that moves and feeds cooperatively in predatory groups. On surfaces, rod-shaped vegetative cells move in search of the prey in a coordinated manner, forming dynamic multicellular groups referred to as swarms. Within the swarms, cells interact with one another and use two separate locomotion systems. Adventurous motility, which drives the movement of individual cells, is associated with the secretion of slime that forms trails at the leading edge of the swarms. It has been proposed that cellular traffic along these trails contributes to M. xanthus social behavior via stigmergic regulation. However, most of the cells travel in groups by using social motility, which is cell contact-dependent and requires a large number of individuals. Exopolysaccharides and the retraction of type IV pili at alternate poles of the cells are the engines associated with social motility. When the swarms encounter prey, the population of M. xanthus lyses and takes up nutrients from nearby cells. This cooperative and highly density-dependent feeding behavior has the advantage that the pool of hydrolytic enzymes and other secondary metabolites secreted by the entire group is shared by the community to optimize the use of the degradation products. This multicellular behavior is especially observed in the absence of nutrients. In this condition, M. xanthus swarms have the ability to organize the gliding movements of 1000s of rods, synchronizing rippling waves of oscillating cells, to form macroscopic fruiting bodies, with three subpopulations of cells showing division of labor. A small fraction of cells either develop into resistant myxospores or remain as peripheral rods, while the majority of cells die, probably to provide nutrients to allow aggregation and spore differentiation. Sporulation within multicellular fruiting bodies has the benefit of enabling survival in hostile environments, and increases germination and growth rates when cells encounter favorable conditions. Herein, we review how these social bacteria cooperate and review the main cell–cell signaling systems used for communication to maintain multicellularity.

Transcription factors, chromatin proteins and the diversification of Hemiptera

Availability of complete genomes provides a means to explore the evolution of enormous developmental, morphological, and behavioral diversity among insects. Hemipterans in particular show great diversity of both morphology and life history within a single order. To better understand the role of transcription regulators in the diversification of hemipterans, using sequence profile searches and hidden Markov models we computationally analyzed transcription factors (TFs) and chromatin proteins (CPs) in the recently available Rhodnius prolixus genome along with 13 other insect and 4 non-insect arthropod genomes. We generated a comprehensive collection of TFs and CPs across arthropods including 303 distinct types of domains in TFs and 139 in CPs. This, along with the availability of two hemipteran genomes, R. prolixus and Acyrthosiphon pisum, helped us identify possible determinants for their dramatic morphological and behavioral divergence. We identified five domain families (i.e. Pipsqueak, SAZ/MADF, THAP, FLYWCH and BED finger) as having undergone differential patterns of lineage-specific expansion in hemipterans or within hemipterans relative to other insects. These expansions appear to be at least in part driven by transposons, with the DNA-binding domains of transposases having provided the raw material for emergence of new TFs. Our analysis suggests that while R. prolixus probably retains a state closer to the ancestral hemipteran, A. pisum represents a highly derived state, with the emergence of asexual reproduction potentially favoring genome duplication and transposon expansion. Both hemipterans are predicted to possess active DNA methylation systems. However, in the course of their divergence, aphids seem to have expanded the ancestral hemipteran DNA methylation along with a distinctive linkage to the histone methylation system, as suggested by expansion of SET domain methylases, including those fused to methylated CpG recognition domains. Thus, differential use of DNA methylation and histone methylation might have played a role in emergence of polyphenism and cyclic parthenogenesis from the ancestral hemipteran.

Evolution of sexes from an ancestral mating-type specification pathway

Male and female sexes have evolved repeatedly in eukaryotes but the origins of dimorphic sexes and their relationship to mating types in unicellular species are not understood. Volvocine algae include isogamous species such as Chlamydomonas reinhardtii, with two equal-sized mating types, and oogamous multicellular species such as Volvox carteri with sperm-producing males and egg-producing females. Theoretical work predicts genetic linkage of a gamete cell-size regulatory gene(s) to an ancestral mating-type locus as a possible step in the evolution of dimorphic gametes, but this idea has not been tested. Here we show that, contrary to predictions, a single conserved mating locus (MT) gene in volvocine algae-MID, which encodes a RWP-RK domain transcription factor-evolved from its ancestral role in C. reinhardtii as a mating-type specifier, to become a determinant of sperm and egg development in V. carteri. Transgenic female V. carteri expressing male MID produced functional sperm packets during sexual development. Transgenic male V. carteri with RNA interference (RNAi)-mediated knockdowns of VcMID produced functional eggs, or self-fertile hermaphrodites. Post-transcriptional controls were found to regulate cell-type-limited expression and nuclear localization of VcMid protein that restricted its activity to nuclei of developing male germ cells and sperm. Crosses with sex-reversed strains uncoupled sex determination from sex chromosome identity and revealed gender-specific roles for male and female mating locus genes in sexual development, gamete fitness and reproductive success. Our data show genetic continuity between the mating-type specification and sex determination pathways of volvocine algae, and reveal evidence for gender-specific adaptations in the male and female mating locus haplotypes of Volvox. These findings will enable a deeper understanding of how a master regulator of mating-type determination in an ancestral unicellular species was reprogrammed to control sexually dimorphic gamete development in a multicellular descendant.

New perspectives on the diversification of the RNA interference system: insights from comparative genomics and small RNA sequencing

Our understanding of the pervasive involvement of small RNAs in regulating diverse biological processes has been greatly augmented by recent application of deep-sequencing technologies to small RNA across diverse eukaryotes. We review the currently known small RNA classes and place them in context of the reconstructed evolutionary history of the RNA interference (RNAi) protein machinery. This synthesis indicates that the earliest versions of eukaryotic RNAi systems likely utilized small RNA processed from three types of precursors: (1) sense-antisense transcriptional products, (2) genome-encoded, imperfectly complementary hairpin sequences, and (3) larger noncoding RNA precursor sequences. Structural dissection of PIWI proteins along with recent discovery of novel families (including Med13 of the Mediator complex) suggest that emergence of a distinct architecture with the N-terminal domains (also occurring separately fused to endoDNases in prokaryotes) formed via duplication of an ancestral unit was key to their recruitment as primary RNAi effectors and use of small RNAs of certain preferred lengths. Prokaryotic PIWI proteins are typically components of several RNA-directed DNA restriction or CRISPR/Cas systems. However, eukaryotic versions appear to have emerged from a subset that evolved RNA-directed RNAi. They were recruited alongside RNaseIII domains and RNA-dependent RNA polymerase (RdRP) domains, also from prokaryotic systems, to form the core eukaryotic RNAi system. Like certain regulatory systems, RNAi diversified into two distinct but linked arms concomitant with eukaryotic nucleocytoplasmic compartmentalization. Subsequent elaboration of RNAi proceeded via diversification of the core protein machinery through lineage-specific expansions and recruitment of new components from prokaryotes (nucleases and small RNA-modifying enzymes), allowing for diversification of associating small RNAs.

Impact of whole-genome and tandem duplications in the expansion and functional diversification of the F-box family in legumes (Fabaceae)

F-box proteins constitute a large gene family that regulates processes from hormone signaling to stress response. F-box proteins are the substrate recognition modules of SCF E3 ubiquitin ligases. Here we report very distinct trends in family size, duplication, synteny and transcription of F-box genes in two nitrogen-fixing legumes, Glycine max (soybean) and Medicago truncatula (alfafa). While the soybean FBX genes emerged mainly through segmental duplications (including whole-genome duplications), M. truncatula genome is dominated by locally-duplicated (tandem) F-box genes. Many of these young FBX genes evolved complex transcriptional patterns, including preferential transcription in different tissues, suggesting that they have probably been recruited to important biochemical pathways (e.g. nodulation and seed development).

Structural and functional characterization of the protein kinase Mps1 in Arabidopsis thaliana

In eukaryotes, protein kinases catalyze the transfer of a gamma-phosphate from ATP (or GTP) to specific amino acids in protein targets. In plants, protein kinases have been shown to participate in signaling cascades driving responses to environmental stimuli and developmental processes. Plant meristems are undifferentiated tissues that provide the major source of cells that will form organs throughout development. However, non-dividing specialized cells can also dedifferentiate and re-initiate cell division if exposed to appropriate conditions. Mps1 (Monopolar spindle) is a dual-specificity protein kinase that plays a critical role in monitoring the accuracy of chromosome segregation in the mitotic checkpoint mechanism. Although Mps1 functions have been clearly demonstrated in animals and fungi, its role in plants is so far unclear. Here, using structural and biochemical analyses here we show that Mps1 has highly similar homologs in many plant genomes across distinct lineages (e.g. AtMps1 in Arabidopsis thaliana). Several structural features (i.e. catalytic site, DFG motif and threonine triad) are clearly conserved in plant Mps1 kinases. Structural and sequence analysis also suggest that AtMps1 interact with other cell cycle proteins, such as Mad2 and MAPK1. By using a very specific Mps1 inhibitor (SP600125) we show that compromised AtMps1 activity hampers the development of A. thaliana seedlings in a dose-dependent manner, especially in secondary roots. Moreover, concomitant administration of the auxin IAA neutralizes the AtMps1 inhibition phenotype, allowing secondary root development. These observations let us to hypothesize that AtMps1 might be a downstream regulator of IAA signaling in the formation of secondary roots. Our results indicate that Mps1 might be a universal component of the Spindle Assembly Checkpoint machinery across very distant lineages of eukaryotes.

Gene flow and biological conflict systems in the origin and evolution of eukaryotes

The endosymbiotic origin of eukaryotes brought together two disparate genomes in the cell. Additionally, eukaryotic natural history has included other endosymbiotic events, phagotrophic consumption of organisms, and intimate interactions with viruses and endoparasites. These phenomena facilitated large-scale lateral gene transfer and biological conflicts. We synthesize information from nearly two decades of genomics to illustrate how the interplay between lateral gene transfer and biological conflicts has impacted the emergence of new adaptations in eukaryotes. Using apicomplexans as example, we illustrate how lateral transfer from animals has contributed to unique parasite-host interfaces comprised of adhesion- and O-linked glycosylation-related domains. Adaptations, emerging due to intense selection for diversity in the molecular participants in organismal and genomic conflicts, being dispersed by lateral transfer, were subsequently exapted for eukaryote-specific innovations. We illustrate this using examples relating to eukaryotic chromatin, RNAi and RNA-processing systems, signaling pathways, apoptosis and immunity. We highlight the major contributions from catalytic domains of bacterial toxin systems to the origin of signaling enzymes (e.g., ADP-ribosylation and small molecule messenger synthesis), mutagenic enzymes for immune receptor diversification and RNA-processing. Similarly, we discuss contributions of bacterial antibiotic/siderophore synthesis systems and intra-genomic and intra-cellular selfish elements (e.g., restriction-modification, mobile elements and lysogenic phages) in the emergence of chromatin remodeling/modifying enzymes and RNA-based regulation. We develop the concept that biological conflict systems served as evolutionary “nurseries” for innovations in the protein world, which were delivered to eukaryotes via lateral gene flow to spur key evolutionary innovations all the way from nucleogenesis to lineage-specific adaptations.

Comparative genomics uncovers novel structural and functional features of the heterotrimeric GTPase signaling system

Though the heterotrimeric G-proteins signaling system is one of the best studied in eukaryotes, its provenance and its prevalence outside of model eukaryotes remains poorly understood. We utilized the wealth of sequence data from recently sequenced eukaryotic genomes to uncover robust G-protein signaling systems in several poorly studied eukaryotic lineages such as the parabasalids, heteroloboseans and stramenopiles. This indicated that the Gα subunit is likely to have separated from the ARF-like GTPases prior to the last eukaryotic common ancestor. We systematically identified the structure and sequence features associated with this divergence and found that most of the neomorphic positions in Gα form a ring of residues centered on the nucleotide binding site, several of which are likely to be critical for interactions with the RGS domain for its GAP function. We also present evidence that in some of the potentially early branching eukaryotic lineages, like Trichomonas, Gα is likely to function independently of the Gβγ subunits. We were able to identify previously unknown Gγ subunits in Naegleria, suggesting that the trimeric version was already present by the time of the divergence of the heteroloboseans from the remaining eukaryotes. Evolution of Gα subunits is dominated by several independent lineage-specific expansions (LSEs). In most of these cases there are concomitant, independent LSEs of RGS proteins along with an extraordinary diversification of their domain architectures. The diversity of RGS domains from Naegleria in particular, which has the largest complement of Gα and RGS proteins for any eukaryote, provides new insights into RGS function and evolution. We uncovered a new class of soluble ligand receptors of bacterial origin with RGS domains and an extraordinary diversity of membrane-linked, redox-associated, adhesion-dependent and small molecule-induced G-protein signaling networks that evolved in early-branching eukaryotes, independently of parallel systems in animals. Furthermore, this newly characterized diversity of RGS domains helps in defining their ancestral conserved interfaces with Gα and also those interfaces that are prone to extensive lineage-specific diversification and are thereby responsible for selectivity in Gα-RGS interactions. Several mushrooms show LSEs of Gαs but not of RGS proteins pointing to the probable differentiation of Gαs in conjunction with mating-type diversity. When combined with the characterization of the 7TM receptors (GPCRs), it becomes apparent that, through much of eukaryotic evolution, cells contained both 7TM receptors that acted as GEFs and those as GAPs (with C-terminal RGS domains) for Gαs. Only in some lineages like animals and stramenopiles the 7TM receptors were restricted to GEF only roles, probably due to selection imposed by the rate-constants of the Gαs that underwent lineage-specific expansion in them. In the alveolate lineage the 7TM receptors occur independently of heterotrimeric G-proteins, suggesting the prevalence of G-protein-independent signaling in these organisms.

Dynamics of a novel centromeric histone variant CenH3 reveals the evolutionary ancestral timing of centromere biogenesis

The centromeric histone H3 variant (CenH3) serves to target the kinetochore to the centromeres and thus ensures correct chromosome segregation during mitosis and meiosis. The Dictyostelium H3-like variant H3v1 was identified as the CenH3 ortholog. Dictyostelium CenH3 has an extended N-terminal domain with no similarity to any other known proteins and a histone fold domain at its C-terminus. Within the histone fold, α-helix 2 (α2) and an extended loop 1 (L1) have been shown to be required for targeting CenH3 to centromeres. Compared to other known and putative CenH3 histones, Dictyostelium CenH3 has a shorter L1, suggesting that the extension is not an obligatory feature. Through ChIP analysis and fluorescence microscopy of live and fixed cells, we provide here the first survey of centromere structure in amoebozoa. The six telocentric centromeres were found to mostly consist of all the DIRS-1 elements and to associate with H3K9me3. During interphase, the centromeres remain attached to the centrosome forming a single CenH3-containing cluster. Loading of Dictyostelium CenH3 onto centromeres occurs at the G2/prophase transition, in contrast to the anaphase/telophase loading of CenH3 observed in metazoans. This suggests that loading during G2/prophase is the ancestral eukaryotic mechanism and that anaphase/telophase loading of CenH3 has evolved more recently after the amoebozoa diverged from the animal linage.