Parent-of-origin effects on quantitative phenotypes in a large Hutterite pedigree

Asymptotically exact fit for linear mixed model in genetic association studies

The linear mixed model (LMM) has become a standard in genetic association studies to account for population stratification and relatedness in the samples to reduce false positives. Much recent progresses in LMM focused on approximate computations. Exact methods remained computationally demanding and without theoretical assurance. The computation is particularly challenging for multiomics studies where tens of thousands of phenotypes are tested for association with millions of genetic markers. We present IDUL and IDUL† that use iterative dispersion updates to fit LMMs, where IDUL† is a modified version of IDUL that guarantees likelihood increase between updates. Practically, IDUL and IDUL† produced identical results, both are markedly more efficient than the state-of-the-art Newton-Raphson method, and in particular, both are highly efficient for additional phenotypes, making them ideal to study genetic determinants of multiomics phenotypes. Theoretically, the LMM likelihood is asymptotically unimodal, and therefore the gradient ascent algorithm IDUL† is asymptotically exact. A software package implementing IDUL and IDUL† for genetic association studies is freely available at https://github.com/haplotype/IDUL.

Intergenerational transmission of the structure of the auditory cortex and reading skills

High-level cognitive skill development relies on genetic and environmental factors, tied to brain structure and function. Inter-individual variability in language and music skills has been repeatedly associated with the structure of the auditory cortex: the shape, size and asymmetry of the transverse temporal gyrus (TTG) or gyri (TTGs). TTG is highly variable in shape and size, some individuals having one single gyrus (also referred to as Heschl's gyrus, HG) while others presenting duplications (with a common stem or fully separated) or higher-order multiplications of TTG. Both genetic and environmental influences on children's cognition, behavior, and brain can to some to degree be traced back to familial and parental factors. In the current study, using a unique MRI dataset of parents and children (135 individuals from 37 families), we ask whether the anatomy of the auditory cortex is related to reading skills, and whether there are intergenerational effects on TTG(s) anatomy. For this, we performed detailed, automatic segmentations of HG and of additional TTG(s), when present, extracting volume, surface area, thickness and shape of the gyri. We tested for relationships between these and reading skill, and assessed their degree of familial similarity and intergenerational transmission effects. We found that volume and area of all identified left TTG(s) combined was positively related to reading scores, both in children and adults. With respect to intergenerational similarities in the structure of the auditory cortex, we identified structural brain similarities for parent-child pairs of the 1st TTG (Heschl's gyrus, HG) (in terms of volume, area and thickness for the right HG, and shape for the left HG) and of the lateralization of all TTG(s) surface area for father-child pairs. Both the HG and TTG-lateralization findings were significantly more likely for parent-child dyads than for unrelated adult-child pairs. Furthermore, we established characteristics of parents' TTG that are related to better reading abilities in children: fathers' small left HG, and a small ratio of HG to planum temporale. Our results suggest intergenerational transmission of specific structural features of the auditory cortex; these may arise from genetics and/or from shared environment.

Detecting genetic effects on phenotype variability to capture gene-by-environment interactions: a systematic method comparison

Genetically associated phenotypic variability has been widely observed across organisms and traits, including in humans. Both gene-gene and gene-environment interactions can lead to an increase in genetically associated phenotypic variability. Therefore, detecting the underlying genetic variants, or variance Quantitative Trait Loci (vQTLs), can provide novel insights into complex traits. Established approaches to detect vQTLs apply different methodologies from variance-only approaches to mean-variance joint tests, but a comprehensive comparison of these methods is lacking. Here, we review available methods to detect vQTLs in humans, carry out a simulation study to assess their performance under different biological scenarios of gene-environment interactions, and apply the optimal approaches for vQTL identification to gene expression data. Overall, with a minor allele frequency (MAF) of less than 0.2, the squared residual value linear model (SVLM) and the deviation regression model (DRM) are optimal when the data follow normal and non-normal distributions, respectively. In addition, the Brown-Forsythe (BF) test is one of the optimal methods when the MAF is 0.2 or larger, irrespective of phenotype distribution. Additionally, a larger sample size and more balanced sample distribution in different exposure categories increase the power of BF, SVLM, and DRM. Our results highlight vQTL detection methods that perform optimally under realistic simulation settings and show that their relative performance depends on the phenotype distribution, allele frequency, sample size, and the type of exposure in the interaction model underlying the vQTL.

Phenome-wide analyses identify an association between the parent-of-origin effects dependent methylome and the rate of aging in humans

BACKGROUND

The variation in the rate at which humans age may be rooted in early events acting through the genomic regions that are influenced by such events and subsequently are related to health phenotypes in later life. The parent-of-origin-effect (POE)-regulated methylome includes regions enriched for genetically controlled imprinting effects (the typical type of POE) and regions influenced by environmental effects associated with parents (the atypical POE). This part of the methylome is heavily influenced by early events, making it a potential route connecting early exposures, the epigenome, and aging. We aim to test the association of POE-CpGs with early and later exposures and subsequently with health-related phenotypes and adult aging.

RESULTS

We perform a phenome-wide association analysis for the POE-influenced methylome using GS:SFHS (Ndiscovery = 5087, Nreplication = 4450). We identify and replicate 92 POE-CpG-phenotype associations. Most of the associations are contributed by the POE-CpGs belonging to the atypical class where the most strongly enriched associations are with aging (DNAmTL acceleration), intelligence, and parental (maternal) smoking exposure phenotypes. A proportion of the atypical POE-CpGs form co-methylation networks (modules) which are associated with these phenotypes, with one of the aging-associated modules displaying increased within-module methylation connectivity with age. The atypical POE-CpGs also display high levels of methylation heterogeneity, fast information loss with age, and a strong correlation with CpGs contained within epigenetic clocks.

CONCLUSIONS

These results identify the association between the atypical POE-influenced methylome and aging and provide new evidence for the "early development of origin" hypothesis for aging in humans.

Phenome-wide analysis identifies parent-of-origin effects on the human methylome associated with changes in the rate of aging

Variation in the rate at which humans age may be rooted in early life events acting through genomic regions that are influenced by such events and subsequently are related to health phenotypes in later life. The parent-of-origin-effect (POE)-regulated methylome includes regions either enriched for genetically controlled imprinting effects (the typical type of POE) or atypical POE introduced by environmental effects associated with parents. This part of the methylome is heavily influenced by early life events, making it a potential route connecting early environmental exposures, the epigenome and the rate of aging. Here, we aim to test the association of POE-influenced methylation of CpG dinucleotides (POE-CpG sites) with early and later environmental exposures and subsequently with health-related phenotypes and adult aging phenotypes. We do this by performing phenome-wide association analyses of the POE-influenced methylome using a large family-based population cohort (GS:SFHS, Ndiscovery=5,087, Nreplication=4,450). At the single CpG level, 92 associations of POE-CpGs with phenotypic variation were identified and replicated. Most of the associations were contributed by POE-CpGs belonging to the atypical class and the most strongly enriched associations were with aging (DNAmTL acceleration), intelligence and parental (maternal) smoking exposure phenotypes. We further found that a proportion of the atypical-POE-CpGs formed co-methylation networks (modules) which are associated with these phenotypes, with one of the aging-associated modules displaying increased internal module connectivity (strength of methylation correlation across constituent CpGs) with age. Atypical POE-CpGs also displayed high levels of methylation heterogeneity and epigenetic drift (i.e. information loss with age) and a strong correlation with CpGs contained within epigenetic clocks. These results identified associations between the atypical-POE-influenced methylome and aging and provided new evidence for the "early development of origin" hypothesis for aging in humans.

Inference of putative cell-type-specific imprinted regulatory elements and genes during human neuronal differentiation

Genomic imprinting results in gene expression bias caused by parental chromosome of origin and occurs in genes with important roles during human brain development. However, the cell-type and temporal specificity of imprinting during human neurogenesis is generally unknown. By detecting within-donor allelic biases in chromatin accessibility and gene expression that are unrelated to cross-donor genotype, we inferred imprinting in both primary human neural progenitor cells and their differentiated neuronal progeny from up to 85 donors. We identified 43/20 putatively imprinted regulatory elements (IREs) in neurons/progenitors, and 133/79 putatively imprinted genes in neurons/progenitors. Although 10 IREs and 42 genes were shared between neurons and progenitors, most putative imprinting was only detected within specific cell types. In addition to well-known imprinted genes and their promoters, we inferred novel putative IREs and imprinted genes. Consistent with both DNA methylation-based and H3K27me3-based regulation of imprinted expression, some putative IREs also overlapped with differentially methylated or histone-marked regions. Finally, we identified a progenitor-specific putatively imprinted gene overlapping with copy number variation that is associated with uniparental disomy-like phenotypes. Our results can therefore be useful in interpreting the function of variants identified in future parent-of-origin association studies.

A family-based study of genetic and epigenetic effects across multiple neurocognitive, motor, social-cognitive and social-behavioral functions

Many psychiatric and neurodevelopmental disorders are known to be heritable, but studies trying to elucidate the genetic architecture of such traits often lag behind studies of somatic traits and diseases. The reasons as to why relatively few genome-wide significant associations have been reported for such traits have to do with the sample sizes needed for the detection of small effects, the difficulty in defining and characterizing the phenotypes, partially due to overlaps in affected underlying domains (which is especially true for cognitive phenotypes), and the complex genetic architectures of the phenotypes, which are not wholly captured in traditional case-control GWAS designs. We aimed to tackle the last two issues by performing GWASs of eight quantitative neurocognitive, motor, social-cognitive and social-behavioral traits, which may be considered endophenotypes for a variety of psychiatric and neurodevelopmental conditions, and for which we employed models capturing both general genetic association and parent-of-origin effects, in a family-based sample comprising 402 children and their parents (mostly family trios). We identified 48 genome-wide significant associations across several traits, of which 3 also survived our strict study-wide quality criteria. We additionally performed a functional annotation of implicated genes, as most of the 48 associations were with variants within protein-coding genes. In total, our study highlighted associations with five genes (TGM3, CACNB4, ANKS1B, CSMD1 and SYNE1) associated with measures of working memory, processing speed and social behavior. Our results thus identify novel associations, including previously unreported parent-of-origin associations with relevant genes, and our top results illustrate new potential gene → endophenotype → disorder pathways.

Detection of Parent-of-Origin Effects for the Variants Associated With Behavioral Disinhibition in the MCTFR Data

Behavioral disinhibition is one of the important characteristics of many mental diseases. It has been reported in literature that serious behavioral disinhibition will affect people's health and greatly reduce people's quality of life. Meanwhile, behavioral disinhibition can easily lead to illegal drug abuse and violent crimes, etc., which will bring great harm to the society. At present, large-scale genome-wide association analysis has identified many loci associated with behavioral disinhibition. However, these studies have not incorporated the parent-of-origin effects (POE) into analysis, which may ignore or underestimate the genetic effects of loci on behavioral disinhibition. Therefore, in this article, we analyzed the five phenotypes related to behavioral disinhibition in the Minnesota Center for Twin and Family Research data (nicotine, alcohol consumption, alcohol dependence, illicit drugs, and non-substance use related behavioral disinhibition), to further explore the POE of variants on behavioral disinhibition. We applied a linear mixed model to test for the POE at a genome-wide scale on five transformed phenotypes, and found nine SNPs with statistically significant POE at the significance level of 5 × 10-8. Among them, SNPs rs4141854, rs9394515, and rs4711553 have been reported to be associated with two neurological disorders (restless legs syndrome and Tourette's syndrome) which are related to behavioral disinhibition; SNPs rs12960235 and rs715351 have been found to be associated with head and neck squamous cell carcinoma, skin cancer and type I diabetes, while both SNPs have not been identified to be related to behavioral disinhibition in literature; SNPs rs704833, rs6837925, rs1863548, and rs11067062 are novel loci identified in this article, and their function annotations have not been reported in literature. Follow-up study in molecular genetics is needed to verify whether they are surely related to behavioral disinhibition.

Parent-of-origin effects propagate through networks to shape metabolic traits

Parent-of-origin effects are unexpectedly common in complex traits, including metabolic and neurological traits. Parent-of-origin effects can be modified by the environment, but the architecture of these gene-by-environmental effects on phenotypes remains to be unraveled. Previously, quantitative trait loci (QTL) showing context-specific parent-of-origin effects on metabolic traits were mapped in the F16 generation of an advanced intercross between LG/J and SM/J inbred mice. However, these QTL were not enriched for known imprinted genes, suggesting another mechanism is needed to explain these parent-of-origin effects phenomena. We propose that non-imprinted genes can generate complex parent-of-origin effects on metabolic traits through interactions with imprinted genes. Here, we employ data from mouse populations at different levels of intercrossing (F0, F1, F2, F16) of the LG/J and SM/J inbred mouse lines to test this hypothesis. Using multiple populations and incorporating genetic, genomic, and physiological data, we leverage orthogonal evidence to identify networks of genes through which parent-of-origin effects propagate. We identify a network comprised of three imprinted and six non-imprinted genes that show parent-of-origin effects. This epistatic network forms a nutritional responsive pathway and the genes comprising it jointly serve cellular functions associated with growth. We focus on two genes, Nnat and F2r, whose interaction associates with serum glucose levels across generations in high-fat-fed females. Single-cell RNAseq reveals that Nnat expression increases and F2r expression decreases in pre-adipocytes along an adipogenic trajectory, a result that is consistent with our observations in bulk white adipose tissue.

Lipid-Associated Variants near ANGPTL3 and LPL Show Parent-of-Origin Specific Effects on Blood Lipid Levels and Obesity

Parent-of-origin effects (POE) and sex-specific parental effects have been reported for plasma lipid levels, and a strong relationship exists between dyslipidemia and obesity. We aim to explore whether genetic variants previously reported to have an association to lipid traits also show POE on blood lipid levels and obesity. Families from the Botnia cohort and the Hungarian Transdanubian Biobank (HTB) were genotyped for 12 SNPs, parental origin of alleles were inferred, and generalized estimating equations were modeled to assess parental-specific associations with lipid traits and obesity. POE were observed for the variants at the TMEM57, DOCK7/ANGPTL3, LPL, and APOA on lipid traits, the latter replicated in HTB. Sex-specific parental effects were also observed; variants at ANGPTL3/DOCK7 showed POE on lipid traits and obesity in daughters only, while those at LPL and TMEM57 showed POE on lipid traits in sons. Variants at LPL and DOCK7/ANGPTL3 showed POE on obesity-related traits in Botnia and HTB, and POE effects on obesity were seen to a higher degree in daughters. This highlights the need to include analysis of POEs in genetic studies of complex traits.

Differential Allelic Expression among Long Non-Coding RNAs

Long non-coding RNAs (lncRNA) comprise a diverse group of non-protein-coding RNAs >200 bp in length that are involved in various normal cellular processes and disease states, and can affect coding gene expression through mechanisms in cis or in trans. Since the discovery of the first functional lncRNAs transcribed by RNA Polymerase II, H19 and Xist, many others have been identified and noted for their unusual transcriptional pattern, whereby expression from one chromosome homolog is strongly favored over the other, also known as mono-allelic or differential allelic expression. lncRNAs with differential allelic expression have been observed to play critical roles in developmental gene regulation, chromosome structure, and disease. Here, we will focus on known examples of differential allelic expression of lncRNAs and highlight recent research describing functional lncRNAs expressed from both imprinted and random mono-allelic expression domains.

Genome-wide survey of parent-of-origin-specific associations across clinical traits derived from electronic health records

Parent-of-origin (PoO) effects refer to the differential phenotypic impacts of genetic variants dependent on their parental inheritance due to imprinting. While PoO effects can influence complex traits, they may be poorly captured by models that do not differentiate the parental origin of the variant. The aim of this study was to conduct a genome-wide screen for PoO effects on a broad range of clinical traits derived from electronic health records (EHR) in the DiscovEHR study enriched with familial relationships. Using pairwise kinship estimates from genetic data and demographic data, we identified 22,051 offspring among 134,049 individuals in the DiscovEHR study. PoO of ~9 million variants was assigned in the offspring by comparing offspring and parental genotypes and haplotypes. We then performed genome-wide PoO association analyses across 154 quantitative and 611 binary traits extracted from EHR. Of the 732 significant PoO associations identified (p < 5 × 10^-8), we attempted to replicate 274 PoO associations in the UK Biobank study with 5,015 offspring and replicated 9 PoO associations (p < 0.05). In summary, our study implements a bioinformatic and statistical approach to examine PoO effects genome-wide in a large population study enriched with familial relationships and systematically characterizes PoO effects on hundreds of clinical traits derived from EHR. Our results suggest that, while the statistical power to detect PoO effects remains modest yet, accurately modeling PoO effects has the potential to find new associations that may have been missed by the standard additive model, further enhancing the mechanistic understanding of genetic influence on complex traits.

Sex-specific differences in peripheral blood leukocyte transcriptional response to LPS are enriched for HLA region and X chromosome genes

Sex-specific differences in prevalence are well documented for many common, complex diseases, especially for immune-mediated diseases, yet the precise mechanisms through which factors associated with biological sex exert their effects throughout life are not well understood. We interrogated sex-specific transcriptional responses of peripheral blood leukocytes (PBLs) to innate immune stimulation by lipopolysaccharide (LPS) in 46 male and 66 female members of the Hutterite community, who practice a communal lifestyle. We identified 1217 autosomal and 54 X-linked genes with sex-specific responses to LPS, as well as 71 autosomal and one X-linked sex-specific expression quantitative trait loci (eQTLs). Despite a similar proportion of the 15 HLA genes responding to LPS compared to all expressed autosomal genes, there was a significant over-representation of genes with sex by treatment interactions among HLA genes. We also observed an enrichment of sex-specific differentially expressed genes in response to LPS for X-linked genes compared to autosomal genes, suggesting that HLA and X-linked genes may disproportionately contribute to sex disparities in risk for immune-mediated diseases.

Evidence for germline non-genetic inheritance of human phenotypes and diseases

It is becoming increasingly apparent that certain phenotypes are inherited across generations independent of the information contained in the DNA sequence, by factors in germ cells that remain largely uncharacterized. As evidence for germline non-genetic inheritance of phenotypes and diseases continues to grow in model organisms, there are fewer reports of this phenomenon in humans, due to a variety of complications in evaluating this mechanism of inheritance in humans. This review summarizes the evidence for germline-based non-genetic inheritance in humans, as well as the significant challenges and important caveats that must be considered when evaluating this process in human populations. Most reports of this process evaluate the association of a lifetime exposure in ancestors with changes in DNA methylation or small RNA expression in germ cells, as well as the association between ancestral experiences and the inheritance of a phenotype in descendants, down to great-grandchildren in some cases. Collectively, these studies provide evidence that phenotypes can be inherited in a DNA-independent manner; the extent to which this process contributes to disease development, as well as the cellular and molecular regulation of this process, remain largely undefined.

Quantitative genome-wide association analyses of receptive language in the Danish High Risk and Resilience Study

BACKGROUND

One of the most basic human traits is language. Linguistic ability, and disability, have been shown to have a strong genetic component in family and twin studies, but molecular genetic studies of language phenotypes are scarce, relative to studies of other cognitive traits and neurodevelopmental phenotypes. Moreover, most genetic studies examining such phenotypes do not incorporate parent-of-origin effects, which could account for some of the heritability of the investigated trait. We performed a genome-wide association study of receptive language, examining both child genetic effects and parent-of-origin effects.

RESULTS

Using a family-based cohort with 400 children with receptive language scores, we found a genome-wide significant paternal parent-of-origin effect with a SNP, rs11787922, on chromosome 9q21.31, whereby the T allele reduced the mean receptive language score by ~ 23, constituting a reduction of more than 1.5 times the population SD (P = 1.04 × 10-8). We further confirmed that this association was not driven by broader neurodevelopmental diagnoses in the child or a family history of psychiatric diagnoses by incorporating covariates for the above and repeating the analysis.

CONCLUSIONS

Our study reports a genome-wide significant association for receptive language skills; to our knowledge, this is the first documented genome-wide significant association for this phenotype. Furthermore, our study illustrates the importance of considering parent-of-origin effects in association studies, particularly in the case of cognitive or neurodevelopmental traits, in which parental genetic data are not always incorporated.

Searching for parent-of-origin effects on cardiometabolic traits in imprinted genomic regions

Cardiometabolic traits pose a major global public health burden. Large-scale genome-wide association studies (GWAS) have identified multiple loci accounting for up to 30% of the genetic variance in complex traits such as cardiometabolic traits. However, the contribution of parent-of-origin effects (POEs) to complex traits has been largely ignored in GWAS. Family-based studies enable the assessment of POEs in genetic association analyses. We investigated POEs on a range of complex traits in 3 family-based studies. The discovery phase was carried out in large pedigrees from the Kibbutzim Family Study (n = 901 individuals) and in 872 parent-offspring trios from the Jerusalem Perinatal Study. Focusing on imprinted genomic regions, we examined parent-specific associations with 12 complex traits (i.e., body-size, blood pressure, lipids), mostly cardiometabolic risk traits. Forty five of the 11,967 SNPs initially found to have POE were evaluated for replication (p value < 1 × 10-4) in Framingham Heart Study families (max n = 8000 individuals). Three common variants yielded evidence of POE in the meta-analysis. Two variants, located on chr6 in the HLA region, showed a paternal effect on height (rs1042136: βpaternal = -0.023, p value = 1.5 × 10-8 and rs1431403: βpaternal = -0.011, p value = 5.4 × 10-6). The corresponding maternally-derived effects were statistically nonsignificant. The variant rs9332053, located on chr13 in RCBTB2 gene, demonstrated a maternal effect on hip circumference (βmaternal = -4.24, p value = 9.6 × 10-6; βpaternal = 1.29, p value = 0.23). These findings provide evidence for the utility of incorporating POEs into association studies of cardiometabolic traits, especially anthropometric traits. The study highlights the benefits of using family-based data for deciphering the genetic architecture of complex traits.

Parent of origin genetic effects on methylation in humans are common and influence complex trait variation

Parent-of-origin effects (POE) exist when there is differential expression of alleles inherited from the two parents. A genome-wide scan for POE on DNA methylation at 639,238 CpGs in 5,101 individuals identifies 733 independent methylation CpGs potentially influenced by POE at a false discovery rate ≤ 0.05 of which 331 had not previously been identified. Cis and trans methylation quantitative trait loci (mQTL) regulate methylation variation through POE at 54% (399/733) of the identified POE-influenced CpGs. The combined results provide strong evidence for previously unidentified POE-influenced CpGs at 171 independent loci. Methylation variation at 14 of the POE-influenced CpGs is associated with multiple metabolic traits. A phenome-wide association analysis using the POE mQTL SNPs identifies a previously unidentified imprinted locus associated with waist circumference. These results provide a high resolution population-level map for POE on DNA methylation sites, their local and distant regulators and potential consequences for complex traits.