ALBERT

Description: Some of the most dangerous pathogens such as Mycobacterium tuberculosis and Yersinia pestis evolve clonally . This means that little or no recombination occurs between strains belonging to these species. Paradoxically, although different members of these species show extreme sequence similarity of orthologous genes, some show considerable intraspecies phenotypic variation, the source of which remains elusive. To examine the possible sources of phenotypic variation within clonal pathogenic bacterial species, we carried out an extensive genomic and pan-genomic analysis of the sources of genetic variation available to a large collection of clonal and nonclonal pathogenic bacterial species. We show that while nonclonal species diversify through a combination of changes to gene sequences, gene loss and gene gain, gene loss completely dominates as a source of genetic variation within clonal species. Indeed, gene loss is so prevalent within clonal species as to lead to levels of gene content variation comparable to those found in some nonclonal species that are much more diverged in their gene sequences and that acquire a substantial number of genes horizontally. Gene loss therefore needs to be taken into account as a potential dominant source of phenotypic variation within clonal bacterial species.

Description: Obligate bacterial symbionts are widespread in many invertebrates, where they are often confined to specialized host cells and are transmitted directly from mother to progeny. Increasing numbers of these bacteria are being characterized but questions remain about their population structure and evolution. Here we take a comparative genomics approach to investigate two prominent bacterial symbionts (BFo1 and BFo2) isolated from geographically separated populations of western flower thrips, Frankliniella occidentalis. Our multifaceted approach to classifying these symbionts includes concatenated multilocus sequence analysis (MLSA) phylogenies, ribosomal multilocus sequence typing (rMLST), construction of whole-genome phylogenies, and in-depth genomic comparisons. We showed that the BFo1 genome clusters more closely to species in the genus Erwinia, and is a putative close relative to Erwinia aphidicola . BFo1 is also likely to have shared a common ancestor with Erwinia pyrifoliae/Erwinia amylovora and the nonpathogenic Erwinia tasmaniensis and genetic traits similar to Erwinia billingiae . The BFo1 genome contained virulence factors found in the genus Erwinia but represented a divergent lineage. In contrast, we showed that BFo2 belongs within the Enterobacteriales but does not group closely with any currently known bacterial species. Concatenated MLSA phylogenies indicate that it may have shared a common ancestor to the Erwinia and Pantoea genera, and based on the clustering of rMLST genes, it was most closely related to Pantoea ananatis but represented a divergent lineage. We reconstructed a core genome of a putative common ancestor of Erwinia and Pantoea and compared this with the genomes of BFo bacteria. BFo2 possessed none of the virulence determinants that were omnipresent in the Erwinia and Pantoea genera. Taken together, these data are consistent with BFo2 representing a highly novel species that maybe related to known Pantoea .

Description: Gene expression evolution occurs through changes in cis - or trans -regulatory elements or both. Interactions between transcription factors (TFs) and their binding sites (TFBSs) constitute one of the most important points where these two regulatory components intersect. In this study, we investigated the evolution of TFBSs in the promoter regions of different Saccharomyces strains and species. We divided the promoter of a gene into the proximal region and the distal region, which are defined, respectively, as the 200-bp region upstream of the transcription starting site and as the 200-bp region upstream of the proximal region. We found that the predicted TFBSs in the proximal promoter regions tend to be evolutionarily more conserved than those in the distal promoter regions. Additionally, Saccharomyces cerevisiae strains used in the fermentation of alcoholic drinks have experienced more TFBS losses than gains compared with strains from other environments (wild strains, laboratory strains, and clinical strains). We also showed that differences in TFBSs correlate with the cis component of gene expression evolution between species (comparing S. cerevisiae and its sister species Saccharomyces paradoxus ) and within species (comparing two closely related S. cerevisiae strains).

Description: Gene duplication is a key factor contributing to phenotype diversity across and within species. Although the availability of complete genomes has led to the extensive study of genomic duplications, the dynamics and variability of gene duplications mediated by retrotransposition are not well understood. Here, we predict mRNA retrotransposition and use comparative genomics to investigate their origin and variability across primates. Analyzing seven anthropoid primate genomes, we found a similar number of mRNA retrotranspositions (~7,500 retrocopies) in Catarrhini (Old Word Monkeys, including humans), but a surprising large number of retrocopies (~10,000) in Platyrrhini (New World Monkeys), which may be a by-product of higher long interspersed nuclear element 1 activity in these genomes. By inferring retrocopy orthology, we dated most of the primate retrocopy origins, and estimated a decrease in the fixation rate in recent primate history, implying a smaller number of species-specific retrocopies. Moreover, using RNA-Seq data, we identified approximately 3,600 expressed retrocopies. As expected, most of these retrocopies are located near or within known genes, present tissue-specific and even species-specific expression patterns, and no expression correlation to their parental genes. Taken together, our results provide further evidence that mRNA retrotransposition is an active mechanism in primate evolution and suggest that retrocopies may not only introduce great genetic variability between lineages but also create a large reservoir of potentially functional new genomic loci in primate genomes.

Description: Viruses rely completely on the hosts’ machinery for translation of viral transcripts. However, for most viruses infecting humans, codon usage preferences (CUPrefs) do not match those of the host. Human papillomaviruses (HPVs) are a showcase to tackle this paradox: they present a large genotypic diversity and a broad range of phenotypic presentations, from asymptomatic infections to productive lesions and cancer. By applying phylogenetic inference and dimensionality reduction methods, we demonstrate first that genes in HPVs are poorly adapted to the average human CUPrefs, the only exception being capsid genes in viruses causing productive lesions. Phylogenetic relationships between HPVs explained only a small proportion of CUPrefs variation. Instead, the most important explanatory factor for viral CUPrefs was infection phenotype, as orthologous genes in viruses with similar clinical presentation displayed similar CUPrefs. Moreover, viral genes with similar spatiotemporal expression patterns also showed similar CUPrefs. Our results suggest that CUPrefs in HPVs reflect either variations in the mutation bias or differential selection pressures depending on the clinical presentation and expression timing. We propose that poor viral CUPrefs may be central to a trade-off between strong viral gene expression and the potential for eliciting protective immune response.

Description: Spliceosomal introns are a hallmark of eukaryotic genes that are hypothesized to play important roles in genome evolution but have poorly understood origins. Although most introns lack sequence homology to each other, new families of spliceosomal introns that are repeated hundreds of times in individual genomes have recently been discovered in a few organisms. The prevalence and conservation of these introner elements (IEs) or introner-like elements in other taxa, as well as their evolutionary relationships to regular spliceosomal introns, are still unknown. Here, we systematically investigate introns in the widespread marine green alga Micromonas and report new families of IEs, numerous intron presence–absence polymorphisms, and potential intron insertion hot-spots. The new families enabled identification of conserved IE secondary structure features and establishment of a novel general model for repetitive intron proliferation across genomes. Despite shared secondary structure, the IE families from each Micromonas lineage bear no obvious sequence similarity to those in the other lineages, suggesting that their appearance is intimately linked with the process of speciation. Two of the new IE families come from an Arctic culture ( Micromonas Clade E2) isolated from a polar region where abundance of this alga is increasing due to climate induced changes. The same two families were detected in metagenomic data from Antarctica—a system where Micromonas has never before been reported. Strikingly high identity between the Arctic isolate and Antarctic coding sequences that flank the IEs suggests connectivity between populations in the two polar systems that we postulate occurs through deep-sea currents. Recovery of Clade E2 sequences in North Atlantic Deep Waters beneath the Gulf Stream supports this hypothesis. Our research illuminates the dynamic relationships between an unusual class of repetitive introns, genome evolution, speciation, and global distribution of this sentinel marine alga.

Description: Cellular Ca 2+ homeostasis is tightly regulated and is pivotal to life. Inositol 1,4,5-trisphosphate receptors (IP 3 Rs) and ryanodine receptors (RyRs) are the major ion channels that regulate Ca 2+ release from intracellular stores. Although these channels have been extensively investigated in multicellular organisms, an appreciation of their evolution and the biology of orthologs in unicellular organisms is largely lacking. Extensive phylogenetic analyses reveal that the IP 3 R gene superfamily is ancient and diverged into two subfamilies, IP 3 R-A and IP 3 R-B/RyR, at the dawn of Opisthokonta. IP 3 R-B/RyR further diversified into IP 3 R-B and RyR at the stem of Filozoa. Subsequent evolution and speciation of Holozoa is associated with duplication of IP 3 R-A and RyR genes, and loss of IP 3 R-B in the vertebrate lineages. To gain insight into the properties of IP 3 R important for the challenges of multicellularity, the IP 3 R-A and IP 3 R-B family orthologs were cloned from Capsaspora owczarzaki, a close unicellular relative to Metazoa (designated as CO.IP 3 R-A and CO.IP 3 R-B). Both proteins were targeted to the endoplasmic reticulum. However, CO.IP 3 R-A, but strikingly not CO.IP 3 R-B, bound IP 3 , exhibited robust Ca 2+ release activity and associated with mammalian IP 3 Rs. These data indicate strongly that CO.IP 3 R-A as an exemplar of ancestral IP 3 R-A orthologs forms bona fide IP 3 -gated channels. Notably, however, CO.IP 3 R-A appears not to be regulated by Ca 2+ , ATP or Protein kinase A-phosphorylation. Collectively, our findings explore the origin, conservation, and diversification of IP 3 R gene families and provide insight into the functionality of ancestral IP 3 Rs and the added specialization of these proteins in Metazoa.

Description: Marriage rules, the community prescriptions that dictate who an individual can or cannot marry, are extremely diverse and universally present in traditional societies. A major focus of research in the early decades of modern anthropology, marriage rules impose social and economic forces that help structure societies and forge connections between them. However, in those early anthropological studies, the biological benefits or disadvantages of marriage rules could not be determined. We revisit this question by applying a novel simulation framework and genome-wide data to explore the effects of Asymmetric Prescriptive Alliance, an elaborate set of marriage rules that has been a focus of research for many anthropologists. Simulations show that strict adherence to these marriage rules reduces genetic diversity on the autosomes, X chromosome and mitochondrial DNA, but relaxed compliance produces genetic diversity similar to random mating. Genome-wide data from the Indonesian community of Rindi, one of the early study populations for Asymmetric Prescriptive Alliance, are more consistent with relaxed compliance than strict adherence. We therefore suggest that, in practice, marriage rules are treated with sufficient flexibility to allow social connectivity without significant degradation of biological diversity.

Description: The environment has profound effects on the expression of many traits and reaction norms describe the expression dynamics of a trait across a broad range of environmental conditions. Here, we analyze gene expression in Drosophila melanogaster across four different developmental temperatures (13–29 °C). Gene expression is highly plastic with 83.3% of the genes being differentially expressed. We distinguished three components of plasticity: 1) Dynamics of gene expression intensity (sum of change), 2) direction of change, and 3) curvature of the reaction norm (linear vs. quadratic). Studying their regulatory architecture we found that all three plasticity components were most strongly affected by the number of different transcription factors (TFs) binding to the target gene. More TFs were found in genes with less expression changes across temperatures. Although the effect of microRNAs was weaker, we consistently noted a trend in the opposite direction. The most plastic genes were regulated by fewer TFs and more microRNAs than less plastic genes. Different patterns of plasticity were also reflected by their functional characterization based on gene ontology. Our results suggest that reaction norms provide an important key to understand the functional requirements of natural populations exposed to variable environmental conditions.

Description: Recent advances in paleogenomic technologies have enabled an increasingly detailed understanding of the evolutionary relationships of now-extinct mammalian taxa. However, a number of enigmatic Quaternary species have never been characterized with molecular data, often because available fossils are rare or are found in environments that are not optimal for DNA preservation. Here, we analyze paleogenomic data extracted from bones attributed to the late Pleistocene western camel, Camelops cf. hesternus, a species that was distributed across central and western North America until its extinction approximately 13,000 years ago. Despite a modal sequence length of only around 35 base pairs, we reconstructed high-coverage complete mitochondrial genomes and low-coverage partial nuclear genomes for each specimen. We find that Camelops is sister to African and Asian bactrian and dromedary camels, to the exclusion of South American camelids (llamas, guanacos, alpacas, and vicuñas). These results contradict previous morphology-based phylogenetic models for Camelops , which suggest instead a closer relationship between Camelops and the South American camelids. The molecular data imply a Late Miocene divergence of the Camelops clade from lineages that separately gave rise to the extant camels of Eurasia. Our results demonstrate the increasing capacity of modern paleogenomic methods to resolve evolutionary relationships among distantly related lineages.