The genomic and cellular basis of biosynthetic innovation in rove beetles.

How evolution at the cellular level potentiates macroevolutionary change is central to understanding biological diversification. The >66,000 rove beetle species (Staphylinidae) form the largest metazoan family. Combining genomic and cell type transcriptomic insights spanning the largest clade, Aleocharinae, we retrace evolution of two cell types comprising a defensive gland-a putative catalyst behind staphylinid megadiversity. We identify molecular evolutionary steps leading to benzoquinone production by one cell type via a mechanism convergent with plant toxin release systems, and synthesis by the second cell type of a solvent that weaponizes the total secretion. This cooperative system has been conserved since the Early Cretaceous as Aleocharinae radiated into tens of thousands of lineages. Reprogramming each cell type yielded biochemical novelties enabling ecological specialization-most dramatically in symbionts that infiltrate social insect colonies via host-manipulating secretions. Our findings uncover cell type evolutionary processes underlying the origin and evolvability of a beetle chemical innovation.

Gene-rich X chromosomes implicate intragenomic conflict in the evolution of bizarre genetic systems.

Haplodiploidy and paternal genome elimination (HD/PGE) are common in invertebrates, having evolved at least two dozen times, all from male heterogamety (i.e., systems with X chromosomes). However, why X chromosomes are important for the evolution of HD/PGE remains debated. The Haploid Viability Hypothesis posits that X-linked genes promote the evolution of male haploidy by facilitating purging recessive deleterious mutations. The Intragenomic Conflict Hypothesis holds that conflict between genes drives genetic system turnover; under this model, X-linked genes could promote the evolution of male haploidy due to conflicts with autosomes over sex ratios and genetic transmission. We studied lineages where we can distinguish these hypotheses: species with germline PGE that retain an XX/X0 sex determination system (gPGE+X). Because evolving PGE in these cases involves changes in transmission without increases in male hemizygosity, a high degree of X linkage in these systems is predicted by the Intragenomic Conflict Hypothesis but not the Haploid Viability Hypothesis. To quantify the degree of X linkage, we sequenced and compared 7 gPGE+X species' genomes with 11 related species with typical XX/XY or XX/X0 genetic systems, representing three transitions to gPGE. We find highly increased X linkage in both modern and ancestral genomes of gPGE+X species compared to non-gPGE relatives and recover a significant positive correlation between percent X linkage and the evolution of gPGE. These empirical results substantiate longstanding proposals for a role for intragenomic conflict in the evolution of genetic systems such as HD/PGE.

Insights into the genome of the 'Loco' Concholepas concholepas (Gastropoda: Muricidae) from low-coverage short-read sequencing: genome size, ploidy, transposable elements, nuclear RNA gene operon, mitochondrial genome, and phylogenetic placement in the family Muricidae.

BACKGROUND: The Peruvian 'chanque' or Chilean 'loco' Concholepas concholepas is an economically, ecologically, and culturally important muricid gastropod heavily exploited by artisanal fisheries in the temperate southeastern Pacific Ocean. In this study, we have profited from a set of bioinformatics tools to recover important biological information of C. concholepas from low-coverage short-read NGS datasets. Specifically, we calculated the size of the nuclear genome, ploidy, and estimated transposable elements content using an in silico k-mer approach, we discovered, annotated, and quantified those transposable elements, we assembled and annotated the 45S rDNA RNA operon and mitochondrial genome, and we confirmed the phylogenetic position of C. concholepas within the muricid subfamily Rapaninae based on translated protein coding genes. RESULTS: Using a k-mer approach, the haploid genome size estimated for the predicted diploid genome of C. concholepas varied between 1.83 Gbp (with kmer = 24) and 2.32 Gbp (with kmer = 36). Between half and two thirds of the nuclear genome of C. concholepas was composed of transposable elements. The most common transposable elements were classified as Long Interspersed Nuclear Elements and Short Interspersed Nuclear Elements, which were more abundant than DNA transposons, simple repeats, and Long Terminal Repeats. Less abundant repeat elements included Helitron mobile elements, 45S rRNA DNA, and Satellite DNA, among a few others.The 45S rRNA DNA operon of C. concholepas that encodes for the ssrRNA, 5.8S rRNA, and lsrRNA genes was assembled into a single contig 8,090 bp long. The assembled mitochondrial genome of C. concholepas is 15,449 bp long and encodes 13 protein coding genes, two ribosomal genes, and 22 transfer RNAs. CONCLUSION: The information gained by this study will inform the assembly of a high quality nuclear genome for C. concholepas and will support bioprospecting and biomonitoring using environmental DNA to advance development of conservation and management plans in this overexploited marine snail.

Whole Snake Genomes from Eighteen Families of Snakes (Serpentes: Caenophidia) and Their Applications to Systematics.

We present genome assemblies for 18 snake species representing 18 families (Serpentes: Caenophidia): Acrochordus granulatus, Aparallactus werneri, Boaedon fuliginosus, Calamaria suluensis, Cerberus rynchops, Grayia smithii, Imantodes cenchoa, Mimophis mahfalensis, Oxyrhabdium leporinum, Pareas carinatus, Psammodynastes pulverulentus, Pseudoxenodon macrops, Pseudoxyrhopus heterurus, Sibynophis collaris, Stegonotus admiraltiensis, Toxicocalamus goodenoughensis, Trimeresurus albolabris, and Tropidonophis doriae. From these new genome assemblies, we extracted thousands of loci commonly used in systematic and phylogenomic studies on snakes, including target-capture datasets composed of UCEs and AHEs, as well as traditional Sanger loci. Phylogenies inferred from the two target-capture loci datasets were identical with each other, and strongly congruent with previously published snake phylogenies. To show additional utility of these non-model genomes for investigative evolutionary research, we mined the genome assemblies of two New Guinea island endemics in our dataset (Stegonotus admiraltiensis and Tropidonophis doriae) for the ATP1a3 gene, a thoroughly researched indicator of resistance to toad toxin ingestion by squamates. We find that both these snakes possess the genotype for toad toxin resistance despite their endemism to New Guinea, a region absent of any toads until the human-mediated introduction of Cane Toads in the 1930s. These species possess identical substitutions that suggest the same bufotoxin resistance as their Australian congenerics (Stegonotus cucullatus and Tropidonophis mairii) which forage on invasive Cane Toads. Herein, we show the utility of short-read high coverage genomes, as well as improving the deficit of available squamate genomes with associated voucher specimens.

'Candidatus Tisiphia' is a widespread Rickettsiaceae symbiont in the mosquito Anopheles plumbeus (Diptera: Culicidae).

Symbiotic bacteria can alter host biology by providing protection from natural enemies, or alter reproduction or vectoral competence. Symbiont-linked control of vector-borne disease in Anopheles has been hampered by a lack of symbioses that can establish stable vertical transmission in the host. Previous screening found the symbiont 'Candidatus Tisiphia' in Anopheles plumbeus, an aggressive biter and potential secondary vector of malaria parasites and West Nile virus. We screened samples collected over 10-years across Germany and used climate databases to assess environmental influence on incidence. We observed a 95% infection rate, and that the frequency of infection did not fluctuate with broad environmental factors. Maternal inheritance is indicated by presence in the ovaries through FISH microscopy. Finally, we assembled a high-quality 1.6 Mbp draft genome of 'Ca. Tisiphia' to explore its phylogeny and potential metabolic competence. The infection is closely related to strains found in Culicoides biting midges and shows similar patterns of metabolism, providing no evidence of the capacity to synthesize B-vitamins. This infection offers avenues for onward research in anopheline mosquito symbioses. Additionally, it provides future opportunity to study the impact of 'Ca. Tisiphia' on natural and transinfected hosts, especially in relation to reproductive fitness and vectorial competence and capacity.

Phylogenomics of Bony-Tongue Fishes (Osteoglossomorpha) Shed Light on the Craniofacial Evolution and Biogeography of the Weakly Electric Clade (Mormyridae).

Bonytongues (Osteoglossomorpha) constitute an ancient clade of teleost fishes distributed in freshwater habitats throughout the world. The group includes well-known species such as arowanas, featherbacks, pirarucus, and the weakly electric fishes in the family Mormyridae. Their disjunct distribution, extreme morphologies, and electrolocating capabilities (Gymnarchidae and Mormyridae) have attracted much scientific interest, but a comprehensive phylogenetic framework for comparative analysis is missing, especially for the species-rich family Mormyridae. Of particular interest are disparate craniofacial morphologies among mormyrids which might constitute an exceptional model system to study convergent evolution. We present a phylogenomic analysis based on 546 exons of 179 species (out of 260), 28 out of 29 genera, and all six families of extant bonytongues. Based on a recent reassessment of the fossil record of osteoglossomorphs, we inferred dates of divergence among transcontinental clades and the major groups. The estimated ages of divergence among extant taxa (e.g., Osteoglossomorpha, Osteoglossiformes, and Mormyroidea) are older than previous reports, but most of the divergence dates obtained for clades on separate continents are too young to be explained by simple vicariance hypotheses. Biogeographic analysis of mormyrids indicates that their high species diversity in the Congo Basin is a consequence of range reductions of previously widespread ancestors and that the highest diversity of craniofacial morphologies among mormyrids originated in this basin. Special emphasis on a taxon-rich representation for mormyrids revealed pervasive misalignment between our phylogenomic results and mormyrid taxonomy due to repeated instances of convergence for extreme craniofacial morphologies. Estimation of ancestral phenotypes revealed contingent evolution of snout elongation and unique projections from the lower jaw to form the distinctive Schnauzenorgan. Synthesis of comparative analyses suggests that the remarkable craniofacial morphologies of mormyrids evolved convergently due to niche partitioning, likely enabled by interactions between their exclusive morphological and electrosensory adaptations. [Africa; ancestral state estimation; diversity; exon capture; freshwater fishes; Phylogenomics.].

Genomic Assessment of the Contribution of the Wolbachia Endosymbiont of Eurosta solidaginis to Gall Induction.

We explored the genome of the Wolbachia strain, wEsol, symbiotic with the plant-gall-inducing fly Eurosta solidaginis with the goal of determining if wEsol contributes to gall induction by its insect host. Gall induction by insects has been hypothesized to involve the secretion of the phytohormones cytokinin and auxin and/or proteinaceous effectors to stimulate cell division and growth in the host plant. We sequenced the metagenome of E. solidaginis and wEsol and assembled and annotated the genome of wEsol. The wEsol genome has an assembled length of 1.66 Mbp and contains 1878 protein-coding genes. The wEsol genome is replete with proteins encoded by mobile genetic elements and shows evidence of seven different prophages. We also detected evidence of multiple small insertions of wEsol genes into the genome of the host insect. Our characterization of the genome of wEsol indicates that it is compromised in the synthesis of dimethylallyl pyrophosphate (DMAPP) and S-adenosyl L-methionine (SAM), which are precursors required for the synthesis of cytokinins and methylthiolated cytokinins. wEsol is also incapable of synthesizing tryptophan, and its genome contains no enzymes in any of the known pathways for the synthesis of indole-3-acetic acid (IAA) from tryptophan. wEsol must steal DMAPP and L-methionine from its host and therefore is unlikely to provide cytokinin and auxin to its insect host for use in gall induction. Furthermore, in spite of its large repertoire of predicted Type IV secreted effector proteins, these effectors are more likely to contribute to the acquisition of nutrients and the manipulation of the host's cellular environment to contribute to growth and reproduction of wEsol than to aid E. solidaginis in manipulating its host plant. Combined with earlier work that shows that wEsol is absent from the salivary glands of E. solidaginis, our results suggest that wEsol does not contribute to gall induction by its host.

The Complete Genome Sequence of Ajuga reptans (Lamiaceae, Lamiales), Bugleweed.

We present the whole genome sequence of Ajuga reptans. Illumina paired-end reads were assembled by a de novo method followed by a finishing step. The raw and assembled data are publicly available via GenBank: Sequence Read Archive (SRR24502601) and assembled genome (JAUEKW000000000).

The Complete Genome Sequences of 7 Hypoplectrus species (Serranidae, Perciformes).

Species in the genus Hypoplectrus are found primarily in coral reefs in the Caribbean Sea and the Gulf of Mexico, particularly around Florida and the Bahamas. They are a popular choice for hobbyist saltwater aquariums. Members of this genus are simultaneous hermaphrodites, as they have both male and female sexual organs at the same time as an adult. We present the genome sequences of 7 species of Hypoplectrus. Raw read data and assembled genomes are available in Genbank.

The Genome Sequences of 118 Taxonomically Diverse Eukaryotes of the Salish Sea.

We present the complete genome sequences of 118 taxonomically diverse eukaryotes from the Salish Sea. Illumina sequencing was performed on genetic material from wild-collected individuals. The reads were assembled using a de novo method followed by a finishing step. The raw and assembled data are publicly available via Genbank.

The complete genome sequences of two species of seventeen-year cicadas: Magicicada septendecim and Magicicada septendecula.

The genus Magicicada (Hemiptera: Cicadidae) includes the periodical cicadas of Eastern North America. Spending the majority of their long lives underground, the adult cicadas emerge every 13 or 17 years to spend 4-6 weeks as adult to mate. We present the whole genome sequences of two species of 17-year cicadas, Magicicada septendecim and Magicicada septendecula. The reads were assembled by a de novo method followed by alignments to related species. Annotation was performed by GeneMark-ES. The raw and assembled data is available via NCBI Short Read Archive and Assembly databases.

Nuclear and mitochondrial genome datasets for spiny lobsters genus Panulirus (Decapoda: Achelata: Palinuridae).

Spiny lobsters (Decapoda: Palinuridae) in the genus Panulirus are targets of lucrative fisheries globally and have relevant ecological functions in tropical and subtropical environments. Only a few, but increasing, number of genetic and genomic resources exist for them. Nuclear and mitochondrial genome assemblies can provide insights into their phylogenetic relationships and support fishery management strategies in species that are heavily exploited. Herein, using Illumina short reads whole genome sequencing, we assembled the nuclear and mitochondrial genomes of a total of 14 species. Genomic DNA was extracted from specimens deposited at Clemson University Crustacean Collection and sequenced in a HiSeq X Ten system. The number of paired-end (PE) reads generated for the different studied species varied between 219,917,346 in P. argus and 70,215,423 in P. cygnus. Nuclear and mitochondrial genomes were 'de novo' assembled. Nuclear genomes ranged between 1,624,400,357 bp in P. guttatus and 935,571,898 bp in P. cygnus with scaffold numbers varying between 466,583 in P. versicolor and 852,228 in P. longipes. Mitochondrial genomes varied between 15,613 bp and 15,768 bp in P. pascuensis and P. versicolor, respectively. The totality of the short reads, nuclear, and mitochondrial genome assemblies are available at NCBI's GenBank.

The Draft Genome Sequences of 63 Species of Catfish (Ariidae, Siluriformes).

We present genome sequences of 63 species of catfish from the Ariidae family. Illumina sequencing was performed on tissue from wild-collected museum specimens. The reads were assembled using a de novo method followed by a finishing step. The raw and assembled data are publicly available via Genbank.

The Complete Genome Sequences of 12 Species of Enteromius (Cyprinidae, Cypriniformes, Actinopteri).

We present the complete genome sequences of 12 species of Enteromius. Illumina sequencing was performed on genetic material from museum specimens. The reads were assembled using a de novo method followed by a finishing step. The raw and assembled data are publicly available via Genbank.

The Complete Genome Sequences of seven species of sea slugs (Nudibranchia, Gastropoda, Mollusca).

We present the complete genome sequences of seven species of sea slugs. Illumina sequencing was performed on tissue from wild-collected museum specimens. The reads were assembled using a de novo method followed by a finishing step. The raw and assembled data are publicly available via Genbank.

The Complete Genome Sequences of 23 Finch Species (Fringillidae, Passeriformes).

We present complete genome sequences of 23 species of finches from 6 genera.

The complete genome sequence of the protogynous hermaphrodite North American annual killifish Millerichthys robustus (Rivulidae, Cyprinodontiformes).

The Mexican Rivulus, Millerichthys robustus, is an enigmatic species of seasonal killifish endemic of the Southeast Mexico that has changed paradigms on the evolution of annualism in killifishes. This species survives in ephemeral environments that experience a period of seasonal drought that causes the death of all adult fish. However, populations persist due to their drought-resistant embryos capable of arresting their development in diapause until the next the rainy season. Millerichthys evolved unique characteristics within annual killifishes as a functional sequential hermaphroditism, in which females change to males (protogynous) under perceived conditions of mate competition. Also, M. robustus express different color phenotypes in both sexes: five-color phenotypes continuously distributed in various perceptual units between yellow and red in males, and different number of ocelli disposed in caudal peduncle in females. The phylogenetic relationships of M. robustus revealed that it is a sister clade to two non-annual species found exclusively in Cuba (Rivulus cylindraceus and Rivulus berovidesi), indicating that the annual life cycle, through the acquisition of embryonic diapause, has evolved independently in this species. Here, we present the complete genome sequences for the North American annual killifish Millerichthys robustus. The raw data and assembled genome are available in GeneBank.

The Complete Genome Sequences of 37 Bird Species from the Ecuadorian Andes.

We present the complete genome sequences of 37 Ecuadorian bird species widespread throughout the tropical Andes.

Phylogenomics of Psammodynastes and Buhoma (Elapoidea: Serpentes), with the description of a new Asian snake family.

Asian mock vipers of the genus Psammodynastes and African forest snakes of the genus Buhoma are two genera belonging to the snake superfamily Elapoidea. The phylogenetic placements of Psammodynastes and Buhoma within Elapoidea has been extremely unstable which has resulted in their uncertain and debated taxonomy. We used ultraconserved elements and traditional nuclear and mitochondrial markers to infer the phylogenetic relationships of these two genera with other elapoids. Psammodynastes, for which a reference genome has been sequenced, were found, with strong branch support, to be a relatively early diverging split within Elapoidea that is sister to a clade consisting of Elapidae, Micrelapidae and Lamprophiidae. Hence, we allocate Psammodynastes to its own family, Psammodynastidae new family. However, the phylogenetic position of Buhoma could not be resolved with a high degree of confidence. Attempts to identify the possible sources of conflict in the rapid radiation of elapoid snakes suggest that both hybridisation/introgression during the rapid diversification, including possible ghost introgression, as well as incomplete lineage sorting likely have had a confounding role. The usual practice of combining mitochondrial loci with nuclear genomic data appears to mislead phylogeny reconstructions in rapid radiation scenarios, especially in the absence of genome scale data.

The Complete Genome Sequence of Ceratonia siliqua L. (Fabaceae, Fabales), Carob.

We present the whole genome sequence of Ceratonia siliqua L. Illumina paired-end reads were assembled by a de novo method followed by a finishing step. The raw and assembled data are publicly available via GenBank: Sequence Read Archive (SRR24502586) and assembled genome (JASKGM000000000).