Investigation of parent-of-origin effects induced by fenofibrate treatment on triglycerides levels

Background

Genome-wide association studies performed on triglycerides (TGs) have not accounted for epigenetic mechanisms that may partially explain trait heritability.

Results

Parent-of-origin (POO) effect association analyses using an agnostic approach or a candidate approach were performed for pretreatment TG levels, posttreatment TG levels, and pre- and posttreatment TG-level differences in the real GAW20 family data set. We detected 22 genetic variants with suggestive POO effects with at least 1 phenotype (P ≤ 10^− 5). We evaluated the association of these 22 significant genetic variants showing POO effects with close DNA methylation probes associated with TGs. A total of 18 DNA methylation probes located in the vicinity of the 22 SNPs were associated with at least 1 phenotype and 6 SNP-probe pairs were associated with DNA methylation probes at the nominal level of P < 0.05, among which 1 pair presented evidence of POO effect. Our analyses identified a paternal effect of SNP rs301621 on the difference between pre- and posttreatment TG levels (P = 1.2 × 10^− 5). This same SNP showed evidence for a maternal effect on methylation levels of a nearby probe (cg10206250; P = 0.01). Using a causal inference test we established that the observed POO effect of rs301621 was not mediated by DNA methylation at cg10206250.

Conclusions

We performed POO effect association analyses of SNPs with TGs, as well as association analyses of SNPs with DNA methylation probes. These analyses, which were followed by a causal inference test, established that the paternal effect at the SNP rs301621 is induced by treatment and is not mediated by methylation level at cg10206250.

     

Asymmetric paternal effect on offspring size linked to parent‐of‐origin expression of an insulin‐like growth factor

Sexual reproduction brings together reproductive partners whose long‐term interests often differ, raising the possibility of conflict over their reproductive investment. Males that enhance maternal investment in their offspring gain fitness benefits, even if this compromises future reproductive investment by iteroparous females. When the conflict occurs at a genomic level, it may be uncovered by crossing divergent populations, as a mismatch in the coevolved patterns of paternal manipulation and maternal resistance may generate asymmetric embryonic growth. We report such an asymmetry in reciprocal crosses between populations of the fish Girardinichthys multiradiatus. We also show that a fragment of a gene which can influence embryonic growth (Insulin‐Like Growth Factor 2; igf2) exhibits a parent‐of‐origin methylation pattern, where the maternally inherited igf2 allele has much more 5′ cytosine methylation than the paternally inherited allele. Our findings suggest that male manipulation of maternal investment may have evolved in fish, while the parent‐of‐origin methylation pattern appears to be a potential candidate mechanism modulating this antagonistic coevolution process. However, disruption of other coadaptive processes cannot be ruled out, as these can lead to similar effects as conflict.     

Genome‐wide screening for DNA variants associated with reading and language traits

Reading and language abilities are heritable traits that are likely to share some genetic influences with each other. To identify pleiotropic genetic variants affecting these traits, we first performed a genome‐wide association scan (GWAS) meta‐analysis using three richly characterized datasets comprising individuals with histories of reading or language problems, and their siblings. GWAS was performed in a total of 1862 participants using the first principal component computed from several quantitative measures of reading‐ and language‐related abilities, both before and after adjustment for performance IQ. We identified novel suggestive associations at the SNPs rs59197085 and rs5995177 (uncorrected P ≈ 10^–7 for each SNP), located respectively at the CCDC136/FLNC and RBFOX2 genes. Each of these SNPs then showed evidence for effects across multiple reading and language traits in univariate association testing against the individual traits. FLNC encodes a structural protein involved in cytoskeleton remodelling, while RBFOX2 is an important regulator of alternative splicing in neurons. The CCDC136/FLNC locus showed association with a comparable reading/language measure in an independent sample of 6434 participants from the general population, although involving distinct alleles of the associated SNP. Our datasets will form an important part of on‐going international efforts to identify genes contributing to reading and language skills.     

Genome‐wide association study of word reading: Overlap with risk genes for neurodevelopmental disorders

Reading disabilities (RD) are the most common neurocognitive disorder, affecting 5% to 17% of children in North America. These children often have comorbid neurodevelopmental/psychiatric disorders, such as attention deficit/hyperactivity disorder (ADHD). The genetics of RD and their overlap with other disorders is incompletely understood. To contribute to this, we performed a genome‐wide association study (GWAS) for word reading. Then, using summary statistics from neurodevelopmental/psychiatric disorders, we computed polygenic risk scores (PRS) and used them to predict reading ability in our samples. This enabled us to test the shared aetiology between RD and other disorders. The GWAS consisted of 5.3 million single nucleotide polymorphisms (SNPs) and two samples; a family‐based sample recruited for reading difficulties in Toronto (n = 624) and a population‐based sample recruited in Philadelphia [Philadelphia Neurodevelopmental Cohort (PNC)] (n = 4430). The Toronto sample SNP‐based analysis identified suggestive SNPs (P ~ 5 × 10^?7) in the ARHGAP23 gene, which is implicated in neuronal migration/axon pathfinding. The PNC gene‐based analysis identified significant associations (P < 2.72 × 10^?6) for LINC00935 and CCNT1, located in the region of the KANSL2/CCNT1/LINC00935/SNORA2B/SNORA34/MIR4701/ADCY6 genes on chromosome 12q, with near significant SNP‐based analysis. PRS identified significant overlap between word reading and intelligence (R² = 0.18, P = 7.25 × 10^?181), word reading and educational attainment (R² = 0.07, P = 4.91 × 10^?48) and word reading and ADHD (R² = 0.02, P = 8.70 × 10^?6; threshold for significance = 7.14 × 10^?3). Overlap was also found between RD and autism spectrum disorder (ASD) as top‐ranked genes were previously implicated in autism by rare and copy number variant analyses. These findings support shared risk between word reading, cognitive measures, educational outcomes and neurodevelopmental disorders, including ASD.     

Genome‐wide association analysis identifies quantitative trait loci for growth in a Landrace purebred population

Growth‐related traits are complex and economically important in the livestock industry. The aim of this study was to identify quantitative trait loci (QTL) and the associated positional candidate genes affecting growth in pigs. A genome‐wide association study (GWAS) was performed using the porcine single‐nucleotide polymorphism (SNP) 60K bead chip. A mixed‐effects model and linear regression approach were used for the GWAS. The data used in the study included 490 purebred Landrace pigs. All experimental animals were genotyped with 39 438 SNPs located throughout the pig autosomes. We identified a strong association between a SNP marker on chromosome 16 and body weight at 71 days of age (ALGA0092396, P = 5.35 × 10^?9, Bonferroni adjusted P < 0.05). The SNP marker was located near the genomic region containing IRX4, which encodes iroquois homeobox 4. This SNP marker could be useful in the selective breeding program after validating its effect on other populations.     

Maternal and paternal effects on the germination time of non‐dormant seeds of a monocarpic perennial species,Aster kantoensis (Compositae)

In plant species producing non‐dormant seeds, the germination time (from the start of imbibition to radicle emergence) is the main factor determining the timing of seedling emergence. We investigated maternal and paternal genetic effects on the germination time of non‐dormant seeds of a monocarpic perennial, Aster kantoensis Kitamura (Compositae). Three sets of reciprocal diallel crosses among five plants were conducted to produce genetic variation in seeds, and the germination time of the progeny of each parent was determined. The effects of the maternal parent and the interaction of maternal and paternal parents on the germination time of progeny were significant in all sets, and the effect of the paternal parent was significant in two of the three sets. This result means that the germination time of the progeny of a maternal or paternal parent can vary with the genotype of its mating partners. Because variation in the emergence time of seedlings contributes to avoiding seedling loss owing to unpredictable environmental changes, genetic variation in the germination time among the progeny of each parent mating with multiple partners could contribute to the establishment of the parent's seedlings in species producing non‐dormant seeds in the field.     

Paternal epigenetic effects of population density on locust phase‐related characteristics associated with heat‐shock protein expression

Many species exhibit transgenerational plasticity by which environmental cues experienced by either parent can be transmitted to their offspring, resulting in phenotypic variants in offspring to match ancestral environments. However, the manner by which paternal experiences affect offspring plasticity through epigenetic inheritance in animals generally remains unclear. In this study, we examined the transgenerational effects of population density on phase‐related traits in the migratory locust Locusta migratoria. Using an experimental design that explicitly controls genetic background, we found that the effects of crowd or isolation rearing on phase plasticity could be inherited to the offspring. The isolation of gregarious locusts resulted in reduced weight in offspring eggs and altered morphometric traits in hatchlings, whereas crowding of solitarious locusts exhibited opposite effects. The consequences of density changes were transmitted by both maternal and paternal inheritance, although the expression of paternal effects was not as pronounced as that of maternal effects. Prominent expression of heat‐shock proteins (Hsps), such as Hsp90, Hsp70 and Hsp20.6, could be triggered by density changes. Hsps were significantly upregulated upon crowding but downregulated upon isolation. The variation in parental Hsp expression was also transmitted to the offspring, in which the pattern of inheritance was consistent with that of phase characteristics. These results revealed a paternal effect on phase polyphenism and Hsp expression induced by population density, and defined a model system that could be used to study the paternal epigenetic inheritance of environmental changes.     

Genome‐wide Linkage and Association Analyses to Identify Genes Influencing Adiponectin Levels: The GEMS Stud

Adiponectin has a variety of metabolic effects on obesity, insulin sensitivity, and atherosclerosis. To identify genes influencing variation in plasma adiponectin levels, we performed genome‐wide linkage and association scans of adiponectin in two cohorts of subjects recruited in the Genetic Epidemiology of Metabolic Syndrome Study. The genome‐wide linkage scan was conducted in families of Turkish and southern European (TSE, n = 789) and Northern and Western European (NWE, N = 2,280) origin. A whole genome association (WGA) analysis (500K Affymetrix platform) was carried out in a set of unrelated NWE subjects consisting of approximately 1,000 subjects with dyslipidemia and 1,000 overweight subjects with normal lipids. Peak evidence for linkage occurred at chromosome 8p23 in NWE subjects (lod = 3.10) and at chromosome 3q28 near ADIPOQ, the adiponectin structural gene, in TSE subjects (lod = 1.70). In the WGA analysis, the single‐nucleotide polymorphisms (SNPs) most strongly associated with adiponectin were rs3774261 and rs6773957 (P < 10^?7). These two SNPs were in high linkage disequilibrium (r² = 0.98) and located within ADIPOQ. Interestingly, our fourth strongest region of association (P < 2 × 10^?5) was to an SNP within CDH13, whose protein product is a newly identified receptor for high‐molecular‐weight species of adiponectin. Through WGA analysis, we confirmed previous studies showing SNPs within ADIPOQ to be strongly associated with variation in adiponectin levels and further observed these to have the strongest effects on adiponectin levels throughout the genome. We additionally identified a second gene (CDH13) possibly influencing variation in adiponectin levels. The impact of these SNPs on health and disease has yet to be determined.     

A genome‐wide association study for reading and language abilities in two population cohorts

Candidate genes have been identified for both reading and language, but most of the heritable variance in these traits remains unexplained. Here, we report a genome‐wide association meta‐analysis of two large cohorts: population samples of Australian twins and siblings aged 12–25 years (n = 1177 from 538 families), and a younger cohort of children of the UK Avon Longitudinal Study of Parents and their Children (aged 8 and 9 years; maximum n = 5472). Suggestive association was indicated for reading measures and non‐word repetition (NWR), with the greatest support found for single nucleotide polymorphisms (SNPs) in the pseudogene, ABCC13 (P = 7.34 × 10^?8), and the gene, DAZAP1 (P = 1.32 × 10^?6). Gene‐based analyses showed significant association (P < 2.8 × 10^?6) for reading and spelling with genes CD2L1, CDC2L2 and RCAN3 in two loci on chromosome 1. Some support was found for the same SNPs having effects on both reading skill and NWR, which is compatible with behavior genetic evidence for influences of reading acquisition on phonological‐task performance. The results implicate novel candidates for study in additional cohorts for reading and language abilities.     

Genome‐wide association of drought‐related and biomass traits with HapMap SNPs in Medicago truncatula

Improving drought tolerance of crop plants is a major goal of plant breeders. In this study, we characterized biomass and drought‐related traits of 220 Medicago truncatula HapMap accessions. Characterized traits included shoot biomass, maximum leaf size, specific leaf weight, stomatal density, trichome density and shoot carbon‐13 isotope discrimination (δ¹³C) of well‐watered M. truncatula plants, and leaf performance in vitro under dehydration stress. Genome‐wide association analyses were carried out using the general linear model (GLM), the standard mixed linear model (MLM) and compressed MLM (CMLM) in TASSEL, which revealed significant overestimation of P‐values by CMLM. For each trait, candidate genes and chromosome regions containing SNP markers were found that are in significant association with the trait. For plant biomass, a 0.5 Mbp region on chromosome 2 harbouring a plasma membrane intrinsic protein, PIP2, was discovered that could potentially be targeted to increase dry matter yield. A protein disulfide isomerase‐like protein was found to be tightly associated with both shoot biomass and leaf size. A glutamate‐cysteine ligase and an aldehyde dehydrogenase family protein with Arabidopsis homologs strongly expressed in the guard cells were two of the top genes identified by stomata density genome‐wide association studies analysis.     

The miR‐379/miR‐410 cluster at the imprinted Dlk1‐Dio3 domain controls neonatal metabolic adaptation

In mammals, birth entails complex metabolic adjustments essential for neonatal survival. Using a mouse knockout model, we identify crucial biological roles for the miR‐379/miR‐410 cluster within the imprinted Dlk1‐Dio3 region during this metabolic transition. The miR‐379/miR‐410 locus, also named C14MC in humans, is the largest known placental mammal‐specific miRNA cluster, whose 39 miRNA genes are expressed only from the maternal allele. We found that heterozygote pups with a maternal—but not paternal—deletion of the miRNA cluster display partially penetrant neonatal lethality with defects in the maintenance of energy homeostasis. This maladaptive metabolic response is caused, at least in part, by profound changes in the activation of the neonatal hepatic gene expression program, pointing to as yet unidentified regulatory pathways that govern this crucial metabolic transition in the newborn's liver. Not only does our study highlight the physiological importance of miRNA genes that recently evolved in placental mammal lineages but it also unveils additional layers of RNA‐mediated gene regulation at the Dlk1‐Dio3 domain that impose parent‐of‐origin effects on metabolic control at birth and have likely contributed to mammal evolution.     

Genome‐wide association study of genetic factors related to confectionery intake: Potential roles of the ADIPOQ gene

Objective: The excessive consumption of confectionery might have adverse effects on human health. To screen genetic factors associated with confectionery‐intake frequency, a genome‐wide association study (GWAS) in Japan was conducted. Design and Methods: For the discovery phase (stage 1), we conducted a GWAS of 939 noncancer patients in a cancer hospital. Additive models were used to test associations between genotypes of approximately 500,000 single‐nucleotide polymorphisms (SNPs) and the confectionery‐intake score (based on intake frequency). We followed‐up association signals with P < 1 × 10^?5 and minor allele frequency >0.01 in stage 1 by genotyping the SNPs of 4,491 participants in a cross‐sectional study within a cohort (replication phase [stage 2]). Results: We identified 12 SNPs in stage 1 that were potentially related to confectionery intake. In stage 2, this association was replicated for one SNP (rs822396; P = 0.049 for stage 2 and 4.2 × 10^?5 for stage 1+2) in intron 1 of the ADIPOQ gene, which encodes the adipokine adiponectin. Conclusions: Given the biological plausibility and previous relevant findings, the association of an SNP in the ADIPOQ gene with a preference for confectionery is worthy of follow‐up and provides a good working hypothesis for experimental testing.     

A genome‐wide search for quantitative trait loci affecting the cortical surface area and thickness of Heschl's gyrus

Heschl's gyrus (HG) is a core region of the auditory cortex whose morphology is highly variable across individuals. This variability has been linked to sound perception ability in both speech and music domains. Previous studies show that variations in morphological features of HG, such as cortical surface area and thickness, are heritable. To identify genetic variants that affect HG morphology, we conducted a genome‐wide association scan (GWAS) meta‐analysis in 3054 healthy individuals using HG surface area and thickness as quantitative traits. None of the single nucleotide polymorphisms (SNPs) showed association P values that would survive correction for multiple testing over the genome. The most significant association was found between right HG area and SNP rs72932726 close to gene DCBLD2 (3q12.1; P = 2.77 × 10^?7). This SNP was also associated with other regions involved in speech processing. The SNP rs333332 within gene KALRN (3q21.2; P = 2.27 × 10^?6) and rs143000161 near gene COBLL1 (2q24.3; P = 2.40 × 10^?6) were associated with the area and thickness of left HG, respectively. Both genes are involved in the development of the nervous system. The SNP rs7062395 close to the X‐linked deafness gene POU3F4 was associated with right HG thickness (Xq21.1; P = 2.38 × 10^?6). This is the first molecular genetic analysis of variability in HG morphology.     

GENETIC COMPATIBILITY AFFECTS DIVISION OF LABOR IN THE ARGENTINE ANT LINEPITHEMA HUMILE

Division of labor is central to the organization of insect societies. Within‐colony comparisons between subfamilies of workers (patrilines or matrilines) revealed genetic effects on division of labor in many social insect species. Although this has been taken as evidence for additive genetic effects on division of labor, it has never been experimentally tested. To determine the relative roles of additive and nonadditive genetic effects (e.g., genetic compatibility, epistasis, and parent‐of‐origin imprinting effects) on worker behavior, we performed controlled crosses using the Argentine ant Linepithema humile. Three of the measured behaviors (the efficiency to collect pupae, the foraging propensity, and the distance between non‐brood‐tenders and brood) were affected by the maternal genetic background and the two others (the efficiency to feed larvae and the distance between brood‐tenders and brood) by the paternal genetic background. Moreover, there were significant interactions between the maternal and paternal genetic backgrounds for three of the five behaviors. These results are most consistent with parent‐of‐origin and genetic compatibility effects on division of labor. The finding of nonadditive genetic effects is in strong contrast with the current view and has important consequences for our understanding of division of labor in insect societies.     

Linkage and Genome‐wide Association Analysis of Obesity‐related Phenotypes: Association of Weight With the MGAT1 Gene

As major risk‐factors for diabetes and cardiovascular diseases, the genetic contribution to obesity‐related traits has been of interest for decades. Recently, a limited number of common genetic variants, which have replicated in different populations, have been identified. One approach to increase the statistical power in genetic mapping studies is to focus on populations with increased levels of linkage disequilibrium (LD) and reduced genetic diversity. We have performed joint linkage and genome‐wide association analyses for weight and BMI in 3,448 (linkage) and 3,925 (association) partly overlapping healthy individuals from five European populations. A total of four chromosomal regions (two for weight and two for BMI) showed suggestive linkage (lod >2.69) either in one of the populations or in the joint data. At the genome‐wide level (nominal P < 1.6 × 10^?7, Bonferroni‐adjusted P < 0.05) one single‐nucleotide polymorphism (SNP) (rs12517906) (nominal P = 7.3 × 10^?8) was associated with weight, whereas none with BMI. The SNP associated with weight is located close to MGAT1. The monoacylglycerol acyltransferase (MGAT) enzyme family is known to be involved in dietary fat absorption. There was no overlap between the linkage regions and the associated SNPs. Our results show that genetic effects influencing weight and BMI are shared across diverse European populations, even though some of these populations have experienced recent population bottlenecks and/or been affected by genetic drift. The analysis enabled us to identify a new candidate gene, MGAT1, associated with weight in women.     

LPHN3 and attention‐deficit/hyperactivity disorder: a susceptibility and pharmacogenetic study

Latrophilin 3 (LPHN3) is a brain‐specific member of the G‐protein coupled receptor family associated to both attention‐deficit/hyperactivity disorder (ADHD) genetic susceptibility and methylphenidate (MPH) pharmacogenetics. Interactions of LPHN3 variants with variants harbored in the 11q chromosome improve the prediction of ADHD development and medication response. The aim of this study was to evaluate the role of LPHN3 variants in childhood ADHD susceptibility and treatment response in a naturalistic clinical cohort. The association between LPHN3 and ADHD was evaluated in 523 children and adolescents with ADHD and 132 controls. In the pharmacogenetic study, 172 children with ADHD were investigated. The primary outcome measure was the parent‐rated Swanson, Nolan and Pelham Scale – version IV applied at baseline, first and third months of treatment with MPH. The results reported herein suggest the CGC haplotype derived from single nucleotide polymorphisms (SNPs) rs6813183, rs1355368 and rs734644 as an ADHD risk haplotype (P = 0.02, OR = 1.46). Although non‐significant after multiple testing correction, its interaction with the 11q chromosome SNP rs965560 slightly increases risk (P = 0.03, OR = 1.55). Homozygous individuals for the CGC haplotype showed faster response to MPH treatment as a significant interaction effect between CGC haplotype and treatment over time was observed (P < 0.001). Homozygous individuals for the GT haplotype derived from SNPs rs6551665 and rs1947275 showed a nominally significant interaction with treatment over time (P = 0.04). Our findings replicate previous findings reporting that LPHN3 confers ADHD susceptibility, and moderates MPH treatment response in children and adolescents with ADHD.     

Effects of advanced paternal age on trajectories of social behavior in offspring

Our study is the first investigation of the effects of advanced paternal age (APA) on the developmental trajectory of social behavior in rodent offspring. Given the strong epidemiological association between APA and sexually dimorphic neurodevelopmental disorders that are characterized by abnormalities in social behavior (autism, schizophrenia), we assessed sociability in male and female inbred mice (C57BL/6J) across postnatal development (N = 104) in relation to paternal age. We found differences in early social behavior in both male and female offspring of older breeders, with differences in this social domain persisting into adulthood in males only. We showed that these social deficits were not present in the fathers of these offspring, confirming a de novo origin of an altered social trajectory in the offspring generation. Our results, highly novel in rodent research, support the epidemiological observations in humans and provide evidence for a causal link between APA, age‐related changes in the paternal sperm DNA and neurodevelopmental disorders in their offspring.