PoPoolation: a toolbox for population genetic analysis of next generation sequencing data from pooled individuals

Exploring the genetics of trotting racing ability in horses using a unique Nordic horse model

Background

Horses have been strongly selected for speed, strength, and endurance-exercise traits since the onset of domestication. As a result, highly specialized horse breeds have developed with many modern horse breeds often representing closed populations with high phenotypic and genetic uniformity. However, a great deal of variation still exists between breeds, making the horse particularly well suited for genetic studies of athleticism. To identify genomic regions associated with athleticism as it pertains to trotting racing ability in the horse, the current study applies a pooled sequence analysis approach using a unique Nordic horse model.

Results

Pooled sequence data from three Nordic horse populations were used for F_ST analysis. After strict filtering, F_ST analysis yielded 580 differentiated regions for trotting racing ability. Candidate regions on equine chromosomes 7 and 11 contained the largest number of SNPs (n = 214 and 147, respectively). GO analyses identified multiple genes related to intelligence, energy metabolism, and skeletal development as potential candidate genes. However, only one candidate region for trotting racing ability overlapped a known racing ability QTL.

Conclusions

Not unexpected for genomic investigations of complex traits, the current study identified hundreds of candidate regions contributing to trotting racing ability in the horse. Likely resulting from the cumulative effects of many variants across the genome, racing ability continues to demonstrate its polygenic nature with candidate regions implicating genes influencing both musculature and neurological development.

Electronic supplementary material

The online version of this article (10.1186/s12864-019-5484-9) contains supplementary material, which is available to authorized users.

Small and genetically highly structured populations in a long-legged bee, Rediviva longimanus, as inferred by pooled RAD-seq

Adaptation to local host plants may impact a pollinator's population genetic structure by reducing gene flow and driving population genetic differentiation, representing an early stage of ecological speciation. South African Rediviva longimanus bees exhibit elongated forelegs, a bizarre adaptation for collecting oil from floral spurs of their Diascia hosts. Furthermore, R. longimanus foreleg length (FLL) differs significantly among populations, which has been hypothesised to result from selection imposed by inter-population variation in Diascia floral spur length. Here, we used a pooled restriction site-associated DNA sequencing (pooled RAD-seq) approach to investigate the population genetic structure of R. longimanus and to test if phenotypic differences in FLL translate into increased genetic differentiation (i) between R. longimanus populations and (ii) between phenotypes across populations. We also inferred the effects of demographic processes on population genetic structure and tested for genetic markers underpinning local adaptation. RESULTS: Populations showed marked genetic differentiation (average FST = 0.165), though differentiation was not statistically associated with differences between populations in FLL. All populations exhibited very low genetic diversity and were inferred to have gone through recent bottleneck events, suggesting extremely low effective population sizes. Genetic differentiation between samples pooled by leg length (short versus long) rather than by population of origin was even higher (FST = 0.260) than between populations, suggesting reduced interbreeding between long and short-legged individuals. Signatures of selection were detected in 1119 (3.8%) of a total of 29,721 SNP markers, CONCLUSIONS: Populations of R. longimanus appear to be small, bottlenecked and isolated. Though we could not detect the effect of local adaptation (FLL in response to floral spurs of host plants) on population genetic differentiation, short and long legged bees appeared to be partially differentiated, suggesting incipient ecological speciation. To test this hypothesis, greater resolution through the use of individual-based whole-genome analyses is now needed to quantify the degree of reproductive isolation between long and short legged bees between and even within populations.

Extreme copy number variation at a tRNA ligase gene affecting phenology and fitness in yellow monkeyflowers

Copy number variation (CNV) is a major part of the genetic diversity segregating within populations, but remains poorly understood relative to single nucleotide variation. Here, we report on a tRNA ligase gene (Migut.N02091; RLG1a) exhibiting unprecedented, and fitness-relevant, CNV within an annual population of the yellow monkeyflower Mimulus guttatus. RLG1a variation was associated with multiple traits in pooled population sequencing (PoolSeq) scans of phenotypic and phenological cohorts. Resequencing of inbred lines revealed intermediate-frequency three-copy variants of RLG1a (trip+; 5/35 = 14%), and trip+ lines exhibited elevated RLG1a expression under multiple conditions. trip+ carriers, in addition to being over-represented in late-flowering and large-flowered PoolSeq populations, flowered later under stressful conditions in a greenhouse experiment (p < 0.05). In wild population samples, we discovered an additional rare RLG1a variant (high+) that carries 250-300 copies of RLG1a totalling ~5.7 Mb (20-40% of a chromosome). In the progeny of a high+ carrier, Mendelian segregation of diagnostic alleles and qPCR-based copy counts indicate that high+ is a single tandem array unlinked to the single-copy RLG1a locus. In the wild, high+ carriers had highest fitness in two particularly dry and/or hot years (2015 and 2017; both p < 0.01), while single-copy individuals were twice as fecund as either CNV type in a lush year (2016: p < 0.005). Our results demonstrate fluctuating selection on CNVs affecting phenological traits in a wild population, suggest that plant tRNA ligases mediate stress-responsive life-history traits, and introduce a novel system for investigating the molecular mechanisms of gene amplification.

Mutation of a serine near the catalytic site of the choline acetyltransferase a gene almost completely abolishes motility of the zebrafish embryo

In zebrafish, the gene choline acetyltransferase a (chata) encodes one of the two ChAT orthologs responsible for the synthesis of acetylcholine. Acetylcholine (ACh) is essential for neuromuscular transmission and its impaired synthesis by ChAT can lead to neuromuscular junction disorders such as congenital myasthenic syndromes in humans. We have identified a novel mutation in the chata gene of zebrafish, chata^tk64, in a collection of uncharacterised ENU-induced mutants. This mutant carries a missense mutation in the codon of a highly conserved serine changing it to an arginine (S102R). This serine is conserved among ChATs from zebrafish, rat, mice and chicken to humans. It resides within the catalytic domain and in the vicinity of the active site of the enzyme. However, it has not been reported so far to be required for enzymatic activity. Modelling of the S102R variant change in the ChAT protein crystal structure suggests that the change affects protein structure and has a direct impact on the catalytic domain of the protein which abolishes embryo motility almost completely.

Mapping and Sequencing of a Significant Quantitative Trait Locus Affecting Resistance to Koi Herpesvirus in Common Carp

Cyprinids are the most highly produced group of fishes globally, with common carp being one of the most valuable species of the group. Koi herpesvirus (KHV) infections can result in high levels of mortality, causing major economic losses, and is listed as a notifiable disease by the World Organization for Animal Health. Selective breeding for host resistance has the potential to reduce morbidity and losses due to KHV. Therefore, improving knowledge about host resistance and methods of incorporating genomic data into breeding for resistance may contribute to a decrease in economic losses in carp farming. In the current study, a population of 1,425 carp juveniles, originating from a factorial cross between 40 sires and 20 dams was challenged with KHV. Mortalities and survivors were recorded and sampled for genotyping by sequencing using Restriction Site-Associated DNA sequencing (RADseq). Genome-wide association analyses were performed to investigate the genetic architecture of resistance to KHV. A genome-wide significant QTL affecting resistance to KHV was identified on linkage group 44, explaining approximately 7% of the additive genetic variance. Pooled whole genome resequencing of a subset of resistant (n = 60) and susceptible animals (n = 60) was performed to characterize QTL regions, including identification of putative candidate genes and functional annotation of associated polymorphisms. The TRIM25 gene was identified as a promising positional and functional candidate within the QTL region of LG 44, and a putative premature stop mutation in this gene was discovered.

Contribution of epigenetic variation to adaptation in Arabidopsis

In plants, transgenerational inheritance of some epialleles has been demonstrated but it remains controversial whether epigenetic variation is subject to selection and contributes to adaptation. Simulating selection in a rapidly changing environment, we compare phenotypic traits and epigenetic variation between Arabidopsis thaliana populations grown for five generations under selection and their genetically nearly identical ancestors. Selected populations of two distinct genotypes show significant differences in flowering time and plant architecture, which are maintained for at least 2–3 generations in the absence of selection. While we cannot detect consistent genetic changes, we observe a reduction of epigenetic diversity and changes in the methylation state of about 50,000 cytosines, some of which are associated with phenotypic changes. Thus, we propose that epigenetic variation is subject to selection and can contribute to rapid adaptive responses, although the extent to which epigenetics plays a role in adaptation is still unclear.

Whether plant epigenetic variation is subject to selection and contributes to adaptation is under debate. Here, the authors compare DNA methylation and phenotypes of Arabidopsis lines subject to simulated selection and their nearly isogenic ancestors and provide evidence that epigenetic variation contributes to adaptive responses.

Disentangling the Causes for Faster-X Evolution in Aphids

The faster evolution of X chromosomes has been documented in several species, and results from the increased efficiency of selection on recessive alleles in hemizygous males and/or from increased drift due to the smaller effective population size of X chromosomes. Aphids are excellent models for evaluating the importance of selection in faster-X evolution because their peculiar life cycle and unusual inheritance of sex chromosomes should generally lead to equivalent effective population sizes for X and autosomes. Because we lack a high-density genetic map for the pea aphid, whose complete genome has been sequenced, we first assigned its entire genome to the X or autosomes based on ratios of sequencing depth in males (X0) to females (XX). Then, we computed nonsynonymous to synonymous substitutions ratios (dN/dS) for the pea aphid gene set and found faster evolution of X-linked genes. Our analyses of substitution rates, together with polymorphism and expression data, showed that relaxed selection is likely to be the greatest contributor to faster-X because a large fraction of X-linked genes are expressed at low rates and thus escape selection. Yet, a minor role for positive selection is also suggested by the difference between substitution rates for X and autosomes for male-biased genes (but not for asexual female-biased genes) and by lower Tajima’s D for X-linked compared with autosomal genes with highly male-biased expression patterns. This study highlights the relevance of organisms displaying alternative chromosomal inheritance to the understanding of forces shaping genome evolution.

GPrank: an R package for detecting dynamic elements from genome-wide time series

Background

Genome-wide high-throughput sequencing (HTS) time series experiments are a powerful tool for monitoring various genomic elements over time. They can be used to monitor, for example, gene or transcript expression with RNA sequencing (RNA-seq), DNA methylation levels with bisulfite sequencing (BS-seq), or abundances of genetic variants in populations with pooled sequencing (Pool-seq). However, because of high experimental costs, the time series data sets often consist of a very limited number of time points with very few or no biological replicates, posing challenges in the data analysis.

Results

Here we present the GPrank R package for modelling genome-wide time series by incorporating variance information obtained during pre-processing of the HTS data using probabilistic quantification methods or from a beta-binomial model using sequencing depth. GPrank is well-suited for analysing both short and irregularly sampled time series. It is based on modelling each time series by two Gaussian process (GP) models, namely, time-dependent and time-independent GP models, and comparing the evidence provided by data under two models by computing their Bayes factor (BF). Genomic elements are then ranked by their BFs, and temporally most dynamic elements can be identified.

Conclusions

Incorporating the variance information helps GPrank avoid false positives without compromising computational efficiency. Fitted models can be easily further explored in a browser. Detection and visualisation of temporally most active dynamic elements in the genome can provide a good starting point for further downstream analyses for increasing our understanding of the studied processes.

A non-coding region near Follistatin controls head colour polymorphism in the Gouldian finch

Discrete colour morphs coexisting within a single population are common in nature. In a broad range of organisms, sympatric colour morphs often display major differences in other traits, including morphology, physiology or behaviour. Despite the repeated occurrence of this phenomenon, our understanding of the genetics that underlie multi-trait differences and the factors that promote the long-term maintenance of phenotypic variability within a freely interbreeding population are incomplete. Here, we investigated the genetic basis of red and black head colour in the Gouldian finch (Erythrura gouldiae), a classic polymorphic system in which naturally occurring colour morphs also display differences in aggressivity and reproductive success. We show that the candidate locus is a small (approx. 70 kb) non-coding region mapping to the Z chromosome near the Follistatin (FST) gene. Unlike recent findings in other systems where phenotypic morphs are explained by large inversions containing hundreds of genes (so-called supergenes), we did not identify any structural rearrangements between the two haplotypes using linked-read sequencing technology. Nucleotide divergence between the red and black alleles was high when compared to the remainder of the Z chromosome, consistent with their maintenance as balanced polymorphisms over several million years. Our results illustrate how pleiotropic phenotypes can arise from simple genetic variation, probably regulatory in nature.

Selection at a genomic region of major effect is responsible for evolution of complex life histories in anadromous steelhead

BACKGROUND

Disparity in the timing of biological events occurs across a variety of systems, yet the understanding of genetic basis underlying diverse phenologies remains limited. Variation in maturation timing occurs in steelhead trout, which has been associated with greb1L, an oestrogen target gene. Previous techniques that identified this gene only accounted for about 0.5-2.0% of the genome and solely investigated coastal populations, leaving uncertainty on the genetic basis of this trait and its prevalence across a larger geographic scale.

RESULTS

We used a three-tiered approach to interrogate the genomic basis of complex phenology in anadromous steelhead. First, fine scale mapping with 5.3 million SNPs from resequencing data covering 68% of the genome confirmed a 309-kb region consisting of four genes on chromosome 28, including greb1L, to be the genomic region of major effect for maturation timing. Second, broad-scale characterization of candidate greb1L genotypes across 59 populations revealed unexpected patterns in maturation phenology for inland fish migrating long distances relative to those in coastal streams. Finally, genotypes from 890 PIT-tag tracked steelhead determined associations with early versus late arrival to spawning grounds that were previously unknown.

CONCLUSIONS

This study clarifies the genetic bases for disparity in phenology observed in steelhead, determining an unanticipated trait association with premature versus mature arrival to spawning grounds and identifying multiple candidate genes potentially contributing to this variation from a single genomic region of major effect. This illustrates how dense genome mapping and detailed phenotypic characterization can clarify genotype to phenotype associations across geographic ranges of species.

The Cytoplasm Affects the Epigenome in Drosophila melanogaster

Cytoplasmic components and their interactions with the nuclear genome may mediate patterns of phenotypic expression to form a joint inheritance system. However, proximate mechanisms underpinning these interactions remain elusive. To independently assess nuclear genetic and epigenetic cytoplasmic effects, we created a full-factorial design in which representative cytoplasms and nuclear backgrounds from each of two geographically disjunct populations of Drosophila melanogaster were matched together in all four possible combinations. To capture slowly-accumulating epimutations in addition to immediately occurring ones, these constructed populations were examined one year later. We found the K4 methylation of histone H3, H3K4me3, an epigenetic marker associated with transcription start-sites had diverged across different cytoplasms. The loci concerned mainly related to metabolism, mitochondrial function, and reproduction. We found little overlap (<8%) in sites that varied genetically and epigenetically, suggesting that epigenetic changes have diverged independently from any cis-regulatory sequence changes. These results are the first to show cytoplasm-specific effects on patterns of nuclear histone methylation. Our results highlight that experimental nuclear-cytoplasm mismatch may be used to provide a platform to identify epigenetic candidate loci to study the molecular mechanisms of cyto-nuclear interactions.

The RNA Complement of Outer Membrane Vesicles From Salmonella enterica Serovar Typhimurium Under Distinct Culture Conditions

Bacterial outer membrane vesicles (OMVs), as well as OMV-associated small RNAs, have been demonstrated to play a role in host-pathogen interactions. The presence of larger RNA transcripts in OMVs has been less studied and their potential role in host-pathogen interactions remains largely unknown. Here we analyze RNA from OMVs secreted by Salmonella enterica serovar Typhimurium (S. Typhimurium) cultured under different conditions, which mimic host-pathogen interactions. S. Typhimurium was grown to exponential and stationary growth phases in minimal growth control medium (phosphate-carbon-nitrogen, PCN), as well as in acidic and phosphate-depleted PCN, comparable to the macrophage environment and inducing therefore the expression of Salmonella pathogenicity island 2 (SPI-2) genes. Moreover, Salmonella pathogenicity island 1 (SPI-1), which is required for virulence during the intestinal phase of infection, was induced by culturing S. Typhimurium to the stationary phase in Lysogeny Broth (LB). For each condition, we identified OMV-associated RNAs that are enriched in the extracellular environment relative to the intracellular space. All RNA classes could be observed, but a vast majority of rRNA was exported in all conditions in variable proportions with a notable decrease in LB SPI-1 inducing media. Several mRNAs and ncRNAs were specifically enriched in/on OMVs dependent on the growth conditions. Important to note is that some RNAs showed identical read coverage profiles intracellularly and extracellularly, whereas distinct coverage patterns were observed for other transcripts, suggesting a specific processing or degradation. Moreover, PCR experiments confirmed that distinct RNAs were present in or on OMVs as full-length transcripts (IsrB-1/2; IsrA; ffs; SsrS; CsrC; pSLT035; 10Sa; rnpB; STM0277; sseB; STM0972; STM2606), whereas others seemed to be rather present in a processed or degraded form. Finally, we show by a digestion protection assay that OMVs are able to prevent enzymatic degradation of given full-length transcripts (SsrS, CsrC, 10Sa, and rnpB). In summary, we show that OMV-associated RNA is clearly different in distinct culture conditions and that at least a fraction of the extracellular RNA is associated as a full-length transcripts with OMVs, indicating that some RNAs are protected by OMVs and thereby leaving open the possibility that those might be functionally active.

A simple genetic basis of adaptation to a novel thermal environment results in complex metabolic rewiring in Drosophila

BACKGROUND

Population genetic theory predicts that rapid adaptation is largely driven by complex traits encoded by many loci of small effect. Because large-effect loci are quickly fixed in natural populations, they should not contribute much to rapid adaptation.

RESULTS

To investigate the genetic architecture of thermal adaptation - a highly complex trait - we performed experimental evolution on a natural Drosophila simulans population. Transcriptome and respiration measurements reveal extensive metabolic rewiring after only approximately 60 generations in a hot environment. Analysis of genome-wide polymorphisms identifies two interacting selection targets, Sestrin and SNF4Aγ, pointing to AMPK, a central metabolic switch, as a key factor for thermal adaptation.

CONCLUSIONS

Our results demonstrate that large-effect loci segregating at intermediate allele frequencies can allow natural populations to rapidly respond to selection. Because SNF4Aγ also exhibits clinal variation in various Drosophila species, we suggest that this large-effect polymorphism is maintained by temporal and spatial temperature variation in natural environments.

Genetic variation underlying local adaptation of diapause induction along a cline in a butterfly

Diapause is a life history strategy allowing individuals to arrest development until favourable conditions return, and it is commonly induced by shortened day length that is latitude specific for local populations. Although understanding the evolutionary dynamics of a threshold trait like diapause induction provides insights into the adaptive process and adaptive potential of populations, the genetic mechanism of variation in photoperiodic induction of diapause is not well understood. Here, we investigate genetic variation underlying latitudinal variation in diapause induction and the selection dynamics acting upon it. Using a genomewide scan for divergent regions between two populations of the butterfly Pararge aegeria that differ strongly in their induction thresholds, we identified and investigated the patterns of variation in those regions. We then tested the association of these regions with diapause induction using between-population crosses, finding significant SNP associations in four genes present in two chromosomal regions, one with the gene period, and the other with the genes kinesin, carnitine O-acetyltransferase and timeless. Patterns of allele frequencies in these two regions in population samples along a latitudinal cline suggest strong selection against heterozygotes at two genes within these loci (period, timeless). Evidence for additional loci modifying the diapause decision was found in patterns of allelic change in relation to induction thresholds over the cline, as well as in backcross analyses. Taken together, population-specific adaptations of diapause induction appear to be due to a combination of alleles of larger and smaller effect size, consistent with an exponential distribution of effect sizes involved in local adaption.

Oak genome reveals facets of long lifespan

Oaks are an important part of our natural and cultural heritage. Not only are they ubiquitous in our most common landscapes¹ but they have also supplied human societies with invaluable services, including food and shelter, since prehistoric times². With 450 species spread throughout Asia, Europe and America³, oaks constitute a critical global renewable resource. The longevity of oaks (several hundred years) probably underlies their emblematic cultural and historical importance. Such long-lived sessile organisms must persist in the face of a wide range of abiotic and biotic threats over their lifespans. We investigated the genomic features associated with such a long lifespan by sequencing, assembling and annotating the oak genome. We then used the growing number of whole-genome sequences for plants (including tree and herbaceous species) to investigate the parallel evolution of genomic characteristics potentially underpinning tree longevity. A further consequence of the long lifespan of trees is their accumulation of somatic mutations during mitotic divisions of stem cells present in the shoot apical meristems. Empirical⁴ and modelling⁵ approaches have shown that intra-organismal genetic heterogeneity can be selected for⁶ and provides direct fitness benefits in the arms race with short-lived pests and pathogens through a patchwork of intra-organismal phenotypes⁷. However, there is no clear proof that large-statured trees consist of a genetic mosaic of clonally distinct cell lineages within and between branches. Through this case study of oak, we demonstrate the accumulation and transmission of somatic mutations and the expansion of disease-resistance gene families in trees.

Pleiotropic effects of regulatory variation in tan result in correlation of two pigmentation traits in Drosophila melanogaster

Traits with a common genetic basis frequently display correlated phenotypic responses to selection or environmental conditions. In Drosophila melanogaster, pigmentation of the abdomen and a trident-shaped region on the thorax are genetically correlated. Here, we used a pooled replicated genomewide association approach (Pool-GWAS) to identify the genetic basis of variation in thoracic trident pigmentation in two Drosophila melanogaster populations. We confirmed the previously reported large effect of ebony and the association of the cosmopolitan inversion In(3R)Payne. For the first time, we identified tan as another major locus contributing to variation in trident pigmentation. Intriguingly, the regulatory variants of tan that were most strongly associated with female abdominal pigmentation also showed a strong association with trident pigmentation. We validated this common genetic basis in transgenic assays and found qualitatively similar effects on trident and abdominal pigmentation. Further work is required to determine whether this genetic correlation is favoured by natural selection or reflects a neutral by-product of a shared regulatory architecture.

Adaptive introgression underlies polymorphic seasonal camouflage in snowshoe hares

Snowshoe hares (Lepus americanus) maintain seasonal camouflage by molting to a white winter coat, but some hares remain brown during the winter in regions with low snow cover. We show that cis-regulatory variation controlling seasonal expression of the Agouti gene underlies this adaptive winter camouflage polymorphism. Genetic variation at Agouti clustered by winter coat color across multiple hare and jackrabbit species, revealing a history of recurrent interspecific gene flow. Brown winter coats in snowshoe hares likely originated from an introgressed black-tailed jackrabbit allele that has swept to high frequency in mild winter environments. These discoveries show that introgression of genetic variants that underlie key ecological traits can seed past and ongoing adaptation to rapidly changing environments.

Extensive exchange of transposable elements in the Drosophila pseudoobscura group

BACKGROUND

As species diverge, so does their transposable element (TE) content. Within a genome, TE families may eventually become dormant due to host-silencing mechanisms, natural selection and the accumulation of inactive copies. The transmission of active copies from a TE families, both vertically and horizontally between species, can allow TEs to escape inactivation if it occurs often enough, as it may allow TEs to temporarily escape silencing in a new host. Thus, the contribution of horizontal exchange to TE persistence has been of increasing interest.

RESULTS

Here, we annotated TEs in five species with sequenced genomes from the D. pseudoobscura species group, and curated a set of TE families found in these species. We found that, compared to host genes, many TE families showed lower neutral divergence between species, consistent with recent transmission of TEs between species. Despite these transfers, there are differences in the TE content between species in the group.

CONCLUSIONS

The TE content is highly dynamic in the D. pseudoobscura species group, frequently transferring between species, keeping TEs active. This result highlights how frequently transposable elements are transmitted between sympatric species and, despite these transfers, how rapidly species TE content can diverge.

Complex signatures of genomic variation of two non-model marine species in a homogeneous environment

BACKGROUND

Genomic tools are increasingly being used on non-model organisms to provide insights into population structure and variability, including signals of selection. However, most studies are carried out in regions with distinct environmental gradients or across large geographical areas, in which local adaptation is expected to occur. Therefore, the focus of this study is to characterize genomic variation and selective signals over short geographic areas within a largely homogeneous region. To assess adaptive signals between microhabitats within the rocky shore, we compared genomic variation between the Cape urchin (Parechinus angulosus), which is a low to mid-shore species, and the Granular limpet (Scutellastra granularis), a high shore specialist.

RESULTS

Using pooled restriction site associated DNA (RAD) sequencing, we described patterns of genomic variation and identified outlier loci in both species. We found relatively low numbers of outlier SNPs within each species, and identified outlier genes associated with different selective pressures than those previously identified in studies conducted over larger environmental gradients. The number of population-specific outlier loci differed between species, likely owing to differential selective pressures within the intertidal environment. Interestingly, the outlier loci were highly differentiated within the two northernmost populations for both species, suggesting that unique evolutionary forces are acting on marine invertebrates within this region.

CONCLUSIONS

Our study provides a background for comparative genomic studies focused on non-model species, as well as a baseline for the adaptive potential of marine invertebrates along the South African west coast. We also discuss the caveats associated with Pool-seq and potential biases of sequencing coverage on downstream genomic metrics. The findings provide evidence of species-specific selective pressures within a homogeneous environment, and suggest that selective forces acting on small scales are just as crucial to acknowledge as those acting on larger scales. As a whole, our findings imply that future population genomic studies should expand from focusing on model organisms and/or studying heterogeneous regions to better understand the evolutionary processes shaping current and future biodiversity patterns, particularly when used in a comparative phylogeographic context.

Genome-wide regulatory deterioration impedes adaptive responses to stress in inbred populations of Drosophila melanogaster

Inbreeding depression is often intensified under environmental stress (i.e., inbreeding-stress interaction). Although the fitness consequences of this phenomenon are well-described, underlying mechanisms such as an increased expression of deleterious alleles under stress, or a lower capacity for adaptive responses to stress with inbreeding, have rarely been investigated. We investigated a fitness component (egg-to-adult viability) and gene-expression patterns using RNA-seq analyses in noninbred control lines and in inbred lines of Drosophila melanogaster exposed to benign temperature or heat stress. We find little support for an increase in the cumulative expression of deleterious alleles under stress. Instead, inbred individuals had a reduced ability to induce an adaptive gene regulatory stress response compared to controls. The decrease in egg-to-adult viability due to stress was most pronounced in the lines with the largest deviation in the adaptive stress response (R2 = 0.48). Thus, we find strong evidence for a lower capacity of inbred individuals to respond by gene regulation to stress and that this is the main driver of inbreeding-stress interactions. In comparison, the altered gene expression due to inbreeding at benign temperature showed no correlation with fitness and was pronounced in genomic regions experiencing the highest increase in homozygosity.

Genomewide signatures of selection in Epichloë reveal candidate genes for host specialization

Host specialization is a key process in ecological divergence and speciation of plant-associated fungi. The underlying determinants of host specialization are generally poorly understood, especially in endophytes, which constitute one of the most abundant components of the plant microbiome. We addressed the genetic basis of host specialization in two sympatric subspecies of grass-endophytic fungi from the Epichloë typhina complex: subsp. typhina and clarkii. The life cycle of these fungi entails unrestricted dispersal of gametes and sexual reproduction before infection of a new host, implying that the host imposes a selective barrier on viability of the progeny. We aimed to detect genes under divergent selection between subspecies, experiencing restricted gene flow due to adaptation to different hosts. Using pooled whole-genome sequencing data, we combined F_ST and D_XY population statistics in genome scans and detected 57 outlier genes showing strong differentiation between the two subspecies. Genomewide analyses of nucleotide diversity (π), Tajima's D and dN/dS ratios indicated that these genes have evolved under positive selection. Genes encoding secreted proteins were enriched among the genes showing evidence of positive selection, suggesting that molecular plant-fungus interactions are strong drivers of endophyte divergence. We focused on five genes encoding secreted proteins, which were further sequenced in 28 additional isolates collected across Europe to assess genetic variation in a larger sample size. Signature of positive selection in these isolates and putative identification of pathogenic function supports our findings that these genes represent strong candidates for host specialization determinants in Epichloë endophytes. Our results highlight the role of secreted proteins as key determinants of host specialization.

G-protein-coupled receptor genes of Dirofilaria immitis

The diversity and uniqueness of nematode heterotrimeric G-protein-coupled receptors (GPCRs) provides impetus for identifying ligands that can be used as therapeutics for treating diseases caused by parasitic nematode infections. In human medicine, GPCRs have represented the largest group of 'drugable' targets exploited in the market today. In the filarial nematode Dirofilaria immitis, which causes heartworm disease, the macrocyclic lactones (ML) have been used as the sole preventatives for more than 25 years and now there is confirmed ML resistance in this parasite. A novel anthelmintic emodepside, with antifilarial activity, can act on a GPCR. In view of the ML resistance, there is an urgent need to identify new drug targets and GPCRs of D. immitis may be promising receptors. Knowledge of polymorphism within the GPCR superfamily is of interest. A total of 127 GPCR genes have been identified, so far, in the genome of D. immitis. Whole genome sequencing data from four ML susceptible and four ML loss of efficacy populations was used to identify 393 polymorphic loci in 35 D. immitis GPCR genes. Out of 57 SNPs in exonic regions, 36 of them caused a change in an amino acid, out of which 2 changed the predicted secondary structure of the protein. Knowledge about GPCR genes and their polymorphism is valuable information for drug design processes. Further studies need to be carried out to more fully understand the implications of each of the SNPs identified by this study.

Genomic selection signatures in sheep from the Western Pyrenees

BACKGROUND

The current large spectrum of sheep phenotypic diversity results from the combined product of sheep selection for different production traits such as wool, milk and meat, and its natural adaptation to new environments. In this study, we scanned the genome of 25 Sasi Ardi and 75 Latxa sheep from the Western Pyrenees for three types of regions under selection: (1) regions underlying local adaptation of Sasi Ardi semi-feral sheep, (2) regions related to a long traditional dairy selection pressure in Latxa sheep, and (3) regions experiencing the specific effect of the modern genetic improvement program established for the Latxa breed during the last three decades.

RESULTS

Thirty-two selected candidate regions including 147 annotated genes were detected by using three statistical parameters: pooled heterozygosity H, Tajima's D, and Wright's fixation index F_st. For Sasi Ardi sheep, chromosomes Ovis aries (OAR)4, 6, and 22 showed the strongest signals and harbored several candidate genes related to energy metabolism and morphology (BBS9, ELOVL3 and LDB1), immunity (NFKB2), and reproduction (H2AFZ). The major genomic difference between Sasi Ardi and Latxa sheep was on OAR6, which is known to affect milk production, with highly selected regions around the ABCG2, SPP1, LAP3, NCAPG, LCORL, and MEPE genes in Latxa sheep. The effect of the modern genetic improvement program on Latxa sheep was also evident on OAR15, on which several olfactory genes are located. We also detected several genes involved in reproduction such as ESR1 and ZNF366 that were affected by this selection program.

CONCLUSIONS

Natural and artificial selection have shaped the genome of both Sasi Ardi and Latxa sheep. Our results suggest that Sasi Ardi traits related to energy metabolism, morphological, reproductive, and immunological features have been under positive selection to adapt this semi-feral sheep to its particular environment. The highly selected Latxa sheep for dairy production showed clear signatures of selection in genomic regions related to milk production. Furthermore, our data indicate that the selection criteria applied in the modern genetic improvement program affect immunity and reproduction traits.

Navigating the Interface Between Landscape Genetics and Landscape Genomics

As next-generation sequencing data become increasingly available for non-model organisms, a shift has occurred in the focus of studies of the geographic distribution of genetic variation. Whereas landscape genetics studies primarily focus on testing the effects of landscape variables on gene flow and genetic population structure, landscape genomics studies focus on detecting candidate genes under selection that indicate possible local adaptation. Navigating the transition between landscape genomics and landscape genetics can be challenging. The number of molecular markers analyzed has shifted from what used to be a few dozen loci to thousands of loci and even full genomes. Although genome scale data can be separated into sets of neutral loci for analyses of gene flow and population structure and putative loci under selection for inference of local adaptation, there are inherent differences in the questions that are addressed in the two study frameworks. We discuss these differences and their implications for study design, marker choice and downstream analysis methods. Similar to the rapid proliferation of analysis methods in the early development of landscape genetics, new analytical methods for detection of selection in landscape genomics studies are burgeoning. We focus on genome scan methods for detection of selection, and in particular, outlier differentiation methods and genetic-environment association tests because they are the most widely used. Use of genome scan methods requires an understanding of the potential mismatches between the biology of a species and assumptions inherent in analytical methods used, which can lead to high false positive rates of detected loci under selection. Key to choosing appropriate genome scan methods is an understanding of the underlying demographic structure of study populations, and such data can be obtained using neutral loci from the generated genome-wide data or prior knowledge of a species' phylogeographic history. To this end, we summarize recent simulation studies that test the power and accuracy of genome scan methods under a variety of demographic scenarios and sampling designs. We conclude with a discussion of additional considerations for future method development, and a summary of methods that show promise for landscape genomics studies but are not yet widely used.

Strong phenotypic plasticity limits potential for evolutionary responses to climate change

Phenotypic plasticity, the expression of multiple phenotypes from one genome, is a widespread adaptation to short-term environmental fluctuations, but whether it facilitates evolutionary adaptation to climate change remains contentious. Here, we investigate seasonal plasticity and adaptive potential in an Afrotropical butterfly expressing distinct phenotypes in dry and wet seasons. We assess the transcriptional architecture of plasticity in a full-factorial analysis of heritable and environmental effects across 72 individuals, and reveal pervasive gene expression differences between the seasonal phenotypes. Strikingly, intra-population genetic variation for plasticity is largely absent, consistent with specialisation to a particular environmental cue reliably predicting seasonal transitions. Under climate change, deteriorating accuracy of predictive cues will likely aggravate maladaptive phenotype-environment mismatches and increase selective pressures on reaction norms. However, the observed paucity of genetic variation for plasticity limits evolutionary responses, potentially weakening prospects for population persistence. Thus, seasonally plastic species may be especially vulnerable to climate change.

Diversity and evolution of the repetitive genomic content in Cannabis sativa

Background

The repetitive content of the genome, once considered to be “junk DNA”, is in fact an essential component of genomic architecture and evolution. In this study, we used the genomes of three varieties of Cannabis sativa, three varieties of Humulus lupulus and one genotype of Morus notabilis to explore their repetitive content using a graph-based clustering method, designed to explore and compare repeat content in genomes that have not been fully assembled.

Results

The repetitive content in the C. sativa genome is mainly composed of the retrotransposons LTR/Copia and LTR/Gypsy (14% and 14.8%, respectively), ribosomal DNA (2%), and low-complexity sequences (29%). We observed a recent copy number expansion in some transposable element families. Simple repeats and low complexity regions of the genome show higher intra and inter species variation.

Conclusions

As with other sequenced genomes, the repetitive content of C. sativa’s genome exhibits a wide range of evolutionary patterns. Some repeat types have patterns of diversity consistent with expansions followed by losses in copy number, while others may have expanded more slowly and reached a steady state. Still, other repetitive sequences, particularly ribosomal DNA (rDNA), show signs of concerted evolution playing a major role in homogenizing sequence variation.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-4494-3) contains supplementary material, which is available to authorized users.

Evolution of cross-resistance to medical triazoles in Aspergillus fumigatus through selection pressure of environmental fungicides

Resistance to medical triazoles in Aspergillus fumigatus is an emerging problem for patients at risk of aspergillus diseases. There are currently two presumed routes for medical triazole-resistance selection: (i) through selection pressure of medical triazoles when treating patients and (ii) through selection pressure from non-medical sterol-biosynthesis-inhibiting (SI) triazole fungicides which are used in the environment. Previous studies have suggested that SI fungicides can induce cross-resistance to medical triazoles. Therefore, to assess the potential of selection of resistance to medical triazoles in the environment, we assessed cross-resistance to three medical triazoles in lineages of A. fumigatus from previous work where we applied an experimental evolution approach with one of five different SI fungicides to select for resistance. In our evolved lines we found widespread cross-resistance indicating that resistance to medical triazoles rapidly arises through selection pressure of SI fungicides. All evolved lineages showed similar evolutionary dynamics to SI fungicides and medical triazoles, which suggests that the mutations inducing resistance to both SI fungicides and medical triazoles are likely to be the same. Whole-genome sequencing revealed that a variety of mutations were putatively involved in the resistance mechanism, some of which are in known target genes.

Quantifying the Survival of Multiple Salmonella enterica Serovars In Vivo via Massively Parallel Whole-Genome Sequencing To Predict Zoonotic Risk

Salmonella enterica is an animal and zoonotic pathogen of worldwide importance. Salmonella serovars that differ in their host and tissue tropisms exist. Cattle are an important reservoir of human nontyphoidal salmonellosis, and contaminated bovine peripheral lymph nodes enter the food chain via ground beef. The relative abilities of different serovars to survive within the bovine lymphatic system are poorly understood and constrain the development of control strategies. This problem was addressed by developing a massively parallel whole-genome sequencing method to study mixed-serovar infections in vivo. Salmonella serovars differ genetically by naturally occurring single nucleotide polymorphisms (SNPs) in certain genes. It was hypothesized that these SNPs could be used as markers to simultaneously identify serovars in mixed populations and quantify the abundance of each member in a population. The performance of the method was validated in vitro using simulated pools containing up to 11 serovars in various proportions. It was then applied to study serovar survival in vivo in cattle challenged orally with the same 11 serovars. All the serovars successfully colonized the bovine lymphatic system, including the peripheral lymph nodes, and thus pose similar risks of zoonosis. This method enables the fates of multiple genetically unmodified strains to be evaluated simultaneously in a single animal. It could be useful in reducing the number of animals required to study mixed-strain infections and in testing the cross-protective efficacy of vaccines and treatments. It also has the potential to be applied to diverse bacterial species which possess shared but polymorphic alleles.

IMPORTANCE While some Salmonella serovars are more frequently isolated from lymph nodes rather than the feces and environment of cattle, the relative abilities of serovars to survive within the lymphatic system of cattle remain ill defined. A sequencing-based method which used available information from sequenced Salmonella genomes to study the dynamics of mixed-serovar infections in vivo was developed. The main advantages of the method include the simultaneous identification and quantification of multiple strains without any genetic modification and minimal animal use. This approach could be used in vaccination trials or in epidemiological surveys where an understanding of the dynamics of closely related strains of a pathogen in mixed populations could inform the prediction of zoonotic risk and the development of intervention strategies.

Elevated Proportions of Deleterious Genetic Variation in Domestic Animals and Plants

A fraction of genetic variants segregating in any population are deleterious, which negatively impacts individual fitness. The domestication of animals and plants is associated with population bottlenecks and artificial selection, which are predicted to increase the proportion of deleterious variants. However, the extent to which this is a general feature of domestic species is unclear. Here, we examine the effects of domestication on the prevalence of deleterious variation using pooled whole-genome resequencing data from five domestic animal species (dog, pig, rabbit, chicken, and silkworm) and two domestic plant species (rice and soybean) compared with their wild ancestors. We find significantly reduced genetic variation and increased proportion of nonsynonymous amino acid changes in all but one of the domestic species. These differences are observable across a range of allele frequencies, both common and rare. We find proportionally more single nucleotide polymorphisms in highly conserved elements in domestic species and a tendency for domestic species to harbor a higher proportion of changes classified as damaging. Our findings most likely reflect an increased incidence of deleterious variants in domestic species, which is most likely attributable to population bottlenecks that lead to a reduction in the efficacy of selection. An exception to this pattern is displayed by European domestic pigs, which do not show traces of a strong population bottleneck and probably continued to exchange genes with wild boar populations after domestication. The results presented here indicate that an elevated proportion of deleterious variants is a common, but not ubiquitous, feature of domestic species.

Discovery of Culex pipiens associated tunisia virus: a new ssRNA(+) virus representing a new insect associated virus family

In the global context of arboviral emergence, deep sequencing unlocks the discovery of new mosquito-borne viruses. Mosquitoes of the species Culex pipiens, C. torrentium, and C. hortensis were sampled from 22 locations worldwide for transcriptomic analyses. A virus discovery pipeline was used to analyze the dataset of 0.7 billion reads comprising 22 individual transcriptomes. Two closely related 6.8 kb viral genomes were identified in C. pipiens and named as Culex pipiens associated tunisia virus (CpATV) strains Ayed and Jedaida. The CpATV genome contained four ORFs. ORF1 possessed helicase and RNA-dependent RNA polymerase (RdRp) domains related to new viral sequences recently found mainly in dipterans. ORF2 and 4 contained a capsid protein domain showing strong homology with Virgaviridae plant viruses. ORF3 displayed similarities with eukaryotic Rhoptry domain and a merozoite surface protein (MSP7) domain only found in mosquito-transmitted Plasmodium, suggesting possible interactions between CpATV and vertebrate cells. Estimation of a strong purifying selection exerted on each ORFs and the presence of a polymorphism maintained in the coding region of ORF3 suggested that both CpATV sequences are genuine functional viruses. CpATV is part of an entirely new and highly diversified group of viruses recently found in insects, and that bears the genomic hallmarks of a new viral family.

Infection via mosquito bite alters Zika virus tissue tropism and replication kinetics in rhesus macaques

Mouse and nonhuman primate models now serve as useful platforms to study Zika virus (ZIKV) pathogenesis, candidate therapies, and vaccines, but they rely on needle inoculation of virus: the effects of mosquito-borne infection on disease outcome have not been explored in these models. Here we show that infection via mosquito bite delays ZIKV replication to peak viral loads in rhesus macaques. Importantly, in mosquito-infected animals ZIKV tissue distribution was limited to hemolymphatic tissues, female reproductive tract tissues, kidney, and liver, potentially emulating key features of human ZIKV infections, most of which are characterized by mild or asymptomatic disease. Furthermore, deep sequencing analysis reveals that ZIKV populations in mosquito-infected monkeys show greater sequence heterogeneity and lower overall diversity than in needle-inoculated animals. This newly developed system will be valuable for studying ZIKV disease because it more closely mimics human infection by mosquito bite than needle-based inoculations.

Vector saliva can affect infectivity and pathogenesis of vector-borne viruses, but this hasn’t been studied for Zika virus infection. Here, Dudley et al. show that mosquito-mediated Zika infection of macaques results in altered replication kinetics and greater sequence heterogeneity.

ReadTools: A universal toolkit for handling sequence data from different sequencing platforms

Sequencing whole genomes has become a standard research tool in many disciplines including Molecular Ecology, but the rapid technological advances in combination with several competing platforms have resulted in a confusing diversity of formats. This lack of standard formats causes several problems, such as undocumented preprocessing steps or the loss of information in downstream software tools, which do not account for the specifics of the different available formats. ReadTools is an open-source Java toolkit designed to standardize and preprocess read data from different platforms. It manages FASTQ- and SAM-formatted inputs while dealing with platform-specific peculiarities and provides a standard SAM compliant output. The code and executable are available at https://github.com/magicDGS/ReadTools.

The dynamic evolution of Drosophila innubila Nudivirus

Viruses coevolve with their hosts to overcome host resistance and gain the upper hand in the evolutionary arms race. Drosophila innubila nudivirus (DiNV) is a double stranded DNA virus, closely related to Oryctes rhinoceros nudivirus (OrNV) and Kallithea virus. DiNV is the first DNA virus found to naturally infect Drosophila and therefore has the potential to be developed as a model for DNA virus immune defense and host/virus coevolution within its well-studied host system. Here we sequence and annotate the genome of DiNV and identify signatures of adaptation, revealing clues for genes involved in host-parasite coevolution. The circular genome is 155,555bp and contains 107 coding open reading frames (ORFs) and a wealth of AT-rich simple sequence repeats. While synteny is highly conserved between DiNV and Kallithea virus, it drops off rapidly as sequences become more divergent, consistent with rampant rearrangements across nudiviruses. Overall, we show that evolution of DiNV is likely due to adaptation of a very few genes coupled with high gene turnover.

Dirofilaria immitis JYD-34 isolate: whole genome analysis

BACKGROUND

Macrocyclic lactone (ML) anthelmintics are used for chemoprophylaxis for heartworm infection in dogs and cats. Cases of dogs becoming infected with heartworms, despite apparent compliance to recommended chemoprophylaxis with approved preventives, has led to such cases being considered as suspected lack of efficacy (LOE). Recently, microfilariae collected from a small number of LOE isolates were used as a source of infection of new host dogs and confirmed to have reduced susceptibility to ML in controlled efficacy studies using L3 challenge in dogs. A specific Dirofilaria immitis laboratory isolate named JYD-34 has also been confirmed to have less than 100% susceptibility to ML-based preventives. For preventive claims against heartworm disease, evidence of 100% efficacy is required by FDA-CVM. It was therefore of interest to determine whether JYD-34 has a genetic profile similar to other documented LOE and confirmed reduced susceptibility isolates or has a genetic profile similar to known ML-susceptible isolates.

METHODS

In this study, the 90Mbp whole genome of the JYD-34 strain was sequenced. This genome was compared using bioinformatics tools to pooled whole genomes of four well-characterized susceptible D. immitis populations, one susceptible Missouri laboratory isolate, as well as the pooled whole genomes of four LOE D. immitis populations. Fixation indexes (F_ST), which allow the genetic structure of each population (isolate) to be compared at the level of single nucleotide polymorphisms (SNP) across the genome, have been calculated. Forty-one previously reported SNP, that appeared to differentiate between susceptible and LOE and confirmed reduced susceptibility isolates, were also investigated in the JYD-34 isolate.

RESULTS

The F_ST analysis, and the analysis of the 41 SNP that appeared to differentiate reduced susceptibility from fully susceptible isolates, confirmed that the JYD-34 isolate has a genome similar to previously investigated LOE isolates, and isolates confirmed to have reduced susceptibility, and to be dissimilar to the susceptible isolates.

CONCLUSIONS

These results provide additional evidence for the link between genotype and the reduced susceptibility phenotype observed in such isolates as JYD-34. Further work on other isolates showing reduced susceptibility to ML is required to demonstrate the value of genetic analysis in predicting the response to ML chemoprophylaxis. The authors suggest that genetic analysis may be useful in helping to interpret the results of in vivo efficacy testing of ML heartworm preventives against D. immitis isolates.

Reproductive isolation and environmental adaptation shape the phylogeography of mountain pine beetle (Dendroctonus ponderosae)

Chromosomal rearrangement can be an important mechanism driving population differentiation and incipient speciation. In the mountain pine beetle (MPB, Dendroctonus ponderosae), deletions on the Y chromosome that are polymorphic among populations are associated with reproductive incompatibility. Here, we used RAD sequencing across the entire MPB range in western North America to reveal the extent of the phylogeographic differences between Y haplotypes compared to autosomal and X-linked loci. Clustering and geneflow analyses revealed three distinct Y haplogroups geographically positioned within and on either side of the Great Basin Desert. Despite close geographic proximity between populations on the boundaries of each Y haplogroup, there was extremely low Y haplogroup mixing among populations, and gene flow on the autosomes was reduced across Y haplogroup boundaries. These results are consistent with a previous study suggesting that independent degradation of a recently evolved neo-Y chromosome in previously isolated populations causes male sterility or inviability among Y haplotype lineages. Phylogeographic results supported historic contraction of MPB into three separate Pleistocene glacial refugia followed by postglacial range expansion and secondary contact. Distinct sets of SNPs were statistically associated with environmental data among the most genetically distinct sets of geographic populations. This finding suggests that the process of adaptation to local climatic conditions is influenced by population genetic structure, with evidence for largely independent evolution in the most genetically isolated Y haplogroup.

Variation in the Intensity of Selection on Codon Bias over Time Causes Contrasting Patterns of Base Composition Evolution in Drosophila

Four-fold degenerate coding sites form a major component of the genome, and are often used to make inferences about selection and demography, so that understanding their evolution is important. Despite previous efforts, many questions regarding the causes of base composition changes at these sites in Drosophila remain unanswered. To shed further light on this issue, we obtained a new whole-genome polymorphism data set from D. simulans. We analyzed samples from the putatively ancestral range of D. simulans, as well as an existing polymorphism data set from an African population of D. melanogaster. By using D. yakuba as an outgroup, we found clear evidence for selection on 4-fold sites along both lineages over a substantial period, with the intensity of selection increasing with GC content. Based on an explicit model of base composition evolution, we suggest that the observed AT-biased substitution pattern in both lineages is probably due to an ancestral reduction in selection intensity, and is unlikely to be the result of an increase in mutational bias towards AT alone. By using two polymorphism-based methods for estimating selection coefficients over different timescales, we show that the selection intensity on codon usage has been rather stable in D. simulans in the recent past, but the long-term estimates in D. melanogaster are much higher than the short-term ones, indicating a continuing decline in selection intensity, to such an extent that the short-term estimates suggest that selection is only active in the most GC-rich parts of the genome. Finally, we provide evidence for complex evolutionary patterns in the putatively neutral short introns, which cannot be explained by the standard GC-biased gene conversion model. These results reveal a dynamic picture of base composition evolution.

Differentiated transcriptional signatures in the maize landraces of Chiapas, Mexico

Background

Landrace farmers are the keepers of crops locally adapted to the environments where they are cultivated. Patterns of diversity across the genome can provide signals of past evolution in the face of abiotic and biotic change. Understanding this rich genetic resource is imperative especially since diversity can provide agricultural security as climate continues to shift.

Results

Here we employ RNA sequencing (RNA-seq) to understand the role that conditions that vary across a landscape may have played in shaping genetic diversity in the maize landraces of Chiapas, Mexico. We collected landraces from three distinct elevational zones and planted them in a midland common garden. Early season leaf tissue was collected for RNA-seq and we performed weighted gene co-expression network analysis (WGCNA). We then used association analysis between landrace co-expression module expression values and environmental parameters of landrace origin to elucidate genes and gene networks potentially shaped by environmental factors along our study gradient. Elevation of landrace origin affected the transcriptome profiles. Two co-expression modules were highly correlated with temperature parameters of landrace origin and queries into their ‘hub’ genes suggested that temperature may have led to differentiation among landraces in hormone biosynthesis/signaling and abiotic and biotic stress responses. We identified several ‘hub’ transcription factors and kinases as candidates for the regulation of these responses.

Conclusions

These findings indicate that natural selection may influence the transcriptomes of crop landraces along an elevational gradient in a major diversity center, and provide a foundation for exploring the genetic basis of local adaptation. While we cannot rule out the role of neutral evolutionary forces in the patterns we have identified, combining whole transcriptome sequencing technologies, established bioinformatics techniques, and common garden experimentation can powerfully elucidate structure of adaptive diversity across a varied landscape. Ultimately, gaining such understanding can facilitate the conservation and strategic utilization of crop genetic diversity in a time of climate change.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-017-4005-y) contains supplementary material, which is available to authorized users.

Unexpected high genetic diversity in small populations suggests maintenance by associative overdominance

The effective population size (N_e ) is a central factor in determining maintenance of genetic variation. The neutral theory predicts that loss of variation depends on N_e , with less genetic drift in larger populations. We monitored genetic drift in 42 Drosophila melanogaster populations of different adult census population sizes (10, 50 or 500) using pooled RAD sequencing. In small populations, variation was lost at a substantially lower rate than expected. This observation was consistent across two ecological relevant thermal regimes, one stable and one with a stressful increase in temperature across generations. Estimated ratios between N_e and adult census size were consistently higher in small than in larger populations. The finding provides evidence for a slower than expected loss of genetic diversity and consequently a higher than expected long-term evolutionary potential in small fragmented populations. More genetic diversity was retained in areas of low recombination, suggesting that associative overdominance, driven by disfavoured homozygosity of recessive deleterious alleles, is responsible for the maintenance of genetic diversity in smaller populations. Consistent with this hypothesis, the X-chromosome, which is largely free of recessive deleterious alleles due to hemizygosity in males, fits neutral expectations even in small populations. Our experiments provide experimental answers to a range of unexpected patterns in natural populations, ranging from variable diversity on X-chromosomes and autosomes to surprisingly high levels of nucleotide diversity in small populations.

Whole-genome sequencing approaches for conservation biology: Advantages, limitations and practical recommendations

Whole-genome resequencing (WGR) is a powerful method for addressing fundamental evolutionary biology questions that have not been fully resolved using traditional methods. WGR includes four approaches: the sequencing of individuals to a high depth of coverage with either unresolved or resolved haplotypes, the sequencing of population genomes to a high depth by mixing equimolar amounts of unlabelled-individual DNA (Pool-seq) and the sequencing of multiple individuals from a population to a low depth (lcWGR). These techniques require the availability of a reference genome. This, along with the still high cost of shotgun sequencing and the large demand for computing resources and storage, has limited their implementation in nonmodel species with scarce genomic resources and in fields such as conservation biology. Our goal here is to describe the various WGR methods, their pros and cons and potential applications in conservation biology. WGR offers an unprecedented marker density and surveys a wide diversity of genetic variations not limited to single nucleotide polymorphisms (e.g., structural variants and mutations in regulatory elements), increasing their power for the detection of signatures of selection and local adaptation as well as for the identification of the genetic basis of phenotypic traits and diseases. Currently, though, no single WGR approach fulfils all requirements of conservation genetics, and each method has its own limitations and sources of potential bias. We discuss proposed ways to minimize such biases. We envision a not distant future where the analysis of whole genomes becomes a routine task in many nonmodel species and fields including conservation biology.

Different Genomic Changes Underlie Adaptive Evolution in Populations of Contrasting History

Experimental evolution is a powerful tool to understand the adaptive potential of populations under environmental change. Here, we study the importance of the historical genetic background in the outcome of evolution at the genome-wide level. Using the natural clinal variation of Drosophila subobscura, we sampled populations from two contrasting latitudes (Adraga, Portugal and Groningen, Netherlands) and introduced them in a new common environment in the laboratory. We characterized the genome-wide temporal changes underlying the evolutionary dynamics of these populations, which had previously shown fast convergence at the phenotypic level, but not at chromosomal inversion frequencies. We found that initially differentiated populations did not converge either at genome-wide level or at candidate SNPs with signs of selection. In contrast, populations from Portugal showed convergence to the control population that derived from the same geographical origin and had been long-established in the laboratory. Candidate SNPs showed a variety of different allele frequency change patterns across generations, indicative of an underlying polygenic basis. We did not detect strong linkage around candidate SNPs, but rather a small but long-ranging effect. In conclusion, we found that history played a major role in genomic variation and evolution, with initially differentiated populations reaching the same adaptive outcome through different genetic routes.

Copy number increases of transposable elements and protein-coding genes in an invasive fish of hybrid origin

Evolutionary dynamics of structural genetic variation in lineages of hybrid origin is not well explored, although structural mutations may increase in controlled hybrid crosses. We therefore tested whether structural variants accumulate in a fish of recent hybrid origin, invasive Cottus, relative to both parental species Cottus rhenanus and Cottus perifretum. Copy-number variation in exons of 10,979 genes was assessed using comparative genome hybridization arrays. Twelve genes showed significantly higher copy numbers in invasive Cottus compared to both parents. This coincided with increased expression for three genes related to vision, detoxification and muscle development, suggesting possible gene dosage effects. Copy number increases of putative transposons were assessed by comparative mapping of genomic DNA reads against a de novo assembly of 1,005 repetitive elements. In contrast to exons, copy number increases of repetitive elements were common (20.7%) in invasive Cottus, whereas decrease was very rare (0.01%). Among the increased repetitive elements, 53.8% occurred at higher numbers in C. perifretum compared to C. rhenanus, while only 1.4% were more abundant in C. rhenanus. This implies a biased mutational process that amplifies genetic material from one ancestor. To assess the frequency of de novo mutations through hybridization, we screened 64 laboratory-bred F2 offspring between the parental species for copy-number changes at five candidate loci. We found no evidence for new structural variants, indicating that they are too rare to be detected given our sampling scheme. Instead, they must have accumulated over more generations than we observed in a controlled cross.

A novel sex-linked mutant affecting tail formation in Hongshan chicken

The Hongshan chicken is a Chinese indigenous breed that has two distinctly different tail types. Some chickens have stunted tails as compared to the normal phenotype, and they are termed rumpless. Rumplessness in other chicken breeds was caused by a reduction in the number of coccygeal vertebrae. However, X-ray examination showed that rumpless Hongshan chickens possess the normal number of coccygeal vertebrae. Our analyses of the main tail feathers and tissue sections led us to speculate that their stunted tail appearance may be the result of abnormal feather development. To investigate the genetic mechanism underlying rumplessness in Hongshan chickens, we analyzed the results of various crosses. The results indicated that rumplessness is a Z-linked dominant character. In addition, we chose some normal and rumpless individuals for pool-sequencing. Nucleotide diversity and Fst were calculated, and a selective sweep was detected on the Z chromosome. These analyses allowed us to reduce the search area to 71.8–72 Mb on the Z chromosome (galGal5.0). A pseudogene LOC431648 located in this region appeared a strong candidate involving in Wnt/β-catenin signaling pathway to regulate feather development in chickens.

Mitogenome Sequencing in the Genus Camelus Reveals Evidence for Purifying Selection and Long-term Divergence between Wild and Domestic Bactrian Camels

The genus Camelus is an interesting model to study adaptive evolution in the mitochondrial genome, as the three extant Old World camel species inhabit hot and low-altitude as well as cold and high-altitude deserts. We sequenced 24 camel mitogenomes and combined them with three previously published sequences to study the role of natural selection under different environmental pressure, and to advance our understanding of the evolutionary history of the genus Camelus. We confirmed the heterogeneity of divergence across different components of the electron transport system. Lineage-specific analysis of mitochondrial protein evolution revealed a significant effect of purifying selection in the concatenated protein-coding genes in domestic Bactrian camels. The estimated dN/dS < 1 in the concatenated protein-coding genes suggested purifying selection as driving force for shaping mitogenome diversity in camels. Additional analyses of the functional divergence in amino acid changes between species-specific lineages indicated fixed substitutions in various genes, with radical effects on the physicochemical properties of the protein products. The evolutionary time estimates revealed a divergence between domestic and wild Bactrian camels around 1.1 [0.58-1.8] million years ago (mya). This has major implications for the conservation and management of the critically endangered wild species, Camelus ferus.

Genome-wide analysis of ivermectin response by Onchocerca volvulus reveals that genetic drift and soft selective sweeps contribute to loss of drug sensitivity

Background

Treatment of onchocerciasis using mass ivermectin administration has reduced morbidity and transmission throughout Africa and Central/South America. Mass drug administration is likely to exert selection pressure on parasites, and phenotypic and genetic changes in several Onchocerca volvulus populations from Cameroon and Ghana—exposed to more than a decade of regular ivermectin treatment—have raised concern that sub-optimal responses to ivermectin's anti-fecundity effect are becoming more frequent and may spread.

Methodology/Principal findings

Pooled next generation sequencing (Pool-seq) was used to characterise genetic diversity within and between 108 adult female worms differing in ivermectin treatment history and response. Genome-wide analyses revealed genetic variation that significantly differentiated good responder (GR) and sub-optimal responder (SOR) parasites. These variants were not randomly distributed but clustered in ~31 quantitative trait loci (QTLs), with little overlap in putative QTL position and gene content between the two countries. Published candidate ivermectin SOR genes were largely absent in these regions; QTLs differentiating GR and SOR worms were enriched for genes in molecular pathways associated with neurotransmission, development, and stress responses. Finally, single worm genotyping demonstrated that geographic isolation and genetic change over time (in the presence of drug exposure) had a significantly greater role in shaping genetic diversity than the evolution of SOR.

Conclusions/Significance

This study is one of the first genome-wide association analyses in a parasitic nematode, and provides insight into the genomics of ivermectin response and population structure of O. volvulus. We argue that ivermectin response is a polygenically-determined quantitative trait (QT) whereby identical or related molecular pathways but not necessarily individual genes are likely to determine the extent of ivermectin response in different parasite populations. Furthermore, we propose that genetic drift rather than genetic selection of SOR is the underlying driver of population differentiation, which has significant implications for the emergence and potential spread of SOR within and between these parasite populations.

Onchocerciasis is a human parasitic disease endemic across large areas of Sub-Saharan Africa, where more than 99% of the estimated 100 million people globally at-risk live. The microfilarial stage of Onchocerca volvulus causes pathologies ranging from mild itching to visual impairment and ultimately, irreversible blindness. Mass administration of ivermectin kills microfilariae and has an anti-fecundity effect on adult worms by temporarily inhibiting the development in utero and/or release into the skin of new microfilariae, thereby reducing morbidity and transmission. Phenotypic and genetic changes in some parasite populations that have undergone multiple ivermectin treatments in Cameroon and Ghana have raised concern that sub-optimal response to ivermectin's anti-fecundity effect may increase in frequency, reducing the impact of ivermectin-based control measures. We used next generation sequencing of small pools of parasites to define genome-wide genetic differences between phenotypically characterised good and sub-optimal responder parasites from Cameroon and Ghana, and identified multiple regions of the genome that differentiated the response types. These regions were largely different between parasites from these two countries but revealed common molecular pathways that might be involved in determining the extent of response to ivermectin's anti-fecundity effect. These data reveal a more complex than previously described pattern of genetic diversity among O. volvulus populations that differ in their geography and response to ivermectin treatment.