Thursday, August 4, 2011

Helpful mutants usually won’t cooperate with each other

How genes interact – the phenomenon termed “epistasis” – is complex, yet must be understood to accurately interpret the contribution of individual genes to the development and behavior of the organism. Two groups -- Chou et al. and Khan et al. -- recently reported in Science the results of their remarkably similar experiments; they isolated individual, beneficial mutations in bacteria then tested how the mutations interact when possessed by the same bacterium.

Previous studies (cited by Chou) suggested that two deleterious mutations in the same pathway are generally less-than-additive but were greater-than-additive when they were in parallel pathways, which seems intuitive. Interactions between mutations in single genes were shown to depend on the background (other genes). However, epistasis among beneficial mutations in different genes was “unexplored”, and the focus of these new experiments. Evolution in laboratory conditions is initially rapid but quickly slows, which fits a model of mutually antagonistic beneficial mutations. Khan states the deceleration is due to either (1) negative epistasis or (2) because the most beneficial would “tend to be incorporated earlier owing to their faster spread and greater success in the face of competing beneficial mutations” (which sounds suspiciously convenient – the best arrive early – but they give a reference: Gerrish 1998 so you can look it up).

Khan grew E. coli with a glucose supplement for 20,000 generations, when they sequenced a clone and identified 45 mutations. They say other beneficial mutations arose but were lost due to “interference” with more-beneficial mutations (24, 26) (which seems to answer their question), or because they were “less able to evolve than the eventual winners” (33). They took the first 5 “that fixed … and whose spread coincided with the period of fastest adaptation” (arbitrary? How many generations was that, 200? 2,000?), which were, in order of appearance: rbs operon, topA, spoT, glmUS promoter, and pykF. Together, these 5 mutations increased fitness ~30%, accounting for ~80% of the fitness increase over the full 20,000 generations.

Then they produced 32 populations of E. coli, one for each possible combination of the 5 genes (Fig. 1 from Khan, copied here, shows the ancestral genotype at the top and the 32 combinations of mutations, with increasing fitness downward). They found that although each combination improved fitness, improvement was less than expected from a “multiplicative null model” (in which the individual fitness effects are multiplied).

Similarly, Chou selected a bacterium to grow with methanol as its sole carbon source and then replaced a key metabolic pathway with a less efficient pathway from a different bacterium. The resulting bacterium grew only one-third as fast as the original. Eight separate populations improved the efficiency over 600 generations. In the fittest strain, 9 mutations were identified; 3 were clearly related to the pathway while 6 others were deemed “unlikely” to contribute to fitness. To determine the interactions among these 3 genes, they constructed 16 strains, one with each combinations of these 3 mutations plus the WT allele. They also observed that combination strains were much less fit than expected if they were acting independently (our old favorite, the multiplicative null model).

These findings should influence our expectations in genetic studies in humans; e.g., two disease-risk alleles may have little additive or even negative effects (ergo protective?), contradicting the simple expectations of most biologists. [Chou distinguished beneficial from detrimental effects.] I wonder if these experiments might be misleading because they first fixed the individual mutants, thereby eliminating mutants that interact positively during selection.


Khan AI, Dinh DM, Schneider D, Lenski RE, Cooper TF."Negative epistasis between beneficial mutations in an evolving bacterial population. " Science. 2011 Jun 3;332(6034):1193-6.

The same issue has a commentary.

Thursday, May 19, 2011

Kidney disease linked to putative autoantigen + HLA

Idiopathic membranous nephropathy is a progressive disease involving the thickening of the basement membranes in glomeruli, which are key blood filtration units in the kidneys. Membranous nephropathy can be caused by exposure to toxins (gold, mercury, some medicines) or autoimmunity, such as lupus. Deposits of antibody-antigen (immune) complexes with complement components can be observed in the glomerulus.

These authors sought genetic associations in 556 biopsy-proven patients (British, French, & Dutch) by comparison of about 300,000 SNPs with matched healthy subjects. They found strong associations with a membrane protein previously implicated in autoimmunity, M-type phospholipase A2 receptor (PLA2R1, p~10E-28), and a histocompatibility gene (HLA-DQA1, p~10E-92). PLA2R1 is normally expressed in human glomeruli, exactly where immune complexes are found in membranous nephropathy patients. PLA2R1 was implicated only recently in autoantibody studies (Beck 2009) and is now known as the major autoantigen in idiopathic membranous nephropathy. The DQA1 association is not surprising, having been discovered by Vaughn et al using the relatively crude restriction fragment length polymorphism (RFLP) analysis and reported way back in 1989. The odds ratio for a single PLA2R1 risk allele is about 2 and for HLA-DQA1 about 6, modest but typical for such association studies. It is therefore astonishing that the risk to individuals possessing homozygous risk alleles at both loci is practically determinate – a 78.5-fold increased risk! – with 42 patients out of the 55 subjects possessing this combination.


How might this happen? The authors’ model is that perhaps the DQA1 molecule binds the PLA2R1 variant peptide and triggers T lymphocytes to help B cells make anti-PLA2R1 autoantibodies that bind to glomerular cells. They could have looked whether any PLA2R1 peptides has anchor residues that determine whether they can fit one of the 35 different allelic forms of DQA1. However, as Segelmark points out in the accompanying review, the strongest SNPs lie within the first introns for both PLA2R1 and DQA1 and (therefore) do not alter the amino acid sequence! The authors discount this, speculating that either the associated SNP is tightly linked to a variant that does change the protein. They could have sequenced the few subjects to identify any rare variant. Segelmark proposes as “more likely” that the SNP changes a regulatory sequence, such as a transcription factor binding site or a microRNA, that increases production of the proteins. A few more facts could help resolve these possibilities.
Risk HLA-DQA1 and PLA(2)R1 alleles in idiopathic membranous nephropathy. Stanescu HC, Arcos-Burgos M, Medlar A, Bockenhauer D, Kottgen A, Dragomirescu L, Voinescu C, Patel N, Pearce K, Hubank M, Stephens HA, Laundy V, Padmanabhan S, Zawadzka A, Hofstra JM, Coenen MJ, den Heijer M, Kiemeney LA, Bacq-Daian D, Stengel B, Powis SH, Brenchley P, Feehally J, Rees AJ, Debiec H, Wetzels JF, Ronco P, Mathieson PW, Kleta R. N Engl J Med. 2011 Feb 17;364(7):616-26.