BioStructures.jl: read, write and manipulate macromolecular structures in Julia.

SUMMARY: Robust, flexible and fast software to read, write and manipulate macromolecular structures is a prerequisite for productively doing structural bioinformatics. We present BioStructures.jl, the first dedicated package in the Julia programming language for dealing with macromolecular structures and the Protein Data Bank. BioStructures.jl builds on the lessons learned with similar packages to provide a large feature set, a flexible object representation and high performance. AVAILABILITY AND IMPLEMENTATION: BioStructures.jl is freely available under the MIT license. Source code and documentation are available at https://github.com/BioJulia/BioStructures.jl. BioStructures.jl is compatible with Julia versions 0.6 and later and is system-independent. CONTACT: j.greener@ucl.ac.uk.

A Sequence Distance Graph framework for genome assembly and analysis.

The Sequence Distance Graph (SDG) framework works with genome assembly graphs and raw data from paired, linked and long reads. It includes a simple deBruijn graph module, and can import graphs using the graphical fragment assembly (GFA) format. It also maps raw reads onto graphs, and provides a Python application programming interface (API) to navigate the graph, access the mapped and raw data and perform interactive or scripted analyses. Its complete workspace can be dumped to and loaded from disk, decoupling mapping from analysis and supporting multi-stage pipelines. We present the design and implementation of the framework, and example analyses scaffolding a short read graph with long reads, and navigating paths in a heterozygous graph for a simulated parent-offspring trio dataset. SDG  is  freely  available  under  the  MIT  license  at https://github.com/bioinfologics/sdg.

Evolutionary genomics of anthroponosis in Cryptosporidium.

Human cryptosporidiosis is the leading protozoan cause of diarrhoeal mortality worldwide, and a preponderance of infections is caused by Cryptosporidium hominis and C. parvum. Both species consist of several subtypes with distinct geographical distributions and host preferences (that is, generalist zoonotic and specialist anthroponotic subtypes). The evolutionary processes that drive the adaptation to the human host and the population structures of Cryptosporidium remain unknown. In this study, we analyse 21 whole-genome sequences to elucidate the evolution of anthroponosis. We show that Cryptosporidium parvum splits into two subclades and that the specialist anthroponotic subtype IIc-a shares a subset of loci with C. hominis that is undergoing rapid convergent evolution driven by positive selection. C. parvum subtype IIc-a also has an elevated level of insertion and deletion mutations in the peri-telomeric genes, which is also a characteristic of other specialist subtypes. Genetic exchange between Cryptosporidium subtypes plays a prominent role throughout the evolution of the genus. Interestingly, recombinant regions are enriched for positively selected genes and potential virulence factors, which indicates adaptive introgression. Analysis of 467 gp60 sequences collected from locations across the world shows that the population genetic structure differs markedly between the main zoonotic subtype (isolation-by-distance) and the anthroponotic subtype (admixed population structure). We also show that introgression between the four anthroponotic Cryptosporidium subtypes and species included in this study has occurred recently, probably within the past millennium.

The ash dieback invasion of Europe was founded by two genetically divergent individuals.

Accelerating international trade and climate change make pathogen spread an increasing concern. Hymenoscyphus fraxineus, the causal agent of ash dieback, is a fungal pathogen that has been moving across continents and hosts from Asian to European ash. Most European common ash trees (Fraxinus excelsior) are highly susceptible to H. fraxineus, although a minority (~5%) have partial resistance to dieback. Here, we assemble and annotate a H. fraxineus draft genome, which approaches chromosome scale. Pathogen genetic diversity across Europe and in Japan, reveals a strong bottleneck in Europe, though a signal of adaptive diversity remains in key host interaction genes. We find that the European population was founded by two divergent haploid individuals. Divergence between these haplotypes represents the ancestral polymorphism within a large source population. Subsequent introduction from this source would greatly increase adaptive potential of the pathogen. Thus, further introgression of H. fraxineus into Europe represents a potential threat and Europe-wide biological security measures are needed to manage this disease.

Evolutionary genetics of immunological supertypes reveals two faces of the Red Queen.

Red Queen host-parasite co-evolution can drive adaptations of immune genes by positive selection that erodes genetic variation (Red Queen arms race) or results in a balanced polymorphism (Red Queen dynamics) and long-term preservation of genetic variation (trans-species polymorphism). These two Red Queen processes are opposite extremes of the co-evolutionary spectrum. Here we show that both Red Queen processes can operate simultaneously by analysing the major histocompatibility complex (MHC) in guppies (Poecilia reticulata and P. obscura) and swamp guppies (Micropoecilia picta). Sub-functionalisation of MHC alleles into 'supertypes' explains how polymorphisms persist during rapid host-parasite co-evolution. Simulations show the maintenance of supertypes as balanced polymorphisms, consistent with Red Queen dynamics, whereas alleles within supertypes are subject to positive selection in a Red Queen arms race. Building on the divergent allele advantage hypothesis, we show that functional aspects of allelic diversity help to elucidate the evolution of polymorphic genes involved in Red Queen co-evolution.

Evolutionary genomics of the cold-adapted diatom Fragilariopsis cylindrus.

The Southern Ocean houses a diverse and productive community of organisms. Unicellular eukaryotic diatoms are the main primary producers in this environment, where photosynthesis is limited by low concentrations of dissolved iron and large seasonal fluctuations in light, temperature and the extent of sea ice. How diatoms have adapted to this extreme environment is largely unknown. Here we present insights into the genome evolution of a cold-adapted diatom from the Southern Ocean, Fragilariopsis cylindrus, based on a comparison with temperate diatoms. We find that approximately 24.7 per cent of the diploid F. cylindrus genome consists of genetic loci with alleles that are highly divergent (15.1 megabases of the total genome size of 61.1 megabases). These divergent alleles were differentially expressed across environmental conditions, including darkness, low iron, freezing, elevated temperature and increased CO2. Alleles with the largest ratio of non-synonymous to synonymous nucleotide substitutions also show the most pronounced condition-dependent expression, suggesting a correlation between diversifying selection and allelic differentiation. Divergent alleles may be involved in adaptation to environmental fluctuations in the Southern Ocean.

Adaptive phenotypic response to climate enabled by epigenetics in a K-strategy species, the fish Leucoraja ocellata (Rajidae).

The relative importance of genetic versus epigenetic changes in adaptive evolution is a hotly debated topic, with studies showing that some species appear to be able to adapt rapidly without significant genetic change. Epigenetic mechanisms may be particularly important for the evolutionary potential of species with long maturation times and low reproductive potential ('K-strategists'), particularly when faced with rapidly changing environmental conditions. Here we study the transcriptome of two populations of the winter skate (Leucoraja ocellata), a typical 'K-strategist', in Atlantic Canada; an endemic population in the southern Gulf of St Lawrence and a large population on the Scotian Shelf. The endemic population has been able to adapt to a 10°C higher water temperature over short evolutionary time (7000 years), dramatically reducing its body size (by 45%) significantly below the minimum maturation size of Scotian Shelf and other populations of winter skate, as well as exhibiting other adaptations in life history and physiology. We demonstrate that the adaptive response to selection has an epigenetic basis, cataloguing 3653 changes in gene expression that may have enabled this species to rapidly respond to the novel environment. We argue that the epigenetic augmentation of species evolutionary potential (its regulation though gene expression) can enable K-strategists to survive and adapt to different environments, and this mechanism may be particularly important for the persistence of sharks, skates and rays in the light of future climate change.

HYBRIDCHECK: software for the rapid detection, visualization and dating of recombinant regions in genome sequence data.

HYBRIDCHECK is a software package to visualize the recombination signal in large DNA sequence data set, and it can be used to analyse recombination, genetic introgression, hybridization and horizontal gene transfer. It can scan large (multiple kb) contigs and whole-genome sequences of three or more individuals. HYBRIDCHECK is written in the r software for OS X, Linux and Windows operating systems, and it has a simple graphical user interface. In addition, the r code can be readily incorporated in scripts and analysis pipelines. HYBRIDCHECK implements several ABBA-BABA tests and visualizes the effects of hybridization and the resulting mosaic-like genome structure in high-density graphics. The package also reports the following: (i) the breakpoint positions, (ii) the number of mutations in each introgressed block, (iii) the probability that the identified region is not caused by recombination and (iv) the estimated age of each recombination event. The divergence times between the donor and recombinant sequence are calculated using a JC, K80, F81, HKY or GTR correction, and the dating algorithm is exceedingly fast. By estimating the coalescence time of introgressed blocks, it is possible to distinguish between hybridization and incomplete lineage sorting. HYBRIDCHECK is libré software and it and its manual are free to download from http://ward9250.github.io/HybridCheck/.

Evidence for suppression of immunity as a driver for genomic introgressions and host range expansion in races of Albugo candida, a generalist parasite.

How generalist parasites with wide host ranges can evolve is a central question in parasite evolution. Albugo candida is an obligate biotrophic parasite that consists of many physiological races that each specialize on distinct Brassicaceae host species. By analyzing genome sequence assemblies of five isolates, we show they represent three races that are genetically diverged by ∼1%. Despite this divergence, their genomes are mosaic-like, with ∼25% being introgressed from other races. Sequential infection experiments show that infection by adapted races enables subsequent infection of hosts by normally non-infecting races. This facilitates introgression and the exchange of effector repertoires, and may enable the evolution of novel races that can undergo clonal population expansion on new hosts. We discuss recent studies on hybridization in other eukaryotes such as yeast, Heliconius butterflies, Darwin's finches, sunflowers and cichlid fishes, and the implications of introgression for pathogen evolution in an agro-ecological environment.