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: Historically, genome-wide and molecular characterization of the genus Listeria has concentrated on the important human pathogen Listeria monocytogenes and a small number of closely related species, together termed Listeria sensu strictu. More recently, a number of genome sequences for more basal, and nonpathogenic, members of the Listeria genus have become available, facilitating a wider perspective on the evolution of pathogenicity and genome level evolutionary dynamics within the entire genus (termed Listeria sensu lato). Here, we have sequenced the genomes of additional Listeria fleischmannii and Listeria newyorkensis isolates and explored the dynamics of genome evolution in Listeria sensu lato. Our analyses suggest that acquisition of genetic material through gene duplication and divergence as well as through lateral gene transfer (mostly from outside Listeria ) is widespread throughout the genus. Novel genetic material is apparently subject to rapid turnover. Multiple lines of evidence point to significant differences in evolutionary dynamics between the most basal Listeria subclade and all other congeners, including both sensu strictu and other sensu lato isolates. Strikingly, these differences are likely attributable to stochastic, population-level processes and contribute to observed variation in genome size across the genus. Notably, our analyses indicate that the common ancestor of Listeria sensu lato lacked flagella, which were acquired by lateral gene transfer by a common ancestor of Listeria grayi and Listeria sensu strictu, whereas a recently functionally characterized pathogenicity island, responsible for the capacity to produce cobalamin and utilize ethanolamine/propane-2-diol, was acquired in an ancestor of Listeria sensu strictu.

Description: Taeniid cestodes (including the human parasites Echinococcus spp. and Taenia solium ) have very few mobile genetic elements (MGEs) in their genome, despite lacking a canonical PIWI pathway. The MGEs of these parasites are virtually unexplored, and nothing is known about their expression and silencing. In this work, we report the discovery of a novel family of small nonautonomous long terminal repeat retrotransposons (also known as terminal-repeat retrotransposons in miniature, TRIMs) which we have named ta-TRIM (taeniid TRIM). ta-TRIM s are only the second family of TRIM elements discovered in animals, and are likely the result of convergent reductive evolution in different taxonomic groups. These elements originated at the base of the taeniid tree and have expanded during taeniid diversification, including after the divergence of closely related species such as Echinococcus multilocularis and Echinococcus granulosus . They are massively expressed in larval stages, from a small proportion of full-length copies and from isolated terminal repeats that show transcriptional read-through into downstream regions, generating novel noncoding RNAs and transcriptional fusions to coding genes. In E. multilocularis , ta-TRIM s are specifically expressed in the germinative cells (the somatic stem cells) during asexual reproduction of metacestode larvae. This would provide a developmental mechanism for insertion of ta-TRIM s into cells that will eventually generate the adult germ line. Future studies of active and inactive ta-TRIM elements could give the first clues on MGE silencing mechanisms in cestodes.

Description: Even though mitochondrial genomes, which characterize eukaryotic cells, were first discovered more than 50 years ago, mitochondrial genomics remains an important topic in molecular biology and genome sciences. The Phylum Alveolata comprises three major groups (ciliates, apicomplexans, and dinoflagellates), the mitochondrial genomes of which have diverged widely. Even though the gene content of dinoflagellate mitochondrial genomes is reportedly comparable to that of apicomplexans, the highly fragmented and rearranged genome structures of dinoflagellates have frustrated whole genomic analysis. Consequently, noncoding sequences and gene arrangements of dinoflagellate mitochondrial genomes have not been well characterized. Here we report that the continuous assembled genome (~326 kb) of the dinoflagellate, Symbiodinium minutum , is AT-rich (~64.3%) and that it contains three protein-coding genes. Based upon in silico analysis, the remaining 99% of the genome comprises transcriptomic noncoding sequences. RNA edited sites and unique, possible start and stop codons clarify conserved regions among dinoflagellates. Our massive transcriptome analysis shows that almost all regions of the genome are transcribed, including 27 possible fragmented ribosomal RNA genes and 12 uncharacterized small RNAs that are similar to mitochondrial RNA genes of the malarial parasite, Plasmodium falciparum . Gene map comparisons show that gene order is only slightly conserved between S. minutu m and P. falciparum . However, small RNAs and intergenic sequences share sequence similarities with P. falciparum , suggesting that the function of noncoding sequences has been preserved despite development of very different genome structures.

Description: Evolutionary studies usually use a two-step process to investigate sequence data. Step one estimates a multiple sequence alignment (MSA) and step two applies phylogenetic methods to ask evolutionary questions of that MSA. Modern phylogenetic methods infer evolutionary parameters using maximum likelihood or Bayesian inference, mediated by a probabilistic substitution model that describes sequence change over a tree. The statistical properties of these methods mean that more data directly translates to an increased confidence in downstream results, providing the substitution model is adequate and the MSA is correct. Many studies have investigated the robustness of phylogenetic methods in the presence of substitution model misspecification, but few have examined the statistical properties of those methods when the MSA is unknown. This simulation study examines the statistical properties of the complete two-step process when inferring sequence divergence and the phylogenetic tree topology. Both nucleotide and amino acid analyses are negatively affected by the alignment step, both through inaccurate guide tree estimates and through overfitting to that guide tree. For many alignment tools these effects become more pronounced when additional sequences are added to the analysis. Nucleotide sequences are particularly susceptible, with MSA errors leading to statistical support for long-branch attraction artifacts, which are usually associated with gross substitution model misspecification. Amino acid MSAs are more robust, but do tend to arbitrarily resolve multifurcations in favor of the guide tree. No inference strategies produce consistently accurate estimates of divergence between sequences, although amino acid MSAs are again more accurate than their nucleotide counterparts. We conclude with some practical suggestions about how to limit the effect of MSA uncertainty on evolutionary inference.

Description: The evolution of mitochondrial information processing pathways, including replication, transcription and translation, is characterized by the gradual replacement of mitochondrial-encoded proteins with nuclear-encoded counterparts of diverse evolutionary origins. Although the ancestral enzymes involved in mitochondrial transcription and replication have been replaced early in eukaryotic evolution, mitochondrial translation is still carried out by an apparatus largely inherited from the α-proteobacterial ancestor. However, variation in the complement of mitochondrial-encoded molecules involved in translation, including transfer RNAs (tRNAs), provides evidence for the ongoing evolution of mitochondrial protein synthesis. Here, we investigate the evolution of the mitochondrial translational machinery using recent genomic and transcriptomic data from animals that have experienced the loss of mt-tRNAs, including phyla Cnidaria and Ctenophora, as well as some representatives of all four classes of Porifera. We focus on four sets of mitochondrial enzymes that directly interact with tRNAs: Aminoacyl-tRNA synthetases, glutamyl-tRNA amidotransferase, tRNA Ile lysidine synthetase, and RNase P. Our results support the observation that the fate of nuclear-encoded mitochondrial proteins is influenced by the evolution of molecules encoded in mitochondrial DNA, but in a more complex manner than appreciated previously. The data also suggest that relaxed selection on mitochondrial translation rather than coevolution between mitochondrial and nuclear subunits is responsible for elevated rates of evolution in mitochondrial translational proteins.

Description: The expansion of DUF1220 domain copy number during human evolution is a dramatic example of rapid and repeated domain duplication. Although patterns of expression, homology, and disease associations suggest a role in cortical development, this hypothesis has not been robustly tested using phylogenetic methods. Here, we estimate DUF1220 domain counts across 12 primate genomes using a nucleotide Hidden Markov Model. We then test a series of hypotheses designed to examine the potential evolutionary significance of DUF1220 copy number expansion. Our results suggest a robust association with brain size, and more specifically neocortex volume. In contradiction to previous hypotheses, we find a strong association with postnatal brain development but not with prenatal brain development. Our results provide further evidence of a conserved association between specific loci and brain size across primates, suggesting that human brain evolution may have occurred through a continuation of existing processes.

Description: Alternative splicing and gene duplication are the two main processes responsible for expanding protein functional diversity. Although gene duplication can generate new genes and alternative splicing can introduce variation through alternative gene products, the interplay between the two processes is complex and poorly understood. Here, we have carried out a study of the evolution of alternatively spliced exons after gene duplication to better understand the interaction between the two processes. We created a manually curated set of 97 human genes with mutually exclusively spliced homologous exons and analyzed the evolution of these exons across five distantly related vertebrates (lamprey, spotted gar, zebrafish, fugu, and coelacanth). Most of these exons had an ancient origin (more than 400 Ma). We found examples supporting two extreme evolutionary models for the behaviour of homologous axons after gene duplication. We observed 11 events in which gene duplication was accompanied by splice isoform separation, that is, each paralog specifically conserved just one distinct ancestral homologous exon. At other extreme, we identified genes in which the homologous exons were always conserved within paralogs, suggesting that the alternative splicing event cannot easily be separated from the function in these genes. That many homologous exons fall in between these two extremes highlights the diversity of biological systems and suggests that the subtle balance between alternative splicing and gene duplication is adjusted to the specific cellular context of each gene.

Description: There is widespread interest today in understanding enhancers, which are regulatory elements typically harboring several transcription factor binding sites and mediating the combinatorial effect of transcription factors on gene expression. The evolution of enhancers poses interesting unanswered questions, for example, the evolutionary time taken for a typical enhancer to emerge or the factors shaping its evolution. Existing approaches to cis -regulatory evolution have often ignored the combinatorial nature and varied biochemical mechanisms of gene regulation encoded in enhancers. We report on our investigation of enhancer evolution through the use of PEBCRES, a framework for evolutionary simulation of enhancers that employs a mechanistic and well-supported sequence-to-expression model to assign fitness to the evolving enhancer genotype. We estimated the time necessary to evolve, from genomic background, enhancers capable of driving complex gene expression patterns similar to those involved in early development in Drosophila. We found the time-to-evolve to range between 0.5 and 10 Myr, and to vary greatly with the target expression pattern, complexity of the real enhancer known to encode that pattern, and the strength of input from specific transcription factors. To our knowledge, this is the first estimate of waiting times for realistic enhancers to evolve. The in silico evolved enhancers had, with a few interesting exceptions, site compositions similar to those seen in real enhancers for the same patterns. Our simulations also revealed that certain features of an enhancer might evolve not due to their biological function but as aids to the evolutionary process itself.

Description: Despite significant progress in the structural and functional characterization of the human genome, understanding of the mechanisms underlying the genetic basis of human phenotypic uniqueness remains limited. Here, I report that transposable element-derived sequences, most notably LTR7/HERV-H, LTR5_Hs, and L1HS, harbor 99.8% of the candidate human-specific regulatory loci (HSRL) with putative transcription factor-binding sites in the genome of human embryonic stem cells (hESC). A total of 4,094 candidate HSRL display selective and site-specific binding of critical regulators (NANOG [Nanog homeobox], POU5F1 [POU class 5 homeobox 1], CCCTC-binding factor [CTCF], Lamin B1), and are preferentially located within the matrix of transcriptionally active DNA segments that are hypermethylated in hESC. hESC-specific NANOG-binding sites are enriched near the protein-coding genes regulating brain size, pluripotency long noncoding RNAs, hESC enhancers, and 5-hydroxymethylcytosine-harboring regions immediately adjacent to binding sites. Sequences of only 4.3% of hESC-specific NANOG-binding sites are present in Neanderthals’ genome, suggesting that a majority of these regulatory elements emerged in Modern Humans. Comparisons of estimated creation rates of novel TF-binding sites revealed that there was 49.7-fold acceleration of creation rates of NANOG-binding sites in genomes of Chimpanzees compared with the mouse genomes and further 5.7-fold acceleration in genomes of Modern Humans compared with the Chimpanzees genomes. Preliminary estimates suggest that emergence of one novel NANOG-binding site detectable in hESC required 466 years of evolution. Pathway analysis of coding genes that have hESC-specific NANOG-binding sites within gene bodies or near gene boundaries revealed their association with physiological development and functions of nervous and cardiovascular systems, embryonic development, behavior, as well as development of a diverse spectrum of pathological conditions such as cancer, diseases of cardiovascular and reproductive systems, metabolic diseases, multiple neurological and psychological disorders. A proximity placement model is proposed explaining how a 33–47% excess of NANOG, CTCF, and POU5F1 proteins immobilized on a DNA scaffold may play a functional role at distal regulatory elements.

Description: Organisms can adapt to local environmental conditions as a plastic response or become adapted through natural selection on genetic variation. The ability to adapt to increased water temperatures will be of paramount importance for many fish species as the climate continues to warm and water resources become limited. Because increased water temperatures will reduce the dissolved oxygen available for fish, we hypothesized that adaptation to low oxygen environments would involve improved respiration through oxidative phosphorylation (OXPHOS). To test this hypothesis, we subjected individuals from two ecologically divergent populations of inland (redband) rainbow trout ( Oncorhynchus mykiss gairdneri ) with historically different temperature regimes (desert and montane) and their F1 progeny to diel cycles of temperature stress and then examined gene expression data for 80 nuclear- and mitochondrial-encoded OXPHOS subunits that participate in respiration. Of the 80 transcripts, 7 showed ≥ 2-fold difference in expression levels in gill tissue from desert fish under heat stress whereas the montane fish had none and the F1 only had one differentially expressed gene. A structural analysis of the proteins encoded by those genes suggests that the response could coordinate the formation of supercomplexes and oligomers. Supercomplexes may increase the efficiency of respiration because complexes I, III, and IV are brought into close proximity and oligomerization of complex V alters the macrostructure of mitochondria to improve respiration. Significant differences in gene expression patterns in response to heat stress in a common environment indicate that the response was not due to plasticity but had a genetic basis.

Description: The enigmatic monocot family Triuridaceae provides a potentially useful model system for studying the effects of an ancient loss of photosynthesis on the plant plastid genome, as all of its members are mycoheterotrophic and achlorophyllous. However, few studies have placed the family in a comparative context, and its phylogenetic placement is only partly resolved. It was also unclear whether any taxa in this family have retained a plastid genome. Here, we used genome survey sequencing to retrieve plastid genome data for Sciaphila densiflora (Triuridaceae) and ten autotrophic relatives in the orders Dioscoreales and Pandanales. We recovered a highly reduced plastome for Sciaphila that is nearly colinear with Carludovica palmata , a photosynthetic relative that belongs to its sister group in Pandanales, Cyclanthaceae–Pandanaceae. This phylogenetic placement is well supported and robust to a broad range of analytical assumptions in maximum-likelihood inference, and is congruent with recent findings based on nuclear and mitochondrial evidence. The 28 genes retained in the S. densiflora plastid genome are involved in translation and other nonphotosynthetic functions, and we demonstrate that nearly all of the 18 protein-coding genes are under strong purifying selection. Our study confirms the utility of whole plastid genome data in phylogenetic studies of highly modified heterotrophic plants, even when they have substantially elevated rates of substitution.

Description: We report the identification of a novel gene family (named MgCRP-I) encoding short secreted cysteine-rich peptides in the Mediterranean mussel Mytilus galloprovincialis . These peptides display a highly conserved pre-pro region and a hypervariable mature peptide comprising six invariant cysteine residues arranged in three intramolecular disulfide bridges. Although their cysteine pattern is similar to cysteines-rich neurotoxic peptides of distantly related protostomes such as cone snails and arachnids, the different organization of the disulfide bridges observed in synthetic peptides and phylogenetic analyses revealed MgCRP-I as a novel protein family. Genome- and transcriptome-wide searches for orthologous sequences in other bivalve species indicated the unique presence of this gene family in Mytilus spp. Like many antimicrobial peptides and neurotoxins, MgCRP-I peptides are produced as pre-propeptides, usually have a net positive charge and likely derive from similar evolutionary mechanisms, that is, gene duplication and positive selection within the mature peptide region; however, synthetic MgCRP-I peptides did not display significant toxicity in cultured mammalian cells, insecticidal, antimicrobial, or antifungal activities. The functional role of MgCRP-I peptides in mussel physiology still remains puzzling.

Description: Most sequenced eukaryotic genomes show a large excess of recent duplicates. As duplicates age, both the population genetic process of failed fixation and the mutation-driven process of nonfunctionalization act to reduce the observed number of duplicates. Understanding the processes generating the age distributions of recent duplicates is important to also understand the role of duplicate genes in the functional divergence of genomes. To date, mechanistic models for duplicate gene retention only account for the mutation-driven nonfunctionalization process. Here, a neutral model for the distribution of synonymous substitutions in duplicated genes which are segregating and expected to never fix in a population is introduced. This model enables differentiation of neutral loss due to failed fixation from loss due to mutation-driven nonfunctionalization. The model has been validated on simulated data and subsequent analysis with the model on genomic data from human and mouse shows that conclusions about the underlying mechanisms for duplicate gene retention can be sensitive to consideration of population genetic processes.

Description: Stenotrophomonas maltophilia , a ubiquitous Gram-negative -proteobacterium, has emerged as an important opportunistic pathogen responsible for nosocomial infections. A major characteristic of clinical isolates is their high intrinsic or acquired antibiotic resistance level. The aim of this study was to decipher the genetic determinism of antibiotic resistance among strains from different origins (i.e., natural environment and clinical origin) showing various antibiotic resistance profiles. To this purpose, we selected three strains isolated from soil collected in France or Burkina Faso that showed contrasting antibiotic resistance profiles. After whole-genome sequencing, the phylogenetic relationships of these 3 strains and 11 strains with available genome sequences were determined. Results showed that a strain’s phylogeny did not match their origin or antibiotic resistance profiles. Numerous antibiotic resistance coding genes and efflux pump operons were revealed by the genome analysis, with 57% of the identified genes not previously described. No major variation in the antibiotic resistance gene content was observed between strains irrespective of their origin and antibiotic resistance profiles. Although environmental strains generally carry as many multidrug resistant (MDR) efflux pumps as clinical strains, the absence of resistance–nodulation–division (RND) pumps (i.e., SmeABC) previously described to be specific to S. maltophilia was revealed in two environmental strains (BurA1 and PierC1). Furthermore the genome analysis of the environmental MDR strain BurA1 showed the absence of SmeABC but the presence of another putative MDR RND efflux pump, named EbyCAB on a genomic island probably acquired through horizontal gene transfer.

Description: The expansion of Bantu-speaking agropastoralist populations had a great impact on the genetic, linguistic, and cultural variation of sub-Saharan Africa. It is generally accepted that Bantu languages originated in an area around the present border between Cameroon and Nigeria approximately 5,000 years ago, from where they spread South and East becoming the largest African linguistic branch. The demic consequences of this event are reflected in the relatively high genetic homogeneity observed across most of sub-Saharan Africa populations. In this work, we explored genome-wide single nucleotide polymorphism data from 28 populations to characterize the genetic components present in sub-Saharan African populations. Combining novel data from four Southern African populations with previously published results, we reject the hypothesis that the "non-Bantu" genetic component reported in South-Eastern Africa (Mozambique) reflects extensive gene flow between incoming agriculturalist and resident hunter-gatherer communities. We alternatively suggest that this novel component is the result of demographic dynamics associated with the Bantu dispersal.

Description: Dietary shifts can drive molecular evolution in mammals and a major transition in human history, the agricultural revolution, favored carbohydrate consumption. We investigated the evolutionary history of nine genes encoding brush-border proteins involved in carbohydrate digestion/absorption. Results indicated widespread adaptive evolution in mammals, with several branches experiencing episodic selection, particularly strong in bats. Many positively selected sites map to functional protein regions (e.g., within glucosidase catalytic crevices), with parallel evolution at SI (sucrase-isomaltase) and MGAM (maltase-glucoamylase). In human populations, five genes were targeted by positive selection acting on noncoding variants within regulatory elements. Analysis of ancient DNA samples indicated that most derived alleles were already present in the Paleolithic. Positively selected variants at SLC2A5 (fructose transporter) were an exception and possibly spread following the domestication of specific fruit crops. We conclude that agriculture determined no major selective event at carbohydrate metabolism genes in humans, with implications for susceptibility to metabolic disorders.

Description: Ferns are one of the few remaining major clades of land plants for which a complete genome sequence is lacking. Knowledge of genome space in ferns will enable broad-scale comparative analyses of land plant genes and genomes, provide insights into genome evolution across green plants, and shed light on genetic and genomic features that characterize ferns, such as their high chromosome numbers and large genome sizes. As part of an initial exploration into fern genome space, we used a whole genome shotgun sequencing approach to obtain low-density coverage (~0.4X to 2X) for six fern species from the Polypodiales ( Ceratopteris , Pteridium , Polypodium , Cystopteris ), Cyatheales ( Plagiogyria ), and Gleicheniales ( Dipteris ). We explore these data to characterize the proportion of the nuclear genome represented by repetitive sequences (including DNA transposons, retrotransposons, ribosomal DNA, and simple repeats) and protein-coding genes, and to extract chloroplast and mitochondrial genome sequences. Such initial sweeps of fern genomes can provide information useful for selecting a promising candidate fern species for whole genome sequencing. We also describe variation of genomic traits across our sample and highlight some differences and similarities in repeat structure between ferns and seed plants.

Description: The evolution of heterogametic sex chromosomes is often—but not always—accompanied by the evolution of dosage compensating mechanisms that mitigate the impact of sex-specific gene dosage on levels of gene expression. One emerging view of this process is that such mechanisms may only evolve in male-heterogametic (XY) species but not in female-heterogametic (ZW) species, which will consequently exhibit "incomplete" sex chromosome dosage compensation. However, recent results suggest that at least some Lepidoptera (moths and butterflies) may prove to be an exception to this prediction. Studies in bombycoid moths indicate the presence of a chromosome-wide epigenetic mechanism that effectively balances Z chromosome gene expression between the sexes by reducing Z-linked expression in males. In contrast, strong sex chromosome dosage effects without any reduction in male Z-linked expression were previously reported in a pyralid moth, suggesting a lack of any such dosage compensating mechanism. Here we report an analysis of sex chromosome dosage compensation in Heliconius butterflies, sampling multiple individuals for several different adult tissues (head, abdomen, leg, mouth, and antennae). Methodologically, we introduce a novel application of linear mixed-effects models to assess dosage compensation, offering a unified statistical framework that can estimate effects specific to chromosome, to sex, and their interactions (i.e., a dosage effect). Our results show substantially reduced Z-linked expression relative to autosomes in both sexes, as previously observed in bombycoid moths. This observation is consistent with an increasing body of evidence that some lepidopteran species possess an epigenetic dosage compensating mechanism that reduces Z chromosome expression in males to levels comparable with females. However, this mechanism appears to be imperfect in Heliconius , resulting in a modest dosage effect that produces an average 5–20% increase in male expression relative to females on the Z chromosome, depending on the tissue. Thus our results in Heliconius reflect a mixture of previous patterns reported for Lepidoptera. In Heliconius, a moderate pattern of incomplete dosage compensation persists apparently despite the presence of an epigenetic dosage compensating mechanism. The chromosomal distributions of sex-biased genes show an excess of male-biased and a dearth of female-biased genes on the Z chromosome relative to autosomes, consistent with predictions of sexually antagonistic evolution.

Description: Mycobacterium avium ( M. a. ) subsp. paratuberculosis (MAP)—the etiologic agent of Johne’s disease—affects cattle, sheep, and other ruminants worldwide. To decipher phenotypic differences among sheep and cattle strains (belonging to MAP-S [Type-I/III], respectively, MAP-C [Type-II]), comparative genome analysis needs data from diverse isolates originating from different geographic regions of the world. This study presents the so far best assembled genome of a MAP-S-strain: Sheep isolate JIII-386 from Germany. One newly sequenced cattle isolate (JII-1961, Germany), four published MAP strains of MAP-C and MAP-S from the United States and Australia, and M. a. subsp. hominissuis (MAH) strain 104 were used for assembly improvement and comparisons. All genomes were annotated by BacProt and results compared with NCBI (National Center for Biotechnology Information) annotation. Corresponding protein-coding sequences (CDSs) were detected, but also CDSs that were exclusively determined by either NCBI or BacProt. A new Shine–Dalgarno sequence motif (5'-AGCTGG-3') was extracted. Novel CDSs including PE-PGRS family protein genes and about 80 noncoding RNAs exhibiting high sequence conservation are presented. Previously found genetic differences between MAP-types are partially revised. Four of ten assumed MAP-S-specific large sequence polymorphism regions (LSP S s) are still present in MAP-C strains; new LSP S s were identified. Independently of the regional origin of the strains, the number of individual CDSs and single nucleotide variants confirms the strong similarity of MAP-C strains and shows higher diversity among MAP-S strains. This study gives ambiguous results regarding the hypothesis that MAP-S is the evolutionary intermediate between MAH and MAP-C, but it clearly shows a higher similarity of MAP to MAH than to Mycobacterium intracellulare .

Description: Eukaryotic organelles depend on nuclear genes to perpetuate their biochemical integrity. This is true for mitochondria in all eukaryotes and plastids in plants and algae. Then how do kleptoplasts, plastids that are sequestered by some sacoglossan sea slugs, survive in the animals’ digestive gland cells in the absence of the algal nucleus encoding the vast majority of organellar proteins? For almost two decades, lateral gene transfer (LGT) from algae to slugs appeared to offer a solution, but RNA-seq analysis, later supported by genome sequencing of slug DNA, failed to find any evidence for such LGT events. Yet, isolated reports continue to be published and are readily discussed by the popular press and social media, making the data on LGT and its support for kleptoplast longevity appear controversial. However, when we take a sober look at the methods used, we realize that caution is warranted in how the results are interpreted. There is no evidence that the evolution of kleptoplasty in sea slugs involves LGT events. Based on what we know about photosystem maintenance in embryophyte plastids, we assume kleptoplasts depend on nuclear genes. However, studies have shown that some isolated algal plastids are, by nature, more robust than those of land plants. The evolution of kleptoplasty in green sea slugs involves many promising and unexplored phenomena, but there is no evidence that any of these require the expression of slug genes of algal origin.

Description: A number of sap-sucking insects harbor endosymbionts, which are thought to play an important role in the development of their hosts. One of the most important rice pests, the brown planthopper (BPH), Nilaparvata lugens (Stål), harbors an obligatory yeast-like symbiont (YLS) that cannot be cultured in vitro. Genomic information on this YLS would be useful to better understand its evolution. In this study, we performed genome sequencing of the YLS using both 454 and Illumina approaches, generating a draft genome that shows a slightly smaller genome size and relatively higher GC content than most ascomycete fungi. A phylogenomic analysis of the YLS supported its close relationship with insect pathogens. We analyzed YLS-specific genes and the categories of genes that are likely to have changed in the YLS during its evolution. The loss of mating type locus demonstrated in the YLS sheds light on the evolution of eukaryotic symbionts. This information about the YLS genome provides a helpful guide for further understanding endosymbiotic associations in hemiptera and the symbiotic replacement of ancient bacteria with a multifunctional YLS seems to have been a successful change.

Description: The RNA-directed DNA methylation (RdDM) pathway can be divided into three phases: 1) small interfering RNA biogenesis, 2) de novo methylation, and 3) chromatin modification. To determine the degree of conservation of this pathway we searched for key genes among land plants. We used OrthoMCL and the OrthoMCL Viridiplantae database to analyze proteomes of species in bryophytes, lycophytes, monilophytes, gymnosperms, and angiosperms. We also analyzed small RNA size categories and, in two gymnosperms, cytosine methylation in ribosomal DNA. Six proteins were restricted to angiosperms, these being NRPD4/NRPE4, RDM1, DMS3 (defective in meristem silencing 3), SHH1 (SAWADEE homeodomain homolog 1), KTF1, and SUVR2, although we failed to find the latter three proteins in Fritillaria persica , a species with a giant genome. Small RNAs of 24 nt in length were abundant only in angiosperms. Phylogenetic analyses of Dicer-like (DCL) proteins showed that DCL2 was restricted to seed plants, although it was absent in Gnetum gnemon and Welwitschia mirabilis . The data suggest that phases (1) and (2) of the RdDM pathway, described for model angiosperms, evolved with angiosperms. The absence of some features of RdDM in F. persica may be associated with its large genome. Phase (3) is probably the most conserved part of the pathway across land plants. DCL2, involved in virus defense and interaction with the canonical RdDM pathway to facilitate methylation of CHH, is absent outside seed plants. Its absence in G. gnemon , and W. mirabilis coupled with distinctive patterns of CHH methylation, suggest a secondary loss of DCL2 following the divergence of Gnetales.

Description: Comparisons of draft genome sequences of three geographically distinct isolates of Fusarium fujikuroi with two recently published genome sequences from the same species suggest diverse profiles of secondary metabolite production within F. fujikuroi . Species- and lineage-specific genes, many of which appear to exhibit expression profiles that are consistent with roles in host–pathogen interactions and adaptation to environmental changes, are concentrated in subtelomeric regions. These genomic compartments also exhibit distinct gene densities and compositional characteristics with respect to other genomic partitions, and likely play a role in the generation of molecular diversity. Our data provide additional evidence that gene duplication, divergence, and differential loss play important roles in F. fujikuroi genome evolution and suggest that hundreds of lineage-specific genes might have been acquired through horizontal gene transfer.

Description: Passerines are the largest avian order, and the 6,000 species comprise more than half of all extant bird species. This successful radiation probably had its origin in the Australasian region, but dating this origin has been difficult due to a scarce fossil record and poor biogeographic assumptions. Many of New Zealand’s endemic passerines fall within the deeper branches of the passerine radiation, and a well resolved phylogeny for the modern New Zealand element in the deeper branches of the oscine lineage will help us understand both oscine and passerine biogeography. To this end we present complete mitochondrial genomes representing all families of New Zealand passerines in a phylogenetic framework of over 100 passerine species. Dating analyses of this robust phylogeny suggest Passeriformes originated in the early Paleocene, with the major lineages of oscines "escaping" from Australasia about 30 Ma, and radiating throughout the world during the Oligocene. This independently derived conclusion is consistent with the passerine fossil record.

Description: Triplophysa fishes are the primary component of the fish fauna on the Tibetan Plateau and are well adapted to the high-altitude environment. Despite the importance of Triplophysa fishes on the plateau, the genetic mechanisms of the adaptations of these fishes to this high-altitude environment remain poorly understood. In this study, we generated the transcriptome sequences for three Triplophysa fishes, that is, Triplophysa siluroides , Triplophysa scleroptera , and Triplophysa dalaica , and used these and the previously available transcriptome and genome sequences from fishes living at low altitudes to identify potential genetic mechanisms for the high-altitude adaptations in Triplophysa fishes. An analysis of 2,269 orthologous genes among cave fish ( Astyanax mexicanus ), zebrafish ( Danio rerio ), large-scale loach ( Paramisgurnus dabryanus ), and Triplophysa fishes revealed that each of the terminal branches of the Triplophysa fishes had a significantly higher ratio of nonsynonymous to synonymous substitutions than that of the branches of the fishes from low altitudes, which provided consistent evidence for genome-wide rapid evolution in the Triplophysa genus. Many of the GO (Gene Ontology) categories associated with energy metabolism and hypoxia response exhibited accelerated evolution in the Triplophysa fishes compared with the large-scale loach. The genes that exhibited signs of positive selection and rapid evolution in the Triplophysa fishes were also significantly enriched in energy metabolism and hypoxia response categories. Our analysis identified widespread Triplophysa -specific nonsynonymous mutations in the fast evolving genes and positively selected genes. Moreover, we detected significant evidence of positive selection in the HIF (hypoxia-inducible factor)-1A and HIF-2B genes in Triplophysa fishes and found that the Triplophysa -specific nonsynonymous mutations in the HIF-1A and HIF-2B genes were associated with functional changes. Overall, our study provides new insights into the adaptations and evolution of fishes in the high-altitude environment of the Tibetan Plateau and complements previous findings on the adaptations of mammals and birds to high altitudes.

Description: Sponges harbor a complex consortium of microbial communities living in symbiotic relationship benefiting each other through the integration of metabolites. The mechanisms influencing a successful microbial association with a sponge partner are yet to be fully understood. Here, we sequenced the genome of Pseudovibrio sp. POLY-S9 strain isolated from the intertidal marine sponge Polymastia penicillus sampled from the Atlantic coast of Portugal to identify the genomic features favoring the symbiotic relationship. The draft genome revealed an exceptionally large genome size of 6.6 Mbp compared with the previously reported genomes of the genus Pseudovibrio isolated from a coral and a sponge larva. Our genomic study detected the presence of several biosynthetic gene clusters—polyketide synthase, nonribosomal peptide synthetase and siderophore—affirming the potential ability of the genus Pseudovibrio to produce a wide variety of metabolic compounds. Moreover, we identified a repertoire of genes encoding adaptive symbioses factors (eukaryotic-like proteins), such as the ankyrin repeats, tetratrico peptide repeats, and Sel1 repeats that improve the attachment to the eukaryotic hosts and the avoidance of the host’s immune response . The genome also harbored a large number of mobile elements (~5%) and gene transfer agents, which explains the massive genome expansion and suggests a possible mechanism of horizontal gene transfer. In conclusion, the genome of POLY-S9 exhibited an increase in size, number of mobile DNA, multiple metabolite gene clusters, and secretion systems, likely to influence the genome diversification and the evolvability.

Description: Gene regulatory networks (GRN) are central to developmental processes. They are composed of transcription factors and signaling molecules orchestrating gene expression modules that tightly regulate the development of organisms. The neural crest (NC) is a multipotent cell population that is considered a key innovation of vertebrates. Its derivatives contribute to shaping the astounding morphological diversity of jaws, teeth, head skeleton, or pigmentation. Here, we study the molecular evolution of the NC GRN by analyzing patterns of molecular divergence for a total of 36 genes in 16 species of bony fishes. Analyses of nonsynonymous to synonymous substitution rate ratios (d N /d S ) support patterns of variable selective pressures among genes deployed at different stages of NC development, consistent with the developmental hourglass model. Model-based clustering techniques of sequence features support the notion of extreme conservation of NC-genes across the entire network. Our data show that most genes are under strong purifying selection that is maintained throughout ray-finned fish evolution. Late NC development genes reveal a pattern of increased constraints in more recent lineages. Additionally, seven of the NC-genes showed signs of relaxation of purifying selection in the famously species-rich lineage of cichlid fishes. This suggests that NC genes might have played a role in the adaptive radiation of cichlids by granting flexibility in the development of NC-derived traits—suggesting an important role for NC network architecture during the diversification in vertebrates.

Description: Cymbomonas tetramitiformis —a marine prasinophyte—is one of only a few green algae that still retain an ancestral particulate-feeding mechanism while harvesting energy through photosynthesis. The genome of the alga is estimated to be 850 Mb–1.2 Gb in size—the bulk of which is filled with repetitive sequences—and is annotated with 37,366 protein-coding gene models. A number of unusual metabolic pathways (for the Chloroplastida) are predicted for C. tetramitiformis , including pathways for Lipid-A and peptidoglycan metabolism. Comparative analyses of the predicted peptides of C. tetramitiformis to sets of other eukaryotes revealed that nonphagocytes are depleted in a number of genes, a proportion of which have known function in feeding. In addition, our analysis suggests that obligatory phagotrophy is associated with the loss of genes that function in biosynthesis of small molecules (e.g., amino acids). Further, C. tetramitiformis and at least one other phago-mixotrophic alga are thus unique, compared with obligatory heterotrophs and nonphagocytes, in that both feeding and small molecule synthesis-related genes are retained in their genomes. These results suggest that early, ancestral host eukaryotes that gave rise to phototrophs had the capacity to assimilate building block molecules from inorganic substances (i.e., prototrophy). The loss of biosynthesis genes, thus, may at least partially explain the apparent lack of instances of permanent incorporation of photosynthetic endosymbionts in later-divergent, auxotrophic eukaryotic lineages, such as metazoans and ciliates.

Description: Chlamydiaceae are obligate intracellular bacteria that cause a diversity of severe infections among humans and livestock on a global scale. Identification of new species since 1989 and emergence of zoonotic infections, including abortion in women, underscore the need for genome sequencing of multiple strains of each species to advance our knowledge of evolutionary dynamics across Chlamydiaceae . Here, we genome sequenced isolates from avian, lower mammalian and human hosts. Based on core gene phylogeny, five isolates previously classified as Chlamydia abortus were identified as members of Chlamydia psittaci and Chlamydia pecorum. Chlamydia abortus is the most recently emerged species and is a highly monomorphic group that lacks the conserved virulence-associated plasmid. Low-level recombination and evidence for adaptation to the placenta echo evolutionary processes seen in recently emerged, highly virulent niche-restricted pathogens, such as Bacillus anthracis . In contrast, gene flow occurred within C. psittaci and other Chlamydiaceae species. The C. psittaci strain RTH, isolated from a red-tailed hawk ( Buteo jamaicensis ), is an outlying strain with admixture of C. abortus , C. psittaci , and its own population markers. An average nucleotide identity of less than 94% compared with other Chlamydiaceae species suggests that RTH belongs to a new species intermediary between C. psittaci and C. abortus . Hawks, as scavengers and predators, have extensive opportunities to acquire multiple species in their intestinal tract. This could facilitate transformation and homologous recombination with the potential for new species emergence. Our findings indicate that incubator hosts such as birds-of-prey likely promote Chlamydiaceae evolution resulting in novel pathogenic lineages.

Description: The causes of the great variation in nucleotide composition of prokaryotic genomes have long been disputed. Here, we use extensive metagenomic and whole-genome data to demonstrate that both phylogeny and the environment shape prokaryotic nucleotide content. We show that across environments, various phyla are characterized by different mean guanine and cytosine (GC) values as well as by the extent of variation on that mean value. At the same time, we show that GC-content varies greatly as a function of environment, in a manner that cannot be entirely explained by disparities in phylogenetic composition. We find environmentally driven differences in nucleotide content not only between highly diverged environments (e.g., soil, vs. aquatic vs. human gut) but also within a single type of environment. More specifically, we demonstrate that some human guts are associated with a microbiome that is consistently more GC-rich across phyla, whereas others are associated with a more AT-rich microbiome. These differences appear to be driven both by variations in phylogenetic composition and by environmental differences—which are independent of these phylogenetic composition differences. Combined, our results demonstrate that both phylogeny and the environment significantly affect nucleotide composition and that the environmental differences affecting nucleotide composition are far subtler than previously appreciated.

Description: Brachiopods are a lineage of invertebrates well known for the breadth and depth of their fossil record. Although the quality of this fossil record attracts the attention of paleontologists, geochemists, and paleoclimatologists, modern day brachiopods are also of interest to evolutionary biologists due to their potential to address a variety of questions ranging from developmental biology to biomineralization. The brachiopod shell is a composite material primarily composed of either calcite or calcium phosphate in close association with proteins and polysaccharides which give these composite structures their material properties. The information content of these biomolecules, sequestered within the shell during its construction, has the potential to inform hypotheses focused on describing how brachiopod shell formation evolved. Here, using high throughput proteomic approaches and next generation sequencing, we have surveyed and characterized the first shell-proteome and shell-forming transcriptome of any brachiopod, the South American Magellania venosa (Rhynchonelliformea: Terebratulida) . We find that the seven most abundant proteins present in the shell are unique to M. venosa , but that these proteins display biochemical features found in other metazoan biomineralization proteins. We can also detect some M. venosa proteins that display significant sequence similarity to other metazoan biomineralization proteins, suggesting that some elements of the brachiopod shell-forming proteome are deeply evolutionarily conserved. We also employed a variety of preparation methods to isolate shell proteins and find that in comparison to the shells of other spiralian invertebrates (such as mollusks) the shell ultrastructure of M. venosa may explain the effects these preparation strategies have on our results.

Description: Free fatty acid receptors (FFAR) belong to a family of five G-protein coupled receptors that are involved in the regulation of lipid metabolism, so that their loss of function increases the risk of obesity. The aim of this study was to determine the expansion of genes encoding paralogs of FFAR2 in the chicken, considered as a model organism for developmental biology and biomedical research. By estimating the gene copy number using quantitative polymerase chain reaction, genomic DNA resequencing, and RNA sequencing data, we showed the existence of 23 ± 1.5 genes encoding FFAR2 paralogs in the chicken genome. The FFAR2 paralogs shared an identity from 87.2% up to 99%. Extensive gene conversion was responsible for this high degree of sequence similarities between these genes, and this concerned especially the four amino acids known to be critical for ligand binding. Moreover, elevated nonsynonymous/synonymous substitution ratios on some amino acids within or in close-vicinity of the ligand-binding groove suggest that positive selection may have reduced the effective rate of gene conversion in this region, thus contributing to diversify the function of some FFAR2 paralogs. All the FFAR2 paralogs were located on a microchromosome in a same linkage group. FFAR2 genes were expressed in different tissues and cells such as spleen, peripheral blood mononuclear cells, abdominal adipose tissue, intestine, and lung, with the highest rate of expression in testis. Further investigations are needed to determine whether these chicken-specific events along evolution are the consequence of domestication and may play a role in regulating lipid metabolism in this species.

Description: Loss of introns in plant mitochondrial genes is commonly explained by retroprocessing. Under this model, an mRNA is reverse transcribed and integrated back into the genome, simultaneously affecting the contents of introns and edited sites. To evaluate the extent to which retroprocessing explains intron loss, we analyzed patterns of intron content and predicted RNA editing for whole mitochondrial genomes of 30 species in the monocot order Alismatales. In this group, we found an unusually high degree of variation in the intron content, even expanding the hitherto known variation among angiosperms. Some species have lost some two-third of the cis -spliced introns. We found a strong correlation between intron content and editing frequency, and detected 27 events in which intron loss is consistent with the presence of nucleotides in an edited state, supporting retroprocessing. However, we also detected seven cases of intron loss not readily being explained by retroprocession. Our analyses are also not consistent with the entire length of a fully processed cDNA copy being integrated into the genome, but instead indicate that retroprocessing usually occurs for only part of the gene. In some cases, several rounds of retroprocessing may explain intron loss in genes completely devoid of introns. A number of taxa retroprocessing seem to be very common and a possibly ongoing process. It affects the entire mitochondrial genome.

Description: The luminous bacterial symbionts of anomalopid flashlight fish are thought to be obligately dependent on their hosts for growth and share several aspects of genome evolution with unrelated obligate symbionts, including genome reduction. However, in contrast to most obligate bacteria, anomalopid symbionts have an active environmental phase that may be important for symbiont transmission. Here we investigated patterns of evolution between anomalopid symbionts compared with patterns in free-living relatives and unrelated obligate symbionts to determine if trends common to obligate symbionts are also found in anomalopid symbionts. Two symbionts, " Candidatus Photodesmus katoptron" and " Candidatus Photodesmus blepharus," have genomes that are highly similar in gene content and order, suggesting genome stasis similar to ancient obligate symbionts present in insect lineages. This genome stasis exists in spite of the symbiont’s inferred ability to recombine, which is frequently lacking in obligate symbionts with stable genomes. Additionally, we used genome comparisons and tests of selection to infer which genes may be particularly important for the symbiont’s ecology compared with relatives. In keeping with obligate dependence, substitution patterns suggest that most symbiont genes are experiencing relaxed purifying selection compared with relatives. However, genes involved in motility and carbon storage, which are likely to be used outside the host, appear to be under increased purifying selection. Two chemoreceptor chemotaxis genes are retained by both species and show high conservation with amino acid sensing genes, suggesting that the bacteria may actively seek out hosts using chemotaxis toward amino acids, which the symbionts are not able to synthesize.

Description: Cynomoriaceae, one of the last unplaced families of flowering plants, comprise one or two species or subspecies of root parasites that occur from the Mediterranean to the Gobi Desert. Using Illumina sequencing, we assembled the mitochondrial and plastid genomes as well as some nuclear genes of a Cynomorium specimen from Italy. Selected genes were also obtained by Sanger sequencing from individuals collected in China and Iran, resulting in matrices of 33 mitochondrial, 6 nuclear, and 14 plastid genes and rDNAs enlarged to include a representative angiosperm taxon sampling based on data available in GenBank. We also compiled a new geographic map to discern possible discontinuities in the parasites’ occurrence. Cynomorium has large genomes of 13.70–13.61 (Italy) to 13.95–13.76 pg (China). Its mitochondrial genome consists of up to 49 circular subgenomes and has an overall gene content similar to that of photosynthetic angiosperms, while its plastome retains only 27 of the normally 116 genes. Nuclear, plastid and mitochondrial phylogenies place Cynomoriaceae in Saxifragales, and we found evidence for several horizontal gene transfers from different hosts, as well as intracellular gene transfers.

Description: The male-specific regions of primate Y-chromosomes (MSY) are enriched for multi-copy genes highly expressed in the testis. These genes are located in large repetitive sequences arranged as palindromes, inverted-, and tandem repeats termed amplicons. In humans, these genes have critical roles in male fertility and are essential for the production of sperm. The structure of human and chimpanzee amplicon sequences show remarkable difference relative to the remainder of the genome, a difference that may be the result of intense selective pressure on male fertility. Four subspecies of common chimpanzees have undergone extended periods of isolation and appear to be in the early process of subspeciation. A recent study found amplicons enriched for testis-expressed genes on the primate X-chromosome the target of hard selective sweeps, and male-fertility genes on the Y-chromosome may also be the targets of selection. However, little is understood about Y-chromosome amplicon diversity within and across chimpanzee populations. Here, we analyze nine common chimpanzee (representing three subspecies: Pan troglodytes schweinfurthii , Pan troglodytes ellioti , and Pan troglodytes verus ) and two bonobo ( Pan paniscus) male whole-genome sequences to assess Y ampliconic copy-number diversity across the Pan genus. We observe that the copy number of Y chromosome amplicons is variable among chimpanzees and bonobos, and identify several lineage-specific patterns, including variable copy number of azoospermia candidates RBMY and DAZ . We detect recurrent switchpoints of copy-number change along the ampliconic tracts across chimpanzee populations, which may be the result of localized genome instability or selective forces.

Description: De novo genes, which originate from ancestral nongenic sequences, are one of the most important sources of protein-coding genes. This origination process is crucial for the adaptation of organisms. However, how de novo genes arise and become fixed in a population or species remains largely unknown. Here, we identified 782 de novo genes from the model plant Arabidopsis thaliana and divided them into three types based on the availability of translational evidence, transcriptional evidence, and neither transcriptional nor translational evidence for their origin. Importantly, by integrating multiple types of omics data, including data from genomes, epigenomes, transcriptomes, and translatomes, we found that epigenetic modifications (DNA methylation and histone modification) play an important role in the origination process of de novo genes. Intriguingly, using the transcriptomes and methylomes from the same population of 84 accessions, we found that de novo genes that are transcribed in approximately half of the total accessions within the population are highly methylated, with lower levels of transcription than those transcribed at other frequencies within the population. We hypothesized that, during the origin of de novo gene alleles, those neutralized to low expression states via DNA methylation have relatively high probabilities of spreading and becoming fixed in a population. Our results highlight the process underlying the origin of de novo genes at the population level, as well as the importance of DNA methylation in this process.

Description: Halophilic microbes living in hypersaline environments must counteract the detrimental effects of low water activity and salt interference. Some halophilic prokaryotes equilibrate their intracellular osmotic strength with the extracellular milieu by importing inorganic solutes, mainly potassium. These "salt-in" organisms characteristically have proteins that are highly enriched with acidic and hydrophilic residues. In contrast, "salt-out" halophiles accumulate large amounts of organic solutes like amino acids, sugars and polyols, and lack a strong signature of halophilicity in the amino acid composition of cytoplasmic proteins. Studies to date have examined halophilic prokaryotes, yeasts, or algae, thus virtually nothing is known about the molecular adaptations of the other eukaryotic microbes, that is, heterotrophic protists (protozoa), that also thrive in hypersaline habitats. We conducted transcriptomic investigations to unravel the molecular adaptations of two obligately halophilic protists, Halocafeteria seosinensis and Pharyngomonas kirbyi . Their predicted cytoplasmic proteomes showed increased hydrophilicity compared with marine protists. Furthermore, analysis of reconstructed ancestral sequences suggested that, relative to mesophiles, proteins in halophilic protists have undergone fewer substitutions from hydrophilic to hydrophobic residues since divergence from their closest relatives. These results suggest that these halophilic protists have a higher intracellular salt content than marine protists. However, absence of the acidic signature of salt-in microbes suggests that Haloc. seosinensis and P. kirbyi utilize organic osmolytes to maintain osmotic equilibrium. We detected increased expression of enzymes involved in synthesis and transport of organic osmolytes, namely hydroxyectoine and myo-inositol, at maximal salt concentration for growth in Haloc. seosinensis , suggesting possible candidates for these inferred organic osmolytes.

Description: MicroRNAs are small (~22 nt) noncoding RNAs that repress translation and therefore regulate the production of proteins from specific target mRNAs. microRNAs have been found to function in diverse aspects of gene regulation within animal development and many other processes. Among invertebrates, both conserved and novel, lineage specific, microRNAs have been extensively studied predominantly in holometabolous insects such as Drosophila melanogaster . However little is known about microRNA repertoires in other arthropod lineages such as the chelicerates. To understand the evolution of microRNAs in this poorly sampled subphylum, we characterized the microRNA repertoire expressed during embryogenesis of the common house spider Parasteatoda tepidariorum . We identified a total of 148 microRNAs in P. tepidariorum representing 66 families. Approximately half of these microRNA families are conserved in other metazoans, while the remainder are specific to this spider. Of the 35 conserved microRNAs families 15 had at least two copies in the P. tepidariorum genome. A BLAST-based approach revealed a similar pattern of duplication in other spiders and a scorpion, but not among other chelicerates and arthropods, with the exception of a horseshoe crab. Among the duplicated microRNAs we found examples of lineage-specific tandem duplications, and the duplication of entire microRNA clusters in three spiders, a scorpion, and in a horseshoe crab. Furthermore, we found that paralogs of many P. tepidariorum microRNA families exhibit arm switching, which suggests that duplication was often followed by sub- or neofunctionalization. Our work shows that understanding the evolution of microRNAs in the chelicerates has great potential to provide insights into the process of microRNA duplication and divergence and the evolution of animal development.

Description: Family Hominidae, which includes humans and great apes, is recognized for unique complex social behavior and intellectual abilities. Despite the increasing genome data, however, the genomic origin of its phenotypic uniqueness has remained elusive. Clade-specific genes and highly conserved noncoding sequences (HCNSs) are among the high-potential evolutionary candidates involved in driving clade-specific characters and phenotypes. On this premise, we analyzed whole genome sequences along with gene orthology data retrieved from major DNA databases to find Hominidae-specific (HS) genes and HCNSs. We discovered that Down syndrome critical region 4 ( DSCR4 ) is the only experimentally verified gene uniquely present in Hominidae. DSCR4 has no structural homology to any known protein and was inferred to have emerged in several steps through LTR/ERV1, LTR/ERVL retrotransposition, and transversion. Using the genomic distance as neutral evolution threshold, we identified 1,658 HS HCNSs. Polymorphism coverage and derived allele frequency analysis of HS HCNSs showed that these HCNSs are under purifying selection, indicating that they may harbor important functions. They are overrepresented in promoters/untranslated regions, in close proximity of genes involved in sensory perception of sound and developmental process, and also showed a significantly lower nucleosome occupancy probability. Interestingly, many ancestral sequences of the HS HCNSs showed very high evolutionary rates. This suggests that new functions emerged through some kind of positive selection, and then purifying selection started to operate to keep these functions.

Description: To progress our understanding of molecular evolution from a collection of well-studied genes toward the level of the cell, we must consider whole systems. Here, we reveal the evolution of an important intracellular signaling system. The calcium-signaling toolkit is made up of different multidomain proteins that have undergone duplication, recombination, sequence divergence, and selection. The picture of evolution, considering the repertoire of proteins in the toolkit of both extant organisms and ancestors, is radically different from that of other systems. In eukaryotes, the repertoire increased in both abundance and diversity at a far greater rate than general genomic expansion. We describe how calcium-based intracellular signaling evolution differs not only in rate but in nature, and how this correlates with the disparity of plants and animals.

Description: Endosymbiosis is an important evolutionary event for organisms, and there is widespread interest in understanding the evolution of endosymbiosis establishment. Hydra is one of the most suitable organisms for studying the evolution of endosymbiosis. Within the genus Hydra , H. viridissima and H. vulgaris show endosymbiosis with green algae. Previous studies suggested that the endosymbiosis in H. vulgaris took place much more recently than that in H. viridissima , noting that the establishment of the interaction between H. vulgaris and its algae is not as stable as in H. viridissima. To investigate the on-going process of endosymbiosis, we first compared growth and tolerance to starvation in symbiotic and aposymbiotic polyps of both species. The results revealed that symbiotic H. viridissima had a higher growth rate and greater tolerance to starvation than aposymbiotic polyps. By contrast, growth of symbiotic H. vulgaris was identical to that of aposymbiotic polyps, and symbiotic H. vulgaris was less tolerant to starvation. Moreover, our gene expression analysis showed a pattern of differential gene expression in H. viridissima similar to that in other endosymbiotically established organisms, and contrary to that observed in H. vulgaris . We also showed that H. viridissima could cope with oxidative stress that caused damage, such as cell death, in H. vulgaris . These observations support the idea that oxidative stress related genes play an important role in the on-going process of endosymbiosis evolution. The different evolutionary stages of endosymbiosis studied here provide a deeper insight into the evolutionary processes occurring toward a stable endosymbiosis.

Description: The plastid genome (plastome) of heterotrophic plants like mycoheterotrophs and parasites shows massive gene losses in consequence to the relaxation of functional constraints on photosynthesis. To understand the patterns of this convergent plastome reduction syndrome in heterotrophic plants, we studied 12 closely related orchids of three different lifeforms from the tribe Neottieae (Orchidaceae). We employ a comparative genomics approach to examine structural and selectional changes in plastomes within Neottieae. Both leafy and leafless heterotrophic species have functionally reduced plastid genome. Our analyses show that genes for the NAD(P)H dehydrogenase complex, the photosystems, and the RNA polymerase have been lost functionally multiple times independently. The physical reduction proceeds in a highly lineage-specific manner, accompanied by structural reconfigurations such as inversions or modifications of the large inverted repeats. Despite significant but minor selectional changes, all retained genes continue to evolve under purifying selection. All leafless Neottia species, including both visibly green and nongreen members, are fully mycoheterotrophic, likely evolved from leafy and partially mycoheterotrophic species. The plastomes of Neottieae span many stages of plastome degradation, including the longest plastome of a mycoheterotroph, providing invaluable insights into the mechanisms of plastome evolution along the transition from autotrophy to full mycoheterotrophy.

Description: In the experimental evolution of microbes such as Escherichia coli , many replicate populations are evolved from a common ancestor. Freezing a population sample supplemented with the cryoprotectant glycerol permits later analysis or restarting of an evolution experiment. Typically, each evolving population, and thus each sample archived in this way, consists of many unique genotypes and phenotypes. The effect of archiving on such a heterogeneous population is unknown. Here, we identified optimal archiving conditions for E. coli. We also used genome sequencing of archived samples to study the effects that archiving has on genomic population diversity. We observed no allele substitutions and mostly small changes in allele frequency. Nevertheless, principal component analysis of genome-scale allelic diversity shows that archiving affects diversity across many loci. We showed that this change in diversity is due to selection rather than drift. In addition, ~1% of rare alleles that occurred at low frequencies were lost after treatment. Our observations imply that archived populations may be used to conduct fitness or other phenotypic assays of populations, in which the loss of a rare allele may have negligible effects. However, caution is appropriate when sequencing populations restarted from glycerol stocks, as well as when using glycerol stocks to restart or replay evolution. This is because the loss of rare alleles can alter the future evolutionary trajectory of a population if the lost alleles were strongly beneficial.

Description: PIWI proteins and PIWI-interacting RNAs (piRNAs) are part of a cellular pathway that has evolved to protect genomes against the proliferation of transposable elements (TEs). PIWIs and piRNAs assemble into complexes that are involved in epigenetic and post-transcriptional repression of TEs. Most of our understanding of the mechanisms of piRNA-mediated TE silencing comes from fruit fly and mouse models. However, even in these well-studied animals it is unclear how piRNA responses relate to variable TE expression and whether the strength of the piRNA response affects TE content over time. Here, we assessed the evolutionary interactions between TE and piRNAs in a statistical framework using three nonmodel laurasiatherian mammals as a study system: dog, horse, and a vesper bat. These three species diverged ~80 million years ago and have distinct genomic TE contents. By comparing species with distinct TE landscapes, we aimed to identify clear relationships among TE content, expression, and piRNAs. We found that the TE subfamilies that are the most transcribed appear to elicit the strongest "ping-pong" response. This was most evident among long interspersed elements, but the relationships between expression and ping-pong pilRNA (piRNA-like) expression were more complex among SINEs. SINE transcripts were equally abundant in the dog and horse yet new SINE insertions were relatively rare in the horse genome, where we identified a stronger piRNA response. Our analyses suggest that the piRNA response can have a strong impact on the TE composition of a genome. However, our results also suggest that the presence of a robust piRNA response is apparently not sufficient to stop TE mobilization and accumulation.

Description: We introduce a gene tree simulator that is designed for use in conjunction with approximate Bayesian computation approaches. We show that it can be used to determine the relative importance of hybrid speciation and introgression compared with incomplete lineage sorting (ILS) in producing patterns of incongruence across gene trees. Important features of the new simulator are (1) a choice of models to capture the decreasing probability of successful hybrid species formation or introgression as a function of genetic distance between potential parent species; (2) the ability for hybrid speciation to result in asymmetrical contributions of genetic material from each parent species; (3) the ability to vary the rates of hybrid speciation, introgression, and divergence speciation in different epochs; and (4) incorporation of the coalescent, so that patterns of incongruence due to ILS can be compared with those due to hybrid evolution. Given a set of gene trees generated by the simulator, we calculate a set of statistics, each measuring in a different way the discordance between the gene trees. We show that these statistics can be used to differentiate whether the gene tree discordance was largely due to hybridization, or only due to lineage sorting.

Description: Beetles represent the largest insect order and they display extreme morphological, ecological and behavioral diversity, which makes them ideal models for evolutionary studies. Here, we present the draft genome of the scarab beetle Oryctes borbonicus , which has a more basal phylogenetic position than the two previously sequenced pest species Tribolium castaneum and Dendroctonus ponderosae providing the potential for sequence polarization. Oryctes borbonicus is endemic to La Réunion, an island located in the Indian Ocean, and is the host of the nematode Pristionchus pacificus , a well-established model organism for integrative evolutionary biology. At 518 Mb, the O. borbonicus genome is substantially larger and encodes more genes than T. castaneum and D. ponderosae . We found that only 25% of the predicted genes of O. borbonicus are conserved as single copy genes across the nine investigated insect genomes, suggesting substantial gene turnover within insects. Even within beetles, up to 21% of genes are restricted to only one species, whereas most other genes have undergone lineage-specific duplications and losses. We illustrate lineage-specific duplications using detailed phylogenetic analysis of two gene families. This study serves as a reference point for insect/coleopteran genomics, although its original motivation was to find evidence for potential horizontal gene transfer (HGT) between O. borbonicus and P. pacificus . The latter was previously shown to be the recipient of multiple horizontally transferred genes including some genes from insect donors. However, our study failed to provide any clear evidence for additional HGTs between the two species.

Description: Collembola (springtails) are detritivorous hexapods that inhabit the soil and its litter layer. The ecology of the springtail Orchesella cincta is extensively studied in the context of adaptation to anthropogenically disturbed areas. Here, we present a draft genome of an O. cincta reference strain with an estimated size of 286.8 Mbp, containing 20,249 genes. In total, 446 gene families are expanded and 1,169 gene families evolved specific to this lineage. Besides these gene families involved in general biological processes, we observe gene clusters participating in xenobiotic biotransformation. Furthermore, we identified 253 cases of horizontal gene transfer (HGT). Although the largest percentage of them originated from bacteria (37.5%), we observe an unusually high percentage (30.4%) of such genes of fungal origin. The majority of foreign genes are involved in carbohydrate metabolism and cellulose degradation. Moreover, some foreign genes (e.g., bacillopeptidases) expanded after HGT. We hypothesize that horizontally transferred genes could be advantageous for food processing in a soil environment that is full of decaying organic material. Finally, we identified several lineage-specific genes, expanded gene families, and horizontally transferred genes, associated with altered gene expression as a consequence of genetic adaptation to metal stress. This suggests that these genome features may be preadaptations allowing natural selection to act on. In conclusion, this genome study provides a solid foundation for further analysis of evolutionary mechanisms of adaptation to environmental stressors.

Description: The question on the patterns and limits of reduction of plastid genomes in nonphotosynthetic plants and the reasons of their conservation is one of the intriguing topics in plant genome evolution. Here, we report sequencing and analysis of plastid genome in nonphotosynthetic orchids Epipogium aphyllum and Epipogium roseum , which, with sizes of 31 and 19 kbp, respectively, represent the smallest plastid genomes characterized by now. Besides drastic reduction, which is expected, we found several unusual features of these "minimal" plastomes: Multiple rearrangements, highly biased nucleotide composition, and unprecedentedly high substitution rate. Only 27 and 29 genes remained intact in the plastomes of E. aphyllum and E. roseum —those encoding ribosomal components, transfer RNAs, and three additional housekeeping genes ( infA , clpP , and accD ). We found no signs of relaxed selection acting on these genes. We hypothesize that the main reason for retention of plastid genomes in Epipogium is the necessity to translate messenger RNAs (mRNAs) of accD and/or clpP proteins which are essential for cell metabolism. However, these genes are absent in plastomes of several plant species; their absence is compensated by the presence of a functional copy arisen by gene transfer from plastid to the nuclear genome. This suggests that there is no single set of plastid-encoded essential genes, but rather different sets for different species and that the retention of a gene in the plastome depends on the interaction between the nucleus and plastids.

Description: Repetitive DNA, including transposable elements (TEs), is found throughout eukaryotic genomes. Annotating and assembling the "repeatome" during genome-wide analysis often poses a challenge. To address this problem, we present dnaPipeTE—a new bioinformatics pipeline that uses a sample of raw genomic reads. It produces precise estimates of repeated DNA content and TE consensus sequences, as well as the relative ages of TE families. We shows that dnaPipeTE performs well using very low coverage sequencing in different genomes, losing accuracy only with old TE families. We applied this pipeline to the genome of the Asian tiger mosquito Aedes albopictu s, an invasive species of human health interest, for which the genome size is estimated to be over 1 Gbp. Using dnaPipeTE, we showed that this species harbors a large (50% of the genome) and potentially active repeatome with an overall TE class and order composition similar to that of Aedes aegypti , the yellow fever mosquito. However, intraorder dynamics show clear distinctions between the two species, with differences at the TE family level. Our pipeline’s ability to manage the repeatome annotation problem will make it helpful for new or ongoing assembly projects, and our results will benefit future genomic studies of A. albopictus.

Description: River buffalo, Bubalus bubalis is a large bovine species frequently used livestock in southern Asia. It is believed that the river buffalo was domesticated from Bubalus arnee , the wild buffalo of mainland Asia, a few thousand years ago, probably during the period of Indus Valley civilization. However, the domestication history of the river buffalo has been the subject of debate for many decades mainly due to the lack of clear archeological evidence and the divisive conclusions of the genetic studies. Therefore, in order to understand the domestication history and genetic relationship among the various river buffalo populations, we analyzed 492-bp region of mitochondrial DNA control region sequences of 414 river buffalo sampled from India, Pakistan, Egypt, and Iran along with the available 403 swamp buffalo sequences. The phylogenetic analyses of our study along with the archaeological evidence suggest that the river buffalo was domesticated in an atypical manner involving continuous introgression of wild animals to the domestic stocks in Indian subcontinent prior to mature phase of Indus Valley civilization (2600–1900 BC). Specifically, our data exclude Mesopotamian region as the place of domestication of the river buffalo.

Description: Although adaptive mutations are often considered to be dominant, it has been recently shown that a substantial proportion of adaptive mutations should display heterozygote advantage. In this work, we take advantage of a recently characterized transposable element insertion mediating oxidative stress response in Drosophila melanogaster to test the dominance effect of an adaptive mutation. The comparison of the survival curves of heterozygous and the two corresponding homozygous flies indicated that the dominance effect of Bari-Jheh depends on the genetic background. Both in homozygous and in heterozygous flies, Bari-Jheh was associated with upregulation of Jheh1 ( Juvenile Hormone Epoxyde Hydrolase 1 ) and/or Jheh2 genes. Our results add to the limited number of studies in which the dominance effect of adaptive mutations has been empirically estimated and highlights the complexity of their inheritance.

Description: The five macromolecular complexes that jointly mediate oxidative phosphorylation (OXPHOS) in mitochondria consist of many more subunits than those of bacteria, yet, it remains unclear by which evolutionary mechanism(s) these novel subunits were recruited. Even less well understood is the structural evolution of mitochondrial ribosomes (mitoribosomes): while it was long thought that their exceptionally high protein content would physically compensate for their uniquely low amount of ribosomal RNA (rRNA), this hypothesis has been refuted by structural studies. Here, we present a cryo-electron microscopy structure of the 73S mitoribosome from Neurospora crassa , together with genomic and proteomic analyses of mitoribosome composition across the eukaryotic domain. Surprisingly, our findings reveal that both structurally and compositionally, mitoribosomes have evolved very similarly to mitochondrial OXPHOS complexes via two distinct phases: A constructive phase that mainly acted early in eukaryote evolution, resulting in the recruitment of altogether approximately 75 novel subunits, and a reductive phase that acted during metazoan evolution, resulting in gradual length-reduction of mitochondrially encoded rRNAs and OXPHOS proteins. Both phases can be well explained by the accumulation of (slightly) deleterious mutations and deletions, respectively, in mitochondrially encoded rRNAs and OXPHOS proteins. We argue that the main role of the newly recruited (nuclear encoded) ribosomal- and OXPHOS proteins is to provide structural compensation to the mutationally destabilized mitochondrially encoded components. While the newly recruited proteins probably provide a selective advantage owing to their compensatory nature, and while their presence may have opened evolutionary pathways toward novel mitochondrion-specific functions, we emphasize that the initial events that resulted in their recruitment was nonadaptive in nature. Our framework is supported by population genetic studies, and it can explain the complete structural evolution of mitochondrial ribosomes and OXPHOS complexes, as well as many observed functions of individual proteins.

Description: Indigenous populations of Malaysia known as Orang Asli (OA) show huge morphological, anthropological, and linguistic diversity. However, the genetic history of these populations remained obscure. We performed a high-density array genotyping using over 2 million single nucleotide polymorphisms in three major groups of Negrito, Senoi, and Proto-Malay. Structural analyses indicated that although all OA groups are genetically closest to East Asian (EA) populations, they are substantially distinct. We identified a genetic affinity between Andamanese and Malaysian Negritos which may suggest an ancient link between these two groups. We also showed that Senoi and Proto-Malay may be admixtures between Negrito and EA populations. Formal admixture tests provided evidence of gene flow between Austro-Asiatic-speaking OAs and populations from Southeast Asia (SEA) and South China which suggest a widespread presence of these people in SEA before Austronesian expansion. Elevated linkage disequilibrium (LD) and enriched homozygosity found in OAs reflect isolation and bottlenecks experienced. Estimates based on N e and LD indicated that these populations diverged from East Asians during the late Pleistocene (14.5 to 8 KYA). The continuum in divergence time from Negritos to Senoi and Proto-Malay in combination with ancestral markers provides evidences of multiple waves of migration into SEA starting with the first Out-of-Africa dispersals followed by Early Train and subsequent Austronesian expansions.

Description: Recent metagenomic studies of environments, such as marine and soil, have significantly enhanced our understanding of the diverse microbial communities living in these habitats and their essential roles in sustaining vast ecosystems. The increase in the number of publications related to soil and marine metagenomics is in sharp contrast to those of air, yet airborne microbes are thought to have significant impacts on many aspects of our lives from their potential roles in atmospheric events such as cloud formation, precipitation, and atmospheric chemistry to their major impact on human health. In this review, we will discuss the current progress in airborne metagenomics, with a special focus on exploring the challenges and opportunities of undertaking such studies. The main challenges of conducting metagenomic studies of airborne microbes are as follows: 1) Low density of microorganisms in the air, 2) efficient retrieval of microorganisms from the air, 3) variability in airborne microbial community composition, 4) the lack of standardized protocols and methodologies, and 5) DNA sequencing and bioinformatics-related challenges. Overcoming these challenges could provide the groundwork for comprehensive analysis of airborne microbes and their potential impact on the atmosphere, global climate, and our health. Metagenomic studies offer a unique opportunity to examine viral and bacterial diversity in the air and monitor their spread locally or across the globe, including threats from pathogenic microorganisms. Airborne metagenomic studies could also lead to discoveries of novel genes and metabolic pathways relevant to meteorological and industrial applications, environmental bioremediation, and biogeochemical cycles.

Description: Within plastid-bearing species, the mutation rate of the plastid genome is often assumed to be greater than that of the mitochondrial genome. This assumption is based on early, pioneering studies of land plant molecular evolution, which uncovered higher rates of synonymous substitution in plastid versus mitochondrial DNAs. However, much of the plastid-containing eukaryotic diversity falls outside of land plants, and the patterns of plastid DNA evolution for embryophytes do not necessarily reflect those of other groups. Recent analyses of plastid and mitochondrial substitution rates in diverse lineages have uncovered very different trends than those recorded for land plants. Here, I explore these new data and argue that for many protists the plastid mutation rate is lower than that of the mitochondrion, including groups with primary or secondary plastids as well as nonphotosynthetic algae. These findings have far-reaching implications for how we view plastid genomes and how their sequences are used for evolutionary analyses, and might ultimately reflect a general tendency toward more efficient DNA repair mechanisms in plastids than in mitochondria.

Description: Understanding how gene regulation evolves is a key area in the current evolutionary field. Gene regulation occurs at various levels. Previous work on the evolution of gene regulation has largely focused on gene transcription. In this study, we used a recently developed ribosomal footprint profiling method to investigate how gene regulation evolves at both the transcription (mRNA abundance) and translation (ribosomal density) levels. By constructing a hybrid between Saccharomyces cerevisiae (Scer) and Saccharomyces bayanus (Sbay), which diverged ~20 Ma, and quantifying transcriptome and translatome in both parental strains and their hybrid, we showed that translation is much more conserved than transcription, mostly due to the buffering effect of translational regulation for the transcriptional divergence. More conservation in translation than transcription is also confirmed by the inheritance mode of transcription and translation between two species. Furthermore, cis and trans effects are widely involved in changes at both transcription and translation levels. Finally, our results showed that genes with certain functions and sequence features might employ specific modes for evolution at these two critical levels of gene regulation. Our results demonstrated that it is essential to investigate the evolution of gene regulation at various levels from different genetic backgrounds to obtain a complete picture of its evolutionary modes in nature.

Description: What makes us human is one of the most interesting and enduring questions in evolutionary biology. To assist in answering this question, we have identified insertions in the human genome which cannot be found in five comparison primate species: Chimpanzee, gorilla, orangutan, gibbon, and macaque. A total of 21,269 nonpolymorphic human-specific insertions were identified, of which only 372 were found in exons. Any function conferred by the remaining 20,897 is likely to be regulatory. Many of these insertions are likely to have been fitness neutral; however, a small number has been identified in genes showing signs of positive selection. Insertions found within positively selected genes show associations to neural phenotypes, which were also enriched in the whole data set. Other phenotypes that are found to be enriched in the data set include dental and sensory perception-related phenotypes, features which are known to differ between humans and other apes. The analysis provides several likely candidates, either genes or regulatory regions, which may be involved in the processes that differentiate humans from other apes.

Description: Bacterial endosymbiosis is an important evolutionary process in insects, which can harbor both obligate and facultative symbionts. The evolution of these symbionts is driven by evolutionary convergence, and they exhibit among the tiniest genomes in prokaryotes. The large host spectrum of facultative symbionts and the high diversity of strategies they use to infect new hosts probably impact the evolution of their genome and explain why they undergo less severe genomic erosion than obligate symbionts. Candidatus Hamiltonella defensa is suitable for the investigation of the genomic evolution of facultative symbionts because the bacteria are engaged in specific relationships in two clades of insects. In aphids, H. defensa is found in several species with an intermediate prevalence and confers protection against parasitoids. In whiteflies, H. defensa is almost fixed in some species of Bemisia tabaci , which suggests an important role of and a transition toward obligate symbiosis. In this study, comparisons of the genome of H. defensa present in two B. tabaci species (Middle East Asia Minor 1 and Mediterranean) and in the aphid Acyrthosiphon pisum revealed that they belong to two distinct clades and underwent specific gene losses. In aphids, it contains highly virulent factors that could allow protection and horizontal transfers. In whiteflies, the genome lost these factors and seems to have a limited ability to acquire genes. However it contains genes that could be involved in the production of essential nutrients, which is consistent with a primordial role for this symbiont. In conclusion, although both lineages of H. defensa have mutualistic interactions with their hosts, their genomes follow distinct evolutionary trajectories that reflect their phenotype and could have important consequences on their evolvability.

Description: By phylogenetic analysis, Mycobacterium kansasii is closely related to Mycobacterium tuberculosis. Yet, although both organisms cause pulmonary disease, M. tuberculosis is a global health menace, whereas M. kansasii is an opportunistic pathogen. To illuminate the differences between these organisms, we have sequenced the genome of M. kansasii ATCC 12478 and its plasmid (pMK12478) and conducted side-by-side in vitro and in vivo investigations of these two organisms. The M. kansasii genome is 6,432,277 bp, more than 2 Mb longer than that of M. tuberculosis H37Rv, and the plasmid contains 144,951 bp. Pairwise comparisons reveal conserved and discordant genes and genomic regions. A notable example of genomic conservation is the virulence locus ESX-1, which is intact and functional in the low-virulence M. kansasii, potentially mediating phagosomal disruption. Differences between these organisms include a decreased predicted metabolic capacity, an increased proportion of toxin–antitoxin genes, and the acquisition of M. tuberculosis -specific genes in the pathogen since their common ancestor. Consistent with their distinct epidemiologic profiles, following infection of C57BL/6 mice, M. kansasii counts increased by less than 10-fold over 6 weeks, whereas M. tuberculosis counts increased by over 10,000-fold in just 3 weeks. Together, these data suggest that M. kansasii can serve as an image of the environmental ancestor of M. tuberculosis before its emergence as a professional pathogen, and can be used as a model organism to study the switch from an environmental opportunistic pathogen to a professional host-restricted pathogen.

Description: Epigenetic control of gene expression has important implications for the regulation of developmental processes, for mediating homeostasis and responses to the external environment, and for transgenerational inheritance of gene expression patterns. Genes that mediate epigenetic control have been well-characterized in Drosophila melanogaster , and we have identified and analyzed an orthologous gene ensemble in Anopheles gambiae that comprises 169 orthologs related to a 215-member epigenetic gene ensemble in D. melanogaster . We find that this ensemble is highly conserved among anopheline mosquitoes, as we identify only seven gene family expansion/contraction events within the ensemble among 12 mosquito species we have studied within the genus Anopheles . Comparative analyses of the epigenetic gene expression across the genera Drosophila and Anopheles reveal distinct tissue-associated expression patterns in the two genera, but similar temporal expression patterns. The A. gambiae complex and D. melanogaster subgroup epigenetic gene ensembles exhibit similar evolutionary rates, as assessed by their respective d N /d S values. These differences in tissue-associated expression patterns, in contrast to similarities in evolutionary rates and temporal expression patterns, may imply that some members of the epigenetic gene ensemble have been redeployed within one or both genera, in comparison to the most recent common ancestor of these two clades. Members of this epigenetic gene ensemble may constitute another set of potential targets for vector control and enable further reductions in the burden of human malaria, by analogy to recent success in development of small molecule antagonists for mammalian epigenetic machinery.

Description: Hirsutella thompsonii ( Ht ) is a fungal pathogen of acarines and the primary cause of epizootics among mites. The draft genomes of two isolates of Ht (MTCC 3556: Ht3 , 34.6 Mb and MTCC 6686: Ht6, 34.7 Mb) are presented and compared with the genomes of Beauveria bassiana ( Bb ) ARSEF 2860 and Ophiocordyceps sinensis ( Os ) CO18. Comparative analysis of carbohydrate active enzymes, pathogen–host interaction genes, metabolism-associated genes, and genes involved in biosynthesis of secondary metabolites in the four genomes was carried out. Reduction in gene family sizes in Ht3 and Os as compared with Ht6 and Bb is observed. Analysis of the mating type genes in Ht reveals the presence of MAT idiomorphs which is suggestive of cryptic sexual traits in Ht . We further identify and classify putative chitinases that may function as virulence factors in fungal entomopathogens due to their role in degradation of arthropod cuticle.

Description: Nucleotide substitutions, small (〈50 bp) insertions or deletions (indels), and large (〉50 bp) deletions are well-known causes of genetic variation within the human genome. We recently reported a previously unrecognized form of polymorphic insertions, termed templated sequence insertion polymorphism (TSIP), in which the inserted sequence was templated from a distant genomic region, and was inserted in the genome through reverse transcription of an RNA intermediate. TSIPs can be grouped into two classes based on nucleotide sequence features at the insertion junctions; class 1 TSIPs show target site duplication, polyadenylation, and preference for insertion at a 5'-TTTT/A-3' sequence, suggesting a LINE-1 based insertion mechanism, whereas class 2 TSIPs show features consistent with repair of a DNA double strand break by nonhomologous end joining. To gain a more complete picture of TSIPs throughout the human population, we evaluated whole-genome sequence from 52 individuals, and identified 171 TSIPs. Most individuals had 25–30 TSIPs, and common (present in 〉20% of individuals) TSIPs were found in individuals throughout the world, whereas rare TSIPs tended to cluster in specific geographic regions. The number of rare TSIPs was greater than the number of common TSIPs, suggesting that TSIP generation is an ongoing process. Intriguingly, mitochondrial sequences were a frequent template for class 2 insertions, used more commonly than any nuclear chromosome. Similar to single nucleotide polymorphisms and indels, we suspect that these TSIPs may be important for the generation of human diversity and genetic diseases, and can be useful in tracking historical migration of populations.

Description: The large number of complete mitochondrial DNA (mtDNA) sequences available for metazoan species makes it a good system for studying genome diversity, although little is known about the mechanisms that promote and/or are correlated with the evolution of this organellar genome. By investigating the molecular evolutionary history of the catalytic and accessory subunits of the mtDNA polymerase, pol , we sought to develop mechanistic insight into its function that might impact genome structure by exploring the relationships between DNA replication and animal mitochondrial genome diversity. We identified three evolutionary patterns among metazoan pol s. First, a trend toward stabilization of both sequence and structure occurred in vertebrates, with both subunits evolving distinctly from those of other animal groups, and acquiring at least four novel structural elements, the most important of which is the HLH-3β (helix-loop-helix, 3 β-sheets) domain that allows the accessory subunit to homodimerize. Second, both subunits of arthropods and tunicates have become shorter and evolved approximately twice as rapidly as their vertebrate homologs. And third, nematodes have lost the gene for the accessory subunit, which was accompanied by the loss of its interacting domain in the catalytic subunit of pol , and they show the highest rate of molecular evolution among all animal taxa. These findings correlate well with the mtDNA genomic features of each group described above, and with their modes of DNA replication, although a substantive amount of biochemical work is needed to draw conclusive links regarding the latter. Describing the parallels between evolution of pol and metazoan mtDNA architecture may also help in understanding the processes that lead to mitochondrial dysfunction and to human disease-related phenotypes.

Description: The ascomycete fungus Ophiostoma novo-ulmi is responsible for the pandemic of Dutch elm disease that has been ravaging Europe and North America for 50 years. We proceeded to annotate the genome of the O. novo-ulmi strain H327 that was sequenced in 2012. The 31.784-Mb nuclear genome (50.1% GC) is organized into 8 chromosomes containing a total of 8,640 protein-coding genes that we validated with RNA sequencing analysis. Approximately 53% of these genes have their closest match to Grosmannia clavigera kw1407, followed by 36% in other close Sordariomycetes, 5% in other Pezizomycotina, and surprisingly few (5%) orphans. A relatively small portion (~3.4%) of the genome is occupied by repeat sequences; however, the mechanism of repeat-induced point mutation appears active in this genome. Approximately 76% of the proteins could be assigned functions using Gene Ontology analysis; we identified 311 carbohydrate-active enzymes, 48 cytochrome P450s, and 1,731 proteins potentially involved in pathogen–host interaction, along with 7 clusters of fungal secondary metabolites. Complementary mating-type locus sequencing, mating tests, and culturing in the presence of elm terpenes were conducted. Our analysis identified a specific genetic arsenal impacting the sexual and vegetative growth, phytopathogenicity, and signaling/plant–defense–degradation relationship between O. novo-ulmi and its elm host and insect vectors.

Description: The vomeronasal organ (VNO) is an olfactory structure for the detection of pheromones. VNO neurons express three groups of unrelated G-protein-coupled receptors. Type-2 vomeronasal receptors (V2Rs) are specifically localized in the basal neurons of the VNO and are believed to sense protein pheromones eliciting specific reproductive behaviors. In murine species, V2Rs are organized into four families. Family-ABD V2Rs are expressed monogenically and coexpress with family-C V2Rs of either subfamily C1 (V2RC1) or subfamily C2 (V2RC2), according to a coordinate temporal diagram. Neurons expressing the phylogenetically ancient V2RC1 coexpress family-BD V2Rs or a specific group of subfamily-A V2Rs (V2RA8-10), whereas a second neuronal subset (V2RC2-positive) coexpresses a recently expanded group of five subfamily-A V2Rs (V2RA1-5) along with vomeronasal-specific Major Histocompatibility Complex molecules (H2-Mv). Through database mining and Sanger sequencing, we have analyzed the onset, diversification, and expansion of the V2R-families throughout the phylogeny of Rodentia. Our results suggest that the separation of V2RC1 and V2RC2 occurred in a Cricetidae ancestor in coincidence with the evolution of the H2-Mv genes; this phylogenetic event did not correspond with the origin of the coexpressing V2RA1-5 genes, which dates back to an ancestral myomorphan lineage. Interestingly, the evolution of receptors within the V2RA1-5 group may be implicated in the origin and diversification of some of the V2R putative cognate ligands, the exocrine secreting peptides. The establishment of V2RC2, which probably reflects the complex expansion and diversification of family-A V2Rs, generated receptors that have probably acquired a more subtle functional specificity.

Description: The increase in phenotypic variability through gene expression noise is proposed to be an evolutionary strategy in selective environments. Differences in promoter-mediated noise between Saccharomyces cerevisiae strains could have been selected for thanks to the benefit conferred by gene expression heterogeneity in the stressful conditions, for instance, those experienced by industrial strains. Here, we used a genome-wide approach to identify promoters conferring high noise levels in the industrial wine strain EC1118. Many promoters of genes related to environmental factors were identified, some of them containing genetic variations compared with their counterpart in the laboratory strain S288c. Each variant of eight promoters has been fused to yeast-Enhanced Green Fluorescent Protein and integrated in the genome of both strains. Some industrial variants conferred higher expression associated, as expected, with lower noise, but other variants either increased or decreased expression without modifying variability, so that they might exhibit different levels of transcriptional-mediated noise at equal mean. At different induction conditions giving similar expression for both variants of the CUP1 promoter, we indeed observed higher noise with the industrial variant. Nevertheless, this difference was only observed in the industrial strain, revealing epistasis in the generation of promoter-mediated noise. Moreover, the increased expression variability conferred by this natural yeast promoter variant provided a clear benefit in the face of an environmental stress. Thus, modulation of gene expression noise by a combination of promoter modifications and trans -influences might be a possible adaptation mechanism in yeast.

Description: The ancestors to the Australian marsupials entered Australia around 60 (54–72) Ma from Antarctica, and radiated into the four living orders Peramelemorphia, Dasyuromorphia, Diprotodontia, and Notoryctemorphia. The relationship between the four Australian marsupial orders has been a long-standing question, because different phylogenetic studies have not been able to consistently reconstruct the same topology. Initial in silico analysis of the Tasmanian devil genome and experimental screening in the seven marsupial orders revealed 20 informative transposable element insertions for resolving the inter- and intraordinal relationships of Australian and South American orders. However, the retrotransposon insertions support three conflicting topologies regarding Peramelemorphia, Dasyuromorphia, and Notoryctemorphia, indicating that the split between the three orders may be best understood as a network. This finding is supported by a phylogenetic reanalysis of nuclear gene sequences, using a consensus network approach that allows depicting hidden phylogenetic conflict, otherwise lost when forcing the data into a bifurcating tree. The consensus network analysis agrees with the transposable element analysis in that all possible topologies regarding Peramelemorphia, Dasyuromorphia, and Notoryctemorphia in a rooted four-taxon topology are equally well supported. In addition, retrotransposon insertion data support the South American order Didelphimorphia being the sistergroup to all other living marsupial orders. The four Australian orders originated within 3 Myr at the Cretaceous–Paleogene boundary. The rapid divergences left conflicting phylogenetic information in the genome possibly generated by incomplete lineage sorting or introgressive hybridization, leaving the relationship among Australian marsupial orders unresolvable as a bifurcating process millions of years later.

Description: Viruses with single-stranded (ss) DNA genomes infect hosts in all three domains of life and include many medically, ecologically, and economically important pathogens. Recently, a new group of ssDNA viruses with chimeric genomes has been discovered through viral metagenomics. These chimeric viruses combine capsid protein genes and replicative protein genes that, respectively, appear to have been inherited from viruses with positive-strand RNA genomes, such as tombusviruses, and ssDNA genomes, such as circoviruses, nanoviruses or geminiviruses. Here, we describe the genome sequence of a new representative of this virus group and reveal an additional layer of chimerism among ssDNA viruses. We show that not only do these viruses encompass genes for capsid proteins and replicative proteins that have distinct evolutionary histories, but also the replicative genes themselves are chimeras of functional domains inherited from viruses of different families. Our results underscore the importance of horizontal gene transfer in the evolution of ssDNA viruses and the role of genetic recombination in the emergence of novel virus groups.

Description: The coat protein complex II (COPII) is responsible for the transport of protein cargoes from the Endoplasmic Reticulum (ER) to the Golgi apparatus. COPII has been functionally characterized extensively in vivo in humans and yeast. This complex shares components with the nuclear pore complex and the Seh1-Associated (SEA) complex, inextricably linking its evolution with that of the nuclear pore and other protocoatomer domain-containing complexes. Importantly, this is one of the last coat complexes to be examined from a comparative genomic and phylogenetic perspective. We use homology searching of eight components across 74 eukaryotic genomes, followed by phylogenetic analyses, to assess both the distribution of the COPII components across eukaryote diversity and to assess its evolutionary history. We report that Sec12, but not Sed4 was present in the Last Eukaryotic Common Ancestor along with Sec16, Sar1, Sec13, Sec31, Sec23, and Sec24. We identify a previously undetected paralog of Sec23 that, at least, predates the archaeplastid clade. We also describe three Sec24 paralogs likely present in the Last Eukaryotic Common Ancestor, including one newly detected that was anciently present but lost from both opisthokonts and excavates. Altogether, we report previously undescribed complexity of the COPII coat in the ancient eukaryotic ancestor and speculate on models for the evolution, not only of the complex, but its relationship to other protocoatomer-derived complexes.

Description: Both hybridization and allopolyploidization generate novel phenotypes by conciliating divergent genomes and regulatory networks in the same cellular context. To understand the rewiring of gene expression in hybrids, the total expression of 21,025 genes and the allele-specific expression of over 11,000 genes were quantified in interspecific hybrids and their parental species, Coffea canephora and Coffea eugenioides using RNA-seq technology. Between parental species, cis - and trans -regulatory divergences affected around 32% and 35% of analyzed genes, respectively, with nearly 17% of them showing both. The relative importance of trans -regulatory divergences between both species could be related to their low genetic divergence and perennial habit. In hybrids, among divergently expressed genes between parental species and hybrids, 77% was expressed like one parent (expression level dominance), including 65% like C. eugenioides . Gene expression was shown to result from the expression of both alleles affected by intertwined parental trans- regulatory factors. A strong impact of C. eugenioides trans -regulatory factors on the upregulation of C. canephora alleles was revealed. The gene expression patterns appeared determined by complex combinations of cis - and trans -regulatory divergences. In particular, the observed biased expression level dominance seemed to be derived from the asymmetric effects of trans -regulatory parental factors on regulation of alleles. More generally, this study illustrates the effects of divergent trans -regulatory parental factors on the gene expression pattern in hybrids. The characteristics of the transcriptional response to hybridization appear to be determined by the compatibility of gene regulatory networks and therefore depend on genetic divergences between the parental species and their evolutionary history.

Description: We study genome-wide nucleotide diversity in three subspecies of extant chimpanzees using exome capture. After strict filtering, Single Nucleotide Polymorphisms and indels were called and genotyped for greater than 50% of exons at a mean coverage of 35 x per individual. Central chimpanzees ( Pan troglodytes troglodytes ) are the most polymorphic (nucleotide diversity, w = 0.0023 per site) followed by Eastern ( P. t. schweinfurthii ) chimpanzees ( w = 0.0016) and Western ( P. t. verus ) chimpanzees ( w = 0.0008). A demographic scenario of divergence without gene flow fits the patterns of autosomal synonymous nucleotide diversity well except for a signal of recent gene flow from Western into Eastern chimpanzees. The striking contrast in X-linked versus autosomal polymorphism and divergence previously reported in Central chimpanzees is also found in Eastern and Western chimpanzees. We show that the direction of selection statistic exhibits a strong nonmonotonic relationship with the strength of purifying selection S , making it inappropriate for estimating S . We instead use counts in synonymous versus nonsynonymous frequency classes to infer the distribution of S coefficients acting on nonsynonymous mutations in each subspecies. The strength of purifying selection we infer is congruent with the differences in effective sizes of each subspecies: Central chimpanzees are undergoing the strongest purifying selection followed by Eastern and Western chimpanzees. Coding indels show stronger selection against indels changing the reading frame than observed in human populations.

Description: Throughout their evolutionary history, genomes acquire new genetic material that facilitates phenotypic innovation and diversification. Developmental processes associated with reproduction are particularly likely to involve novel genes. Abundant gene creation impacts the evolution of chromosomal gene content and general regulatory mechanisms such as dosage compensation. Numerous studies in model organisms have found complex and, at times contradictory, relationships among these genomic attributes highlighting the need to examine these patterns in other systems characterized by abundant sexual selection. Therefore, we examined the association among novel gene creation, tissue-specific gene expression, and chromosomal gene content within stalk-eyed flies. Flies in this family are characterized by strong sexual selection and the presence of a newly evolved X chromosome. We generated RNA-seq transcriptome data from the testes for three species within the family and from seven additional tissues in the highly dimorphic species, Teleopsis dalmanni . Analysis of dipteran gene orthology reveals dramatic testes-specific gene creation in stalk-eyed flies, involving numerous gene families that are highly conserved in other insect groups. Identification of X-linked genes for the three species indicates that the X chromosome arose prior to the diversification of the family. The most striking feature of this X chromosome is that it is highly masculinized, containing nearly twice as many testes-specific genes as expected based on its size. All the major processes that may drive differential sex chromosome gene content—creation of genes with male-specific expression, development of male-specific expression from pre-existing genes, and movement of genes with male-specific expression—are elevated on the X chromosome of T. dalmanni . This masculinization occurs despite evidence that testes expressed genes do not achieve the same levels of gene expression on the X chromosome as they do on the autosomes.

Description: For 200 years, zoologists have relied on phenotypes to learn about the evolution of animals. A glance at the genotype, even through several gene markers, revolutionized our understanding of animal phylogeny. Recent advances in sequencing techniques allow researchers to study speciation mechanisms and the link between genotype and phenotype using complete genomes. We sequenced and assembled a complete genome of the Cloudless Sulphur ( Phoebis sennae ) from a single wild-caught specimen. This genome was used as reference to compare genomes of six specimens, three from the eastern populations (Oklahoma and north Texas), referred to as a subspecies Phoebis sennae eubule , and three from the southwestern populations (south Texas) known as a subspecies Phoebis sennae marcellina . While the two subspecies differ only subtly in phenotype and mitochondrial DNA, comparison of their complete genomes revealed consistent and significant differences, which are more prominent than those between tiger swallowtails Pterourus canadensis and Pterourus glaucus . The two sulphur taxa differed in histone methylation regulators, chromatin-associated proteins, circadian clock, and early development proteins. Despite being well separated on the whole-genome level, the two taxa show introgression, with gene flow mainly from P. s. marcellina to P. s. eubule . Functional analysis of introgressed genes reveals enrichment in transmembrane transporters. Many transporters are responsible for nutrient uptake, and their introgression may be of selective advantage for caterpillars to feed on more diverse food resources. Phylogenetically, complete genomes place family Pieridae away from Papilionidae, which is consistent with previous analyses based on several gene markers.

Description: In eukaryotes, genes are nonrandomly organized into short gene-dense regions or "gene-clusters" interspersed by long gene-poor regions. How these gene-clusters have evolved is not entirely clear. Gene duplication may not account for all the gene-clusters since the genes in most of the clusters do not exhibit significant sequence similarity. In this study, using genome-wide data sets from budding yeast, fruit-fly, and human, we show that: 1) long-range evolutionary repositioning of genes strongly associate with their spatial proximity in the nucleus; 2) presence of evolutionary DNA break-points at involved loci hints at their susceptibility to undergo long-range genomic rearrangements; and 3) correlated epigenetic and transcriptional states of engaged genes highlight the underlying evolutionary constraints. The significance of observation 1, 2, and 3 are particularly stronger for the instances of inferred evolutionary gain, as compared with loss, of linear gene-clustering. These observations suggest that the long-range genomic rearrangements guided through 3D genome organization might have contributed to the evolution of gene order. We further hypothesize that the evolution of linear gene-clusters in eukaryotic genomes might have been mediated through spatial interactions among distant loci in order to optimize co-ordinated regulation of genes. We model this hypothesis through a heuristic model of gene-order evolution.

Description: In recent years a number of studies have brought attention to the role of positive selection during the evolution of antigenic escape by influenza strains. Particularly, the identification of positively selected sites within antigenic domains of viral surface proteins has been used to suggest that the evolution of viral–host receptor binding specificity is driven by selection. Here we show that, following the 1918 outbreak, the antigenic sites of the hemagglutinin (HA) viral surface protein and the stalk region of neuraminidase became substitution hotspots. The hotspots show similar patterns of nucleotide substitution bias at synonymous and nonsynonymous sites. Such bias imposes directionality in amino acid replacements that can influence signals of selection at antigenic sites. Our results suggest that the high accumulation of substitutions within the antigenic sites of HA can explain not only cases of antigenic escape by antigenic drift but also lead to occasional episodes of viral extinction.

Description: Mycobacterium phlei , a nontuberculosis mycobacterial species, was first described in 1898–1899. We present the complete genome sequence for the M. phlei CCUG21000 T type strain and the draft genomes for four additional strains. The genome size for all fiveis 5.3 Mb with 69.4% Guanine-Cytosine content. This is 0.35 Mbp smaller than the previously reported M. phlei RIVM draft genome. The size difference is attributed partly to large bacteriophage sequence fragments in the M. phlei RIVM genome. Comparative analysis revealed the following: 1) A CRISPR system similar to Type 1E ( cas 3) in M. phlei RIVM; 2) genes involved in polyamine metabolism and transport ( potAD , potF ) that are absent in other mycobacteria, and 3) strain-specific variations in the number of -factor genes. Moreover, M. phlei has as many as 82 mce (mammalian cell entry) homologs and many of the horizontally acquired genes in M. phlei are present in other environmental bacteria including mycobacteria that share similar habitat. Phylogenetic analysis based on 693 Mycobacterium core genes present in all complete mycobacterial genomes suggested that its closest neighbor is Mycobacterium smegmatis JS623 and Mycobacterium rhodesiae NBB3, while it is more distant to M. smegmatis mc2 155.

Description: For both historical and technical reasons, 16S ribosomal RNA has been the most common molecular marker used to analyze the contents of microbial communities. However, its slow rate of evolution hinders the resolution of closely related bacteria—individual 16S-phylotypes, particularly when clustered at 97% sequence identity, conceal vast amounts of species- and strain-level variation. Protein-coding genes, which evolve more quickly, are useful for differentiating among more recently diverged lineages, but their application is complicated by difficulties in designing low-redundancy primers that amplify homologous regions from distantly related taxa. Given the now-common practice of multiplexing hundreds of samples, adopting new genes usually entails the synthesis of large sets of barcoded primers. To circumvent problems associated with use of protein-coding genes to survey microbial communities, we develop an approach—termed phyloTAGs —that offers an automatic solution for primer design and can be easily adapted to target different taxonomic groups and/or different protein-coding regions. We applied this method to analyze diversity within the gorilla gut microbiome and recovered hundreds of strains that went undetected after deep-sequencing of 16S rDNA amplicons. PhyloTAGs provides a powerful way to recover the fine-level diversity within microbial communities and to study stability and dynamics of bacterial populations.

Description: Here we analyzed whole-genome sequences of 3,969 influenza A(H1N1)pdm09 and 4,774 A(H3N2) strains that circulated during 2009–2015 in the world. The analysis revealed changes at 481 and 533 amino acid sites in proteins of influenza A(H1N1)pdm09 and A(H3N2) strains, respectively. Many of these changes were introduced as a result of random drift. However, there were 61 and 68 changes that were present in relatively large number of A(H1N1)pdm09 and A(H3N2) strains, respectively, that circulated during relatively long time. We named these amino acid substitutions evolutionary markers, as they seemed to contain valuable information regarding the viral evolution. Interestingly, influenza A(H1N1)pdm09 and A(H3N2) viruses acquired non-overlapping sets of evolutionary markers. We next analyzed these characteristic markers in vaccine strains recommended by the World Health Organization for the past five years. Our analysis revealed that vaccine strains carried only few evolutionary markers at antigenic sites of viral hemagglutinin (HA) and neuraminidase (NA). The absence of these markers at antigenic sites could affect the recognition of HA and NA by human antibodies generated in response to vaccinations. This could, in part, explain moderate efficacy of influenza vaccines during 2009–2014. Finally, we identified influenza A(H1N1)pdm09 and A(H3N2) strains, which contain all the evolutionary markers of influenza A strains circulated in 2015, and which could be used as vaccine candidates for the 2015/2016 season. Thus, genome-wide analysis of evolutionary markers of influenza A(H1N1)pdm09 and A(H3N2) viruses may guide selection of vaccine strain candidates.

Description: Genome-wide studies of African populations have the potential to reveal powerful insights into the evolution of our species, as these diverse populations have been exposed to intense selective pressures imposed by infectious diseases, diet, and environmental factors. Within Africa, the Sahel Belt extensively overlaps the geographical center of several endemic infections such as malaria, trypanosomiasis, meningitis, and hemorrhagic fevers. We screened 2.5 million single nucleotide polymorphisms in 161 individuals from 13 Sahelian populations, which together with published data cover Western, Central, and Eastern Sahel, and include both nomadic and sedentary groups. We confirmed the role of this Belt as a main corridor for human migrations across the continent. Strong admixture was observed in both Central and Eastern Sahelian populations, with North Africans and Near Eastern/Arabians, respectively, but it was inexistent in Western Sahelian populations. Genome-wide local ancestry inference in admixed Sahelian populations revealed several candidate regions that were significantly enriched for non-autochthonous haplotypes, and many showed to be under positive selection. The DARC gene region in Arabs and Nubians was enriched for African ancestry, whereas the RAB3GAP1 / LCT / MCM6 region in Oromo, the TAS2R gene family in Fulani, and the ALMS1 / NAT8 in Turkana and Samburu were enriched for non-African ancestry. Signals of positive selection varied in terms of geographic amplitude. Some genomic regions were selected across the Belt, the most striking example being the malaria-related DARC gene. Others were Western-specific (oxytocin, calcium, and heart pathways), Eastern-specific (lipid pathways), or even population-restricted ( TAS2R genes in Fulani, which may reflect sexual selection).

Description: The comparative genomics revolution of the past decade has enabled the discovery of functional elements in the human genome via sequence comparison. While that is so, an important class of elements, those specific to humans, is entirely missed by searching for sequence conservation across species. Here we present an analysis based on variation data among human genomes that utilizes a supervised machine learning approach for the identification of human-specific purifying selection in the genome. Using only allele frequency information from the complete low-coverage 1000 Genomes Project data set in conjunction with a support vector machine trained from known functional and nonfunctional portions of the genome, we are able to accurately identify portions of the genome constrained by purifying selection. Our method identifies previously known human-specific gains or losses of function and uncovers many novel candidates. Candidate targets for gain and loss of function along the human lineage include numerous putative regulatory regions of genes essential for normal development of the central nervous system, including a significant enrichment of gain of function events near neurotransmitter receptor genes. These results are consistent with regulatory turnover being a key mechanism in the evolution of human-specific characteristics of brain development. Finally, we show that the majority of the genome is unconstrained by natural selection currently, in agreement with what has been estimated from phylogenetic methods but in sharp contrast to estimates based on transcriptomics or other high-throughput functional methods.

Description: Previous studies focused on psychrophilic adaptation generally have demonstrated that multiple mechanisms work together to increase protein flexibility and activity, as well as to decrease the thermostability of proteins. However, the relationship between high and low temperature adaptations remains unclear. To investigate this issue, we collected the available predicted whole proteome sequences of species with different optimal growth temperatures, and analyzed amino acid variations and substitutional asymmetry in pairs of homologous proteins from related species. We found that changes in amino acid composition associated with low temperature adaptation did not exhibit a coherent opposite trend when compared with changes in amino acid composition associated with high temperature adaptation. This result indicates that during their evolutionary histories the proteome-scale evolutionary patterns associated with prokaryotes exposed to low temperature environments were distinct from the proteome-scale evolutionary patterns associated with prokaryotes exposed to high temperature environments in terms of changes in amino acid composition of the proteins.

Description: The genome sequences of the plastid and mitochondrion of white spruce ( Picea glauca ) were assembled from whole-genome shotgun sequencing data using ABySS. The sequencing data contained reads from both the nuclear and organellar genomes, and reads of the organellar genomes were abundant in the data as each cell harbors hundreds of mitochondria and plastids. Hence, assembly of the 123-kb plastid and 5.9-Mb mitochondrial genomes were accomplished by analyzing data sets primarily representing low coverage of the nuclear genome. The assembled organellar genomes were annotated for their coding genes, ribosomal RNA, and transfer RNA. Transcript abundances of the mitochondrial genes were quantified in three developmental tissues and five mature tissues using data from RNA-seq experiments. C-to-U RNA editing was observed in the majority of mitochondrial genes, and in four genes, editing events were noted to modify ACG codons to create cryptic AUG start codons. The informatics methodology presented in this study should prove useful to assemble organellar genomes of other plant species using whole-genome shotgun sequencing data.

Description: Evolution and maintenance of genetic recombination and its relation to the mutational process is a long-standing, fundamental problem in evolutionary biology that is linked to the general problem of evolution of evolvability. We explored a stochastic model of the evolution of recombination using additive fitness and infinite allele assumptions but no assumptions on the sign or magnitude of the epistasis and the distribution of mutation effects. In this model, fluctuating negative epistasis and predominantly deleterious mutations arise naturally as a consequence of the additive fitness and a reservoir from which new alleles arrive with a fixed distribution of fitness effects. Analysis of the model revealed a nonmonotonic effect of recombination intensity on fitness, with an optimal recombination rate value which maximized fitness in steady state. The optimal recombination rate depended on the mutation rate and was evolvable, that is, subject to selection. The predictions of the model were compatible with the observations on the dependence between genome rearrangement rate and gene flux in microbial genomes.

Description: Olfactory receptors (ORs) are membrane proteins that mediate the detection of odorants in the environment, and are the largest vertebrate gene family. Comparative studies of mammalian genomes indicate that OR repertoires vary widely, even between closely related lineages, as a consequence of frequent OR gains and losses. Several studies also suggest that mammalian OR repertoires are influenced by life history traits. Sauropsida is a diverse group of vertebrates group that is the sister group to mammals, and includes birds, testudines, squamates, and crocodilians, and represents a natural system to explore predictions derived from mammalian studies. In this study, we analyzed olfactory receptor (OR) repertoire variation among several representative species and found that the number of intact OR genes in sauropsid genomes analyzed ranged over an order of magnitude, from 108 in the green anole to over 1,000 in turtles. Our results suggest that different sauropsid lineages have highly divergent OR repertoire composition that derive from lineage-specific combinations of gene expansions, losses, and retentions of ancestral OR genes. These differences also suggest that varying degrees of adaption related to life history have shaped the unique OR repertoires observed across sauropsid lineages.

Description: In clonal systems, interpreting driver genes in terms of molecular networks helps understanding how these drivers elicit an adaptive phenotype. Obtaining such a network-based understanding depends on the correct identification of driver genes. In clonal systems, independent evolved lines can acquire a similar adaptive phenotype by affecting the same molecular pathways, a phenomenon referred to as parallelism at the molecular pathway level. This implies that successful driver identification depends on interpreting mutated genes in terms of molecular networks. Driver identification and obtaining a network-based understanding of the adaptive phenotype are thus confounded problems that ideally should be solved simultaneously. In this study, a network-based eQTL method is presented that solves both the driver identification and the network-based interpretation problem. As input the method uses coupled genotype-expression phenotype data (eQTL data) of independently evolved lines with similar adaptive phenotypes and an organism-specific genome-wide interaction network. The search for mutational consistency at pathway level is defined as a subnetwork inference problem, which consists of inferring a subnetwork from the genome-wide interaction network that best connects the genes containing mutations to differentially expressed genes. Based on their connectivity with the differentially expressed genes, mutated genes are prioritized as driver genes. Based on semisynthetic data and two publicly available data sets, we illustrate the potential of the network-based eQTL method to prioritize driver genes and to gain insights in the molecular mechanisms underlying an adaptive phenotype. The method is available at http://bioinformatics.intec.ugent.be/phenetic_eqtl/index.html

Description: Facilitating the evolution of new gene functions, gene duplication is a major mechanism driving evolutionary innovation. Gene family expansions relevant to host/symbiont interactions are increasingly being discovered in eukaryotes that host endosymbiotic microbes. Such discoveries entice speculation that gene duplication facilitates the evolution of novel, endosymbiotic relationships. Here, using a comparative transcriptomic approach combined with differential gene expression analysis, we investigate the importance of endosymbiosis in retention of amino acid transporter paralogs in aphid genomes. To pinpoint the timing of amino acid transporter duplications we inferred gene phylogenies for five aphid species and three outgroups. We found that while some duplications arose in the aphid common ancestor concurrent with endosymbiont acquisition, others predate aphid divergence from related insects without intracellular symbionts, and still others appeared during aphid diversification. Interestingly, several aphid-specific paralogs have conserved enriched expression in bacteriocytes, the insect cells that host primary symbionts. Conserved bacteriocyte enrichment suggests that the transporters were recruited to the aphid/endosymbiont interface in the aphid common ancestor, consistent with a role for gene duplication in facilitating the evolution of endosymbiosis in aphids. In contrast, the temporal variability of amino acid transporter duplication indicates that endosymbiosis is not the only trait driving selection for retention of amino acid transporter paralogs in sap-feeding insects. This study cautions against simplistic interpretations of the role of gene family expansion in the evolution of novel host/symbiont interactions by further highlighting that multiple complex factors maintain gene family paralogs in the genomes of eukaryotes that host endosymbiotic microbes.

Description: Tibetan chicken, unlike their lowland counterparts, exhibit specific adaptations to high-altitude conditions. The genetic mechanisms of such adaptations in highland chickens were determined by resequencing the genomes of four highland (Tibetan and Lhasa White) and four lowland (White Leghorn, Lindian, and Chahua) chicken populations. Our results showed an evident genetic admixture in Tibetan chickens, suggesting a history of introgression from lowland gene pools. Genes showing positive selection in highland populations were related to cardiovascular and respiratory system development, DNA repair, response to radiation, inflammation, and immune responses, indicating a strong adaptation to oxygen scarcity and high-intensity solar radiation. The distribution of allele frequencies of nonsynonymous single nucleotide polymorphisms between highland and lowland populations was analyzed using chi-square test, which showed that several differentially distributed genes with missense mutations were enriched in several functional categories, especially in blood vessel development and adaptations to hypoxia and intense radiation. RNA sequencing revealed that several differentially expressed genes were enriched in gene ontology terms related to blood vessel and respiratory system development. Several candidate genes involved in the development of cardiorespiratory system ( FGFR1 , CTGF , ADAM9 , JPH2 , SATB1 , BMP4 , LOX , LPR , ANGPTL4 , and HYAL1 ), inflammation and immune responses ( AIRE , MYO1F , ZAP70 , DDX60 , CCL19 , CD47 , JSC , and FAS ), DNA repair, and responses to radiation ( VCP , ASH2L , and FANCG ) were identified to play key roles in the adaptation to high-altitude conditions. Our data provide new insights into the unique adaptations of highland animals to extreme environments.

Description: We performed a detailed whole-genome sequence analysis of Providencia rettgeri H1736, a multidrug-resistant clinical pathogen isolated in Israel in 2011. The objective was to describe the genomic flexibility of this bacterium that has greatly contributed to its pathogenicity. The genome has a chromosome size of 4,609,352 bp with 40.22% GC content. Five plasmids were predicted, as well as other mobile genetic elements (MGEs) including phages, genomic islands, and integrative and conjugative elements. The resistome consisted of a total of 27 different antibiotic resistance genes including bla NDM-1 , mostly located on MGEs. Phenotypically, the bacteria displayed resistance to a total of ten different antimicrobial classes. Various features such as metabolic operons (including a novel carbapenem biosynthesis operon) and virulence genes were also borne on the MGEs, making P. rettgeri H1736 significantly different from other P. rettgeri isolates. A large quantity of the genetic diversity that exists in P. rettgeri H1736 was due to extensive horizontal gene transfer events, leading to an enormous presence of MGEs in its genome. Most of these changes contributed toward the pathogenic evolution of this bacterium.

Description: Aspergillus ochraceoroseus and Aspergillus rambellii were isolated from soil detritus in Taï National Park, Ivory Coast, Africa. The Type strain for each species happens to be the only representative ever sampled. Both species secrete copious amounts of aflatoxin B 1 and sterigmatocystin, because each of their genomes contains clustered genes for biosynthesis of these mycotoxins. We sequenced their genomes using a personal genome machine and found them to be smaller in size ( A. ochraceoroseus = 23.9 Mb and A. rambellii = 26.1 Mb), as well as in numbers of predicted genes (7,837 and 7,807, respectively), compared to other sequenced Aspergilli. Our findings also showed that the A. ochraceoroseus Type strain contains a single MAT1-1 gene, while the Type strain of A. rambellii contains a single MAT1-2 gene, indicating that these species are heterothallic (self-infertile). These draft genomes will be useful for understanding the genes and pathways necessary for the cosynthesis of these two toxic secondary metabolites as well as the evolution of these pathways in aflatoxigenic fungi.

Description: Escherichia coli strain MRE600 was originally identified for its low RNase I activity and has therefore been widely adopted by the biomedical research community as a preferred source for the expression and purification of transfer RNAs and ribosomes. Despite its widespread use, surprisingly little information about its genome or genetic content exists. Here, we present the first de novo assembly and description of the MRE600 genome and epigenome. To provide context to these studies of MRE600, we include comparative analyses with E. coli K-12 MG1655 (K12). Pacific Biosciences Single Molecule, Real-Time sequencing reads were assembled into one large chromosome (4.83 Mb) and three smaller plasmids (89.1, 56.9, and 7.1 kb). Interestingly, the 7.1-kb plasmid possesses genes encoding a colicin E1 protein and its associated immunity protein. The MRE600 genome has a G + C content of 50.8% and contains a total of 5,181 genes, including 4,913 protein-encoding genes and 268 RNA genes. We identified 41,469 modified DNA bases (0.83% of total) and found that MRE600 lacks the gene for type I methyltransferase, EcoKI. Phylogenetic, taxonomic, and genetic analyses demonstrate that MRE600 is a divergent E. coli strain that displays features of the closely related genus, Shigella . Nevertheless, comparative analyses between MRE600 and E. coli K12 show that these two strains exhibit nearly identical ribosomal proteins, ribosomal RNAs, and highly homologous tRNA species. Substantiating prior suggestions that MRE600 lacks RNase I activity, the RNase I-encoding gene, rna , contains a single premature stop codon early in its open-reading frame.