Excess Heritability Contribution of Alcohol Consumption Variants in the “Missing Heritability” of Type 2 Diabetes Mellitus

We aim to compare the relative heritability contributed by variants of behavior-related environmental phenotypes and elucidate the role of these factors in the conundrum of “missing heritability” of type 2 diabetes. Methods: We used Linkage-Disequilibrium Adjusted Kinships (LDAK) and LDAK-Thin models to calculate the relative heritability of each variant and compare the relative heritability for each phenotype. Biological analysis was carried out for the phenotype whose variants made a significant contribution. Potential hub genes were prioritized based on topological parameters of the protein-–protein interaction network. We included 16 behavior-related phenotypes and 2607 valid variants. In the LDAK model, we found the variants of alcohol consumption and caffeine intake were identified as contributing higher relative heritability than that of the random variants. Compared with the relative expected heritability contributed by the variants associated with type 2 diabetes, the relative expected heritability contributed by the variants associated with these two phenotypes was higher. In the LDAK-Thin model, the relative heritability of variants of 11 phenotypes was statistically higher than random variants. Biological function analysis showed the same distributions among type 2 diabetes and alcohol consumption. We eventually screened out 31 hub genes interacting intensively, four of which were validated and showed the upregulated expression pattern in blood samples seen in type 2 diabetes cases. Conclusion: We found that alcohol consumption contributed higher relative heritability. Hub genes may influence the onset of type 2 diabetes by a mediating effect or a pleiotropic effect. Our results provide new insight to reveal the role of behavior-related factors in the conundrum of “missing heritability” of type 2 diabetes.

Structural Variants Are a Major Component of the Missing Heritability of Autism Spectrum Disorder

The heritability of autism spectrum disorder (ASD), based on 680,000 families and five countries, is estimated to be nearly 80%, yet we lack genetic markers that adequately explain it. It is increasingly clear that genomic structural variants (SVs) are a major component of the missing heritability for many complex phenotypes. Here we use a novel method to identify SVs based on non-Mendelian inheritance patterns in pedigrees using parent-child genotypes from ASD families and demonstrate that the genes that the ASD-specific SVs overlap recapitulate the known molecular biology of ASD including dendritic spinogenesis, axon guidance, and chromatin modification. We further define fine-grained biological pathways that strongly implicate aberrant early development of the cerebellum. Importantly, using these previously excluded variants, we identify the ACMSD gene in the kynurenine pathway as significantly associated with non-verbal cases of ASD and we then use an explainable artificial intelligence approach to define subgroups for future diagnosis and deployment of personalized medicine.

Exploring the missing heritability in subjects with hearing loss, enlarged vestibular aqueducts, and a single or no pathogenic SLC26A4 variant

Pathogenic variants in SLC26A4 have been associated with autosomal recessive hearing loss (arHL) and a unilateral or bilateral enlarged vestibular aqueduct (EVA). SLC26A4 is the second most frequently mutated gene in arHL. Despite the strong genotype–phenotype correlation, a significant part of cases remains genetically unresolved. In this study, we investigated a cohort of 28 Dutch index cases diagnosed with HL in combination with an EVA but without (M0) or with a single (M1) pathogenic variant in SLC26A4. To explore the missing heritability, we first determined the presence of the previously described EVA-associated haplotype (Caucasian EVA (CEVA)), characterized by 12 single nucleotide variants located upstream of SLC26A4. We found this haplotype and a delimited V1-CEVA haplotype to be significantly enriched in our M1 patient cohort (10/16 cases). The CEVA haplotype was also present in two M0 cases (2/12). Short- and long-read whole genome sequencing and optical genome mapping could not prioritize any of the variants present within the CEVA haplotype as the likely pathogenic defect. Short-read whole-genome sequencing of the six M1 cases without this haplotype and the two M0/CEVA cases only revealed previously overlooked or misinterpreted splice-altering SLC26A4 variants in two cases, who are now genetically explained. No deep-intronic or structural variants were identified in any of the M1 subjects. With this study, we have provided important insights that will pave the way for elucidating the missing heritability in M0 and M1 SLC26A4 cases. For pinpointing the pathogenic effect of the CEVA haplotype, additional analyses are required addressing defect(s) at the RNA, protein, or epigenetic level.

Decoding ‘Unnecessary Complexity’: A Law of Complexity and a Concept of Hidden Variation Behind “Missing Heritability” in Precision Medicine

The high hopes for the Human Genome Project and personalized medicine were not met because the relationship between genotypes and phenotypes turned out to be more complex than expected. In a previous study we laid the foundation of a theory of complexity and showed that because of the blind nature of evolution, and molecular and historical contingency, cells have accumulated unnecessary complexity, complexity beyond what is necessary and sufficient to describe an organism. Here we provide empirical evidence and show that unnecessary complexity has become integrated into the genome in the form of redundancy and is relevant to molecular evolution of phenotypic complexity. Unnecessary complexity creates uncertainty between molecular and phenotypic complexity, such that phenotypic complexity (CP) is higher than molecular complexity (CM), which is higher than DNA complexity (CD). The qualitative inequality in complexity is based on the following hierarchy: CP > CM > CD. This law-like relationship holds true for all complex traits, including complex diseases. We present a hypothesis of two types of variation, namely open and closed (hidden) systems, show that hidden variation provides a hitherto undiscovered “third source” of phenotypic variation, beside genotype and environment, and argue that “missing heritability” for some complex diseases is likely to be a case of “diluted heritability”. There is a need for radically new ways of thinking about the principles of genotype–phenotype relationship. Understanding how cells use hidden, pathway variation to respond to stress can shed light on why two individuals who share the same risk factors may not develop the same disease, or how cancer cells escape death.

The contribution of rare whole genome sequencing variants to plasma protein levels and to the missing heritability

Despite the success in identifying effects of common genetic variants, using genome-wide association studies (GWAS), much of the genetic contribution to complex traits remains unexplained. Here, we analysed high coverage whole-genome sequencing (WGS) data, to evaluate the contribution of rare genetic variants to 414 plasma proteins. The frequency distribution of genetic variants was skewed towards the rare spectrum, and damaging variants were more often rare. However, only 2.24% of the heritability was estimated to be explained by rare variants. A gene-based approach, developed to also capture the effect of rare variants, identified associations for 249 of the proteins, which was 25% more as compared to a GWAS. Out of those, 24 associations were driven by rare variants, clearly highlighting the capacity of aggregated tests and WGS data. We conclude that, while many rare variants have considerable phenotypic effects, their contribution to the missing heritability is limited by their low frequencies.

Exploring the Missing Heritability in Subjects With Hearing Loss, Enlarged Vestibular Aqueducts, and A Single or No Pathogenic SLC26A4 Variant

Pathogenic variants in SLC26A4 have been associated with autosomal recessive hearing loss (arHL) and a unilateral or bilateral enlarged vestibular aqueduct (EVA). SLC26A4 is the second most frequently mutated gene in arHL. Despite the strong genotype-phenotype correlation, a significant part of SLC26A4 cases remains genetically unresolved. In this study, we investigated a cohort of 28 Dutch index cases diagnosed with HL in combination with an EVA but without (M0) or with a single (M1) pathogenic variant in SLC26A4. To explore the missing heritability, short- and long-read whole genome sequencing and optical genome mapping were performed. We found a previously described EVA-associated haplotype (Caucasian EVA (CEVA)) to be significantly enriched in our M1 patient cohort. The haplotype was also present in two M0 cases. Despite extensive genetic analyses, we were not able to prioritize any of the variants present within the haplotype as the likely pathogenic defect, and therefore additional analyses addressing the defect(s) at the RNA, protein, or epigenetic level are required. Whole genome sequencing also revealed splice-altering SLC26A4 variants in two M1 cases, which are now genetically explained, but no deep-intronic or copy number variants. With these findings, we have provided important insights that will pave the way for elucidating the missing heritability in M0 and M1 SLC26A4 cases.