Chromosomal inversions from an initial ecotypic divergence drive a gradual repeated radiation of Galápagos beetles.

Island faunas exhibit some of the most iconic examples where similar forms repeatedly evolve within different islands. Yet, whether these deterministic evolutionary trajectories within islands are driven by an initial, singular divergence and the subsequent exchange of individuals and adaptive genetic variation between islands remains unclear. Here, we study a gradual, repeated evolution of low-dispersive highland ecotypes from a dispersive lowland ecotype of Calosoma beetles along the island progression of the Galápagos. We show that repeated highland adaptation involved selection on multiple shared alleles within extensive chromosomal inversions that originated from an initial adaptation event on the oldest island. These highland inversions first spread through dispersal of highland individuals. Subsequent admixture with the lowland ecotype resulted in polymorphic dispersive populations from which the highland populations evolved on the youngest islands. Our findings emphasize the significance of an ancient divergence in driving repeated evolution and highlight how a mixed contribution of inter-island colonization and within-island evolution can shape parallel species communities.

Identification of a non-canonical ciliate nuclear genetic code where UAA and UAG code for different amino acids.

The genetic code is one of the most highly conserved features across life. Only a few lineages have deviated from the "universal" genetic code. Amongst the few variants of the genetic code reported to date, the codons UAA and UAG virtually always have the same translation, suggesting that their evolution is coupled. Here, we report the genome and transcriptome sequencing of a novel uncultured ciliate, belonging to the Oligohymenophorea class, where the translation of the UAA and UAG stop codons have changed to specify different amino acids. Genomic and transcriptomic analyses revealed that UAA has been reassigned to encode lysine, while UAG has been reassigned to encode glutamic acid. We identified multiple suppressor tRNA genes with anticodons complementary to the reassigned codons. We show that the retained UGA stop codon is enriched in the 3'UTR immediately downstream of the coding region of genes, suggesting that there is functional drive to maintain tandem stop codons. Using a phylogenomics approach, we reconstructed the ciliate phylogeny and mapped genetic code changes, highlighting the remarkable number of independent genetic code changes within the Ciliophora group of protists. According to our knowledge, this is the first report of a genetic code variant where UAA and UAG encode different amino acids.

Different waves of postglacial recolonisation and genomic structure of bank vole populations in NE Poland.

Previous studies indicated that in some species phylogeographic patterns obtained in the analysis of nuclear and mitochondrial DNA (mtDNA) markers can be different. Such mitonuclear discordance can have important evolutionary and ecological consequences. In the present study, we aimed to check whether there was any discordance between mtDNA and nuclear DNA in the bank vole population in the contact zone of its two mtDNA lineages. We analysed the population genetic structure of bank voles using genome-wide genetic data (SNPs) and diversity of sequenced heart transcriptomes obtained from selected individuals from three populations inhabiting areas outside the contact zone. The SNP genetic structure of the populations confirmed the presence of at least two genetic clusters, and such division was concordant with the patterns obtained in the analysis of other genetic markers and functional genes. However, genome-wide SNP analyses revealed the more detailed structure of the studied population, consistent with more than two bank vole recolonisation waves, as recognised previously in the study area. We did not find any significant differences between individuals representing two separate mtDNA lineages of the species in functional genes coding for protein-forming complexes, which are involved in the process of cell respiration in mitochondria. We concluded that the contemporary genetic structure of the populations and the width of the contact zone were shaped by climatic and environmental factors rather than by genetic barriers. The studied populations were likely isolated in separate Last Glacial Maximum refugia for insufficient amount of time to develop significant genetic differentiation.

Fumigant Toxicity of Essential Oils against Frankliniella occidentalis and F. insularis (Thysanoptera: Thripidae) as Affected by Polymer Release and Adjuvants.

Frankliniella occidentalis is among the most economically significant pests of greenhouse crops, whose resistance to conventional insecticides has created demand for biopesticides such as essential oils. We assessed the fumigant toxicity of linalool against F. occidentalis, F. insularis, and Solanum lycopersicum. Thrips were fumigated with polyacrylamide hydrogels containing either (R)-linalool, (S)-linalool, racemic linalool, or a binary mixture of (R)-linalool with one of twelve adjuvants (i.e., peppermint, cedarwood, neem, clove, coconut, jojoba, soybean, olive, α-terpineol, 1,8-cineole, trans-anethole, or (R)-pulegone). Solanum lycopersicum seedlings were exposed to (R)-linalool or a mixture of (R)-linalool and peppermint oil via conditioned hydrogels or foliar spray. For F. insularis, (R)-linalool was more toxic than (S)-linalool, with LC50 values of 11.7 mg/L air and 16.7 mg/L air, respectively. Similarly for F. occidentalis, (R)-linalool was more toxic than (S)-linalool, with LC50 values of 29.0 mg/L air and 34.9 mg/L air, respectively. Peppermint oil and α-terpineol were the only synergists, while the other adjuvants exhibited varying degrees of antagonism. All seedling treatments demonstrated phytotoxicity, but symptoms were most severe for foliar sprays and mixtures containing peppermint oil. While hydrogels conditioned in linalool may be a favorable substitute to conventional insecticides, the cross-resistance demonstrated herein indicates that expectations should be metered.

Long read sequencing reveals novel isoforms and insights into splicing regulation during cell state changes.

BACKGROUND: Alternative splicing is a key mechanism underlying cellular differentiation and a driver of complexity in mammalian neuronal tissues. However, understanding of which isoforms are differentially used or expressed and how this affects cellular differentiation remains unclear. Long read sequencing allows full-length transcript recovery and quantification, enabling transcript-level analysis of alternative splicing processes and how these change with cell state. Here, we utilise Oxford Nanopore Technologies sequencing to produce a custom annotation of a well-studied human neuroblastoma cell line SH-SY5Y, and to characterise isoform expression and usage across differentiation. RESULTS: We identify many previously unannotated features, including a novel transcript of the voltage-gated calcium channel subunit gene, CACNA2D2. We show differential expression and usage of transcripts during differentiation identifying candidates for future research into state change regulation. CONCLUSIONS: Our work highlights the potential of long read sequencing to uncover previously unknown transcript diversity and mechanisms influencing alternative splicing.

Extensive microbial diversity within the chicken gut microbiome revealed by metagenomics and culture.

BACKGROUND: The chicken is the most abundant food animal in the world. However, despite its importance, the chicken gut microbiome remains largely undefined. Here, we exploit culture-independent and culture-dependent approaches to reveal extensive taxonomic diversity within this complex microbial community. RESULTS: We performed metagenomic sequencing of fifty chicken faecal samples from two breeds and analysed these, alongside all (n = 582) relevant publicly available chicken metagenomes, to cluster over 20 million non-redundant genes and to construct over 5,500 metagenome-assembled bacterial genomes. In addition, we recovered nearly 600 bacteriophage genomes. This represents the most comprehensive view of taxonomic diversity within the chicken gut microbiome to date, encompassing hundreds of novel candidate bacterial genera and species. To provide a stable, clear and memorable nomenclature for novel species, we devised a scalable combinatorial system for the creation of hundreds of well-formed Latin binomials. We cultured and genome-sequenced bacterial isolates from chicken faeces, documenting over forty novel species, together with three species from the genus Escherichia, including the newly named species Escherichia whittamii. CONCLUSIONS: Our metagenomic and culture-based analyses provide new insights into the bacterial, archaeal and bacteriophage components of the chicken gut microbiome. The resulting datasets expand the known diversity of the chicken gut microbiome and provide a key resource for future high-resolution taxonomic and functional studies on the chicken gut microbiome.

MolluscDB: a genome and transcriptome database for molluscs.

As sequencing becomes more accessible and affordable, the analysis of genomic and transcriptomic data has become a cornerstone of many research initiatives. Communities with a focus on particular taxa or ecosystems need solutions capable of aggregating genomic resources and serving them in a standardized and analysis-friendly manner. Taxon-focussed resources can be more flexible in addressing the needs of a research community than can universal or general databases. Here, we present MolluscDB, a genome and transcriptome database for molluscs. MolluscDB offers a rich ecosystem of tools, including an Ensembl browser, a BLAST server for homology searches and an HTTP server from which any dataset present in the database can be downloaded. To demonstrate the utility of the database and verify the quality of its data, we imported data from assembled genomes and transcriptomes of 22 species, estimated the phylogeny of Mollusca using single-copy orthologues, explored patterns of gene family size change and interrogated the data for biomineralization-associated enzymes and shell matrix proteins. MolluscDB provides an easy-to-use and openly accessible data resource for the research community. This article is part of the Theo Murphy meeting issue 'Molluscan genomics: broad insights and future directions for a neglected phylum'.

Deciphering mollusc shell production: the roles of genetic mechanisms through to ecology, aquaculture and biomimetics.

Most molluscs possess shells, constructed from a vast array of microstructures and architectures. The fully formed shell is composed of calcite or aragonite. These CaCO3 crystals form complex biocomposites with proteins, which although typically less than 5% of total shell mass, play significant roles in determining shell microstructure. Despite much research effort, large knowledge gaps remain in how molluscs construct and maintain their shells, and how they produce such a great diversity of forms. Here we synthesize results on how shell shape, microstructure, composition and organic content vary among, and within, species in response to numerous biotic and abiotic factors. At the local level, temperature, food supply and predation cues significantly affect shell morphology, whilst salinity has a much stronger influence across latitudes. Moreover, we emphasize how advances in genomic technologies [e.g. restriction site-associated DNA sequencing (RAD-Seq) and epigenetics] allow detailed examinations of whether morphological changes result from phenotypic plasticity or genetic adaptation, or a combination of these. RAD-Seq has already identified single nucleotide polymorphisms associated with temperature and aquaculture practices, whilst epigenetic processes have been shown significantly to modify shell construction to local conditions in, for example, Antarctica and New Zealand. We also synthesize results on the costs of shell construction and explore how these affect energetic trade-offs in animal metabolism. The cellular costs are still debated, with CaCO3 precipitation estimates ranging from 1-2 J/mg to 17-55 J/mg depending on experimental and environmental conditions. However, organic components are more expensive (~29 J/mg) and recent data indicate transmembrane calcium ion transporters can involve considerable costs. This review emphasizes the role that molecular analyses have played in demonstrating multiple evolutionary origins of biomineralization genes. Although these are characterized by lineage-specific proteins and unique combinations of co-opted genes, a small set of protein domains have been identified as a conserved biomineralization tool box. We further highlight the use of sequence data sets in providing candidate genes for in situ localization and protein function studies. The former has elucidated gene expression modularity in mantle tissue, improving understanding of the diversity of shell morphology synthesis. RNA interference (RNAi) and clustered regularly interspersed short palindromic repeats - CRISPR-associated protein 9 (CRISPR-Cas9) experiments have provided proof of concept for use in the functional investigation of mollusc gene sequences, showing for example that Pif (aragonite-binding) protein plays a significant role in structured nacre crystal growth and that the Lsdia1 gene sets shell chirality in Lymnaea stagnalis. Much research has focused on the impacts of ocean acidification on molluscs. Initial studies were predominantly pessimistic for future molluscan biodiversity. However, more sophisticated experiments incorporating selective breeding and multiple generations are identifying subtle effects and that variability within mollusc genomes has potential for adaption to future conditions. Furthermore, we highlight recent historical studies based on museum collections that demonstrate a greater resilience of molluscs to climate change compared with experimental data. The future of mollusc research lies not solely with ecological investigations into biodiversity, and this review synthesizes knowledge across disciplines to understand biomineralization. It spans research ranging from evolution and development, through predictions of biodiversity prospects and future-proofing of aquaculture to identifying new biomimetic opportunities and societal benefits from recycling shell products.

A newly developed genetic sex marker and its application to understanding chemically induced feminisation in roach (Rutilus rutilus).

Oestrogenic wastewater treatment works (WwTW) effluents discharged into UK rivers have been shown to affect sexual development, including inducing intersex, in wild roach (Rutilus rutilus). This can result in a reduced breeding capability with potential population level impacts. In the absence of a sex probe for roach it has not been possible to confirm whether intersex fish in the wild arise from genetic males or females, or whether sex reversal occurs in the wild, as this condition can be induced experimentally in controlled exposures to WwTW effluents and a steroidal oestrogen. Using restriction site-associated DNA sequencing (RAD-seq), we identified a candidate for a genetic sex marker and validated this marker as a sex probe through PCR analyses of samples from wild roach populations from nonpolluted rivers. We also applied the sex marker to samples from roach exposed experimentally to oestrogen and oestrogenic effluents to confirm suspected phenotypic sex reversal from males to females in some treatments, and also that sex-reversed males are able to breed as females. We then show, unequivocally, that intersex in wild roach populations results from feminisation of males, but find no strong evidence for complete sex reversal in wild roach at river sites contaminated with oestrogens. The discovered marker has utility for studies in roach on chemical effects, wild stock assessments, and reducing the number of fish used where only one sex is required for experimentation. Furthermore, we show that the marker can be applied nondestructively using a fin clip or skin swab, with animal welfare benefits.

Correction to: Assessing genotype-phenotype associations in three dorsal colour morphs in the meadow spittlebug Philaenus spumarius (L.) (Hemiptera: Aphrophoridae) using genomic and transcriptomic resources.

Following publication of the original article [1], it has been brought to the authors' attention that in their paper (Rodrigues et al. 2016) they reported the genome size based on 2C values (diploid genome) when it is more common to present it as 1C value.

Field cricket genome reveals the footprint of recent, abrupt adaptation in the wild.

Evolutionary adaptation is generally thought to occur through incremental mutational steps, but large mutational leaps can occur during its early stages. These are challenging to study in nature due to the difficulty of observing new genetic variants as they arise and spread, but characterizing their genomic dynamics is important for understanding factors favoring rapid adaptation. Here, we report genomic consequences of recent, adaptive song loss in a Hawaiian population of field crickets (Teleogryllus oceanicus). A discrete genetic variant, flatwing, appeared and spread approximately 15 years ago. Flatwing erases sound-producing veins on male wings. These silent flatwing males are protected from a lethal, eavesdropping parasitoid fly. We sequenced, assembled and annotated the cricket genome, produced a linkage map, and identified a flatwing quantitative trait locus covering a large region of the X chromosome. Gene expression profiling showed that flatwing is associated with extensive genome-wide effects on embryonic gene expression. We found that flatwing male crickets express feminized chemical pheromones. This male feminizing effect, on a different sexual signaling modality, is genetically associated with the flatwing genotype. Our findings suggest that the early stages of evolutionary adaptation to extreme pressures can be accompanied by greater genomic and phenotypic disruption than previously appreciated, and highlight how abrupt adaptation might involve suites of traits that arise through pleiotropy or genomic hitchhiking.

The Laboratory Domestication of Zebrafish: From Diverse Populations to Inbred Substrains.

We know from human genetic studies that practically all aspects of biology are strongly influenced by the genetic background, as reflected in the advent of "personalized medicine." Yet, with few exceptions, this is not taken into account when using laboratory populations as animal model systems for research in these fields. Laboratory strains of zebrafish (Danio rerio) are widely used for research in vertebrate developmental biology, behavior, and physiology, for modeling diseases, and for testing pharmaceutic compounds in vivo. However, all of these strains are derived from artificial bottleneck events and therefore are likely to represent only a fraction of the genetic diversity present within the species. Here, we use restriction site-associated DNA sequencing to genetically characterize wild populations of zebrafish from India, Nepal, and Bangladesh, and to compare them to previously published data on four common laboratory strains. We measured nucleotide diversity, heterozygosity, and allele frequency spectra, and find that wild zebrafish are much more diverse than laboratory strains. Further, in wild zebrafish, there is a clear signal of GC-biased gene conversion that is missing in laboratory strains. We also find that zebrafish populations in Nepal and Bangladesh are most distinct from all other strains studied, making them an attractive subject for future studies of zebrafish population genetics and molecular ecology. Finally, isolates of the same strains kept in different laboratories show a pattern of ongoing differentiation into genetically distinct substrains. Together, our findings broaden the basis for future genetic, physiological, pharmaceutic, and evolutionary studies in Danio rerio.

Detailed insights into pan-European population structure and inbreeding in wild and hatchery Pacific oysters (Crassostrea gigas) revealed by genome-wide SNP data.

Cultivated bivalves are important not only because of their economic value, but also due to their impacts on natural ecosystems. The Pacific oyster (Crassostrea gigas) is the world's most heavily cultivated shellfish species and has been introduced to all continents except Antarctica for aquaculture. We therefore used a medium-density single nucleotide polymorphism (SNP) array to investigate the genetic structure of this species in Europe, where it was introduced during the 1960s and has since become a prolific invader of coastal ecosystems across the continent. We analyzed 21,499 polymorphic SNPs in 232 individuals from 23 localities spanning a latitudinal cline from Portugal to Norway and including the source populations of Japan and Canada. We confirmed the results of previous studies by finding clear support for a southern and a northern group, with the former being indistinguishable from the source populations indicating the absence of a pronounced founder effect. We furthermore conducted a large-scale comparison of oysters sampled from the wild and from hatcheries to reveal substantial genetic differences including significantly higher levels of inbreeding in some but not all of the sampled hatchery cohorts. These findings were confirmed by a smaller but representative SNP dataset generated using restriction site-associated DNA sequencing. We therefore conclude that genomic approaches can generate increasingly detailed insights into the genetics of wild and hatchery produced Pacific oysters.

Preliminary insights into the genetics of bank vole tolerance to Puumala hantavirus in Sweden.

Natural reservoirs of zoonotic pathogens generally seem to be capable of tolerating infections. Tolerance and its underlying mechanisms remain difficult to assess using experiments or wildlife surveys. High-throughput sequencing technologies give the opportunity to investigate the genetic bases of tolerance, and the variability of its mechanisms in natural populations. In particular, population genomics may provide preliminary insights into the genes shaping tolerance and potentially influencing epidemiological dynamics. Here, we addressed these questions in the bank vole Myodes glareolus, the specific asymptomatic reservoir host of Puumala hantavirus (PUUV), which causes nephropathia epidemica (NE) in humans. Despite the continuous spatial distribution of M. glareolus in Sweden, NE is endemic to the northern part of the country. Northern bank vole populations in Sweden might exhibit tolerance strategies as a result of coadaptation with PUUV. This may favor the circulation and maintenance of PUUV and lead to high spatial risk of NE in northern Sweden. We performed a genome-scan study to detect signatures of selection potentially correlated with spatial variations in tolerance to PUUV. We analyzed six bank vole populations from Sweden, sampled from northern NE-endemic to southern NE-free areas. We combined candidate gene analyses (Tlr4, Tlr7, and Mx2 genes) and high-throughput sequencing of restriction site-associated DNA (RAD) markers. Outlier loci showed high levels of genetic differentiation and significant associations with environmental data including variations in the regional number of NE human cases. Among the 108 outliers that matched to mouse protein-coding genes, 14 corresponded to immune-related genes. The main biological pathways found to be significantly enriched corresponded to immune processes and responses to hantavirus, including the regulation of cytokine productions, TLR cascades, and IL-7, VEGF, and JAK-STAT signaling. In the future, genome-scan replicates and functional experimentations should enable to assess the role of these biological pathways in M. glareolus tolerance to PUUV.

The Genomic Basis of Color Pattern Polymorphism in the Harlequin Ladybird.

Many animal species comprise discrete phenotypic forms. A common example in natural populations of insects is the occurrence of different color patterns, which has motivated a rich body of ecological and genetic research [1-6]. The occurrence of dark, i.e., melanic, forms displaying discrete color patterns is found across multiple taxa, but the underlying genomic basis remains poorly characterized. In numerous ladybird species (Coccinellidae), the spatial arrangement of black and red patches on adult elytra varies wildly within species, forming strikingly different complex color patterns [7, 8]. In the harlequin ladybird, Harmonia axyridis, more than 200 distinct color forms have been described, which classic genetic studies suggest result from allelic variation at a single, unknown, locus [9, 10]. Here, we combined whole-genome sequencing, population-based genome-wide association studies, gene expression, and functional analyses to establish that the transcription factor Pannier controls melanic pattern polymorphism in H. axyridis. We show that pannier is necessary for the formation of melanic elements on the elytra. Allelic variation in pannier leads to protein expression in distinct domains on the elytra and thus determines the distinct color patterns in H. axyridis. Recombination between pannier alleles may be reduced by a highly divergent sequence of ∼170 kb in the cis-regulatory regions of pannier, with a 50 kb inversion between color forms. This most likely helps maintain the distinct alleles found in natural populations. Thus, we propose that highly variable discrete color forms can arise in natural populations through cis-regulatory allelic variation of a single gene.

Inter and Intraspecific Genomic Divergence in Drosophila montana Shows Evidence for Cold Adaptation.

The genomes of species that are ecological specialists will likely contain signatures of genomic adaptation to their niche. However, distinguishing genes related to ecological specialism from other sources of selection and more random changes is a challenge. Here, we describe the genome of Drosophila montana, which is the most extremely cold-adapted Drosophila species known. We use branch tests to identify genes showing accelerated divergence in contrasts between cold- and warm-adapted species and identify about 250 genes that show differences, possibly driven by a lower synonymous substitution rate in cold-adapted species. We also look for evidence of accelerated divergence between D. montana and D. virilis, a previously sequenced relative, but do not find strong evidence for divergent selection on coding sequence variation. Divergent genes are involved in a variety of functions, including cuticular and olfactory processes. Finally, we also resequenced three populations of D. montana from across its ecological and geographic range. Outlier loci were more likely to be found on the X chromosome and there was a greater than expected overlap between population outliers and those genes implicated in cold adaptation between Drosophila species, implying some continuity of selective process at these different evolutionary scales.

A simple and robust real-time qPCR method for the detection of PIK3CA mutations.

PIK3CA mutations are seemingly the most common driver mutations in breast cancer with H1047R and E545K being the most common of these, accounting together for around 60% of all PIK3CA mutations and have promising therapeutic implications. Given the low sensitivity and the high cost of current genotyping methods we sought to develop fast, simple and inexpensive assays for PIK3CA H1047R and E545K mutation screening in clinical material. The methods we describe are based on a real-time PCR including a mutation specific primer combined with a non-productive oligonucleotide which inhibits wild-type amplification and a parallel internal control reaction. We demonstrate consistent detection of PIK3CA H1047R mutant DNA in genomic DNA extracted from frozen breast cancer biopsies, FFPE material or cancer cell lines with a detection sensitivity of approximately 5% mutant allele fraction and validate these results using both Sanger sequencing and deep next generation sequencing methods. The detection sensitivity for PIK3CA E545K mutation was approximately 10%. We propose these methods as simple, fast and inexpensive diagnostic tools to determine PIK3CA mutation status.

Maintaining their genetic distance: Little evidence for introgression between widely hybridizing species of Geum with contrasting mating systems.

Within the plant kingdom, many genera contain sister lineages with contrasting outcrossing and inbreeding mating systems that are known to hybridize. The evolutionary fate of these sister lineages is likely to be influenced by the extent to which they exchange genes. We measured gene flow between outcrossing Geum rivale and selfing Geum urbanum, sister species that hybridize in contemporary populations. We generated and used a draft genome of G. urbanum to develop dd-RAD data scorable in both species. Coalescent analysis of RAD data from allopatric populations indicated that the species diverged 2-3 Mya, and that historical gene flow between them was extremely low (1 migrant every 25 generations). Comparison of genetic divergence between species in sympatry and allopatry, together with an analysis of allele frequencies in potential parental and hybrid populations, provided no evidence of contemporary introgression in sympatric populations. Cluster- and species-specific marker analyses revealed that, apart from four early-generation hybrids, individuals in sympatric populations fell into two genetically distinct groups that corresponded exactly to their morphological species classification with maximum individual admixture estimates of only 1-3%. However, we did observe joint segregation of four putatively introgressed SNPs across two scaffolds in the G. urbanum population that was associated with significant morphological variation, interpreted as tentative evidence for rare, recent interspecific gene flow. Overall, our results indicate that despite the presence of hybrids in contemporary populations, genetic exchange between G. rivale and G. urbanum has been extremely limited throughout their evolutionary history.