Adaptive substitutions underlying cardiac glycoside insensitivity in insects exhibit epistasis in vivo

  1. Andrew M Taverner
  2. Lu Yang
  3. Zachary J Barile
  4. Becky Lin
  5. Julie Peng
  6. Ana P Pinharanda
  7. Arya S Rao
  8. Bartholomew P Roland
  9. Aaron D Talsma
  10. Daniel Wei
  11. Georg Petschenka
  12. Michael J Palladino  Is a corresponding author
  13. Peter Andolfatto  Is a corresponding author
  1. Princeton University, United States
  2. University of Pittsburgh, United States
  3. University of Pittsburgh School of Medicine, United States
  4. Columbia University, United States
  5. Justus-Liebig-Universität Gießen, Germany
4 figures, 1 table and 5 additional files

Figures

Figure 1 with 4 supplements
Deleterious effects of substitutions at positions 111 and 122 are ameliorated by the permissive substitution A119S.

(A) Patterns of amino acid substitution at sites implicated in CG sensitivity near the H1-H2 transmembrane domain of ATPα1 in representative species from six insect orders (see Supplementary file 1

https://doi.org/10.7554/eLife.48224.002
Figure 1—figure supplement 1
Engineering strategy for generating amino acid substitution lines.

A targeted replacement of exons 2-6b with a mini-white gene was generated by ends-out homologous recombination. The mini-white gene was removed using Cre-lox site-specific recombination to create …

https://doi.org/10.7554/eLife.48224.003
Figure 1—figure supplement 2
Phylogeny showing relationships of the sampled species.

In red are species known specialists on CG-containing plants or inferred to be CG-insensitive.

https://doi.org/10.7554/eLife.48224.005
Figure 1—figure supplement 3
Variant sites in ATPα1 most strongly correlated with substitutions with site 119.

Shown is the distribution of P-values for strength of phylogenetic correlation between site 119 and 270 non-singleton amino acid variants in an alignment including 174 ATPα1 sequences representing …

https://doi.org/10.7554/eLife.48224.006
Figure 1—figure supplement 4
Crystal structure of NKA bound to the cardiac glycoside ouabain (PDB:4HYT).

Estimated co-crystal structure of Drosophila melanogaster NKA bound to the cardiac glycoside (CG) ouabain. A homology model was constructed with Swiss-Model (Waterhouse et al., 2018; https://swissmod…

https://doi.org/10.7554/eLife.48224.007
Figure 2 with 3 supplements
A119S ameliorates detrimental effects of CG-insensitivity substitutions on NKA function.

Plotted are relative activity as a function of increasing concentrations of the cardiac glycoside (CG) ouabain. (A) Inhibition curves are plotted for heterozygous individuals to allow comparison of …

https://doi.org/10.7554/eLife.48224.008
Figure 2—figure supplement 1
Biphasic inhibition curve fitting allows the estimation of relative activities of two enzyme forms in a mixture.

(A) Shown in gray is the monophasic curve fit corresponding to observed data (gray points) for the D. melanogaster line one homozygous wild- type (+/+) construct. In black is a hypothetical biphasic …

https://doi.org/10.7554/eLife.48224.009
Figure 2—figure supplement 2
Allele-specific expression (ASE) in heterozygous lines.

Allele-specific expression (ASE) in heterozygous lines. Estimates and 95% confidence intervals are based on Cochran-Mantel-Haenzel framework analysis – that is plotted is ‘relative risk’ of the …

https://doi.org/10.7554/eLife.48224.010
Figure 2—figure supplement 3
Enzyme inhibition curves for Q111V (V) with and without A119S (S).

Following Figure 2, inhibition curves are plotted for heterozygous individuals to allow comparison of the effects of these substitutions in the presence and absence of A119S. Points represent …

https://doi.org/10.7554/eLife.48224.011
Figure 3 with 2 supplements
A119S ameliorates detrimental effects of CG-insensitivity substitutions on neural function.

Plotted are recovery times for individuals (open circles), and means with approximate 95% confidence bounds (solid circles with whiskers), following mechanical over-stimulation (aka, the ‘bang …

https://doi.org/10.7554/eLife.48224.012
Figure 3—figure supplement 1
Bang sensitivity phenotypes of lines with heterozygous substitutions at sites 104, 111 and 122.

Bang sensitivity phenotypes of lines with heterozygous substitutions at sites 104, 111 and 122. Δ2-6b is a loss of function deletion of ATPα1. Red letters indicate substitutions at position 111; …

https://doi.org/10.7554/eLife.48224.013
Figure 3—figure supplement 2
Bang sensitivity phenotypes of substitutions at sites 111 and 122 on the background of A119S.

Bang sensitivity phenotypes of substitutions at sites 111 and 122 on the background of A119S. Δ2-6b is a loss of function deletion of ATPα1. Red letters indicate substitutions at position 111; green …

https://doi.org/10.7554/eLife.48224.014
Figure 4 with 1 supplement
Adult survival upon 7 day exposure to the CG ouabain.

Plotted is the log relative risk for treatments (5, 10 or 20 mM ouabain) relative to no treatment controls (no ouabain) for (A) heterozygous strains and (B) homozygous strains. Estimates (points) …

https://doi.org/10.7554/eLife.48224.015
Figure 4—figure supplement 1
Adult survival of all homozygous strains with A119S upon 7 day exposure to CGs.

Adult survival of homozygous strains with A119S upon 7 day exposure to CGs. Plotted is the log relative risk of treatment, for 5 mM (blue) or 10 mM (red) ouabain, relative to no treatment controls …

https://doi.org/10.7554/eLife.48224.016

Tables

Key resources table
Reagent type
(species) or
resource
DesignationSource or
reference
IdentifiersAdditional
information
Gene (Drosophila melanogaster)ATPalpha1NAFLYB:FBgn0027548
Strain, strain background (D. melanogaster)ATPalpha1 △2-6b founder line, w[1118];;ATPa△2-6b attP/TM6B,Tb[1]this paperN/Aavailable on request
Strain, strain background (D. melanogaster)15 engineered variants of ATPalpha1this paperN/ASee Supplementary file 2, available on request
Genetic reagent (D. melanogaster)y[1], w[67c23], P{y[+mDint2]=Crey}1b;;D[*]/TM3, Sb[1]Bloomington Drosophila Stock CenterBDSC:851 FLYB:FBst0000851; RRID:BDSC_851Source of Cre protein
Genetic reagent (D. melanogaster)w[*];;ry[506] Dr[1]/TM6B, P{w[+mC]=Dfd-EYFP}3, Sb[1],Tb[1],ca[1]Bloomington Drosophila Stock CenterBDSC:8704 FLYB:FBst0008704; RRID:BDSC_8704Balancer, chr 3
Sequence-based reagentPCR/sequencing primers, Sanger-based sequencingthis paperN/AATPa1 H1-H2 Forward: AAAACTTGGAGCGCGATGGT; ATPa1 H1-H2 Reverse: ATACGGTCGCCGAACTTCAC
Sequence-based reagentPCR/sequencing primers for Illumina-based sequencingthis paperN/AATPa1 H1-H2 Forward 2834: GTCTCGTGGGCTCGGCGCTTTCAGACACATCCCGA; ATPa1 H1-H2 Forward 3185: TCGTCGGCAGCGTCTGAGAAAATGCCCGTCACGA
Sequence-based reagentPCR-added linker primers, Illumina-based sequencingthis paperN/AATPa1 H1-H2 i5: AATGATACGGCGACCACCGAGATCTACACnnnnnnnnTCGTCGGCAGCGTCA; ATPa1 H1-H2 i7: CAAGCAGAAGACGGCATACGAGATnnnnnnnnGTCTCGTGGGCTCGG; where ‘nnnnnnnn’=any sequence index
Chemical compound, drugOuabainSigmaCat#: O3125
Chemical compound, drugPorcine sodium potassium ATPase, cerebral cortexMillipore-SigmaCat#: A7510
Recombinant DNA reagentpGX-attB-ATPa2-6bthis paperN/Aavailable on request
Software, algorithmPAML v4.8PMID: 17483113RRID:SCR_014932http://abacus.gene.ucl.ac.uk/software/paml.html
Software, algorithmBayesTraits v3.01PMID: 16685633RRID:SCR_014487http://www.evolution.rdg.ac.uk/BayesTraitsV3.0.1/BayesTraitsV3.0.1.html
Othermultiple sequence alignment (protein)this paperN/ASupplementary file 3
Othermultiple sequence alignment (nucleotide)this paperN/ASupplementary file 4

Additional files

Supplementary file 1

Summary of amino acid variation at sites implicated in CG-sensitivity of ATPαone in surveyed species.

Columns correspond to 35 sites implicated in cardenolide-sensitivity with the addition of six additional sites of interest (112, 114, 119, 787, 874, 898). Following convention, position is standardized relative to the sheep (Ovis aries) sequence NM_001009360 minus five amino acid residues from 5'end. The reference sequence refers to the consensus sequence among non-specialist species. A dot indicates identity with the reference sequence and dashes indicate missing data. Species highlighted in green are specialists on cardiac glycoside (CG) containing plants.

https://doi.org/10.7554/eLife.48224.004
Supplementary file 2

Includes sources of sequence data used in this study, references for phylogenetic relationships a list of transgenic strains generated in this study.

https://doi.org/10.7554/eLife.48224.017
Supplementary file 3

Fasta nucleotide sequence alignment of all sequences used in this study.

https://doi.org/10.7554/eLife.48224.018
Supplementary file 4

Fasta amino acid sequence alignment of all sequences used in this study.

https://doi.org/10.7554/eLife.48224.019
Transparent reporting form
https://doi.org/10.7554/eLife.48224.020

Download links