Use of whole genome sequencing to estimate the mutation rate of Mycobacterium tuberculosis during latent infection

Tuberculosis poses a global health emergency, which has been compounded by the emergence of drug-resistant Mycobacterium tuberculosis (Mtb) strains. We used whole-genome sequencing to compare the accumulation of mutations in Mtb isolated from cynomolgus macaques with active, latent or reactivated disease. We sequenced 33 Mtb isolates from nine macaques with an average genome coverage of 93% and an average read depth of 117×. Based on the distribution of SNPs observed, we calculated the mutation rates for these disease states. We found a similar mutation rate during latency as during active disease or in a logarithmically growing culture over the same period of time. The pattern of polymorphisms suggests that the mutational burden in vivo is because of oxidative DNA damage. We show that Mtb continues to acquire mutations during disease latency, which may explain why isoniazid monotherapy for latent tuberculosis is a risk factor for the emergence of isoniazid resistance.     

Trinucleotide expansion in haploid germ cells by gap repair

Huntington disease (HD) is one of eight progressive neurodegenerative disorders in which the underlying mutation is a CAG expansion encoding a polyglutamine tract. The mechanism of trinucleotide expansion is poorly understood. Expansion is mediated by misaligned pairing of repeats and the inappropriate formation of DNA secondary structure as the duplex unpairs. It has never been clear, however, whether duplex unpairing occurs during mitotic replication or during strand-break repair. In simple organisms, trinucleotide expansion arises by replication slippage on either the leading or the lagging strand, homologous recombination, gene conversion, double-strand break repair and base excision repair; it is not clear which of these mechanisms is used in mammalian cells in vivo. We have followed heritable changes in CAG length in male transgenic mice. In germ cells, expansion is limited to the post-meiotic, haploid cell and therefore cannot involve mitotic replication or recombination between a homologous chromosome or a sister chromatid. Our data support a model in which expansion in the germ cells arises by gap repair and depends on a complex containing Msh2. Expansion occurs during gap-filling synthesis when DNA loops comprising the CAG trinucleotide repeats are sealed into the DNA strand.     

Unraveling adaptive evolution: how a single point mutation affects the protein coregulation network

Understanding the mechanisms of evolution requires identification of the molecular basis of the multiple (pleiotropic) effects of specific adaptive mutations. We have characterized the pleiotropic effects on protein levels of an adaptive single-base pair substitution in the coding sequence of a signaling pathway gene in the bacterium Pseudomonas fluorescens SBW25. We find 52 proteomic changes, corresponding to 46 identified proteins. None of these proteins is required for the adaptive phenotype. Instead, many are found within specific metabolic pathways associated with fitness-reducing (that is, antagonistic) effects of the mutation. The affected proteins fall within a single coregulatory network. The mutation 'rewires' this network by drawing particular proteins into tighter coregulating relationships. Although these changes are specific to the mutation studied, the quantitatively altered proteins are also affected in a coordinated way in other examples of evolution to the same niche.     

The MLH1 D132H variant is associated with susceptibility to sporadic colorectal cancer

Most susceptibility to colorectal cancer (CRC) is not accounted for by known risk factors. Because MLH1, MSH2 and MSH6 mutations underlie high-penetrance CRC susceptibility in hereditary nonpolyposis colon cancer (HNPCC), we hypothesized that attenuated alleles might also underlie susceptibility to sporadic CRC. We looked for gene variants associated with HNPCC in Israeli probands with familial CRC unstratified with respect to the microsatellite instability (MSI) phenotype. Association studies identified a new MLH1 variant (415G-->C, resulting in the amino acid substitution D132H) in approximately 1.3% of Israeli individuals with CRC self-described as Jewish, Christian and Muslim. MLH1 415C confers clinically significant susceptibility to CRC. In contrast to classic HNPCC, CRCs associated with MLH1 415C usually do not have the MSI defect, which is important for clinical mutation screening. Structural and functional analyses showed that the normal ATPase function of MLH1 is attenuated, but not eliminated, by the MLH1 415G-->C mutation. The new MLH1 variant confers a high risk of CRC and identifies a previously unrecognized mechanism in microsatellite-stable tumors. These studies suggest that variants of mismatch repair proteins with attenuated function may account for a higher proportion of susceptibility to sporadic microsatellite-stable CRC than previously assumed.     

CpG island methylator phenotype underlies sporadic microsatellite instability and is tightly associated with BRAF mutation in colorectal cancer

Aberrant DNA methylation of CpG islands has been widely observed in human colorectal tumors and is associated with gene silencing when it occurs in promoter areas. A subset of colorectal tumors has an exceptionally high frequency of methylation of some CpG islands, leading to the suggestion of a distinct trait referred to as 'CpG island methylator phenotype', or 'CIMP'. However, the existence of CIMP has been challenged. To resolve this continuing controversy, we conducted a systematic, stepwise screen of 195 CpG island methylation markers using MethyLight technology, involving 295 primary human colorectal tumors and 16,785 separate quantitative analyses. We found that CIMP-positive (CIMP+) tumors convincingly represent a distinct subset, encompassing almost all cases of tumors with BRAF mutation (odds ratio = 203). Sporadic cases of mismatch repair deficiency occur almost exclusively as a consequence of CIMP-associated methylation of MLH1 . We propose a robust new marker panel to classify CIMP+ tumors.     

High frequency of BRAF mutations in nevi

To evaluate the timing of mutations in BRAF (v-raf murine sarcoma viral oncogene homolog B1) during melanocytic neoplasia, we carried out mutation analysis on microdissected melanoma and nevi samples. We observed mutations resulting in the V599E amino-acid substitution in 41 of 60 (68%) melanoma metastases, 4 of 5 (80%) primary melanomas and, unexpectedly, in 63 of 77 (82%) nevi. These data suggest that mutational activation of the RAS/RAF/MAPK pathway in nevi is a critical step in the initiation of melanocytic neoplasia but alone is insufficient for melanoma tumorigenesis.     

Reciprocal crossover asymmetry and meiotic drive in a human recombination hot spot

Human DNA diversity arises ultimately from germline mutation that creates new haplotypes that can be reshuffled by meiotic recombination. Reciprocal crossover generates recombinant haplotypes but should not influence the frequencies of alleles in a population. We demonstrate crossover asymmetry at a recombination hot spot in the major histocompatibility complex, whereby reciprocal exchanges in sperm map to different locations in the hot spot. We identify a single-nucleotide polymorphism at the center of the hot spot and show that, when heterozygous, it seems sufficient to cause this asymmetry, apparently by influencing the efficiency of highly localized crossover initiation. As a consequence, crossovers in heterozygotes are accompanied by biased gene conversion, most likely occurring by gap repair, that can also affect nearby polymorphisms through repair of an extended gap. The result is substantial over-transmission of the recombination-suppressing allele and neighboring markers to crossover products. Computer simulations show that this meiotic drive, although weak at the population level, is sufficient to favor eventual fixation of the recombination-suppressing variant. These findings provide an explanation for the relatively uniform widths of human crossover hot spots and suggest that hot spots may be generally prone to extinction by meiotic drive.     

Inherited variants of MYH associated with somatic G:C-->T:A mutations in colorectal tumors

Inherited defects of base excision repair have not been associated with any human genetic disorder, although mutations of the genes mutM and mutY, which function in Escherichia coli base excision repair, lead to increased transversions of G:C to T:A. We have studied family N, which is affected with multiple colorectal adenomas and carcinoma but lacks an inherited mutation of the adenomatous polyposis coli gene (APC) that is associated with familial adenomatous polyposis. Here we show that 11 tumors from 3 affected siblings contain 18 somatic inactivating mutations of APC and that 15 of these mutations are G:C-->A transversions--a significantly greater proportion than is found in sporadic tumors or in tumors associated with familial adenomatous polyposis. Analysis of the human homolog of mutY, MYH, showed that the siblings were compound heterozygotes for the nonconservative missense variants Tyr165Cys and Gly382Asp. These mutations affect residues that are conserved in mutY of E. coli (Tyr82 and Gly253). Tyrosine 82 is located in the pseudo-helix-hairpin-helix (HhH) motif and is predicted to function in mismatch specificity. Assays of adenine glycosylase activity of the Tyr82Cys and Gly253Asp mutant proteins with 8-oxoG:A and G:A substrates show that their activity is reduced significantly. Our findings link the inherited variants in MYH to the pattern of somatic APC mutation in family N and implicate defective base excision repair in predisposition to tumors in humans.     

Identical point mutations of PMP-22 in Trembler-J mouse and Charcot-Marie-Tooth disease type 1A.

We have investigated the peripheral myelin protein gene, PMP-22, in a family with Charcot-Marie-Tooth disease type 1A (CMT1A). The DNA duplication commonly found in CMT1A was absent in this family, but strong linkage existed between the disease and the CMT1A marker VAW409R3 on chromosome 17p11.2. We found a point mutation in PMP-22 which was completely linked with the disease. The mutation, a proline for leucine substitution in the first putative transmembrane domain, is identical to that recently found in the Trembler-J mouse. The presence of this PMP-22 defect in this CMT1A family and the location of PMP-22 within the DNA duplication associated with CMT1A suggest that both structural alteration and overexpression of PMP-22 may lead to the disease.     

A canine DNM1 mutation is highly associated with the syndrome of exercise-induced collapse

Labrador retrievers are the most common dog breed in the world, with over 200,000 new kennel club registrations per year. The syndrome of exercise-induced collapse (EIC) in this breed is manifested by muscle weakness, incoordination and life-threatening collapse after intense exercise. Using a genome-wide microsatellite marker scan for linkage in pedigrees, we mapped the EIC locus to canine chromosome 9. We then used SNP association and haplotype analysis to fine map the locus, and identified a mutation in the dynamin 1 gene (DNM1) that causes an R256L substitution in a highly conserved region of the protein. This first documented mammalian DNM1 mutation is present at a high frequency in the breed and is a compelling candidate causal mutation for EIC, as the dynamin 1 protein has an essential role in neurotransmission and synaptic vesicle endocytosis.     

Gelsolin-derived familial amyloidosis caused by asparagine or tyrosine substitution for aspartic acid at residue 187.

Dominantly inherited familial amyloidosis, Finnish type (FAF) is caused by the accumulation of a 71-amino acid amyloidogenic fragment of mutant gelsolin (GSN). FAF is common in Finland but is very rare elsewhere. In Finland and in two American families, the mutation is a G654A transition leading to an Asp to Asn substitution at residue 187. We found the same mutation in a Dutch family but a Danish FAF family had a G654T mutation, predicting Asp to Tyr at residue 187. We also found the G654T transversion in a Czech family. Using GSN polymorphisms, different haplotypes were found in the Danish and Czech families. We conclude that substitution of the uncharged Asn or Tyr for the acidic Asp at residue 187 creates a conformation that may be preferentially amyloidogenic for GSN.     

Defective colour vision associated with a missense mutation in the human green visual pigment gene.

All red/green colour vision defects described so far have been associated with gross rearrangements within the red/green opsin gene array (Xq28). We now describe a male with severe deuteranomaly without such a rearrangement. A substitution of a highly conserved cysteine by arginine at position 203 in the green opsins presumably accounted for his colour vision defect. Surprisingly, this mutation was fairly common (2%) in the population but apparently was not always expressed. In analogy with nonexpression of some 5'green-red hybrid genes in persons with normal colour vision, we suggest that failure of manifestation occurs when the mutant gene is located at a distal (3') position among several green opsin genes. This mutation might also predispose to certain X-linked retinal dystrophies.     

Identification of SLC7A7, encoding y+LAT-1, as the lysinuric protein intolerance gene

Lysinuric protein intolerance (LPI; OMIM 222700) is a rare, recessive disorder with a worldwide distribution, but with a high prevalence in the Finnish population; symptoms include failure to thrive, growth retardation, muscle hypotonia and hepatosplenomegaly. A defect in the plasma membrane transport of dibasic amino acids has been demonstrated at the baso-lateral membrane of epithelial cells in small intestine and in renal tubules and in plasma membrane of cultured skin fibroblasts from LPI patients. The gene causing LPI has been assigned by linkage analysis to 14q11-13. Here we report mutations in SLC7A7 cDNA (encoding y+L amino acid transporter-1, y+LAT-1), which expresses dibasic amino-acid transport activity and is located in the LPI region, in 31 Finnish LPI patients and 1 Spanish patient. The Finnish patients are homozygous for a founder missense mutation leading to a premature stop codon. The Spanish patient is a compound heterozygote with a missense mutation in one allele and a frameshift mutation in the other. The frameshift mutation generates a premature stop codon, eliminating the last one-third of the protein. The missense mutation abolishes y+LAT-1 amino-acid transport activity when co-expressed with the heavy chain of the cell-surface antigen 4F2 (4F2hc, also known as CD98) in Xenopus laevis oocytes. Our data establish that mutations in SLC7A7 cause LPI.     

Mutations in a delta 8-delta 7 sterol isomerase in the tattered mouse and X-linked dominant chondrodysplasia punctata. jderry@immunex.com.

Tattered (Td) is an X-linked, semi-dominant mouse mutation associated with prenatal male lethality. Heterozygous females are small and at 4-5 days of age develop patches of hyperkeratotic skin where no hair grows, resulting in a striping of the coat in adults. Craniofacial anomalies and twisted toes have also been observed in some affected females. A potential second allele of Td has also been described. The phenotype of Td is similar to that seen in heterozygous females with human X-linked dominant chondrodysplasia punctata (CDPX2, alternatively known as X-linked dominant Conradi-Hünermann-Happle syndrome) as well as another X-linked, semi-dominant mouse mutation, bare patches (Bpa). The Bpa gene has recently been identified and encodes a protein with homology to 3beta-hydroxysteroid dehydrogenases that functions in one of the later steps of cholesterol biosynthesis. CDPX2 patients display skin defects including linear or whorled atrophic and pigmentary lesions, striated hyperkeratosis, coarse lusterless hair and alopecia, cataracts and skeletal abnormalities including short stature, rhizomelic shortening of the limbs, epiphyseal stippling and craniofacial defects (MIM 302960). We have now identified the defect in Td mice as a single amino acid substitution in the delta8-delta7 sterol isomerase emopamil binding protein (Ebp; encoded by Ebp in mouse) and identified alterations in human EBP in seven unrelated CDPX2 patients.     

Epistatic buffering of fitness loss in yeast double deletion strains

Interactions between deleterious mutations have been insufficiently studied, despite the fact that their strength and direction are critical for understanding the evolution of genetic recombination and the buildup of mutational load in populations. We compiled a list of 758 yeast gene deletions causing growth defects (from the Munich Information Center for Protein Sequences database and ref. 7). Using BY4741 and BY4742 single-deletion strains, we carried out 639 random crosses and assayed growth curves of the resulting progeny. We show that the maximum growth rate averaged over strains lacking deletions and those with double deletions is higher than that of strains with single deletions, indicating a positive epistatic effect. This tendency is shared by genes belonging to a variety of functional classes. Based on our data and former theoretical work, we suggest that epistasis is likely to diminish the negative effects of mutations when the ability to produce biomass at high rates contributes significantly to fitness.     

The peripheral myelin gene PMP-22/GAS-3 is duplicated in Charcot-Marie-Tooth disease type 1A.

Charcot-Marie-Tooth disease type 1A (CMT1A) is associated with a DNA duplication at chromosome 17p11.2. In view of the point mutation in the gene for peripheral myelin protein pmp-22/gas-3 in Trembler mice, a murine model for CMT1A, we have analysed whether this gene is altered in CMT1A. Here we show that the human homologue of the murine pmp-22 gene is located within the CMT1A DNA duplication, which is a direct repeat and does not interrupt the coding region of PMP-22. Expression of PMP-22 in CMT1A fibroblasts is similar to expression in control fibroblasts. Increased gene dosage or altered PMP-22 expression in the peripheral nervous system are therefore possible mechanisms by which PMP-22 is involved in CMT1A.