Phylogenomic analyses of the genus Drosophila reveals genomic signals of climate adaptation.

Many Drosophila species differ widely in their distributions and climate niches, making them excellent subjects for evolutionary genomic studies. Here, we have developed a database of high-quality assemblies for 46 Drosophila species and one closely related Zaprionus. Fifteen of the genomes were newly sequenced, and 20 were improved with additional sequencing. New or improved annotations were generated for all 47 species, assisted by new transcriptomes for 19. Phylogenomic analyses of these data resolved several previously ambiguous relationships, especially in the melanogaster species group. However, it also revealed significant phylogenetic incongruence among genes, mainly in the form of incomplete lineage sorting in the subgenus Sophophora but also including asymmetric introgression in the subgenus Drosophila. Using the phylogeny as a framework and taking into account these incongruences, we then screened the data for genome-wide signals of adaptation to different climatic niches. First, phylostratigraphy revealed relatively high rates of recent novel gene gain in three temperate pseudoobscura and five desert-adapted cactophilic mulleri subgroup species. Second, we found differing ratios of nonsynonymous to synonymous substitutions in several hundred orthologues between climate generalists and specialists, with trends for significantly higher ratios for those in tropical and lower ratios for those in temperate-continental specialists respectively than those in the climate generalists. Finally, resequencing natural populations of 13 species revealed tropics-restricted species generally had smaller population sizes, lower genome diversity and more deleterious mutations than the more widespread species. We conclude that adaptation to different climates in the genus Drosophila has been associated with large-scale and multifaceted genomic changes.

Application of a bioinformatics training delivery method for reaching dispersed and distant trainees.

Many initiatives have addressed the global need to upskill biologists in bioinformatics tools and techniques. Australia is not unique in its requirement for such training, but due to its large size and relatively small and geographically dispersed population, Australia faces specific challenges. A combined training approach was implemented by the authors to overcome these challenges. The "hybrid" method combines guidance from experienced trainers with the benefits of both webinar-style delivery and concurrent face-to-face hands-on practical exercises in classrooms. Since 2017, the hybrid method has been used to conduct 9 hands-on bioinformatics training sessions at international scale in which over 800 researchers have been trained in diverse topics on a range of software platforms. The method has become a key tool to ensure scalable and more equitable delivery of short-course bioinformatics training across Australia and can be easily adapted to other locations, topics, or settings.

Tracking genetic invasions: Genome-wide single nucleotide polymorphisms reveal the source of pyrethroid-resistant Aedes aegypti (yellow fever mosquito) incursions at international ports.

Biological invasions are increasing globally in number and extent despite efforts to restrict their spread. Knowledge of incursion pathways is necessary to prevent new invasions and to design effective biosecurity protocols at source and recipient locations. This study uses genome-wide single nucleotide polymorphisms (SNPs) to determine the origin of 115 incursive Aedes aegypti(yellow fever mosquito) detected at international ports in Australia and New Zealand. We also genotyped mosquitoes at three point mutations in the voltage-sensitive sodium channel (Vssc) gene: V1016G, F1534C and S989P. These mutations confer knockdown resistance to synthetic pyrethroid insecticides, widely used for controlling invertebrate pests. We first delineated reference populations using Ae. aegypti sampled from 15 locations in Asia, South America, Australia and the Pacific Islands. Incursives were assigned to these populations using discriminant analysis of principal components (DAPC) and an assignment test with a support vector machine predictive model. Bali, Indonesia, was the most common origin of Ae. aegypti detected in Australia, while Ae. aegypti detected in New Zealand originated from Pacific Islands such as Fiji. Most incursives had the same allelic genotype across the three Vsscgene point mutations, which confers strong resistance to synthetic pyrethroids, the only insecticide class used in current, widely implemented aircraft disinsection protocols endorsed by the World Health Organization (WHO). Additionally, all internationally assigned Ae. aegypti had Vssc point mutations linked to pyrethroid resistance that are not found in Australian populations. These findings demonstrate that protocols for preventing introductions of invertebrates must consider insecticide resistance, and highlight the usefulness of genomic data sets for managing global biosecurity objectives.

A Re-Evaluation of Chironomid Deformities as an Environmental Stress Response: Avoiding Survivorship Bias and Testing Noncontaminant Biological Factors.

Larval deformities have been observed in chironomids, and are thought to be associated with aquatic contaminant exposure. However, in laboratory assays, deformities have not been linked with contaminants in the absence of potential confounding variables including mortality, which introduces a survivorship bias. There is also a paucity of data on noncontaminant causes. In addition, power analyses are rarely undertaken, meaning that effect sizes detectable are usually uncertain. We therefore aimed to clarify factors associated with deformities, by running survivorship bias-free (i.e., sublethal) assays, assessing contaminant (copper and imidacloprid) and noncontaminant (malnutrition) stressors, and considering natural biological (metamorphosis) factors in Chironomus tepperi. We included a posteriori power analyses for all tests. Our assays found no significant association between tested factors and deformity rate. Power analyses indicated that the stressor experiment had moderate power to detect deformity effects. The metamorphosis assay had relatively lower power (due to an unexpectedly high control deformity rate), highlighting the importance of power tests in these types of evaluations. These results, in conjunction with others recently published, raise doubts as to the causal effects of environmental stressors on deformity incidence. By avoiding survivorship bias, and by testing noncontaminant factors and statistical power, we present a more holistic methodology, to resolve ongoing uncertainty in this area. We also discuss possible future directions for chironomid deformity research, and concerns regarding survivorship bias in ecotoxicology. Environ Toxicol Chem 2019;38:1658-1667. © 2019 SETAC.

Genomic changes associated with adaptation to arid environments in cactophilic Drosophila species.

BACKGROUND: Insights into the genetic capacities of species to adapt to future climate change can be gained by using comparative genomic and transcriptomic data to reconstruct the genetic changes associated with such adaptations in the past. Here we investigate the genetic changes associated with adaptation to arid environments, specifically climatic extremes and new cactus hosts, through such an analysis of five repleta group Drosophila species. RESULTS: We find disproportionately high rates of gene gains in internal branches in the species' phylogeny where cactus use and subsequently cactus specialisation and high heat and desiccation tolerance evolved. The terminal branch leading to the most heat and desiccation resistant species, Drosophila aldrichi, also shows disproportionately high rates of both gene gains and positive selection. Several Gene Ontology terms related to metabolism were enriched in gene gain events in lineages where cactus use was evolving, while some regulatory and developmental genes were strongly selected in the Drosophila aldrichi branch. Transcriptomic analysis of flies subjected to sublethal heat shocks showed many more downregulation responses to the stress in a heat sensitive versus heat resistant species, confirming the existence of widespread regulatory as well as structural changes in the species' differing adaptations. Gene Ontology terms related to metabolism were enriched in the differentially expressed genes in the resistant species while terms related to stress response were over-represented in the sensitive one. CONCLUSION: Adaptations to new cactus hosts and hot desiccating environments were associated with periods of accelerated evolutionary change in diverse biochemistries. The hundreds of genes involved suggest adaptations of this sort would be difficult to achieve in the timeframes projected for anthropogenic climate change.

Establishing a distributed national research infrastructure providing bioinformatics support to life science researchers in Australia.

EMBL Australia Bioinformatics Resource (EMBL-ABR) is a developing national research infrastructure, providing bioinformatics resources and support to life science and biomedical researchers in Australia. EMBL-ABR comprises 10 geographically distributed national nodes with one coordinating hub, with current funding provided through Bioplatforms Australia and the University of Melbourne for its initial 2-year development phase. The EMBL-ABR mission is to: (1) increase Australia's capacity in bioinformatics and data sciences; (2) contribute to the development of training in bioinformatics skills; (3) showcase Australian data sets at an international level and (4) enable engagement in international programs. The activities of EMBL-ABR are focussed in six key areas, aligning with comparable international initiatives such as ELIXIR, CyVerse and NIH Commons. These key areas-Tools, Data, Standards, Platforms, Compute and Training-are described in this article.

A Wolbachia infection from Drosophila that causes cytoplasmic incompatibility despite low prevalence and densities in males.

Wolbachia bacteria are common insect endosymbionts transmitted maternally and capable of spreading through insect populations by cytoplasmic incompatibility (CI) when infected males cause embryo death after mating with uninfected females. Selection in the Wolbachia endosymbiont occurs on female hosts and is expected to favour strong maternal transmission to female offspring, even at the cost of reduced CI. With maternal leakage, nuclear genes are expected to be selected to suppress cytoplasmic incompatibility caused by males while also reducing any deleterious effects associated with the infection. Here we describe a new type of Wolbachia strain from Drosophila pseudotakahashii likely to have arisen from evolutionary processes on host and/or Wolbachia genomes. This strain is often absent from adult male offspring, but always transmitted to females. It leads to males with low or non-detectable Wolbachia that nevertheless show CI. When detected in adult males, the infection has a low density relative to that in females, a phenomenon not previously seen in Wolbachia infections of Drosophila. This Wolbachia strain is common in natural populations, and shows reduced CI when older (infected) males are crossed. These patterns highlight that endosymbionts can have strong sex-specific effects and that high frequency Wolbachia strains persist through effects on female reproduction. Female-limited Wolbachia infections may be of applied interest if the low level of Wolbachia in males reduces deleterious fitness effects on the host.

Testing the niche-breadth-range-size hypothesis: habitat specialization vs. performance in Australian alpine daisies.

Relatively common species within a clade are expected to perform well across a wider range of conditions than their rarer relatives, yet experimental tests of this "niche-breadth-range-size" hypothesis remain surprisingly scarce. Rarity may arise due to trade-offs between specialization and performance across a wide range of environments. Here we use common garden and reciprocal transplant experiments to test the niche-breadth-range-size hypothesis, focusing on four common and three rare endemic alpine daisies (Brachyscome spp.) from the Australian Alps. We used three experimental contexts: (1) alpine reciprocal seedling experiment, a test of seedling survival and growth in three alpine habitat types differing in environmental quality and species diversity; (2) warm environment common garden, a test of whether common daisy species have higher growth rates and phenotypic plasticity, assessed in a common garden in a warmer climate and run simultaneously with experiment 1; and (3) alpine reciprocal seed experiment, a test of seed germination capacity and viability in the same three alpine habitat types as in experiment 1. In the alpine reciprocal seedling experiment, survival of all species was highest in the open heathland habitat where overall plant diversity is high, suggesting a general, positive response to a relatively productive, low-stress environment. We found only partial support for higher survival of rare species in their habitats of origin. In the warm environment common garden, three common daisies exhibited greater growth and biomass than two rare species, but the other rare species performed as well as the common species. In the alpine reciprocal seed experiment, common daisies exhibited higher germination across most habitats, but rare species maintained a higher proportion of viable seed in all conditions, suggesting different life history strategies. These results indicate that some but not all rare, alpine endemics exhibit stress tolerance at the cost of reduced growth rates in low-stress environments compared to common species. Finally, these findings suggest the seed stage is important in the persistence of rare species, and they provide only weak support at the seedling stage for the niche-breadth-range-size hypothesis.

Four simple recommendations to encourage best practices in research software.

Scientific research relies on computer software, yet software is not always developed following practices that ensure its quality and sustainability. This manuscript does not aim to propose new software development best practices, but rather to provide simple recommendations that encourage the adoption of existing best practices. Software development best practices promote better quality software, and better quality software improves the reproducibility and reusability of research. These recommendations are designed around Open Source values, and provide practical suggestions that contribute to making research software and its source code more discoverable, reusable and transparent. This manuscript is aimed at developers, but also at organisations, projects, journals and funders that can increase the quality and sustainability of research software by encouraging the adoption of these recommendations.

Best practice data life cycle approaches for the life sciences.

Throughout history, the life sciences have been revolutionised by technological advances; in our era this is manifested by advances in instrumentation for data generation, and consequently researchers now routinely handle large amounts of heterogeneous data in digital formats. The simultaneous transitions towards biology as a data science and towards a 'life cycle' view of research data pose new challenges. Researchers face a bewildering landscape of data management requirements, recommendations and regulations, without necessarily being able to access data management training or possessing a clear understanding of practical approaches that can assist in data management in their particular research domain. Here we provide an overview of best practice data life cycle approaches for researchers in the life sciences/bioinformatics space with a particular focus on 'omics' datasets and computer-based data processing and analysis. We discuss the different stages of the data life cycle and provide practical suggestions for useful tools and resources to improve data management practices.

Genomic Trajectories to Desiccation Resistance: Convergence and Divergence Among Replicate Selected Drosophila Lines.

Adaptation to environmental stress is critical for long-term species persistence. With climate change and other anthropogenic stressors compounding natural selective pressures, understanding the nature of adaptation is as important as ever in evolutionary biology. In particular, the number of alternative molecular trajectories available for an organism to reach the same adaptive phenotype remains poorly understood. Here, we investigate this issue in a set of replicated Drosophila melanogaster lines selected for increased desiccation resistance-a classical physiological trait that has been closely linked to Drosophila species distributions. We used pooled whole-genome sequencing (Pool-Seq) to compare the genetic basis of their selection responses, using a matching set of replicated control lines for characterizing laboratory (lab-)adaptation, as well as the original base population. The ratio of effective population size to census size was high over the 21 generations of the experiment at 0.52-0.88 for all selected and control lines. While selected SNPs in replicates of the same treatment (desiccation-selection or lab-adaptation) tended to change frequency in the same direction, suggesting some commonality in the selection response, candidate SNP and gene lists often differed among replicates. Three of the five desiccation-selection replicates showed significant overlap at the gene and network level. All five replicates showed enrichment for ovary-expressed genes, suggesting maternal effects on the selected trait. Divergence between pairs of replicate lines for desiccation-candidate SNPs was greater than between pairs of control lines. This difference also far exceeded the divergence between pairs of replicate lines for neutral SNPs. Overall, while there was overlap in the direction of allele frequency changes and the network and functional categories affected by desiccation selection, replicates showed unique responses at all levels, likely reflecting hitchhiking effects, and highlighting the challenges in identifying candidate genes from these types of experiments when traits are likely to be polygenic.

Tropical Drosophila pandora carry Wolbachia infections causing cytoplasmic incompatibility or male killing.

Wolbachia infections have been described in several Drosophila species, but relatively few have been assessed for phenotypic effects. Cytoplasmic incompatibility (CI) is the most common phenotypic effect that has been detected, while some infections cause male killing or feminization, and many Wolbachia infections have few host effects. Here, we describe two new infections in a recently described species, Drosophila pandora, one of which causes near-complete CI and near-perfect maternal transmission (the "CI" strain). The other infection is a male killer (the "MK" strain), which we confirm by observing reinitiation of male production following tetracycline treatment. No incompatibility was detected in crosses between CI strain males and MK strain females, and rare MK males do not cause CI. Molecular analyses indicate that the CI and MK infections are distantly related and the CI infection is closely related to the wRi infection of Drosophila simulans. Two population surveys indicate that all individuals are infected with Wolbachia, but the MK infection is uncommon. Given patterns of incompatibility among the strains, the infection dynamics is expected to be governed by the relative fitness of the females, suggesting that the CI infection should have a higher fitness. This was evidenced by changes in infection frequencies and sex ratios in population cages initiated at different starting frequencies of the infections.

Cross-Study Comparison Reveals Common Genomic, Network, and Functional Signatures of Desiccation Resistance in Drosophila melanogaster.

Repeated attempts to map the genomic basis of complex traits often yield different outcomes because of the influence of genetic background, gene-by-environment interactions, and/or statistical limitations. However, where repeatability is low at the level of individual genes, overlap often occurs in gene ontology categories, genetic pathways, and interaction networks. Here we report on the genomic overlap for natural desiccation resistance from a Pool-genome-wide association study experiment and a selection experiment in flies collected from the same region in southeastern Australia in different years. We identified over 600 single nucleotide polymorphisms associated with desiccation resistance in flies derived from almost 1,000 wild-caught genotypes, a similar number of loci to that observed in our previous genomic study of selected lines, demonstrating the genetic complexity of this ecologically important trait. By harnessing the power of cross-study comparison, we narrowed the candidates from almost 400 genes in each study to a core set of 45 genes, enriched for stimulus, stress, and defense responses. In addition to gene-level overlap, there was higher order congruence at the network and functional levels, suggesting genetic redundancy in key stress sensing, stress response, immunity, signaling, and gene expression pathways. We also identified variants linked to different molecular aspects of desiccation physiology previously verified from functional experiments. Our approach provides insight into the genomic basis of a complex and ecologically important trait and predicts candidate genetic pathways to explore in multiple genetic backgrounds and related species within a functional framework.

A collection of Australian Drosophila datasets on climate adaptation and species distributions.

The Australian Drosophila Ecology and Evolution Resource (ADEER) collates Australian datasets on drosophilid flies, which are aimed at investigating questions around climate adaptation, species distribution limits and population genetics. Australian drosophilid species are diverse in climatic tolerance, geographic distribution and behaviour. Many species are restricted to the tropics, a few are temperate specialists, and some have broad distributions across climatic regions. Whereas some species show adaptability to climate changes through genetic and plastic changes, other species have limited adaptive capacity. This knowledge has been used to identify traits and genetic polymorphisms involved in climate change adaptation and build predictive models of responses to climate change. ADEER brings together 103 datasets from 39 studies published between 1982-2013 in a single online resource. All datasets can be downloaded freely in full, along with maps and other visualisations. These historical datasets are preserved for future studies, which will be especially useful for assessing climate-related changes over time.

Validation of Pooled Whole-Genome Re-Sequencing in Arabidopsis lyrata.

Sequencing pooled DNA of multiple individuals from a population instead of sequencing individuals separately has become popular due to its cost-effectiveness and simple wet-lab protocol, although some criticism of this approach remains. Here we validated a protocol for pooled whole-genome re-sequencing (Pool-seq) of Arabidopsis lyrata libraries prepared with low amounts of DNA (1.6 ng per individual). The validation was based on comparing single nucleotide polymorphism (SNP) frequencies obtained by pooling with those obtained by individual-based Genotyping By Sequencing (GBS). Furthermore, we investigated the effect of sample number, sequencing depth per individual and variant caller on population SNP frequency estimates. For Pool-seq data, we compared frequency estimates from two SNP callers, VarScan and Snape; the former employs a frequentist SNP calling approach while the latter uses a Bayesian approach. Results revealed concordance correlation coefficients well above 0.8, confirming that Pool-seq is a valid method for acquiring population-level SNP frequency data. Higher accuracy was achieved by pooling more samples (25 compared to 14) and working with higher sequencing depth (4.1× per individual compared to 1.4× per individual), which increased the concordance correlation coefficient to 0.955. The Bayesian-based SNP caller produced somewhat higher concordance correlation coefficients, particularly at low sequencing depth. We recommend pooling at least 25 individuals combined with sequencing at a depth of 100× to produce satisfactory frequency estimates for common SNPs (minor allele frequency above 0.05).

Biological invasions, climate change and genomics.

The rate of biological invasions is expected to increase as the effects of climate change on biological communities become widespread. Climate change enhances habitat disturbance which facilitates the establishment of invasive species, which in turn provides opportunities for hybridization and introgression. These effects influence local biodiversity that can be tracked through genetic and genomic approaches. Metabarcoding and metagenomic approaches provide a way of monitoring some types of communities under climate change for the appearance of invasives. Introgression and hybridization can be followed by the analysis of entire genomes so that rapidly changing areas of the genome are identified and instances of genetic pollution monitored. Genomic markers enable accurate tracking of invasive species' geographic origin well beyond what was previously possible. New genomic tools are promoting fresh insights into classic questions about invading organisms under climate change, such as the role of genetic variation, local adaptation and climate pre-adaptation in successful invasions. These tools are providing managers with often more effective means to identify potential threats, improve surveillance and assess impacts on communities. We provide a framework for the application of genomic techniques within a management context and also indicate some important limitations in what can be achieved.

Limited genetic divergence among Australian alpine Poa tussock grasses coupled with regional structuring points to ongoing gene flow and taxonomic challenges.

BACKGROUND AND AIMS: While molecular approaches can often accurately reconstruct species relationships, taxa that are incompletely differentiated pose a challenge even with extensive data. Such taxa are functionally differentiated, but may be genetically differentiated only at small and/or patchy regions of the genome. This issue is considered here in Poa tussock grass species that dominate grassland and herbfields in the Australian alpine zone. METHODS: Previously reported tetraploidy was confirmed in all species by sequencing seven nuclear regions and five microsatellite markers. A Bayesian approach was used to co-estimate nuclear and chloroplast gene trees with an overall dated species tree. The resulting species tree was used to examine species structure and recent hybridization, and intertaxon fertility was tested by experimental crosses. KEY RESULTS: Species tree estimation revealed Poa gunnii, a Tasmanian endemic species, as sister to the rest of the Australian alpine Poa. The taxa have radiated in the last 0·5-1·2 million years and the non-gunnii taxa are not supported as genetically distinct. Recent hybridization following past species divergence was also not supported. Ongoing gene flow is suggested, with some broad-scale geographic structure within the group. CONCLUSIONS: The Australian alpine Poa species are not genetically distinct despite being distinguishable phenotypically, suggesting recent adaptive divergence with ongoing intertaxon gene flow. This highlights challenges in using conventional molecular taxonomy to infer species relationships in recent, rapid radiations.

A next-generation sequencing method for overcoming the multiple gene copy problem in polyploid phylogenetics, applied to Poa grasses.

BACKGROUND: Polyploidy is important from a phylogenetic perspective because of its immense past impact on evolution and its potential future impact on diversification, survival and adaptation, especially in plants. Molecular population genetics studies of polyploid organisms have been difficult because of problems in sequencing multiple-copy nuclear genes using Sanger sequencing. This paper describes a method for sequencing a barcoded mixture of targeted gene regions using next-generation sequencing methods to overcome these problems. RESULTS: Using 64 3-bp barcodes, we successfully sequenced three chloroplast and two nuclear gene regions (each of which contained two gene copies with up to two alleles per individual) in a total of 60 individuals across 11 species of Australian Poa grasses. This method had high replicability, a low sequencing error rate (after appropriate quality control) and a low rate of missing data. Eighty-eight percent of the 320 gene/individual combinations produced sequence reads, and >80% of individuals produced sufficient reads to detect all four possible nuclear alleles of the homeologous nuclear loci with 95% probability.We applied this method to a group of sympatric Australian alpine Poa species, which we discovered to share an allopolyploid ancestor with a group of American Poa species. All markers revealed extensive allele sharing among the Australian species and so we recommend that the current taxonomy be re-examined. We also detected hypermutation in the trnH-psbA marker, suggesting it should not be used as a land plant barcode region. Some markers indicated differentiation between Tasmanian and mainland samples. Significant positive spatial genetic structure was detected at <100 km with chloroplast but not nuclear markers, which may be a result of restricted seed flow and long-distance pollen flow in this wind-pollinated group. CONCLUSIONS: Our results demonstrate that 454 sequencing of barcoded amplicon mixtures can be used to reliably sample all alleles of homeologous loci in polyploid species and successfully investigate phylogenetic relationships among species, as well as to investigate phylogeographic hypotheses. This next-generation sequencing method is more affordable than and at least as reliable as bacterial cloning. It could be applied to any experiment involving sequencing of amplicon mixtures.