Inter-chromosomal level of genome organization and longevity-related phenotypes in humans

Evidence of epistasis in regions of long-range linkage disequilibrium across five complex diseases in the UK Biobank and eMERGE datasets

Leveraging linkage disequilibrium (LD) patterns as representative of population substructure enables the discovery of additive association signals in genome-wide association studies (GWASs). Standard GWASs are well-powered to interrogate additive models; however, new approaches are required for invesigating other modes of inheritance such as dominance and epistasis. Epistasis, or non-additive interaction between genes, exists across the genome but often goes undetected because of a lack of statistical power. Furthermore, the adoption of LD pruning as customary in standard GWASs excludes detection of sites that are in LD but might underlie the genetic architecture of complex traits. We hypothesize that uncovering long-range interactions between loci with strong LD due to epistatic selection can elucidate genetic mechanisms underlying common diseases. To investigate this hypothesis, we tested for associations between 23 common diseases and 5,625,845 epistatic SNP-SNP pairs (determined by Ohta's D statistics) in long-range LD (>0.25 cM). Across five disease phenotypes, we identified one significant and four near-significant associations that replicated in two large genotype-phenotype datasets (UK Biobank and eMERGE). The genes that were most likely involved in the replicated associations were (1) members of highly conserved gene families with complex roles in multiple pathways, (2) essential genes, and/or (3) genes that were associated in the literature with complex traits that display variable expressivity. These results support the highly pleiotropic and conserved nature of variants in long-range LD under epistatic selection. Our work supports the hypothesis that epistatic interactions regulate diverse clinical mechanisms and might especially be driving factors in conditions with a wide range of phenotypic outcomes.

Hierarchical dinucleotide distribution in genome along evolution and its effect on chromatin packing

It describes how hierarchical CpG distribution on the genomes of species change in the evolution and the correlation with the chromatin structure.

Dinucleotide densities and their distribution patterns vary significantly among species. Previous studies revealed that CpG is susceptible to methylation, enriched at topologically associating domain boundaries and its distribution along the genome correlates with chromatin compartmentalization. However, the multi-scale organizations of CpG in the linear genome, their role in chromatin organization, and how they change along the evolution are only partially understood. By comparing the CpG distribution at different genomic length scales, we quantify the difference between the CpG distributions of different species and evaluate how the hierarchical uneven CpG distribution appears in evolution. The clustering of species based on the CpG distribution is consistent with the phylogenetic tree. Interestingly, we found the CpG distribution and chromatin structure to be correlated in many different length scales, especially for mammals and avians, consistent with the mosaic CpG distribution in the genomes of these species.

Coordinated Action of Biological Processes during Embryogenesis Can Cause Genome-Wide Linkage Disequilibrium in the Human Genome and Influence Age-Related Phenotypes

A role of non-Mendelian inheritance in genetics of complex, age-related traits is becoming increasingly recognized. Recently, we reported on two inheritable clusters of SNPs in extensive genome-wide linkage disequilibrium (LD) in the Framingham Heart Study (FHS), which were associated with the phenotype of premature death. Here we address biologically-related properties of these two clusters. These clusters have been unlikely selected randomly because they are functionally and structurally different from matched sets of randomly selected SNPs. For example, SNPs in LD from each cluster are highly significantly enriched in genes (p=7.1×10-22 and p=5.8×10-18), in general, and in short genes (p=1.4×10-47 and p=4.6×10-7), in particular. Mapping of SNPs in LD to genes resulted in two, partly overlapping, networks of 1764 and 4806 genes. Both these networks were gene enriched in developmental processes and in biological processes tightly linked with development including biological adhesion, cellular component organization, locomotion, localization, signaling, (p<10-4, q<10-4 for each category). Thorough analysis suggests connections of these genetic networks with different stages of embryogenesis and highlights biological interlink of specific processes enriched for genes from these networks. The results suggest that coordinated action of biological processes during embryogenesis may generate genome-wide networks of genetic variants, which may influence complex age-related phenotypes characterizing health span and lifespan.