Special features of RAD Sequencing data:implications for genotyping

Restriction site-associated DNA Sequencing (RAD-Seq) is an economical and efficient method for SNP discovery and genotyping. As with other sequencing-by-synthesis methods, RAD-Seq produces stochastic count data and requires sensitive analysis to develop or genotype markers accurately. We show that there are several sources of bias specific to RAD-Seq that are not explicitly addressed by current genotyping tools, namely restriction fragment bias, restriction site heterozygosity and PCR GC content bias. We explore the performance of existing analysis tools given these biases and discuss approaches to limiting or handling biases in RAD-Seq data. While these biases need to be taken seriously, we believe RAD loci affected by them can be excluded or processed with relative ease in most cases and that most RAD loci will be accurately genotyped by existing tools

Recombination rate variation shapes barriers to introgression across butterfly genomes.

Hybridisation and introgression can dramatically alter the relationships among groups of species, leading to phylogenetic discordance across the genome and between populations. Introgression can also erode species differences over time, but selection against introgression at certain loci acts to maintain postmating species barriers. Theory predicts that species barriers made up of many loci throughout the genome should lead to a broad correlation between introgression and recombination rate, which determines the extent to which selection on deleterious foreign alleles will affect neutral alleles at physically linked loci. Here, we describe the variation in genealogical relationships across the genome among three species of Heliconius butterflies: H. melpomene (mel), H. cydno (cyd), and H. timareta (tim), using whole genomes of 92 individuals, and ask whether this variation can be explained by heterogeneous barriers to introgression. We find that species relationships vary predictably at the chromosomal scale. By quantifying recombination rate and admixture proportions, we then show that rates of introgression are predicted by variation in recombination rate. This implies that species barriers are highly polygenic, with selection acting against introgressed alleles across most of the genome. In addition, long chromosomes, which have lower recombination rates, produce stronger barriers on average than short chromosomes. Finally, we find a consistent difference between two species pairs on either side of the Andes, which suggests differences in the architecture of the species barriers. Our findings illustrate how the combined effects of hybridisation, recombination, and natural selection, acting at multitudes of loci over long periods, can dramatically sculpt the phylogenetic relationships among species

A conserved set of maternal genes? Insights from a molluscan transcriptome

The early animal embryo is entirely reliant on maternal gene products for a ‘jump-start’ that transforms a transcriptionally inactive embryo into a fully functioning zygote. Despite extensive work on model species, it has not been possible to perform a comprehensive comparison of maternally-provisioned transcripts across the Bilateria because of the absence of a suitable dataset from the Lophotrochozoa. As part of an ongoing effort to identify the maternal gene that determines left-right asymmetry in snails, we have generated transcriptome data from 1 to 2-cell and ~32-cell pond snail (Lymnaea stagnalis) embryos. Here, we compare these data to maternal transcript datasets from other bilaterian metazoan groups, including representatives of the Ecydysozoa and Deuterostomia. We found that between 5 and 10% of all L. stagnalis maternal transcripts (~300-400 genes) are also present in the equivalent arthropod (Drosophila melanogaster), nematode (Caenorhabditis elegans), urochordate (Ciona intestinalis) and chordate (Homo sapiens, Mus musculus, Danio rerio) datasets. While the majority of these conserved maternal transcripts (“COMATs”) have housekeeping gene functions, they are a non-random subset of all housekeeping genes, with an overrepresentation of functions associated with nucleotide binding, protein degradation and activities associated with the cell cycle. We conclude that a conserved set of maternal transcripts and their associated functions may be a necessary starting point of early development in the Bilateria. For the wider community interested in discovering conservation of gene expression in early bilaterian development, the list of putative COMATs may be useful resource

How do predators generalize warning signals in simple and complex prey communities? Insights from a videogame

The persistence of distinct warning signals within and between sympatric mimetic communities is a puzzling evolutionary question because selection favours convergence of colour patterns among toxic species. Such convergence is partly shaped by predators' reaction to similar but not identical stimulus (i.e. generalization behaviour), and generalization by predators is likely to be shaped by the diversity of local prey. However, studying generalization behaviour is generally limited to simple variations of prey colour patterns. Here, we used a computer game played by humans as surrogate predators to investigate generalization behaviours in simple (4 morphs) and complex (10 morphs) communities of unprofitable (associated with a penalty) and profitable butterflies. Colour patterns used in the game are observed in the natural populations of unprofitable butterfly species such as Heliconius numata. Analyses of 449 game participants' behaviours show that players avoided unprofitable prey more readily in simple than in complex communities. However, generalization was observed only in players that faced complex communities, enhancing the protection of profitable prey that looked similar to at least one unprofitable morph. Additionally, similarity among unprofitable prey also reduced attack rates only in complex communities. These results are consistent with previous studies using avian predators but artificial colour patterns and suggest that mimicry is more likely to evolve in complex communities where increases in similarity are more likely to be advantageous

Ancestral Admixture Is the Main Determinant of Global Biodiversity in Fission Yeast

Mutation and recombination are key evolutionary processes governing phenotypic variation and reproductive isolation. We here demonstrate that biodiversity within all globally known strains of Schizosaccharomyces pombe arose through admixture between two divergent ancestral lineages. Initial hybridization was inferred to have occurred ∼20-60 sexual outcrossing generations ago consistent with recent, human-induced migration at the onset of intensified transcontinental trade. Species-wide heritable phenotypic variation was explained near-exclusively by strain-specific arrangements of alternating ancestry components with evidence for transgressive segregation. Reproductive compatibility between strains was likewise predicted by the degree of shared ancestry. To assess the genetic determinants of ancestry block distribution across the genome, we characterized the type, frequency, and position of structural genomic variation using nanopore and single-molecule real-time sequencing. Despite being associated with double-strand break initiation points, over 800 segregating structural variants exerted overall little influence on the introgression landscape or on reproductive compatibility between strains. In contrast, we found strong ancestry disequilibrium consistent with negative epistatic selection shaping genomic ancestry combinations during the course of hybridization. This study provides a detailed, experimentally tractable example that genomes of natural populations are mosaics reflecting different evolutionary histories. Exploiting genome-wide heterogeneity in the history of ancestral recombination and lineage-specific mutations sheds new light on the population history of S. pombe and highlights the importance of hybridization as a creative force in generating biodiversity

Genome-wide genetic marker discovery and genotyping using next-generation sequencing,”

Abstract | The advent of next-generation sequencing (NGS) has revolutionized genomic and transcriptomic approaches to biology. These new sequencing tools are also valuable for the discovery, validation and assessment of genetic markers in populations. Here we review and discuss best practices for several NGS methods for genome-wide genetic marker development and genotyping that use restriction enzyme digestion of target genomes to reduce the complexity of the target. These new methods -which include reduced-representation sequencing using reduced-representation libraries (RRLs) or complexity reduction of polymorphic sequences (CRoPS), restriction-site-associated DNA sequencing (RAD-seq) and low coverage genotyping -are applicable to both model organisms with high-quality reference genome sequences and, excitingly, to non-model species with no existing genomic data

Frequent use of paracetamol and risk of allergic disease among women in an Ethiopian population

Introduction The hypothesis that paracetamol might increase the risk of asthma and other allergic diseases have gained support from a range of independent studies. However, in studies based in developed countries, the possibility that paracetamol and asthma are associated through aspirin avoidance is difficult to exclude. Objectives To explore this hypothesis among women in a developing country, where we have previously reported aspirin avoidance to be rare. Methods In 2005/6 a population based cohort of 1065 pregnant women was established in Butajira, Ethiopia and baseline demographic data collected. At 3 years post birth, an interview-based questionnaire administered to 945 (94%) of these women collected data on asthma, eczema, and hay fever in the past 12 month, frequency of paracetamol use and potential confounders. Allergen skin tests to Dermatophagoides pteronyssinus and cockroach were also performed. The independent effects of paracetamol use on allergic outcomes were determined using multiple logistic regression analysis. Findings The prevalence of asthma, eczema and hay fever was 1.7%, 0.9% and 3.8% respectively; of any one of these conditions 5.5%, and of allergen sensitization 7.8%. Paracetamol use in the past month was reported by 29%, and associations of borderline significance were seen for eczema (adjusted OR (95% CI) = 8.51 (1.68 to 43.19) for 1–3 tablets and 2.19 (0.36 to 13.38) for ≥4 tablets, compared to no tablets in the past month; overall p = 0.055) and for ‘any allergic condition’ (adjusted OR (95% CI) = 2.73 (1.22 to 6.11) for 1–3 tablets and 1.35 (0.67 to 2.70) for ≥4 tablets compared to 0 in the past month; overall p = 0.071). Conclusions This study provides further cross-sectional evidence that paracetamol use increases the risk of allergic disease

Major improvements to the Heliconius melpomene genome assembly used to confirm 10 chromosome fusion events in 6 million years of butterfly evolution

The Heliconius butterflies are a widely studied adaptive radiation of 46 species spread across Central and South America, several of which are known to hybridize in the wild. Here, we present a substantially improved assembly of the Heliconius melpomene genome, developed using novel methods that should be applicable to improving other genome assemblies produced using short read sequencing. First, we whole-genome-sequenced a pedigree to produce a linkage map incorporating 99% of the genome. Second, we incorporated haplotype scaffolds extensively to produce a more complete haploid version of the draft genome. Third, we incorporated ~20x coverage of Pacific Biosciences sequencing, and scaffolded the haploid genome using an assembly of this long-read sequence. These improvements result in a genome of 795 scaffolds, 275 Mb in length, with an N50 length of 2.1 Mb, an N50 number of 34, and with 99% of the genome placed, and 84% anchored on chromosomes. We use the new genome assembly to confirm that the Heliconius genome underwent 10 chromosome fusions since the split with its sister genus Eueides, over a period of about 6 million yr

Variation in the SLC23A1 gene does not influence cardiometabolic outcomes to the extent expected given its association with L-ascorbic acid

Background: Observational studies showed that circulating l-ascorbic acid (vitamin C) is inversely associated with cardiometabolic traits. However, these studies were susceptible to confounding and reverse causation. Objectives:We assessed the relation between l-ascorbic acid and 10 cardiometabolic traits by using a single nucleotide polymorphism in the solute carrier family 23 member 1 (SLC23A1) gene (rs33972313) associated with circulating l-ascorbic acid concentrations. The observed association between rs33972313 and cardiometabolic outcomes was compared with that expected given the rs33972313-l-ascorbic acid and l-ascorbic acid–outcome associations. Design: A meta-analysis was performed in the following 5 independent studies: the British Women's Heart and Health Study (n = 1833), the MIDSPAN study (n = 1138), the Ten Towns study (n = 1324), the British Regional Heart Study (n = 2521), and the European Prospective Investigation into Cancer (n = 3737). Results: With the use of a meta-analysis of observational estimates, inverse associations were shown between l-ascorbic acid and systolic blood pressure, triglycerides, and the waist-hip ratio [the strongest of which was the waist-hip ratio (−0.13-SD change; 95% CI: −0.20-, −0.07-SD change; P = 0.0001) per SD increase in l-ascorbic acid], and a positive association was shown with high-density lipoprotein (HDL) cholesterol. The variation at rs33972313 was associated with a 0.18-SD (95% CI: 0.10-, 0.25-SD; P = 3.34 × 10−6) increase in l-ascorbic acid per effect allele. There was no evidence of a relation between the variation at rs33972313 and any cardiometabolic outcome. Although observed estimates were not statistically different from expected associations between rs33972313 and cardiometabolic outcomes, estimates for low-density lipoprotein cholesterol, HDL cholesterol, triglycerides, glucose, and body mass index were in the opposite direction to those expected. Conclusions: The nature of the genetic association exploited in this study led to limited statistical application, but despite this, when all cardiometabolic traits were assessed, there was no evidence of any trend supporting a protective role of l-ascorbic acid. In the context of existing work, these results add to the suggestion that observational relations between l-ascorbic acid and cardiometabolic health may be attributable to confounding and reverse causation

The Far-Ultraviolet Spectrum and Short Timescale Variability of AM Herculis from Observations with the Hopkins Ultraviolet Telescope

Using the Hopkins Ultraviolet Telescope (HUT), we have obtained 850-1850 angstrom spectra of the magnetic cataclysmic variable star AM Her in the high state. These observations provide high time resolution spectra of AM Her in the FUV and sample much of the orbital period of the system. The spectra are not well-modelled in terms of simple white dwarf (WD) atmospheres, especially at wavelengths shortward of Lyman alpha. The continuum flux changes by a factor of 2 near the Lyman limit as a function of orbital phase; the peak fluxes are observed near magnetic phase 0.6 when the accreting pole of the WD is most clearly visible. The spectrum of the hotspot can be modelled in terms of a 100 000 K WD atmosphere covering 2% of the WD surface. The high time resolution of the HUT data allows an analysis of the short term variability and shows the UV luminosity to change by as much as 50% on timescales as short as 10 s. This rapid variability is shown to be inconsistent with the clumpy accretion model proposed to account for the soft X-ray excess in polars. We see an increase in narrow line emission during these flares when the heated face of the secondary is in view. The He II narrow line flux is partially eclipsed at secondary conjunction, implying that the inclination of the system is greater than 45 degrees. We also present results from models of the heated face of the secondary. These models show that reprocessing on the face of the secondary star of X-ray/EUV emission from the accretion region near the WD can account for the intensities and kinematics of most of the narrow line components observed.Comment: 19 pp., 12 fig., 3 tbl. To appear in The Astrophysical Journal. Also available at http://greeley.pha.jhu.edu/papers/amherpp.ps.g