ALBERT

Description: The rapid and widespread evolution of fungicide resistance remains a challenge for crop disease management. The demethylation inhibitor (DMI) class of fungicides is a widely used chemistry for managing disease, but there has been a gradual decline in efficacy in many crop pathosystems. Reliance on DMI fungicides has increased resistance in populations of the plant pathogenic fungus Cercospora beticola worldwide. To better understand the genetic and evolutionary basis for DMI resistance in C. beticola, a genome-wide association study (GWAS) and selective sweep analysis were conducted for the first time in this species. We performed whole-genome resequencing of 190 C. beticola isolates infecting sugar beet (Beta vulgaris ssp. vulgaris). All isolates were phenotyped for sensitivity to the DMI tetraconazole. Intragenic markers on chromosomes 1, 4, and 9 were significantly associated with DMI fungicide resistance, including a polyketide synthase gene and the gene encoding the DMI target CbCYP51. Haplotype analysis of CbCYP51 identified a synonymous mutation (E170) and nonsynonymous mutations (L144F, I387M, and Y464S) associated with DMI resistance. Genome-wide scans of selection showed that several of the GWAS mutations for fungicide resistance resided in regions that have recently undergone a selective sweep. Using radial plate growth on selected media as a fitness proxy, we did not find a trade-off associated with DMI fungicide resistance. Taken together, we show that population genomic data from a crop pathogen can allow the identification of mutations conferring fungicide resistance and inform about their origins in the pathogen population.

Description: A growing body of evidence has underscored the role of horizontal gene transfer (HGT) in animal evolution. Previously, we discovered the horizontal transfer of the gene encoding the eukaryotic genotoxin cytolethal distending toxin B (cdtB) from the pea aphid Acyrthosiphon pisum secondary endosymbiont (APSE) phages to drosophilid and aphid nuclear genomes. Here, we report cdtB in the nuclear genome of the gall-forming “swede midge” Contarinia nasturtii (Diptera: Cecidomyiidae) via HGT. We searched all available gall midge genome sequences for evidence of APSE-to-insect HGT events and found five toxin genes (aip56, cdtB, lysozyme, rhs, and sltxB) transferred horizontally to cecidomyiid nuclear genomes. Surprisingly, phylogenetic analyses of HGT candidates indicated APSE phages were often not the ancestral donor lineage of the toxin gene to cecidomyiids. We used a phylogenetic signal statistic to test a transfer-by-proximity hypothesis for animal HGT, which suggested that microbe-to-insect HGT was more likely between taxa that share environments than those from different environments. Many of the toxins we found in midge genomes target eukaryotic cells, and catalytic residues important for toxin function are conserved in insect copies. This class of horizontally transferred, eukaryotic cell-targeting genes is potentially important in insect adaptation.

Description: Authentic DNA sequences are crucial for reliable evolutionary inference. Concerns about the identification of DNA sequences have been voiced several times in the past but few quantitative studies exist. Mitogenomes play important roles in phylogenetics, phylogeography, population genetics, and DNA identification. However, the large number of mitogenomes being published routinely, often in brief data papers, has raised questions about their authenticity. In this study, we quantify problematic mitogenomes of birds and their reusage in other papers. Of 1,876 complete or partial mitogenomes of birds published until January 1, 2020, the authenticity of 1,559 could be assessed with sequences of conspecifics. Of these, 78 (5.0%) were found to be problematic, including 45 curated reference sequences. Problems were due to misidentification (33), chimeras of two or three species (23), sequencing errors/numts (18), incorrect sequence assembly (1), mislabeling at GenBank but not in the final paper (2), or vice versa (1). The number of problematic mitogenomes has increased sharply since 2012. Worryingly, these problematic sequences have been reused 436 times in other papers, including 385 times in phylogenies. No less than 53% of all mitogenomic phylogenies/networks published until January 1, 2020 included at least one problematic mitogenome. Problematic mitogenomes have resulted in incorrect phylogenetic hypotheses and proposals for unwarranted taxonomic revision, and may have compromised comparative analyses and measurements of divergence times. Our results indicate that a major upgrade of quality control measures is warranted. We propose a comprehensive set of measures that may serve as a new standard for publishing mitogenome sequences.

Description: Nucleotides are an important class of odorants for aquatic vertebrates such as frogs and fishes, but also have manifold signaling roles in other cellular processes. Recently, an adenosine receptor believed to belong to the adora2 clade has been identified as an olfactory receptor in zebrafish. Here, we set out to elucidate the evolutionary history of both this gene and its olfactory function. We have performed a thorough phylogenetic study in vertebrates, chordates and their sister group, ambulacraria, and show that the origin of the zebrafish olfactory receptor gene can be traced back to the most recent common ancestor of all three groups as a segregate sister clade (adorb) to the adora gene family. Eel, carp, and clawed frog all express adorb in a sparse and distributed pattern within their olfactory epithelium very similar to the pattern observed for zebrafish that is, consistent with a function as olfactory receptor. In sharp contrast, lamprey adorb-expressing cells are absent from the sensory region of the lamprey nose, but form a contiguous domain directly adjacent to the sensory region. Double-labeling experiments confirmed the expression of lamprey adorb in nonneuronal cells and are consistent with an expression in neuronal progenitor cells. Thus, adorb may have undergone a switch of function in the jawed lineage of vertebrates towards a role as olfactory receptor.

Description: Microbial studies typically involve the sequencing and assembly of draft genomes for individual microbes or whole microbiomes. Given a draft genome, one first task is to determine its phylogenetic context, that is, to place it relative to the set of related reference genomes. We provide a new interactive graphical tool that addresses this task using Mash sketches to compare against all bacterial and archaeal representative genomes in the Genome Taxonomy Database taxonomy, all within the framework of SplitsTree5. The phylogenetic context of the query sequences is then displayed as a phylogenetic outline, a new type of phylogenetic network that is more general than a phylogenetic tree, but significantly less complex than other types of phylogenetic networks. We propose to use such networks, rather than trees, to represent phylogenetic context, because they can express uncertainty in the placement of taxa, whereas a tree must always commit to a specific branching pattern. We illustrate the new method using a number of draft genomes of different assembly quality.

Description: The giant black tiger shrimp (Penaeus monodon) is native to the Indo-Pacific and is the second most farmed penaeid shrimp species globally. Understanding genetic structure, connectivity, and local adaptation among Indo-Pacific black tiger shrimp populations is important for informing sustainable fisheries management and aquaculture breeding programs. Population genetic and outlier detection analyses were undertaken using 10,593 genome-wide single nucleotide polymorphisms (SNPs) from 16 geographically disparate Indo-Pacific P. monodon populations. Levels of genetic diversity were highest for Southeast Asian populations and were lowest for Western Indian Ocean (WIO) populations. Both neutral (n = 9,930) and outlier (n = 663) loci datasets revealed a pattern of strong genetic structure of P. monodon corresponding with broad geographical regions and clear genetic breaks among samples within regions. Neutral loci revealed seven genetic clusters and the separation of Fiji and WIO clusters from all other clusters, whereas outlier loci revealed six genetic clusters and high genetic differentiation among populations. The neutral loci dataset estimated five migration events that indicated migration to Southeast Asia from the WIO, with partial connectivity to populations in both oceans. We also identified 26 putatively adaptive SNPs that exhibited significant Pearson correlation (P 〈 0.05) between minor allele frequency and maximum or minimum sea surface temperature. Matched transcriptome contig annotations suggest putatively adaptive SNPs involvement in cellular and metabolic processes, pigmentation, immune response, and currently unknown functions. This study provides novel genome-level insights that have direct implications for P. monodon aquaculture and fishery management practices.

Description: Many animals including birds, reptiles, insects, and teleost fishes can see ultraviolet (UV) light (shorter than 400 nm), which has functional importance for foraging and communication. For coral reef fishes, shallow reef environments transmit a broad spectrum of light, rich in UV, driving the evolution of diverse spectral sensitivities. However, the identities and sites of the specific visual genes that underly vision in reef fishes remain elusive and are useful in determining how evolution has tuned vision to suit life on the reef. We investigated the visual systems of 11 anemonefish (Amphiprioninae) species, specifically probing for the molecular pathways that facilitate UV-sensitivity. Searching the genomes of anemonefishes, we identified a total of eight functional opsin genes from all five vertebrate visual opsin subfamilies. We found rare instances of teleost UV-sensitive SWS1 opsin gene duplications that produced two functionally coding paralogs (SWS1α and SWS1β) and a pseudogene. We also found separate green sensitive RH2A opsin gene duplicates not yet reported in the family Pomacentridae. Transcriptome analysis revealed false clown anemonefish (Amphiprion ocellaris) expressed one rod opsin (RH1) and six cone opsins (SWS1β, SWS2B, RH2B, RH2A-1, RH2A-2, LWS) in the retina. Fluorescent in situ hybridization highlighted the (co-)expression of SWS1β with SWS2B in single cones, and either RH2B, RH2A, or RH2A together with LWS in different members of double cone photoreceptors (two single cones fused together). Our study provides the first in-depth characterization of visual opsin genes found in anemonefishes and provides a useful basis for the further study of UV-vision in reef fishes.

Description: Species of infraorder Gryllidea, or crickets, are useful invertebrate models for studying developmental biology and neuroscience. They have also attracted attention as alternative protein sources for human food and animal feed. Mitochondrial genomic information on related invertebrates, such as katydids, and locusts, has recently become available in attempt to clarify the controversial classification schemes, although robust phylogenetic relationships with emphasis on crickets remain elusive. Here, we report newly sequenced complete mitochondrial genomes of crickets to study their phylogeny, genomic rearrangements, and adaptive evolution. First, we conducted de novo assembly of mitochondrial genomes from eight cricket species and annotated protein-coding genes and transfer and ribosomal RNAs using automatic annotations and manual curation. Next, by combining newly described protein-coding genes with public data of the complete Gryllidea genomes and gene annotations, we performed phylogenetic analysis and found gene order rearrangements in several branches. We further analyzed genetic signatures of selection in ant-loving crickets (Myrmecophilidae), which are small wingless crickets that inhabit ant nests. Three distinct approaches revealed two positively selected sites in the cox1 gene in these crickets. Protein 3D structural analyses suggested that these selected sites could influence the interaction of respiratory complex proteins, conferring benefits to ant-loving crickets with a unique ecological niche and morphology. These findings enhance our understanding of the genetic basis of cricket evolution without relying on estimates based on a limited number of molecular markers.

Description: Assessing the compatibility between gene family phylogenies is a crucial and often computationally demanding step in many phylogenomic analyses. Here, we describe the Evolutionary Similarity Index (IES), a means to assess shared evolution between gene families using a weighted orthogonal distance regression model applied to sequence distances. The utilization of pairwise distance matrices circumvents comparisons between gene tree topologies, which are inherently uncertain and sensitive to evolutionary model choice, phylogenetic reconstruction artifacts, and other sources of error. Furthermore, IES enables the many-to-many pairing of multiple copies between similarly evolving gene families. This is done by selecting non-overlapping pairs of copies, one from each assessed family, and yielding the least sum of squared residuals. Analyses of simulated gene family data sets show that IES’s accuracy is on par with popular tree-based methods while also less susceptible to noise introduced by sequence alignment and evolutionary model fitting. Applying IES to an empirical data set of 1,322 genes from 42 archaeal genomes identified eight major clusters of gene families with compatible evolutionary trends. The most cohesive cluster consisted of 62 genes with compatible evolutionary signal, which occur as both single-copy and multiple homologs per genome; phylogenetic analysis of concatenated alignments from this cluster produced a tree closely matching previously published species trees for Archaea. Four other clusters are mainly composed of accessory genes with limited distribution among Archaea and enriched toward specific metabolic functions. Pairwise evolutionary distances obtained from these accessory gene clusters suggest patterns of interphyla horizontal gene transfer. An IES implementation is available at https://github.com/lthiberiol/evolSimIndex.

Description: There are more than 100 species of American didelphid marsupials (opossums and mouse opossums). Limited genomic resources for didelphids exists, with only two publicly available genome assemblies compared with dozens in the case of their Australasian counterparts. This discrepancy impedes evolutionary and ecological research. To address this gap, we assembled a high-quality chromosome-level genome of the agile gracile mouse opossum (Gracilinanus agilis) using a combination of stLFR sequencing, polishing with mate-pair data, and anchoring onto pseudochromosomes using Hi-C. This species employs a rare life-history strategy, semelparity, and all G. agilis males and most females die at the end of their first breeding season after succumbing to stress and exhaustion. The 3.7-Gb chromosome-level assembly, with 92.6% anchored onto pseudochromosomes, has a scaffold N50 of 683.5 Mb and a contig N50 of 56.9 kb. The genome assembly shows high completeness, with a mammalian BUSCO score of 88.1%. Around 49.7% of the genome contains repetitive elements. Gene annotation yielded 24,425 genes, of which 83.9% were functionally annotated. The G. agilis genome is an important resource for future studies of marsupial biology, evolution, and conservation.

Description: Wolbachia is a widespread, vertically transmitted bacterial endosymbiont known for manipulating arthropod reproduction. Its most common form of reproductive manipulation is cytoplasmic incompatibility (CI), observed when a modification in the male sperm leads to embryonic lethality unless a compatible rescue factor is present in the female egg. CI attracts scientific attention due to its implications for host speciation and in the use of Wolbachia for controlling vector-borne diseases. However, our understanding of CI is complicated by the complexity of the phenotype, whose expression depends on both symbiont and host factors. In the present study, we perform a comparative analysis of nine complete Wolbachia genomes with known CI properties in the same genetic host background, Drosophila simulans STC. We describe genetic differences between closely related strains and uncover evidence that phages and other mobile elements contribute to the rapid evolution of both genomes and phenotypes of Wolbachia. Additionally, we identify both known and novel genes associated with the modification and rescue functions of CI. We combine our observations with published phenotypic information and discuss how variability in cif genes, novel CI-associated genes, and Wolbachia titer might contribute to poorly understood aspects of CI such as strength and bidirectional incompatibility. We speculate that high titer CI strains could be better at invading new hosts already infected with a CI Wolbachia, due to a higher rescue potential, and suggest that titer might thus be a relevant parameter to consider for future strategies using CI Wolbachia in biological control.

Description: Adenosine Deaminases that Act on RNA (ADARs) are RNA editing enzymes that play a dynamic and nuanced role in regulating transcriptome and proteome diversity. This editing can be highly selective, affecting a specific site within a transcript, or nonselective, resulting in hyperediting. ADAR editing is important for regulating neural functions and autoimmunity, and has a key role in the innate immune response to viral infections, where editing can have a range of pro- or antiviral effects and can contribute to viral evolution. Here we examine the role of ADAR editing across a broad range of viral groups. We propose that the effect of ADAR editing on viral replication, whether pro- or antiviral, is better viewed as an axis rather than a binary, and that the specific position of a given virus on this axis is highly dependent on virus- and host-specific factors, and can change over the course of infection. However, more research needs to be devoted to understanding these dynamic factors and how they affect virus-ADAR interactions and viral evolution. Another area which warrants significant attention is the effect of virus-ADAR interactions on host-ADAR interactions, particularly in light of the crucial role of ADAR in regulating neural functions. Answering these questions will be essential to developing our understanding of the relationship between ADAR editing and viral infection. In turn, this will further our understanding of the effects of viruses such as SARS-CoV-2, as well as many others, and thereby influence our approach to treating these deadly diseases.

Description: Many adaptive traits are polygenic and frequently more loci contributing to the phenotype are segregating than needed to express the phenotypic optimum. Experimental evolution with replicated populations adapting to a new controlled environment provides a powerful approach to study polygenic adaptation. Since genetic redundancy often results in non-parallel selection responses among replicates, we propose a modified Evolve and Resequence (E&R) design that maximizes the similarity among replicates. Rather than starting from many founders, we only use two inbred Drosophila melanogaster strains and expose them to a very extreme, hot temperature environment (29 °C). After 20 generations, we detect many genomic regions with a strong, highly parallel selection response in 10 evolved replicates. The X chromosome has a more pronounced selection response than the autosomes, which may be attributed to dominance effects. Furthermore, we find that the median selection coefficient for all chromosomes is higher in our two-genotype experiment than in classic E&R studies. Since two random genomes harbor sufficient variation for adaptive responses, we propose that this approach is particularly well-suited for the analysis of polygenic adaptation.

Description: Conservation genomics has made dramatic improvements over the past decade, leveraging the power of genomes to infer diverse parameters central to conservation management questions. However, much of this effort has focused upon vertebrate species, despite insects providing similar flagship status with the added benefit of smaller genomes, shorter generation times and extensive historical collections in museums. Here we present the genome of the Apollo butterfly (Parnassius apollo, Papilionidae), an iconic endangered butterfly, which like many species in this genus, needs conservation genomic attention yet lacks a genome. Using 68.7 Gb of long-read data (N50 = 15.2 kb) we assembled a 1.4 Gb genome for the Apollo butterfly, making this the largest sequenced Lepidopteran genome to date. The assembly was highly contiguous (N50 = 7.1 Mb) and complete (97% of Lepidopteran BUSCOs were single-copy and complete) and consisted of 1,707 contigs. Using RNAseq data and Arthropoda proteins, we annotated 28.3K genes. Alignment with the closest-related chromosome-level assembly, Papilio bianor, reveals a highly conserved chromosomal organization, albeit genome size is highly expanded in the Apollo butterfly, due primarily to a dramatic increase in repetitive element content. Using this alignment for superscaffolding places the P. apollo genome in to 31 chromosomal scaffolds, and together with our functional annotation, provides an essential resource for advancing conservation genomics in a flagship species for insect conservation.

Description: The green rice leafhopper, Nephotettix cincticeps (Uhler), is an important rice pest and a vector of the rice dwarf virus in Asia. Here, we produced a high-quality chromosome-level genome assembly of 753.23 Mb using PacBio (∼110×) and Hi-C data (∼94×). It contained 163 scaffolds and 950 contigs, whose scaffold/contig N50 lengths reached 85.36/2.57 Mb. And 731.19 Mb (97.07%) of the assembly was anchored into eight pseudo-chromosomes. Genome completeness was attained to 97.0% according to the insect reference Benchmarking Universal Single-Copy Orthologs (BUSCO) gene set (n = 1,367). We masked 347.10 Mb (46.08%) of the genome as repetitive elements. Nine hundred sixty-two non-coding RNAs were identified and 14,337 protein-coding genes were predicted. We also assigned GO term and KEGG pathway annotations for 10,049 and 9,251 genes, respectively. Significantly expanded gene families were primarily involved in immunity, cuticle, digestion, detoxification, and embryonic development. This study provided a crucial genomic resource for better understanding on the biology and evolution in family Cicadellidae.

Description: Museum collections contain enormous quantities of insect specimens collected over the past century, covering a period of increased and varied insecticide usage. These historic collections are therefore incredibly valuable as genomic snapshots of organisms before, during, and after exposure to novel selective pressures. However, these samples come with their own challenges compared with present-day collections, as they are fragile and retrievable DNA is low yield and fragmented. In this article, we tested several DNA extraction procedures across pinned historic Diptera specimens from four disease vector genera: Anopheles, Aedes, Culex, and Glossina. We identify an approach that minimizes morphological damage while maximizing DNA retrieval for Illumina library preparation and sequencing that can accommodate the fragmented and low yield nature of historic DNA. We identify several key points in retrieving sufficient DNA while keeping morphological damage to a minimum: an initial rehydration step, a short incubation without agitation in a modified low salt Proteinase K buffer (referred to as “lysis buffer C” throughout), and critical point drying of samples post-extraction to prevent tissue collapse caused by air drying. The suggested method presented here provides a solid foundation for exploring the genomes and morphology of historic Diptera collections.

Description: The contribution of gene duplications to the evolution of eukaryotic genomes is well studied. By contrast, studies of gene duplications in prokaryotes are scarce and generally limited to a handful of genes or careful analysis of a few prokaryotic lineages. Systematic broad-scale studies of prokaryotic genomes that sample available data are lacking, leaving gaps in our understanding of the contribution of gene duplications as a source of genetic novelty in the prokaryotic world. Here, we report conservative and robust estimates for the frequency of recent gene duplications within prokaryotic genomes relative to recent lateral gene transfer (LGT), as mechanisms to generate multiple copies of related sequences in the same genome. We obtain our estimates by focusing on evolutionarily recent events among 5,655 prokaryotic genomes, thereby avoiding vagaries of deep phylogenetic inference and confounding effects of ancient events and differential loss. We find that recent, genome-specific gene duplications are at least 50 times less frequent and probably 100 times less frequent than recent, genome-specific, gene acquisitions via LGT. The frequency of gene duplications varies across lineages and functional categories. The findings improve our understanding of genome evolution in prokaryotes and have far-reaching implications for evolutionary models that entail LGT to gene duplications ratio as a parameter.

Description: Members of the gram-positive bacterial genus Staphylococcus have historically been classified into coagulase-positive Staphylococcus (CoPS) and coagulase-negative Staphylococcus (CoNS) based on the diagnostic presentation of the coagulase protein. Previous studies have noted the importance of horizontal gene transfer (HGT) and recombination in the more well-known CoPS species Staphylococcus aureus, yet little is known of the contributions of these processes in CoNS evolution. In this study, we aimed to elucidate the phylogenetic relationships, genomic characteristics, and frequencies of HGT in CoNS, which are now being recognized as major opportunistic pathogens of humans. We compiled a data set of 1,876 publicly available named CoNS genomes. These can be delineated into 55 species based on allele differences in 462 core genes and variation in accessory gene content. CoNS species are a reservoir of transferrable genes associated with resistance to diverse classes of antimicrobials. We also identified nine types of the mobile genetic element SCCmec, which carries the methicillin resistance determinant mecA. Other frequently transferred genes included those associated with resistance to heavy metals, surface-associated proteins related to virulence and biofilm formation, type VII secretion system, iron capture, recombination, and metabolic enzymes. The highest frequencies of receipt and donation of recombined DNA fragments were observed in Staphylococcus capitis, Staphylococcus caprae, Staphylococcus hominis, Staphylococcus haemolyticus, and members of the Saprophyticus species group. The variable rates of recombination and biases in transfer partners imply that certain CoNS species function as hubs of gene flow and major reservoir of genetic diversity for the entire genus.

Description: The reed warbler (Acrocephalus scirpaceus) is a long-distance migrant passerine with a wide distribution across Eurasia. This species has fascinated researchers for decades, especially its role as host of a brood parasite, and its capacity for rapid phenotypic change in the face of climate change. Currently, it is expanding its range northwards in Europe, and is altering its migratory behavior in certain areas. Thus, there is great potential to discover signs of recent evolution and its impact on the genomic composition of the reed warbler. Here, we present a high-quality reference genome for the reed warbler, based on PacBio, 10×, and Hi-C sequencing. The genome has an assembly size of 1,075,083,815 bp with a scaffold N50 of 74,438,198 bp and a contig N50 of 12,742,779 bp. BUSCO analysis using aves_odb10 as a model showed that 95.7% of BUSCO genes were complete. We found unequivocal evidence of two separate macrochromosomal fusions in the reed warbler genome, in addition to the previously identified fusion between chromosome Z and a part of chromosome 4A in the Sylvioidea superfamily. We annotated 14,645 protein-coding genes, and a BUSCO analysis of the protein sequences indicated 97.5% completeness. This reference genome will serve as an important resource, and will provide new insights into the genomic effects of evolutionary drivers such as coevolution, range expansion, and adaptations to climate change, as well as chromosomal rearrangements in birds.

Description: The Himalayan giant honeybee, Apis laboriosa, is the largest individual honeybee with major ecological and economic importance in high-latitude environments. However, our understanding of its environmental adaptations is circumscribed by the paucity of genomic data for this species. Here, we provide a draft genome of wild A. laboriosa, along with a comparison to its closely related species, Apis dorsata. The draft genome of A. laboriosa based on the de novo assembly is 226.1 Mbp in length with a scaffold N50 size of 3.34 Mbp, a GC content of 32.2%, a repeat content of 6.86%, and a gene family number of 8,404. Comparative genomics analysis revealed that the genes in A. laboriosa genome have undergone stronger positive selection (2.5 times more genes) and more recent duplication/loss events (6.1 times more events) than those in the A. dorsata genome. Our study implies the potential molecular mechanisms underlying the high-altitude adaptation of A. laboriosa and will catalyze future comparative studies to understand the environmental adaptation of modern honeybees.

Description: The bluntnose knifefish Brachyhypopomus occidentalis is a primary freshwater fish from north-western South America and Lower Central America. Like other Gymnotiformes, it has an electric organ that generates electric discharges used for both communication and electrolocation. We assembled a high-quality reference genome sequence of B. occidentalis by combining Oxford Nanopore and 10X Genomics linked-reads technologies. We also describe its demographic history in the context of the rise of the Isthmus of Panama. The size of the assembled genome is 540.3 Mb with an N50 scaffold length of 5.4 Mb, which includes 93.8% complete, 0.7% fragmented, and 5.5% of missing vertebrate/Actinoterigie Benchmarking Universal Single-Copy Orthologs. Repetitive elements account for 11.04% of the genome, and 34,347 protein-coding genes were predicted, of which 23,935 have been functionally annotated. Demographic analysis suggests a rapid effective population expansion between 3 and 5 Myr, corresponding to the final closure of the Isthmus of Panama (2.8–3.5 Myr). This event was followed by a sudden and constant population decline during the last 1 Myr, likely associated with strong shifts in both precipitation and sea level during the Pleistocene glacial-interglacial cycles. The de novo genome assembly of B. occidentalis will provide novel insights into the molecular basis of both electric signal productions and detection and will be fundamental for understanding the processes that have shaped the diversity of Neotropical freshwater environments.

Description: Genome sizes of eukaryotic organisms vary substantially, with whole-genome duplications (WGD) and transposable element expansion acting as main drivers for rapid genome size increase. The two North American mudminnows, Umbra limi and Umbra pygmaea, feature genomes about twice the size of their sister lineage Esocidae (e.g., pikes and pickerels). However, it is unknown whether all Umbra species share this genome expansion and which causal mechanisms drive this expansion. Using flow cytometry, we find that the genome of the European mudminnow is expanded similarly to both North American species, ranging between 4.5 and 5.4 pg per diploid nucleus. Observed blocks of interstitially located telomeric repeats in U. limi suggest frequent Robertsonian rearrangements in its history. Comparative analyses of transcriptome and genome assemblies show that the genome expansion in Umbra is driven by the expansion of DNA transposon and unclassified repeat sequences without WGD. Furthermore, we find a substantial ongoing expansion of repeat sequences in the Alaska blackfish Dallia pectoralis, the closest relative to the family Umbridae, which might mark the beginning of a similar genome expansion. Our study suggests that the genome expansion in mudminnows, driven mainly by transposon expansion, but not WGD, occurred before the separation into the American and European lineage.

Description: We provide a new, annotated genome assembly of Neomicropteryx cornuta, a species of the so-called mandibulate archaic moths (Lepidoptera: Micropterigidae). These moths belong to a lineage that is thought to have split from all other Lepidoptera more than 300 Ma and are consequently vital to understanding the early evolution of superorder Amphiesmenoptera, which contains the order Lepidoptera (butterflies and moths) and its sister order Trichoptera (caddisflies). Using PacBio HiFi sequencing reads, we assembled a highly contiguous genome with a contig N50 of nearly 17 Mb. The assembled genome length of 541,115,538 bp is about half the length of the largest published Amphiesmenoptera genome (Limnephilus lunatus, Trichoptera) and double the length of the smallest (Papilio polytes, Lepidoptera). We find high recovery of universal single copy orthologs with 98.1% of BUSCO genes present and provide a genome annotation of 15,643 genes aided by resolved isoforms from PacBio IsoSeq data. This high-quality genome assembly provides an important resource for studying ecological and evolutionary transitions in the early evolution of Amphiesmenoptera.

Description: Although sex is now accepted as a ubiquitous and ancestral feature of eukaryotes, direct observation of sex is still lacking in most unicellular eukaryotic lineages. Evidence of sex is frequently indirect and inferred from the identification of genes involved in meiosis from whole genome data and/or the detection of recombination signatures from genetic diversity in natural populations. In haploid unicellular eukaryotes, sex-related chromosomes are named mating-type (MTs) chromosomes and generally carry large genomic regions where recombination is suppressed. These regions have been characterized in Fungi and Chlorophyta and determine gamete compatibility and fusion. Two candidate MT+ and MT− alleles, spanning 450–650 kb, have recently been described in Ostreococcus tauri, a marine phytoplanktonic alga from the Mamiellophyceae class, an early diverging branch in the green lineage. Here, we investigate the architecture and evolution of these candidate MT+ and MT− alleles. We analyzed the phylogenetic profile and GC content of MT gene families in eight different genomes whose divergence has been previously estimated at up to 640 Myr, and found evidence that the divergence of the two MT alleles predates speciation in the Ostreococcus genus. Phylogenetic profiles of MT trans-specific polymorphisms in gametologs disclosed candidate MTs in two additional species, and possibly a third. These Mamiellales MT candidates are likely to be the oldest mating-type loci described to date, which makes them fascinating models to investigate the evolutionary mechanisms of haploid sex determination in eukaryotes.

Description: The invasive Japanese stiltgrass (Microstegium vimineum) affects a wide range of ecosystems and threatens biodiversity across the eastern USA. However, the mechanisms underlying rapid adaptation, plasticity, and epigenetics in the invasive range are largely unknown. We present a chromosome-level assembly for M. vimineum to investigate genome dynamics, evolution, adaptation, and the genomics of phenotypic plasticity. We generated a 1.12 Gb genome with scaffold N50 length of 53.44 Mb respectively, taking a de novo assembly approach that combined PacBio and Dovetail Genomics Omni-C sequencing. The assembly contains 23 pseudochromosomes, representing 99.96% of the genome. BUSCO assessment indicated that 80.3% of Poales gene groups are present in the assembly. The genome is predicted to contain 39,604 protein-coding genes, of which 26,288 are functionally annotated. Furthermore, 66.68% of the genome is repetitive, of which unclassified (35.63%) and long terminal repeat (LTR) retrotransposons (26.90%) are predominant. Similar to other grasses, Gypsy (41.07%) and Copia (32%) are the most abundant LTR-retrotransposon families. The majority of LTR-retrotransposons are derived from a significant expansion in the past 1-2 million years, suggesting the presence of relatively young LTR-retrotransposon lineages. We find corroborating evidence from Ks plots for a stiltgrass-specific duplication event, distinct from the more ancient grass-specific duplication event. The assembly and annotation of M. vimineum will serve as an essential genomic resource facilitating studies of the invasion process, the history and consequences of polyploidy in grasses, and provides a crucial tool for natural resource managers.

Description: As humans populated the world, they adapted to many varying environmental factors, including climate, diet, and pathogens. Because many of these adaptations were mediated by multiple non-coding variants with small effects on gene regulation, it has been difficult to link genomic signals of selection to specific genes, and to describe the regulatory response to selection. To overcome this challenge, we adapted PrediXcan, a machine learning method for imputing gene regulation from genotype data, to analyze low-coverage ancient human DNA (aDNA). First, we used simulated genomes to benchmark strategies for adapting PrediXcan to increase robustness to incomplete data. Applying the resulting models to 490 ancient Eurasians, we found that genes with the strongest divergent regulation among ancient populations with hunter-gatherer, pastoralist, and agricultural lifestyles are enriched for metabolic and immune functions. Next, we explored the contribution of divergent gene regulation to two traits with strong evidence of recent adaptation: dietary metabolism and skin pigmentation. We found enrichment for divergent regulation among genes proposed to be involved in diet-related local adaptation, and the predicted effects on regulation often suggest explanations for known signals of selection, e.g. at FADS1, GPX1, and LEPR. In contrast, skin pigmentation genes show little regulatory change over a 38,000-year time series of 2999 ancient Europeans, suggesting that adaptation mainly involved large-effect coding variants. This work demonstrates that combining aDNA with present-day genomes is informative about the biological differences among ancient populations, the role of gene regulation in adaptation, and the relationship between genetic diversity and complex traits.

Description: Birds in the clade Palaeognathae, excluding Tinamiformes, have morphologically conserved karyotypes and less differentiated ZW sex chromosomes compared with those of other birds. In particular, the sex chromosomes of the ostrich and emu have exceptionally large recombining pseudoautosomal regions (PARs), while non-PARs are classified into two strata according to the date of their origins: stratum 0 (S0) and stratum 1 (S1). However, the construction and analysis of the genome sequences in these regions in the clade Palaeognathae can be challenging because assembling the S1 region is difficult owing to low sequence diversity between gametologs (Z-linked and W-linked sequences). We addressed this issue by applying the Platanus-allee assembler and successfully constructed the haplotype-resolved (phased) assembly for female emu, cassowary, and ostrich using only sequence read data derived from the Illumina platform. Comparative genomic and phylogenetic analyses based on assembled Z-linked and W-linked sequences confirmed that the S1 region of emu and cassowary formed in their common ancestor. Moreover, the interspersed repetitive sequence landscapes in the S1 regions of female emu showed an expansion of younger repetitive elements in the W-linked S1 region, suggesting an interruption in homologous recombination in the S1 region. These results provide novel insights into the trajectory of sex chromosome evolution in the clade Palaeognathae and suggest that the Illumina-based phased assembly method is an effective approach for elucidating the evolutionary process underlying the transition from homomorphic to differentiated sex chromosomes.

Description: The true bug, Aspongopus chinensis Dallas, 1851 (Hemiptera: Dinidoridae), is a fascinating insect with prolonged diapause and medicinal properties but also a notorious pest. However, because of the lack of genomic resources, an in-depth understanding of its biological characteristics is lacking. Here, we report the first genome assembly of A. chinensis anchored to 10 pseudochromosomes, which was achieved by combining PacBio long reads and Hi-C sequencing data. This chromosome-level genome assembly was 1.55 Gb in size with a scaffold N50 of 156 Mb. The benchmarking universal single-copy ortholog (BUSCO) analysis of the assembly captured 96.6% of the BUSCO genes. A total of 686,888,052 bp of repeat sequences, 18,511 protein-coding genes, and 1,749 noncoding RNAs were annotated. By comparing the A. chinensis genome with that of 8 homologous insects and 2 model organisms, 213 rapidly evolving gene families were identified, including 83 expanded and 130 contracted gene families. The functional enrichment of Gene Ontology and KEGG pathways showed that the significantly expanded gene families were primarily involved in metabolism, immunity, detoxification, and DNA/RNA replication associated with stress responses. The data reported here shed light on the ecological adaptation of A. chinensis and further expanded our understanding of true bug evolution in general.

Description: Amino acid substitution models are commonly used for phylogenetic inference, for ancestral sequence reconstruction, and for the inference of positive selection. All commonly used models explicitly assume that each site evolves independently, an assumption that is violated by both linkage and protein structural and functional constraints. We introduce two new models for amino acid substitution which incorporate linkage between sites, each based on the (population-genetic) Moran model. The first model is a generalized population process tracking arbitrarily many sites which undergo mutation, with individuals replaced according to their fitnesses. This model provides a reasonably complete framework for simulations but is numerically and analytically intractable. We also introduce a second model which includes several simplifying assumptions but for which some theoretical results can be derived. We analyze the simplified model to determine conditions where linkage is likely to have meaningful effects on sitewise substitution probabilities, as well as conditions under which the effects are likely to be negligible. These findings are an important step in the generation of tractable phylogenetic models that parameterize selective coefficients for amino acid substitution while accounting for linkage of sites leading to both hitchhiking and background selection.

Description: Advances in phylogenomics contribute toward resolving long-standing evolutionary questions. Notwithstanding, genetic diversity contained within more than a billion biological specimens deposited in natural history museums remains recalcitrant to analysis owing to challenges posed by its intrinsically degraded nature. Yet that tantalizing resource could be critical in overcoming taxon sampling constraints hindering our ability to address major evolutionary questions. We addressed this impediment by developing phyloHyRAD, a new bioinformatic pipeline enabling locus recovery at a broad evolutionary scale from HyRAD-X exome capture of museum specimens of low DNA integrity using a benchtop RAD-derived exome-complexity-reduction probe set developed from high DNA integrity specimens. Our new pipeline can also successfully align raw RNAseq transcriptomic and ultraconserved element reads with the RAD-derived probe catalog. Using this method, we generated a robust timetree for Carabinae beetles, the lack of which had precluded study of macroevolutionary trends pertaining to their biogeography and wing-morphology evolution. We successfully recovered up to 2,945 loci with a mean of 1,788 loci across the exome of specimens of varying age. Coverage was not significantly linked to specimen age, demonstrating the wide exploitability of museum specimens. We also recovered fragmentary mitogenomes compatible with Sanger-sequenced mtDNA. Our phylogenomic timetree revealed a Lower Cretaceous origin for crown group Carabinae, with the extinct Aplothorax Waterhouse, 1841 nested within the genus Calosoma Weber, 1801 demonstrating the junior synonymy of Aplothorax syn. nov., resulting in the new combination Calosomaburchellii (Waterhouse, 1841) comb. nov. This study compellingly illustrates that HyRAD-X and phyloHyRAD efficiently provide genomic-level data sets informative at deep evolutionary scales.

Description: Ceratonova shasta is an important myxozoan pathogen affecting the health of salmonid fishes in the Pacific Northwest of North America. Ceratonova shasta exists as a complex of host-specific genotypes, some with low to moderate virulence, and one that causes a profound, lethal infection in susceptible hosts. High throughput sequencing methods are powerful tools for discovering the genetic basis of these host/virulence differences, but deep sequencing of myxozoans has been challenging due to extremely fast molecular evolution of this group, yielding strongly divergent sequences that are difficult to identify, and unavoidable host contamination. We designed and optimized different bioinformatic pipelines to address these challenges. We obtained a unique set of comprehensive, host-free myxozoan RNA-seq data from C. shasta genotypes of varying virulence from different salmonid hosts. Analyses of transcriptome-wide genetic distances and maximum likelihood multigene phylogenies elucidated the evolutionary relationship between lineages and demonstrated the limited resolution of the established Internal Transcribed Spacer marker for C. shasta genotype identification, as this marker fails to differentiate between biologically distinct genotype II lineages from coho salmon and rainbow trout. We further analyzed the data sets based on polymorphisms in two gene groups related to virulence: cell migration and proteolytic enzymes including their inhibitors. The developed single-nucleotide polymorphism-calling pipeline identified polymorphisms between genotypes and demonstrated that variations in both motility and protease genes were associated with different levels of virulence of C. shasta in its salmonid hosts. The prospective use of proteolytic enzymes as promising candidates for targeted interventions against myxozoans in aquaculture is discussed. We developed host-free transcriptomes of a myxozoan model organism from strains that exhibited different degrees of virulence, as a unique source of data that will foster functional gene analyses and serve as a base for the development of potential therapeutics for efficient control of these parasites.

Description: The fraction of the human genome that is functional is a question of both evolutionary and practical importance. Studies of sequence divergence have suggested that the functional fraction of the human genome is likely to be no more than ∼15%. In contrast, the ENCODE project, a systematic effort to map regions of transcription, transcription factor association, chromatin structure, and histone modification, assigned function to 80% of the human genome. In this article, we examine whether and how an analysis based on mutational load might set a limit on the functional fraction. In order to do so, we characterize the distribution of fitness of a large, finite, diploid population at mutation-selection equilibrium. In particular, if mean fitness is ∼1, the fitness of the fittest individual likely to occur cannot be unreasonably high. We find that at equilibrium, the distribution of log fitness has variance nus, where u is the per-base deleterious mutation rate, n is the number of functional sites (and hence incorporates the functional fraction f), and s is the selection coefficient of deleterious mutations. In a large (N=109) reproducing population, the fitness of the fittest individual likely to exist is ∼e5nus. These results apply to both additive and recessive fitness schemes. Our approach is different from previous work that compared mean fitness at mutation-selection equilibrium with the fitness of an individual who has no deleterious mutations; we show that such an individual is exceedingly unlikely to exist. We find that the functional fraction is not very likely to be limited substantially by mutational load, and that any such limit, if it exists, depends strongly on the selection coefficients of new deleterious mutations.

Description: Animals that use venom to feed on a wide diversity of prey may evolve a complex mixture of toxins to target a variety of physiological processes and prey-defense mechanisms. Blarina brevicauda, the northern short-tailed shrew, is one of few venomous mammals, and is also known to eat evolutionarily divergent prey. Despite their complex diet, earlier proteomic and transcriptomic studies of this shrew’s venom have only identified two venom proteins. Here, we investigated with comprehensive molecular approaches whether B. brevicauda venom is more complex than previously understood. We generated de novo assemblies of a B. brevicauda genome and submaxillary-gland transcriptome, as well as sequenced the salivary proteome. Our findings show that B. brevicauda’s venom composition is simple relative to their broad diet and is likely limited to seven proteins from six gene families. Additionally, we explored expression levels and rate of evolution of these venom genes and the origins of key duplications that led to toxin neofunctionalization. We also found three proteins that may be involved in endogenous self-defense. The possible synergism of the toxins suggests that vertebrate prey may be the main target of the venom. Further functional assays for all venom proteins on both vertebrate and invertebrate prey would provide further insight into the ecological relevance of venom in this species.

Description: In 1957, an invasive and highly defensive honey bee began to spread across Brazil. In the previous year, Brazilian researchers hoped to produce a subtropical-adapted honey bee by crossing local commercial honey bees (of European origin) with a South African honey bee subspecies (Apis mellifera scutellata; an A-lineage honey bee subspecies). The resulting cross—African hybrid honey bees (AHBs)—escaped from their enclosure and spread through the Americas. Today, AHB is the most common honey bee from Northern Argentina to the Southern United States. AHBs are much more likely to sting nest intruders than managed European-derived honey bee colonies. Previous studies have explored how genetic variation contributes to differences in defense response between European-derived honey bee and AHB. Although this work demonstrated very strong genetic effects on defense response, they have yet to pinpoint which genes influence variation in defense response within AHBs, specifically. We quantified defense response for 116 colonies in Brazil and performed pooled sequencing on the most phenotypically divergent samples. We identified 65 loci containing 322 genes that were significantly associated with defense response. Loci were strongly associated with metabolic function, consistent with previous functional genomic analyses of this phenotype. Additionally, defense-associated loci had nonrandom and unexpected patterns of admixture. Defense response was not simply the product of more A-lineage honey bee ancestry as previously assumed, but rather an interaction between A-lineage and European alleles. Our results suggest that a combination of A-lineage and European alleles play roles in defensive behavior in AHBs.

Description: The water skinks Eulamprus tympanum and Eulamprus heatwolei show thermally induced sex determination where elevated temperatures give rise to male offspring. Paradoxically, Eulamprus species reproduce in temperatures of 12–15 °C making them outliers when compared with reptiles that use temperature as a cue for sex determination. Moreover, these two species are among the very few viviparous reptiles reported to have thermally induced sex determination. Thus, we tested whether these skinks possess undetected sex chromosomes with thermal override. We produced transcriptome and genome data for E. heatwolei. We found that E. heatwolei presents XY chromosomes that include 14 gametologs with regulatory functions. The Y chromosomal region is 79–116 Myr old and shared between water and spotted skinks. Our work provides clear evidence that climate could be useful to predict the type of sex determination systems in reptiles and it also indicates that viviparity is strictly associated with sex chromosomes.

Description: Schmidingerella arcuata is an ecologically important tintinnid ciliate that has long served as a model species in plankton trophic ecology. We present a partial micronuclear genome and macronuclear transcriptome resource for S. arcuata, acquired using single-cell techniques, and we report on pilot analyses including functional annotation and genome architecture. Our analysis shows major fragmentation, elimination, and scrambling in the micronuclear genome of S. arcuata. This work introduces a new nonmodel genome resource for the study of ciliate ecology and genomic biology and provides a detailed functional counterpart to ecological research on S. arcuata.

Description: The solute carrier (SLC) transporter superfamily comprises an ancient and ubiquitous group of proteins capable of translocating a range of nutrients, endogenous molecules, and xenobiotics. Although the group has been the subject of intense investigation in both bacteria and mammals, its systematic identification in arthropods has not yet been undertaken. Here, we present a genome-wide identification of all 66 human SLC families in 174 arthropod species. A pipeline (SLC_id) was constructed to identify and group SLCs using a combination of hidden Markov model and BLAST searches followed by filtering based on polypeptide length and the number of transmembrane domains. Comparative analysis of the number of transporters in each family across diverse arthropod lineages was accomplished using one-way analysis of variance (ANOVA) and the Computational Analysis of gene Family Evolution (CAFE). These results suggested that many SLC families have undergone expansions or contractions in particular evolutionary lineages. Notably, the sugar transporting SLC2 family was significantly larger in insects compared with arachnids. This difference may have been complemented by a rapid expansion of the SLC60 family in arachnids which also acts on dietary sugars. Furthermore, the SLC33 family underwent a recent and drastic expansion in aphids, although the biological relevance of this expansion was not possible to infer. Information on specific SLC transporter families across arthropod species can be accessed through an R shiny web application at http://chrysalida.imbb.forth.gr : 3838/Arthropod_SLC_Database/. The present study greatly facilitates further investigation of the diverse group of SLC transporters in arthropods.

Description: Despite the extreme and varying environmental conditions prevalent in the Arabian Peninsula, it has experienced several waves of human migrations following the out-of-Africa diaspora. Eventually, the inhabitants of the peninsula region adapted to the hot and dry environment. The adaptation and natural selection that shaped the extant human populations of the Arabian Peninsula region have been scarcely studied. In an attempt to explore natural selection in the region, we analyzed 662,750 variants in 583 Kuwaiti individuals. We searched for regions in the genome that display signatures of positive selection in the Kuwaiti population using an integrative approach in a conservative manner. We highlight a haplotype overlapping TNKS that showed strong signals of positive selection based on the results of the multiple selection tests conducted (integrated Haplotype Score, Cross Population Extended Haplotype Homozygosity, Population Branch Statistics, and log-likelihood ratio scores). Notably, the TNKS haplotype under selection potentially conferred a fitness advantage to the Kuwaiti ancestors for surviving in the harsh environment while posing a major health risk to present-day Kuwaitis.

Description: The ribosome’s common core, comprised of ribosomal RNA (rRNA) and universal ribosomal proteins, connects all life back to a common ancestor and serves as a window to relationships among organisms. The rRNA of the common core is similar to rRNA of extant bacteria. In eukaryotes, the rRNA of the common core is decorated by expansion segments (ESs) that vastly increase its size. Supersized ESs have not been observed previously in Archaea, and the origin of eukaryotic ESs remains enigmatic. We discovered that the large ribosomal subunit (LSU) rRNA of two Asgard phyla, Lokiarchaeota and Heimdallarchaeota, considered to be the closest modern archaeal cell lineages to Eukarya, bridge the gap in size between prokaryotic and eukaryotic LSU rRNAs. The elongated LSU rRNAs in Lokiarchaeota and Heimdallarchaeota stem from two supersized ESs, called ES9 and ES39. We applied chemical footprinting experiments to study the structure of Lokiarchaeota ES39. Furthermore, we used covariation and sequence analysis to study the evolution of Asgard ES39s and ES9s. By defining the common eukaryotic ES39 signature fold, we found that Asgard ES39s have more and longer helices than eukaryotic ES39s. Although Asgard ES39s have sequences and structures distinct from eukaryotic ES39s, we found overall conservation of a three-way junction across the Asgard species that matches eukaryotic ES39 topology, a result consistent with the accretion model of ribosomal evolution.

Description: We examine the genetic history and population status of Hawaiian hoary bats (Lasiurus semotus), the most isolated bats on Earth, and their relationship to northern hoary bats (Lasiurus cinereus), through whole-genome analysis of single-nucleotide polymorphisms mapped to a de novo-assembled reference genome. Profiles of genomic diversity and divergence indicate that Hawaiian hoary bats are distinct from northern hoary bats, and form a monophyletic group, indicating a single ancestral colonization event 1.34 Ma, followed by substantial divergence between islands beginning 0.51 Ma. Phylogenetic analysis indicates Maui is central to the radiation across the archipelago, with the southward expansion to Hawai‘i and westward to O‘ahu and Kaua‘i. Because this endangered species is of conservation concern, a clearer understanding of the population genetic structure of this bat in the Hawaiian Islands is of timely importance.

Description: Evolve and resequencing (E&R) studies investigate the genomic responses of adaptation during experimental evolution. Because replicate populations evolve in the same controlled environment, consistent responses to selection across replicates are frequently used to identify reliable candidate regions that underlie adaptation to a new environment. However, recent work demonstrated that selection signatures can be restricted to one or a few replicate(s) only. These selection signatures frequently have weak statistical support, and given the difficulties of functional validation, additional evidence is needed before considering them as candidates for functional analysis. Here, we introduce an experimental procedure to validate candidate loci with weak or replicate-specific selection signature(s). Crossing an evolved population from a primary E&R experiment to the ancestral founder population reduces the frequency of candidate alleles that have reached a high frequency. We hypothesize that genuine selection targets will experience a repeatable frequency increase after the mixing with the ancestral founders if they are exposed to the same environment (secondary E&R experiment). Using this approach, we successfully validate two overlapping selection targets, which showed a mutually exclusive selection signature in a primary E&R experiment of Drosophila simulans adapting to a novel temperature regime. We conclude that secondary E&R experiments provide a reliable confirmation of selection signatures that either are not replicated or show only a low statistical significance in a primary E&R experiment unless epistatic interactions predominate. Such experiments are particularly helpful to prioritize candidate loci for time-consuming functional follow-up investigations.

Description: The application of Wolbachia in insect pest and vector control requires the establishment of genotypically stable host associations. The cytoplasmic incompatibility (CI) inducing Wolbachia strain wCer2 naturally occurs in the cherry fruit fly Rhagoletis cerasi as co-infection with other strains and was transferred to other fruit fly species by embryonic microinjections. We obtained wCer2 genome data from its native and three novel hosts, Drosophila simulans, Drosophila melanogaster, and Ceratitis capitata and assessed its genome stability, characteristics, and CI factor (cif) genes. De novo assembly was successful from Wolbachia cell-enriched singly infected D. simulans embryos, with minimal host and other bacterial genome traces. The low yield of Wolbachia sequence reads from total genomic extracts of one multiply infected R. cerasi pupa and one singly infected C. capitata adult limited de novo assemblies but was sufficient for comparative analyses. Across hosts wCer2 was stable in genome synteny and content. Polymorphic nucleotide sites were found in wCer2 of each host; however, only one nucleotide was different between R. cerasi and C. capitata, and none between replicated D. simulans lines. The wCer2 genome is highly similar to wAu (D. simulans), wMel (D. melanogaster), and wRec (Drosophila recens). In contrast to wMel and wRec (each with one cif gene pair) and wAu (without any cif genes), wCer2 has three pairs of Type I cif genes, and one Type V cifB gene without a cifA complement. This may explain previously reported CI patterns of wCer2, including incomplete rescue of its own CI modification in three novel host species.

Description: Corals are the ecosystem engineers of coral reefs, one of the most biodiverse marine ecosystems. The ability of corals to form reefs depends on the precipitation of calcium carbonate (CaCO3) under biological control. However, several mechanisms underlying coral biomineralization remain elusive, for example, whether corals employ different molecular machineries to deposit different CaCO3 polymorphs (i.e., aragonite or calcite). Here, we used tandem mass spectrometry (MS/MS) to compare the proteins occluded in the skeleton of three octocoral and one scleractinian species: Tubipora musica and Sinularia cf. cruciata (calcite sclerites), the blue coral Heliopora coerulea (aragonitic skeleton), and the scleractinian aragonitic Montipora digitata. Reciprocal Blast analysis revealed extremely low overlap between aragonitic and calcitic species, while a core set of proteins is shared between octocorals producing calcite sclerites. However, the carbonic anhydrase CruCA4 is present in the skeletons of both polymorphs. Phylogenetic analysis highlighted several possible instances of protein co-option in octocorals. These include acidic proteins and scleritin, which appear to have been secondarily recruited for calcification and likely derive from proteins playing different functions. Similarities between octocorals and scleractinians included presence of a galaxin-related protein, carbonic anhydrases, and one hephaestin-like protein. Although the first two appear to have been independently recruited, the third appear to share a common origin. This work represents the first attempt to identify and compare proteins associated with coral skeleton polymorph diversity, providing several new research targets and enabling both future functional and evolutionary studies aimed at elucidating the origin and evolution of coral biomineralization.

Description: Human-specific pseudogenization of the CMAH gene eliminated the mammalian sialic acid (Sia) Neu5Gc (generating an excess of its precursor Neu5Ac), thus changing ubiquitous cell surface “self-associated molecular patterns” that modulate innate immunity via engagement of CD33-related-Siglec receptors. The Alu-fusion-mediated loss-of-function of CMAH fixed ∼2–3 Ma, possibly contributing to the origins of the genus Homo. The mutation likely altered human self-associated molecular patterns, triggering multiple events, including emergence of human-adapted pathogens with strong preference for Neu5Ac recognition and/or presenting Neu5Ac-containing molecular mimics of human glycans, which can suppress immune responses via CD33-related-Siglec engagement. Human-specific alterations reported in some gene-encoding Sia-sensing proteins suggested a “hotspot” in hominin evolution. The availability of more hominid genomes including those of two extinct hominins now allows full reanalysis and evolutionary timing. Functional changes occur in 8/13 members of the human genomic cluster encoding CD33-related Siglecs, all predating the human common ancestor. Comparisons with great ape genomes indicate that these changes are unique to hominins. We found no evidence for strong selection after the Human–Neanderthal/Denisovan common ancestor, and these extinct hominin genomes include almost all major changes found in humans, indicating that these changes in hominin sialobiology predate the Neanderthal–human divergence ∼0.6 Ma. Multiple changes in this genomic cluster may also explain human-specific expression of CD33rSiglecs in unexpected locations such as amnion, placental trophoblast, pancreatic islets, ovarian fibroblasts, microglia, Natural Killer(NK) cells, and epithelia. Taken together, our data suggest that innate immune interactions with pathogens markedly altered hominin Siglec biology between 0.6 and 2 Ma, potentially affecting human evolution.

Description: Transposable elements (TEs) are repetitive regions of DNA that are able to self-replicate and reinsert themselves throughout host genomes. Since the discovery of TEs, a prevalent question has been whether increasing TE copy number has an effect on the fitness of their hosts. P-elements (PEs) in Drosophila are a well-studied TE that has strong phenotypic effects. When a female without PEs (M) is crossed to a male with them (P), the resulting females are often sterile, a phenomenon called hybrid dysgenesis (HD). Here, we used short- and long-read sequencing to infer the number of PEs in the genomes of dozens of isofemale lines from two Drosophila species and measured whether the magnitude of HD was correlated with the number of PEs in the paternal genome. Consistent with previous reports, we find evidence for a positive correlation between the paternal PE copy number and the magnitude of HD in progeny from ♀M × ♂ P crosses for both species. Other crosses are not affected by the number of PE copies. We also find that the correlation between the strength of HD and PE copy number differs between species, which suggests that there are genetic differences that might make some genomes more resilient to the potentially deleterious effects of TEs. Our results suggest that PE copy number interacts with other factors in the genome and the environment to cause HD and that the importance of these interactions is species specific.

Description: Recombination confers a major evolutionary advantage by breaking up linkage disequilibrium between harmful and beneficial mutations, thereby facilitating selection. However, in species that are only periodically sexual, such as many microbial eukaryotes, the realized rate of recombination is also affected by the frequency of sex, meaning that infrequent sex can increase the effects of selection at linked sites despite high recombination rates. Despite this, the rate of sex of most facultatively sexual species is unknown. Here, we use genomewide patterns of linkage disequilibrium to infer fine-scale recombination rate variation in the genome of the facultatively sexual green alga Chlamydomonas reinhardtii. We observe recombination rate variation of up to two orders of magnitude and find evidence of recombination hotspots across the genome. Recombination rate is highest flanking genes, consistent with trends observed in other nonmammalian organisms, though intergenic recombination rates vary by intergenic tract length. We also find a positive relationship between nucleotide diversity and physical recombination rate, suggesting a widespread influence of selection at linked sites in the genome. Finally, we use estimates of the effective rate of recombination to calculate the rate of sex that occurs in natural populations, estimating a sexual cycle roughly every 840 generations. We argue that the relatively infrequent rate of sex and large effective population size creates a population genetic environment that increases the influence of selection on linked sites across the genome.

Description: Members of the Roseobacter group are known for their different ecologically relevant metabolic traits and high abundance in many marine environments. This includes traits like carbon monoxide oxidation, sulfur oxidation, nitrogen oxidation, DMSP demethylation, denitrification, and production of bioactive compounds. Nevertheless, their role in the marine biogeochemical cycles remains to be elucidated. Roseobacter ponti DSM 106830T, also designated strain MM-7T (=KCTC 52469T =NBRC 112431T), is a novel type strain of the Roseobacter group, which was proposed as new Roseobacter species. It was isolated from seawater of the Yellow Sea in South Korea. We report the complete genome sequence of R. ponti DSM 106830T, which belongs to the family Rhodobacteraceae. The genome of R. ponti DSM 106830T comprises a single circular chromosome (3,861,689 bp) with a GC content of 60.52% and an additional circular plasmid (p1) of 100,942 bp with a GC content of 61.51%. The genome encodes 3,812 putative genes, including 3 rRNA, 42 tRNA, 1 tmRNA, and 3 ncRNA. The genome information was used to perform a phylogenetic analysis, which confirmed that the strain represents a new species. Moreover, the genome sequence enabled the investigation of the metabolic capabilities and versatility of R. ponti DSM 106830T. Finally, it provided insight into the high niche adaptation potential of Roseobacter group members.

Description: Rhizobia are soil bacteria capable of forming symbiotic nitrogen-fixing nodules associated with leguminous plants. In fast-growing legume-nodulating rhizobia, such as the species in the family Rhizobiaceae, the symbiotic plasmid is the main genetic basis for nitrogen-fixing symbiosis, and is susceptible to horizontal gene transfer. To further understand the symbioses evolution in Rhizobiaceae, we analysed the pan-genome of this family based on 92 genomes of type/reference strains and reconstructed its phylogeny using a phylogenomics approach. Intriguingly, although the genetic expansion that occurred in chromosomal regions was the main reason for the high proportion of low-frequency flexible gene families in the pan-genome, gene gain events associated with accessory plasmids introduced more genes into the genomes of nitrogen-fixing species. For symbiotic plasmids, although horizontal gene transfer frequently occurred, transfer may be impeded by, such as, the host’s physical isolation and soil conditions, even among phylogenetically close species. During coevolution with leguminous hosts, the plasmid system, including accessory and symbiotic plasmids, may have evolved over a time span, and provided rhizobial species with the ability to adapt to various environmental conditions and helped them achieve nitrogen fixation. These findings provide new insights into the phylogeny of Rhizobiaceae and advance our understanding of the evolution of symbiotic nitrogen fixation.

Description: We previously introduced a numerical quantity called the stability (Ps) of an inferred tree and showed that for the tree to be reliable this stability as well as the reliability of the tree, which is usually computed as the bootstrap probability (Pb), must be high. However, if genome duplication occurs in a species, a gene family of the genome also duplicates, and for this reason alone some Ps values can be high in a tree of the duplicated gene families. In addition, the topology of the duplicated gene family can be similar to that of the original gene family if such gene families are identifiable. After genome duplication, however, the gene families are often partially deleted or partially duplicated, and the duplicated gene family may not show the same topology as that of the original family. It is therefore necessary to compute the similarity of the topologies of the duplicated and the original gene families. In this paper, we introduce another quantity called the reproducibility (Pr) for measuring the similarity of the two gene families. To show how to compute the Pr values, we first compute the Pb and Ps values for each of the MHC class II α and β chain gene families, which were apparently generated by genome duplication. We then compute the Pr values for the MHC class II α and β chain gene families. The Pr values for the α and β chain gene families are now low, and this suggests that the diploidization of gene segregation has occurred after the genome duplication. Currently higher animals, defined as animals with complex phenotypic characters, generally have a higher genome size, and this increase in genome size appears to have been caused by genome duplication and diploidization of gene segregation after genome duplication.

Description: Salmonella enterica serovar Infantis (S. Infantis) is one of the dominant serovars of the bacterial pathogen S. enterica. In recent years, the number of human infections caused by S. Infantis has been increasing in many countries, and often the emerging population harbors a unique virulence-resistant megaplasmid called plasmid of emerging S. Infantis (pESI). Here, we report the complete gap-free genome sequence of the S. Infantis Israeli emerging clone and compare its chromosome and pESI sequences with other complete S. Infantis genomes. We show a conserved presence of the Salmonella pathogenicity islands 1–6, 9, 11, 12, and CS54 and a common integration of five bacteriophages in the S. Infantis chromosome. In contrast, we found variable presence of additionally three chromosomally integrated phages and eight modular regions in pESI, which contribute to the genetic and phenotypic diversity (including antimicrobial resistance) of this ubiquitous foodborne pathogen.

Description: This study uses population genomic data to estimate demographic and selection parameters in two sister lineages of the wild yeast Saccharomyces paradoxus and compare their evolution. We first estimate nucleotide and recombinational diversities in each of the two lineages to infer their population size and frequency of sex and then analyze the rate of mutation accumulation since divergence from their inferred common ancestor to estimate the generation time and efficacy of selection. We find that one of the lineages has significantly higher silent nucleotide diversity and lower linkage disequilibrium, indicating a larger population with more frequent sexual generations. The same lineage also shows shorter generation time and higher efficacy of purifying selection, the latter consistent with the finding of larger population size and more frequent sex. Similar analyses are also performed on the ancestries of individual strains within lineages and we find significant differences between strains implying variation in rates of mitotic cell divisions. Our sample includes some strains originating in the Chernobyl nuclear-accident exclusion zone, which has been subjected to high levels of radiation for nearly 30 years now. We find no evidence, however, for increased rates of mutation. Finally, there is a positive correlation between rates of mutation accumulation and length of growing period, as measured by latitude of the place of origin of strains. Our study illustrates the power of genomic analyses in estimating population and life history parameters and testing predictions based on population genetic theory.

Description: Pseudomonas aeruginosa filamentous (Pf) bacteriophages are important factors contributing to the pathogenicity of this opportunistic bacterium, including biofilm formation and suppression of bacterial phagocytosis by macrophages. In addition, the capacity of Pf phages to form liquid crystal structures and their high negative charge density makes them potent sequesters of cationic antibacterial agents, such as aminoglycoside antibiotics or host antimicrobial peptides. Therefore, Pf phages have been proposed as a potential biomarker for risk of antibiotic resistance development. The majority of studies describing biological functions of Pf viruses have been performed with only three of them: Pf1, Pf4, and Pf5. However, our analysis revealed that Pf phages exist as two evolutionary lineages (I and II), characterized by substantially different structural/morphogenesis properties, despite sharing the same integration sites in the host chromosomes. All aforementioned model Pf phages are members of the lineage I. Hence, it is reasonable to speculate that their interactions with P. aeruginosa and impact on its pathogenicity may be not completely extrapolated to the lineage II members. Furthermore, in order to organize the present numerical nomenclature of Pf phages, we propose a more informative approach based on the insertion sites, that is, Pf-tRNA-Gly, -Met, -Sec, -tmRNA, and -DR (direct repeats), which are fully compatible with one of five types of tyrosine integrases/recombinases XerC/D carried by these viruses. Finally, we discuss possible evolutionary mechanisms behind this division and consequences from the perspective of virus–virus, virus–bacterium, and virus–human interactions.

Description: Apoptosis is a fundamental feature of multicellular animals and is best understood in mammals, flies, and nematodes, with the invertebrate models being thought to represent a condition of ancestral simplicity. However, the existence of a leukemia-like cancer in the softshell clam Mya arenaria provides an opportunity to re-evaluate the evolution of the genetic machinery of apoptosis. Here, we report the whole-genome sequence for M. arenaria which we leverage with existing data to test evolutionary hypotheses on the origins of apoptosis in animals. We show that the ancestral bilaterian p53 locus, a master regulator of apoptosis, possessed a complex domain structure, in contrast to that of extant ecdysozoan p53s. Further, ecdysozoan taxa, but not chordates or lophotrochozoans like M. arenaria, show a widespread reduction in apoptosis gene copy number. Finally, phylogenetic exploration of apoptosis gene copy number reveals a striking linkage with p53 domain complexity across species. Our results challenge the current understanding of the evolution of apoptosis and highlight the ancestral complexity of the bilaterian apoptotic tool kit and its subsequent dismantlement during the ecdysozoan radiation.

Description: Cold seeps, characterized by the methane, hydrogen sulfide, and other hydrocarbon chemicals, foster one of the most widespread chemosynthetic ecosystems in deep sea that are densely populated by specialized benthos. However, scarce genomic resources severely limit our knowledge about the origin and adaptation of life in this unique ecosystem. Here, we present a genome of a deep-sea limpet Bathyacmaea lactea, a common species associated with the dominant mussel beds in cold seeps. We yielded 54.6 gigabases (Gb) of Nanopore reads and 77.9-Gb BGI-seq raw reads, respectively. Assembly harvested a 754.3-Mb genome for B. lactea, with 3,720 contigs and a contig N50 of 1.57 Mb, covering 94.3% of metazoan Benchmarking Universal Single-Copy Orthologs. In total, 23,574 protein-coding genes and 463.4 Mb of repetitive elements were identified. We analyzed the phylogenetic position, substitution rate, demographic history, and TE activity of B. lactea. We also identified 80 expanded gene families and 87 rapidly evolving Gene Ontology categories in the B. lactea genome. Many of these genes were associated with heterocyclic compound metabolism, membrane-bounded organelle, metal ion binding, and nitrogen and phosphorus metabolism. The high-quality assembly and in-depth characterization suggest the B. lactea genome will serve as an essential resource for understanding the origin and adaptation of life in the cold seeps.

Description: Nematodes are highly abundant animals with diverse habitats and lifestyles. Some are free living whereas others parasitize animals or plants, and among the latter, infection abilities change across developmental stages to infect hosts and complete life cycles. To determine the relationship between transcriptome evolution and morphological divergences among nematodes, we compared 48 transcriptomes of different developmental stages across eight nematode species. The transcriptomes were clustered broadly into embryo, larva, and adult stages, with the developmental plastic stages were separated from common larval stages within the larval branch. This suggests that development was the major determining factor after lifestyle changes, such as parasitism, during transcriptome evolution. Such patterns were partly accounted for by tissue-specific genes—such as those in oocytes and the hypodermis—being expressed at different proportions. Although nematodes typically have 3–5 larval stages, the transcriptomes for these stages were found to be highly correlated within each species, suggesting high similarity among larval stages across species. For the Caenorhabditis elegans–Caenorhabditis briggsae and Strongyloides stercoralis–Strongyloides venezuelensis comparisons, we found that ∼50% of genes were expressed at multiple stages, whereas half of their orthologs were also expressed in multiple but different stages. Such frequent changes in expression have resulted in concerted transcriptome evolution across adjacent stages, thus generating species-specific transcriptomes over the course of nematode evolution. Our study provides a first insight into the evolution of nematode transcriptomes beyond embryonic development.

Description: The origin of plastids (chloroplasts) by endosymbiosis stands as one of the most important events in the history of eukaryotic life. The genetic, biochemical, and cell biological integration of a cyanobacterial endosymbiont into a heterotrophic host eukaryote approximately a billion years ago paved the way for the evolution of diverse algal groups in a wide range of aquatic and, eventually, terrestrial environments. Plastids have on multiple occasions also moved horizontally from eukaryote to eukaryote by secondary and tertiary endosymbiotic events. The overall picture of extant photosynthetic diversity can best be described as “patchy”: Plastid-bearing lineages are spread far and wide across the eukaryotic tree of life, nested within heterotrophic groups. The algae do not constitute a monophyletic entity, and understanding how, and how often, plastids have moved from branch to branch on the eukaryotic tree remains one of the most fundamental unsolved problems in the field of cell evolution. In this review, we provide an overview of recent advances in our understanding of the origin and spread of plastids from the perspective of comparative genomics. Recent years have seen significant improvements in genomic sampling from photosynthetic and nonphotosynthetic lineages, both of which have added important pieces to the puzzle of plastid evolution. Comparative genomics has also allowed us to better understand how endosymbionts become organelles.

Description: Introgression events and population admixture occurred among Sus species across the Eurasian mainland in the Middle Pleistocene, which reflects the local adaption of different populations and contributes to evolutionary novelty. Previous findings on these population introgressions were largely based on extensive genome-wide single nucleotide polymorphism (SNP) information, ignoring structural variants (SVs) as an important alternative resource of genetic variations. Here, we profiled the genome-wide SVs and explored the formation of pattern-related SVs, indicating that PRE1-SS is a recently active subfamily that was strongly associated with introgression events in multiple Asian and European pig populations. As reflected by the three different combination haplotypes from two specific patterns and known phylogenetic relationships in Eurasian boars, we identified the Asian Northern wild pigs as having experienced introgression from European wild boars around 0.5 to 0.2 Mya and having received latitude-related selection. During further exploration of the influence of pattern-related SVs on gene functions, we found substantial sequence changes in 199 intron regions of 54 genes and 3 exon regions of 3 genes (HDX, TRO and SMIM1), implying that the pattern-related SVs were highly related to positive selection and adaption of pigs. Our findings revealed novel introgression events in Eurasian wild boars, providing a timeline of population admixture and divergence across the Eurasian mainland in the Middle Pleistocene.

Description: Long-term balancing selection results in a build-up of alleles at similar frequencies and a deficit of substitutions when compared with an outgroup at a locus. The previously published β(1) statistics detect balancing selection using only polymorphism data. We now propose the β(2) statistic which detects balancing selection using both polymorphism and substitution data. In addition, we derive the variance of all β statistics, allowing for their standardization and thereby reducing the influence of parameters which can confound other selection tests. The standardized β statistics outperform existing summary statistics in simulations, indicating β is a well-powered and widely applicable approach for detecting balancing selection. We apply the β(2) statistic to 1000 Genomes data and report two missense mutations with high β scores in the ACSBG2 gene. An implementation of all β statistics and their standardization are available in the BetaScan2 software package at https://github.com/ksiewert/BetaScan.

Description: The discovery of new members of the Tc1/mariner superfamily of transposons is expected based on the increasing availability of genome sequencing data. Here, we identified a new DD35E family termed Traveler (TR). Phylogenetic analyses of its DDE domain and full-length transposase showed that, although TR formed a monophyletic clade, it exhibited the highest sequence identity and closest phylogenetic relationship with DD34E/Tc1. This family displayed a very restricted taxonomic distribution in the animal kingdom and was only detected in ray-finned fish, anura, and squamata, including 91 vertebrate species. The structural organization of TRs was highly conserved across different classes of animals. Most intact TR transposons had a length of ∼1.5 kb (range 1,072–2,191 bp) and harbored a single open reading frame encoding a transposase of ∼340 aa (range 304–350 aa) flanked by two short-terminal inverted repeats (13–68 bp). Several conserved motifs, including two helix-turn-helix motifs, a GRPR motif, a nuclear localization sequence, and a DDE domain, were also identified in TR transposases. This study also demonstrated the presence of horizontal transfer events of TRs in vertebrates, whereas the average sequence identities and the evolutionary dynamics of TR elements across species and clusters strongly indicated that the TR family invaded the vertebrate lineage very recently and that some of these elements may be currently active, combining the intact TR copies in multiple lineages of vertebrates. These data will contribute to the understanding of the evolutionary history of Tc1/mariner transposons and that of their hosts.

Description: Tunicates, the closest living relatives of vertebrates, have served as a foundational model of early embryonic development for decades. Comparative studies of tunicate phylogeny and genome evolution provide a critical framework for analyzing chordate diversification and the emergence of vertebrates. Toward this goal, we sequenced the genome of Corella inflata (Ascidiacea, Phlebobranchia), so named for the capacity to brood self-fertilized embryos in a modified, “inflated” atrial chamber. Combining the new genome sequence for Co. inflata with publicly available tunicate data, we estimated a tunicate species phylogeny, reconstructed the ancestral Hox gene cluster at important nodes in the tunicate tree, and compared patterns of gene loss between Co. inflata and Ciona robusta, the prevailing tunicate model species. Our maximum-likelihood and Bayesian trees estimated from a concatenated 210-gene matrix were largely concordant and showed that Aplousobranchia was nested within a paraphyletic Phlebobranchia. We demonstrated that this relationship is not an artifact due to compositional heterogeneity, as had been suggested by previous studies. In addition, within Thaliacea, we recovered Doliolida as sister to the clade containing Salpida and Pyrosomatida. The Co. inflata genome provides increased resolution of the ancestral Hox clusters of key tunicate nodes, therefore expanding our understanding of the evolution of this cluster and its potential impact on tunicate morphological diversity. Our analyses of other gene families revealed that several cardiovascular associated genes (e.g., BMP10, SCL2A12, and PDE2a) absent from Ci. robusta, are present in Co. inflata. Taken together, our results help clarify tunicate relationships and the genomic content of key ancestral nodes within this phylogeny, providing critical insights into tunicate evolution.

Description: The germlines of metazoans contain transposable elements (TEs) causing genetic instability and affecting fitness. To protect the germline from TE activity, gonads of metazoans produce TE-derived PIWI-interacting RNAs (piRNAs) that silence TE expression. In Drosophila, our understanding of piRNA biogenesis is mainly based on studies of the Drosophila melanogaster female germline. However, it is not known whether piRNA functions are also important in the male germline or whether and how piRNAs are affected by the global genomic context. To address these questions, we compared genome sequences, transcriptomes, and small RNA libraries extracted from entire testes and ovaries of two sister species: D. melanogaster and Drosophila simulans. We found that most TE-derived piRNAs were produced in ovaries and that piRNA pathway genes were strongly overexpressed in ovaries compared with testes, indicating that the silencing of TEs by the piRNA pathway mainly took place in the female germline. To study the relationship between host piRNAs and TE landscape, we analyzed TE genomic features and how they correlate with piRNA production in the two species. In D. melanogaster, we found that TE-derived piRNAs target recently active TEs. In contrast, although Drosophila simulans TEs do not display any features of recent activity, the host still intensively produced silencing piRNAs targeting old TE relics. Together, our results show that the piRNA silencing response mainly takes place in Drosophila ovaries and indicate that the host piRNA response is implemented following a burst of TE activity and could persist long after the extinction of active TE families.

Description: The evolution of the Tumor Necrosis Factor superfamily (TNFSF) in early vertebrates is inferred by comparing the TNFSF genes found in humans and nine fishes: three agnathans, two chondrichthyans, three actinopterygians and the sarcopterygian Latimeria chalumnae. By combining phylogenetic and synteny analyses, the TNFSF sequences detected are classified into five clusters of genes and 24 orthology groups. A model for their evolution since the origin of vertebrates is proposed. Fifteen TNFSF genes emerged from just three progenitors due to the whole genome duplications (WGDs) that occurred before the agnathan/gnathostome split. Later, gnathostomes not only kept most of the genes emerged in the WGDs but soon added several tandem duplicates. More recently, complex, lineage-specific patterns of duplications and losses occurred in different gnathostome lineages. In agnathan species only 7-8 TNFSF genes are detected, because this lineage soon lost six of the genes emerged in the ancestral WGDs and additional losses in both hagfishes and lampreys later occurred. The orthologs of many of these lost genes are, in mammals, ligands of death domain-containing TNFSF receptors, indicating that the extrinsic apoptotic pathway became simplified in the agnathan lineage. From the patterns of emergence of these genes, it is deduced that both the regulation of apoptosis and the control of the NF-κB pathway that depends in modern mammals on TNFSF members emerged before the ancestral vertebrate WGDs.

Description: The fungal-interactive (fungiphilic) strains BS001, BS007, BS110, and BS437 have previously been preliminarily assigned to the species Paraburkholderia terrae. However, in the (novel) genus Paraburkholderia, an as-yet unresolved subgroup exists, that clusters around Paraburkholderia hospita (containing the species P. terrae, P. hospita, and Paraburkholderia caribensis). To shed light on the precise relationships across the respective type strains and the novel fungiphiles, we here compare their genomic and ecophysiological features. To reach this goal, the genomes of the three type strains, with sizes ranging from 9.0 to 11.5 Mb, were de novo sequenced and the high-quality genomes analyzed. Using whole-genome, ribosomal RNA and marker-gene-concatenate analyses, close relationships between P. hospita DSM 17164T and P. terrae DSM 17804T, versus more remote relationships to P. caribensis DSM 13236T, were found. All four fungiphilic strains clustered closely to the two-species cluster. Analyses of average nucleotide identities (ANIm) and tetranucleotide frequencies (TETRA) confirmed the close relationships between P. hospita DSM 17164T and P. terrae DSM 17804T (ANIm = 95.42; TETRA = 0.99784), as compared with the similarities of each one of these strains to P. caribensis DSM 13236T. A species cluster was thus proposed. Furthermore, high similarities of the fungiphilic strains BS001, BS007, BS110, and BS437 with this cluster were found, indicating that these strains also make part of it, being closely linked to P. hospita DSM 17164T (ANIm = 99%; TETRA = 0.99). We propose to coin this cluster the P. hospita species cluster (containing P. hospita DSM 17164T, P. terrae DSM 17804T, and strains BS001, BS007, BS110, and BS437), being clearly divergent from the closely related species P. caribensis (type strain DSM 13236T). Moreover, given their close relatedness to P. hospita DSM 17164T within the cluster, we propose to rename the four fungiphilic strains as members of P. hospita. Analysis of migratory behavior along with fungal growth through soil revealed both P. terrae DSM 17804T and P. hospita DSM 17164T (next to the four fungiphilic strains) to be migration-proficient, whereas P. caribensis DSM 13236T was a relatively poor migrator. Examination of predicted functions across the genomes of the seven investigated strains, next to several selected additional ones, revealed the common presence of features in the P. hospita cluster strains that are potentially important in interactions with soil fungi. Thus, genes encoding specific metabolic functions, biofilm formation (pelABCDEFG, pgaABCD, alginate-related genes), motility/chemotaxis, type-4 pili, and diverse secretion systems were found.

Description: Gene duplication of green (RH2) opsin genes and their spectral differentiation are well documented in many teleost fish. However, their evolutionary divergence or conservation patterns among phylogenetically close but ecologically diverse species is not well explored. Medaka fish (genus Oryzias) are broadly distributed in fresh and brackish waters of Asia, with many species being laboratory-housed and feasible for genetic studies. We previously showed that a Japan strain (HNI) of medaka (Oryzias latipes) possessed three RH2 opsin genes (RH2-A, RH2-B, and RH2-C) encoding spectrally divergent photopigments. Here, we examined the three RH2 opsin genes from six Oryzias species representing three species groups: the latipes, the celebensis, and the javanicus. Photopigment reconstitution revealed that the peak absorption spectra (λmax) of RH2-A were divergent among the species (447–469 nm), whereas those of RH2-B and RH2-C were conservative (516–519 and 486–493 nm, respectively). For the RH2-A opsins, the largest spectral shift was detected in the phylogenetic branch leading to the latipes group. A single amino acid replacement T94C explained most of the spectral shift. For RH2-B and -C opsins, we detected tracts of gene conversion between the two genes homogenizing them. Nevertheless, several amino acid differences were maintained. We showed that the spectral difference between the two opsins was attributed to largely the E/Q amino acid difference at the site 122 and to several sites with individually small spectral effects. These results depict dynamism of spectral divergence of orthologous and paralogous green opsin genes in phylogenetically close but ecologically diverse species exemplified by medaka.

Description: Eukaryotes are typically depicted as descendants of archaea, but their genomes are evolutionary chimeras with genes stemming from archaea and bacteria. Which prokaryotic heritage predominates? Here, we have clustered 19,050,992 protein sequences from 5,443 bacteria and 212 archaea with 3,420,731 protein sequences from 150 eukaryotes spanning six eukaryotic supergroups. By downsampling, we obtain estimates for the bacterial and archaeal proportions. Eukaryotic genomes possess a bacterial majority of genes. On average, the majority of bacterial genes is 56% overall, 53% in eukaryotes that never possessed plastids, and 61% in photosynthetic eukaryotic lineages, where the cyanobacterial ancestor of plastids contributed additional genes to the eukaryotic lineage. Intracellular parasites, which undergo reductive evolution in adaptation to the nutrient rich environment of the cells that they infect, relinquish bacterial genes for metabolic processes. Such adaptive gene loss is most pronounced in the human parasite Encephalitozoon intestinalis with 86% archaeal and 14% bacterial derived genes. The most bacterial eukaryote genome sampled is rice, with 67% bacterial and 33% archaeal genes. The functional dichotomy, initially described for yeast, of archaeal genes being involved in genetic information processing and bacterial genes being involved in metabolic processes is conserved across all eukaryotic supergroups.

Description: Baboons are one of the most abundant large nonhuman primates and are widely studied in biomedical, behavioral, and anthropological research. Despite this, our knowledge of their evolutionary and demographic history remains incomplete. Here, we report a 0.9-fold coverage genome sequence from a 5800-year-old baboon from the site of Ha Makotoko in Lesotho. The ancient baboon is closely related to present-day Papio ursinus individuals from southern Africa—indicating a high degree of continuity in the southern African baboon population. This level of population continuity is rare in recent human populations but may provide a good model for the evolution of Homo and other large primates over similar timespans in structured populations throughout Africa.

Description: Our ability to correctly reconstruct a phylogenetic tree is strongly affected by both systematic errors and the amount of phylogenetic signal in the data. Current approaches to tackle tree reconstruction artifacts, such as the use of parameter-rich models, do not translate readily to single-gene alignments. This, coupled with the limited amount of phylogenetic information contained in single-gene alignments, makes gene trees particularly difficult to reconstruct. Opsin phylogeny illustrates this problem clearly. Opsins are G-protein coupled receptors utilized in photoreceptive processes across Metazoa and their protein sequences are roughly 300 amino acids long. A number of incongruent opsin phylogenies have been published and opsin evolution remains poorly understood. Here, we present a novel approach, the canary sequence approach, to investigate and potentially circumvent errors in single-gene phylogenies. First, we demonstrate our approach using two well-understood cases of long-branch attraction in single-gene data sets, and simulations. After that, we apply our approach to a large collection of well-characterized opsins to clarify the relationships of the three main opsin subfamilies.

Description: Introgressive hybridization results in the transfer of genetic material between species, often with fitness implications for the recipient species. The development of statistical methods for detecting the signatures of historical introgression in whole-genome data has been a major area of focus. Although existing techniques are able to identify the taxa that exchanged genes during introgression using a four-taxon system, most methods do not explicitly distinguish which taxon served as donor and which as recipient during introgression (i.e., polarization of introgression directionality). Existing methods that do polarize introgression are often only able to do so when there is a fifth taxon available and that taxon is sister to one of the taxa involved in introgression. Here, we present divergence-based introgression polarization (DIP), a method for polarizing introgression using patterns of sequence divergence across whole genomes, which operates in a four-taxon context. Thus, DIP can be applied to infer the directionality of introgression when additional taxa are not available. We use simulations to show that DIP can polarize introgression and identify potential sources of bias in the assignment of directionality, and we apply DIP to a well-described hominin introgression event.

Description: Animals depend on the sequential oxidation of organic molecules to survive; thus, oxygen-carrying/transporting proteins play a fundamental role in aerobic metabolism. Globins are the most common and widespread group of respiratory proteins. They can be divided into three types: circulating intracellular, noncirculating intracellular, and extracellular, all of which have been reported in annelids. The diversity of oxygen transport proteins has been underestimated across metazoans. We probed 250 annelid transcriptomes in search of globin diversity in order to elucidate the evolutionary history of this gene family within this phylum. We report two new globin types in annelids, namely androglobins and cytoglobins. Although cytoglobins and myoglobins from vertebrates and from invertebrates are referred to by the same name, our data show they are not genuine orthologs. Our phylogenetic analyses show that extracellular globins from annelids are more closely related to extracellular globins from other metazoans than to the intracellular globins of annelids. Broadly, our findings indicate that multiple gene duplication and neo-functionalization events shaped the evolutionary history of the globin family.

Description: Toll-like receptors (TLRs) are a complex family of innate immune genes that are well characterized in mammals and birds but less well understood in nonavian sauropsids (reptiles). The advent of highly contiguous draft genomes of nonmodel organisms enables study of such gene families through analysis of synteny and sequence identity. Here, we analyze TLR genes from the genomes of 22 tetrapod species. Findings reveal a TLR8 gene expansion in crocodilians and turtles (TLR8B), and a second duplication (TLR8C) specifically within turtles, followed by pseudogenization of that gene in the nonfreshwater species (desert tortoise and green sea turtle). Additionally, the Mojave desert tortoise (Gopherus agassizii) has a stop codon in TLR8B (TLR8-1) that is polymorphic among conspecifics. Revised orthology further reveals a new TLR homolog, TLR21-like, which is exclusive to lizards, snakes, turtles, and crocodilians. These analyses were made possible by a new draft genome assembly of the desert tortoise (gopAga2.0), which used chromatin-based assembly to yield draft chromosomal scaffolds (L50 = 26 scaffolds, N50 = 28.36 Mb, longest scaffold = 107 Mb) and an enhanced de novo genome annotation with 25,469 genes. Our three-step approach to orthology curation and comparative analysis of TLR genes shows what new insights are possible using genome assemblies with chromosome-scale scaffolds that permit integration of synteny conservation data.

Description: Here we present the draft genome sequence of Setaria digitata, a parasitic nematode affecting cattle. Due to its similarity to Wuchereria bancrofti, the parasitic nematode that causes lymphatic filariasis in humans, S. digitata has been used as a model organism at the genomic level to find drug targets which can be used for the development of novel drugs and/or vaccines for human filariasis. Setaria digitata causes cerebrospinal nematodiasis in goats, sheep, and horses posing a serious threat to livestock in developing countries. The genome sequence of S. digitata will assist in finding candidate genes to use as drug targets in both S. digitata and W. bancrofti. The assembled draft genome is ∼90 Mb long and contains 8,974 genomic scaffolds with a G+C content of 31.73%.

Description: Positive selection and purifying selection reduce levels of variation at linked neutral loci. One consequence of these processes is that the amount of neutral diversity and the meiotic recombination rate are predicted to be positively correlated across the genome—a prediction met in some species but not others. To better document the prevalence of selection at linked sites, we used new and published whole-genome sequences to survey nucleotide variation in population samples of the western European house mouse (Mus musculus domesticus) from Germany, France, and Gough Island, a remote volcanic island in the south Atlantic. Correlations between sequence variation and recombination rates estimated independently from dense linkage maps were consistently very weak (ρ ≤ 0.06), though they exceeded conventional significance thresholds. This pattern persisted in comparisons between genomic regions with the highest and lowest recombination rates, as well as in models incorporating the density of transcribed sites, the density of CpG dinucleotides, and divergence between mouse and rat as covariates. We conclude that natural selection affects linked neutral variation in a restricted manner in the western European house mouse.

Description: Members of genus Martes provide early warning signals about forest ecosystem health and are designated as a Management Indicator Species. As one of the most widespread members in Martes, the sable (Martes zibellina) is a circumboreal small predator found throughout all taiga zoogeographical zones of Eurasia and shows distinct population differentiation and morphological variations. To support further studies on striking local adaptation and population evolution, we present the first sable genome, assembled de novo from an individual originating in the Great Khingan Mountains (China). The assembled genome is 2.42 Gb, consisting of 15,814 scaffolds with a scaffold N50 of 5.20 Mb. Searches for complete Mammalia BUSCO (Benchmarking Universal Single-Copy Ortholog) gene groups found that 95.15% of the curated single-copy orthologs were assembled as complete, suggesting a high level of completeness of the genome. We totally predicted 19,413 protein-coding genes, and 0.82 Gb of repeat sequences was annotated. We also detected 1,257 olfactory receptor genes and found more functional olfactory receptor genes in sable than in other Mustelidae species, which provide a possible genetic explanation for the acute sense of smell of the sable for searching the preys under deep snow. Phylogenetic analyses revealed that the ferret (Mustela putorius furo) and sea otter (Enhydra lutris) form a clade that is sister to the sable, which was dated ∼16.4 Ma. Overall, our study provided the first reference genome for research in a broad range of areas including local adaptations, population evolution, conservation, and management for sable.

Description: The candidate phyla radiation (CPR) is a proposed subdivision within the bacterial domain comprising several candidate phyla. CPR organisms are united by small genome and physical sizes, lack several metabolic enzymes, and populate deep branches within the bacterial subtree of life. These features raise intriguing questions regarding their origin and mode of evolution. In this study, we performed a comparative and phylogenomic analysis to investigate CPR origin and evolution. Unlike previous gene/protein sequence-based reports of CPR evolution, we used protein domain superfamilies classified by protein structure databases to resolve the evolutionary relationships of CPR with non-CPR bacteria, Archaea, Eukarya, and viruses. Across all supergroups, CPR shared maximum superfamilies with non-CPR bacteria and were placed as deep branching bacteria in most phylogenomic trees. CPR contributed 1.22% of new superfamilies to bacteria including the ribosomal protein L19e and encoded four core superfamilies that are likely involved in cell-to-cell interaction and establishing episymbiotic lifestyles. Although CPR and non-CPR bacterial proteomes gained common superfamilies over the course of evolution, CPR and Archaea had more common losses. These losses mostly involved metabolic superfamilies. In fact, phylogenies built from only metabolic protein superfamilies separated CPR and non-CPR bacteria. These findings indicate that CPR are bacterial organisms that have probably evolved in an Archaea-like manner via the early loss of metabolic functions. We also discovered that phylogenies built from metabolic and informational superfamilies gave contrasting views of the groupings among Archaea, Bacteria, and Eukarya, which add to the current debate on the evolutionary relationships among superkingdoms.

Description: Much of the enormous phenotypic variation observed across human populations is thought to have arisen from events experienced as our ancestors peopled different regions of the world. However, little is known about the genes involved in these population-specific adaptations. Here, we explore this problem by simultaneously examining population-specific genetic and expression differentiation in four human populations. In particular, we derive a branch-based estimator of population-specific differentiation in four populations, and apply this statistic to single-nucleotide polymorphism and RNA-seq data from Italian, British, Finish, and Yoruban populations. As expected, genome-wide estimates of genetic and expression differentiation each independently recapitulate the known relationships among these four human populations, highlighting the utility of our statistic for identifying putative targets of population-specific adaptations. Moreover, genes with large copy number variations display elevated levels of population-specific genetic and expression differentiation, consistent with the hypothesis that gene duplication and deletion events are key reservoirs of adaptive variation. Further, many top-scoring genes are well-known targets of adaptation in Europeans, including those involved in lactase persistence and vitamin D absorption, and a handful of novel candidates represent promising avenues for future research. Together, these analyses reveal that our statistic can aid in uncovering genes involved in population-specific genetic and expression differentiation, and that such genes often play important roles in a diversity of adaptive and disease-related phenotypes in humans.

Description: Biodegradation of the phenylurea herbicide linuron appears a specialization within a specific clade of the Variovorax genus. The linuron catabolic ability is likely acquired by horizontal gene transfer but the mechanisms involved are not known. The full-genome sequences of six linuron-degrading Variovorax strains isolated from geographically distant locations were analyzed to acquire insight into the mechanisms of genetic adaptation toward linuron metabolism. Whole-genome sequence analysis confirmed the phylogenetic position of the linuron degraders in a separate clade within Variovorax and indicated that they unlikely originate from a common ancestral linuron degrader. The linuron degraders differentiated from Variovorax strains that do not degrade linuron by the presence of multiple plasmids of 20–839 kb, including plasmids of unknown plasmid groups. The linuron catabolic gene clusters showed 1) high conservation and synteny and 2) strain-dependent distribution among the different plasmids. Most of them were bordered by IS1071 elements forming composite transposon structures, often in a multimeric array configuration, appointing IS1071 as a key element in the recruitment of linuron catabolic genes in Variovorax. Most of the strains carried at least one (catabolic) broad host range plasmid that might have been a second instrument for catabolic gene acquisition. We conclude that clade 1 Variovorax strains, despite their different geographical origin, made use of a limited genetic repertoire regarding both catabolic functions and vehicles to acquire linuron biodegradation.

Description: Through its long history of artificial selection, the rock pigeon (Columba livia Gmelin 1789) was forged into a large number of domestic breeds. The incredible amount of phenotypic diversity exhibited in these breeds has long held the fascination of scholars, particularly those interested in biological inheritance and evolution. However, exploiting them as a model system is challenging, as unlike with many other domestic species, few reliable records exist about the origins of, and relationships between, each of the breeds. Therefore, in order to broaden our understanding of the complex evolutionary relationships among pigeon breeds, we generated genome-wide data by performing the genotyping-by-sequencing (GBS) method on close to 200 domestic individuals representing over 60 breeds. We analyzed these GBS data alongside previously published whole-genome sequencing data, and this combined analysis allowed us to conduct the most extensive phylogenetic analysis of the group, including two feral pigeons and one outgroup. We improve previous phylogenies, find considerable population structure across the different breeds, and identify unreported interbreed admixture events. Despite the reduced number of loci relative to whole-genome sequencing, we demonstrate that GBS data provide sufficient analytical power to investigate intertwined evolutionary relationships, such as those that are characteristic of animal domestic breeds. Thus, we argue that future studies should consider sequencing methods akin to the GBS approach as an optimal cost-effective approach for addressing complex phylogenies.

Description: Many regulatory networks appear to involve partially redundant enhancers. Traditionally, such enhancers have been hypothesized to originate mainly by sequence duplication. An alternative model postulates that they arise independently, through convergent evolution. This mechanism appears to be counterintuitive to natural selection: Redundant sequences are expected to either diverge and acquire new functions or accumulate mutations and become nonfunctional. Nevertheless, we show that at least 31% of the redundant enhancer pairs in the human genome (and 17% in the mouse genome) indeed originated in this manner. Specifically, for virtually all transposon-derived redundant enhancer pairs, both enhancer partners have evolved independently, from the exaptation of two different transposons. In addition to conferring robustness to the system, redundant enhancers could provide an evolutionary advantage by fine-tuning gene expression. Consistent with this hypothesis, we observed that the target genes of redundant enhancers exhibit higher expression levels and tissue specificity as compared with other genes. Finally, we found that although enhancer redundancy appears to be an intrinsic property of certain mammalian regulatory networks, the corresponding enhancers are largely species-specific. In other words, the redundancy in these networks is most likely a result of convergent evolution.

Description: Comparative genomics has begun to elucidate the genomic basis of social life in insects, but insight into the genomic basis of spider sociality has lagged behind. To begin, to characterize genomic signatures associated with the evolution of social life in spiders, we performed one of the first spider comparative genomics studies including five solitary species and two social species, representing two independent origins of sociality in the genus Stegodyphus. We found that the two social spider species had a large expansion of gene families associated with transport and metabolic processes and an elevated genome-wide rate of molecular evolution compared with the five solitary spider species. Genes that were rapidly evolving in the two social species relative to the five solitary species were enriched for transport, behavior, and immune functions, whereas genes that were rapidly evolving in the solitary species were enriched for energy metabolism processes. Most rapidly evolving genes in the social species Stegodyphus dumicola were broadly expressed across four tissues and enriched for transport functions, but 12 rapidly evolving genes showed brain-specific expression and were enriched for social behavioral processes. Altogether, our study identifies putative genomic signatures and potential candidate genes associated with spider sociality. These results indicate that future spider comparative genomic studies, including broader sampling and additional independent origins of sociality, can further clarify the genomic causes and consequences of social life.

Description: It has long been challenging to uncover the molecular mechanisms behind striking morphological innovations such as mammalian pregnancy. We studied the power of a robust comparative orthology pipeline based on gene synteny to address such problems. We inferred orthology relations between human genes and genes from each of 43 other vertebrate genomes, resulting in ∼18,000 orthologous pairs for each genome comparison. By identifying genes that first appear coincident with origin of the placental mammals, we hypothesized that we would define a subset of the genome enriched for genes that played a role in placental evolution. We thus pinpointed orthologs that appeared before and after the divergence of eutherian mammals from marsupials. Reinforcing previous work, we found instead that much of the genetic toolkit of mammalian pregnancy evolved through the repurposing of preexisting genes to new roles. These genes acquired regulatory controls for their novel roles from a group of regulatory genes, many of which did in fact originate at the appearance of the eutherians. Thus, orthologs appearing at the origin of the eutherians are enriched in functions such as transcriptional regulation by Krüppel-associated box-zinc-finger proteins, innate immune responses, keratinization, and the melanoma-associated antigen protein class. Because the cellular mechanisms of invasive placentae are similar to those of metastatic cancers, we then used our orthology inferences to explore the association between placenta invasion and cancer metastasis. Again echoing previous work, we find that genes that are phylogenetically older are more likely to be implicated in cancer development.

Description: Bacterial genes are sometimes found to be inactivated by mutation. This inactivation may be observable simply because selection for function is intermittent or too weak to eliminate inactive alleles quickly. Here, I investigate cases in Salmonella enterica where inactivation is instead positively selected. These are identified by a rate of introduction of premature stop codons to a gene that is higher than expected under selective neutrality, as assessed by comparison to the rate of synonymous changes. I identify 84 genes that meet this criterion at a 10% false discovery rate. Many of these genes are involved in virulence, motility and chemotaxis, biofilm formation, and resistance to antibiotics or other toxic substances. It is hypothesized that most of these genes are subject to an ongoing process in which inactivation is favored under rare conditions, but the inactivated allele is deleterious under most other conditions and is subsequently driven to extinction by purifying selection.

Description: Though large multilocus genomic data sets have led to overall improvements in phylogenetic inference, they have posed the new challenge of addressing conflicting signals across the genome. In particular, ancestral population structure, which has been uncovered in a number of diverse species, can skew gene tree frequencies, thereby hindering the performance of species tree estimators. Here we develop a novel maximum likelihood method, termed TASTI (Taxa with Ancestral structure Species Tree Inference), that can infer phylogenies under such scenarios, and find that it has increasing accuracy with increasing numbers of input gene trees, contrasting with the relatively poor performances of methods not tailored for ancestral structure. Moreover, we propose a supertree approach that allows TASTI to scale computationally with increasing numbers of input taxa. We use genetic simulations to assess TASTI’s performance in the three- and four-taxon settings and demonstrate the application of TASTI on a six-species Afrotropical mosquito data set. Finally, we have implemented TASTI in an open-source software package for ease of use by the scientific community.

Description: The great gerbil (Rhombomys opimus) is a social rodent living in permanent, complex burrow systems distributed throughout Central Asia, where it serves as the main host of several important vector-borne infectious pathogens including the well-known plague bacterium (Yersinia pestis). Here, we present a continuous annotated genome assembly of the great gerbil, covering over 96% of the estimated 2.47-Gb genome. Taking advantage of the recent genome assemblies of the sand rat (Psammomys obesus) and the Mongolian gerbil (Meriones unguiculatus), comparative immunogenomic analyses reveal shared gene losses within TLR gene families (i.e., TLR8, TLR10, and the entire TLR11-subfamily) for Gerbillinae, accompanied with signs of diversifying selection of TLR7 and TLR9. Most notably, we find a great gerbil-specific duplication of the MHCII DRB locus. In silico analyses suggest that the duplicated gene provides high peptide binding affinity for Yersiniae epitopes as well as Leishmania and Leptospira epitopes, putatively leading to increased capability to withstand infections by these pathogens. Our study demonstrates the power of whole-genome sequencing combined with comparative genomic analyses to gain deeper insight into the immunogenomic landscape of the great gerbil and its close relatives.

Description: Phages can fundamentally alter the physiology and metabolism of their hosts. Although these phages are ubiquitous in the bacterial world, they have seldom been described among endosymbiotic bacteria. One notable exception is the APSE phage that is found associated with the gammaproteobacterial Hamiltonella defensa, hosted by several insect species. This secondary facultative endosymbiont is not necessary for the survival of its hosts but can infect certain individuals or even whole populations. Its infection in aphids is often associated with protection against parasitoid wasps. This protective phenotype has actually been linked to the infection of the symbiont strain with an APSE, which carries a toxin cassette that varies among so-called “types.” In the present work, we seek to expand our understanding of the diversity of APSE phages as well as the relations of their Hamiltonella hosts. For this, we assembled and annotated the full genomes of 16 APSE phages infecting Hamiltonella symbionts across ten insect species. Molecular and phylogenetic analyses suggest that recombination has occurred repeatedly among lineages. Comparative genomics of the phage genomes revealed two variable regions that are useful for phage typing. Additionally, we find that mobile elements could play a role in the acquisition of new genes in the toxin cassette. Altogether, we provide an unprecedented view of APSE diversity and their genome evolution across aphids. This genomic investigation will provide a valuable resource for the design and interpretation of experiments aiming at understanding the protective phenotype these phages confer to their insect hosts.

Description: The Oryza officinalis complex is the largest species group in Oryza, with more than nine species from four continents, and is a tertiary gene pool that can be exploited in breeding programs for the improvement of cultivated rice. Most diploid and tetraploid members of this group have a C genome. Using a new reference C genome for the diploid species Oryza officinalis, and draft genomes for two other C genome diploid species O. eichingeri and O. rhizomatis, we examine the influence of transposable elements on genome structure and provide a detailed phylogeny and evolutionary history of the Oryza C genomes. The O. officinalis genome is 1.6 times larger than the A genome of cultivated O. sativa, mostly due to proliferation of Gypsy type long-terminal repeat (LTR) transposable elements, but overall syntenic relationships are maintained with other Oryza genomes (A, B and F). Draft genome assemblies of the two other C genome diploid species, O. eichingeri and O. rhizomatis, and short-read resequencing of a series of other C genome species and accessions reveal that after the divergence of the C genome progenitor, there was still a substantial degree of variation within the C genome species through proliferation and loss of both DNA and LTR transposable elements. We provide a detailed phylogeny and evolutionary history of the Oryza C genomes, and a genomic resource for the exploitation of the Oryza tertiary gene pool.

Description: Major phenotypic innovations in social amoeba evolution occurred at the transition between the Polysphondylia and group 4 Dictyostelia, which comprise the model organism Dictyostelium discoideum, such as the formation of a new structure, the basal disk. Basal disk differentiation and robust stalk formation require the morphogen DIF-1, synthesized by the polyketide synthase StlB, the des-methyl-DIF-1 methyltransferase DmtA, and the chlorinase ChlA, which are conserved throughout Dictyostelia. To understand how the basal disk and other innovations evolved in group 4, we sequenced and annotated the Polysphondylium violaceum (Pvio) genome, performed cell type-specific transcriptomics to identify cell-type marker genes, and developed transformation and gene knock-out procedures for Pvio. We used the novel methods to delete the Pvio stlB gene. The Pvio stlB− mutants formed misshapen curly sorogens with thick and irregular stalks. As fruiting body formation continued, the upper stalks became more regular, but structures contained 40% less spores. The stlB− sorogens overexpressed a stalk gene and underexpressed a (pre)spore gene. Normal fruiting body formation and sporulation were restored in Pvio stlB− by including DIF-1 in the supporting agar. These data indicate that, although conserved, stlB and its product(s) acquired both a novel role in the group 4 Dictyostelia and a role opposite to that in its sister group.

Description: Relatively little is known about the processes shaping population structure in cooperatively breeding insect species, despite the long-hypothesized importance of population structure in shaping patterns of cooperative breeding. Polistes paper wasps are primitively eusocial insects, with a cooperative breeding system in which females often found nests in cooperative associations. Prior mark-recapture studies of Polistes have documented extreme female philopatry, although genetic studies frequently fail to detect the strong population structure expected for highly philopatric species. Together these findings have led to lack of consensus on the degree of dispersal and population structure in these species. This study assessed population structure of female Polistes fuscatus wasps at three scales: within a single site, throughout Central New York, and across the Northeastern United States. Patterns of spatial genetic clustering and isolation by distance were observed in nuclear and mitochondrial genomes at the continental scale. Remarkably, population structure was evident even at fine spatial scales within a single collection site. However, P. fuscatus had low levels of genetic differentiation across long distances. These results suggest that P. fuscatus wasps may employ multiple dispersal strategies, including extreme natal philopatry as well as longer-distance dispersal. We observed greater genetic differentiation in mitochondrial genes than in the nuclear genome, indicative of increased dispersal distances in males. Our findings support the hypothesis that limited female dispersal contributes toward population structure in paper wasps.

Description: The pervasiveness of sex despite its well-known costs is a long-standing puzzle in evolutionary biology. Current explanations for the success of sex in nature largely rely on the adaptive significance of the new or rare genotypes that sex may generate. Less explored is the possibility that sex-underlying molecular mechanisms can enhance fitness and convey benefits to the individuals that bear the immediate costs of sex. Here, we show that the molecular environment associated with self-fertilization can increase stress resistance in the ciliate Paramecium tetraurelia. This advantage is independent of new genetic variation, coupled with a reduced nutritional input, and offers fresh insights into the mechanistic origin of sex. In addition to providing evidence that the molecular underpinnings of sexual reproduction and the stress response are linked in P. tetraurelia, these findings supply an integrative explanation for the persistence of self-fertilization in this ciliate.

Description: Koreans are thought to be an ethnic group of admixed northern and southern subgroups. However, the exact genetic origins of these two remain unclear. In addition, the past admixture is presumed to have taken place on the Korean peninsula, but there is no genomic scale analysis exploring the origin, composition, admixture, or the past migration of Koreans. Here, 88 Korean genomes compared with 91 other present-day populations showed two major genetic components of East Siberia and Southeast Asia. Additional paleogenomic analysis with 115 ancient genomes from Pleistocene hunter-gatherers to Iron Age farmers showed a gradual admixture of Tianyuan (40 ka) and Devil’s gate (8 ka) ancestries throughout East Asia and East Siberia up until the Neolithic era. Afterward, the current genetic foundation of Koreans may have been established through a rapid admixture with ancient Southern Chinese populations associated with Iron Age Cambodians. We speculate that this admixing trend initially occurred mostly outside the Korean peninsula followed by continuous spread and localization in Korea, corresponding to the general admixture trend of East Asia. Over 70% of extant Korean genetic diversity is explained to be derived from such a recent population expansion and admixture from the South.

Description: Retrotransposons, DNA sequences capable of creating copies of themselves, compose about half of the human genome and played a central role in the evolution of mammals. Their current position in the host genome is the result of the retrotranscription process and of the following host genome evolution. We apply a model from statistical physics to show that the genomic distribution of the two most populated classes of retrotransposons in human deviates from random placement, and that this deviation increases with time. The time dependence suggests a major role of the host genome dynamics in shaping the current retrotransposon distributions. Focusing on a neutral scenario, we show that a simple model based on random placement followed by genome expansion and sequence duplications can reproduce the empirical retrotransposon distributions, even though more complex and possibly selective mechanisms can have contributed. Besides the inherent interest in understanding the origin of current retrotransposon distributions, this work sets a general analytical framework to analyze quantitatively the effects of genome evolutionary dynamics on the distribution of genomic elements.

Description: Genome-wide nucleotide composition varies widely among species. Despite extensive research, the source of genome-wide nucleotide composition diversity remains elusive. Yeast mitochondrial genomes (mitogenomes) are highly A+T rich, and they provide a unique opportunity to study the evolution of AT-biased landscape. In this study, we sequenced ten complete mitogenomes of the Saccharomycodes ludwigii yeast with 8% G+C content, the lowest genome-wide %(G+C) in all published genomes to date. The Saccharomycodes ludwigii mitogenomes have high densities of short tandem repeats but severely underrepresented mononucleotide repeats. Comparative population genomics of these record-setting A+T rich genomes show dynamic indel mutations and strong mutation bias towards A/T. Indel mutations play a greater role in genomic variation among very closely related strains than nucleotide substitutions. Indels have resulted in presence-absence polymorphism of tRNAArg (ACG) among Saccharomycodes ludwigii mitogenomes. Interestingly, these mitogenomes have undergone recombination, a genetic process that can increase G+C content by GC-biased gene conversion. Finally, the expected equilibrium G+C content under mutation pressure alone is higher than observed G+C content, suggesting existence of mechanisms other than AT-biased mutation operating to increase A/T. Together, our findings shed new lights on mechanisms driving extremely AT-rich genomes.

Description: Ancestral sequence reconstruction (ASR) uses an alignment of extant protein sequences, a phylogeny describing the history of the protein family and a model of the molecular-evolutionary process to infer the sequences of ancient proteins, allowing researchers to directly investigate the impact of sequence evolution on protein structure and function. Like all statistical inferences, ASR can be sensitive to violations of its underlying assumptions. Previous studies have shown that, whereas phylogenetic uncertainty has only a very weak impact on ASR accuracy, uncertainty in the protein sequence alignment can more strongly affect inferred ancestral sequences. Here, we show that errors in sequence alignment can produce errors in ASR across a range of realistic and simplified evolutionary scenarios. Importantly, sequence reconstruction errors can lead to errors in estimates of structural and functional properties of ancestral proteins, potentially undermining the reliability of analyses relying on ASR. We introduce an alignment-integrated ASR approach that combines information from many different sequence alignments. We show that integrating alignment uncertainty improves ASR accuracy and the accuracy of downstream structural and functional inferences, often performing as well as highly accurate structure-guided alignment. Given the growing evidence that sequence alignment errors can impact the reliability of ASR studies, we recommend that future studies incorporate approaches to mitigate the impact of alignment uncertainty. Probabilistic modeling of insertion and deletion events has the potential to radically improve ASR accuracy when the model reflects the true underlying evolutionary history, but further studies are required to thoroughly evaluate the reliability of these approaches under realistic conditions.

Description: Gene losses in plastid genomes (plastomes) are often accompanied by functional transfer to the nucleus or substitution of an alternative nuclear-encoded gene. Despite the highly conserved gene content in plastomes of photosynthetic land plants, recent gene loss events have been documented in several disparate angiosperm clades. Among these lineages, Passiflora lacks several essential ribosomal genes, rps7, rps16, rpl20, rpl22, and rpl32, the two largest plastid genes, ycf1 and ycf2, and has a highly divergent rpoA. Comparative transcriptome analyses were performed to determine the fate of the missing genes in Passiflora. Putative functional transfers of rps7, rpl22, and rpl32 to nucleus were detected, with the nuclear transfer of rps7, representing a novel event in angiosperms. Plastid-encoded rps7 was transferred into the intron of a nuclear-encoded plastid-targeted thioredoxin m-type gene, acquiring its plastid transit peptide (TP). Plastid rpl20 likely experienced a novel substitution by a duplicated, nuclear-encoded mitochondrial-targeted rpl20 that has a similar gene structure. Additionally, among rosids, evidence for a third independent transfer of rpl22 in Passiflora was detected that gained a TP from a nuclear gene containing an organelle RNA recognition motif. Nuclear transcripts representing rpoA, ycf1, and ycf2 were not detected. Further analyses suggest that the divergent rpoA remains functional and that the gene is under positive or purifying selection in different clades. Comparative analyses indicate that alternative translocon and motor protein complexes may have substituted for the loss of ycf1 and ycf2 in Passiflora.