Parallel shift of DNA methylation and gene expression toward the mean in mouse spleen with aging.

Age-associated DNA-methylation drift (AMD) manifests itself in two ways in mammals: global decrease (hypomethylation) and local increase of DNA methylation (hypermethylation). To comprehend the principle behind this bidirectional AMD, we studied methylation states of spatially clustered CpG dinucleotides in mouse splenic DNA using reduced-representation-bisulfite-sequencing (RRBS). The mean methylation levels of whole CpGs declined with age. Promoter-resident CpGs, generally weakly methylated (<5%) in young mice, became hypermethylated in old mice, whereas CpGs in gene-body and intergenic regions, initially moderately (~33%) and extensively (>80%) methylated, respectively, were hypomethylated in the old. Chromosome-wise analysis of methylation revealed that inter-individual heterogeneities increase with age. The density of nearby CpGs was used to classify individual CpGs, which found hypermethylation in CpG-rich regions and hypomethylation in CpG-poor regions. When genomic regions were grouped by methylation level, high-methylation regions tended to become hypomethylated whereas low-methylation regions tended to become hypermethylated, regardless of genomic structure/function. Data analysis revealed that while methylation level and CpG density were interdependent, methylation level was a better predictor of the AMD pattern representing a shift toward the mean. Further analysis of gene-expression data showed a decrease in the expression of highly-expressed genes and an increase in the expression of lowly-expressed genes with age. This shift towards the mean in gene-expression changes was correlated with that of methylation changes, indicating a potential link between the two age-associated changes. Our findings suggest that age-associated hyper- and hypomethylation events are stochastic and attributed to malfunctioning intrinsic mechanisms for methylation maintenance in low- and high-methylation regions, respectively.

K120R mutation inactivates p53 by creating an aberrant splice site leading to nonsense-mediated mRNA decay.

The point mutation that substitutes lysine with arginine at position 120 of human p53 has been characterized as a missense mutation. The K120R mutation renders the p53 protein disabled for acetylation and, as a result, defective for apoptotic function, which provides a mechanistic link between the missense mutation and tumorigenesis. However, we noticed the failures of tumorigenesis in mice with the mutation, and of the related studies to notice that it has arbitrarily reflected in amino acid change through a sequence modification (AGA) of the original tumor mutation (AGG) by codon degeneracy. Unlike this modified version, we also discovered a novel splicing site the original mutation, TP53 c.359A>G, may induce. Using a human induced pluripotent stem cell line that was engineered to be homozygous for the original mutation, we here identified that the accidental splicing site generates a defective transcript variant with a frame-shifted premature termination codon which is subjected to nonsense-mediated mRNA decay. The authentic splicing still occurs but in extremely low amounts. Taken together, this mutation causes depletion of cellular p53 via defective mRNA, suggesting a new link to tumorigenesis.

Generation of gene-corrected iPSC line from Parkinson's disease patient iPSC line with alpha-SNCA A53T mutation.

Parkinson's disease (PD) is the second most common age-related neurodegenerative disorder. PD can result from a mutation of alpha-synuclein (α-SNCA), such as α-SNCA A53T. Using episomal vectors, induced pluripotent stem cells (iPSCs) were generated from skin fibroblasts with the α-SNCA A53T mutation. A huge bacterial artificial chromosome (BAC) harboring the normal α-SNCA gene successfully corrected the α-SNCA A53T-mutant iPSCs. Melting curve analysis for allelic composition indicated that the BAC DNA was precisely targeted to the α-SNCA A53T mutation allele, without random integration. The corrected PD-iPSCs displayed the normal karyotype and pluripotency, with the capability to differentiate to any cell type.

Age-associated bimodal transcriptional drift reduces intergenic disparities in transcription.

This study addressed the question of how well the quantitative transcriptome structure established in early life is maintained and how consistently it appears with increasing age, and if there is age-associated alteration of gene expression (A3GE), how much influence the Huntington's disease (HD) genotype exerts on it. We examined 285 exonic sequences of 175 genes using targeted PCR sequencing in skeletal muscle, brain, and splenic CD4+ T cells of wild-type and HD mice. In contrast to the muscle and brain, T cells exhibited large A3GE, suggesting a strong contribution to functional decline of the organism. This A3GE was markedly intensified in age-matched HD T cells, which exhibited accelerated aging as determined by reduced telomere length. Regression analysis suggested that gene expression levels change at a rate of approximately 3% per month with age. We found a bimodal relationship in A3GE in T cells in that weakly expressed genes in young mice were increasingly transcribed in older animals whereas highly expressed genes in the young were decreasingly expressed with age. This bimodal transcriptional drift in the T cell transcriptome data causes the differences in transcription rate between genes to progressively reduce with age.

Age-associated chromatin relaxation is enhanced in Huntington's disease mice.

Expansion of polyglutamine stretch in the huntingtin (HTT) protein is a major cause of Huntington's disease (HD). The polyglutamine part in HTT interacts with various proteins implicated in epigenetic regulation of genes, suggesting that mutant HTT may disturb the integrity of the epigenetic system. Here, we used a PCRseq-based method to examine expression profile of 395 exonic segments from 260 "epi-driver" genes in splenic T lymphocytes from aged HD mice. We identified 67 exonic segments differentially expressed between young and aged HD mice, most of them upregulated in the aged. Polycomb-repressive complex (PRC)-regulated genes (PRGs) were markedly upregulated in aged HD mice, consistent with downregulation of PRC genes. Epi-driver gene categories of lysine-methylation, lysine-demethylation, arginine-methylation, and PRG showed differential age-associated changes between HD and control. Analyzing the pattern of change in epi-driver gene expressions hinted at an enhanced shift in HD chromatin to a more accessible state with age, which was experimentally demonstrated by DNase-I-hypersensitivity sequencing showing increased chromatin accessibility in HD cells compared to control. We suggest the global change can potentially relieve chromatin-induced repression of many genes, and the unintended expressions of some detrimental proteins could alter T cell function to a greater degree in aged HD mice.

Normalization of human RNA-seq experiments using chimpanzee RNA as a spike-in standard.

Normalization of human RNA-seq experiments employing chimpanzee RNA as a spike-in standard is reported. Human and chimpanzee RNAs exhibit single nucleotide variations (SNVs) in average 210-bp intervals. Spike-in chimpanzee RNA would behave the same as the human counterparts during the whole NGS procedures owing to the high sequence similarity. After discrimination of species origins of the NGS reads based on SNVs, the chimpanzee reads were used to read-by-read normalize biases and variations of human reads. By this approach, as many as 10,119 transcripts were simultaneously normalized for the entire NGS procedures leading to accurate and reproducible quantification of differential gene expression. In addition, incomparable data sets from different in-process degradations or from different library preparation methods were made well comparable by the normalization. Based on these results, we expect that the normalization approaches using near neighbor genomes as internal standards could be employed as a standard protocol, which will improve both accuracy and comparability of NGS results across different sample batches, laboratories and NGS platforms.

Concurrent progress of reprogramming and gene correction to overcome therapeutic limitation of mutant ALK2-iPSC.

Fibrodysplasia ossificans progressiva (FOP) syndrome is caused by mutation of the gene ACVR1, encoding a constitutive active bone morphogenetic protein type I receptor (also called ALK2) to induce heterotopic ossification in the patient. To genetically correct it, we attempted to generate the mutant ALK2-iPSCs (mALK2-iPSCs) from FOP-human dermal fibroblasts. However, the mALK2 leads to inhibitory pluripotency maintenance, or impaired clonogenic potential after single-cell dissociation as an inevitable step, which applies gene-correction tools to induced pluripotent stem cells (iPSCs). Thus, current iPSC-based gene therapy approach reveals a limitation that is not readily applicable to iPSCs with ALK2 mutation. Here we developed a simplified one-step procedure by simultaneously introducing reprogramming and gene-editing components into human fibroblasts derived from patient with FOP syndrome, and genetically treated it. The mixtures of reprogramming and gene-editing components are composed of reprogramming episomal vectors, CRISPR/Cas9-expressing vectors and single-stranded oligodeoxynucleotide harboring normal base to correct ALK2 c.617G>A. The one-step-mediated ALK2 gene-corrected iPSCs restored global gene expression pattern, as well as mineralization to the extent of normal iPSCs. This procedure not only helps save time, labor and costs but also opens up a new paradigm that is beyond the current application of gene-editing methodologies, which is hampered by inhibitory pluripotency-maintenance requirements, or vulnerability of single-cell-dissociated iPSCs.

H19ICR mediated transcriptional silencing does not require target promoter methylation.

Transcription of the reciprocally imprinted genes Insulin-like growth factor 2 (Igf2) and H19 is orchestrated by the 2.4-kb H19 Imprinting Control Region (H19ICR) located upstream of H19. Three known functions are associated with the H19ICR: (1) it is a germline differentially methylated region, (2) it is a transcriptional insulator, and (3) it is a transcriptional silencer. The molecular mechanisms of the DMR and insulator functions have been well characterized but the basis for the ICR's silencer function is less well understood. In order to study the role the H19ICR intrinsically plays in gene silencing, we transferred the 2.4-kb H19ICR to a heterologous non-imprinted location on chromosome 5, upstream of the alpha fetoprotein (Afp) promoter. Independent of its orientation, the 2.4-kb H19ICR silences transcription from the paternal Afp promoter. Thus silencing is a function intrinsic to this DNA element. Further, ICR mediated silencing is a developmental process that, unexpectedly, does not occur through DNA methylation at the target promoter.

SiNG-PCRseq: Accurate inter-sequence quantification achieved by spiking-in a neighbor genome for competitive PCR amplicon sequencing.

Despite the recent technological advances in DNA quantitation by sequencing, accurate delineation of the quantitative relationship among different DNA sequences is yet to be elaborated due to difficulties in correcting the sequence-specific quantitation biases. We here developed a novel DNA quantitation method via spiking-in a neighbor genome for competitive PCR amplicon sequencing (SiNG-PCRseq). This method utilizes genome-wide chemically equivalent but easily discriminable homologous sequences with a known copy arrangement in the neighbor genome. By comparing the amounts of selected human DNA sequences simultaneously to those of matched sequences in the orangutan genome, we could accurately draw the quantitative relationships for those sequences in the human genome (root-mean-square deviations <0.05). Technical replications of cDNA quantitation performed using different reagents at different time points also resulted in excellent correlations (R(2) > 0.95). The cDNA quantitation using SiNG-PCRseq was highly concordant with the RNA-seq-derived version in inter-sample comparisons (R(2) = 0.88), but relatively discordant in inter-sequence quantitation (R(2) < 0.44), indicating considerable level of sequence-dependent quantitative biases in RNA-seq. Considering the measurement structure explicitly relating the amount of different sequences within a sample, SiNG-PCRseq will facilitate sharing and comparing the quantitation data generated under different spatio-temporal settings.

Quantifying difference in gene expression profile between bovine blastocysts derived by in vitro fertilization and somatic cell nuclear transfer.

Epigenetic reprogramming intensely occurs in somatic-cell nuclear transfer (SCNT) embryos, which highlights the importance of proper expressions of reprogramming-related genes in SCNT embryos. We here assessed gene expression profiles (GEPs) difference between bovine blastocyst groups derived by in-vitro fertilization (IVF) or SCNT; in SCNT, cumulus cells and ear skin fibroblasts were used for cSCNT and fSCNT blastocysts, respectively. We obtained GEPs of 15 reprogramming-related genes in single blastocysts using multiplex PCR and found a broad range of variations in their GEPs. Weighted root-mean-square deviation (wRMSD) analysis, which calculates the deviation of SCNT blastocysts' GEPs from IVF blastocysts' mean GEP, found a significant difference between IVF and fSCNT and between cSCNT and fSCNT blastocysts (p < 0.001) but not between IVF and cSCNT. Since the fibroblasts' GEP was more distant from the IVF blastocysts' than the cumulus cells', it might partly explain the less similarity of fSCNT blastocysts' GEPs to the IVF's mean GEP. Our wRMSD method succeeds in expressing in figures how different two comparable embryo groups of different derivations are in GEP, which would be useful to select a better embryo derivation protocol among the candidates prior to field applications.

Assessment of difference in gene expression profile between embryos of different derivations.

Researchers have exerted sustained efforts to improve the viability of somatic cell nuclear transfer (SCNT) embryos, testing their experimental designs and probing the resultant embryos. However, the lack of a reliable method to estimate the efficacy of these experimental attempts is a chief hindrance to tackling the low-viability problem in SCNT. Here, we introduce a procedure that assesses the degree of difference in gene expression profiles (GEPs) of blastocysts from each other as a representative control of good quality. We first adapted a multiplex reverse transcription-polymerase chain reaction strategy to obtain GEPs for 15 reprogramming-related genes from single mouse blastocysts. GEPs of individual blastocysts displayed a broad range of variations, the extent of which was calculated using a weighted root mean square deviation (wRMSD). wRMSD-based quantitation of GEP difference (qGEP) found that GEP difference between in vivo-derived blastocysts (in vivo) and SCNT blastocysts was greater than the difference between in vivo blastocysts and in vitro-produced (IVP) blastocysts, demonstrating that the SCNT group was more distantly related to the in vivo group than the IVP group. Our qGEP approach for grading individual blastocysts would be useful for selecting a better protocol to derive embryos of better quality prior to field applications.

Identification of cis-regulatory variations in the IL6R gene through the inheritance assessment of allelic transcription.

BACKGROUND: The level of circulating interleukin-6 receptor in human blood varies depending on the genetic and/or physiological causes, and has been implicated in the development of chronic inflammatory diseases. METHOD: The cis-regulatory effects of genetic variations on the transcription of interleukin-6 receptor gene, IL6R,were studied by assessing allelic transcriptions in the immortalized lymphocytes derived from unrelated and familial samples. RESULTS: The assays for allelic transcription in the cells from unrelated subjects demonstrated an extensive and variable range of allelic transcriptional imbalances, suggesting an operation of multiple cis-regulations with varying degrees on the locus. Analysis of the familial samples illustrated the Mendelian inheritance of allelic transcriptions,enabling us to assign each haplotype allele into one of the 3 transcriptional strengths. A comparison of the allele structures based on the transcriptional attributes highlighted 2 SNP variations, rs952146 and rs4845617, as being associated with higher allelic transcription. Consistently, lymphocytes that were homozygous for the 2SNPs exhibited differences in their transcript levels depending on the haplotypes. CONCLUSION: Inheritance assessment of allelic transcription of IL6R identified 2 SNPs that are associated with transcriptional variation in cis. GENERAL SIGNIFICANCE: Our results not only demonstrate genetic variations that are associated with IL6R transcription in cis but also demonstrate an effective genetic approach for isolating cis-regulatory variations.

Accurate measurement of the relative abundance of different DNA species in complex DNA mixtures.

A molecular tool that can compare the abundances of different DNA sequences is necessary for comparing intergenic or interspecific gene expression. We devised and verified such a tool using a quantitative competitive polymerase chain reaction approach. For this approach, we adapted a competitor array, an artificially made plasmid DNA in which all the competitor templates for the target DNAs are arranged with a defined ratio, and melting analysis for allele quantitation for accurate quantitation of the fractional ratios of competitively amplified DNAs. Assays on two sets of DNA mixtures with explicitly known compositional structures of the test sequences were performed. The resultant average relative errors of 0.059 and 0.021 emphasize the highly accurate nature of this method. Furthermore, the method's capability of obtaining biological data is demonstrated by the fact that it can illustrate the tissue-specific quantitative expression signatures of the three housekeeping genes G6pdx, Ubc, and Rps27 by using the forms of the relative abundances of their transcripts, and the differential preferences of Igf2 enhancers for each of the multiple Igf2 promoters for the transcription.

Inactivation of the MSLNL gene encoding mesothelin-like protein during African great ape evolution.

Loss of gene function is implicated in the emergence of novel phenotypes during organism evolution. Here, we report the inactivation of the MSLNL gene encoding mesothelin-like protein in African great ape evolution. Human MSLNL has a nonsense mutation in exon 10 and two polymorphic mutations: a frameshift in exon 3 and a nonsense codon in exon 8. The gorilla gene also shows multiple deleterious mutations, including a premature stop codon, a deletion, and a splice site mutation. Molecular evolutionary analysis indicated relaxed selection pressure on MSLNL in African great ape lineages, which suggested that MSLNL might have become inactivated before the divergence of human, chimpanzee and gorilla. The mouse Mslnl gene is highly expressed in olfactory epithelium and moderately expressed in several other tissues. We propose that the loss of MSLNL may be associated with the evolution of the olfactory system in African great apes including human.

SNP genotyping through the melting analysis of unlabelled oligonucleotide applied on dilute PCR amplicon.

Adaptation of DNA melting analysis for polymorphic single nucleotides (SNPs) genotyping using an unlabeled oligonucleotide probe for polymorphic DNAs under the presence of fluorescent DNA binding dye necessitates a reaction condition where the probe efficiently associates with a target strand that is PCR amplified. We present experimental evidence that application of an unlabeled probe to a dilute PCR amplicon provides a condition such that the fluorescent signals gained subsequently by probe melting are sufficient to discriminate allelic identities. This approach is best exploited by adapting the multiplexing PCR technique in order to cover multiple SNPs for given samples. 3'-end modification of the probe is unnecessary as the amplicon dilution step provides a way of inactivating the polymerase through divalent cation chelation. With the use of low-cost reagents and ordinary laboratory equipment, this method offers a rapid, simple and cost-efficient way of SNP genotyping.

Genetic Effects of FTO and MC4R Polymorphisms on Body Mass in Constitutional Types.

Sasang constitutional medicine (SCM), a Korean tailored medicine, categorizes human beings into four types through states of physiological imbalances and responsiveness to herbal medicine. One SCM type susceptible to obesity seems sensitive to energy intake due to an imbalance toward preserving energy. Common variants of fat mass and obesity associated (FTO) and melanocortin 4 receptor (MC4R) genes have been associated with increased body mass index (BMI) by affecting energy intake. Here, we statistically examined the association of FTO and MC4R polymorphisms with BMI in two populations with 1370 Koreans before and after SCM typing, and with the lowering of BMI in 538 individuals who underwent a 1-month lifestyle intervention. The increased BMI replicated the association with FTO haplotypes (effect size ≃ 1.1 kg/m(2)) and MC4R variants (effect size ≃ 0.64 kg/m(2)). After the lifestyle intervention, the carriers of the haplotype represented by the minor allele of rs1075440 had a tendency to lose more waist-to-hip ratio (0.76%) than non-carriers. The constitutional discrepancy for the accumulation of body mass by the effects of FTO and/or MC4R variants seemed to reflect the physique differences shown in each group of SCM constitutional types. In conclusion, FTO and MC4R polymorphisms appear to play an important role in weight gain, while only FTO variants play a role in weight loss after lifestyle intervention. Different trends were observed among individuals of SCM types, especially for weight gain. Therefore, classification of individuals based on physiological imbalance would offer a good genetic stratification system in assessing the effects of obesity genes.

Reciprocal regulation of gene expression by Ephedra herba in mouse brain.

Although Ephedra herba has been used to treat a variety of diseases, the molecular mechanism has not been defined clearly. To identify the molecular events, a microarray analysis was conducted using the brain tissue of mice orally administered aqueous extracts of Ephedra herba, in which the expression level of many genes was altered depending on the duration of treatment time. In particular, temporal changes in the expressions of genes can be classified into two major reciprocal patterns in which 175 and 282 genes were included in pattern 1 and 2, respectively. Classification based on gene ontology revealed different functional implications of genes between temporal pattern 1 and 2. Through promoter analysis of temporally co-expressed genes, putative co-regulated genes could be segregated based on the similarity of the transcription factor binding sites (TFBS). Interestingly, two temporal patterns were associated with different sub-classes of TFBS similarity, which implies the presence of dual regulatory mechanisms. Finally, through interaction analysis, core node genes that could interact with many other nodes were identified. By integrating the interaction information with promoter clustering, temporally co-regulated key components could be selected in mouse brain after administration of Ephedra herba.

Genetic approach to elucidation of sasang constitutional medicine.

Sasang Constitutional Medicine (SCM) offers a medical principle that classifies humans into four constitution groups and guides their treatment with constitution-matched medical assistance. The principle of this traditional medicine, although requires significant scientific support, appears to suggest a genetic influence on constitution type. The relative frequency of constitution types in a population, for instance, has remained relatively constant since Jema Lee first described them from his observations. In addition, the body compartment concept of SCM appears to be related to the anterio-posterior patterning of the embryonic gut and associated internal organs. This study describes the attributes of the constitution concept of SCM that can be interpreted in the language of genetics and current approaches to identity the genetic factors that make up the constitution. These efforts should make it possible to interpret the principle of this traditional medicine scientifically. Considering the recent trend in medicine that pursues individualized or tailored medical offerings, once SCM is proven to be explainable with scientific evidence, it will be able to contribute to and take a place in the rapidly evolving medicine environment.

Relative quantitation of restriction fragment length polymorphic DNAs via DNA melting analysis provides an effective way to determine allele frequencies.

To accurately and precisely estimate the allele frequencies in DNA pools for a cost-effective approach to correlate genetic variations to phenotypic traits, we exploited differential melting kinetics between restriction fragment length polymorphic DNAs. The allele frequencies of two SNPs in a series of DNA mixtures with known allelic compositions of the SNPs were determined by analyzing the meltings of restricted PCR amplicons, yielding a result with a root mean square error (RMSE) of 0.014 relative to the expected values and a standard deviation (SD) of 0.018 from triplicate measurements. This method was then applied in the measurement of genotype frequencies in DNA pools in which varying numbers of genomic DNAs were intermingled while maintaining uniform quantitative contribution. Analyses of 10 SNPs demonstrated the feasibility of this method in an economical and highly accurate manner as the results yielded an RMSE value of 0.027 and a SD of 0.019.

Inactivation of MOXD2 and S100A15A by exon deletion during human evolution.

We devised a bioinformatics method for systematic identification of putative human-specific exon-deletion mutations that occurred after the divergence of human and chimpanzee and experimentally verified 2 of the predicted mutations in MOXD2 and S100A15A genes. MOXD2 gene encodes a monooxygenase that is highly conserved in mammals and is mostly expressed in the olfactory epithelium in mouse. The presence of a deletion of the last 2 exons and a polymorphic nonsense mutation in exon 6 suggests that MOXD2 gene is inactive in humans. S100A15A is a member of the S100 family of calcium-binding proteins, the mouse ortholog of which is expressed during epidermal maturation. Human S100A15A gene is likely to be inactive because the start codon-bearing exon is deleted in human. We propose that modification or inactivation of MOXD2 and S100A15A genes have contributed to the loss of certain smell sense in humans and to the development of human skin.