The genomic causes and effects of divergent ecological selection during speciation are still poorly understood. Here, we report the discovery and detailed characterization of early-stage adaptive divergence of two cichlid fish ecomorphs in a small (700m diameter) isolated crater lake in Tanzania. The ecomorphs differ in depth preference, male breeding color, body shape, diet and trophic morphology. With whole genome sequences of 146 fish, we identify 98 clearly demarcated genomic 'islands' of high differentiation and demonstrate association of genotypes across these islands to divergent mate preferences. The islands contain candidate adaptive genes enriched for functions in sensory perception (including rhodopsin and other twilight vision associated genes), hormone signaling and morphogenesis. Our study suggests mechanisms and genomic regions that may play a role in the closely related mega-radiation of Lake Malawi. * Publisher's Disclaimer: This manuscript has been accepted for publication in Science. This version has not undergone final editing.Please refer to the complete version of record at http://www.sciencemag.org/. The manuscript may not be reproduced or used in any manner that does not fall within the fair use provisions of the

We used 20 de novo genome assemblies to probe the speciation history and architecture of gene flow in rapidly radiating Heliconius butterflies. Our tests to distinguish incomplete lineage sorting from introgression indicate that gene flow has obscured several ancient phylogenetic relationships in this group over large swathes of the genome. Introgressed loci are underrepresented in low-recombination and gene-rich regions, consistent with the purging of foreign alleles more tightly linked to incompatibility loci. Here, we identify a hitherto unknown inversion that traps a color pattern switch locus. We infer that this inversion was transferred between lineages by introgression and is convergent with a similar rearrangement in another part of the genus. These multiple de novo genome sequences enable improved understanding of the importance of introgression and selective processes in adaptive radiation.

Reconstruction of target genomes from sequence data produced by instruments that are agnostic as to the species-of-origin may be confounded by contaminant DNA. Whether introduced during sample processing or through co-extraction alongside the target DNA, if insufficient care is taken during the assembly process, the final assembled genome may be a mixture of data from several species. Such assemblies can confound sequence-based biological inference and, when deposited in public databases, may be included in downstream analyses by users unaware of underlying problems. We present BlobToolKit, a software suite to aid researchers in identifying and isolating non-target data in draft and publicly available genome assemblies. BlobToolKit can be used to process assembly, read and analysis files for fully reproducible interactive exploration in the browser-based Viewer. BlobToolKit can be used during assembly to filter non-target DNA, helping researchers produce assemblies with high biological credibility. We have been running an automated BlobToolKit pipeline on eukaryotic assemblies publicly available in the International Nucleotide Sequence Data Collaboration and are making the results available through a public instance of the Viewer at https://blobtoolkit.genomehubs.org/view. We aim to complete analysis of all publicly available genomes and then maintain currency with the flow of new genomes. We have worked to embed these views into the presentation of genome assemblies at the European Nucleotide Archive, providing an indication of assembly quality alongside the public record with links out to allow full exploration in the Viewer.

The Plasmodium falciparum variant erythrocyte surface antigens known as PfEMP1, encoded by the var gene family, are thought to play a crucial role in malaria pathogenesis because they mediate adhesion to host cells and immuno-modulation. Var genes have been divided into three major groups (A, B and C) and two intermediate groups (B/A and B/C) on the basis of their genomic location and upstream sequence. We analysed expressed sequence tags of the var gene DBLα domain to investigate var gene transcription in relation to disease severity in Malian children. We found that P. falciparum isolates from children with cerebral malaria (unrousable coma) predominantly transcribe var genes with DBLα1-like domains that are characteristic of Group A or B/A var genes. In contrast, isolates from children with equally high parasite burdens but no symptoms or signs of severe malaria (hyperparasitaemia patients) predominantly transcribe var genes with DBLα0-like domains that are characteristic of the B and C-related var gene groups. These results suggest that var genes with DBLα1-like domains (Group A or B/A) may be implicated in the pathogenesis of cerebral malaria, while var genes with DBLα0-like domains promote less virulent malaria infections.

Comparing newly obtained and previously known nucleotide and amino-acid sequences underpins modern biological research. BLAST is a well-established tool for such comparisons but is challenging to use on new data sets. We combined a user-centric design philosophy with sustainable software development approaches to create Sequenceserver, a tool for running BLAST and visually inspecting BLAST results for biological interpretation. Sequenceserver uses simple algorithms to prevent potential analysis errors and provides flexible text-based and visual outputs to support researcher productivity. Our software can be rapidly installed for use by individuals or on shared servers.

Previous studies of highly branched mutants in pea (rms1–rms5), Arabidopsis thaliana (max1–max4), petunia (dad1–dad3), and rice (d3, d10, htd1/d17, d14, d27) identified strigolactones or their derivates (SLs), as shoot branching inhibitors. This recent discovery offers the possibility of using SLs to regulate branching commercially, for example, in chrysanthemum, an important cut flower crop. To investigate this option, SL physiology and molecular biology were studied in chrysanthemum (Dendranthema grandiflorum), focusing on the CCD8/MAX4/DAD1/RMS1/D10 gene. Our results suggest that, as has been proposed for Arabidopsis, the ability of SLs to inhibit bud activity depends on the presence of a competing auxin source. The chrysanthemum SL biosynthesis gene, CCD8 was cloned, and found to be regulated in a similar, but not identical way to known CCD8s. Expression analyses revealed that DgCCD8 is predominantly expressed in roots and stems, and is up-regulated by exogenous auxin. Exogenous SL can down-regulate DgCCD8 expression, but this effect can be overridden by apical auxin application. This study provides evidence that SLs are promising candidates to alter the shoot branching habit of chrysanthemum.

Human dispersal of organisms is an important process modifying natural patterns of biodiversity. Such dispersal generates new patterns of genetic diversity that overlie natural phylogeographical signatures, allowing discrimination between alternative dispersal mechanisms. Here we use allele frequency and DNA sequence data to distinguish between alternative scenarios (unassisted range expansion and long range introduction) for the colonization of northern Europe by an oak-feeding gallwasp, Andricus kollari. Native to Mediterranean latitudes from Portugal to Iran, this species became established in northern Europe following human introduction of a host plant, the Turkey oak Quercus cerris. Colonization of northern Europe is possible through three alternative routes: (i) unassisted range expansion from natural populations in the Iberian Peninsula; (ii) unassisted range expansion from natural populations in Italy and Hungary; or (iii) descent from populations imported to the UK as trade goods from the eastern Mediterranean in the 1830s. We show that while populations in France were colonized from sources in Italy and Hungary, populations in the UK and neighbouring parts of coastal northern Europe encompass allozyme and sequence variation absent from the known native range. Further, these populations show demographic signatures expected for large stable populations, rather than signatures of rapid population growth from small numbers of founders. The extent and spatial distribution of genetic diversity in the UK suggests that these A. kollari populations are derived from introductions of large numbers of individuals from each of two genetically divergent centres of diversity in the eastern Mediterranean. The strong spatial patterning in genetic diversity observed between different regions of northern Europe, and between sites in the UK, is compatible with leptokurtic models of population establishment.

Strigolactones (SLs) are carotenoid-derived phytohormones with diverse roles. They are secreted from roots as attractants for arbuscular mycorrhizal fungi and have a wide range of endogenous functions, such as regulation of root and shoot system architecture. To date, six genes associated with SL synthesis and signaling have been molecularly identified using the shootbranching mutants more axillary growth (max) of Arabidopsis (Arabidopsis thaliana) and dwarf (d) of rice (Oryza sativa). Here, we present a phylogenetic analysis of the MAX/D genes to clarify the relationships of each gene with its wider family and to allow the correlation of events in the evolution of the genes with the evolution of SL function. Our analysis suggests that the notion of a distinct SL pathway is inappropriate. Instead, there may be a diversity of SL-like compounds, the response to which requires a D14/D14-like protein. This ancestral system could have been refined toward distinct ligand-specific pathways channeled through MAX2, the most downstream known component of SL signaling. MAX2 is tightly conserved among land plants and is more diverged from its nearest sister clade than any other SL-related gene, suggesting a pivotal role in the evolution of SL signaling. By contrast, the evidence suggests much greater flexibility upstream of MAX2. The MAX1 gene is a particularly strong candidate for contributing to diversification of inputs upstream of MAX2. Our functional analysis of the MAX1 family demonstrates the early origin of its catalytic function and both redundancy and functional diversification associated with its duplication in angiosperm lineages.