The direction of microsatellite mutations is dependent upon allele length

Microsatellites, comprising tandemly repeated short nucleotide sequences, are ubiquitous in eukaryotic genomes. Mutations within microsatellites are frequent, altering their overall length by insertion or deletion of a small number of repeat units, with a rate as high as 10(-3) in humans. Despite their high mutability, stable allele frequency distributions are typically observed for microsatellites in humans as well as other primates, although the mechanism maintaining these stable distributions remains unclear. Previous studies have suggested that microsatellite mutations occur more frequently in longer alleles and favour expansion. Generalizing these results has been hindered because the sample sizes were small, only a small subset of alleles for any marker was studied and the direction of mutation (expansion or contraction) was not rigorously determined. Here we examine 236 mutations at 122 tetranucleotide repeat markers and find that the rate of contraction mutations increases exponentially with allele size, whereas the rate of expansion mutations is constant across the entire allele distribution. The overall rate of expansion mutations does not differ from that of contractions. Our findings offer an explanation for the stationary allele distribution of microsatellites.     

Episodic evolution of pyrin in primates: human mutations recapitulate ancestral amino acid states

Familial Mediterranean fever (FMF; MIM 249100) is an autosomal recessive disease characterized by recurrent attacks of fever with synovial, pleural or peritoneal inflammation. The disease is caused by mutations in the gene encoding the pyrin protein. Human population studies have revealed extremely high allele frequencies for several different pyrin mutations, leading to the conclusion that the mutant alleles confer a selective advantage. Here we examine the ret finger protein (rfp) domain (which contains most of the disease-causing mutations) of pyrin during primate evolution. Amino acids that cause human disease are often present as wild type in other species. This is true at positions 653 (a novel mutation), 680, 681, 726, 744 and 761. For several of these human mutations, the mutant represents the reappearance of an ancestral amino acid state. Examination of lineage-specific dN/dS ratios revealed a pattern consistent with the signature of episodic positive selection. Our data, together with previous human population studies, indicate that selective pressures may have caused functional evolution of pyrin in humans and other primates.     

Mutations in the 70K peroxisomal membrane protein gene in Zellweger syndrome.

The peroxisomal membrane protein, with a relative molecular mass of 70,000 (M(r) 70K) (PMP70), is an important component of peroxisomal membranes and an ATP-binding cassette protein. To investigate its possible involvement in Zellweger syndrome (ZS), an inborn error of peroxisome assembly, we cloned and sequenced cDNAs for human PMP70 and mapped the gene to chromosome 1. Amongst 32 probands with ZS or related disorders, we found two mutant PMP70 alleles in single ZS probands from the same complementation group. One allele has a donor splice site mutation and the second a missense mutation. Our results suggest that PMP70 plays an important role in peroxisome biogenesis and that mutations in PMP70 may be responsible for a subset of ZS patients.     

Conserved noncoding sequences are selectively constrained and not mutation cold spots

Noncoding genetic variants are likely to influence human biology and disease, but recognizing functional noncoding variants is difficult. Approximately 3% of noncoding sequence is conserved among distantly related mammals, suggesting that these evolutionarily conserved noncoding regions (CNCs) are selectively constrained and contain functional variation. However, CNCs could also merely represent regions with lower local mutation rates. Here we address this issue and show that CNCs are selectively constrained in humans by analyzing HapMap genotype data. Specifically, new (derived) alleles of SNPs within CNCs are rarer than new alleles in nonconserved regions (P = 3 x 10(-18)), indicating that evolutionary pressure has suppressed CNC-derived allele frequencies. Intronic CNCs and CNCs near genes show greater allele frequency shifts, with magnitudes comparable to those for missense variants. Thus, conserved noncoding variants are more likely to be functional. Allele frequency distributions highlight selectively constrained genomic regions that should be intensively surveyed for functionally important variation.     

Extensive allelic variation and ultrashort telomeres in senescent human cells

By imposing a limit on the proliferative lifespan of most somatic cells, telomere erosion represents an innate mechanism for tumor suppression and may contribute to age-related disease. A detailed understanding of the pathways that link shortened telomeres to replicative senescence has been severely hindered by the inability of current methods to analyze telomere dynamics in detail. Here we describe single telomere length analysis (STELA), a PCR-based approach that accurately measures the full spectrum of telomere lengths from individual chromosomes. STELA analysis of human XpYp telomeres in fibroblasts identifies several features of telomere biology. We observe bimodal distributions of telomeres in normal fibroblasts; these distributions result from inter-allelic differences of up to 6.5 kb, indicating that unexpectedly large-scale differences in zygotic telomere length are maintained throughout development. Most telomeres shorten in a gradual fashion consistent with simple losses through end replication, and the rates of erosion are independent of allele size. Superimposed on this are occasional, more substantial changes in length, which may be the consequence of additional mutational mechanisms. Notably, some alleles show almost complete loss of TTAGGG repeats at senescence.     

Mutant chromatin remodeling protein SMARCAL1 causes Schimke immuno-osseous dysplasia.

Schimke immuno-osseous dysplasia (SIOD, MIM 242900) is an autosomal-recessive pleiotropic disorder with the diagnostic features of spondyloepiphyseal dysplasia, renal dysfunction and T-cell immunodeficiency. Using genome-wide linkage mapping and a positional candidate approach, we determined that mutations in SMARCAL1 (SWI/SNF2-related, matrix-associated, actin-dependent regulator of chromatin, subfamily a-like 1), are responsible for SIOD. Through analysis of data from persons with SIOD in 26 unrelated families, we observed that affected individuals from 13 of 23 families with severe disease had two alleles with nonsense, frameshift or splicing mutations, whereas affected individuals from 3 of 3 families with milder disease had a missense mutation on each allele. These observations indicate that some missense mutations allow retention of partial SMARCAL1 function and thus cause milder disease.     

Genotype-based screen for ENU-induced mutations in mouse embryonic stem cells

The ability to generate mutations is a prerequisite to functional genetic analysis. Despite a long history of using mice as a model system for genetic analysis, the scientific community has not generated a comprehensive collection of multiple alleles for most mouse genes. The chemical mutagen of choice for mouse has been N-ethyl-N-nitrosourea (ENU), an alkylating agent that mainly causes base substitutions in DNA, and therefore allows for recovery of complete and partial loss-, as well as gain-, of-function alleles . Specific locus tests designed to detect recessive mutations showed that ENU is the most efficient mutagen in mouse with an approximate mutation rate of 1 in 1,000 gametes. In fact, several genome-wide and region-specific screens based on phenotypes have been carried out. The anticipation of the completion of the human and mouse genome projects, however, now emphasizes genotype-driven genetics--from sequence to mutants. To take advantage of the mutagenicity of ENU and its ability to create allelic series of mutations, we have developed a complementary approach to generating mutations using mouse embryonic stem (ES) cells. We show that a high mutation frequency can be achieved and that modulating DNA-repair activities can enhance this frequency. The treated cells retain germline competency, thereby rendering this approach applicable for efficient generation of an allelic series of mutations pivotal to a fine-tuned dissection of biological pathways.     

Antagonism between DNA hypermethylation and enhancer-blocking activity at the H19 DMD is uncovered by CpG mutations

Imprinted expression at the H19-Igf2 locus depends on a differentially methylated domain (DMD) that acts both as a maternal-specific, methylation-sensitive insulator and as a paternal-specific site of hypermethylation. Four repeats in the DMD bind CCCTC-binding factor (CTCF) on the maternal allele and have been proposed to attract methylation on the paternal allele. We introduced point mutations into the DMD to deplete the repeats of CpGs while retaining CTCF-binding and enhancer-blocking activity. Maternal inheritance of the mutations left H19 expression and Igf2 imprinting intact, consistent with the idea that the DMD acts as an insulator. Conversely, paternal inheritance of these mutations disrupted maintenance of DMD methylation, resulting in biallelic H19 expression. Furthermore, an insulator was established on the paternally inherited mutated allele in vivo, reducing Igf2 expression and resulting in a 40% reduction in size of newborn offspring. Thus, the nine CpG mutations in the DMD showed that the two parental-specific roles of the H19 DMD, methylation maintenance and insulator assembly, are antagonistic.     

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.     

Germline mutations in the ribonuclease L gene in families showing linkage with HPC1.

Although prostate cancer is the most common non-cutaneous malignancy diagnosed in men in the United States, little is known about inherited factors that influence its genetic predisposition. Here we report that germline mutations in the gene encoding 2'-5'-oligoadenylate(2-5A)-dependent RNase L (RNASEL) segregate in prostate cancer families that show linkage to the HPC1 (hereditary prostate cancer 1) region at 1q24-25 (ref. 9). We identified RNASEL by a positional cloning/candidate gene method, and show that a nonsense mutation and a mutation in an initiation codon of RNASEL segregate independently in two HPC1-linked families. Inactive RNASEL alleles are present at a low frequency in the general population. RNASEL regulates cell proliferation and apoptosis through the interferon-regulated 2-5A pathway and has been suggested to be a candidate tumor suppressor gene. We found that microdissected tumors with a germline mutation showed loss of heterozygosity and loss of RNase L protein, and that RNASEL activity was reduced in lymphoblasts from heterozyogous individuals compared with family members who were homozygous with respect to the wildtype allele. Thus, germline mutations in RNASEL may be of diagnostic value, and the 2-5A pathway might provide opportunities for developing therapies for those with prostate cancer.     

A de novo paradigm for mental retardation

The per-generation mutation rate in humans is high. De novo mutations may compensate for allele loss due to severely reduced fecundity in common neurodevelopmental and psychiatric diseases, explaining a major paradox in evolutionary genetic theory. Here we used a family based exome sequencing approach to test this de novo mutation hypothesis in ten individuals with unexplained mental retardation. We identified and validated unique non-synonymous de novo mutations in nine genes. Six of these, identified in six different individuals, are likely to be pathogenic based on gene function, evolutionary conservation and mutation impact. Our findings provide strong experimental support for a de novo paradigm for mental retardation. Together with de novo copy number variation, de novo point mutations of large effect could explain the majority of all mental retardation cases in the population.     

MLH3: a DNA mismatch repair gene associated with mammalian microsatellite instability

DNA mismatch repair is important because of its role in maintaining genomic integrity and its association with hereditary non-polyposis colon cancer (HNPCC). To identify new human mismatch repair proteins, we probed nuclear extracts with the conserved carboxy-terminal MLH1 interaction domain. Here we describe the cloning and complete genomic sequence of MLH3, which encodes a new DNA mismatch repair protein that interacts with MLH1. MLH3 is more similar to mismatch repair proteins from yeast, plants, worms and bacteria than to any known mammalian protein, suggesting that its conserved sequence may confer unique functions in mice and humans. Cells in culture stably expressing a dominant-negative MLH3 protein exhibit microsatellite instability. Mlh3 is highly expressed in gastrointestinal epithelium and physically maps to the mouse complex trait locus colon cancer susceptibility I (Ccs1). Although we were unable to identify a mutation in the protein-coding region of Mlh3 in the susceptible mouse strain, colon tumours from congenic Ccs1 mice exhibit microsatellite instability. Functional redundancy among Mlh3, Pms1 and Pms2 may explain why neither Pms1 nor Pms2 mutant mice develop colon cancer, and why PMS1 and PMS2 mutations are only rarely found in HNPCC families.     

Characterization of single-nucleotide polymorphisms in coding regions of human genes.

A major goal in human genetics is to understand the role of common genetic variants in susceptibility to common diseases. This will require characterizing the nature of gene variation in human populations, assembling an extensive catalogue of single-nucleotide polymorphisms (SNPs) in candidate genes and performing association studies for particular diseases. At present, our knowledge of human gene variation remains rudimentary. Here we describe a systematic survey of SNPs in the coding regions of human genes. We identified SNPs in 106 genes relevant to cardiovascular disease, endocrinology and neuropsychiatry by screening an average of 114 independent alleles using 2 independent screening methods. To ensure high accuracy, all reported SNPs were confirmed by DNA sequencing. We identified 560 SNPs, including 392 coding-region SNPs (cSNPs) divided roughly equally between those causing synonymous and non-synonymous changes. We observed different rates of polymorphism among classes of sites within genes (non-coding, degenerate and non-degenerate) as well as between genes. The cSNPs most likely to influence disease, those that alter the amino acid sequence of the encoded protein, are found at a lower rate and with lower allele frequencies than silent substitutions. This likely reflects selection acting against deleterious alleles during human evolution. The lower allele frequency of missense cSNPs has implications for the compilation of a comprehensive catalogue, as well as for the subsequent application to disease association.     

Germline loss-of-function mutations in SPRED1 cause a neurofibromatosis 1-like phenotype

We report germline loss-of-function mutations in SPRED1 in a newly identified autosomal dominant human disorder. SPRED1 is a member of the SPROUTY/SPRED family of proteins that act as negative regulators of RAS->RAF interaction and mitogen-activated protein kinase (MAPK) signaling. The clinical features of the reported disorder resemble those of neurofibromatosis type 1 and consist of multiple café-au-lait spots, axillary freckling and macrocephaly. Melanocytes from a café-au-lait spot showed, in addition to the germline SPRED1 mutation, an acquired somatic mutation in the wild-type SPRED1 allele, indicating that complete SPRED1 inactivation is needed to generate a café-au-lait spot in this syndrome. This disorder is yet another member of the recently characterized group of phenotypically overlapping syndromes caused by mutations in the genes encoding key components of the RAS-MAPK pathway. To our knowledge, this is the first report of mutations in the SPRY (SPROUTY)/SPRED family of genes in human disease.     

A natural allele of Nxf1 suppresses retrovirus insertional mutations

Endogenous retroviruses have shaped the evolution of mammalian genomes. Host genes that control the effects of retrovirus insertions are therefore of great interest. The modifier-of-vibrator-1 locus (Mvb1) controls levels of correctly processed mRNA from genes mutated by endogenous retrovirus insertions into introns, including the Pitpn(vb) tremor mutation and the Eya1(BOR) model of human branchiootorenal syndrome. Positional complementation cloning identifies Mvb1 as the nuclear export factor Nxf1, providing an unexpected link between the mRNA export receptor and pre-mRNA processing. Population structure of the suppressive allele in wild Mus musculus castaneus suggests selective advantage. A congenic Mvb1(CAST) allele is a useful tool for modifying gene expression from existing mutations and could be used to manipulate engineered mutations containing retroviral elements.     

Exome sequencing identifies ACSF3 as a cause of combined malonic and methylmalonic aciduria

We used exome sequencing to identify the genetic basis of combined malonic and methylmalonic aciduria (CMAMMA). We sequenced the exome of an individual with CMAMMA and followed up with sequencing of eight additional affected individuals (cases). This included one individual who was identified and diagnosed by searching an exome database. We identify mutations in ACSF3, encoding a putative methylmalonyl-CoA and malonyl-CoA synthetase as a cause of CMAMMA. We also examined a canine model of CMAMMA, which showed pathogenic mutations in a predicted ACSF3 ortholog. ACSF3 mutant alleles occur with a minor allele frequency of 0.0058 in ～1,000 control individuals, predicting a CMAMMA population incidence of ～1:30,000. ACSF3 deficiency is the first human disorder identified as caused by mutations in a gene encoding a member of the acyl-CoA synthetase family, a diverse group of evolutionarily conserved proteins, and may emerge as one of the more common human metabolic disorders.     

Large expansion of the ATTCT pentanucleotide repeat in spinocerebellar ataxia type 10

Spinocerebellar ataxia type 10 (SCA10; MIM 603516; refs 1,2) is an autosomal dominant disorder characterized by cerebellar ataxia and seizures. The gene SCA10 maps to a 3.8-cM interval on human chromosome 22q13-qter (refs 1,2). Because several other SCA subtypes show trinucleotide repeat expansions, we examined microsatellites in this region. We found an expansion of a pentanucleotide (ATTCT) repeat in intron 9 of SCA10 in all patients in five Mexican SCA10 families. There was an inverse correlation between the expansion size, up to 22.5 kb larger than the normal allele, and the age of onset (r2=0.34, P=0.018). Analysis of 562 chromosomes from unaffected individuals of various ethnic origins (including 242 chromosomes from Mexican persons) showed a range of 10 to 22 ATTCT repeats with no evidence of expansions. Our data indicate that the new SCA10 intronic ATTCT pentanucleotide repeat in SCA10 patients is unstable and represents the largest microsatellite expansion found so far in the human genome.