Mark F. Richardson

The Drosophila melanogaster Genetic Reference Panel

A major challenge of biology is understanding the relationship between molecular genetic variation and variation in quantitative traits, including fitness. This relationship determines our ability to predict phenotypes from genotypes and to understand how evolutionary forces shape variation within and between species. Previous efforts to dissect the genotype–phenotype map were based on incomplete genotypic information. Here, we describe the Drosophila melanogaster Genetic Reference Panel (DGRP), a community resource for analysis of population genomics and quantitative traits. The DGRP consists of fully sequenced inbred lines derived from a natural population. Population genomic analyses reveal reduced polymorphism in centromeric autosomal regions and the X chromosome, evidence for positive and negative selection, and rapid evolution of the X chromosome. Many variants in novel genes, most at low frequency, are associated with quantitative traits and explain a large fraction of the phenotypic variance. The DGRP facilitates genotype–phenotype mapping using the power of Drosophila genetics.

Space-time wiring specificity supports direction selectivity in the retina

How does the mammalian retina detect motion? This classic problem in visual neuroscience has remained unsolved for 50 years. In search of clues, here we reconstruct Off-type starburst amacrine cells (SACs) and bipolar cells (BCs) in serial electron microscopic images with help from EyeWire, an online community of ‘citizen neuroscientists’. On the basis of quantitative analyses of contact area and branch depth in the retina, we find evidence that one BC type prefers to wire with a SAC dendrite near the SAC soma, whereas another BC type prefers to wire far from the soma. The near type is known to lag the far type in time of visual response. A mathematical model shows how such ‘space–time wiring specificity’ could endow SAC dendrites with receptive fields that are oriented in space–time and therefore respond selectively to stimuli that move in the outward direction from the soma.

Population Genomics of the Wolbachia Endosymbiont in Drosophila melanogaster

Wolbachia are maternally inherited symbiotic bacteria, commonly found in arthropods, which are able to manipulate the reproduction of their host in order to maximise their transmission. The evolutionary history of endosymbionts like Wolbachia can be revealed by integrating information on infection status in natural populations with patterns of sequence variation in Wolbachia and host mitochondrial genomes. Here we use whole-genome resequencing data from 290 lines of Drosophila melanogaster from North America, Europe, and Africa to predict Wolbachia infection status, estimate relative cytoplasmic genome copy number, and reconstruct Wolbachia and mitochondrial genome sequences. Overall, 63% of Drosophila strains were predicted to be infected with Wolbachia by our in silico analysis pipeline, which shows 99% concordance with infection status determined by diagnostic PCR. Complete Wolbachia and mitochondrial genomes show congruent phylogenies, consistent with strict vertical transmission through the maternal cytoplasm and imperfect transmission of Wolbachia. Bayesian phylogenetic analysis reveals that the most recent common ancestor of all Wolbachia and mitochondrial genomes in D. melanogaster dates to around 8,000 years ago. We find evidence for a recent global replacement of ancestral Wolbachia and mtDNA lineages, but our data suggest that the derived wMel lineage arose several thousand years ago, not in the 20th century as previously proposed. Our data also provide evidence that this global replacement event is incomplete and is likely to be one of several similar incomplete replacement events that have occurred since the out-of-Africa migration that allowed D. melanogaster to colonize worldwide habitats. This study provides a complete genomic analysis of the evolutionary mode and temporal dynamics of the D. melanogaster–Wolbachia symbiosis, as well as important resources for further analyses of the impact of Wolbachia on host biology.

A genetic perspective on rapid evolution in cane toads (Rhinella marina)

The process of biological invasion exposes a species to novel pressures, in terms of both the environments it encounters and the evolutionary consequences of range expansion. Several invaders have been shown to exhibit rapid evolutionary changes in response to those pressures, thus providing robust opportunities to clarify the processes at work during rapid phenotypic transitions. The accelerating pace of invasion of cane toads (Rhinella marina) in tropical Australia during its 80‐year history has been well characterized at the phenotypic level, including common‐garden experiments that demonstrate heritability of several dispersal‐relevant traits. Individuals from the invasion front (and their progeny) show distinctive changes in morphology, physiology and behaviour that, in combination, result in far more rapid dispersal than is true of conspecifics from long‐colonized areas. The extensive body of work on cane toad ecology enables us to place into context studies of the genetic basis of these traits. Our analyses of differential gene expression from toads from both ends of this invasion‐history transect reveal substantial upregulation of many genes, notably those involved in metabolism and cellular repair. Clearly, then, the dramatically rapid phenotypic evolution of cane toads in Australia has been accompanied by substantial shifts in gene expression, suggesting that this system is well suited to investigating the genetic underpinnings of invasiveness.

Load carrying during locomotion in the barnacle goose (Branta leucopsis): The effect of load placement and size☆

Load carrying has been used to study the energetics and mechanics of locomotion in a range of taxa. Here we investigated the energetic and kinematic effects of trunk and limb loading in walking barnacle geese (Branta leucopsis). A directly proportional relationship between increasing back-load mass and metabolic rate was established, indicating that the barnacle goose can carry back loads (up to 20% of body mass) more economically than the majority of mammals. The increased cost of supporting and propelling the body during locomotion is likely to account for a major proportion of the extra metabolic cost. Sternal loads up to 15% of body mass were approximately twice as expensive to carry as back loads. Given the key role in dorso-ventral movement of the sternum during respiration we suggest that moving this extra mass may account for the elevated metabolic rate. Loading the distal limb with 5% extra mass incurred the greatest proportional rise in metabolism, and also caused increases in stride length, swing duration and stride frequency during locomotion. The increased work required to move the loaded limb may explain the high cost of walking.

De Novo Assembly and Characterization of the Invasive Northern Pacific Seastar Transcriptome.

Invasive species are a major threat to global biodiversity but can also serve as valuable model systems to examine important evolutionary processes. While the ecological aspects of invasions have been well documented, the genetic basis of adaptive change during the invasion process has been hampered by a lack of genomic resources for the majority of invasive species. Here we report the first larval transcriptomic resource for the Northern Pacific Seastar, Asterias amurensis, an invasive marine predator in Australia. Approximately 117.5 million 100 base-pair (bp) paired-end reads were sequenced from a single RNA-Seq library from a pooled set of full-sibling A. amurensis bipinnaria larvae. We evaluated the efficacy of a pre-assembly error correction pipeline on subsequent de novo assembly. Error correction resulted in small but important improvements to the final assembly in terms of mapping statistics and core eukaryotic genes representation. The error-corrected de novo assembly resulted in 115,654 contigs after redundancy clustering. 41,667 assembled contigs were homologous to sequences from NCBI’s non-redundant protein and UniProt databases. We assigned Gene Ontology, KEGG Orthology, Pfam protein domain terms and predicted protein-coding sequences to > 36,000 contigs. The final transcriptome dataset generated here provides functional information for 18,319 unique proteins, comprising at least 11,355 expressed genes. Furthermore, we identified 9,739 orthologs to P. miniata proteins, evaluated our annotation pipeline and generated a list of 150 candidate genes for responses to several environmental stressors that may be important for adaptation of A. amurensis in the invasive range. Our study has produced a large set of A. amurensis RNA contigs with functional annotations that can serve as a resource for future comparisons to other echinoderm transcriptomes and gene expression studies. Our data can be used to study the genetic basis of adaptive change and other important evolutionary processes during a successful invasion.

Multiple dispersal vectors drive range expansion in an invasive marine species

The establishment and subsequent spread of invasive species is widely recognized as one of the most threatening processes contributing to global biodiversity loss. This is especially true for marine and estuarine ecosystems, which have experienced significant increases in the number of invasive species with the increase in global maritime trade. Understanding the rate and mechanisms of range expansion is therefore of significant interest to ecologists and conservation managers alike. Using a combination of population genetic surveys, environmental DNA (eDNA) plankton sampling and hydrodynamic modelling, we examined the patterns of introduction of the predatory Northern Pacific seastar (Asterias amurensis) and pathways of secondary spread within southeast Australia. Genetic surveys across the invasive range reveal some genetic divergence between the two main invasive regions and no evidence of ongoing gene flow, a pattern that is consistent with the establishment of the second invasive region via a human‐mediated translocation event. In contrast, hydrodynamic modelling combined with eDNA plankton sampling demonstrated that the establishment of range expansion populations within a region is consistent with natural larval dispersal and recruitment. Our results suggest that both anthropogenic and natural dispersal vectors have played an important role in the range expansion of this species in Australia. The multiple modes of spread combined with high levels of fecundity and a long larval duration in A. amurensis suggests it is likely to continue its range expansion and significantly impact Australian marine ecosystems.

propr: An R-package for Identifying Proportionally Abundant Features Using Compositional Data Analysis

In the life sciences, many assays measure only the relative abundances of components in each sample. Such data, called compositional data, require special treatment to avoid misleading conclusions. Awareness of the need for caution in analyzing compositional data is growing, including the understanding that correlation is not appropriate for relative data. Recently, researchers have proposed proportionality as a valid alternative to correlation for calculating pairwise association in relative data. Although the question of how to best measure proportionality remains open, we present here a computationally efficient R package that implements three measures of proportionality. In an effort to advance the understanding and application of proportionality analysis, we review the mathematics behind proportionality, demonstrate its application to genomic data, and discuss some ongoing challenges in the analysis of relative abundance data.

Development of novel microsatellite markers for the invasive Northern Pacific seastar, Asterias amurensis

The Northern Pacific seastar, Asterias amurensis, is a benthic marine predator, which has recently established several invasive populations in Australian waters. To investigate population structure, diversity and patterns of connectivity, we isolated and characterised 27 microsatellite loci and tested their polymorphism based on 46 individuals from two invasive populations. The mean allelic richness was 4.33; observed heterozygosity was 0.42, while the percentage of polymorphic loci was 92.6%. The polymorphic markers will prove useful in the assessment of population genetic parameters, in both invasive and native A. amurensis populations.

Understanding sequencing data as compositions: an outlook and review.

Abstract Motivation Although seldom acknowledged explicitly, count data generated by sequencing platforms exist as compositions for which the abundance of each component (e.g. gene or transcript) is only coherently interpretable relative to other components within that sample. This property arises from the assay technology itself, whereby the number of counts recorded for each sample is constrained by an arbitrary total sum (i.e. library size). Consequently, sequencing data, as compositional data, exist in a non-Euclidean space that, without normalization or transformation, renders invalid many conventional analyses, including distance measures, correlation coefficients and multivariate statistical models. Results The purpose of this review is to summarize the principles of compositional data analysis (CoDA), provide evidence for why sequencing data are compositional, discuss compositionally valid methods available for analyzing sequencing data, and highlight future directions with regard to this field of study. Supplementary information Supplementary data are available at Bioinformatics online.