Ecological gradients drive insect wing loss and speciation: The role of the alpine treeline.

Alpine ecosystems are frequently characterized by an abundance of wing-reduced insect species, but the drivers of this biodiversity remain poorly understood. Insect wing reduction in these environments has variously been attributed to altitude, temperature, isolation, habitat stability or decreased habitat size. We used fine-scale ecotypic and genomic analyses, along with broad-scale distributional analyses of ecotypes, to unravel the ecological drivers of wing reduction in the wing-dimorphic stonefly Zelandoperla fenestrata complex. Altitudinal transects within populations revealed dramatic wing reduction over very fine spatial scales, tightly linked to the alpine treeline. Broad biogeographical analyses confirm that the treeline has a much stronger effect on these ecotype distributions than altitude per se. Molecular analyses revealed parallel genomic divergence between vestigial-winged (high altitude) and full-winged (low altitude) ecotypes across distinct streams. These data thus highlight the role of the alpine treeline as a key driver of rapid speciation, providing a new model for ecological diversification along exposure gradients.

Genomics detects population structure within and between ocean basins in a circumpolar seabird: The white-chinned petrel.

The Southern Ocean represents a continuous stretch of circumpolar marine habitat, but the potential physical and ecological drivers of evolutionary genetic differentiation across this vast ecosystem remain unclear. We tested for genetic structure across the full circumpolar range of the white-chinned petrel (Procellaria aequinoctialis) to unravel the potential drivers of population differentiation and test alternative population differentiation hypotheses. Following range-wide comprehensive sampling, we applied genomic (genotyping-by-sequencing or GBS; 60,709 loci) and standard mitochondrial-marker approaches (cytochrome b and first domain of control region) to quantify genetic diversity within and among island populations, test for isolation by distance, and quantify the number of genetic clusters using neutral and outlier (non-neutral) loci. Our results supported the multi-region hypothesis, with a range of analyses showing clear three-region genetic population structure, split by ocean basin, within two evolutionary units. The most significant differentiation between these regions confirmed previous work distinguishing New Zealand and nominate subspecies. Although there was little evidence of structure within the island groups of the Indian or Atlantic oceans, a small set of highly-discriminatory outlier loci could assign petrels to ocean basin and potentially to island group, though the latter needs further verification. Genomic data hold the key to revealing substantial regional genetic structure within wide-ranging circumpolar species previously assumed to be panmictic.

Comparative transcriptomic analysis of a wing-dimorphic stonefly reveals candidate wing loss genes.

BACKGROUND: The genetic basis of wing development has been well characterised for model insect species, but remains poorly understood in phylogenetically divergent, non-model taxa. Wing-polymorphic insect species potentially provide ideal systems for unravelling the genetic basis of secondary wing reduction. Stoneflies (Plecoptera) represent an anciently derived insect assemblage for which the genetic basis of wing polymorphism remains unclear. We undertake quantitative RNA-seq of sympatric full-winged versus vestigial-winged nymphs of a widespread wing-dimorphic New Zealand stonefly, Zelandoperla fenestrata, to identify genes potentially involved in wing development and secondary wing loss. RESULTS: Our analysis reveals substantial differential expression of wing-development genes between full-winged versus vestigial-winged stonefly ecotypes. Specifically, of 23 clusters showing significant similarity to Drosophila wing development-related genes and their pea aphid orthologues, nine were significantly upregulated in full-winged stonefly ecotypes, whereas only one cluster (teashirt) was substantially upregulated in the vestigial-winged ecotype. CONCLUSIONS: These findings suggest remarkable conservation of key wing-development pathways throughout 400 Ma of insect evolution. The finding that two Juvenile Hormone pathway clusters were significantly upregulated in vestigial-winged Zelandoperla supports the hypothesis that Juvenile Hormone may play a key role in modulating insect wing polymorphism, as has previously been suggested for other insect lineages.

The genomic ancestry, landscape genetics and invasion history of introduced mice in New Zealand.

The house mouse (Mus musculus) provides a fascinating system for studying both the genomic basis of reproductive isolation, and the patterns of human-mediated dispersal. New Zealand has a complex history of mouse invasions, and the living descendants of these invaders have genetic ancestry from all three subspecies, although most are primarily descended from M. m. domesticus. We used the GigaMUGA genotyping array (approximately 135 000 loci) to describe the genomic ancestry of 161 mice, sampled from 34 locations from across New Zealand (and one Australian city-Sydney). Of these, two populations, one in the south of the South Island, and one on Chatham Island, showed complete mitochondrial lineage capture, featuring two different lineages of M. m. castaneus mitochondrial DNA but with only M. m. domesticus nuclear ancestry detectable. Mice in the northern and southern parts of the North Island had small traces (approx. 2-3%) of M. m. castaneus nuclear ancestry, and mice in the upper South Island had approximately 7-8% M. m. musculus nuclear ancestry including some Y-chromosomal ancestry-though no detectable M. m. musculus mitochondrial ancestry. This is the most thorough genomic study of introduced populations of house mice yet conducted, and will have relevance to studies of the isolation mechanisms separating subspecies of mice.

Genotyping-by-sequencing supports a genetic basis for wing reduction in an alpine New Zealand stonefly.

Wing polymorphism is a prominent feature of numerous insect groups, but the genomic basis for this diversity remains poorly understood. Wing reduction is a commonly observed trait in many species of stoneflies, particularly in cold or alpine environments. The widespread New Zealand stonefly Zelandoperla fenestrata species group (Z. fenestrata, Z. tillyardi, Z. pennulata) contains populations ranging from fully winged (macropterous) to vestigial-winged (micropterous), with the latter phenotype typically associated with high altitudes. The presence of flightless forms on numerous mountain ranges, separated by lowland fully winged populations, suggests wing reduction has occurred multiple times. We use Genotyping by Sequencing (GBS) to test for genetic differentiation between fully winged (n = 62) and vestigial-winged (n = 34) individuals, sampled from a sympatric population of distinct wing morphotypes, to test for a genetic basis for wing morphology. While we found no population genetic differentiation between these two morphotypes across 6,843 SNP loci, we did detect several outlier loci that strongly differentiated morphotypes across independent tests. These findings indicate that small regions of the genome are likely to be highly differentiated between morphotypes, suggesting a genetic basis for wing reduction. Our results provide a clear basis for ongoing genomic analysis to elucidate critical regulatory pathways for wing development in Pterygota.

Genomic Changes Associated with Reproductive and Migratory Ecotypes in Sockeye Salmon (Oncorhynchus nerka).

Mechanisms underlying adaptive evolution can best be explored using paired populations displaying similar phenotypic divergence, illuminating the genomic changes associated with specific life history traits. Here, we used paired migratory [anadromous vs. resident (kokanee)] and reproductive [shore- vs. stream-spawning] ecotypes of sockeye salmon (Oncorhynchus nerka) sampled from seven lakes and two rivers spanning three catchments (Columbia, Fraser, and Skeena) in British Columbia, Canada to investigate the patterns and processes underlying their divergence. Restriction-site associated DNA sequencing was used to genotype this sampling at 7,347 single nucleotide polymorphisms, 334 of which were identified as outlier loci and candidates for divergent selection within at least one ecotype comparison. Sixty-eight of these outliers were present in two or more comparisons, with 33 detected across multiple catchments. Of particular note, one locus was detected as the most significant outlier between shore and stream-spawning ecotypes in multiple comparisons and across catchments (Columbia, Fraser, and Snake). We also detected several genomic islands of divergence, some shared among comparisons, potentially showing linked signals of differential selection. The single nucleotide polymorphisms and genomic regions identified in our study offer a range of mechanistic hypotheses associated with the genetic basis of O. nerka life history variation and provide novel tools for informing fisheries management.

An ancient selective sweep linked to reproductive life history evolution in sockeye salmon.

Study of parallel (or convergent) phenotypic evolution can provide important insights into processes driving sympatric, ecologically-mediated divergence and speciation, as ecotype pairs may provide a biological replicate of the underlying signals and mechanisms. Here, we provide evidence for a selective sweep creating an island of divergence associated with reproductive behavior in sockeye salmon (Oncorhynchus nerka), identifying a series of linked single nucleotide polymorphisms across a ~22,733 basepair region spanning the leucine-rich repeat-containing protein 9 gene exhibiting signatures of divergent selection associated with stream- and shore-spawning in both anadromous and resident forms across their pan-Pacific distribution. This divergence likely occurred ~3.8 Mya (95% HPD = 2.1-6.03 Mya), after sockeye separated from pink (O. gorbuscha) and chum (O. keta) salmon, but prior to the Pleistocene glaciations. Our results suggest recurrent evolution of reproductive ecotypes across the native range of O. nerka is at least partially associated with divergent selection of pre-existing genetic variation within or linked to this region. As sockeye salmon are unique among Pacific salmonids in their flexibility to spawn in lake-shore benthic environments, this region provides great promise for continued investigation of the genomic basis of O. nerka life history evolution, and, more broadly, for increasing our understanding of the heritable basis of adaptation of complex traits in novel environments.

Sockeye salmon repatriation leads to population re-establishment and rapid introgression with native kokanee.

Re-establishing salmonid populations to areas historically occupied has the substantial potential for conservation gains; however, such interventions also risk negatively impacting native resident stocks. Here, we assessed the success of the hatchery-assisted reintroduction of anadromous sockeye salmon (Oncorhynchus nerka) into Skaha Lake, British Columbia, Canada, and evaluated the genetic consequences for native kokanee, a freshwater-obligate ecotype, using single nucleotide polymorphism genotypic data collected from the reference samples of spawning Okanagan River sockeye and Skaha Lake kokanee presockeye reintroduction, along with annual trawl survey and angler-caught samples obtained over an eight-year period. Significant differentiation was detected between sockeye and kokanee reference samples, with >99% stock assignment. Low proportions of sockeye and hybrids were detected within 2008 and 2010 age-0 trawl samples; however, by 2012, 28% were sockeye, rising to 41% in 2014. The number of hybrids detected rose proportionally with the increase in sockeye and exhibited an intermediate phenotype. Our results indicate that the reintroduction of anadromous sockeye to Skaha Lake is succeeding, with large numbers returning to spawn. However, hybridization with native kokanee is of concern due to the potential for demographic or genetic swamping, with ongoing genetic monitoring necessary to assess the long-term effects of introgression and to support interactive fisheries management.

The complete mitochondrial genomes of two chiton species (Sypharochiton pelliserpentis and Sypharochiton sinclairi) obtained using Illumina next generation sequencing.

Using an Illumina platform, we shot-gun sequenced the complete mitochondrial genomes of two sister chiton species (Sypharochiton pelliserpentis and Sypharochiton sinclairi) to an average coverage of 172× and 60×, respectively. We performed a de novo assembly using SOAPdenovo2 and determined the total mitogenome lengths to be 15,048 and 15,028 bps, respectively. The gene organization was similar to that of other chitons, with 13 protein-coding genes, 24 transfer RNAs and 2 ribosomal RNAs. These data will contribute for resolving the taxonomy and population genetic structures of these species.