ALBERT

Description: Admixture mapping has been enormously resourceful in identifying genetic variations linked to phenotypes, adaptation, and diseases. In this study through analysis of copy number variable regions (CNVRs), we report extensive restructuring in the genomes of the recently admixed African-Indian population (OG-W-IP) that inhabits a highly saline environment in Western India. The study included subjects from OG-W-IP (OG), five different Indian and three HapMap populations that were genotyped using Affymetrix version 6.0 arrays. Copy number variations (CNVs) detected using Birdsuite were used to define CNVRs. Population structure with respect to CNVRs was delineated using random forest approach. OG genomes have a surprising excess of CNVs in comparison to other studied populations. Individual ancestry proportions computed using STRUCTURE also reveals a unique genetic component in OGs. Population structure analysis with CNV genotypes indicates OG to be distant from both the African and Indian ancestral populations. Interestingly, it shows genetic proximity with respect to CNVs to only one Indian population IE-W-LP4, which also happens to reside in the same geographical region. We also observe a significant enrichment of molecular processes related to ion binding and receptor activity in genes encompassing OG-specific CNVRs. Our results suggest that retention of CNVRs from ancestral natives and de novo acquisition of CNVRs could accelerate the process of adaptation especially in an extreme environment. Additionally, this population would be enormously useful for dissecting genes and delineating the involvement of CNVs in salt adaptation.

Description: As an important subtype of structural variations, chromosomal translocation is associated with various diseases, especially cancers, by disrupting gene structures and functions. Traditional methods for identifying translocations are time consuming and have limited resolutions. Recently, a few studies have employed next-generation sequencing (NGS) technology for characterizing chromosomal translocations on human genome, obtaining high-throughput results with high resolutions. However, these studies are mainly focused on mechanism-specific or site-specific translocation mapping. In this study, we conducted a comprehensive genome-wide analysis on the characterization of human chromosomal material exchange with regard to the chromosome translocations. Using NGS data of 1,481 subjects from the 1000 Genomes Project, we identified 15,349,092 translocated DNA fragment pairs, ranging from 65 to 1,886 bp and with an average size of approximately 102 bp. On average, each individual genome carried about 10,364 pairs, covering approximately 0.069% of the genome. We identified 16 translocation hot regions, among which two regions did not contain repetitive fragments. Results of our study overlapped with a majority of previous results, containing approximately 79% of approximately 2,340 translocations characterized in three available translocation databases. In addition, our study identified five novel potential recurrent chromosomal material exchange regions with greater than 20% detection rates. Our results will be helpful for an accurate characterization of translocations in human genomes, and contribute as a resource for future studies of the roles of translocations in human disease etiology and mechanisms.

Description: Angiosperm mitochondrial genomes exhibit many unusual properties, including heterogeneous nucleotide composition and exceptionally large and variable genome sizes. Determining the role of nonadaptive mechanisms such as mutation bias in shaping the molecular evolution of these unique genomes has proven challenging because their dynamic structures generally prevent identification of homologous intergenic sequences for comparative analyses. Here, we report an analysis of angiosperm mitochondrial DNA sequences that are derived from inserted plastid DNA ( mtpts ). The availability of numerous completely sequenced plastid genomes allows us to infer the evolutionary history of these insertions, including the specific nucleotide substitutions and indels that have occurred because their incorporation into the mitochondrial genome. Our analysis confirmed that many mtpts have a complex history, including frequent gene conversion and multiple examples of horizontal transfer between divergent angiosperm lineages. Nevertheless, it is clear that the majority of extant mtpt sequence in angiosperms is the product of recent transfer (or gene conversion) and is subject to rapid loss/deterioration, suggesting that most mtpts are evolving relatively free from functional constraint. The evolution of mtpt sequences reveals a pattern of biased mutational input in angiosperm mitochondrial genomes, including an excess of small deletions over insertions and a skew toward nucleotide substitutions that increase AT content. However, these mutation biases are far weaker than have been observed in many other cellular genomes, providing insight into some of the notable features of angiosperm mitochondrial architecture, including the retention of large intergenic regions and the relatively neutral GC content found in these regions.

Description: Neotropical primates (NP) are presently distributed in the New World from Mexico to northern Argentina, comprising three large families, Cebidae, Atelidae, and Pitheciidae, consequently to their diversification following their separation from Old World anthropoids near the Eocene/Oligocene boundary, some 40 Ma. The evolution of NP has been intensively investigated in the last decade by studies focusing on their phylogeny and timescale. However, despite major efforts, the phylogenetic relationship between these three major clades and the age of their last common ancestor are still controversial because these inferences were based on limited numbers of loci and dating analyses that did not consider the evolutionary variation associated with the distribution of gene trees within the proposed phylogenies. We show, by multispecies coalescent analyses of selected genome segments, spanning along 92,496,904 bp that the early diversification of extant NP was marked by a 2-fold increase of their effective population size and that Atelids and Cebids are more closely related respective to Pitheciids. The molecular phylogeny of NP has been difficult to solve because of population-level phenomena at the early evolution of the lineage. The association of evolutionary variation with the distribution of gene trees within proposed phylogenies is crucial for distinguishing the mean genetic divergence between species (the mean coalescent time between loci) from speciation time. This approach, based on extensive genomic data provided by new generation DNA sequencing, provides more accurate reconstructions of phylogenies and timescales for all organisms.

Description: Dogs shared a much closer relationship with humans than any other domesticated animals, probably due to their unique social cognitive capabilities, which were hypothesized to be a by-product of selection for tameness toward humans. Here, we demonstrate that genes involved in glutamate metabolism, which account partially for fear response, indeed show the greatest population differentiation by whole-genome comparison of dogs and wolves. However, the changing direction of their expression supports a role in increasing excitatory synaptic plasticity in dogs rather than reducing fear response. Because synaptic plasticity are widely believed to be cellular correlates of learning and memory, this change may alter the learning and memory abilities of ancient scavenging wolves, weaken the fear reaction toward humans, and prompt the initial interspecific contact.

Description: Hirsutella minnesotensis [Ophiocordycipitaceae (Hypocreales, Ascomycota)] is a dominant endoparasitic fungus by using conidia that adhere to and penetrate the secondary stage juveniles of soybean cyst nematode. Its genome was de novo sequenced and compared with five entomopathogenic fungi in the Hypocreales and three nematode-trapping fungi in the Orbiliales (Ascomycota). The genome of H. minnesotensis is 51.4 Mb and encodes 12,702 genes enriched with transposable elements up to 32%. Phylogenomic analysis revealed that H. minnesotensis was diverged from entomopathogenic fungi in Hypocreales. Genome of H. minnesotensis is similar to those of entomopathogenic fungi to have fewer genes encoding lectins for adhesion and glycoside hydrolases for cellulose degradation, but is different from those of nematode-trapping fungi to possess more genes for protein degradation, signal transduction, and secondary metabolism. Those results indicate that H. minnesotensis has evolved different mechanism for nematode endoparasitism compared with nematode-trapping fungi. Transcriptomics analyses for the time-scale parasitism revealed the upregulations of lectins, secreted proteases and the genes for biosynthesis of secondary metabolites that could be putatively involved in host surface adhesion, cuticle degradation, and host manipulation. Genome and transcriptome analyses provided comprehensive understanding of the evolution and lifestyle of nematode endoparasitism.

Description: Anopheles gambiae is a major mosquito vector of malaria in Africa. Although increased use of insecticide-based vector control tools has decreased malaria transmission, elimination is likely to require novel genetic control strategies. It can be argued that the absence of an A. gambiae inbred line has slowed progress toward genetic vector control . In order to empower genetic studies and enable precise and reproducible experimentation, we set out to create an inbred line of this species. We found that amenability to inbreeding varied between populations of A. gambiae . After full-sib inbreeding for ten generations, we genotyped 112 individuals—56 saved prior to inbreeding and 56 collected after inbreeding—at a genome-wide panel of single nucleotide polymorphisms (SNPs). Although inbreeding dramatically reduced diversity across much of the genome, we discovered numerous, discrete genomic blocks that maintained high heterozygosity. For one large genomic region, we were able to definitively show that high diversity is due to the persistent polymorphism of a chromosomal inversion. Inbred lines in other eukaryotes often exhibit a qualitatively similar retention of polymorphism when typed at a small number of markers. Our whole-genome SNP data provide the first strong, empirical evidence supporting associative overdominance as the mechanism maintaining higher than expected diversity in inbred lines. Although creation of A. gambiae lines devoid of nearly all polymorphism may not be feasible, our results provide critical insights into how more fully isogenic lines can be created.

Description: Two-pronged bristletails (Diplura) are traditionally classified into three major superfamilies: Campodeoidea, Projapygoidea, and Japygoidea. The interrelationships of these three superfamilies and the monophyly of Diplura have been much debated. Few previous studies included Projapygoidea in their phylogenetic considerations, and its position within Diplura still is a puzzle from both morphological and molecular points of view. Until now, no mitochondrial genome has been sequenced for any projapygoid species. To fill in this gap, we determined and annotated the complete mitochondrial genome of Octostigma sinensis (Octostigmatidae, Projapygoidea), and of three more dipluran species, one each from the Campodeidae, Parajapygidae, and Japygidae. All four newly sequenced dipluran mtDNAs encode the same set of genes in the same gene order as shared by most crustaceans and hexapods. Secondary structure truncations have occurred in trnR , trnC , trnS1 , and trnS2 , and the reduction of transfer RNA D-arms was found to be taxonomically correlated, with Campodeoidea having experienced the most reduction. Partitioned phylogenetic analyses, based on both amino acids and nucleotides of the protein-coding genes plus the ribosomal RNA genes, retrieve significant support for a monophyletic Diplura within Pancrustacea, with Projapygoidea more closely related to Campodeoidea than to Japygoidea. Another key finding is that monophyly of Diplura cannot be recovered unless Projapygoidea is included in the phylogenetic analyses; this explains the dipluran polyphyly found by past mitogenomic studies. Including Projapygoidea increased the sampling density within Diplura and probably helped by breaking up a long-branch-attraction artifact. This finding provides an example of how proper sampling is significant for phylogenetic inference.

Description: Major unresolved questions in evolutionary genetics include determining the contributions of different mutational sources to the total pool of genetic variation in a species, and understanding how these different forms of genetic variation interact with natural selection. Recent work has shown that structural variants (SVs) (insertions, deletions, inversions, and transpositions) are a major source of genetic variation, often outnumbering single nucleotide variants in terms of total bases affected. Despite the near ubiquity of SVs, major questions about their interaction with natural selection remain. For example, how does the allele frequency spectrum of SVs differ when compared with single nucleotide variants? How often do SVs affect genes, and what are the consequences? To begin to address these questions, we have systematically identified and characterized a large set of submicroscopic insertion and deletion (indel) variants (between 1 and 200 kb in length) among ten inbred lines from a single natural population of the plant species Mimulus guttatus . After extensive computational filtering, we focused on a set of 4,142 high-confidence indels that showed an experimental validation rate of 73%. All but one of these indels were less than 200 kb. Although the largest were generally at lower frequencies in the population, a surprising number of large indels are at intermediate frequencies. Although indels overlapping with genes were much rarer than expected by chance, approximately 600 genes were affected by an indel. Nucleotide-binding site leucine-rich repeat (NBS–LRR) defense response genes were the most enriched among the gene families affected. Most indels associated with genes were rare and appeared to be under purifying selection, though we do find four high-frequency derived insertion alleles that show signatures of recent positive selection.

Description: Phylogenetic inference is widely used to investigate the relationships between homologous sequences. RNA molecules have played a key role in these studies because they are present throughout life and tend to evolve slowly. Phylogenetic inference has been shown to be dependent on the substitution model used. A wide range of models have been developed to describe RNA evolution, either with 16 states describing all possible canonical base pairs or with 7 states where the 10 mismatched nucleotides are reduced to a single state. Formal model selection has become a standard practice for choosing an inferential model and works well for comparing models of a specific type, such as comparisons within nucleotide models or within amino acid models. Model selection cannot function across different sized state spaces because the likelihoods are conditioned on different data. Here, we introduce statistical state-space projection methods that allow the direct comparison of likelihoods between nucleotide models and 7-state and 16-state RNA models. To demonstrate the general applicability of our new methods, we extract 287 RNA families from genomic alignments and perform model selection. We find that in 281/287 families, RNA models are selected in preference to nucleotide models, with simple 7-state RNA models selected for more conserved families with shorter stems and more complex 16-state RNA models selected for more divergent families with longer stems. Other factors, such as the function of the RNA molecule or the GC-content, have limited impact on model selection. Our models and model selection methods are freely available in the open-source PHASE 3.0 software.

Description: Symbiotic associations between animals and microbes are ubiquitous in nature, with an estimated 15% of all insect species harboring intracellular bacterial symbionts. Most bacterial symbionts share many genomic features including small genomes, nucleotide composition bias, high coding density, and a paucity of mobile DNA, consistent with long-term host association. In this study, we focus on the early stages of genome degeneration in a recently derived insect-bacterial mutualistic intracellular association. We present the complete genome sequence and annotation of Sitophilus oryzae primary endosymbiont (SOPE). We also present the finished genome sequence and annotation of strain HS, a close free-living relative of SOPE and other insect symbionts of the Sodalis -allied clade, whose gene inventory is expected to closely resemble the putative ancestor of this group. Structural, functional, and evolutionary analyses indicate that SOPE has undergone extensive adaptation toward an insect-associated lifestyle in a very short time period. The genome of SOPE is large in size when compared with many ancient bacterial symbionts; however, almost half of the protein-coding genes in SOPE are pseudogenes. There is also evidence for relaxed selection on the remaining intact protein-coding genes. Comparative analyses of the whole-genome sequence of strain HS and SOPE highlight numerous genomic rearrangements, duplications, and deletions facilitated by a recent expansion of insertions sequence elements, some of which appear to have catalyzed adaptive changes. Functional metabolic predictions suggest that SOPE has lost the ability to synthesize several essential amino acids and vitamins. Analyses of the bacterial cell envelope and genes encoding secretion systems suggest that these structures and elements have become simplified in the transition to a mutualistic association.

Description: Sea ice is a highly dynamic and productive environment that includes a diverse array of psychrophilic prokaryotic and eukaryotic taxa distinct from the underlying water column. Because sea ice has only been extensive on Earth since the mid-Eocene, it has been hypothesized that bacteria highly adapted to inhabit sea ice have traits that have been acquired through horizontal gene transfer (HGT). Here we compared the genomes of the psychrophilic bacterium Psychroflexus torquis ATCC 700755 T , associated with both Antarctic and Arctic sea ice, and its closely related nonpsychrophilic sister species, P. gondwanensis ACAM 44 T . Results show that HGT has occurred much more extensively in P. torquis in comparison to P. gondwanensis . Genetic features that can be linked to the psychrophilic and sea ice-specific lifestyle of P. torquis include genes for exopolysaccharide (EPS) and polyunsaturated fatty acid (PUFA) biosynthesis, numerous specific modes of nutrient acquisition, and proteins putatively associated with ice-binding, light-sensing (bacteriophytochromes), and programmed cell death (metacaspases). Proteomic analysis showed that several genes associated with these traits are highly translated, especially those involved with EPS and PUFA production. Because most of the genes relating to the ability of P. torquis to dwell in sea-ice ecosystems occur on genomic islands that are absent in closely related P. gondwanensis , its adaptation to the sea-ice environment appears driven mainly by HGT. The genomic islands are rich in pseudogenes, insertional elements, and addiction modules, suggesting that gene acquisition is being followed by a process of genome reduction potentially indicative of evolving ecosystem specialism.

Description: Papillomaviruses (PVs) are widespread pathogens. However, the extent of PV infections in bats remains largely unknown. This work represents the first comprehensive study of PVs in Iberian bats. We identified four novel PVs in the mucosa of free-ranging Eptesicus serotinus (EserPV1, EserPV2, and EserPV3) and Rhinolophus ferrumequinum (RferPV1) individuals and analyzed their phylogenetic relationships within the viral family. We further assessed their prevalence in different populations of E. serotinus and its close relative E. isabellinus . Although it is frequent to read that PVs co-evolve with their host, that PVs are highly species-specific, and that PVs do not usually recombine, our results suggest otherwise. First, strict virus–host co-evolution is rejected by the existence of five, distantly related bat PV lineages and by the lack of congruence between bats and bat PVs phylogenies. Second, the ability of EserPV2 and EserPV3 to infect two different bat species ( E. serotinus and E. isabellinus ) argues against strict host specificity. Finally, the description of a second noncoding region in the RferPV1 genome reinforces the view of an increased susceptibility to recombination in the E2-L2 genomic region. These findings prompt the question of whether the prevailing paradigms regarding PVs evolution should be reconsidered.

Description: Sex chromosome dosage compensation balances homogametic sex chromosome expression with autosomal expression in the heterogametic sex, leading to sex chromosome expression parity between the sexes. If compensation is incomplete, this can lead to expression imbalance and sex-biased gene expression. Recent work has uncovered an intriguing and variable pattern of dosage compensation across species that includes a lack of complete dosage compensation in ZW species compared with XY species. This has led to the hypothesis that ZW species do not require complete compensation or that complete compensation would negatively affect their fitness. To date, only one study, a study of the moth Bombyx mori, has discovered evidence for complete dosage compensation in a ZW species. We examined another moth species, Manduca sexta , using high-throughput sequencing to survey gene expression in the head tissue of males and females. We found dosage compensation to be complete in M. sexta with average expression between the Z chromosome in males and females being equal. When genes expressed at very low levels are removed by filtering, we found that average autosome expression was highly similar to average Z expression, suggesting that the majority of genes in M. sexta are completely dosage compensated. Further, this compensation was accompanied by sex-specific gene expression associated with important sexually dimorphic traits. We suggest that complete dosage compensation in ZW species might be more common than previously appreciated and linked to additional selective processes, such as sexual selection. More ZW and lepidopteran species should now be examined in a phylogenetic framework, to understand the evolution of dosage compensation.

Description: While many vertebrate transcription factor (TF) families are conserved, the C2H2 zinc finger (ZNF) family stands out as a notable exception. In particular, novel ZNF gene types have arisen, duplicated, and diverged independently throughout evolution to yield many lineage-specific TF genes. This evolutionary dynamic not only raises many intriguing questions but also severely complicates identification of those ZNF genes that remain functionally conserved. To address this problem, we searched for vertebrate "DNA binding orthologs" by mining ZNF loci from eight sequenced genomes and then aligning the patterns of DNA-binding amino acids, or "fingerprints," extracted from the encoded ZNF motifs. Using this approach, we found hundreds of lineage-specific genes in each species and also hundreds of orthologous groups. Most groups of orthologs displayed some degree of fingerprint divergence between species, but 174 groups showed fingerprint patterns that have been very rigidly conserved. Focusing on the dynamic KRAB-ZNF subfamily—including nearly 400 human genes thought to possess potent KRAB-mediated epigenetic silencing activities—we found only three genes conserved between mammals and nonmammalian groups. These three genes, members of an ancient familial cluster, encode an unusual KRAB domain that functions as a transcriptional activator. Evolutionary analysis confirms the ancient provenance of this activating KRAB and reveals the independent expansion of KRAB-ZNFs in every vertebrate lineage. Most human ZNF genes, from the most deeply conserved to the primate-specific genes, are highly expressed in immune and reproductive tissues, indicating that they have been enlisted to regulate evolutionarily divergent biological traits.

Description: Toll-like receptors (TLRs) are directly involved in host–pathogen interactions. Polymorphisms in these genes are associated with susceptibility to infectious diseases. To understand the influence of environment and pathogen diversity on the evolution of TLR genes, we have undertaken a large-scale population-genetic study. Our study included two hunter–gatherer tribal populations and one urbanized nontribal population from India with distinct ethnicities ( n = 266) and 14 populations inhabiting four different continents ( n = 1,092). From the data on DNA sequences of cell-surface TLR genes, we observed an excess of rare variants and a large number of low frequency haplotypes in each gene. Nonsynonymous changes were few in every population and the commonly used statistical tests for detecting natural selection provided evidence of purifying selection. The evidence of purifying selection acting on the cell-surface TLRs of the innate immune system is not consistent with Haldane’s theory of coevolution of immunity genes, at least of innate immunity genes, with pathogens. Our study provides evidence that genes of the cell-surface TLRs, that is, TLR2 and TLR4 , have been so optimized to defend the host against microbial infections that new mutations in these genes are quickly eliminated.

Description: Nonretroviral integrated RNA viruses (NIRVs) are genes of nonretroviral RNA viruses found in the genomes of many eukaryotic organisms. NIRVs are thought to sometimes confer virus resistance, meaning that they could impact spread of the virus in the host population. However, a NIRV that is expressed may also impact the evolution of virus populations within host organisms. Here, we experimentally addressed the evolution of a virus in a host expressing a NIRV using Tobacco etch virus (TEV), a plant RNA virus, and transgenic tobacco plants expressing its replicase, NIb. We found that a virus missing the NIb gene, TEV- NIb , which is incapable of autonomous replication in wild-type plants, had a higher fitness than the full-length TEV in the transgenic plants. Moreover, when the full-length TEV was evolved by serial passages in transgenic plants, we observed genomic deletions within NIb —and in some cases the adjacent cistrons—starting from the first passage. When we passaged TEV and TEV- NIb in transgenic plants, we found mutations in proteolytic sites, but these only occurred in TEV- NIb lineages, suggesting the adaptation of polyprotein processing to altered NIb expression. These results raise the possibility that NIRV expression can indeed induce the deletion of the corresponding genes in the viral genome, resulting in the formation of viruses that are replication defective in hosts that do not express the same NIRV. Moreover, virus genome evolution was contingent upon the deletion of the viral replicase, suggesting NIRV expression could also alter patterns of virus evolution.

Description: Whole genome duplication (WGD) is widespread in flowering plants and is a driving force in angiosperm diversification. The redundancy introduced by WGD allows the evolution of novel gene interactions and functions, although the patterns and processes of diversification are poorly understood. We identified ~2,000 pairs of paralogous genes in Gossypium raimondii (cotton) resulting from an approximately 60 My old 5- to 6-fold ploidy increase. Gene expression analyses revealed that, in G. raimondii , 99.4% of the gene pairs exhibit differential expression in at least one of the three tissues (petal, leaf, and seed), with 93% to 94% exhibiting differential expression on a per-tissue basis. For 1,666 (85%) pairs, differential expression was observed in all tissues. These observations were mirrored in a time series of G. raimondii seed, and separately in leaf, petal, and seed of G. arboreum , indicating expression level diversification before species divergence. A generalized linear model revealed 92.4% of the paralog pairs exhibited expression divergence, with most exhibiting significant gene and tissue interactions indicating complementary expression patterns in different tissues. These data indicate massive, near-complete expression level neo- and/or subfunctionalization among ancient gene duplicates, suggesting these processes are essential in their maintenance over ~ 60 Ma.

Description: The histone modification H3K27me3 is involved in repression of transcription and plays a crucial role in developmental transitions in both animals and plants. It is deposited by PRC2 (Polycomb repressive complex 2), a conserved protein complex. In Arabidopsis thaliana , H3K27me3 is found at 15% of all genes. These tend to encode transcription factors and other regulators important for development. However, it is not known how PRC2 is recruited to target loci nor how this set of target genes arose during Arabidopsis evolution. To resolve the latter, we integrated A. thaliana gene families with five independent genome-wide H3K27me3 data sets. Gene families were either significantly enriched or depleted of H3K27me3, showing a strong impact of shared ancestry to H3K27me3 distribution. To quantify this, we performed ancestral state reconstruction of H3K27me3 on phylogenetic trees of gene families. The set of H3K27me3-marked genes changed less than expected by chance, suggesting that H3K27me3 was retained after gene duplication. This retention suggests that the PRC2-recruiting signal could be encoded in the DNA and also conserved among certain duplicated genes. Indeed, H3K27me3-marked genes were overrepresented among paralogs sharing conserved noncoding sequences (CNSs) that are enriched with transcription factor binding sites. The association of upstream CNSs with H3K27me3-marked genes represents the first genome-wide connection between H3K27me3 and potential regulatory elements in plants. Thus, we propose that CNSs likely function as part of the PRC2 recruitment in plants.

Description: Most land plant plastomes contain two copies of a large inverted repeat (IR) that promote high-frequency homologous recombination to generate isomeric genomic forms. Among conifer plastomes, this canonical IR is highly reduced in Pinaceae and completely lost from cupressophytes. However, both lineages have acquired short, novel IRs, some of which also exhibit recombinational activity to generate genomic structural diversity. This diversity has been shown to exist between, and occasionally within, cupressophyte species, but it is not known whether multiple genomic forms coexist within individual plants. To examine the recombinational potential of the novel cupressophyte IRs within individuals and between species, we sequenced the plastomes of four closely related species of Juniperus . The four plastomes have identical gene content and genome organization except for a large 36 kb inversion between approximately 250 bp IR containing trnQ-UUG . Southern blotting showed that different isomeric versions of the plastome predominate among individual junipers, whereas polymerase chain reaction and high-throughput read-pair mapping revealed the substoichiometric presence of the alternative isomeric form within each individual plant. Furthermore, our comparative genomic studies demonstrate that the predominant and substoichiometric arrangements of this IR have changed several times in other cupressophytes as well. These results provide compelling evidence for substoichiometric shifting of plastomic forms during cupressophyte evolution and suggest that substoichiometric shifting activity in plastid genomes may be adaptive.

Description: Chlorarachniophyte and cryptophyte algae have complex plastids that were acquired by the uptake of a green or red algal endosymbiont via secondary endosymbiosis. The plastid is surrounded by four membranes, and a relict nucleus, called the nucleomorph, remains in the periplastidal compartment that is the remnant cytoplasm of the endosymbiont. Thus, these two algae possess four different genomes in a cell: Nuclear, nucleomorph, plastid, and mitochondrial. Recently, sequencing of the nuclear genomes of the chlorarachniophyte Bigelowiella natans and the cryptophyte Guillardia theta has been completed, and all four genomes have been made available. However, the copy number of each genome has never been investigated. It is important to know the actual DNA content of each genome, especially the highly reduced nucleomorph genome, for studies on genome evolution. In this study, we calculated genomic copy numbers in B. natans and G. theta using a real-time quantitative polymerase chain reaction approach. The nuclear genomes were haploid in both species, whereas the nucleomorph genomes were estimated to be diploid and tetraploid, respectively. Mitochondria and plastids contained a large copy number of genomic DNA in each cell. In the secondary endosymbioses of chlorarachniophytes and cryptophytes, the endosymbiont nuclear genomes were highly reduced in size and in the number of coding genes, whereas the chromosomal copy number was increased, as in bacterial endosymbiont genomes. This suggests that polyploidization is a general characteristic of highly reduced genomes in broad prokaryotic and eukaryotic endosymbionts.

Description: The Arabidopsis DA1 genes appear to have multiple functions in regulating organ size and abiotic stress response, but the biological roles of its closely related genes remain unknown. Evolutionary analyses might provide some clues to aid in an understanding of their functional diversification. In this work, we characterized the molecular evolution and expressional diversification of DA1 -like genes. Surveying 354 sequenced genomes revealed 142 DA1 -like genes only in plants, indicating plant-specificity of these genes. The DA1-like protein modular structure was composed of two UIMs (ubiquitin interaction motifs), one LIM-domain (from lin-11 , isl-1 , and mec-3 ), and a conserved C-terminal, and was distinguishable from the already defined three groups of LIM-domain proteins. We further found that the DA1 -like genes diverged into Classes I and II at the ancestor of seed plants and acquired 13 clade-specific residues during their evolutionary history. Moreover, diverse intron size evolution was noted following the transition from size-expandable introns to minimal ones, accompanying the emergence and diversification of angiosperms. Functional diversification as it relates to gene expression was further investigated in soybean. Glycine max DA1 genes showed diverse tissues expression patterns during development and had substantially varied abiotic stress response expression. Thus, variations in the coding regions, intron size, and gene expression contributed to the functional diversification of this gene family. Our data suggest that the evolution of the DA1 -like genes facilitated the development of diverse molecular and functional diversification patterns to accompany the successful radiation of plants into diverse environments during evolution.

Description: Many insects harbor inherited bacterial endosymbionts. Although some of them are not strictly essential and are considered facultative, they can be a key to host survival under specific environmental conditions, such as parasitoid attacks, climate changes, or insecticide pressures. The whitefly Bemisia tabaci is at the top of the list of organisms inflicting agricultural damage and outbreaks, and changes in its distribution may be associated to global warming. In this work, we have sequenced and analyzed the genome of Cardinium cBtQ1, a facultative bacterial endosymbiont of B. tabaci and propose that it belongs to a new taxonomic family, which also includes Candidatus Amoebophilus asiaticus and Cardinium cEper1, endosymbionts of amoeba and wasps, respectively. Reconstruction of their last common ancestors’ gene contents revealed an initial massive gene loss from the free-living ancestor. This was followed in Cardinium by smaller losses, associated with settlement in arthropods. Some of these losses, affecting cofactor and amino acid biosynthetic encoding genes, took place in Cardinium cBtQ1 after its divergence from the Cardinium cEper1 lineage and were related to its settlement in the whitefly and its endosymbionts. Furthermore, the Cardinium cBtQ1 genome displays a large proportion of transposable elements, which have recently inactivated genes and produced chromosomal rearrangements. The genome also contains a chromosomal duplication and a multicopy plasmid, which harbors several genes putatively associated with gliding motility, as well as two other genes encoding proteins with potential insecticidal activity. As gene amplification is very rare in endosymbionts, an important function of these genes cannot be ruled out.

Description: Mammalian genomes are replete with millions of polymorphic sites, among which those genetic variants that are colocated on the same chromosome and exist close to one another form blocks of closely linked mutations known as haplotypes. The linkage within haplotypes is constantly disrupted due to meiotic recombination events. Whole ensembles of such numerous haplotypes are subjected to evolutionary pressure, where mutations influence each other and should be considered as a whole entity—a gigantic matrix, unique for each individual genome. This idea was implemented into a computational approach, named Genome Evolution by Matrix Algorithms (GEMA) to model genomic changes taking into account all mutations in a population. GEMA has been tested for modeling of entire human chromosomes. The program can precisely mimic real biological processes that have influence on genome evolution such as: 1) Authentic arrangements of genes and functional genomic elements, 2) frequencies of various types of mutations in different nucleotide contexts, and 3) nonrandom distribution of meiotic recombination events along chromosomes. Computer modeling with GEMA has demonstrated that the number of meiotic recombination events per gamete is among the most crucial factors influencing population fitness. In humans, these recombinations create a gamete genome consisting on an average of 48 pieces of corresponding parental chromosomes. Such highly mosaic gamete structure allows preserving fitness of population under the intense influx of novel mutations (40 per individual) even when the number of mutations with deleterious effects is up to ten times more abundant than those with beneficial effects.

Description: Teleost fish underwent whole-genome duplication around 450 Ma followed by diploidization and loss of 80–85% of the duplicated genes. To identify a deep signature of this teleost-specific whole-genome duplication (TSGD), we searched for duplicated genes that were systematically and uniquely retained in one or other of the superorders Ostariophysi and Acanthopterygii. TSGD paralogs comprised 17–21% of total gene content. Some 2.6% (510) of TSGD paralogs were present as pairs in the Ostariophysi genomes of Danio rerio (Cypriniformes) and Astyanax mexicanus (Characiformes) but not in species from four orders of Acanthopterygii (Gasterosteiformes, Gasterosteus aculeatus ; Tetraodontiformes, Tetraodon nigroviridis ; Perciformes, Oreochromis niloticus ; and Beloniformes, Oryzias latipes ) where a single copy was identified. Similarly, 1.3% (418) of total gene number represented cases where TSGD paralogs pairs were systematically retained in the Acanthopterygian but conserved as a single copy in Ostariophysi genomes. We confirmed the generality of these results by phylogenetic and synteny analysis of 40 randomly selected linage-specific paralogs (LSPs) from each superorder and completed with the transcriptomes of three additional Ostariophysi species ( Ictalurus punctatus [Siluriformes], Sinocyclocheilus species [Cypriniformes], and Piaractus mesopotamicus [Characiformes]). No chromosome bias was detected in TSGD paralog retention. Gene ontology (GO) analysis revealed significant enrichment of GO terms relative to the human GO SLIM database for "growth," "Cell differentiation," and "Embryo development" in Ostariophysi and for "Transport," "Signal Transduction," and "Vesicle mediated transport" in Acanthopterygii. The observed patterns of paralog retention are consistent with different diploidization outcomes having contributed to the evolution/diversification of each superorder.

Description: To study the male and female lineages of East Asian and European humans, we have sequenced 25 short tandem repeat markers on 453 Y-chromosomes and collected sequences of 72 complete mitochondrial genomes to construct independent phylogenetic trees for male and female lineages. The results indicate that East Asian individuals fall into two clades, one that includes East Asian individuals only and a second that contains East Asian and European individuals. Surprisingly, the European individuals did not form an independent clade, but branched within in the East Asians. We then estimated the divergence time of the root of the European clade as ~41,000 years ago. These data indicate that, contrary to traditional views, Europeans diverged from East Asians around that time. We also address the origin of the Ainu lineage in northern Japan.

Description: Capsaspora owczarzaki , a protistan symbiont of the pulmonate snail Biomphalaria glabrata , is the centre of much interest in evolutionary biology due to its close relationship to Metazoa. The whole genome sequence of this protist has revealed new insights into the ancestral genome composition of Metazoa, in particular with regard to gene families involved in the evolution of multicellularity. The draft genome revealed the presence of 23 families of transposable element, made up from DNA transposon as well as long terminal repeat (LTR) and non-LTR retrotransposon families. The phylogenetic analyses presented here show that all of the transposable elements identified in the C. owczarzaki genome have orthologous families in Metazoa, indicating that the ancestral metazoan also had a rich diversity of elements. Molecular evolutionary analyses also show that the majority of families has recently been active within the Capsaspora genome. One family now appears to be inactive and a further five families show no evidence of current transposition. Most individual element copies are evolutionarily young; however, a small proportion of inserts appear to have persisted for longer in the genome. The families present in the genome show contrasting population histories and appear to be in different stages of their life cycles. Transcriptome data have been analyzed from multiple stages in the C. owczarzaki life cycle. Expression levels vary greatly both between families and between different stages of the life cycle, suggesting an unexpectedly complex level of transposable element regulation in a single celled organism.

Description: With their complex life cycle and highly derived body plan, echinoderms are unique among bilaterians. Although early development has been intensively studied, the molecular mechanisms underlying development of the adult echinoderm and its unusual radial body plan are largely unknown. To investigate the evolution of developmental changes in gene expression underlying radial body plan development and metamorphosis, we assembled a reference transcriptome de novo and used RNA-seq to measure gene expression profiles across larval, metamorphic, and postmetamorphic life cycle phases in the sea urchin Heliocidaris erythrogramma . Our results present a high-resolution view of gene expression dynamics during the complex transition from pre- to postmetamorphic development and suggest that distinct sets of regulatory and effector proteins are used during different life history phases. These analyses provide an important foundation for more detailed analyses of the evolution of the radial adult body of echinoderms.

Description: Current hypotheses about the history of cellular life are mainly based on analyses of cultivated organisms, but these represent only a small fraction of extant biodiversity. The sequencing of new environmental lineages therefore provides an opportunity to test, revise, or reject existing ideas about the tree of life and the origin of eukaryotes. According to the textbook three domains hypothesis, the eukaryotes emerge as the sister group to a monophyletic Archaea. However, recent analyses incorporating better phylogenetic models and an improved sampling of the archaeal domain have generally supported the competing eocyte hypothesis, in which core genes of eukaryotic cells originated from within the Archaea, with important implications for eukaryogenesis. Given this trend, it was surprising that a recent analysis incorporating new genomes from uncultivated Archaea recovered a strongly supported three domains tree. Here, we show that this result was due in part to the use of a poorly fitting phylogenetic model and also to the inclusion by an automated pipeline of genes of putative bacterial origin rather than nucleocytosolic versions for some of the eukaryotes analyzed. When these issues were resolved, analyses including the new archaeal lineages placed core eukaryotic genes within the Archaea. These results are consistent with a number of recent studies in which improved archaeal sampling and better phylogenetic models agree in supporting the eocyte tree over the three domains hypothesis.

Description: Ascidians are a fascinating group of filter-feeding marine chordates characterized by rapid evolution of both sequences and structure of their nuclear and mitochondrial genomes. Moreover, they include several model organisms used to investigate complex biological processes in chordates. To study the evolutionary dynamics of ascidians at short phylogenetic distances, we sequenced 13 new mitogenomes and analyzed them, together with 15 other available mitogenomes, using a novel approach involving detailed whole-mitogenome comparisons of conspecific and congeneric pairs. The evolutionary rate was quite homogeneous at both intraspecific and congeneric level, and the lowest congeneric rates were found in cryptic (morphologically undistinguishable) and in morphologically very similar species pairs. Moreover, congeneric nonsynonymous rates (d N ) were up to two orders of magnitude higher than in intraspecies pairs. Overall, a clear-cut gap sets apart conspecific from congeneric pairs. These evolutionary peculiarities allowed easily identifying an extraordinary intraspecific variability in the model ascidian Botryllus schlosseri , where most pairs show a d N value between that observed at intraspecies and congeneric level, yet consistently lower than that of the Ciona intestinalis cryptic species pair. These data suggest ongoing speciation events producing genetically distinct B. schlosseri entities. Remarkably, these ongoing speciation events were undetectable by the cox1 barcode fragment, demonstrating that, at low phylogenetic distances, the whole mitogenome has a higher resolving power than cox1 . Our study shows that whole-mitogenome comparative analyses, performed on a suitable sample of congeneric and intraspecies pairs, may allow detecting not only cryptic species but also ongoing speciation events.

Description: Studies of domesticated animals have led to the suggestion that domestication could have significant effects on patterns of molecular evolution. In particular, analyses of mitochondrial genome sequences from domestic dogs and yaks have yielded higher ratios of non-synonymous to synonymous substitutions in the domesticated lineages than in their wild relatives. These results are important because they imply that changes to selection or population size operating over a short timescale can cause significant changes to the patterns of mitochondrial molecular evolution. In this study, our aim is to test whether the impact on mitochondrial genome evolution is a general feature of domestication or whether it is specific to particular examples. We test whether domesticated mammals and birds have consistently different patterns of molecular evolution than their wild relatives for 16 phylogenetically independent comparisons of mitochondrial genome sequences. We find no consistent difference in branch lengths or d N / d S between domesticated and wild lineages. We also find no evidence that our failure to detect a consistent pattern is due to the short timescales involved or low genetic distance between domesticated lineages and their wild relatives. However, removing comparisons where the wild relative may also have undergone a bottleneck does reveal a pattern consistent with reduced effective population size in domesticated lineages. Our results suggest that, although some domesticated lineages may have undergone changes to selective regime or effective population size that could have affected mitochondrial evolution, it is not possible to generalize these patterns over all domesticated mammals and birds.

Description: Mitochondrial ORFans (open reading frames having no detectable homology and with unknown function) were discovered in bivalve molluscs with doubly uniparental inheritance (DUI) of mitochondria. In these animals, two mitochondrial lineages are present, one transmitted through eggs (F-type), the other through sperm (M-type), each showing a specific ORFan. In this study, we used in situ hybridization and immunocytochemistry to provide evidence for the expression of Ruditapes philippinarum male-specific ORFan ( orf21 ): both the transcript and the protein (RPHM21) were localized in spermatogenic cells and mature spermatozoa; the protein was localized in sperm mitochondria and nuclei, and in early embryos. Also, in silico analyses of orf21 flanking region and RPHM21 structure supported its derivation from viral sequence endogenization. We propose that RPHM21 prevents the recognition of M-type mitochondria by the degradation machinery, allowing their survival in the zygote. The process might involve a mechanism similar to that of Modulators of Immune Recognition, viral proteins involved in the immune recognition pathway, to which RPHM21 showed structural similarities. A viral origin of RPHM21 may also support a developmental role, because some integrated viral elements are involved in development and sperm differentiation of their host. Mitochondrial ORFans could be responsible for or participate in the DUI mechanism and their viral origin could explain the acquired capability of M-type mitochondria to avoid degradation and invade the germ line, that is what viruses do best: to elude host immune system and proliferate.

Description: Because of methodological breakthroughs and the availability of an increasing amount of whole-genome sequence data, horizontal transfers (HTs) in eukaryotes have received much attention recently. Contrary to similar analyses in prokaryotes, most studies in eukaryotes usually investigate particular sequences corresponding to transposable elements (TEs), neglecting the other components of the genome. We present a new methodological framework for the genome-wide detection of all putative horizontally transferred sequences between two species that requires no prior knowledge of the transferred sequences. This method provides a broader picture of HTs in eukaryotes by fully exploiting complete-genome sequence data. In contrast to previous genome-wide approaches, we used a well-defined statistical framework to control for the number of false positives in the results, and we propose two new validation procedures to control for confounding factors. The first validation procedure relies on a comparative analysis with other species of the phylogeny to validate HTs for the nonrepeated sequences detected, whereas the second one built upon the study of the dynamics of the detected TEs. We applied our method to two closely related Drosophila species, Drosophila melanogaster and D. simulans , in which we discovered 10 new HTs in addition to all the HTs previously detected in different studies, which underscores our method’s high sensitivity and specificity. Our results favor the hypothesis of multiple independent HTs of TEs while unraveling a small portion of the network of HTs in the Drosophila phylogeny.

Description: From their origin as an early alpha proteobacterial endosymbiont to their current state as cellular organelles, large-scale genomic reorganization has taken place in the mitochondria of all main eukaryotic lineages. So far, most studies have focused on plant and animal mitochondrial (mt) genomes (mtDNA), but fungi provide new opportunities to study highly differentiated mtDNAs. Here, we analyzed 38 complete fungal mt genomes to investigate the evolution of mtDNA gene order among fungi. In particular, we looked for evidence of nonhomologous intrachromosomal recombination and investigated the dynamics of gene rearrangements. We investigated the effect that introns, intronic open reading frames (ORFs), and repeats may have on gene order. Additionally, we asked whether the distribution of transfer RNAs (tRNAs) evolves independently to that of mt protein-coding genes. We found that fungal mt genomes display remarkable variation between and within the major fungal phyla in terms of gene order, genome size, composition of intergenic regions, and presence of repeats, introns, and associated ORFs. Our results support previous evidence for the presence of mt recombination in all fungal phyla, a process conspicuously lacking in most Metazoa. Overall, the patterns of rearrangements may be explained by the combined influences of recombination (i.e., most likely nonhomologous and intrachromosomal), accumulated repeats, especially at intergenic regions, and to a lesser extent, mobile element dynamics.

Description: Polar bears ( Ursus maritimus ) face extremely cold temperatures and periods of fasting, which might result in more severe energetic challenges than those experienced by their sister species, the brown bear ( U. arctos ). We have examined the mitochondrial and nuclear genomes of polar and brown bears to investigate whether polar bears demonstrate lineage-specific signals of molecular adaptation in genes associated with cellular respiration/energy production. We observed increased evolutionary rates in the mitochondrial cytochrome c oxidase I gene in polar but not brown bears. An amino acid substitution occurred near the interaction site with a nuclear-encoded subunit of the cytochrome c oxidase complex and was predicted to lead to a functional change, although the significance of this remains unclear. The nuclear genomes of brown and polar bears demonstrate different adaptations related to cellular respiration. Analyses of the genomes of brown bears exhibited substitutions that may alter the function of proteins that regulate glucose uptake, which could be beneficial when feeding on carbohydrate-dominated diets during hyperphagia, followed by fasting during hibernation. In polar bears, genes demonstrating signatures of functional divergence and those potentially under positive selection were enriched in functions related to production of nitric oxide (NO), which can regulate energy production in several different ways. This suggests that polar bears may be able to fine-tune intracellular levels of NO as an adaptive response to control trade-offs between energy production in the form of adenosine triphosphate versus generation of heat (thermogenesis).

Description: The G-protein-coupled receptor (GPCR) signaling system is one of the main signaling pathways in eukaryotes. Here, we analyze the evolutionary history of all its components, from receptors to regulators, to gain a broad picture of its system-level evolution. Using eukaryotic genomes covering most lineages sampled to date, we find that the various components of the GPCR signaling pathway evolved independently, highlighting the modular nature of this system. Our data show that some GPCR families, G proteins, and regulators of G proteins diversified through lineage-specific diversifications and recurrent domain shuffling. Moreover, most of the gene families involved in the GPCR signaling system were already present in the last common ancestor of eukaryotes. Furthermore, we show that the unicellular ancestor of Metazoa already had most of the cytoplasmic components of the GPCR signaling system, including, remarkably, all the G protein alpha subunits, which are typical of metazoans. Thus, we show how the transition to multicellularity involved conservation of the signaling transduction machinery, as well as a burst of receptor diversification to cope with the new multicellular necessities.

Description: Chlorophylls (Chls) play pivotal roles in energy absorption and transduction and also in charge separation in reaction centers in all photosynthetic organisms. In Chl biosynthesis steps, only a step for the enzymatic reduction of protochlorophyllide (Pchlide) to chlorophyllide (Chlide) is mediated by both nuclear- and chloroplast-encoded genes in land plants. Many plants encode the genes for light-dependent Pchlide reductase (LPOR) and light-independent Pchlide reductase (DPOR) in the nucleus and chloroplast genome, respectively. During the diversification of land plants, the reduction step of Pchlide to Chlide has become solely dependent on LPOR, and the genes for DPOR have been lost from chloroplast genome. It remains unclear why DPOR persists in some land plants, how they were eliminated from chloroplast genomes during the diversification of land plants, and under what environmental conditions DPOR was required. We demonstrate that DPOR is functional in liverwort ( Marchantia polymorpha L.) and plays an important role in Chl biosynthesis. Having established a plastid transformation system in liverwort, we disrupted chl B, which encodes a subunit of DPOR in the M . polymorpha chloroplast genome. Morphological and Chl content analysis of a chl B mutant grown under different photoperiods revealed that DPOR is particularly required for Chl biosynthesis under short-day conditions. Our findings suggest that an environmental condition in the form of photoperiod is an important factor that determines the loss or retention of chloroplast-encoded genes mediating Pchlide reduction to Chlide.

Description: Previous studies have found that DNA-flanking low-complexity regions (LCRs) have an increased substitution rate. Here, the substitution rate was confirmed to increase in the vicinity of LCRs in several primate species, including humans. This effect was also found among human sequences from the 1000 Genomes Project. A strong correlation was found between average substitution rate per site and distance from the LCR, as well as the proportion of genes with gaps in the alignment at each site and distance from the LCR. Along with substitution rates, d N / d S ratios were also determined for each site, and the proportion of sites undergoing negative selection was found to have a negative relationship with distance from the LCR.

Description: The 20 canonical amino acids of the genetic code have been invariant over 3 billion years of biological evolution. Although various aminoacyl-tRNA synthetases can charge their cognate tRNAs with amino acid analogs, there has been no known displacement of any canonical amino acid from the code. Experimental departure from this universal protein alphabet comprising the canonical amino acids was first achieved in the mutants of the Bacillus subtilis QB928 strain, which after serial selection and mutagenesis led to the HR23 strain that could use 4-fluorotryptophan (4FTrp) but not canonical tryptophan (Trp) for propagation. To gain insight into this displacement of Trp from the genetic code by 4FTrp, genome sequencing was performed on LC33 (a precursor strain of HR23), HR23, and TR7 (a revertant of HR23 that regained the capacity to propagate on Trp). Compared with QB928, the negative regulator mtrB of Trp transport was found to be knocked out in LC33, HR23, and TR7, and sigma factor sigB was mutated in HR23 and TR7. Moreover, rpoBC encoding RNA polymerase subunits were mutated in three independent isolates of TR7 relative to HR23. Increased expression of sigB was also observed in HR23 and in TR7 growing under 4FTrp. These findings indicated that stabilization of the genetic code can be provided by just a small number of analog-sensitive proteins, forming an oligogenic barrier that safeguards the canonical amino acids throughout biological evolution.

Description: Studies of Hordeum vulgare subsp. spontaneum , the wild progenitor of cultivated barley, have mostly relied on materials collected decades ago and maintained since then ex situ in germplasm repositories. We analyzed spatial genetic variation in wild barley populations collected rather recently, exploring sequence variations at seven single-copy nuclear loci, and inferred the relationships among these populations and toward the genepool of the crop. The wild barley collection covers the whole natural distribution area from the Mediterranean to Middle Asia. In contrast to earlier studies, Bayesian assignment analyses revealed three population clusters, in the Levant, Turkey, and east of Turkey, respectively. Genetic diversity was exceptionally high in the Levant, while eastern populations were depleted of private alleles. Species distribution modeling based on climate parameters and extant occurrence points of the taxon inferred suitable habitat conditions during the ice-age, particularly in the Levant and Turkey. Together with the ecologically wide range of habitats, they might contribute to structured but long-term stable populations in this region and their high genetic diversity. For recently collected individuals, Bayesian assignment to geographic clusters was generally unambiguous, but materials from genebanks often showed accessions that were not placed according to their assumed geographic origin or showed traces of introgression from cultivated barley. We assign this to gene flow among accessions during ex situ maintenance. Evolutionary studies based on such materials might therefore result in wrong conclusions regarding the history of the species or the origin and mode of domestication of the crop, depending on the accessions included.

Description: Photosynthesis by diatoms accounts for roughly one-fifth of global primary production, but despite this, relatively little is known about their plastid genomes. We report the completely sequenced plastid genomes for eight phylogenetically diverse diatoms and show them to be variable in size, gene and foreign sequence content, and gene order. The genomes contain a core set of 122 protein-coding genes, with 15 additional genes exhibiting complex patterns of 1) gene losses at varying phylogenetic scales, 2) functional transfers to the nucleus, 3) gene duplication, divergence, and differential retention of paralogs, and 4) acquisitions of putatively functional recombinase genes from resident plasmids. The newly sequenced genomes also contain several previously unreported genes, highlighting how poorly characterized diatom plastid genomes are overall. Genome size variation reflects major expansions of the inverted repeat region in some cases but, more commonly, large-scale expansions of intergenic regions, many of which contain unique open reading frames of likely foreign origin. Although many gene clusters are conserved across species, rearrangements appear to be frequent in most lineages.

Description: The discovery of Chromera velia , a free-living photosynthetic relative of apicomplexan pathogens, has provided an unexpected opportunity to study the algal ancestry of malaria parasites. In this work, we compared the molecular footprints of a eukaryote-to-eukaryote endosymbiosis in C. velia to their equivalents in peridinin-containing dinoflagellates (PCD) to reevaluate recent claims in favor of a common ancestry of their plastids. To this end, we established the draft genome and a set of full-length cDNA sequences from C. velia via next-generation sequencing. We documented the presence of a single coxI gene in the mitochondrial genome, which thus represents the genetically most reduced aerobic organelle identified so far, but focused our analyses on five "lucky genes" of the Calvin cycle. These were selected because of their known support for a common origin of complex plastids from cryptophytes, alveolates (represented by PCDs), stramenopiles, and haptophytes (CASH) via a single secondary endosymbiosis with a red alga. As expected, our broadly sampled phylogenies of the nuclear-encoded Calvin cycle markers support a rhodophycean origin for the complex plastid of Chromera . However, they also suggest an independent origin of apicomplexan and dinophycean (PCD) plastids via two eukaryote-to-eukaryote endosymbioses. Although at odds with the current view of a common photosynthetic ancestry for alveolates, this conclusion is nonetheless in line with the deviant plastome architecture in dinoflagellates and the morphological paradox of four versus three plastid membranes in the respective lineages. Further support for independent endosymbioses is provided by analysis of five additional markers, four of them involved in the plastid protein import machinery. Finally, we introduce the "rhodoplex hypothesis" as a convenient way to designate evolutionary scenarios where CASH plastids are ultimately the product of a single secondary endosymbiosis with a red alga but were subsequently horizontally spread via higher-order eukaryote-to-eukaryote endosymbioses.

Description: Several bacterial families are known to be highly abundant within the human microbiome, but their ecological roles and evolutionary histories have yet to be investigated in depth. One such family, Lachnospiraceae (phylum Firmicutes, class Clostridia) is abundant in the digestive tracts of many mammals and relatively rare elsewhere. Members of this family have been linked to obesity and protection from colon cancer in humans, mainly due to the association of many species within the group with the production of butyric acid, a substance that is important for both microbial and host epithelial cell growth. We examined the genomes of 30 Lachnospiraceae isolates to better understand the origin of butyric acid capabilities and other ecological adaptations within this group. Butyric acid production-related genes were detected in fewer than half of the examined genomes with the distribution of this function likely arising in part from lateral gene transfer (LGT). An investigation of environment-specific functional signatures indicated that human gut-associated Lachnospiraceae possess genes for endospore formation, whereas other members of this family lack key sporulation-associated genes, an observation supported by analysis of metagenomes from the human gut, oral cavity, and bovine rumen. Our analysis demonstrates that adaptation to an ecological niche and acquisition of defining functional roles within a microbiome can arise through a combination of both habitat-specific gene loss and LGT.

Description: Scale insects (Hemiptera: Coccoidae) constitute a very diverse group of sap-feeding insects with a large diversity of symbiotic associations with bacteria. Here, we present the complete genome sequence, metabolic reconstruction, and comparative genomics of the flavobacterial endosymbiont of the giant scale insect Llaveia axin axin. The gene repertoire of its 309,299 bp genome was similar to that of other flavobacterial insect endosymbionts though not syntenic. According to its genetic content, essential amino acid biosynthesis is likely to be the flavobacterial endosymbiont's principal contribution to the symbiotic association with its insect host. We also report the presence of a -proteobacterial symbiont that may be involved in waste nitrogen recycling and also has amino acid biosynthetic capabilities that may provide metabolic precursors to the flavobacterial endosymbiont. We propose " Candidatus Walczuchella monophlebidarum" as the name of the flavobacterial endosymbiont of insects from the Monophlebidae family.

Description: Insertion sequences (ISs) are small transposable elements widespread in bacterial genomes, where they play an essential role in chromosome evolution by stimulating recombination and genetic flow. Despite their ubiquity, it is unclear how ISs interact with the host. Here, we report a survey of the orientation patterns of ISs in bacterial chromosomes with the objective of gaining insight into the interplay between ISs and host chromosomal functions. We find that a significant fraction of IS families present a consistent and family-specific orientation bias with respect to chromosomal DNA replication, especially in Firmicutes. Additionally, we find that the transposases of up to nine different IS families with different transposition pathways interact with the β sliding clamp, an essential replication factor, suggesting that this is a widespread mechanism of interaction with the host. Although we find evidence that the interaction with the β sliding clamp is common to all bacterial phyla, it also could explain the observed strong orientation bias found in Firmicutes, because in this group β is asymmetrically distributed during synthesis of the leading or lagging strands. Besides the interaction with the β sliding clamp, other asymmetries also play a role in the biased orientation of some IS families. The utilization of the highly conserved replication sliding clamps suggests a mechanism for host regulation of IS proliferation and also a universal platform for IS dispersal and transmission within bacterial populations and among phylogenetically distant species.

Description: The genus Streptococcus comprises important pathogens that have a severe impact on human health and are responsible for substantial economic losses to agriculture. Here, we utilize 46 Streptococcus genome sequences (44 species), including eight species sequenced here, to provide the first genomic level insight into the evolutionary history and genetic basis underlying the functional diversity of all major groups of this genus. Gene gain/loss analysis revealed a dynamic pattern of genome evolution characterized by an initial period of gene gain followed by a period of loss, as the major groups within the genus diversified. This was followed by a period of genome expansion associated with the origins of the present extant species. The pattern is concordant with an emerging view that genomes evolve through a dynamic process of expansion and streamlining. A large proportion of the pan-genome has experienced lateral gene transfer (LGT) with causative factors, such as relatedness and shared environment, operating over different evolutionary scales. Multiple gene ontology terms were significantly enriched for each group, and mapping terms onto the phylogeny showed that those corresponding to genes born on branches leading to the major groups represented approximately one-fifth of those enriched. Furthermore, despite the extensive LGT, several biochemical characteristics have been retained since group formation, suggesting genomic cohesiveness through time, and that these characteristics may be fundamental to each group. For example, proteolysis: mitis group; urea metabolism: salivarius group; carbohydrate metabolism: pyogenic group; and transcription regulation: bovis group.

Description: Recent studies have demonstrated a high genomic plasticity in Acinetobacter baumannii , which may explain its high capacity to acquire multiple antibiotic resistance determinants and to survive in the hospital environment. Acinetobacter baumannii strain A118 ( Ab A118) was isolated in the year 1995 from a blood culture of an intensive care unit patient. As this particular strain showed some peculiar characteristic such as being naturally competent and susceptible to numerous antibiotics, we performed whole-genome comparison (WGC) studies to gain insights into the nature and extent of the genomic differences. The Ab A118 genome is approximately 3,824 kb long with a 38.4% GC content and contains 3,520 coding sequences. WGC studies showed that the Ab A118 genome has 98% average nucleotide identity with that of A. baumannii ATCC 17978, and 96% average nucleotide identity with that of strains AYE and ACICU. At least 12 inversions, 275 insertions, and 626 deletions were identified when the Ab A118 genome was compared with those of strains ATCC 17978, AYE, and ACICU using MAUVE WGC. Multiple gene order arrangements were observed among the analyzed strains. MAUVE WGC analysis identified 19 conserved segments, known as locally colinear blocks. The number of single nucleotide polymorphisms found when comparing the Ab A118 genome with that of strains ATCC 17978, AYE, and ACICU was 43,784 (1.1496%), 44,130 (1.158%), and 43,914 (1.153%), respectively. Genes comEA, pilQ, pilD, pilF, comL, pilA, comEC, pilI, pilH, pilO, pilN, pilY1(comC), pilE, pilR, and comM , potentially involved in natural competence were found in the Ab A118 genome. In particular, unlike in most strains where comM is interrupted by an insertion of a resistance island (AbaR), in strain Ab A118 it is uninterrupted.

Description: Discoba (Excavata) is an ancient group of eukaryotes with great morphological and ecological diversity. Unlike the other major divisions of Discoba (Jakobida and Euglenozoa), little is known about the mitochondrial DNAs (mtDNAs) of Heterolobosea. We have assembled a complete mtDNA genome from the aggregating heterolobosean amoeba, Acrasis kona , which consists of a single circular highly AT-rich (83.3%) molecule of 51.5 kb. Unexpectedly, A. kona mtDNA is missing roughly 40% of the protein-coding genes and nearly half of the transfer RNAs found in the only other sequenced heterolobosean mtDNAs, those of Naegleria spp. Instead, over a quarter of A. kona mtDNA consists of novel open reading frames. Eleven of the 16 protein-coding genes missing from A. kona mtDNA were identified in its nuclear DNA and polyA RNA, and phylogenetic analyses indicate that at least 10 of these 11 putative nuclear-encoded mitochondrial (NcMt) proteins arose by direct transfer from the mitochondrion. Acrasis kona mtDNA also employs C-to-U type RNA editing, and 12 homologs of DYW-type pentatricopeptide repeat (PPR) proteins implicated in plant organellar RNA editing are found in A. kona nuclear DNA. A mapping of mitochondrial gene content onto a consensus phylogeny reveals a sporadic pattern of relative stasis and rampant gene loss in Discoba. Rampant loss occurred independently in the unique common lineage leading to Heterolobosea + Tsukubamonadida and later in the unique lineage leading to Acrasis. Meanwhile, mtDNA gene content appears to be remarkably stable in the Acrasis sister lineage leading to Naegleria and in their distant relatives Jakobida .

Description: Feathers are hallmark avian integument appendages, although they were also present on theropods. They are composed of flexible corneous materials made of α- and β-keratins, but their genomic organization and their functional roles in feathers have not been well studied. First, we made an exhaustive search of α- and β-keratin genes in the new chicken genome assembly (Galgal4). Then, using transcriptomic analysis, we studied α- and β-keratin gene expression patterns in five types of feather epidermis. The expression patterns of β-keratin genes were different in different feather types, whereas those of α-keratin genes were less variable. In addition, we obtained extensive α- and β-keratin mRNA in situ hybridization data, showing that α-keratins and β-keratins are preferentially expressed in different parts of the feather components. Together, our data suggest that feather morphological and structural diversity can largely be attributed to differential combinations of α- and β-keratin genes in different intrafeather regions and/or feather types from different body parts. The expression profiles provide new insights into the evolutionary origin and diversification of feathers. Finally, functional analysis using mutant chicken keratin forms based on those found in the human α-keratin mutation database led to abnormal phenotypes. This demonstrates that the chicken can be a convenient model for studying the molecular biology of human keratin-based diseases.

Description: NAD + is an essential molecule for life, present in each living cell. It can function as an electron carrier or cofactor in redox biochemistry and energetics, and serves as substrate to generate the secondary messenger cyclic ADP ribose and nicotinic acid adenine dinucleotide phosphate. Although de novo NAD + biosynthesis is essential, different metabolic pathways exist in different eukaryotic clades. The kynurenine pathway starting with tryptophan was most likely present in the last common ancestor of all eukaryotes, and is active in fungi and animals. The aspartate pathway, detected in most photosynthetic eukaryotes, was probably acquired from the cyanobacterial endosymbiont that gave rise to chloroplasts. An evolutionary analysis of enzymes catalyzing de novo NAD + biosynthesis resulted in evolutionary trees incongruent with established organismal phylogeny, indicating numerous gene transfers. Endosymbiotic gene transfers probably introduced the aspartate pathway into eukaryotes and may have distributed it among different photosynthetic clades. In addition, several horizontal gene transfers substituted eukaryotic genes with bacterial orthologs. Although horizontal gene transfer is accepted as a key mechanism in prokaryotic evolution, it is supposed to be rare in eukaryotic evolution. The essential metabolic pathway of de novo NAD + biosynthesis in eukaryotes was shaped by numerous gene transfers.

Description: Many studies on conserved noncoding sequences (CNSs) have found that CNSs are enriched significantly in regulatory sequence elements. We conducted whole-genome analysis on plant CNSs to identify lineage-specific CNSs in eudicots, monocots, angiosperms, and vascular plants based on the premise that lineage-specific CNSs define lineage-specific characters and functions in groups of organisms. We identified 27 eudicot, 204 monocot, 6,536 grass, 19 angiosperm, and 2 vascular plant lineage-specific CNSs (lengths range from 16 to 1,517 bp) that presumably originated in their respective common ancestors. A stronger constraint on the CNSs located in the untranslated regions was observed. The CNSs were often flanked by genes involved in transcription regulation. A drop of A+T content near the border of CNSs was observed and CNS regions showed a higher nucleosome occupancy probability. These CNSs are candidate regulatory elements, which are expected to define lineage-specific features of various plant groups.

Description: Males and females of the same species differ largely in gene expression, which accounts for most of the morphological and physiological differences and sex-specific phenotypes. Here, we analyzed sex-specific gene expression in the brain and the gonads of cichlid fishes from Lake Tanganyika belonging to four different lineages, so-called tribes (Eretmodini, Ectodini, Haplochromini, and Lamprologini), using the outgroup Nile tilapia ( Oreochromis niloticus ) as reference. The comparison between male and female brains revealed few differences between the sexes, consistent in all investigated species. The gonads, on the other hand, showed a large fraction of differentially expressed transcripts with the majority of them showing the same direction of expression in all four species. All here-studied cichlids, especially the three investigated mouth-breeding species, showed a trend toward more male- than female-biased transcripts. Transcripts, which were female-biased in expression in all four species, were overrepresented on linkage group (LG)1 in the reference genome and common male-biased transcripts showed accumulation on LG23, the presumable sex chromosomes of the Nile tilapia. Sex-specific transcripts contained candidate genes for sex determination and differentiation in fishes, especially members of the transforming growth factor-β-superfamily and the Wnt-pathway and also prominent members of the sox-, dm-domain-, and high mobility group-box families. We further confirmed our previous finding on species/lineage-specific gene expression shifts in the sex steroid pathway, including synthesizing enzymes as the aromatase cyp19a1 and estrogen and androgen receptors.

Description: Lignocellulosic plant material is a viable source of biomass to produce alternative energy including ethanol and other biofuels. However, several factors—including toxic byproducts from biomass pretreatment and poor fermentation of xylose and other pentose sugars—currently limit the efficiency of microbial biofuel production. To begin to understand the genetic basis of desirable traits, we characterized three strains of Saccharomyces cerevisiae with robust growth in a pretreated lignocellulosic hydrolysate or tolerance to stress conditions relevant to industrial biofuel production, through genome and transcriptome sequencing analysis. All stress resistant strains were highly mosaic, suggesting that genetic admixture may contribute to novel allele combinations underlying these phenotypes. Strain-specific gene sets not found in the lab strain were functionally linked to the tolerances of particular strains. Furthermore, genes with signatures of evolutionary selection were enriched for functional categories important for stress resistance and included stress-responsive signaling factors. Comparison of the strains’ transcriptomic responses to heat and ethanol treatment—two stresses relevant to industrial bioethanol production—pointed to physiological processes that were related to particular stress resistance profiles. Many of the genotype-by-environment expression responses occurred at targets of transcription factors with signatures of positive selection, suggesting that these strains have undergone positive selection for stress tolerance. Our results generate new insights into potential mechanisms of tolerance to stresses relevant to biofuel production, including ethanol and heat, present a backdrop for further engineering, and provide glimpses into the natural variation of stress tolerance in wild yeast strains.

Description: The family Actinomycetaceae comprises several important pathogens that impose serious threat to human health and cause substantial infections of economically important animals. However, the phylogeny and evolutionary dynamic of this family are poorly characterized. Here, we provide detailed description of the genome characteristics of Trueperella pyogenes , a prevalent opportunistic bacterium that belongs to the family Actinomycetaceae , and the results of comparative genomics analyses suggested that T. pyogenes was a more versatile pathogen than Arcanobacterium haemolyticum in adapting various environments. We then performed phylogenetic analyses at the genomic level and showed that, on the whole, the established members of the family Actinomycetaceae were clearly separated with high bootstrap values but confused with the dominant genus Actinomyces , because the species of genus Actinomyces were divided into three main groups with different G+C content. Although T. pyogenes and A. haemolyticum were found to share the same branch as previously determined, our results of single nucleotide polymorphism tree and genome clustering as well as predicted intercellular metabolic analyses provide evidence that they are phylogenetic neighbors. Finally, we found that the gene gain/loss events occurring in each species may play an important role during the evolution of Actinomycetaceae from free-living to a specific lifestyle.

Description: Introns in protein-coding genes are very rare in hemiascomycetous yeast genomes. It has been suggested that these species have experienced extensive intron loss during their evolution from the postulated intron-rich fungal ancestor. However, no intron-devoid yeast species have been identified and some of the introns remaining within the genomes of intron-poor species, such as Saccharomyces cerevisiae , appear to be beneficial during growth under stress conditions. In order to reveal the pattern of intron retention within intron-poor yeast species and better understand the mechanisms of intron evolution, we generated a comprehensive set of 250 orthologous introns in the 20 species that comprise the Saccharomycetaceae, by analyzing RNA deep-sequencing data and alignments of intron-containing genes. Analysis of these intron sets shows that intron loss is at least two orders of magnitude more frequent than intron gain. Fine mapping of intron positions shows that intron sliding is rare, and that introns are almost always removed without changing the primary sequence of the encoded protein. The latter finding is consistent with the prevailing view that homologous recombination between reverse-transcribed mature mRNAs and the corresponding genomic locus is the primary mechanism of intron loss. However, we also find evidence that loss of a small number of introns is mediated by micro-homology, and that the number of intron losses is diminished in yeast species that have lost the microhomology end joining and nonhomologous end joining machinery.

Description: The increasing availability of mitochondrial (mt) sequence data from various yeasts provides a tool to study genomic evolution within and between different species. While the genomes from a range of lineages are available, there is a lack of information concerning intraspecific mtDNA diversity. Here, we analyzed the mt genomes of 50 strains from Lachancea thermotolerans , a protoploid yeast species that has been isolated from several locations (Europe, Asia, Australia, South Africa, and North / South America) and ecological sources (fruit, tree exudate, plant material, and grape and agave fermentations). Protein-coding genes from the mtDNA were used to construct a phylogeny, which reflected a similar, yet less resolved topology than the phylogenetic tree of 50 nuclear genes. In comparison to its sister species Lachancea kluyveri , L. thermotolerans has a smaller mt genome. This is due to shorter intergenic regions and fewer introns, of which the latter are only found in COX1 . We revealed that L. kluyveri and L. thermotolerans share similar levels of intraspecific divergence concerning the nuclear genomes. However, L. thermotolerans has a more highly conserved mt genome with the coding regions characterized by low rates of nonsynonymous substitution. Thus, in the mt genomes of L. thermotolerans, stronger purifying selection and lower mutation rates potentially shape genome diversity in contract to what was found for L. kluyveri , demonstrating that the factors driving mt genome evolution are different even between closely related species.

Description: Rolling-circle transposons ( Helitrons ) are a newly discovered group of mobile DNA widespread in plant and invertebrate genomes but limited to the bat family Vespertilionidae among mammals. Little is known about the long-term impact of Helitron activity because the genomes where Helitron activity has been extensively studied are predominated by young families. Here, we report a comprehensive catalog of vetted Helitrons from the 7 x Myotis lucifugus genome assembly. To estimate the timing of transposition, we scored presence/absence across related vespertilionid genome sequences with estimated divergence times. This analysis revealed that the Helibat family has been a persistent source of genomic innovation throughout the vespertilionid diversification from approximately 30–36 Ma to as recently as approximately 1.8–6 Ma. This is the first report of persistent Helitron transposition over an extended evolutionary timeframe. These findings illustrate that the pattern of Helitron activity is akin to the vertical persistence of LINE retrotransposons in primates and other mammalian lineages. Like retrotransposition in primates, rolling-circle transposition has generated lineage-specific variation and accounts for approximately 110 Mb, approximately 6% of the genome of M. lucifugus . The Helitrons carry a heterogeneous assortment of host sequence including retroposed messenger RNAs, retrotransposons, DNA transposons, as well as introns, exons and regulatory regions (promoters, 5'-untranslated regions [UTRs], and 3'-UTRs) of which some are evolving in a pattern suggestive of purifying selection. Evidence that Helitrons have contributed putative promoters, exons, splice sites, polyadenylation sites, and microRNA-binding sites to transcripts otherwise conserved across mammals is presented, and the implication of Helitron activity to innovation in these unique mammals is discussed.

Description: Gene duplication leads to paralogy, which complicates the de novo assembly of genotyping-by-sequencing (GBS) data. The issue of paralogous genes is exacerbated in plants, because they are particularly prone to gene duplication events. Paralogs are normally filtered from GBS data before undertaking population genomics or phylogenetic analyses. However, gene duplication plays an important role in the functional diversification of genes and it can also lead to the formation of postzygotic barriers. Using populations and closely related species of a tropical mountain shrub, we examine 1) the genomic differentiation produced by putative orthologs, and 2) the distribution of recent gene duplication among lineages and geography. We find high differentiation among populations from isolated mountain peaks and species-level differentiation within what is morphologically described as a single species. The inferred distribution of paralogs among populations is congruent with taxonomy and shows that GBS could be used to examine recent gene duplication as a source of genomic differentiation of nonmodel species.

Description: A significant limitation when testing the putative single origin of primary plastids and the monophyly of the Archaeplastida supergroup, comprised of the red algae, viridiplants, and glaucophytes, is the scarce nuclear and organellar genome data available from the latter lineage. The Glaucophyta are a key algal group when investigating the origin and early diversification of photosynthetic eukaryotes. However, so far only the plastid and mitochondrial genomes of the glaucophytes Cyanophora paradoxa (strain CCMP 329) and Glaucocystis nostochinearum (strain UTEX 64) have been completely sequenced. Here, we present the complete mitochondrial genomes of Gloeochaete wittrockiana SAG 46.84 (36.05 kb; 33 protein-coding genes, 6 unidentified open reading frames [ORFs], and 28 transfer RNAs [tRNAs]) and Cyanoptyche gloeocystis SAG 4.97 (33.24 kb; 33 protein-coding genes, 6 unidentified ORFs, and 26 tRNAs), which represent two genera distantly related to the "well-known" Cyanophora and Glaucocystis. The mitochondrial gene repertoire of the four glaucophyte species is highly conserved, whereas the gene order shows considerable variation. Phylogenetic analyses of 14 mitochondrial genes from representative taxa from the major eukaryotic supergroups, here including novel sequences from the glaucophytes Cyanophora tetracyanea (strain NIES-764) and Cyanophora biloba (strain UTEX LB 2766), recover a clade uniting the three Archaeplastida lineages; this recovery is dependent on our novel glaucophyte data, demonstrating the importance of greater taxon sampling within the glaucophytes.

Description: Many mosquito species serve as vectors of diseases such as malaria and yellow fever, wherein pathogen transmission is tightly associated with the reproductive requirement of taking vertebrate blood meals. Toxorhynchites is one of only three known mosquito genera that does not host-seek and initiates egg development in the absence of a blood-derived protein bolus. These remarkable differences make Toxorhynchites an attractive comparative reference for understanding mosquito chemosensation as it pertains to host-seeking. We performed deep transcriptome profiling of adult female Toxorhynchites amboinensis bodies, antennae and maxillary palps, and identified 25,084 protein-coding "genes" in the de novo assembly. Phylogenomic analysis of 4,266 single-copy "genes" from T. amboinensis , Aedes aegypti , Anopheles gambiae, and Culex quinquefasciatus robustly supported Ae. aegypti as the closest relative of T. amboinensis, with the two species diverged approximately 40 Ma. We identified a large number of T. amboinensis chemosensory "genes," the majority of which have orthologs in other mosquitoes. Finally, cross-species expression analyses indicated that patterns of chemoreceptor transcript abundance were very similar for chemoreceptors that are conserved between T. amboinensis and Ae. aegypti , whereas T. amboinensis appeared deficient in the variety of expressed, lineage-specific chemoreceptors. Our transcriptome assembly of T. amboinensis represents the first comprehensive genomic resource for a nonblood-feeding mosquito and establishes a foundation for future comparative studies of blood-feeding and nonblood-feeding mosquitoes. We hypothesize that chemosensory genes that display discrete patterns of evolution and abundance between T. amboinensis and blood-feeding mosquitoes are likely to play critical roles in host-seeking and hence the vectorial capacity.

Description: Recombination clusters nonuniformly across mammalian genomes at discrete genomic loci referred to as recombination hotspots. Despite their ubiquitous presence, individual hotspots rapidly lose their activities, and the molecular and evolutionary mechanisms underlying such frequent hotspot turnovers (the so-called "recombination hotspot paradox") remain unresolved. Even though some sequence motifs are significantly associated with hotspots, multiple lines of evidence indicate that factors other than underlying sequences, such as epigenetic modifications, may affect the evolution of recombination hotspots. Thus, identifying epigenetic factors that covary with recombination at fine-scale is a promising step for this important research area. It was previously reported that recombination rates correlate with indirect measures of DNA methylation in the human genome. Here, we analyze experimentally determined DNA methylation and histone modification of human sperms, and show that the correlation between DNA methylation and recombination in long-range windows does not hold with respect to the spatial and temporal variation of recombination at hotspots. On the other hand, two histone modifications (H3K4me3 and H3K27me3) overlap extensively with recombination hotspots. Similar trends were observed in mice. These results indicate that specific histone modifications rather than DNA methylation are associated with the rapid evolution of recombination hotspots. Furthermore, many human recombination hotspots occupy "bivalent" chromatin regions that harbor both active (H3K4me3) and repressive (H3K27me3) marks. This may explain why human recombination hotspots tend to occur in nongenic regions, in contrast to yeast and Arabidopsis hotspots that are characterized by generally active chromatins. Our results highlight the dynamic epigenetic underpinnings of recombination hotspot evolution.

Description: Under certain circumstances, X-linked loci are expected to experience more adaptive substitutions than similar autosomal loci. To look for evidence of faster-X evolution, we analyzed the evolutionary rates of coding sequences in two sets of Drosophila species, the melanogaster and pseudoobscura clades, using whole-genome sequences. One of these, the pseudoobscura clade, contains a centric fusion between the ancestral X chromosome and the autosomal arm homologous to 3L in D. melanogaster . This offers an opportunity to study the same loci in both an X-linked and an autosomal context, and to compare these loci with those that are only X-linked or only autosomal. We therefore investigated these clades for evidence of faster-X evolution with respect to nonsynonymous substitutions, finding mixed results. Overall, there was consistent evidence for a faster-X effect in the melanogaster clade, but not in the pseudoobscura clade, except for the comparison between D. pseudoobscura and its close relative, Drosophila persimilis . An analysis of polymorphism data on a set of genes from D. pseudoobscura that evolve rapidly with respect to their protein sequences revealed no evidence for a faster-X effect with respect to adaptive protein sequence evolution; their rapid evolution is instead largely attributable to lower selective constraints. Faster-X evolution in the melanogaster clade was not related to male-biased gene expression; surprisingly, however, female-biased genes showed evidence for faster-X effects, perhaps due to their sexually antagonistic effects in males.

Description: New genes, which provide material for evolutionary innovation, have been extensively studied for many years in animals where it is observed that they commonly show an expression bias for the testis. Thus, the testis is a major source for the generation of new genes in animals. The source tissue for new genes in plants is unclear. Here, we find that new genes in plants show a bias in expression to mature pollen, and are also enriched in a gene coexpression module that correlates with mature pollen in Arabidopsis thaliana . Transposable elements are significantly enriched in the new genes, and the high activity of transposable elements in the vegetative nucleus, compared with the germ cells, suggests that new genes are most easily generated in the vegetative nucleus in the mature pollen. We propose an "out of pollen" hypothesis for the origin of new genes in flowering plants.

Description: The tricarboxylic acid (TCA) cycle, a crucial component of respiratory metabolism, is composed of a set of eight enzymes present in the mitochondrial matrix. However, most of the TCA cycle enzymes are encoded in the nucleus in higher eukaryotes. In addition, evidence has accumulated demonstrating that nuclear genes were acquired from the mitochondrial genome during the course of evolution. For this reason, we here analyzed the evolutionary history of all TCA cycle enzymes in attempt to better understand the origin of these nuclear-encoded proteins. Our results indicate that prior to endosymbiotic events the TCA cycle seemed to operate only as isolated steps in both the host (eubacterial cell) and mitochondria (alphaproteobacteria). The origin of isoforms present in different cell compartments might be associated either with gene-transfer events which did not result in proper targeting of the protein to mitochondrion or with duplication events. Further in silico analyses allow us to suggest new insights into the possible roles of TCA cycle enzymes in different tissues. Finally, we performed coexpression analysis using mitochondrial TCA cycle genes revealing close connections among these genes most likely related to the higher efficiency of oxidative phosphorylation in this specialized organelle. Moreover, these analyses allowed us to identify further candidate genes which might be used for metabolic engineering purposes given the importance of the TCA cycle during development and/or stress situations.

Description: Bacterial genomics has greatly expanded our understanding of microdiversification patterns within a species, but analyses at higher taxonomical levels are necessary to understand and predict the independent rise of pathogens in a genus. We have sampled, sequenced, and assessed the diversity of genomes of validly named and tentative species of the Acinetobacter genus, a clade including major nosocomial pathogens and biotechnologically important species. We inferred a robust global phylogeny and delimited several new putative species. The genus is very ancient and extremely diverse: Genomes of highly divergent species share more orthologs than certain strains within a species. We systematically characterized elements and mechanisms driving genome diversification, such as conjugative elements, insertion sequences, and natural transformation. We found many error-prone polymerases that may play a role in resistance to toxins, antibiotics, and in the generation of genetic variation. Surprisingly, temperate phages, poorly studied in Acinetobacter , were found to account for a significant fraction of most genomes. Accordingly, many genomes encode clustered regularly interspaced short palindromic repeats (CRISPR)-Cas systems with some of the largest CRISPR-arrays found so far in bacteria. Integrons are strongly overrepresented in Acinetobacter baumannii , which correlates with its frequent resistance to antibiotics. Our data suggest that A. baumannii arose from an ancient population bottleneck followed by population expansion under strong purifying selection. The outstanding diversification of the species occurred largely by horizontal transfer, including some allelic recombination, at specific hotspots preferentially located close to the replication terminus. Our work sets a quantitative basis to understand the diversification of Acinetobacter into emerging resistant and versatile pathogens.

Description: The extracellular matrix (ECM) is a defining characteristic of metazoans and consists of a meshwork of self-assembling, fibrous proteins, and their functionally related neighbours. Previous studies, focusing on a limited number of gene families, suggest that vertebrate complexity predominantly arose through the duplication and subsequent modification of retained, preexisting ECM genes. These genes provided the structural underpinnings to support a variety of specialized tissues, as well as a platform for the organization of spatio-temporal signaling and cell migration. However, the relative contributions of ancient versus novel domains to ECM evolution have not been quantified across the full range of ECM proteins. Here, utilizing a high quality list comprising 324 ECM genes, we reveal general and clade-specific domain combinations, identifying domains of eukaryotic and metazoan origin recruited into new roles in approximately two-third of the ECM proteins in humans representing novel vertebrate proteins. We show that, rather than acquiring new domains, sampling of new domain combinations has been key to the innovation of paralogous ECM genes during vertebrate evolution. Applying a novel framework for identifying potentially important, noncontiguous, conserved arrangements of domains, we find that the distinct biological characteristics of the ECM have arisen through unique evolutionary processes. These include the preferential recruitment of novel domains to existing architectures and the utilization of high promiscuity domains in organizing the ECM network around a connected array of structural hubs. Our focus on ECM proteins reveals that distinct types of proteins and/or the biological systems in which they operate have influenced the types of evolutionary forces that drive protein innovation. This emphasizes the need for rigorously defined systems to address questions of evolution that focus on specific systems of interacting proteins.

Description: With the goal to contribute for the understanding of satellite DNA evolution and its genomic involvement, in this work it was isolated and characterized the first satellite DNA (PSUcentSat) from Phodopus sungorus (Cricetidae). Physical mapping of this sequence in P. sungorus showed large PSUcentSat arrays located at the heterochromatic (peri)centromeric region of five autosomal pairs and Y-chromosome. The presence of orthologous PSUcentSat sequences in the genomes of other Cricetidae and Muridae rodents was also verified, presenting however, an interspersed chromosomal distribution. This distribution pattern suggests a PSUcentSat-scattered location in an ancestor of Muridae/Cricetidae families, that assumed afterwards, in the descendant genome of P. sungorus a restricted localization to few chromosomes in the (peri)centromeric region. We believe that after the divergence of the studied species, PSUcentSat was most probably highly amplified in the (peri)centromeric region of some chromosome pairs of this hamster by recombinational mechanisms. The bouquet chromosome configuration (prophase I) possibly displays an important role in this selective amplification, providing physical proximity of centromeric regions between chromosomes with similar size and/or morphology. This seems particularly evident for the acrocentric chromosomes of P. sungorus (including the Y-chromosome), all presenting large PSUcentSat arrays at the (peri)centromeric region. The conservation of this sequence in the studied genomes and its (peri)centromeric amplification in P. sungorus strongly suggests functional significance, possibly displaying this satellite family different functions in the different genomes. The verification of PSUcentSat transcriptional activity in normal proliferative cells suggests that its transcription is not stage-limited, as described for some other satellites.

Description: Monotremes (platypus and echidna) are the descendants of the oldest ancestor of all extant mammals distinguished from other mammals by mode of reproduction. Monotremes lay eggs following a short gestation period and after an even briefer incubation period, altricial hatchlings are nourished over a long lactation period with milk secreted by nipple-less mammary patches located on the female’s abdomen. Milk is the sole source of nutrition and immune protection for the developing young until weaning. Using transcriptome and mass spectrometry analysis of milk cells and milk proteins, respectively, a novel Monotreme Lactation Protein ( MLP ) was identified as a major secreted protein in milk. We show that platypus and short-beaked echidna MLP genes show significant homology and are unique to monotremes. The MLP transcript was shown to be expressed in a variety of tissues; however, highest expression was observed in milk cells and was expressed constitutively from early to late lactation. Analysis of recombinant MLP showed that it is an N-linked glycosylated protein and biophysical studies predicted that MLP is an amphipathic, α-helical protein, a typical feature of antimicrobial proteins. Functional analysis revealed MLP antibacterial activity against both opportunistic pathogenic Staphylococcus aureus and commensal Enterococcus faecalis bacteria but showed no effect on Escherichia coli, Pseudomonas aeruginosa, Staphylococcus epidermidis, and Salmonella enterica . Our data suggest that MLP is an evolutionarily ancient component of milk-mediated innate immunity absent in other mammals. We propose that MLP evolved specifically in the monotreme lineage supporting the evolution of lactation in these species to provide bacterial protection, at a time when mammals lacked nipples.

Description: The Trypanosoma brucei complex contains a number of subspecies with exceptionally variable life histories, including zoonotic subspecies, which are causative agents of human African trypanosomiasis (HAT) in sub-Saharan Africa. Paradoxically, genomic variation between taxa is extremely low. We analyzed the whole-genome sequences of 39 isolates across the T. brucei complex from diverse hosts and regions, identifying 608,501 single nucleotide polymorphisms that represent 2.33% of the nuclear genome. We show that human pathogenicity occurs across a wide range of parasite genotypes, and taxonomic designation does not reflect genetic variation across the group, as previous studies have suggested based on a small number of genes. This genome-wide study allowed the identification of significant host and geographic location associations. Strong purifying selection was detected in genomic regions associated with cytoskeleton structure, and regulatory genes associated with antigenic variation, suggesting conservation of these regions in African trypanosomes. In agreement with expectations drawn from meiotic reciprocal recombination, differences in average linkage disequilibrium between chromosomes in T. brucei correlate positively with chromosome size. In addition to insights into the life history of a diverse group of eukaryotic parasites, the documentation of genomic variation across the T. brucei complex and its association with specific hosts and geographic localities will aid in the development of comprehensive monitoring tools crucial to the proposed elimination of HAT by 2020, and on a shorter term, for monitoring the feared merger between the two human infective parasites, T. brucei rhodesiense and T. b. gambiense, in northern Uganda.

Description: Transposable elements (TEs) are one of the most important features of genome architecture, so their evolution and relationship with host defense mechanisms have been topics of intense study, especially in model systems such as Drosophila melanogaster . Recently, a novel small RNA-based defense mechanism in animals called the Piwi-interacting RNA (piRNA) pathway was discovered to form an adaptive defense mechanism against TEs. To investigate the relationship between piRNA and TE content between strains of a species, we sequenced piRNAs from 16 inbred lines of D. melanogaster from the Drosophila Genetic Reference Panel. Instead of a global correlation of piRNA expression and TE content, we found evidence for a host response through de novo piRNA production from novel TE insertions. Although approximately 20% of novel TE insertions induced de novo piRNA production, the abundance of de novo piRNAs was low and did not markedly affect the global pool of ovarian piRNAs. Our results provide new insights into the evolution of TEs and the piRNA system in an important model organism.

Description: Haloarchaeal genomes are generally composed of multiple replicons, and each replicon has a single or multiple replication origin(s). The comparative genomic analysis of replication origins from closely related species can be used to reveal the evolutionary mechanisms that account for the development of multiple origin systems. Multiple replication origins have been in silico and experimentally investigated in Haloarcula hispanica , which raise the possibility for comparisons of multiple replication origins in Haloarcula species. Thus, we performed a comparison of H. hispanica replication origins with those from five additional Haloarcula species. We demonstrated that the multiple replication origins in the chromosome were evolved independently multiple times from the oriC1 -dependent ancestral chromosome. Particularly, the two origins oriC1 and oriC2 were conserved in location, and both of them were adjacent to an rRNA operon, suggestive of correlations in replication and expression of surrounding genes that may promote the conservation of these two origins. Some chromosomal variable regions were used as hotspots for origin evolution in which replication origins were continually being acquired, lost, and disrupted. Furthermore, we demonstrated that autonomously replicating sequence plasmids with H. hispanica minichromosomal replication origins were extremely unstable. Because both organization and replication origins of minichromosomes were not conserved, we proposed an association between the evolution of extrachromosomal replicons and origin variation. Taken together, we provided insights into the evolutionary history of multiple replication origins in Haloarcula species, and proposed a general model of association between the dynamics of multiple replication origins and the evolution of multireplicon genome architecture in haloarchaea.

Description: Hsp90 is one of the most abundant and conserved proteins in the cell. Reduced levels or activity of Hsp90 causes defects in many cellular processes and also reveals genetic and nongenetic variation within a population. Despite information about Hsp90 protein–protein interactions, a global view of the Hsp90-regulated proteome in yeast is unavailable. To investigate the degree of dependency of individual yeast proteins on Hsp90, we used the "stable isotope labeling by amino acids in cell culture" method coupled with mass spectrometry to quantify around 4,000 proteins in low-Hsp90 cells. We observed that 904 proteins changed in their abundance by more than 1.5-fold. When compared with the transcriptome of the same population of cells, two-thirds of the misregulated proteins were observed to be affected posttranscriptionally, of which the majority were downregulated. Further analyses indicated that the downregulated proteins are highly conserved and assume central roles in cellular networks with a high number of protein interacting partners, suggesting that Hsp90 buffers genetic and nongenetic variation through regulating protein network hubs. The downregulated proteins were enriched for essential proteins previously not known to be Hsp90-dependent. Finally, we observed that downregulation of transcription factors and mating pathway components by attenuating Hsp90 function led to decreased target gene expression and pheromone response, respectively, providing a direct link between observed proteome regulation and cellular phenotypes.

Description: Streptococcus sinensis is a recently discovered human pathogen isolated from blood cultures of patients with infective endocarditis. Its phylogenetic position, as well as those of its closely related species, remains inconclusive when single genes were used for phylogenetic analysis. For example, S. sinensis branched out from members of the anginosus, mitis, and sanguinis groups in the 16S ribosomal RNA gene phylogenetic tree, but it was clustered with members of the anginosus and sanguinis groups when groEL gene sequences used for analysis. In this study, we sequenced the draft genome of S. sinensis and used a polyphasic approach, including concatenated genes, whole genomes, and matrix-assisted laser desorption ionization-time of flight mass spectrometry to analyze the phylogeny of S. sinensis . The size of the S. sinensis draft genome is 2.06 Mb, with GC content of 42.2%. Phylogenetic analysis using 50 concatenated genes or whole genomes revealed that S. sinensis formed a distinct cluster with Streptococcus oligofermentans and Streptococcus cristatus , and these three streptococci were clustered with the "sanguinis group." As for phylogenetic analysis using hierarchical cluster analysis of the mass spectra of streptococci, S. sinensis also formed a distinct cluster with S. oligofermentans and S. cristatus , but these three streptococci were clustered with the "mitis group." On the basis of the findings, we propose a novel group, named "sinensis group," to include S. sinensis , S. oligofermentans, and S. cristatus , in the Streptococcus genus. Our study also illustrates the power of phylogenomic analyses for resolving ambiguities in bacterial taxonomy.

Description: Understanding the relative contributions of various evolutionary processes—purifying selection, neutral drift, and adaptation—is fundamental to evolutionary biology. A common metric to distinguish these processes is the ratio of nonsynonymous to synonymous substitutions (i.e., d N /d S ) interpreted from the neutral theory as a null model. However, from biophysical considerations, mutations have non-negligible effects on the biophysical properties of proteins such as folding stability. In this work, we investigated how stability affects the rate of protein evolution in phylogenetic trees by using simulations that combine explicit protein sequences with associated stability changes. We first simulated myoglobin evolution in phylogenetic trees with a biophysically realistic approach that accounts for 3D structural information and estimates of changes in stability upon mutation. We then compared evolutionary rates inferred directly from simulation to those estimated using maximum-likelihood (ML) methods. We found that the d N /d S estimated by ML methods ( ML ) is highly predictive of the per gene d N /d S inferred from the simulated phylogenetic trees. This agreement is strong in the regime of high stability where protein evolution is neutral. At low folding stabilities and under mutation-selection balance, we observe deviations from neutrality (per gene d N /d S 〉 1 and d N /d S 〈 1). We showed that although per gene d N /d S is robust to these deviations, ML tests for positive selection detect statistically significant per site d N /d S 〉 1. Altogether, we show how protein biophysics affects the dN/dS estimations and its subsequent interpretation. These results are important for improving the current approaches for detecting positive selection.

Description: The albuminoid gene family comprises vitamin D-binding protein ( GC ), alpha-fetoprotein ( AFP ), afamin ( AFM ), and albumin ( ALB ). Albumin is the most abundant human serum protein, and, as the other family members, acts as a transporter of endogenous and exogenous substances including thyroxine, fatty acids, and drugs. Instead, the major cargo of GC is 25-hydroxyvitamin D. We performed an evolutionary study of albuminoid genes and we show that ALB evolved adaptively in mammals. Most positively selected sites are located within albumin-binding sites for fatty acids and thyroxine, as well as at the contact surface with neonatal Fc receptor. Positive selection was also detected for residues forming the prostaglandin-binding pocket. Adaptation to hibernation/torpor might explain the signatures of episodic positive selection we detected for few mammalian lineages. Application of a population genetics–phylogenetics approach showed that purifying selection represented a major force acting on albuminoid genes in both humans and chimpanzees, with the strongest constraint observed for human GC . Population genetic analysis revealed that GC was also the target of locally exerted selective pressure, which drove the frequency increase of different haplotypes in distinct human populations. A search for known variants that modulate GC and 25-hydroxyvitamin D concentrations revealed linkage disequilibrium with positively selected variants, although European and Asian major GC haplotypes carry alleles with reported opposite effect on GC concentration. Data herein indicate that albumin, an extremely abundant housekeeping protein, was the target of pervasive and episodic selection in mammals, whereas GC represented a selection target during the recent evolution of human populations.

Description: DNA methylation in the genome plays a fundamental role in the regulation of gene expression and is widespread in the genome of eukaryotic species. For example, in higher vertebrates, there is a "global" methylation pattern involving complete methylation of CpG sites genome-wide, except in promoter regions that are typically enriched for CpG dinucleotides, or so called "CpG islands." Here, we comprehensively examined and compared the distribution of CpG sites within ten model eukaryotic species and linked the observed patterns to the role of DNA methylation in controlling gene transcription. The analysis revealed two distinct but conserved methylation patterns for gene promoters in human and mouse genomes, involving genes with distinct distributions of promoter CpGs and gene expression patterns. Comparative analysis with four other higher vertebrates revealed that the primary regulatory role of the DNA methylation system is highly conserved in higher vertebrates.

Description: The bread wheat ( Triticum aestivum L.) genotype "Chinese Spring" ("CS") is the reference base in wheat genetics and genomics. Pericentric rearrangements in this genotype were systematically assessed by analyzing homoeoloci for a set of nonredundant genes from Brachypodium distachyon , Triticum urartu , and Aegilops tauschii in the CS chromosome shotgun sequence obtained from individual chromosome arms flow-sorted from CS aneuploid lines. Based on patterns of their homoeologous arm locations, 551 genes indicated the presence of pericentric inversions in at least 10 of the 21 chromosomes. Available data from deletion bin-mapped expressed sequence tags and genetic mapping in wheat indicated that all inversions had breakpoints in the low-recombinant gene-poor pericentromeric regions. The large number of putative intrachromosomal rearrangements suggests the presence of extensive structural differences among the three subgenomes, at least some of which likely occurred during the production of the aneuploid lines of this hexaploid wheat genotype. These differences could have significant implications in wheat genome research where comparative approaches are used such as in ordering and orientating sequence contigs and in gene cloning.

Description: β-defensins are a family of important peptides of innate immunity, involved in host defense, immunomodulation, reproduction, and pigmentation. Genes encoding β-defensins show evidence of birth-and-death evolution, adaptation by amino acid sequence changes, and extensive copy number variation (CNV) within humans and other species. The role of CNV in the adaptation of β-defensins to new functions remains unclear, as does the adaptive role of CNV in general. Here, we fine-map CNV of a cluster of β-defensins in humans and rhesus macaques. Remarkably, we found that the structure of the CNV is different between primates, with distinct mutational origins and CNV boundaries defined by retroviral long terminal repeat elements. Although the human β-defensin CNV region is 322 kb and encompasses several genes, including β-defensins, a long noncoding RNA gene, and testes-specific zinc-finger transcription factors, the orthologous region in the rhesus macaque shows CNV of a 20-kb region, containing only a single gene, the ortholog of the human β-defensin-2 gene. Despite its independent origins, the range of gene copy numbers in the rhesus macaque is similar to humans. In addition, the rhesus macaque gene has been subject to divergent positive selection at the amino acid level following its initial duplication event between 3 and 9.5 Ma, suggesting adaptation of this gene as the macaque successfully colonized novel environments outside Africa. Therefore, the molecular phenotype of β-defensin-2 CNV has undergone convergent evolution, and this gene shows evidence of adaptation at the amino acid level in rhesus macaques.

Description: Translation elongation factor-1 alpha (EF1A) and the related GTPase EF-like (EFL) are two proteins with a complex mutually exclusive distribution across the tree of eukaryotes. Recent surveys revealed that the distribution of the two GTPases in even closely related taxa is frequently at odds with their phylogenetic relationships. Here, we investigate the distribution of EF1A and EFL in the alveolate supergroup. Alveolates comprise three major lineages: ciliates and apicomplexans encode EF1A, whereas dinoflagellates encode EFL. We searched transcriptome databases for seven early-diverging alveolate taxa that do not belong to any of these groups: colpodellids, chromerids, and colponemids. Current data suggest all seven are expected to encode EF1A, but we find three genera encode EFL: Colpodella , Voromonas , and the photosynthetic Chromera . Comparing this distribution with the phylogeny of alveolates suggests that EF1A and EFL evolution in alveolates cannot be explained by a simple horizontal gene transfer event or lineage sorting.

Description: The high-altitude hypoxic environment represents one of the most extreme challenges for mammals. Previous studies of humans on the Tibetan plateau and in the Andes Mountains have identified statistical signatures of selection in different sets of loci. Here, we first measured the hemoglobin levels in village dogs from Tibet and those from Chinese lowlands. We found that the hemoglobin levels are very similar between the two groups, suggesting that Tibetan dogs might share similar adaptive strategies as the Tibetan people. Through a whole-genome sequencing approach, we have identified EPAS1 and HBB as candidate genes for the hypoxic adaptation on the Tibetan plateau. The population genetic analysis shows a significant convergence between humans and dogs in Tibet. The similarities in the sets of loci that exhibit putative signatures of selection and the hemoglobin levels between humans and dogs of the same environment, but not between human populations in different regions, suggests an extraordinary landscape of convergent evolution between human beings and their best friend on the Tibetan plateau.

Description: Recent studies in paleovirology have uncovered myriads of endogenous viral elements (EVEs) integrated in the genome of their eukaryotic hosts. These fragments result from endogenization, that is, integration of the viral genome into the host germline genome followed by vertical inheritance. So far, most studies have used a virus-centered approach, whereby endogenous copies of a particular group of viruses were searched in all available sequenced genomes. Here, we follow a host-centered approach whereby the genome of a given species is comprehensively screened for the presence of EVEs using all available complete viral genomes as queries. Our analyses revealed that 54 EVEs corresponding to 10 different viral lineages belonging to 5 viral families ( Bunyaviridae , Circoviridae , Parvoviridae , and Totiviridae ) and one viral order ( Mononegavirales ) became endogenized in the genome of the isopod crustacean Armadillidium vulgare . We show that viral endogenization occurred recurrently during the evolution of isopods and that A. vulgare viral lineages were involved in multiple host switches that took place between widely divergent taxa. Furthermore, 30 A. vulgare EVEs have uninterrupted open reading frames, suggesting they result from recent endogenization of viruses likely to be currently infecting isopod populations. Overall, our work shows that isopods have been and are still infected by a large variety of viruses. It also extends the host range of several families of viruses and brings new insights into their evolution. More generally, our results underline the power of paleovirology in characterizing the viral diversity currently infecting eukaryotic taxa.

Description: Tetrapyrroles such as heme and chlorophyll are essential for biological processes, including oxygenation, respiration, and photosynthesis. In the tetrapyrrole biosynthesis pathway, protoporphyrinogen IX oxidase (Protox) catalyzes the formation of protoporphyrin IX, the last common intermediate for the biosynthesis of heme and chlorophyll. Three nonhomologous isofunctional enzymes, HemG, HemJ, and HemY, for Protox have been identified. To reveal the distribution and evolution of the three Protox enzymes, we identified homologs of each along with other heme biosynthetic enzymes by whole-genome clustering across three domains of life. Most organisms possess only one of the three Protox types, with some exceptions. Detailed phylogenetic analysis revealed that HemG is mostly limited to -Proteobacteria whereas HemJ may have originated within α-Proteobacteria and transferred to other Proteobacteria and Cyanobacteria. In contrast, HemY is ubiquitous in prokaryotes and is the only Protox in eukaryotes, so this type may be the ancestral Protox. Land plants have a unique HemY homolog that is also shared by Chloroflexus species, in addition to the main HemY homolog originating from Cyanobacteria. Meanwhile, organisms missing any Protox can be classified into two groups; those lacking most heme synthetic genes, which necessarily depend on external heme supply, and those lacking only genes involved in the conversion of uroporphyrinogen III into heme, which would use a precorrin2-dependent alternative pathway. However, hemN encoding coproporphyrinogen IX oxidase was frequently found in organisms lacking Protox enzyme, which suggests a unique role of this gene other than in heme biosynthesis.

Description: Bartonellae are mammalian pathogens vectored by blood-feeding arthropods. Although of increasing medical importance, little is known about their ecological past, and host associations are underexplored. Previous studies suggest an influence of horizontal gene transfers in ecological niche colonization by acquisition of host pathogenicity genes. We here expand these analyses to metabolic pathways of 28 Bartonella genomes, and experimentally explore the distribution of bartonellae in 21 species of blood-feeding arthropods. Across genomes, repeated gene losses and horizontal gains in the phospholipid pathway were found. The evolutionary timing of these patterns suggests functional consequences likely leading to an early intracellular lifestyle for stem bartonellae. Comparative phylogenomic analyses discover three independent lineage-specific reacquisitions of a core metabolic gene—NAD(P)H-dependent glycerol-3-phosphate dehydrogenase ( gpsA )—from Gammaproteobacteria and Epsilonproteobacteria. Transferred genes are significantly closely related to invertebrate Arsenophonus -, and Serratia -like endosymbionts , and mammalian Helicobacter -like pathogens, supporting a cellular association with arthropods and mammals at the base of extant Bartonella spp. Our studies suggest that the horizontal reacquisitions had a key impact on bartonellae lineage specific ecological and functional evolution.

Description: The majority of our knowledge about mitochondrial genomes of Viridiplantae comes from land plants, but much less is known about their green algal relatives. In the green algal order Sphaeropleales (Chlorophyta), only one representative mitochondrial genome is currently available—that of Acutodesmus obliquus . Our study adds nine completely sequenced and three partially sequenced mitochondrial genomes spanning the phylogenetic diversity of Sphaeropleales. We show not only a size range of 25–53 kb and variation in intron content (0–11) and gene order but also conservation of 13 core respiratory genes and fragmented ribosomal RNA genes. We also report an unusual case of gene arrangement convergence in Neochloris aquatica , where the two rns fragments were secondarily placed in close proximity. Finally, we report the unprecedented usage of UCG as stop codon in Pseudomuriella schumacherensis . In addition, phylogenetic analyses of the mitochondrial protein-coding genes yield a fully resolved, well-supported phylogeny, showing promise for addressing systematic challenges in green algae.

Description: Distributive conjugal transfer (DCT) is a newly described mechanism of lateral gene transfer (LGT) that results in a mosaic transconjugant structure, similar to the products of meiosis. We have tested popular LGT detection methods on whole-genome sequence data from experimental DCT transconjugants and used the best performing methods to compare genomic signatures of DCT with those of LGT through natural transformation, conjugative plasmids, and mobile genetic elements (MGE). We found that DCT results in transfer of larger chromosomal segments, that these segments are distributed more broadly around the chromosome, and that a greater proportion of the chromosome is affected by DCT than by other mechanisms of LGT. We used the best performing methods to characterize LGT in Mycobacterium canettii , the mycobacterial species most closely related to Mycobacterium tuberculosis . Patterns of LGT among M. canettii were highly distinctive. Gene flow appeared unidirectional, from lineages with minimal evidence of LGT to isolates with a substantial proportion (6–13%) of sites identified as recombinant. Among M. canettii isolates with evidence of LGT, recombinant fragments were larger and more evenly distributed relative to bacteria that undergo LGT through natural transformation, conjugative plasmids, and MGE. Spatial bias in M. canettii was also unusual in that patterns of recombinant fragment sharing mirrored overall phylogenetic structure. Based on the proportion of recombinant sites, the size of recombinant fragments, their spatial distribution and lack of association with MGE, as well as unidirectionality of DNA transfer, we conclude that DCT is the predominant mechanism of LGT among M. canettii .

Description: The spatiotemporal control of gene expression is crucial for the successful completion of animal development. The evolutionary constraints on development are particularly strong for the mid-embryonic stage when body segments are specified, as evidenced by a high degree of morphological and protein-coding gene conservation during this period—a phenomenon known as the developmental hourglass. The discovery of microRNA-mediated gene control revealed an entirely new layer of complexity of the molecular networks that orchestrate development. However, the constraints on microRNA developmental expression and evolution, and the implications for animal evolution are less well understood. To systematically explore the conservation of microRNAs during development, we carried out a genome-wide comparative study of microRNA expression levels throughout the ontogenesis of two divergent fruit flies, Drosophila melanogaster and D. virilis . We show that orthologous microRNAs display highly similar temporal profiles regardless of their mutation rates, suggesting that the timely expression of microRNA genes can be more constrained than their sequence. Furthermore, transitions between key developmental events in the different species are accompanied by conserved shifts in microRNA expression profiles, with the mid-embryonic period between gastrulation and segmentation characterized by the highest similarity of microRNA expression. The conservation of microRNA expression therefore displays an hourglass pattern similar to that observed for protein-coding genes.

Description: Rice wine fermentation represents a unique environment for the evolution of the budding yeast, Saccharomyces cerevisiae . To understand how the selection pressure shaped the yeast genome and gene regulation, we determined the genome sequence and transcriptome of a S. cerevisiae strain YHJ7 isolated from Chinese rice wine (Huangjiu), a popular traditional alcoholic beverage in China. By comparing the genome of YHJ7 to the lab strain S288c, a Japanese sake strain K7, and a Chinese industrial bioethanol strain YJSH1, we identified many genomic sequence and structural variations in YHJ7, which are mainly located in subtelomeric regions, suggesting that these regions play an important role in genomic evolution between strains. In addition, our comparative transcriptome analysis between YHJ7 and S288c revealed a set of differentially expressed genes, including those involved in glucose transport (e.g., HXT 2, HXT7 ) and oxidoredutase activity (e.g., AAD 10, ADH 7). Interestingly, many of these genomic and transcriptional variations are directly or indirectly associated with the adaptation of YHJ7 strain to its specific niches. Our molecular evolution analysis suggested that Japanese sake strains (K7/UC5) were derived from Chinese rice wine strains (YHJ7) at least approximately 2,300 years ago, providing the first molecular evidence elucidating the origin of Japanese sake strains. Our results depict interesting insights regarding the evolution of yeast during rice wine fermentation, and provided a valuable resource for genetic engineering to improve industrial wine-making strains.

Description: Campylobacter is the most common bacterial cause of gastroenteritis in the European Union with over 200,000 laboratory-confirmed cases reported annually. This is the first study to describe findings related to comparative genomics analyses of the sequence type (ST)-677 clonal complex (CC), a Campylobacter jejuni lineage associated with bacteremia cases in humans. We performed whole-genome sequencing, using Illumina HiSeq sequencing technology, on five related ST-677 CC isolates from two chicken farms to identify microevolution taking place at the farms. Our further aim was to identify novel putative virulence determinants from the ST-677 CC genomes. For this purpose, clinical isolates of the same CC were included in comparative genomic analyses against well-known reference strains of C. jejuni . Overall, the ST-677 CC was recognized as a highly clonal lineage with relatively small differences between the genomes. Among the farm isolates differences were identified mainly in the lengths of the homopolymeric tracts in genes related to the capsule, lipo-oligosaccharide, and flagella. We identified genomic features shared with C. jejuni subsp. doylei , which has also been shown to be associated with bacteremia in humans. These included the degradation of the cytolethal distending toxin operon and similarities between the capsular polysaccharide biosynthesis loci. The phase-variable GDP-mannose 4,6-dehydratase (EC 4.2.1.47) ( wcbK , CAMP1649), associated with the capsular polysaccharide biosynthesis locus, may play a central role in ST-677 CC conferring acid and serum resistance during different stages of infection. Homology-based searches revealed several additional novel features and characteristics, including two putative type Vb secretion systems and a novel restriction modification/methyltransferase gene cluster, putatively associated with pathogenesis and niche adaptation.

Description: Within the phylum Nematoda, plant-parasitism is hypothesized to have arisen independently on at least four occasions. The most economically damaging plant-parasitic nematode species, and consequently the most widely studied, are those that feed as they migrate destructively through host roots causing necrotic lesions (migratory endoparasites) and those that modify host root tissue to create a nutrient sink from which they feed (sedentary endoparasites). The false root-knot nematode Nacobbus aberrans is the only known species to have both migratory endoparasitic and sedentary endoparasitic stages within its life cycle. Moreover, its sedentary stage appears to have characteristics of both the root-knot and the cyst nematodes. We present the first large-scale genetic resource of any false-root knot nematode species. We use RNAseq to describe relative abundance changes in all expressed genes across the life cycle to provide interesting insights into the biology of this nematode as it transitions between modes of parasitism. A multigene phylogenetic analysis of N. aberrans with respect to plant-parasitic nematodes of all groups confirms its proximity to both cyst and root-knot nematodes. We present a transcriptome-wide analysis of both lateral gene transfer events and the effector complement. Comparing parasitism genes of typical root-knot and cyst nematodes to those of N. aberrans has revealed interesting similarities. Importantly, genes that were believed to be either cyst nematode, or root-knot nematode, "specific" have both been identified in N. aberrans . Our results provide insights into the characteristics of a common ancestor and the evolution of sedentary endoparasitism of plants by nematodes.

Description: Complex nongradual evolutionary processes such as gene remodeling are difficult to model, to visualize, and to investigate systematically. Despite these challenges, the creation of composite (or mosaic) genes by combination of genetic segments from unrelated gene families was established as an important adaptive phenomena in eukaryotic genomes. In contrast, almost no general studies have been conducted to quantify composite genes in viruses. Although viral genome mosaicism has been well-described, the extent of gene mosaicism and its rules of emergence remain largely unexplored. Applying methods from graph theory to inclusive similarity networks, and using data from more than 3,000 complete viral genomes, we provide the first demonstration that composite genes in viruses are 1) functionally biased, 2) involved in key aspects of the arm race between cells and viruses, and 3) can be classified into two distinct types of composite genes in all viral classes. Beyond the quantification of the widespread recombination of genes among different viruses of the same class, we also report a striking sharing of genetic information between viruses of different classes and with different nucleic acid types. This latter discovery provides novel evidence for the existence of a large and complex mobilome network, which appears partly bound by the sharing of genetic information and by the formation of composite genes between mobile entities with different genetic material. Considering that there are around 10E31 viruses on the planet, gene remodeling appears as a hugely significant way of generating and moving novel sequences between different kinds of organisms on Earth.

Description: The alpha subunits of voltage-gated calcium channels (Ca v s) are large transmembrane proteins responsible for crucial physiological processes in excitable cells. They are assisted by three auxiliary subunits that can modulate their electrical behavior. Little is known about the evolution and roles of the various subunits of Ca v s in nonbilaterian animals and in nonanimal lineages. For this reason, we mapped the phyletic distribution of the four channel subunits and reconstructed their phylogeny. Although alpha subunits have deep evolutionary roots as ancient as the split between plants and opistokonths, beta subunits appeared in the last common ancestor of animals and their close-relatives choanoflagellates, gamma subunits are a bilaterian novelty and alpha2/delta subunits appeared in the lineage of Placozoa, Cnidaria, and Bilateria. We note that gene losses were extremely common in the evolution of Ca v s, with noticeable losses in multiple clades of subfamilies and also of whole Ca v families. As in vertebrates, but not protostomes, Ca v channel genes duplicated in Cnidaria. We characterized by in situ hybridization the tissue distribution of alpha subunits in the sea anemone Nematostella vectensis , a nonbilaterian animal possessing all three Ca v subfamilies common to Bilateria. We find that some of the alpha subunit subtypes exhibit distinct spatiotemporal expression patterns. Further, all six sea anemone alpha subunit subtypes are conserved in stony corals, which separated from anemones 500 MA. This unexpected conservation together with the expression patterns strongly supports the notion that these subtypes carry unique functional roles.

Description: Micro-RNA (miRNA) genes encode abundant small regulatory RNAs that play key roles during development and in homeostasis by fine tuning and buffering gene expression. This layer of regulatory control over transcriptional networks is preserved by selection across deep evolutionary time, yet selection pressures on individual miRNA genes in contemporary populations remain poorly characterized in any organism. Here, we quantify nucleotide variability for 129 miRNAs in the genome of the nematode Caenorhabditis remanei to understand the microevolution of this important class of regulatory genes. Our analysis of three population samples and C. remanei ’s sister species revealed ongoing natural selection that constrains evolution of all sequence domains within miRNA hairpins. We also show that new miRNAs evolve faster than older miRNAs but that selection nevertheless favors their persistence. Despite the ongoing importance of purging of new mutations, we discover a trove of 〉400 natural miRNA sequence variants that include single nucleotide polymorphisms in seed motifs, indels that ablate miRNA functional domains, and origination of new miRNAs by duplication. Moreover, we demonstrate substantial nucleotide divergence of pre-miRNA hairpin alleles between populations and sister species. These findings from the first global survey of miRNA microevolution in Caenorhabditis support the idea that changes in gene expression, mediated through divergence in miRNA regulation, can contribute to phenotypic novelty and adaptation to specific environments in the present day as well as the distant past.

Description: Respiratory complexes are encoded by two genomes (mitochondrial DNA [mtDNA] and nuclear DNA [nDNA]). Although the importance of intergenomic coadaptation is acknowledged, the forces and constraints shaping such coevolution are largely unknown. Previous works using cytochrome c oxidase (COX) as a model enzyme have led to the so-called "optimizing interaction" hypothesis. According to this view, mtDNA-encoded residues close to nDNA-encoded residues evolve faster than the rest of positions, favoring the optimization of protein–protein interfaces. Herein, using evolutionary data in combination with structural information of COX, we show that failing to discern the effects of interaction from other structural and functional effects can lead to deceptive conclusions such as the "optimizing hypothesis." Once spurious factors have been accounted for, data analysis shows that mtDNA-encoded residues engaged in contacts are, in general, more constrained than their noncontact counterparts. Nevertheless, noncontact residues from the surface of COX I subunit are a remarkable exception, being subjected to an exceptionally high purifying selection that may be related to the maintenance of a suitable heme environment. We also report that mtDNA-encoded residues involved in contacts with other mtDNA-encoded subunits are more constrained than mtDNA-encoded residues interacting with nDNA-encoded polypeptides. This differential behavior cannot be explained on the basis of predicted thermodynamic stability, as interactions between mtDNA-encoded subunits contribute more weakly to the complex stability than those interactions between subunits encoded by different genomes. Therefore, the higher conservation observed among mtDNA-encoded residues involved in intragenome interactions is likely due to factors other than structural stability.

Description: Nucleotide-binding oligomerization domain (NOD)-like receptors (NLRs) are intracellular receptors that control innate immunity and other biotic interactions in animals and plants. NLRs have been characterized in plant and animal lineages, but in fungi, this gene family has not been systematically described. There is however previous indications of the involvement of NLR-like genes in nonself recognition and programmed cell death in fungi. We have analyzed 198 fungal genomes for the presence of NLRs and have annotated a total of 5,616 NLR candidates. We describe their phylogenetic distribution, domain organization, and evolution. Fungal NLRs are characterized by a great diversity of domain organizations, suggesting frequently occurring combinatorial assortments of different effector, NOD and repeat domains. The repeat domains are of the WD, ANK, and TPR type; no LRR motifs were found. As previously documented for WD-repeat domains of fungal NLRs, TPR, and ANK repeats evolve under positive selection and show highly conserved repeats and repeat length polymorphism, suggesting the possibility of concerted evolution of these repeats. We identify novel effector domains not previously found associated with NLRs, whereas others are related to effector domains of plant or animals NLRs. In particular, we show that the HET domain found in fungal NLRs may be related to Toll/interleukin-1 receptor domains found in animal and plant immune receptors. This description of fungal NLR repertoires reveals both similarities and differences with plant and animals NLR collections, highlights the importance of domain reassortment and repeat evolution and provides a novel entry point to explore the evolution of NLRs in eukaryotes.

Description: Mongolians have played a significant role in modern human evolution, especially after the rise of Genghis Khan (1162[?]–1227). Although the social cultural impacts of Genghis Khan and the Mongolian population have been well documented, explorations of their genome structure and genetic imprints on other human populations have been lacking. We here present the genome of a Mongolian male individual. The genome was de novo assembled using a total of 130.8-fold genomic data produced from massively parallel whole-genome sequencing. We identified high-confidence variation sets, including 3.7 million single nucleotide polymorphisms (SNPs) and 756,234 short insertions and deletions. Functional SNP analysis predicted that the individual has a pathogenic risk for carnitine deficiency. We located the patrilineal inheritance of the Mongolian genome to the lineage D3a through Y haplogroup analysis and inferred that the individual has a common patrilineal ancestor with Tibeto-Burman populations and is likely to be the progeny of the earliest settlers in East Asia. We finally investigated the genetic imprints of Mongolians on other human populations using different approaches. We found varying degrees of gene flows between Mongolians and populations living in Europe, South/Central Asia, and the Indian subcontinent. The analyses demonstrate that the genetic impacts of Mongolians likely resulted from the expansion of the Mongolian Empire in the 13th century. The genome will be of great help in further explorations of modern human evolution and genetic causes of diseases/traits specific to Mongolians.

Description: The symptoms of Clostridium difficile infection are caused by two closely related toxins, TcdA and TcdB, which are encoded by the 19.6 kb Pathogenicity Locus (PaLoc). The PaLoc is variably present among strains, and in this respect it resembles a mobile genetic element. The C. difficile population structure consists mainly of five phylogenetic clades designated 1–5. Certain genotypes of clade 5 are associated with recently emergent highly pathogenic strains causing human disease and animal infections. The aim of this study was to explore the evolutionary history of the PaLoc in C. difficile clade 5. Phylogenetic analyses and annotation of clade 5 PaLoc variants and adjoining genomic regions were undertaken using a representative collection of toxigenic and nontoxigenic strains. Comparison of the core genome and PaLoc phylogenies obtained for clade 5 and representatives of the other clades identified two distinct PaLoc acquisition events, one involving a toxin A + B + PaLoc variant and the other an A – B + variant. Although the exact mechanism of each PaLoc acquisition is unclear, evidence of possible homologous recombination with other clades and between clade 5 lineages was found within the PaLoc and adjacent regions. The generation of nontoxigenic variants by PaLoc loss via homologous recombination with PaLoc-negative members of other clades was suggested by analysis of cdu2 , although none is likely to have occurred recently. A variant of the putative holin gene present in the clade 5 A – B + PaLoc was likely acquired via allelic exchange with an unknown element. Fine-scale phylogenetic analysis of C. difficile clade 5 revealed the extent of its genetic diversity, consistent with ancient evolutionary origins and a complex evolutionary history for the PaLoc.