Adaptive Evolution and the Birth of CTCF Binding Sites in the Drosophila Genome

Changes in the physical interaction between cis-regulatory DNA sequences and proteins drive the evolution of gene expression. However, it has proven difficult to accurately quantify evolutionary rates of such binding change or to estimate the relative effects of selection and drift in shaping the binding evolution. Here we examine the genome-wide binding of CTCF in four species of Drosophila separated by between ~2.5 and 25 million years. CTCF is a highly conserved protein known to be associated with insulator sequences in the genomes of human and Drosophila. Although the binding preference for CTCF is highly conserved, we find that CTCF binding itself is highly evolutionarily dynamic and has adaptively evolved. Between species, binding divergence increased linearly with evolutionary distance, and CTCF binding profiles are diverging rapidly at the rate of 2.22% per million years (Myr). At least 89 new CTCF binding sites have originated in the Drosophila melanogaster genome since the most recent common ancestor with Drosophila simulans. Comparing these data to genome sequence data from 37 different strains of Drosophila melanogaster, we detected signatures of selection in both newly gained and evolutionarily conserved binding sites. Newly evolved CTCF binding sites show a significantly stronger signature for positive selection than older sites. Comparative gene expression profiling revealed that expression divergence of genes adjacent to CTCF binding site is significantly associated with the gain and loss of CTCF binding. Further, the birth of new genes is associated with the birth of new CTCF binding sites. Our data indicate that binding of Drosophila CTCF protein has evolved under natural selection, and CTCF binding evolution has shaped both the evolution of gene expression and genome evolution during the birth of new genes

Insertion of Horizontally Transferred Genes within Conserved Syntenic Regions of Yeast Genomes

Horizontal gene transfer has been occasionally mentioned in eukaryotic genomes, but such events appear much less numerous than in prokaryotes, where they play important functional and evolutionary roles. In yeasts, few independent cases have been described, some of which corresponding to major metabolic functions, but no systematic screening of horizontally transferred genes has been attempted so far. Taking advantage of the synteny conservation among five newly sequenced and annotated genomes of Saccharomycetaceae, we carried out a systematic search for HGT candidates amidst genes present in only one species within conserved synteny blocks. Out of 255 species-specific genes, we discovered 11 candidates for HGT, based on their similarity with bacterial proteins and on reconstructed phylogenies. This corresponds to a minimum of six transfer events because some horizontally acquired genes appear to rapidly duplicate in yeast genomes (e.g. YwqG genes in Kluyveromyces thermotolerans and serine recombinase genes of the IS607 family in Saccharomyces kluyveri). We show that the resulting copies are submitted to a strong functional selective pressure. The mechanisms of DNA transfer and integration are discussed, in relation with the generally small size of HGT candidates. Our results on a limited set of species expand by 50% the number of previously published HGT cases in hemiascomycetous yeasts, suggesting that this type of event is more frequent than usually thought. Our restrictive method does not exclude the possibility that additional HGT events exist. Actually, ancestral events common to several yeast species must have been overlooked, and the absence of homologs in present databases leaves open the question of the origin of the 244 remaining species-specific genes inserted within conserved synteny blocks

Widespread Site-Dependent Buffering of Human Regulatory Polymorphism

The average individual is expected to harbor thousands of variants within non-coding genomic regions involved in gene regulation. However, it is currently not possible to interpret reliably the functional consequences of genetic variation within any given transcription factor recognition sequence. To address this, we comprehensively analyzed heritable genome-wide binding patterns of a major sequence-specific regulator (CTCF) in relation to genetic variability in binding site sequences across a multi-generational pedigree. We localized and quantified CTCF occupancy by ChIP-seq in 12 related and unrelated individuals spanning three generations, followed by comprehensive targeted resequencing of the entire CTCF–binding landscape across all individuals. We identified hundreds of variants with reproducible quantitative effects on CTCF occupancy (both positive and negative). While these effects paralleled protein–DNA recognition energetics when averaged, they were extensively buffered by striking local context dependencies. In the significant majority of cases buffering was complete, resulting in silent variants spanning every position within the DNA recognition interface irrespective of level of binding energy or evolutionary constraint. The prevalence of complex partial or complete buffering effects severely constrained the ability to predict reliably the impact of variation within any given binding site instance. Surprisingly, 40% of variants that increased CTCF occupancy occurred at positions of human–chimp divergence, challenging the expectation that the vast majority of functional regulatory variants should be deleterious. Our results suggest that, even in the presence of “perfect” genetic information afforded by resequencing and parallel studies in multiple related individuals, genomic site-specific prediction of the consequences of individual variation in regulatory DNA will require systematic coupling with empirical functional genomic measurements

Divergence of the Yeast Transcription Factor FZF1 Affects Sulfite Resistance

Changes in gene expression are commonly observed during evolution. However, the phenotypic consequences of expression divergence are frequently unknown and difficult to measure. Transcriptional regulators provide a mechanism by which phenotypic divergence can occur through multiple, coordinated changes in gene expression during development or in response to environmental changes. Yet, some changes in transcriptional regulators may be constrained by their pleiotropic effects on gene expression. Here, we use a genome-wide screen for promoters that are likely to have diverged in function and identify a yeast transcription factor, FZF1, that has evolved substantial differences in its ability to confer resistance to sulfites. Chimeric alleles from four Saccharomyces species show that divergence in FZF1 activity is due to changes in both its coding and upstream noncoding sequence. Between the two closest species, noncoding changes affect the expression of FZF1, whereas coding changes affect the expression of SSU1, a sulfite efflux pump activated by FZF1. Both coding and noncoding changes also affect the expression of many other genes. Our results show how divergence in the coding and promoter region of a transcription factor alters the response to an environmental stress

Measuring the Evolutionary Rewiring of Biological Networks

We have accumulated a large amount of biological network data and expect even more to come. Soon, we anticipate being able to compare many different biological networks as we commonly do for molecular sequences. It has long been believed that many of these networks change, or “rewire”, at different rates. It is therefore important to develop a framework to quantify the differences between networks in a unified fashion. We developed such a formalism based on analogy to simple models of sequence evolution, and used it to conduct a systematic study of network rewiring on all the currently available biological networks. We found that, similar to sequences, biological networks show a decreased rate of change at large time divergences, because of saturation in potential substitutions. However, different types of biological networks consistently rewire at different rates. Using comparative genomics and proteomics data, we found a consistent ordering of the rewiring rates: transcription regulatory, phosphorylation regulatory, genetic interaction, miRNA regulatory, protein interaction, and metabolic pathway network, from fast to slow. This ordering was found in all comparisons we did of matched networks between organisms. To gain further intuition on network rewiring, we compared our observed rewirings with those obtained from simulation. We also investigated how readily our formalism could be mapped to other network contexts; in particular, we showed how it could be applied to analyze changes in a range of “commonplace” networks such as family trees, co-authorships and linux-kernel function dependencies

Preserving and Using Germplasm and Dissociated Embryonic Cells for Conserving Caribbean and Pacific Coral

Coral reefs are experiencing unprecedented degradation due to human activities, and protecting specific reef habitats may not stop this decline, because the most serious threats are global (i.e., climate change), not local. However, ex situ preservation practices can provide safeguards for coral reef conservation. Specifically, modern advances in cryobiology and genome banking could secure existing species and genetic diversity until genotypes can be introduced into rehabilitated habitats. We assessed the feasibility of recovering viable sperm and embryonic cells post-thaw from two coral species, Acropora palmata and Fungia scutaria that have diffferent evolutionary histories, ecological niches and reproductive strategies. In vitro fertilization (IVF) of conspecific eggs using fresh (control) spermatozoa revealed high levels of fertilization (>90% in A. palmata; >84% in F. scutaria; P>0.05) that were unaffected by tested sperm concentrations. A solution of 10% dimethyl sulfoxide (DMSO) at cooling rates of 20 to 30°C/min most successfully cryopreserved both A. palmata and F. scutaria spermatozoa and allowed producing developing larvae in vitro. IVF success under these conditions was 65% in A. palmata and 53% in F. scutaria on particular nights; however, on subsequent nights, the same process resulted in little or no IVF success. Thus, the window for optimal freezing of high quality spermatozoa was short (∼5 h for one night each spawning cycle). Additionally, cryopreserved F. scutaria embryonic cells had∼50% post-thaw viability as measured by intact membranes. Thus, despite some differences between species, coral spermatozoa and embryonic cells are viable after low temperature (−196°C) storage, preservation and thawing. Based on these results, we have begun systematically banking coral spermatozoa and embryonic cells on a large-scale as a support approach for preserving existing bio- and genetic diversity found in reef systems

Identification and Genome-Wide Prediction of DNA Binding Specificities for the ApiAP2 Family of Regulators from the Malaria Parasite

The molecular mechanisms underlying transcriptional regulation in apicomplexan parasites remain poorly understood. Recently, the Apicomplexan AP2 (ApiAP2) family of DNA binding proteins was identified as a major class of transcriptional regulators that are found across all Apicomplexa. To gain insight into the regulatory role of these proteins in the malaria parasite, we have comprehensively surveyed the DNA-binding specificities of all 27 members of the ApiAP2 protein family from Plasmodium falciparum revealing unique binding preferences for the majority of these DNA binding proteins. In addition to high affinity primary motif interactions, we also observe interactions with secondary motifs. The ability of a number of ApiAP2 proteins to bind multiple, distinct motifs significantly increases the potential complexity of the transcriptional regulatory networks governed by the ApiAP2 family. Using these newly identified sequence motifs, we infer the trans-factors associated with previously reported plasmodial cis-elements and provide evidence that ApiAP2 proteins modulate key regulatory decisions at all stages of parasite development. Our results offer a detailed view of ApiAP2 DNA binding specificity and take the first step toward inferring comprehensive gene regulatory networks for P. falciparum

Transcription Factors Mat2 and Znf2 Operate Cellular Circuits Orchestrating Opposite- and Same-Sex Mating in Cryptococcus neoformans

Cryptococcus neoformans is a human fungal pathogen that undergoes a dimorphic transition from a unicellular yeast to multicellular hyphae during opposite sex (mating) and unisexual reproduction (same-sex mating). Opposite- and same-sex mating are induced by similar environmental conditions and involve many shared components, including the conserved pheromone sensing Cpk1 MAPK signal transduction cascade that governs the dimorphic switch in C. neoformans. However, the homeodomain cell identity proteins Sxi1α/Sxi2a encoded by the mating type locus that are essential for completion of sexual reproduction following cell–cell fusion during opposite-sex mating are dispensable for same-sex mating. Therefore, identification of downstream targets of the Cpk1 MAPK pathway holds the key to understanding molecular mechanisms governing the two distinct developmental fates. Thus far, homology-based approaches failed to identify downstream transcription factors which may therefore be species-specific. Here, we applied insertional mutagenesis via Agrobacterium-mediated transformation and transcription analysis using whole genome microarrays to identify factors involved in C. neoformans differentiation. Two transcription factors, Mat2 and Znf2, were identified as key regulators of hyphal growth during same- and opposite-sex mating. Mat2 is an HMG domain factor, and Znf2 is a zinc finger protein; neither is encoded by the mating type locus. Genetic, phenotypic, and transcriptional analyses of Mat2 and Znf2 provide evidence that Mat2 is a downstream transcription factor of the Cpk1 MAPK pathway whereas Znf2 functions as a more terminal hyphal morphogenesis determinant. Although the components of the MAPK pathway including Mat2 are not required for virulence in animal models, Znf2, as a hyphal morphology determinant, is a negative regulator of virulence. Further characterization of these elements and their target circuits will reveal genes controlling biological processes central to fungal development and virulence

Mineral Type and Solution Chemistry Affect the Structure and Composition of Actively Growing Bacterial Communities as Revealed by Bromodeoxyuridine Immunocapture and 16S rRNA Pyrosequencing

© 2016, Springer Science+Business Media New York. Understanding how minerals affect bacterial communities and their in situ activities in relation to environmental conditions are central issues in soil microbial ecology, as minerals represent essential reservoirs of inorganic nutrients for the biosphere. To determine the impact of mineral type and solution chemistry on soil bacterial communities, we compared the diversity, composition, and functional abilities of a soil bacterial community incubated in presence/absence of different mineral types (apatite, biotite, obsidian). Microcosms were prepared containing different liquid culture media devoid of particular essential nutrients, the nutrients provided only in the introduced minerals and therefore only available to the microbial community through mineral dissolution by biotic and/or abiotic processes. By combining functional screening of bacterial isolates and community analysis by bromodeoxyuridine DNA immunocapture and 16S rRNA gene pyrosequencing, we demonstrated that bacterial communities were mainly impacted by the solution chemistry at the taxonomic level and by the mineral type at the functional level. Metabolically active bacterial communities varied with solution chemistry and mineral type. Burkholderia were significantly enriched in the obsidian treatment compared to the biotite treatment and were the most effective isolates at solubilizing phosphorous or mobilizing iron, in all the treatments. A detailed analysis revealed that the 16S rRNA gene sequences of the OTUs or isolated strains assigned as Burkholderia in our study showed high homology with effective mineral-weathering bacteria previously recovered from the same experimental site