The CpG dinucleotide and human genetic disease

Age at onset in genetic prion disease and the design of preventive clinical trials

OBJECTIVE

To determine whether preventive trials in genetic prion disease could be designed to follow presymptomatic mutation carriers to onset of disease.

METHODS

We assembled age at onset or death data from 1,094 individuals with high penetrance mutations in the prion protein gene (PRNP) in order to generate survival and hazard curves and test for genetic modifiers of age at onset. We used formulae and simulations to estimate statistical power for clinical trials.

RESULTS

Genetic prion disease age at onset varies over several decades for the most common mutations and neither sex, parent's age at onset, nor PRNP codon 129 genotype provided additional explanatory power to stratify trials. Randomized preventive trials would require hundreds or thousands of at-risk individuals in order to be statistically powered for an endpoint of clinical onset, posing prohibitive cost and delay and likely exceeding the number of individuals available for such trials.

CONCLUSION

The characterization of biomarkers suitable to serve as surrogate endpoints will be essential for the prevention of genetic prion disease. Parameters such as longer trial duration, increased enrollment, and the use of historical controls in a postmarketing study could provide opportunities for subsequent determination of clinical benefit.

ALPK1 missense pathogenic variant in five families leads to ROSAH syndrome, an ocular multisystem autosomal dominant disorder

Purpose

To identify the molecular cause in five unrelated families with a distinct autosomal dominant ocular systemic disorder we called ROSAH syndrome due to clinical features of retinal dystrophy, optic nerve edema, splenomegaly, anhidrosis, and migraine headache.

Methods

Independent discovery exome and genome sequencing in families 1, 2, and 3, and confirmation in families 4 and 5. Expression of wild-type messenger RNA and protein in human and mouse tissues and cell lines. Ciliary assays in fibroblasts from affected and unaffected family members.

Results

We found the heterozygous missense variant in the ɑ-kinase gene, ALPK1, (c.710C>T, [p.Thr237Met]), segregated with disease in all five families. All patients shared the ROSAH phenotype with additional low-grade ocular inflammation, pancytopenia, recurrent infections, and mild renal impairment in some. ALPK1 was notably expressed in retina, retinal pigment epithelium, and optic nerve, with immunofluorescence indicating localization to the basal body of the connecting cilium of the photoreceptors, and presence in the sweat glands. Immunocytofluorescence revealed expression at the centrioles and spindle poles during metaphase, and at the base of the primary cilium. Affected family member fibroblasts demonstrated defective ciliogenesis.

Conclusion

Heterozygosity for ALPK1, p.Thr237Met leads to ROSAH syndrome, an autosomal dominant ocular systemic disorder.

Mitochondrial stress triggers a pro-survival response through epigenetic modifications of nuclear DNA

Mitochondrial dysfunction represents an important cellular stressor and when intense and persistent cells must unleash an adaptive response to prevent their extinction. Furthermore, mitochondria can induce nuclear transcriptional changes and DNA methylation can modulate cellular responses to stress. We hypothesized that mitochondrial dysfunction could trigger an epigenetically mediated adaptive response through a distinct DNA methylation patterning. We studied cellular stress responses (i.e., apoptosis and autophagy) in mitochondrial dysfunction models. In addition, we explored nuclear DNA methylation in response to this stressor and its relevance in cell survival. Experiments in cultured human myoblasts revealed that intense mitochondrial dysfunction triggered a methylation-dependent pro-survival response. Assays done on mitochondrial disease patient tissues showed increased autophagy and enhanced DNA methylation of tumor suppressor genes and pathways involved in cell survival regulation. In conclusion, mitochondrial dysfunction leads to a "pro-survival" adaptive state that seems to be triggered by the differential methylation of nuclear genes.

Defining the Role of Nucleotide Flipping in Enzyme Specificity Using 19F NMR

A broad range of proteins employ nucleotide flipping to recognize specific sites in nucleic acids, including DNA glycosylases, which remove modified nucleobases to initiate base excision repair. Deamination, a pervasive mode of damage, typically generates lesions that are recognized by glycosylases as being foreign to DNA. However, deamination of 5-methylcytosine (mC) generates thymine, a canonical DNA base, presenting a challenge for damage recognition. Nevertheless, repair of mC deamination is important because the resulting G·T mispairs cause C → T transition mutations, and mC is abundant in all three domains of life. Countering this threat are three types of glycosylases that excise thymine from G·T mispairs, including thymine DNA glycosylase (TDG). These enzymes must minimize excision of thymine that is not generated by mC deamination, in A·T pairs and in polymerase-generated G·T mispairs. TDG preferentially removes thymine from DNA contexts in which cytosine methylation is prevalent, including CG and one non-CG site. This remarkable context specificity could be attained through modulation of nucleotide flipping, a reversible step that precedes base excision. We tested this idea using fluorine NMR and DNA containing 2'-fluoro-substituted nucleotides. We find that dT nucleotide flipping depends on DNA context and is efficient only in contexts known to feature cytosine methylation. We also show that a conserved Ala residue limits thymine excision by hindering nucleotide flipping. A linear free energy correlation reveals that TDG attains context specificity for thymine excision through modulation of nucleotide flipping. Our results provide a framework for characterizing nucleotide flipping in nucleic acids using ¹⁹F NMR.

Estimating the number of genetic mutations (hits) required for carcinogenesis based on the distribution of somatic mutations

Individual instances of cancer are primarily a result of a combination of a small number of genetic mutations (hits). Knowing the number of such mutations is a prerequisite for identifying specific combinations of carcinogenic mutations and understanding the etiology of cancer. We present a mathematical model for estimating the number of hits based on the distribution of somatic mutations. The model is fundamentally different from previous approaches, which are based on cancer incidence by age. Our somatic mutation based model is likely to be more robust than age-based models since it does not require knowing or accounting for the highly variable mutation rate, which can vary by over three orders of magnitude. In fact, we find that the number of somatic mutations at diagnosis is weakly correlated with age at cancer diagnosis, most likely due to the extreme variability in mutation rates between individuals. Comparing the distribution of somatic mutations predicted by our model to the actual distribution from 6904 tumor samples we estimate the number of hits required for carcinogenesis for 17 cancer types. We find that different cancer types exhibit distinct somatic mutational profiles corresponding to different numbers of hits. Why might different cancer types require different numbers of hits for carcinogenesis? The answer may provide insight into the unique etiology of different cancer types.

Cancer is primarily a result of genetic mutations. Each individual instance of cancer is initiated by a specific combination of a small number of mutations (hits). In trying to identify these combinations of mutations, it is important to know how many hits to look for. However, there are conflicting estimates for the number of hits. We present a fundamentally different model for estimating the number of hits. We found that the number hits ranges from two-eight depending on cancer type. These findings may provide insight into the unique characteristics of different cancer types.

Heat stress upregulates G-T mismatch binding activities in zebrafish (Danio rerio) embryos preexposed and nonexposed to a sublethal level of cadmium (Cd)

G-T mispair frequently appears in eukaryotic DNA due to the spontaneous deamination of 5-methylcytosine paired with guanine and is therefore an important target for DNA mismatch repair (MMR). Our earlier studies showed the downregulation of G-T binding activities in cadmium (Cd)-exposed (Danio rerio) embryos. Since elevation of water temperature was reported to increase Cd toxicity in zebrafish, this study explored whether heat stress affected zebrafish mismatch binding capacity in the absence or presence of Cd. Heat stress (37 °C for 30 min) induced heat shock protein 70 mRNA expression in embryos at 10 and 24 h post fertilization (hpf). Heat stress weakly upregulated normal G-T sensing machinery and inhibited G-T recognition activity in embryos preexposed to 3 μM Cd for 9 h. Either heat shock or a 23-h Cd treatment alone caused a 1.7-fold stimulation of G-T binding capacity in 24 hpf embryos and heat stress of Cd-preexposed embryos further enhanced G-T binding activity to 2.5 fold of control. Normal and Cd-downregulated loop binding activities in 10 and 24 hpf embryos were almost unreactive to heat shock. Heat stress-upregulated G-T sensing in nonexposed, but not in Cd-preexposed, 24 hpf embryos correlated with stronger gene activities encoding MMR-linked mismatch detecting factors MutS homolog 2 and 6 plus a higher DNA binding activity of the transcription factor Sp1 that regulates msh2/msh6 expression. Our results suggested the importance of heat shock response in facilitating the correction of G-T mismatch in developing zebrafish even under Cd exposure.

Methylation-Induced Hypermutation in Natural Populations of Bacteria

Methylation of DNA at the C-5 position of cytosine occurs in diverse organisms. This modification can increase the rate of C→T transitions at the methylated position. In Escherichia coli and related enteric bacteria, the inner C residues of the sequence CCWGG (W is A or T) are methylated by the Dcm enzyme. These sites are hot spots of mutation during rapid growth in the laboratory but not in nondividing cells, in which repair by the Vsr protein is effective. It has been suggested that hypermutation at these sites is a laboratory artifact and does not occur in nature. Many other methyltransferases, with a variety of specificities, can be found in bacteria, usually associated with restriction enzymes and confined to a subset of the population. Their methylation targets are also possible sites of hypermutation. Here, I show using whole-genome sequence data for thousands of isolates that there is indeed considerable hypermutation at Dcm sites in natural populations: their transition rate is approximately eight times the average. I also demonstrate hypermutability of targets of restriction-associated methyltransferases in several distantly related bacteria: methylation increases the transition rate by a factor ranging from 12 to 58. In addition, I demonstrate how patterns of hypermutability inferred from massive sequence data can be used to determine previously unknown methylation patterns and methyltransferase specificities.IMPORTANCE A common type of DNA modification, addition of a methyl group to cytosine (C) at carbon atom C-5, can greatly increase the rate of mutation of the C to a T. In mammals, methylation of CG sequences increases the rate of CG→TG mutations. It is unknown whether cytosine C-5 methylation increases the mutation rate in bacteria under natural conditions. I show that sites methylated by the Dcm enzyme exhibit an 8-fold increase in mutation rate in natural bacterial populations. I also show that modifications at other sites in various bacteria also increase the mutation rate, in some cases by a factor of forty or more. Finally, I demonstrate how this phenomenon can be used to infer sequence specificities of methylation enzymes.

Genome-wide epigenetic landscape of pig lincRNAs and their evolution during porcine domestication

AIM

We aimed to identify previously unreported long intergenic noncoding RNAs (lincRNAs) in the porcine liver, an important metabolic tissue, and further illustrate the epigenomic landscapes and the evolution of lincRNAs.

MATERIALS & METHODS

We used porcine omics data and comprehensively analyzed and identified lincRNAs and their methylation, expression and evolutionary patterns during pig domestication.

RESULTS

LincRNAs exhibit highly methylated promoter and downstream regions, as well as lower expression levels and higher tissue specificity than protein-coding genes. We identified a batch of lincRNAs with selection signals that are associated with pig domestication, which are more highly expressed in the liver than in other tissues (19:10/8/6/3/2/1/1). Interestingly, the lincRNA linc-sscg1779 and its target gene C6, which is crucial in liver metabolism, are differentially expressed during pig domestication.

CONCLUSION

Although they may originate from noisy transcripts, lincRNAs may be subjected to artificial selection. This phenomenon implies the functional importance of lincRNAs in pig domestication.

MBD4 guards against methylation damage and germ line deficiency predisposes to clonal hematopoiesis and early-onset AML

The tendency of 5-methylcytosine (5mC) to undergo spontaneous deamination has had a major role in shaping the human genome, and this methylation damage remains the primary source of somatic mutations that accumulate with age. How 5mC deamination contributes to cancer risk in different tissues remains unclear. Genomic profiling of 3 early-onset acute myeloid leukemias (AMLs) identified germ line loss of MBD4 as an initiator of 5mC-dependent hypermutation. MBD4-deficient AMLs display a 33-fold higher mutation burden than AML generally, with >95% being C>T in the context of a CG dinucleotide. This distinctive signature was also observed in sporadic cancers that acquired biallelic mutations in MBD4 and in Mbd4 knockout mice. Sequential sampling of germ line cases demonstrated repeated expansion of blood cell progenitors with pathogenic mutations in DNMT3A, a key driver gene for both clonal hematopoiesis and AML. Our findings reveal genetic and epigenetic factors that shape the mutagenic influence of 5mC. Within blood cells, this links methylation damage to the driver landscape of clonal hematopoiesis and reveals a conserved path to leukemia. Germ line MBD4 deficiency enhances cancer susceptibility and predisposes to AML.

Novel splice site IDUA gene mutation in Tunisian pedigrees with hurler syndrome

BACKGROUND

The mucopolysaccharidosis type I (MPS I) is a lysosomal storage disease resulting from the defective activity of the enzyme α-L-iduronidase (IDUA). The disease has three major clinical subtypes (severe Hurler syndrome, intermediate Hurler-Scheie syndrome and attenuated Scheie syndrome). We aim to identify the genetic variants in MPS I patients and to investigate the effect of the novel splice site mutation on splicing of IDUA- mRNA variability using bioinformatics tools.

METHODS

The IDUA mutations were determined in four MPS I patients from four families from Northern Tunisia, by amplifying and sequencing each of the IDUA exons and intron-exon junctions.

RESULTS

One novel splice site IDUA mutation, c.1650 + 1G > T in intron 11 and two previously reported mutations, p.A75T and p.R555H, were detected. The patients in families 1 and 2 who have the Hurler phenotype were homozygotes for the novel splice site mutation c.1650 + 1G > T. The patient in family 3, who also had the Hurler phenotype, was a compound heterozygote for the novel splice site mutation c.1650 + 1G > T and for the previously reported missense mutation p.A75T. The patient in family 4 who had the Hurler-Scheie phenotype was a compound heterozygote for the novel splice site mutation c.1650 + 1G > T and for the previously reported missense mutation p.R555H. In addition, four known IDUA polymorphisms were identified. Bioinformatics tools allowed us to associate the variant c.1650 + 1G > T with the severe clinical phenotype of MPS I. This variant affects the essential nucleotide + 1 (G to T) of the donor splice site of IDUA intron 11. The G > T in intron 11 leads to wild type donor site broken with minus 19.97% value compared to normal value with 0%, hence the new splice site acceptor has plus 5.59%.

CONCLUSIONS

The present findings indicate that the identified mutations facilitate the accurate carrier detection (genetic counseling of at-risk relatives) and the molecular prenatal diagnosis in Tunisia.

KLK10 exon 3 unmethylated PCR product concentration: a new potential early diagnostic marker in ovarian cancer? - A pilot study

BACKGROUND

KLK10 exon 3 hypermethylation correlated to tumor-specific lack of KLK10 expression in cancer cell lines and primary tumors. In the present study we investigate the possible role of KLK10 exon 3 methylation in ovarian tumor diagnosis and prognosis.

RESULTS

Qualitative methylation-specific PCR (MSP) results did not show statistically significant differences in patient group samples (normal and tumor) where all samples were positive only for the unmethylated-specific PCR except for two malignant samples that were either doubly positive (serous carcinoma) or doubly negative (Sertoli-Leydig cell tumor) for the two MSP tests. However, KLK10 exon 3 unmethylated PCR product concentration (ng/μl) showed statistically significant differences in benign and malignant patient group samples; mean ± SD (n): tumor: 0.077 ± 0.035 (14) and 0.047 ± 0.021 (15), respectively, p-value = 0.011; and normal: 0.094 ± 0.039 (7) and 0.046 ± 0.027 (6), respectively, p-value = 0.031. Moreover, ROC curve analysis of KLK10 exon 3 unmethylated PCR product concentration in overall patient group samples showed good diagnostic ability (AUC = 0.778; p-value = 0.002). Patient survival (living and died) showed statistically significant difference according to preoperative serum CA125 concentration (U/ml); median (n): 101.25 (10) and 1252 (5), respectively, p-value = 0.037, but not KLK10 exon 3 unmethylated PCR product concentration (ng/μl) in overall malignant patient samples; mean ± SD (n): 0.042 ± 0.015 (14) and 0.055 ± 0.032 (7), p-value = 0.228.

CONCLUSION

To the best of our knowledge, this is the first report on KLK10 exon 3 unmethylated PCR product concentration as potential early epigenetic diagnostic marker in primary ovarian tumors. Taken into account the limitations in our study (small sample size and semi-quantitative PCR product analysis) further studies are strongly recommended.

Missense mutations in the gene encoding prothrombin corresponding to Arg596 cause antithrombin resistance and thrombomodulin resistance

Antithrombin (AT) and thrombomodulin (TM) play important roles in the process of natural anticoagulation in vivo. Recently, we reported that the prothrombin Yukuhashi mutation (p.Arg596Leu) was associated with AT and TM resistance-related thrombophilia. To assess the AT and TM resistances associated with other missense mutations by single base substitution in the Arg596 codon, we generated recombinant variants (596Gln, 596Trp, 596Gly, and 596Pro) and investigated the effects on AT and TM anticoagulant functions. All variants except 596Pro were secreted in amounts comparable to that of the wild-type but exhibited variable procoagulant activities. After a 30-minute inactivation by AT, the relative residual activity of wild-type thrombin decreased to 15 ± 4.0%, in contrast to values of all variants were maintained at above 80%. The thrombin–AT complex formation, as determined by enzyme-linked immunosorbent assay, was reduced with all tested variants in the presence and absence of heparin. In the presence of soluble TM (sTM), the relative fibrinogen clotting activity of wild-type thrombin decreased to 16 ± 0.12%, whereas that of tested variants was 37%–56%. In a surface plasmon resonance assay, missense Arg596 mutations reduced thrombin–TM affinity to an extent similar to the reduction of fibrinogen clotting inhibition. In the presence of sTM or cultured endothelial-like cells, APC generation was enhanced differently by variant thrombins in a thrombin–TM affinity- dependent manner. These data indicate that prothrombin Arg596 missense mutations lead to AT and TM resistance in the variant thrombins and suggest that prothrombin Arg596 is important for AT- and TM- mediated anticoagulation.

Presence of a consensus DNA motif at nearby DNA sequence of the mutation susceptible CG nucleotides

Complexity in tissues affected by cancer arises from somatic mutations and epigenetic modifications in the genome. The mutation susceptible hotspots present within the genome indicate a non-random nature and/or a position specific selection of mutation. An association exists between the occurrence of mutations and epigenetic DNA methylation. This study is primarily aimed at determining mutation status, and identifying a signature for predicting mutation prone zones of tumor suppressor (TS) genes. Nearby sequences from the top five positions having a higher mutation frequency in each gene of 42 TS genes were selected from a cosmic database and were considered as mutation prone zones. The conserved motifs present in the mutation prone DNA fragments were identified. Molecular docking studies were done to determine putative interactions between the identified conserved motifs and enzyme methyltransferase DNMT1. Collective analysis of 42 TS genes found GC as the most commonly replaced and AT as the most commonly formed residues after mutation. Analysis of the top 5 mutated positions of each gene (210 DNA segments for 42 TS genes) identified that CG nucleotides of the amino acid codons (e.g., Arginine) are most susceptible to mutation, and found a consensus DNA "T/AGC/GAGGA/TG" sequence present in these mutation prone DNA segments. Similar to TS genes, analysis of 54 oncogenes not only found CG nucleotides of the amino acid Arg as the most susceptible to mutation, but also identified the presence of similar consensus DNA motifs in the mutation prone DNA fragments (270 DNA segments for 54 oncogenes) of oncogenes. Docking studies depicted that, upon binding of DNMT1 methylates to this consensus DNA motif (C residues of CpG islands), mutation was likely to occur. Thus, this study proposes that DNMT1 mediated methylation in chromosomal DNA may decrease if a foreign DNA segment containing this consensus sequence along with CG nucleotides is exogenously introduced to dividing cancer cells.

The population genetics of human disease: The case of recessive, lethal mutations

Do the frequencies of disease mutations in human populations reflect a simple balance between mutation and purifying selection? What other factors shape the prevalence of disease mutations? To begin to answer these questions, we focused on one of the simplest cases: recessive mutations that alone cause lethal diseases or complete sterility. To this end, we generated a hand-curated set of 417 Mendelian mutations in 32 genes reported to cause a recessive, lethal Mendelian disease. We then considered analytic models of mutation-selection balance in infinite and finite populations of constant sizes and simulations of purifying selection in a more realistic demographic setting, and tested how well these models fit allele frequencies estimated from 33,370 individuals of European ancestry. In doing so, we distinguished between CpG transitions, which occur at a substantially elevated rate, and three other mutation types. Intriguingly, the observed frequency for CpG transitions is slightly higher than expectation but close, whereas the frequencies observed for the three other mutation types are an order of magnitude higher than expected, with a bigger deviation from expectation seen for less mutable types. This discrepancy is even larger when subtle fitness effects in heterozygotes or lethal compound heterozygotes are taken into account. In principle, higher than expected frequencies of disease mutations could be due to widespread errors in reporting causal variants, compensation by other mutations, or balancing selection. It is unclear why these factors would have a greater impact on disease mutations that occur at lower rates, however. We argue instead that the unexpectedly high frequency of disease mutations and the relationship to the mutation rate likely reflect an ascertainment bias: of all the mutations that cause recessive lethal diseases, those that by chance have reached higher frequencies are more likely to have been identified and thus to have been included in this study. Beyond the specific application, this study highlights the parameters likely to be important in shaping the frequencies of Mendelian disease alleles.

What determines the frequencies of disease mutations in human populations? To begin to answer this question, we focus on one of the simplest cases: mutations that cause completely recessive, lethal Mendelian diseases. We first review theory about what to expect from mutation and selection in a population of finite size and generate predictions based on simulations using a plausible demographic scenario of recent human evolution. For a highly mutable type of mutation, transitions at CpG sites, we find that the predictions are close to the observed frequencies of recessive lethal disease mutations. For less mutable types, however, predictions substantially under-estimate the observed frequency. We discuss possible explanations for the discrepancy and point to a complication that, to our knowledge, is not widely appreciated: that there exists ascertainment bias in disease mutation discovery. Specifically, we suggest that alleles that have been identified to date are likely the ones that by chance have reached higher frequencies and are thus more likely to have been mapped. More generally, our study highlights the factors that influence the frequencies of Mendelian disease alleles.

Epigenetic variation between urban and rural populations of Darwin's finches

Background

The molecular basis of evolutionary change is assumed to be genetic variation. However, growing evidence suggests that epigenetic mechanisms, such as DNA methylation, may also be involved in rapid adaptation to new environments. An important first step in evaluating this hypothesis is to test for the presence of epigenetic variation between natural populations living under different environmental conditions.

Results

In the current study we explored variation between populations of Darwin’s finches, which comprise one of the best-studied examples of adaptive radiation. We tested for morphological, genetic, and epigenetic differences between adjacent “urban” and “rural” populations of each of two species of ground finches, Geospiza fortis and G. fuliginosa, on Santa Cruz Island in the Galápagos. Using data collected from more than 1000 birds, we found significant morphological differences between populations of G. fortis, but not G. fuliginosa. We did not find large size copy number variation (CNV) genetic differences between populations of either species. However, other genetic variants were not investigated. In contrast, we did find dramatic epigenetic differences between the urban and rural populations of both species, based on DNA methylation analysis. We explored genomic features and gene associations of the differentially DNA methylated regions (DMR), as well as their possible functional significance.

Conclusions

In summary, our study documents local population epigenetic variation within each of two species of Darwin’s finches.

Electronic supplementary material

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

Biological function in the twilight zone of sequence conservation

Strong DNA conservation among divergent species is an indicator of enduring functionality. With weaker sequence conservation we enter a vast ‘twilight zone’ in which sequence subject to transient or lower constraint cannot be distinguished easily from neutrally evolving, non-functional sequence. Twilight zone functional sequence is illuminated instead by principles of selective constraint and positive selection using genomic data acquired from within a species’ population. Application of these principles reveals that despite being biochemically active, most twilight zone sequence is not functional.

Methyl-dependent and spatial-specific DNA recognition by the orthologous transcription factors human AP-1 and Epstein-Barr virus Zta

Activator protein 1 (AP-1) is a transcription factor that recognizes two versions of a 7-base pair response element, either 5΄-TGAGTCA-3΄ or 5΄-MGAGTCA-3΄ (where M = 5-methylcytosine). These two elements share the feature that 5-methylcytosine and thymine both have a methyl group in the same position, 5-carbon of the pyrimidine, so each of them has two methyl groups at nucleotide positions 1 and 5 from the 5΄ end, resulting in four methyl groups symmetrically positioned in duplex DNA. Epstein-Barr Virus Zta is a key transcriptional regulator of the viral lytic cycle that is homologous to AP-1. Zta recognizes several methylated Zta-response elements, including meZRE1 (5΄-TGAGMCA-3΄) and meZRE2 (5΄-TGAGMGA-3΄), where a methylated cytosine occupies one of the inner thymine residues corresponding to the AP-1 element, resulting in the four spatially equivalent methyl groups. Here, we study how AP-1 and Zta recognize these methyl groups within their cognate response elements. These methyl groups are in van der Waals contact with a conserved di-alanine in AP-1 dimer (Ala265 and Ala266 in Jun), or with the corresponding Zta residues Ala185 and Ser186 (via its side chain carbon Cβ atom). Furthermore, the two ZRE elements differ at base pair 6 (C:G versus G:C), forming a pseudo-symmetric sequence (meZRE1) or an asymmetric sequence (meZRE2). In vitro DNA binding assays suggest that Zta has high affinity for all four sequences examined, whereas AP-1 has considerably reduced affinity for the asymmetric sequence (meZRE2). We ascribe this difference to Zta Ser186 (a unique residue for Zta) whose side chain hydroxyl oxygen atom interacts with the two half sites differently, whereas the corresponding Ala266 of AP-1 Jun protein lacks such flexibility. Our analyses demonstrate a novel mechanism of 5mC/T recognition in a methylation-dependent, spatial and sequence-specific approach by basic leucine-zipper transcriptional factors.

DNA repair and mutations during quiescence in yeast

Life is maintained through alternating phases of cell division and quiescence. The causes and consequences of spontaneous mutations have been extensively explored in proliferating cells, and the major sources include errors of DNA replication and DNA repair. The foremost consequences are genetic variations within a cell population that can lead to heritable diseases and drive evolution. While most of our knowledge on DNA damage response and repair has been gained through cells actively dividing, it remains essential to also understand how DNA damage is metabolized in cells which are not dividing. In this review, we summarize the current knowledge concerning the type of lesions that arise in non-dividing budding and fission yeast cells, as well as the pathways used to repair them. We discuss the contribution of these models to our current understanding of age-related pathologies.

UV-induced damage to DNA: effect of cytosine methylation on pyrimidine dimerization

Methylation/demethylation of cytosine plays an important role in epigenetic signaling, the reversibility of epigenetic modifications offering important opportunities for targeted therapies. Actually, methylated sites have been correlated with mutational hotspots detected in skin cancers. The present brief review discusses the physicochemical parameters underlying the specific ultraviolet-induced reactivity of methylated cytosine. It focuses on dimerization reactions giving rise to cyclobutane pyrimidine dimers and pyrimidine (6–4) pyrimidone adducts. According to recent studies, four conformational and electronic factors that are affected by cytosine methylation may control these reactions: the red-shift of the absorption spectrum, the lengthening of the excited state lifetime, changes in the sugar puckering modifying the stacking between reactive pyrimidines and an increase in the rigidity of duplexes favoring excitation energy transfer toward methylated pyrimidines.

Mechanisms of DNA damage, repair, and mutagenesis

Living organisms are continuously exposed to a myriad of DNA damaging agents that can impact health and modulate disease-states. However, robust DNA repair and damage-bypass mechanisms faithfully protect the DNA by either removing or tolerating the damage to ensure an overall survival. Deviations in this fine-tuning are known to destabilize cellular metabolic homeostasis, as exemplified in diverse cancers where disruption or deregulation of DNA repair pathways results in genome instability. Because routinely used biological, physical and chemical agents impact human health, testing their genotoxicity and regulating their use have become important. In this introductory review, we will delineate mechanisms of DNA damage and the counteracting repair/tolerance pathways to provide insights into the molecular basis of genotoxicity in cells that lays the foundation for subsequent articles in this issue. Environ. Mol. Mutagen. 58:235-263, 2017. © 2017 Wiley Periodicals, Inc.

Identification of a missense HOXD13 mutation in a Chinese family with syndactyly type I-c using exome sequencing

Syndactyly is one of the most common hereditary limb malformations, and is characterized by the fusion of specific fingers and/or toes. Syndactyly type I‑c is associated with bilateral cutaneous or bony webbing of the third and fourth fingers and occasionally of the third to fifth fingers, with normal feet. The aim of the present study was to identify the genetic basis of syndactyly type I‑c in four generations of a Chinese Han family by exome sequencing. Exome sequencing was conducted in the proband of the family, followed by direct sequencing of other family members of the same ancestry, as well as 100 ethnically‑matched, unrelated normal controls. A missense mutation, c.917G>A (p.R306Q), was identified in the homeobox D13 gene (HOXD13). Sanger sequencing verified the presence of this mutation in all of the affected family members. By contrast, this mutation was absent in the unaffected family members and the 100 ethnically‑matched normal controls. The results suggest that the c.917G>A (p.R306Q) mutation in the HOXD13 gene, may be responsible for syndactyly type I‑c in this family. Exome sequencing may therefore be a powerful tool for identifying mutations associated with syndactyly, which is a disorder with high genetic and clinical heterogeneity. The results provide novel insights into the etiology and diagnosis of syndactyly, and may influence genetic counseling and the clinical management of the disease.

Mapping genotype-phenotype associations of nsSNPs in coiled-coil oligomerization domains of the human proteome

We assessed the impact of disease mutations (DMs) versus polymorphisms (PYs) in coiled-coil (CC) domains in UniProt by modeling the structural and functional impact of variants in silico with the CC prediction program Multicoil. The structural impact of variants was evaluated with respect to three main metrics: the oligomerization score-to determine whether the variant is stabilizing or destabilizing-the oligomerization state, and the register-specific score. The functional impact was queried indirectly in several ways. First, we examined marginally stable CCs that were either stabilized or destabilized by the variant. Second, we looked for variants that altered the register of the wild-type CC near wild-type irregularities of likely functional importance, such as skips and stammers. Third, we searched for variants that altered the oligomerization state of the CC. DMs tended to be more destabilizing than PYs; but interestingly, PYs were more frequently associated with predicted changes in the oligomerization state. The functional impact was also queried by testing the association of CC variants with multiple phenotypes, that is, pleiotropy. Mutations in CC regions of proteins cause 155 different phenotypes and are more frequently associated with pleiotropy than proteins in general. Importantly, the CC region itself often encodes the pleiotropy.

The axon degeneration gene SARM1 is evolutionarily distinct from other TIR domain-containing proteins

Many forms of neurodegenerative disease are characterized by Wallerian degeneration, an active program of axonal destruction. Recently, the important player which enacts Wallerian degeneration was discovered, the multidomain protein SARM1. Since the SARM1 protein has classically been thought of as an innate immune molecule, its role in Wallerian degeneration has raised questions on the evolutionary forces acting on it. Here, we synthesize a picture of SARM1's evolution through various organisms by examining the molecular and genetic changes of SARM1 and the genes around it. Using proteins that possess domains homologous to SARM1, we established distances and Ka/Ks values through 5671 pairwise species-species comparisons. We demonstrate that SARM1 diverged across species in a pattern similar to other SAM domain-containing proteins. This is surprising, because it was expected that SARM1 would behave more like its TIR domain relatives. Going along with this divorce from TIR, we also noted that SARM1's TIR is under stronger purifying selection than the rest of the TIR domain-containing proteins (remaining highly conserved). In addition, SARM1's synteny analysis reveals that the surrounding gene cluster is highly conserved, functioning as a potential nexus of gene functionality across species. Taken together, SARM1 demonstrates a unique evolutionary pattern, separate from the TIR domain protein family.

Differences in the rare variant spectrum among human populations

Mutations occur at vastly different rates across the genome, and populations, leading to differences in the spectrum of segregating polymorphisms. Here, we investigate variation in the rare variant spectrum in a sample of human genomes representing all major world populations. We find at least two distinct signatures of variation. One, consistent with a previously reported signature is characterized by an increased rate of TCC>TTC mutations in people from Western Eurasia and South Asia, likely related to differences in the rate, or efficiency of repair, of damage due to deamination of methylated guanine. We describe the geographic extent of this signature and show that it is detectable in the genomes of ancient, but not archaic humans. The second signature is private to certain Native American populations, and is concentrated at CpG sites. We show that this signature is not driven by differences in the CpG mutation rate, but is a result of the fact that highly mutable CpG sites are more likely to undergo multiple independent mutations across human populations, and the spectrum of such mutations is highly sensitive to recent demography. Both of these effects dramatically affect the spectrum of rare variants across human populations, and should be taken into account when using mutational clocks to make inference about demography.

Transthyretin V122I (pV142I)* cardiac amyloidosis: an age-dependent autosomal dominant cardiomyopathy too common to be overlooked as a cause of significant heart disease in elderly African Americans

Since the identification of a valine-to-isoleucine substitution at position 122 (TTR V122I; pV142I) in the transthyretin (TTR)-derived fibrils extracted from the heart of a patient with late-onset cardiac amyloidosis, it has become clear that the amyloidogenic mutation and the disease occur almost exclusively in individuals of identifiable African descent. In the United States, the amyloidogenic allele frequency is 0.0173 and is carried by 3.5% of community-dwelling African Americans. Genotyping across Africa indicates that the origin of the allele is in the West African countries that were the major source of the slave trade to North America. At autopsy, the allele was found to be associated with cardiac TTR amyloid deposition in all the carriers after age 65 years; however, the clinical penetrance varies, resulting in substantial heart disease in some carriers and few symptoms in others. The allele has been found in 10% of African Americans older than age 65 with severe congestive heart failure. At this time there are potential forms of therapy in clinical trials. The combination of a highly accurate genetic test and the potential for specific therapy demands a greater awareness of this autosomal dominant, age-dependent cardiac disease in the cardiology community.Genet Med advance online publication 19 January 2017.

Statistical Methods for Identifying Sequence Motifs Affecting Point Mutations

Mutation processes differ between types of point mutation, genomic locations, cells, and biological species. For some point mutations, specific neighboring bases are known to be mechanistically influential. Beyond these cases, numerous questions remain unresolved, including: what are the sequence motifs that affect point mutations? How large are the motifs? Are they strand symmetric? And, do they vary between samples? We present new log-linear models that allow explicit examination of these questions, along with sequence logo style visualization to enable identifying specific motifs. We demonstrate the performance of these methods by analyzing mutation processes in human germline and malignant melanoma. We recapitulate the known CpG effect, and identify novel motifs, including a highly significant motif associated with A[Formula: see text]G mutations. We show that major effects of neighbors on germline mutation lie within [Formula: see text] of the mutating base. Models are also presented for contrasting the entire mutation spectra (the distribution of the different point mutations). We show the spectra vary significantly between autosomes and X-chromosome, with a difference in T[Formula: see text]C transition dominating. Analyses of malignant melanoma confirmed reported characteristic features of this cancer, including statistically significant strand asymmetry, and markedly different neighboring influences. The methods we present are made freely available as a Python library https://bitbucket.org/pycogent3/mutationmotif.

Activation induced deaminase mutational signature overlaps with CpG methylation sites in follicular lymphoma and other cancers

Follicular lymphoma (FL) is an uncurable cancer characterized by progressive severity of relapses. We analyzed sequence context specificity of mutations in the B cells from a large cohort of FL patients. We revealed substantial excess of mutations within a novel hybrid nucleotide motif: the signature of somatic hypermutation (SHM) enzyme, Activation Induced Deaminase (AID), which overlaps the CpG methylation site. This finding implies that in FL the SHM machinery acts at genomic sites containing methylated cytosine. We identified the prevalence of this hybrid mutational signature in many other types of human cancer, suggesting that AID-mediated, CpG-methylation dependent mutagenesis is a common feature of tumorigenesis.

The Impact of Linked Selection in Chimpanzees: A Comparative Study

Levels of nucleotide diversity vary greatly across the genomes of most species owing to multiple factors. These include variation in the underlying mutation rates, as well as the effects of both direct and linked selection. Fundamental to interpreting the relative importance of these forces is the common observation of a strong positive correlation between nucleotide diversity and recombination rate. While indeed observed in humans, the interpretation of this pattern has been difficult in the absence of high-quality polymorphism data and recombination maps in closely related species. Here, we characterize genetic features driving nucleotide diversity in Western chimpanzees using a recently generated whole genome polymorphism data set. Our results suggest that recombination rate is the primary predictor of nucleotide variation with a strongly positive correlation. In addition, telomeric distance, regional GC-content, and regional CpG-island content are strongly negatively correlated with variation. These results are compared with humans, with both similarities and differences interpreted in the light of the estimated effective population sizes of the two species as well as their strongly differing recent demographic histories.

Acute Intermittent Porphyria: Predicted Pathogenicity of HMBS Variants Indicates Extremely Low Penetrance of the Autosomal Dominant Disease

Acute intermittent porphyria results from hydroxymethylbilane synthase (HMBS) mutations that markedly decrease HMBS enzymatic activity. This dominant disease is diagnosed when heterozygotes have life-threatening acute attacks, while most heterozygotes remain asymptomatic and undiagnosed. Although >400 HMBS mutations have been reported, the prevalence of pathogenic HMBS mutations in genomic/exomic databases, and the actual disease penetrance are unknown. Thus, we interrogated genomic/exomic databases, identified non-synonymous variants (NSVs) and consensus splice-site variants (CSSVs) in various demographic/racial groups, and determined the NSV's pathogenicity by prediction algorithms and in vitro expression assays. Caucasians had the most: 58 NSVs and two CSSVs among ∼92,000 alleles, a 0.00575 combined allele frequency. In silico algorithms predicted 14 out of 58 NSVs as "likely-pathogenic." In vitro expression identified 10 out of 58 NSVs as likely-pathogenic (seven predicted in silico), which together with two CSSVs had a combined allele frequency of 0.00056. Notably, six presumably pathogenic mutations/NSVs in the Human Gene Mutation Database were benign. Compared with the recent prevalence estimate of symptomatic European heterozygotes (∼0.000005), the prevalence of likely-pathogenic HMBS mutations among Caucasians was >100 times more frequent. Thus, the estimated penetrance of acute attacks was ∼1% of heterozygotes with likely-pathogenic mutations, highlighting the importance of predisposing/protective genes and environmental modifiers that precipitate/prevent the attacks.

Analysis of protein-coding genetic variation in 60,706 humans

Large-scale reference data sets of human genetic variation are critical for the medical and functional interpretation of DNA sequence changes. Here we describe the aggregation and analysis of high-quality exome (protein-coding region) DNA sequence data for 60,706 individuals of diverse ancestries generated as part of the Exome Aggregation Consortium (ExAC). This catalogue of human genetic diversity contains an average of one variant every eight bases of the exome, and provides direct evidence for the presence of widespread mutational recurrence. We have used this catalogue to calculate objective metrics of pathogenicity for sequence variants, and to identify genes subject to strong selection against various classes of mutation; identifying 3,230 genes with near-complete depletion of predicted protein-truncating variants, with 72% of these genes having no currently established human disease phenotype. Finally, we demonstrate that these data can be used for the efficient filtering of candidate disease-causing variants, and for the discovery of human 'knockout' variants in protein-coding genes.

The prevalence and distribution of the amyloidogenic transthyretin (TTR) V122I allele in Africa

Background

Transthyretin (TTR) pV142I (rs76992529‐A) is one of the 113 variants in the human TTR gene associated with systemic amyloidosis. It results from a G to A transition at a CG dinucleotide in the codon for amino acid 122 of the mature protein (TTR V122I). The allele frequency is 0.0173 in African Americans.

Methods

PCR‐based assays to genotype 2767 DNA samples obtained from participants in genetic studies from various African populations supplemented with sequencing data from 529 samples within the 1000 Genomes Project.

Results

The rs76992529‐A variant allele was most prevalent (allele frequency 0.0253) in the contiguous West African countries of Sierra Leone, Guinea, Ivory Coast, Burkina Faso, Ghana, and Nigeria. In other African countries, the mean allele frequency was 0.011.

Conclusions

Our data are consistent with a small number of founder carriers of the amyloidogenic TTR V122I (p.Val142Ile) allele in southern West Africa, with no apparent advantage or disadvantage of an allele carrying newborn reaching adulthood. In U.S. African Americans, the allele represents a significant risk for congestive heart failure late in life. If clinical penetrance is similar in African countries with high allele frequencies, then cardiac amyloidosis could also represent a significant cause of heart disease in the elderly in those populations.

Principles of genetic variations and molecular diseases: applications in hemophilia A

DNA structure alterations are the ultimate source of genetic variations. Without them, evolution would be impossible. While they are essential for DNA diversity, defects in DNA synthesis can lead to numerous genetic diseases. Due to increasingly innovative technologies, our knowledge of the human genome and genetic diseases has grown considerably over the last few years, allowing us to detect another class of variants affecting the chromosomal structure. DNA sequence can be altered in multiple ways: DNA sequence changes by substitution, deletion, or duplication of some nucleotides; chromosomal structure alterations by deletion, duplication, translocation, and inversion, ranging in size from kilobases to mega bases; changes in the cell's genome size. If the alteration is located within a gene and sufficiently deleterious, it can cause genetic disorders. Due to the F8 gene's high rate of new small mutations and its location at the tip of X chromosome, containing high repetitive sequences, a wide variety of genetic variants has been described as the cause of hemophilia A (HA). In addition to the F8 intron 22 repeat inversion, HA can also result from point mutations, other inversions, complex rearrangements, such as duplications or deletions, and transposon insertions causing phenotypes of variable severity characterized by complete or partial deficiency of circulating FVIII. This review aims to present the origins, mechanisms, and consequences of F8 alterations. A sound understanding of the multiple genetic mechanisms responsible for HA is essential to determine the appropriate strategy for molecular diagnosis and detected each type of genetic variant.

Nucleotide composition bias and codon usage trends of gene populations in Mycoplasma capricolum subsp. capricolum and M. Agalactiae

Because of the low GC content of the gene population, amino acids of the two mycoplasmas tend to be encoded by synonymous codons with an A or T end. Compared with the codon usage of ovine, Mycoplasma capricolum and M. agalactiae tend to select optimal codons, which are rare codons in ovine. Due to codon usage pattern caused by genes with key biological functions, the overall codon usage trends represent a certain evolutionary direction in the life cycle of the two mycoplasmas. The overall codon usage trends of a gene population of M. capricolum subsp. capricolum can be obviously separated from other mycoplasmas, and the overall codon usage trends of M. agalactiae are highly similar to those of M. bovis. These results partly indicate the independent evolution of the two mycoplasmas without the limits of the host cell's environment. The GC and AT skews estimate nucleotide composition bias at different positions of nucleotide triplets and the protein consideration caused by the nucleotide composition bias at codon positions 1 and 2 largely take part in synonymous codon usage patterns of the two mycoplasmas. The correlation between the codon adaptation index and codon usage variation indicates that the effect of codon usage on gene expression in M. capricolum subsp. capricolum is opposite to that of M. agalactiae, further suggesting independence of the evolutionary process influencing the overall codon usage trends of gene populations of mycoplasmas.

Epigenetics: A primer for clinicians

With recent advances in cellular biology, we now appreciate that modifications to DNA and histones can have a profound impact on transcription and function, even in the absence of changes to DNA sequence. These modifications, now commonly referred to as "epigenetic" alterations, have changed how we understand cell behavior, reprogramming and differentiation and have provided significant insight into the mechanisms underlying carcinogenesis. Epigenetic alterations, to this point, are largely identified by changes in DNA methylation and hydroxymethylation as well as methylation, acetylation, and phosphorylation of histone tails. These modifications enable significant flexibility in gene expression, rather than just turning genes "ON" or "OFF." Herein we describe the epigenetic landscape in the regulation of gene expression with a particular focus on interrogating DNA methylation in myeloid malignancy.

Assessment of PAX6 alleles in 66 families with aniridia

We report on PAX6 alleles associated with a clinical diagnosis of classical aniridia in 81 affected individuals representing 66 families. Allelic variants expected to affect PAX6 function were identified in 61 families (76 individuals). Ten cases of sporadic aniridia (10 families) had complete (8 cases) or partial (2 cases) deletion of the PAX6 gene. Sequence changes that introduced a premature termination codon into the open reading frame of PAX6 occurred in 47 families (62 individuals). Three individuals with sporadic aniridia (three families) had sequence changes (one deletion, two run-on mutations) expected to result in a C-terminal extension. An intronic deletion of unknown functional significance was detected in one case of sporadic aniridia (one family), but not in unaffected relatives. Within these 61 families, single nucleotide substitutions accounted for 30/61 (49%), indels for 23/61 (38%), and complete deletion of the PAX6 locus for 8/61 (13%). In five cases of sporadic aniridia (five families), no disease-causing mutation in the coding region was detected. In total, 23 unique variants were identified that have not been reported in the Leiden Open Variation Database (LOVD) database. Within the group assessed, 92% had sequence changes expected to reduce PAX6 function, confirming the primacy of PAX6 haploinsufficiency as causal for aniridia.