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.

Widespread somatic L1 retrotransposition in normal colorectal epithelium.

Throughout an individual's lifetime, genomic alterations accumulate in somatic cells1-11. However, the mutational landscape induced by retrotransposition of long interspersed nuclear element-1 (L1), a widespread mobile element in the human genome12-14, is poorly understood in normal cells. Here we explored the whole-genome sequences of 899 single-cell clones established from three different cell types collected from 28 individuals. We identified 1,708 somatic L1 retrotransposition events that were enriched in colorectal epithelium and showed a positive relationship with age. Fingerprinting of source elements showed 34 retrotransposition-competent L1s. Multidimensional analysis demonstrated that (1) somatic L1 retrotranspositions occur from early embryogenesis at a substantial rate, (2) epigenetic on/off of a source element is preferentially determined in the early organogenesis stage, (3) retrotransposition-competent L1s with a lower population allele frequency have higher retrotransposition activity and (4) only a small fraction of L1 transcripts in the cytoplasm are finally retrotransposed in somatic cells. Analysis of matched cancers further suggested that somatic L1 retrotransposition rate is substantially increased during colorectal tumourigenesis. In summary, this study illustrates L1 retrotransposition-induced somatic mosaicism in normal cells and provides insights into the genomic and epigenomic regulation of transposable elements over the human lifetime.

Clonal dynamics in early human embryogenesis inferred from somatic mutation.

Cellular dynamics and fate decision in early human embryogenesis remain largely unknown owing to the challenges of performing studies in human embryos1. Here, we explored whole-genomes of 334 single-cell colonies and targeted deep sequences of 379 bulk tissues obtained from various anatomical locations of seven recently deceased adult human donors. Using somatic mutations as an intrinsic barcode, we reconstructed early cellular phylogenies that demonstrate (1) an endogenous mutational rate that is higher in the first cell division but decreases to approximately one per cell per cell division later in life; (2) universal unequal contribution of early cells to embryo proper, resulting from early cellular bottlenecks that stochastically set aside epiblast cells within the embryo; (3) examples of varying degrees of early clonal imbalances between tissues on the left and right sides of the body, different germ layers and specific anatomical parts and organs; (4) emergence of a few ancestral cells that will substantially contribute to adult cell pools in blood and liver; and (5) presence of mitochondrial DNA heteroplasmy in the fertilized egg. Our approach also provides insights into the age-related mutational processes and loss of sex chromosomes in normal somatic cells. In sum, this study provides a foundation for future studies to complete cellular phylogenies in human embryogenesis.

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.

The obesity-risk variant of FTO is inversely related with the So-Eum constitutional type: genome-wide association and replication analyses.

BACKGROUND: Body constitutional types described in the traditional Korean medicine system, Sasang constitutional medicine, are heritable, as has been revealed by twin and family studies. Thus, individuals with the same constitution type usually have similar pathophysiological and psychological traits. In several recent genome-wide association (GWA) analyses performed to identify constitution-associated variants, the association signals were not replicated due to small sample size and dissimilar, non-objective methods for classification of the constitutional types. METHODS: We conducted GWA analysis and followed replication analysis in two large populations (5,490 subjects: 3,810 subjects at discovery stage and 1,680 subjects at replication stage) to identify the replicable constitution-associated variants, wherein subjects with the highest tertile of constitution probability values versus the reference with the lowest tertile of the values obtained from a recently developed constitution analysis tool were compared. RESULTS: We found that the obesity-risk variant in intron 1 of the fat mass and obesity-associated (FTO) gene was replicably inversely associated with the So-Eum (SE) type, characterized by reduced appetite, slim body, and cautious personality (rs7193144 in combined samples: odds ratio = 0.729, p = 1.47 × 10(-7)), and substantial association signal remained after controlling for body mass index (BMI). In contrast, the association of the variant with the Tae-Eum type, characterized by high body mass, disappeared after controlling BMI. CONCLUSIONS: In summary, the obesity-risk variant in FTO intron 1 was inversely associated with the SE type, independent of BMI, which corresponded well with the characteristics of the SE type, such as the lowest body mass and lowest susceptibility to metabolic disorders among the constitutional types. Therefore, the obesity-risk variant of FTO associated with body mass increase might be involved in the determination of body constitution type.

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.

Association of apolipoprotein A5 gene -1131T>C polymorphism with the risk of metabolic syndrome in Korean subjects.

We assessed the associations between the APOA5 -1131T>C polymorphism and lipid parameters and other risk factors of the metabolic syndrome in Korean subjects. A total of 2,901 participants from 20 oriental medical hospitals in Korea were enrolled between 2006 and 2011. According to the modified National Cholesterol Education Program Adult Treatment Panel III definitions, subjects were classified into the metabolic syndrome group and control group. The APOA5 -1131T>C genotype was significantly associated with serum high-density lipoprotein cholesterol levels (effect = - 1.700 mg/dL, P=6.550-E07) in the total study population after adjustment for differences in age and gender. The association of the APOA5 -1131T>C genotype with serum log-transformed triglyceride was also significant in an additive genetic model (effect = 0.056 mg/dL, P=2.286E-19). After adjustment for age and gender, we determined that the odds ratio for the occurrence of the metabolic syndrome was 1.322 for C-allele carriers in the additive model (95% CI = [1.165 - 1.501], P=1.48E-05). In the current study, we demonstrated that the APOA5 -1131T>C polymorphism is associated with the metabolic syndrome because of its remarkable effect on serum triglyceride levels in Korean subjects.