Clinical utility of pharmacogenetics in a psychiatric and primary care population

This study evaluated the timing, use, and clinical outcomes of the GeneFolio® Pharmacogenomic Panel in a healthcare setting with patients managed by primary care providers or by psychiatrists. Participants were randomized to receive a pharmacogenetics report at four weeks or 12 weeks. After DNA collection and genetic analysis, pharmacists produced a recommendation report which was given to providers at the randomization week. The four-week group decreased depression severity (PHQ-9 and BDI) faster than the 12-week group (p = 0.0196), and psychiatrists' patients decreased their depression severity faster than primary care patients (PHQ-9 p = 0.0005, BDI p = 0.0218). Mean mental quality of life increased over time (p < 0.0001), but it increased slower for patients taking drugs in the Significant drug-drug-gene interaction category (p = 0.0012). Mental quality of life, depression severity, and clinical outcomes were improved by GeneFolio® pharmacogenomic testing regardless of provider type, with earlier testing improving outcomes sooner.

A genome-wide association study of total child psychiatric problems scores

Substantial genetic correlations have been reported across psychiatric disorders and numerous cross-disorder genetic variants have been detected. To identify the genetic variants underlying general psychopathology in childhood, we performed a genome-wide association study using a total psychiatric problem score. We analyzed 6,844,199 common SNPs in 38,418 school-aged children from 20 population-based cohorts participating in the EAGLE consortium. The SNP heritability of total psychiatric problems was 5.4% (SE = 0.01) and two loci reached genome-wide significance: rs10767094 and rs202005905. We also observed an association of SBF2, a gene associated with neuroticism in previous GWAS, with total psychiatric problems. The genetic effects underlying the total score were shared with common psychiatric disorders only (attention-deficit/hyperactivity disorder, anxiety, depression, insomnia) (rG > 0.49), but not with autism or the less common adult disorders (schizophrenia, bipolar disorder, or eating disorders) (rG < 0.01). Importantly, the total psychiatric problem score also showed at least a moderate genetic correlation with intelligence, educational attainment, wellbeing, smoking, and body fat (rG > 0.29). The results suggest that many common genetic variants are associated with childhood psychiatric symptoms and related phenotypes in general instead of with specific symptoms. Further research is needed to establish causality and pleiotropic mechanisms between related traits.

Genomic and phenotypic insights from an atlas of genetic effects on DNA methylation

Characterizing genetic influences on DNA methylation (DNAm) provides an opportunity to understand mechanisms underpinning gene regulation and disease. In the present study, we describe results of DNAm quantitative trait locus (mQTL) analyses on 32,851 participants, identifying genetic variants associated with DNAm at 420,509 DNAm sites in blood. We present a database of >270,000 independent mQTLs, of which 8.5% comprise long-range (trans) associations. Identified mQTL associations explain 15-17% of the additive genetic variance of DNAm. We show that the genetic architecture of DNAm levels is highly polygenic. Using shared genetic control between distal DNAm sites, we constructed networks, identifying 405 discrete genomic communities enriched for genomic annotations and complex traits. Shared genetic variants are associated with both DNAm levels and complex diseases, but only in a minority of cases do these associations reflect causal relationships from DNAm to trait or vice versa, indicating a more complex genotype-phenotype map than previously anticipated.

Genetic association study of childhood aggression across raters, instruments, and age

Childhood aggressive behavior (AGG) has a substantial heritability of around 50%. Here we present a genome-wide association meta-analysis (GWAMA) of childhood AGG, in which all phenotype measures across childhood ages from multiple assessors were included. We analyzed phenotype assessments for a total of 328 935 observations from 87 485 children aged between 1.5 and 18 years, while accounting for sample overlap. We also meta-analyzed within subsets of the data, i.e., within rater, instrument and age. SNP-heritability for the overall meta-analysis (AGGoverall) was 3.31% (SE = 0.0038). We found no genome-wide significant SNPs for AGGoverall. The gene-based analysis returned three significant genes: ST3GAL3 (P = 1.6E-06), PCDH7 (P = 2.0E-06), and IPO13 (P = 2.5E-06). All three genes have previously been associated with educational traits. Polygenic scores based on our GWAMA significantly predicted aggression in a holdout sample of children (variance explained = 0.44%) and in retrospectively assessed childhood aggression (variance explained = 0.20%). Genetic correlations (rg) among rater-specific assessment of AGG ranged from rg = 0.46 between self- and teacher-assessment to rg = 0.81 between mother- and teacher-assessment. We obtained moderate-to-strong rgs with selected phenotypes from multiple domains, but hardly with any of the classical biomarkers thought to be associated with AGG. Significant genetic correlations were observed with most psychiatric and psychological traits (range [Formula: see text]: 0.19-1.00), except for obsessive-compulsive disorder. Aggression had a negative genetic correlation (rg = ~-0.5) with cognitive traits and age at first birth. Aggression was strongly genetically correlated with smoking phenotypes (range [Formula: see text]: 0.46-0.60). The genetic correlations between aggression and psychiatric disorders were weaker for teacher-reported AGG than for mother- and self-reported AGG. The current GWAMA of childhood aggression provides a powerful tool to interrogate the rater-specific genetic etiology of AGG.

DNA methylation signatures of aggression and closely related constructs: A meta-analysis of epigenome-wide studies across the lifespan

DNA methylation profiles of aggressive behavior may capture lifetime cumulative effects of genetic, stochastic, and environmental influences associated with aggression. Here, we report the first large meta-analysis of epigenome-wide association studies (EWAS) of aggressive behavior (N = 15,324 participants). In peripheral blood samples of 14,434 participants from 18 cohorts with mean ages ranging from 7 to 68 years, 13 methylation sites were significantly associated with aggression (alpha = 1.2 × 10-7; Bonferroni correction). In cord blood samples of 2425 children from five cohorts with aggression assessed at mean ages ranging from 4 to 7 years, 83% of these sites showed the same direction of association with childhood aggression (r = 0.74, p = 0.006) but no epigenome-wide significant sites were found. Top-sites (48 at a false discovery rate of 5% in the peripheral blood meta-analysis or in a combined meta-analysis of peripheral blood and cord blood) have been associated with chemical exposures, smoking, cognition, metabolic traits, and genetic variation (mQTLs). Three genes whose expression levels were associated with top-sites were previously linked to schizophrenia and general risk tolerance. At six CpGs, DNA methylation variation in blood mirrors variation in the brain. On average 44% (range = 3-82%) of the aggression-methylation association was explained by current and former smoking and BMI. These findings point at loci that are sensitive to chemical exposures with potential implications for neuronal functions. We hope these results to be a starting point for studies leading to applications as peripheral biomarkers and to reveal causal relationships with aggression and related traits.

Large-scale association analyses identify host factors influencing human gut microbiome composition

To study the effect of host genetics on gut microbiome composition, the MiBioGen consortium curated and analyzed genome-wide genotypes and 16S fecal microbiome data from 18,340 individuals (24 cohorts). Microbial composition showed high variability across cohorts: only 9 of 410 genera were detected in more than 95% of samples. A genome-wide association study of host genetic variation regarding microbial taxa identified 31 loci affecting the microbiome at a genome-wide significant (P < 5 × 10-8) threshold. One locus, the lactase (LCT) gene locus, reached study-wide significance (genome-wide association study signal: P = 1.28 × 10-20), and it showed an age-dependent association with Bifidobacterium abundance. Other associations were suggestive (1.95 × 10-10 < P < 5 × 10-8) but enriched for taxa showing high heritability and for genes expressed in the intestine and brain. A phenome-wide association study and Mendelian randomization identified enrichment of microbiome trait loci in the metabolic, nutrition and environment domains and suggested the microbiome might have causal effects in ulcerative colitis and rheumatoid arthritis.

Sub30: Protocol for the Sub30 feasibility study of a pre-hospital Extracorporeal membrane oxygenation (ECMO) capable advanced resuscitation team at achieving blood flow within 30 ​min in patients with refractory out-of-hospital cardiac arrest

Background

Out-of-hospital cardiac arrest carries a poor prognosis with survival less than 10% in many patient cohorts. Survival is inversely associated with duration of resuscitation as external chest compressions do not provide sufficient blood flow to prevent irreversible organ damage during a prolonged resuscitation. Extracorporeal membrane oxygenation (ECMO) instituted during cardiac arrest can provide normal physiological blood flows and is termed Extracorporeal Cardio-Pulmonary Resuscitation (ECPR). ECPR may improve survival when used with in-hospital cardiac arrests. This possible survival benefit has not been replicated in trials of out-of-hospital cardiac arrests, possibly because of the additional time it takes to transport the patient to hospital and initiate ECPR. Pre-hospital ECPR may shorten the time between cardiac arrest and physiological blood flows, potentially improving survival. It may also mitigate some of the neurological injury that many survivors suffer.

Methods

Sub30 is a prospective six patient feasibility study. The primary aim is to test whether it is possible to institute ECPR within 30 ​min of collapse in adult patients with refractory out of hospital cardiac arrest (OHCA). The secondary aims are to gather preliminary data on clinical outcomes, resource utilisation, and health economics associated with rapid ECPR delivery in order to plan any subsequent clinical investigation or clinical service. On study days a dedicated fast-response vehicle with ECPR capability will be tasked to out-of-hospital cardiac arrests in an area of London served by Barts Heart Centre. If patients suffer a cardiac arrest refractory to standard advanced resuscitation and meet eligibility criteria, ECPR will be started in the pre-hospital environment.

Discussion

Delivering pre-hospital ECPR within 30 ​min of an out-of-hospital cardiac arrest presents significant ethical, clinical, governance and logistical challenges. Prior to conducting an efficacy study of ECPR the feasibility of timely and safe application must be demonstrated first. Extensive planning, multiple high-fidelity multiagency simulations and a unique collaboration between pre-hospital and in-hospital institutions will allow us to test the feasibility of this intervention in London. The study has been reviewed, refined and endorsed by the International ECMO Network (ECMONet).

Trial registration

Clinicaltrials. gov NCT03700125, prospectively registered October 9, 2018.

Low prevalence of male microchimerism in women with Mayer-Rokitansky-Küster-Hauser syndrome

STUDY QUESTION

Is there an increased prevalence of male microchimerism in women with Mayer–Rokitansky–Küster–Hauser (MRKH) syndrome, as evidence of fetal exposure to blood and anti-Müllerian hormone (AMH) from a (vanished) male co-twin resulting in regression of the Müllerian duct derivatives?

SUMMARY ANSWER

Predominant absence of male microchimerism in adult women with MRKH syndrome does not support our hypothesis that intrauterine blood exchange with a (vanished) male co-twin is the pathophysiological mechanism.

WHAT IS KNOWN ALREADY

The etiology of MRKH is unclear. Research on the phenotype analogous condition in cattle (freemartinism) has yielded the hypothesis that Müllerian duct development is inhibited by exposure to AMH in utero. In cattle, the male co-twin has been identified as the source for AMH, which is transferred via placental blood exchange. In human twins, a similar exchange of cellular material has been documented by detection of chimerism, but it is unknown whether this has clinical consequences.

STUDY DESIGN, SIZE, DURATION

An observational case–control study was performed to compare the presence of male microchimerism in women with MRKH syndrome and control women. Through recruitment via the Dutch patients’ association of women with MRKH (comprising 300 members who were informed by email or regular mail), we enrolled 96 patients between January 2017 and July 2017. The control group consisted of 100 women who reported never having been pregnant.

PARTICIPANTS/MATERIALS, SETTING, METHODS

After written informed consent, peripheral blood samples were obtained by venipuncture, and genomic DNA was extracted. Male microchimerism was detected by Y-chromosome–specific real-time quantitative PCR, with use of DYS14 marker. Possible other sources for microchimerism, for example older brothers, were evaluated using questionnaire data.

MAIN RESULTS AND THE ROLE OF CHANCE

The final analysis included 194 women: 95 women with MRKH syndrome with a mean age of 40.9 years and 99 control women with a mean age of 30.2 years. In total, 54 women (56.8%) were identified as having typical MRKH syndrome, and 41 women (43.2%) were identified as having atypical MRKH syndrome (when extra-genital malformations were present). The prevalence of male microchimerism was significantly higher in the control group than in the MRKH group (17.2% versus 5.3%, P = 0.009). After correcting for age, women in the control group were 5.8 times more likely to have male microchimerism (odds ratio 5.84 (CI 1.59–21.47), P = 0.008). The mean concentration of male microchimerism in the positive samples was 56.0 male genome equivalent per 1 000 000 cells. The prevalence of male microchimerism was similar in women with typical MRKH syndrome and atypical MRKH syndrome (5.6% versus 4.9%, P = 0.884). There were no differences between women with or without microchimerism in occurrence of alternative sources of XY cells, such as older brothers, previous blood transfusion, or history of sexual intercourse.

LIMITATIONS, REASON FOR CAUTION

We are not able to draw definitive conclusions regarding the occurrence of AMH exchange during embryologic development in women with MRKH syndrome. Our subject population includes all adult women and therefore is reliant on long-term prevalence of microchimerism. Moreover, we have only tested blood, and, theoretically, the cells may have grafted anywhere in the body during development. It must also be considered that the exchange of AMH may occur without the transfusion of XY cells and therefore cannot be discovered by chimerism detection.

WIDER IMPLICATIONS OF THE FINDINGS

This is the first study to test the theory that freemartinism causes the MRKH syndrome in humans. The study aimed to test the presence of male microchimerism in women with MRKH syndrome as a reflection of early fetal exposure to blood and AMH from a male (vanished) co-twin. We found that male microchimerism was only present in 5.3% of the women with MRKH syndrome, a significantly lower percentage than in the control group (17.2%). Our results do not provide evidence for an increased male microchimerism in adult women with MRKH as a product of intrauterine blood exchange. However, the significant difference in favor of the control group is of interest to the ongoing discussion on microchimeric cell transfer and the possible sources of XY cells.

STUDY FUNDING/COMPETING INTEREST(S)

None.

TRIAL REGISTRATION NUMBER

Dutch trial register, NTR5961.

Intergenerational Transmission of Education and ADHD: Effects of Parental Genotypes

It remains a challenge to determine whether children resemble their parents due to nature, nurture, or a mixture of both. Here we used a design that exploits the distinction between transmitted and non-transmitted alleles in genetic transmission from parent to offspring. Two separate polygenic scores (PGS) were calculated on the basis of the transmitted and non-transmitted alleles. The effect of the non-transmitted PGS is necessarily mediated by parental phenotypes, insofar as they contribute to the rearing environment of the offspring (genetic nurturing). We calculated transmitted and non-transmitted PGSs associated with adult educational attainment (EA) and PGSs associated with childhood ADHD in a general population sample of trios, i.e. child or adult offspring and their parents (N = 1120-2518). We tested if the EA and ADHD (non-)transmitted PGSs were associated with childhood academic achievement and ADHD in offspring. Based on the earlier findings for shared environment, we hypothesized to find genetic nurturing for academic achievement, but not for ADHD. In adults, both transmitted (R2 = 7.6%) and non-transmitted (R2 = 1.7%) EA PGSs were associated with offspring EA, evidencing genetic nurturing. In children around age 12, academic achievement was associated with the transmitted EA PGSs (R2 = 5.7%), but we found no support for genetic nurturing (R2 ~ 0.1%). The ADHD PGSs were not significantly associated with academic achievement (R2 ~ 0.6%). ADHD symptoms in children were only associated with transmitted EA PGSs and ADHD PGSs (R2 = 1-2%). Based on these results, we conclude that the associations between parent characteristics and offspring outcomes in childhood are mainly to be attributable to the effects of genes that are shared by parents and children.

Urinary Amine and Organic Acid Metabolites Evaluated as Markers for Childhood Aggression: The ACTION Biomarker Study

Biomarkers are of interest as potential diagnostic and predictive instruments in personalized medicine. We present the first urinary metabolomics biomarker study of childhood aggression. We aim to examine the association of urinary metabolites and neurotransmitter ratios involved in key metabolic and neurotransmitter pathways in a large cohort of twins (N = 1,347) and clinic-referred children (N = 183) with an average age of 9.7 years. This study is part of ACTION (Aggression in Children: Unraveling gene-environment interplay to inform Treatment and InterventiON strategies), in which we developed a standardized protocol for large-scale collection of urine samples in children. Our analytical design consisted of three phases: a discovery phase in twins scoring low or high on aggression (N = 783); a replication phase in twin pairs discordant for aggression (N = 378); and a validation phase in clinical cases and matched twin controls (N = 367). In the discovery phase, 6 biomarkers were significantly associated with childhood aggression, of which the association of O-phosphoserine (β = 0.36; SE = 0.09; p = 0.004), and gamma-L-glutamyl-L-alanine (β = 0.32; SE = 0.09; p = 0.01) remained significant after multiple testing. Although non-significant, the directions of effect were congruent between the discovery and replication analyses for six biomarkers and two neurotransmitter ratios and the concentrations of 6 amines differed between low and high aggressive twins. In the validation analyses, the top biomarkers and neurotransmitter ratios, with congruent directions of effect, showed no significant associations with childhood aggression. We find suggestive evidence for associations of childhood aggression with metabolic dysregulation of neurotransmission, oxidative stress, and energy metabolism. Although replication is required, our findings provide starting points to investigate causal and pleiotropic effects of these dysregulations on childhood aggression.

Chimerism in health and potential implications on behavior: A systematic review

In this review, we focus on the phenomenon of chimerism and especially microchimerism as one of the currently underexplored explanations for differences in health and behavior. Chimerism is an amalgamation of cells from two or more unique zygotes within a single organism, with microchimerism defined by a minor cell population of <1%. This article first presents an overview of the primary techniques employed to detect and quantify the presence of microchimerism and then reviews empirical studies of chimerism in mammals including primates and humans. In women, male microchimerism, a condition suggested to be the result of fetomaternal exchange in utero, is relatively easily detected by polymerase chain reaction molecular techniques targeting Y-chromosomal markers. Consequently, studies of chimerism in human diseases have largely focused on diseases with a predilection for females including autoimmune diseases, and female cancers. We detail studies of chimerism in human diseases and also discuss some potential implications in behavior. Understanding the prevalence of chimerism and the associated health outcomes will provide invaluable knowledge of human biology and guide novel approaches for treating diseases.

The genetic architecture of the human cerebral cortex

The cerebral cortex underlies our complex cognitive capabilities, yet little is known about the specific genetic loci that influence human cortical structure. To identify genetic variants that affect cortical structure, we conducted a genome-wide association meta-analysis of brain magnetic resonance imaging data from 51,665 individuals. We analyzed the surface area and average thickness of the whole cortex and 34 regions with known functional specializations. We identified 199 significant loci and found significant enrichment for loci influencing total surface area within regulatory elements that are active during prenatal cortical development, supporting the radial unit hypothesis. Loci that affect regional surface area cluster near genes in Wnt signaling pathways, which influence progenitor expansion and areal identity. Variation in cortical structure is genetically correlated with cognitive function, Parkinson's disease, insomnia, depression, neuroticism, and attention deficit hyperactivity disorder.

Utilization of archived neonatal dried blood spots for genome-wide genotyping

Introduction

Heel pricks are performed on newborns for diagnostic screenings of various pre-symptomatic metabolic and genetic diseases. Excess blood is spotted on Guthrie cards and archived by many states in biobanks for follow-up diagnoses and public health research. However, storage environment may vary across biobanks and across time within biobanks. With increased applications of DNA extracted from spots for genetic studies, identifying factors associated with genotyping success is critical to maximize DNA quality for future studies.

Method

We evaluated 399 blood spots, which were part of a genome-wide association study of childhood leukemia risk in children with Down syndrome, archived at the Michigan Neonatal Biobank between 1992 and 2008. High quality DNA was defined as having post-quality control call rate ≥ 99.0% based on the Illumina GenomeStudio 2.0 GenCall algorithm after processing the samples on the Illumina Infinium Global Screening Array. Bivariate analyses and multivariable logistic regression models were applied to evaluate effects of storage environment and storage duration on DNA genotyping quality.

Results

Both storage environment and duration were associated with sample genotyping call rates (p-values < 0.001). Sample call rates were associated with storage duration independent of storage environment (p-trend = 0.006 for DBS archived in an uncontrolled environment and p-trend = 0.002 in a controlled environment). However, 95% of the total sample had high genotyping quality with a call rate ≥ 95.0%, a standard threshold for acceptable sample quality in many genetic studies.

Conclusion

Blood spot DNA quality was lower in samples archived in uncontrolled storage environments and for samples archived for longer durations. Still, regardless of storage environment or duration, neonatal biobanks including the Michigan Neonatal Biobanks can provide access to large collections of spots with DNA quality acceptable for most genotyping studies.

Abstract P4-05-20: Case report - 16 year life history and genomic evolution of an ER+ HER2- breast cancer

Introduction

Metastasis is the primary cause of cancer-related death, with genomic evolution between primary and metastatic breast cancer only recently being studied in smaller numbers. We reconstructed the tumor evolution of a female patient originally diagnosed with localized ER+ HER2- disease to investigate the molecular changes that triggered relapse after long term aromatase inhibitor (AI) therapy and tracked them during subsequent treatments. The patient was diagnosed in 2002, relapsed in 2012 (regional lymph node), progressed in 2015 (liver metastasis), and passed away in 2018. We performed DNA- and RNA-sequencing on FFPE tumor samples and multiple liquid biopsies. The patient underwent eight lines of therapy including genome guided targeted therapy.

Materials and Methods

Whole exome and transcriptome sequencing was performed on the patient's primary tumor site at diagnosis (T1) and a metastatic lymph node (M1). Targeted DNA sequencing was performed on a metastatic liver sample (M2, 315 genes) and on liquid biopsies (n=11, 70 genes) every 2-3 months from 2015-2018.

Results

The mutational landscape between timepoints T1, M1 and M2 revealed low mutational burden at diagnosis which ultimately doubled at progression with active mutational processes (APOBEC, POLE) showing concordance between timepoints. DNA-based analysis revealed multiple mutations in the PI3K-Akt signaling pathway. Two mutations were shared between the samples: a clonal PIK3CA Q546R and a TP53 E180K mutation detected as clonal in M1/M2 and subclonal in T1. Individual samples carried private driver mutations, ranging from n=4 in T1, n=7 in M1 to n=16 in M2. While T1 harbored a single clonal driver mutation, with additional driver mutations being subclonal, M1 and M2 harbored multiple clonal driver mutations and a smaller percentage of subclonal driver mutations. A phylogenetic tree showed branching off between T1 and M1/M2. Liquid biopsies, while initially negative, started to detect mutations found among T1, M1 and M2, with increasing allele frequency (AF) throughout time. By the end of 2016, an ESR1 E380Q mutation not previously detected was found at a low AF that had rapidly increased by the end of 2017. RNA-Seq analysis revealed 908 (209 up, 699 down) genes to be differentially expressed between T1 and M1. Most significantly downregulated genes in M1 were TFF1 and PGR indicating loss-of-sensitivity / resistance to AI therapy. Most upregulated genes in M1 included PTHLH, S100P, and SOX2 which have been implicated in promoting tumor growth and metastasis. Pathway analysis between T1 and M1 showed enrichment for the hallmark gene sets ‘Estrogen Response early and late’ and ‘Epithelial mesenchymal transition’.

Conclusions

This phylogenetic reconstruction of the life history of a single patient's cancer revealed an increased mutational rate over time between the patient's primary tumor at diagnosis and samples taken at relapse; and a high level of heterogeneity with metastatic sites sharing only one driver mutation with the primary tumor at diagnosis and having acquired multiple de-novo driver mutations. Transcriptome profiling aided in uncovering the mechanisms that lead to the patient's initial relapse, indicating expression profile changes that promote tumor progression and metastasis as well as reduced sensitivity / resistance to the patients AI based therapy. Monitoring tumor progression through frequent liquid biopsies allowed for the detection of mutations from multiple metastatic sites, in addition to the detection and growth of an AI resistant clone harboring an ESR1 E380Q mutation that was not present in any of the solid tissue samples.

Citation Format: Bing Xu, Anu Amallraja, Padmapriya Swaminathan, Rachel Elsey, Christel Davis, Stephanie Theel, Sarah Viet, Jason Petersen, Amy Krie, Gareth E Davies, Casey B Williams, Erik Ehli, Tobias Meissner. Case report - 16 year life history and genomic evolution of an ER+ HER2- breast cancer [abstract]. In: Proceedings of the 2019 San Antonio Breast Cancer Symposium; 2019 Dec 10-14; San Antonio, TX. Philadelphia (PA): AACR; Cancer Res 2020;80(4 Suppl):Abstract nr P4-05-20.

DNA Methylation Signatures of Breastfeeding in Buccal Cells Collected in Mid-Childhood

Breastfeeding has long-term benefits for children that may be mediated via the epigenome. This pathway has been hypothesized, but the number of empirical studies in humans is small and mostly done by using peripheral blood as the DNA source. We performed an epigenome-wide association study (EWAS) in buccal cells collected around age nine (mean = 9.5) from 1006 twins recruited by the Netherlands Twin Register (NTR). An age-stratified analysis examined if effects attenuate with age (median split at 10 years; n<10 = 517, mean age = 7.9; n>10 = 489, mean age = 11.2). We performed replication analyses in two independent cohorts from the NTR (buccal cells) and the Avon Longitudinal Study of Parents and Children (ALSPAC) (peripheral blood), and we tested loci previously associated with breastfeeding in epigenetic studies. Genome-wide DNA methylation was assessed with the Illumina Infinium MethylationEPIC BeadChip (Illumina, San Diego, CA, USA) in the NTR and with the HumanMethylation450 Bead Chip in the ALSPAC. The duration of breastfeeding was dichotomized ('never' vs. 'ever'). In the total sample, no robustly associated epigenome-wide significant CpGs were identified (α = 6.34 × 10-8). In the sub-group of children younger than 10 years, four significant CpGs were associated with breastfeeding after adjusting for child and maternal characteristics. In children older than 10 years, methylation differences at these CpGs were smaller and non-significant. The findings did not replicate in the NTR sample (n = 98; mean age = 7.5 years), and no nearby sites were associated with breastfeeding in the ALSPAC study (n = 938; mean age = 7.4). Of the CpG sites previously reported in the literature, three were associated with breastfeeding in children younger than 10 years, thus showing that these CpGs are associated with breastfeeding in buccal and blood cells. Our study is the first to show that breastfeeding is associated with epigenetic variation in buccal cells in children. Further studies are needed to investigate if methylation differences at these loci are caused by breastfeeding or by other unmeasured confounders, as well as what mechanism drives changes in associations with age.

Cohabitation is associated with a greater resemblance in gut microbiota which can impact cardiometabolic and inflammatory risk

BACKGROUND

The gut microbiota composition is known to be influenced by a myriad of factors including the host genetic profile and a number of environmental influences. Here, we focus on the environmental influence of cohabitation on the gut microbiota as well as whether these environmentally influenced microorganisms are associated with cardiometabolic and inflammatory burden. We perform this by investigating the gut microbiota composition of various groups of related individuals including cohabitating monozygotic (MZ) twin pairs, non-cohabitating MZ twin pairs and spouse pairs.

RESULTS

A stronger correlation between alpha diversity was found in cohabitating MZ twins (45 pairs, r = 0.64, p = 2.21 × 10^- 06) than in non-cohabitating MZ twin pairs (121 pairs, r = 0.42, p = 1.35 × 10^- 06). Although the correlation of alpha diversity did not attain significance between spouse pairs (42 pairs, r = 0.23, p = 0.15), the correlation was still higher than those in the 209 unrelated pairs (r = - 0.015, p = 0.832). Bray-Curtis (BC) dissimilarity metrics showed cohabitating MZ twin pairs had the most similar gut microbiota communities which were more similar than the BC values of non-cohabitating MZ twins (empirical p-value = 0.0103), cohabitating spouses (empirical p-value = 0.0194), and pairs of unrelated non-cohabitating individuals (empirical p-value< 0.00001). There was also a significant difference between the BC measures from the spouse pairs and those from the unrelated non-cohabitating individuals (empirical p-value< 0.00001). Intraclass correlation coefficients were calculated between the various groups of interest and the results indicate the presence of OTUs with an environmental influence and one OTU that appeared to demonstrate genetic influences. One of the OTUs (Otu0190) was observed to have a significant association with both the cardiometabolic and inflammatory burden scores (p's < 0.05).

CONCLUSIONS

Through the comparison of the microbiota contents of MZ twins with varying cohabitation status and spousal pairs, we showed evidence of environmentally influenced OTUs, one of which had a significant association with cardiometabolic and inflammatory burden scores.

APOBEC mutation drives early-onset squamous cell carcinomas in recessive dystrophic epidermolysis bullosa

Recessive dystrophic epidermolysis bullosa (RDEB) is a rare inherited skin and mucous membrane fragility disorder complicated by early-onset, highly malignant cutaneous squamous cell carcinomas (SCCs). The molecular etiology of RDEB SCC, which arises at sites of sustained tissue damage, is unknown. We performed detailed molecular analysis using whole-exome, whole-genome, and RNA sequencing of 27 RDEB SCC tumors, including multiple tumors from the same patient and multiple regions from five individual tumors. We report that driver mutations were shared with spontaneous, ultraviolet (UV) light-induced cutaneous SCC (UV SCC) and head and neck SCC (HNSCC) and did not explain the early presentation or aggressive nature of RDEB SCC. Instead, endogenous mutation processes associated with apolipoprotein B mRNA-editing enzyme catalytic polypeptide-like (APOBEC) deaminases dominated RDEB SCC. APOBEC mutation signatures were enhanced throughout RDEB SCC tumor evolution, relative to spontaneous UV SCC and HNSCC mutation profiles. Sixty-seven percent of RDEB SCC driver mutations was found to emerge as a result of APOBEC and other endogenous mutational processes previously associated with age, potentially explaining a >1000-fold increased incidence and the early onset of these SCCs. Human papillomavirus-negative basal and mesenchymal subtypes of HNSCC harbored enhanced APOBEC mutational signatures and transcriptomes similar to those of RDEB SCC, suggesting that APOBEC deaminases drive other subtypes of SCC. Collectively, these data establish specific mutagenic mechanisms associated with chronic tissue damage. Our findings reveal a cause for cancers arising at sites of persistent inflammation and identify potential therapeutic avenues to treat RDEB SCC.

Inherited genetic susceptibility to acute lymphoblastic leukemia in Down syndrome

Children with Down syndrome (DS) have a 20-fold increased risk of acute lymphoblastic leukemia (ALL) and distinct somatic features, including CRLF2 rearrangement in ∼50% of cases; however, the role of inherited genetic variation in DS-ALL susceptibility is unknown. We report the first genome-wide association study of DS-ALL, comprising a meta-analysis of 4 independent studies, with 542 DS-ALL cases and 1192 DS controls. We identified 4 susceptibility loci at genome-wide significance: rs58923657 near IKZF1 (odds ratio [OR], 2.02; Pmeta = 5.32 × 10-15), rs3731249 in CDKN2A (OR, 3.63; Pmeta = 3.91 × 10-10), rs7090445 in ARID5B (OR, 1.60; Pmeta = 8.44 × 10-9), and rs3781093 in GATA3 (OR, 1.73; Pmeta = 2.89 × 10-8). We performed DS-ALL vs non-DS ALL case-case analyses, comparing risk allele frequencies at these and other established susceptibility loci (BMI1, PIP4K2A, and CEBPE) and found significant association with DS status for CDKN2A (OR, 1.58; Pmeta = 4.1 × 10-4). This association was maintained in separate regression models, both adjusting for and stratifying on CRLF2 overexpression and other molecular subgroups, indicating an increased penetrance of CDKN2A risk alleles in children with DS. Finally, we investigated functional significance of the IKZF1 risk locus, and demonstrated mapping to a B-cell super-enhancer, and risk allele association with decreased enhancer activity and differential protein binding. IKZF1 knockdown resulted in significantly higher proliferation in DS than non-DS lymphoblastoid cell lines. Our findings demonstrate a higher penetrance of the CDKN2A risk locus in DS and serve as a basis for further biological insights into DS-ALL etiology.

Driver Mutations in Normal Airway Epithelium Elucidate Spatiotemporal Resolution of Lung Cancer

Rationale: Uninvolved normal-appearing airway epithelium has been shown to exhibit specific mutations characteristic of nearby non-small cell lung cancers (NSCLCs). Yet, its somatic mutational landscape in patients with early-stage NSCLC is unknown.Objectives: To comprehensively survey the somatic mutational architecture of the normal airway epithelium in patients with early-stage NSCLC.Methods: Multiregion normal airways, comprising tumor-adjacent small airways, tumor-distant large airways, nasal epithelium and uninvolved normal lung (collectively airway field), matched NSCLCs, and blood cells (n = 498) from 48 patients were interrogated for somatic single-nucleotide variants by deep-targeted DNA sequencing and for chromosomal allelic imbalance events by genome-wide genotype array profiling. Spatiotemporal relationships between the airway field and NSCLCs were assessed by phylogenetic analysis.Measurements and Main Results: Genomic airway field carcinogenesis was observed in 25 cases (52%). The airway field epithelium exhibited a total of 269 somatic mutations in most patients (n = 36) including key drivers that were shared with the NSCLCs. Allele frequencies of these acquired variants were overall higher in NSCLCs. Integrative analysis of single-nucleotide variants and allelic imbalance events revealed driver genes with shared "two-hit" alterations in the airway field (e.g., TP53, KRAS, KEAP1, STK11, and CDKN2A) and those with single hits progressing to two in the NSCLCs (e.g., PIK3CA and NOTCH1).Conclusions: Tumor-adjacent and tumor-distant normal-appearing airway epithelia exhibit somatic driver alterations that undergo selection-driven clonal expansion in NSCLC. These events offer spatiotemporal insights into the development of NSCLC and, thus, potential targets for early treatment.

Testing Familial Transmission of Smoking With Two Different Research Designs

Introduction

Classical twin studies show that smoking is heritable. To determine if shared family environment plays a role in addition to genetic factors, and if they interact (G×E), we use a children-of-twins design. In a second sample, we measure genetic influence with polygenic risk scores (PRS) and environmental influence with a question on exposure to smoking during childhood.

Methods

Data on smoking initiation were available for 723 children of 712 twins from the Netherlands Twin Register (64.9% female, median birth year 1985). Children were grouped in ascending order of risk, based on smoking status and zygosity of their twin-parent and his/her co-twin: never smoking twin-parent with a never smoking co-twin; never smoking twin-parent with a smoking dizygotic co-twin; never smoking twin-parent with a smoking monozygotic co-twin; and smoking twin-parent with a smoking or never smoking co-twin. For 4072 participants from the Netherlands Twin Register (67.3% female, median birth year 1973), PRS for smoking were computed and smoking initiation, smoking heaviness, and exposure to smoking during childhood were available.

Results

Patterns of smoking initiation in the four group children-of-twins design suggested shared familial influences in addition to genetic factors. PRS for ever smoking were associated with smoking initiation in all individuals. PRS for smoking heaviness were associated with smoking heaviness in individuals exposed to smoking during childhood, but not in non-exposed individuals.

Conclusions

Shared family environment influences smoking, over and above genetic factors. Genetic risk of smoking heaviness was only important for individuals exposed to smoking during childhood, versus those not exposed (G×E).

Implications

This study adds to the very few existing children-of-twins (CoT) studies on smoking and combines a CoT design with a second research design that utilizes polygenic risk scores and data on exposure to smoking during childhood. The results show that shared family environment affects smoking behavior over and above genetic factors. There was also evidence for gene-environment interaction (G×E) such that genetic risk of heavy versus light smoking was only important for individuals who were also exposed to (second-hand) smoking during childhood. Together, these findings give additional incentive to recommending parents not to expose their children to cigarette smoking.

The role of genetic variation in DGKK on moderate and severe hypospadias

BACKGROUND

Recent genome-wide association studies of hypospadias have implicated the role of genetic variants in or near the diacylglycerol kinase kappa (DGKK) gene. However, these variants are largely identified among samples of mild and moderate hypospadias cases. Therefore, we evaluated previously identified DGKK variants among second- and third-degree hypospadias cases and controls recruited in Arkansas, a state characterized by a high birth prevalence of hypospadias.

METHODS

Second- and third-degree hypospadias non-Hispanic white cases (n = 36 and n = 9, respectively) and controls (n = 45) were recruited at Arkansas Children's Hospital. Preputial tissue was collected on cases and controls between 2013 and 2017. Cases and controls were genotyped using the Illumina Infinium Global Screening Array. We used logistic regression models to assess the association of genotyped and imputed genetic variants mapped to the DGKK region with second- and third-degree hypospadias.

RESULTS

All families self-reported as non-Hispanic white and genetic principal component analyses did not demonstrate evidence of population stratification. Five DGKK variants previously reported as associated with hypospadias were identified in the genotype data. None of the variants were associated with second- or third-degree hypospadias (range of odds ratios = 0.7-0.9, all p > .05).

CONCLUSIONS

In our analyses, genetic variation in DGKK does not play a role in the development of moderate and severe hypospadias. Our findings provide support to the etiologic heterogeneity of hypospadias by all classifications of severity.

Association studies of up to 1.2 million individuals yield new insights into the genetic etiology of tobacco and alcohol use

Tobacco and alcohol use are leading causes of mortality that influence risk for many complex diseases and disorders1. They are heritable2,3 and etiologically related4,5 behaviors that have been resistant to gene discovery efforts6-11. In sample sizes up to 1.2 million individuals, we discovered 566 genetic variants in 406 loci associated with multiple stages of tobacco use (initiation, cessation, and heaviness) as well as alcohol use, with 150 loci evidencing pleiotropic association. Smoking phenotypes were positively genetically correlated with many health conditions, whereas alcohol use was negatively correlated with these conditions, such that increased genetic risk for alcohol use is associated with lower disease risk. We report evidence for the involvement of many systems in tobacco and alcohol use, including genes involved in nicotinic, dopaminergic, and glutamatergic neurotransmission. The results provide a solid starting point to evaluate the effects of these loci in model organisms and more precise substance use measures.

Dopaminergic Genetic Variants and Voluntary Externally Paced Exercise Behavior

PURPOSE

Most candidate gene studies on the neurobiology of voluntary exercise behavior have focused on the dopaminergic signaling pathway and its role in the mesolimbic reward system. We hypothesized that dopaminergic candidate genes may influence exercise behavior through additional effects on executive functioning and that these effects are only detected when the types of exercise activity are taken into account.

METHODS

Data on voluntary exercise behavior and at least one single-nucleotide polymorphism/variable number of tandem repeat (VNTR) were available for 12,929 participants of the Netherlands Twin Registry. Exercise activity was classified as externally paced if a high level of executive function skill was required. The total volume of voluntary exercise (minutes per week) as well as the volume specifically spent on externally paced activities were tested for association with nine functional dopaminergic polymorphisms (DRD1: rs265981, DRD2/ANKK1: rs1800497, DRD3: rs6280, DRD4: VNTR 48 bp, DRD5: VNTR 130-166 bp, DBH: rs2519152, DAT1: VNTR 40 bp, COMT: rs4680, MAOA: VNTR 30 bp), a polygenic score (PGS) based on nine alleles leading to lower dopamine responsiveness, and a PGS based on three alleles associated with both higher reward sensitivity and better executive functioning (DRD2/ANKK1: "G" allele, COMT: Met allele, DAT1: 440-bp allele).

RESULTS

No association with total exercise volume or externally paced exercise volume was found for individual alleles or the nine-allele PGS. The volume of externally paced exercise behavior was significantly associated with the reward and executive function congruent PGS. This association was driven by the DAT1 440-bp and COMT Met allele, which acted as increaser alleles for externally paced exercise behavior.

CONCLUSIONS

Taking into account the types of exercise activity may increase the success of identifying genetic variants and unraveling the neurobiology of voluntary exercise behavior.

The moderating role of SES on genetic differences in educational achievement in the Netherlands

Parental socioeconomic status (SES) is a strong predictor of children's educational achievement (EA), with an increasing effect throughout development. Inequality in educational outcomes between children from different SES backgrounds exists in all Western countries. It has been proposed that a cause of this inequality lies in the interplay between genetic effects and SES on EA, which might depend on society and the equality of the education system. This study adopted two approaches, a classical twin design and polygenic score (PGS) approach, to address the effect of parental SES on EA in a large sample of 12-year-old Dutch twin pairs (2479 MZ and 4450 DZ twin pairs with PGSs for educational attainment available in 2335 children) from the Netherlands Twin Register (NTR). The findings of this study indicated that average EA increased with increasing parental SES. The difference in EA between boys and girls became smaller in the higher SES groups. The classical twin design analyses based on genetic covariance structure modeling pointed to lower genetic, environmental, and thus phenotypic variation in EA at higher SES. Independent from a child's PGS, parental SES predicted EA. However, the strength of the association between PGS and EA did not depend on parental SES. In a within-family design, the twin with a higher PGS scored higher on EA than the co-twin, demonstrating that the effect of the PGS on EA was at least partly independent from parental SES. To conclude, EA depended on SES both directly and indirectly, and SES moderated the additive genetic and environmental components of EA. Adding information from PGS, in addition to parental SES, improved the prediction of children's EA.

Genome-wide association meta-analysis of age at first cannabis use

BACKGROUND AND AIMS

Cannabis is one of the most commonly used substances among adolescents and young adults. Earlier age at cannabis initiation is linked to adverse life outcomes, including multi-substance use and dependence. This study estimated the heritability of age at first cannabis use and identified associations with genetic variants.

METHODS

A twin-based heritability analysis using 8055 twins from three cohorts was performed. We then carried out a genome-wide association meta-analysis of age at first cannabis use in a discovery sample of 24 953 individuals from nine European, North American and Australian cohorts, and a replication sample of 3735 individuals.

RESULTS

The twin-based heritability for age at first cannabis use was 38% [95% confidence interval (CI) = 19-60%]. Shared and unique environmental factors explained 39% (95% CI = 20-56%) and 22% (95% CI = 16-29%). The genome-wide association meta-analysis identified five single nucleotide polymorphisms (SNPs) on chromosome 16 within the calcium-transporting ATPase gene (ATP2C2) at P < 5E-08. All five SNPs are in high linkage disequilibrium (LD) (r2  > 0.8), with the strongest association at the intronic variant rs1574587 (P = 4.09E-09). Gene-based tests of association identified the ATP2C2 gene on 16q24.1 (P = 1.33e-06). Although the five SNPs and ATP2C2 did not replicate, ATP2C2 has been associated with cocaine dependence in a previous study. ATP2B2, which is a member of the same calcium signalling pathway, has been associated previously with opioid dependence. SNP-based heritability for age at first cannabis use was non-significant.

CONCLUSION

Age at cannabis initiation appears to be moderately heritable in western countries, and individual differences in onset can be explained by separate but correlated genetic liabilities. The significant association between age of initiation and ATP2C2 is consistent with the role of calcium signalling mechanisms in substance use disorders.

Genome-wide analysis of DNA methylation in buccal cells: a study of monozygotic twins and mQTLs

BACKGROUND

DNA methylation arrays are widely used in epigenome-wide association studies and methylation quantitative trait locus (mQTL) studies. Here, we performed the first genome-wide analysis of monozygotic (MZ) twin correlations and mQTLs on data obtained with the Illumina MethylationEPIC BeadChip (EPIC array) and compared the performance of the EPIC array to the Illumina HumanMethylation450 BeadChip (HM450 array) for buccal-derived DNA.

RESULTS

Good-quality EPIC data were obtained for 102 buccal-derived DNA samples from 49 MZ twin pairs (mean age = 7.5 years, range = 1-10). Differences between MZ twins in the cellular content of buccal swabs were a major driver for differences in their DNA methylation profiles, highlighting the importance to adjust for cellular composition in DNA methylation studies of buccal-derived DNA. After adjusting for cellular composition, the genome-wide mean correlation (r) between MZ twins was 0.21 for the EPIC array, and cis mQTL analysis in 84 twins identified 1,296,323 significant associations (FDR 5%), encompassing 33,749 methylation sites and 616,029 genetic variants. MZ twin correlations were slightly larger (p < 2.2 × 10-16) for novel EPIC probes (N = 383,066, mean r = 0.22) compared to probes that are also present on HM450 (N = 406,822, mean r = 0.20). In line with this observation, a larger percentage of novel EPIC probes was associated with genetic variants (novel EPIC probes with significant mQTL 4.7%, HM450 probes with mQTL 3.9%, p < 2.2 × 10-16). Methylation sites with a large MZ correlation and sites associated with mQTLs were most strongly enriched in epithelial cell DNase I hypersensitive sites (DHSs), enhancers, and histone mark H3K4me3.

CONCLUSIONS

We conclude that the contribution of familial factors to individual differences in DNA methylation and the effect of mQTLs are larger for novel EPIC probes, especially those within regulatory elements connected to active regions specific to the investigated tissue.

Associations between loneliness and personality are mostly driven by a genetic association with Neuroticism

OBJECTIVE

Loneliness is an aversive response to a discrepancy between desired and actual social relationships and correlates with personality. We investigate the relationship of loneliness and personality in twin family and molecular genetic data.

METHOD

Phenotypic correlations between loneliness and the Big Five personality traits were estimated in 29,625 adults, and in a group with genome-wide genotype data (N = 4,222), genetic correlations were obtained. We explored whether genetic correlations may reflect causal relationships by investigating within monozygotic twin pair differences (N_pairs  = 2,662), by longitudinal within-subject changes in personality and loneliness (N = 4,260-9,238 longitudinal comparisons), and by longitudinal cross-lagged panel analyses (N = 15,628). Finally, we tested whether genetic correlations were due to cross-trait assortative mating (N_{spouse pairs}  = 4,436).

RESULTS

The strongest correlations with loneliness were observed for Neuroticism (r = .55) and Extraversion (r = -.33). Only Neuroticism showed a high correlation with loneliness independent of other personality traits (r = .50), so follow-up analyses focused on Neuroticism. The genetic correlation between loneliness and Neuroticism from genotyped variants was .71; a significant reciprocal causal relationship and nonsignificant cross-trait assortative mating imply that this is at least partly due to mediated pleiotropy.

CONCLUSIONS

We show that the relationship between loneliness and personality is largely explained by its relationship with Neuroticism, which is substantially genetic in nature.

Abstract 222: Genome-wide association study of acute lymphoblastic leukemia in children with Down syndrome

Purpose: Children with Down syndrome (DS) have a 20-fold increased risk of acute lymphoblastic leukemia (ALL) compared to children without DS. While genome-wide association studies (GWAS) have identified several susceptibility loci in childhood ALL, studies of ALL in children with DS are lacking. Therefore, we conducted the first GWAS of DS-ALL.

Methods: We analyzed independent cohorts of: 1) 226 newly diagnosed DS-ALL cases from Children's Oncology Group (COG) ALL trials (2000-2013) and 436 DS controls from the National Down Syndrome Project (NDSP), 2) 124 additional COG ALL cases (2011-2015) and 336 additional NDSP DS controls, 3) 20 DS-ALL cases and 275 DS controls from Michigan neonatal bloodspots, and 4) 157 DS-ALL cases and 145 DS controls largely from neonatal bloodspots from California and Washington. Genotyping was performed with Affymetrix or Illumina single nucleotide polymorphism (SNP) arrays. STRUCTURE software was used to define European (372 cases, 1,056 controls), Hispanic (140 cases, 136 controls), and African (15 cases, 62 controls) genetic ancestry. After genome-wide imputation and quality control, ancestry- and cohort-specific associations were evaluated at &gt;6,000,000 autosomal SNPs with minor allele frequency ≥1%. Associations were meta-analyzed across cohort and ancestry groups, assuming additive allelic effects.

Results: Genome-wide significant (p&lt;5x10-8) association signals were identified for known ALL susceptibility loci, including rs58923657 near IKZF1 (Odds Ratio [OR]=2.02, p=5.32x10-15), CDKN2A missense mutation rs3731249 (OR=3.63, p=3.91x10-10), rs3781093 near GATA3 (OR=1.73, p=2.89x10-8), and rs7090445 near ARID5B (OR=1.57, p=2.93x10-8). A novel potential risk locus was identified at chromosome 20q11.21 (rs78019519, OR=3.17, p=5.11x10-7) with consistent effects observed across each cohort and ancestry group. This SNP is in the promoter region of the oncogene TPX2 and is also associated with expression of HM13 in whole blood in the Genotype-Tissue Expression (GTEx) Portal.

Conclusion: We confirmed that known ALL susceptibility loci in children without DS, including IKZF1, CDKN2A, GATA3, PIP4K2A and ARID5B, also confer risk of ALL in children with DS, with CDKN2A showing the largest effect size. We also identified a potentially novel locus associated with ALL susceptibility in DS at chromosome 20q11.21. Additional investigation of these loci is ongoing and may advance our understanding of DS-ALL etiology and biology.

Citation Format: Austin L. Brown, Adam J. de Smith, Michael E. Scheurer, Noah A. Kallsen, Shanna A. Peyton, Gareth E. Davies, Erik A. Ehli, Michael E. Zwick, Naomi Winick, Kelly Maloney, Anne L. Angiolillo, Reuven Schore, MIchael M. Burke, Wanda L. Salzer, Nyla A. Heerema, Andrew J. Carroll, Michael J. Borowitz, Brent L. Wood, William L. Carroll, Elizabeth A. Raetz, Elanor Feingold, Stephanie L. Sherman, Wenjian Yang, Meenakshi Devidas, Kyle Walsh, Andrew T. DeWan, Maria S. Pombo-de-Oliveira, Jeffrey W. Taub, Daniel Sinnett, Jasmine Healy, Jillian M. Birch, Lisa F. Barcellos, Helen Hansen, Ivan Smirnov, Charles G. Mullighan, Stephen P. Hunger, Ching-Hon Pui, Mignon Loh, Joe L. Wiemels, Xiaomei Ma, Catherine Metayer, Beth A. Mueller, Mary V. Relling, Jun J. Yang, Philip J. Lupo, Karen R. Rabin. Genome-wide association study of acute lymphoblastic leukemia in children with Down syndrome [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2018; 2018 Apr 14-18; Chicago, IL. Philadelphia (PA): AACR; Cancer Res 2018;78(13 Suppl):Abstract nr 222.

Predicting loneliness with polygenic scores of social, psychological and psychiatric traits

Loneliness is a heritable trait that accompanies multiple disorders. The association between loneliness and mental health indices may partly be due to inherited biological factors. We constructed polygenic scores for 27 traits related to behavior, cognition and mental health and tested their prediction for self-reported loneliness in a population-based sample of 8798 Dutch individuals. Polygenic scores for major depressive disorder (MDD), schizophrenia and bipolar disorder were significantly associated with loneliness. Of the Big Five personality dimensions, polygenic scores for neuroticism and conscientiousness also significantly predicted loneliness, as did the polygenic scores for subjective well-being, tiredness and self-rated health. When including all polygenic scores simultaneously into one model, only 2 major depression polygenic scores remained as significant predictors of loneliness. When controlling only for these 2 MDD polygenic scores, only neuroticism and schizophrenia remain significant. The total variation explained by all polygenic scores collectively was 1.7%. The association between the propensity to feel lonely and the susceptibility to psychiatric disorders thus pointed to a shared genetic etiology. The predictive power of polygenic scores will increase as the power of the genome-wide association studies on which they are based increases and may lead to clinically useful polygenic scores that can inform on the genetic predisposition to loneliness and mental health.

An Extended Twin-Pedigree Study of Neuroticism in the Netherlands Twin Register

For the participants in the Netherlands Twin Register (NTR) we constructed the extended pedigrees which specify all relations among nuclear and larger twin families in the register. A total of 253,015 subjects from 58,645 families were linked to each other, to the degree that we had information on the relations among participants. We describe the algorithm that was applied to construct the pedigrees. For > 30,000 adolescent and adult NTR participants data were available on harmonized neuroticism scores. We analyzed these data in the Mendel software package (Lange et al., Bioinformatics 29(12):1568-1570, 2013) to estimate the contributions of additive and non-additive genetic factors. In contrast to much of the earlier work based on twin data rather than on extended pedigrees, we could also estimate the contribution of shared household effects in the presence of non-additive genetic factors. The estimated broad-sense heritability of neuroticism was 47%, with almost equal contributions of additive and non-additive (dominance) genetic factors. A shared household effect explained 13% and unique environmental factors explained the remaining 40% of the variance in neuroticism.

Abstract 2463: Tumor profiling of separated carcinomatous and sarcomatous components from uterine carcinosarcoma biopsies provides insights into their development

Uterine carcinosarcoma (UCS) is a rare and aggressive form of uterine cancer. It is bi-phasic, exhibiting histological features of both malignant epithelial (carcinomatous) and mesenchymal (sarcomatous) elements. Studies have indicated that UCS arises from sarcomatous differentiation of high-grade carcinoma while others have suggested a bi-clonal nature. Given these differences, we sought to separate the carcinoma and sarcoma elements of UCS to try to understand their molecular differences and gain further insights into how these tumors develop. We macrodissected carcinomatous, sarcomatous, and normal cells from formalin fixed paraffin embedded (FFPE) uterine samples of 10 UCS patients. DNA and RNA were isolated and extracted using the Qiagen AllPrep DNA/RNA FFPE kit. Whole-genome SNP microarrays and deep sequencing of 26 cancer genes was performed, using the Illumina Infinium OmniExpressExome array and the TruSight Tumor panel respectively. Illumina HiSeq mRNA sequencing was also performed to quantify gene expression. The genomic allelic imbalance (AI) profiling, called from the SNP data by hapLOH, showed that sarcoma samples were more aberrant than their carcinoma counterparts (abstract 131, AACR 2016). From the targeted sequencing, the Illumina Amplicon-DS Somatic Variant Caller was employed to call somatic mutations. Mutations were identified in TP53 in both the sarcoma and carcinoma samples of all 10 patients. Frequently mutated genes included APC, EGFR, MET and MSH6 which were found in 60-80% of the patients. Genes mutated in less than 50% of the patients included PTEN, KRAS, KIT, FBXW7, PIK3CA, FGFR2, and CTNNB1. Current results showed no association of a mutated gene to either the sarcoma or carcinoma component of UCS. RSEM, STAR and EBSeq were applied to the RNA-seq data for gene expression quantification. Approximately 2500 genes were identified as being differentially expressed (DE) between normal and carcinoma samples. Just over 4000 genes were identified as being DE between normal and sarcoma samples. 75% of the DE genes in the carcinoma were also identified in the sarcoma. Using DAVID functional annotation tool, we characterized these gene sets with KEGG pathways. Deregulated pathways identified in both carcinoma and sarcoma include: cell cycle, transcriptional regulation, Ras and p53 signaling. Some additional pathways are putatively associated with sarcoma only, including MAPK and PI3K-Akt signaling. We report here the differences between sarcoma and carcinoma components of UCS from multiple molecular perspectives. From the genomic AI and DE analysis, the carcinoma aberrations appear to be mostly a subset of the sarcoma tumor profiles, where sarcoma samples appear to be more highly aberrant compared to the carcinoma samples. One possible inference from this is that the sarcoma originated and evolved from the carcinoma cells.

Citation Format: Yihua Liu, Zachary Weber, F. Anthony San Lucas, Aditya Deshpande, Raed Sulaiman, Mary Fagerness, Natasha Flier, Joseph Sulaiman, Christel M. Davis, Jerry Fowler, Gareth E. Davies, David Starks, Luis Rojas-Espaillat, Paul Scheet, Erik A. Ehli. Tumor profiling of separated carcinomatous and sarcomatous components from uterine carcinosarcoma biopsies provides insights into their development [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2017; 2017 Apr 1-5; Washington, DC. Philadelphia (PA): AACR; Cancer Res 2017;77(13 Suppl):Abstract nr 2463. doi:10.1158/1538-7445.AM2017-2463

Mutational Heterogeneity in APC and KRAS Arises at the Crypt Level and Leads to Polyclonality in Early Colorectal Tumorigenesis

Purpose: The majority of genomic alterations causing intratumoral heterogeneity (ITH) in colorectal cancer are thought to arise during early stages of carcinogenesis as a burst but only after truncal mutations in APC have expanded a single founder clone. We have investigated if the initial source of ITH is consequent to multiple independent lineages derived from different crypts harboring distinct truncal APC and driver KRAS mutations, thus challenging the prevailing monoclonal monocryptal model.Experimental Design: High-depth next-generation sequencing and SNP arrays were performed in whole-lesion extracts of 37 familial adenomatous polyposis colorectal adenomas. Also, ultrasensitive genotyping of hotspot mutations of APC and KRAS was performed using nanofluidic PCRs in matched bulk biopsies (n = 59) and crypts (n = 591) from 18 adenomas and seven carcinomas and adjacent normal tissues.Results: Multiple co-occurring truncal APC and driver KRAS alterations were uncovered in whole-lesion extracts from adenomas and subsequently confirmed to belong to multiple clones. Ultrasensitive genotyping of bulk biopsies and crypts revealed novel undetected APC mutations that were prominent among carcinomas, whereas abundant wild-type APC crypts were detected in adenomas. KRAS mutational heterogeneity within crypts was evident in both adenomas and carcinomas with a higher degree of concordance between biopsy and crypt genotyping in carcinomas. Nonrandom heterogeneity among crypts was also observed.Conclusions: The striking degree of nonrandom intercrypt heterogeneity in truncal and driver gene mutations observed in adenomas and carcinomas is consistent with a polycryptal model derived from multiple independent initiation linages as the source of early ITH in colorectal carcinogenesis. Clin Cancer Res; 23(19); 5936-47. ©2017 AACR.

Abstract A02: High-depth sequencing reveals the presence of multi-clonality in colorectal premalignancy

Clonality in carcinomas is thought to arise from the accumulation of mutations and changes in selective pressure during carcinogenesis and further stages of tumor development. However, the identification of clonality through different molecular and systems biology tools has challenged the current understanding of carcinogenesis and progression. Familial Adenomatous Polyposis (FAP) is a genetic condition secondary to germline mutations in APC and characterized by the development of hundreds of polyps in the colorectum. Therefore, FAP is an ideal model to study the presence of clonality during early stages of colorectal carcinogenesis in humans. Second hits in APC are a valuable and intuitive method to establish the presence of clonality in FAP, as additional somatic APC mutations will reflect the presence of additional clones. Here, we present a targeted high-depth sequencing analysis of APC and KRAS in conjunction with high density SNP array data in 37 colorectal adenomas and matched adjacent uninvolved mucosa samples. Ampliseq sequencing was performed using the next-generation sequencing platform IT Personal Genome Machine (PGM; Life Technologies) using the Ion PGM 200 Sequencing Kit on an Ion 318 Chip Kit (Life Technologies, average of 70% bases at 500x). Ion Torrent Variant Caller v4.2 was run in the somatic low stringency proton mode to detect variants in APC and KRAS against hg19. Then, normal variants were subtracted from matched adenoma variants to create a list of somatic candidates for each adenoma and normal pair. SNP array analysis was performed using the HumanOmniExpressExome-8 v1.2 BeadChip (Illumina, San Diego, CA). To validate the presence of double somatic hits in APC we performed a colony genotyping assay from one adenoma sample. First, we identified APC somatic mutations or loss of heterogeneity (LOH) in 5q (second hits) in 81% (30/37) and KRAS somatic mutations in 35% (16/37) of adenomas analyzed. Six adenomas (16%) presented double somatic APC mutations with different allelic frequencies with three of them showing a significant difference in allelic frequencies among the somatic APC hits (&gt;30%); in addition, one adenoma presented two KRAS mutations. The APC clonality was validated in one adenoma using colony analysis. The region that contains the two somatic mutations (c.4348C&gt;T and c.4267_4280del) and the germline mutation (c.2894delA) was amplified by Long-Range PCR and cloned into the pGEM-T plasmid. DNA from one hundred colonies was extracted and genotyped using Taqman assay probes. All samples were genotyped in triplicate and positive and negative controls for all mutations were included in every single plate. A total of 25% of the colonies harbored a wild-type allele, 36% the known germline APC mutation, 15% the APC somatic c.4348C&gt;T mutation and 24% the APC somatic c.4267_4280del mutation. In conclusion, our analysis showed that there are multiple mutated APC subclones driving initiation of colorectal carcinogenesis, thus supporting experimentally a finding of multi-clonality in premalignant colonic lesions. These findings support a model of colorectal carcinogenesis that is based on polycryptal polyclonal adenomas and calls for revisiting the “Big-Bang” model based on the presence of multiple initiator clones. Further analyses aimed to understand the functional impact and the hierarchical relations among the different identified subclones are warranted.

Citation Format: Ester Borras, Kyle Chang, Anthony San Lucas, Gita Masand, Axel Delgado Amador, Jerry Fowler, Gareth E. Davies, Erik A. Ehli, Y. Nancy You, Patrick M. Lynch, Ernest T. Hawk, Gabriel Capella, Paul A. Scheet, Eduardo Vilar. High-depth sequencing reveals the presence of multi-clonality in colorectal premalignancy. [abstract]. In: Proceedings of the AACR Special Conference on Colorectal Cancer: From Initiation to Outcomes; 2016 Sep 17-20; Tampa, FL. Philadelphia (PA): AACR; Cancer Res 2017;77(3 Suppl):Abstract nr A02.

Novel genetic loci associated with hippocampal volume

The hippocampal formation is a brain structure integrally involved in episodic memory, spatial navigation, cognition and stress responsiveness. Structural abnormalities in hippocampal volume and shape are found in several common neuropsychiatric disorders. To identify the genetic underpinnings of hippocampal structure here we perform a genome-wide association study (GWAS) of 33,536 individuals and discover six independent loci significantly associated with hippocampal volume, four of them novel. Of the novel loci, three lie within genes (ASTN2, DPP4 and MAST4) and one is found 200 kb upstream of SHH. A hippocampal subfield analysis shows that a locus within the MSRB3 gene shows evidence of a localized effect along the dentate gyrus, subiculum, CA1 and fissure. Further, we show that genetic variants associated with decreased hippocampal volume are also associated with increased risk for Alzheimer's disease (rg=-0.155). Our findings suggest novel biological pathways through which human genetic variation influences hippocampal volume and risk for neuropsychiatric illness.

Novel genetic loci underlying human intracranial volume identified through genome-wide association

Intracranial volume reflects the maximally attained brain size during development, and remains stable with loss of tissue in late life. It is highly heritable, but the underlying genes remain largely undetermined. In a genome-wide association study of 32,438 adults, we discovered five previously unknown loci for intracranial volume and confirmed two known signals. Four of the loci were also associated with adult human stature, but these remained associated with intracranial volume after adjusting for height. We found a high genetic correlation with child head circumference (ρgenetic = 0.748), which indicates a similar genetic background and allowed us to identify four additional loci through meta-analysis (Ncombined = 37,345). Variants for intracranial volume were also related to childhood and adult cognitive function, and Parkinson's disease, and were enriched near genes involved in growth pathways, including PI3K-AKT signaling. These findings identify the biological underpinnings of intracranial volume and their link to physiological and pathological traits.

A Genome-Wide Association Meta-Analysis of Attention-Deficit/Hyperactivity Disorder Symptoms in Population-Based Pediatric Cohorts

OBJECTIVE

The aims of this study were to elucidate the influence of common genetic variants on childhood attention-deficit/hyperactivity disorder (ADHD) symptoms, to identify genetic variants that explain its high heritability, and to investigate the genetic overlap of ADHD symptom scores with ADHD diagnosis.

METHOD

Within the EArly Genetics and Lifecourse Epidemiology (EAGLE) consortium, genome-wide single nucleotide polymorphisms (SNPs) and ADHD symptom scores were available for 17,666 children (<13 years of age) from nine population-based cohorts. SNP-based heritability was estimated in data from the three largest cohorts. Meta-analysis based on genome-wide association (GWA) analyses with SNPs was followed by gene-based association tests, and the overlap in results with a meta-analysis in the Psychiatric Genomics Consortium (PGC) case-control ADHD study was investigated.

RESULTS

SNP-based heritability ranged from 5% to 34%, indicating that variation in common genetic variants influences ADHD symptom scores. The meta-analysis did not detect genome-wide significant SNPs, but three genes, lying close to each other with SNPs in high linkage disequilibrium (LD), showed a gene-wide significant association (p values between 1.46 × 10(-6) and 2.66 × 10(-6)). One gene, WASL, is involved in neuronal development. Both SNP- and gene-based analyses indicated overlap with the PGC meta-analysis results with the genetic correlation estimated at 0.96.

CONCLUSION

The SNP-based heritability for ADHD symptom scores indicates a polygenic architecture, and genes involved in neurite outgrowth are possibly involved. Continuous and dichotomous measures of ADHD appear to assess a genetically common phenotype. A next step is to combine data from population-based and case-control cohorts in genetic association studies to increase sample size and to improve statistical power for identifying genetic variants.

Abstract 131: Allelic imbalance analysis of uterine carcinosarcoma: An inquiry into the dual nature of the neoplasm

Uterine carcinosarcoma (UCS), also known as Malignant Mixed Mullerian Tumor, is a neoplasm of the female genital tract that shows histological features of both carcinoma and sarcoma. Some studies have indicated that UCS arises from sarcomatous differentiation of high-grade carcinoma while others have confirmed a bi-clonal nature of the tumor. Given these differences, further investigation in the origin of this tumor with newer approaches and technologies is warranted. In this study we determined the allelic imbalance (AI) profiles of the two components of UCS and compared the AI events that were shared by, or differed between, them. Samples were obtained from 10 patients diagnosed with UCS and were preserved in formalin fixed paraffin embedded (FFPE) blocks. The two components (carcinoma & sarcoma) were identified and micro-dissected, and whole-genome DNA micro-arrays were applied to the isolated DNA. We then applied a sensitive computational method that combines haplotype information with SNP array data to identify somatic segmental copy number variations and copy-neutral loss of heterozygosity (cnLOH). After we obtained the AI events we filtered small events (&lt;2 MB) to rule out germline duplications. We then obtained events that are common to both components by using a criterion of 80% reciprocal overlap. All such events were classified as a gain, loss or cnLOH. We found that, on an average, 80% of the events found in the carcinoma component are common with the sarcoma component. Conversely, only 64% of events seen in the sarcoma component are found in the carcinoma. This may be in part due to the fact that the sarcoma component usually showed greater number of events and most of these events were larger than the corresponding carcinoma component. Also the sarcoma component showed an increased degree of allelic imbalance than the counterpart. We then looked into particular regions of known tumor suppressor and oncogenes and identified loss in the TP53 gene region (60% of tumors), gains in PIK3CA and PTEN gene regions (40% of tumors), among other events. Our analysis showed that the sarcomatous component shared major portions of the AI profile with the carcinomatous component, but also showed additional events and a greater degree of AI. One explanation for this observation can be that the sarcomatous component retains the AI profile of the carcinomatous component and also gathers additional imbalance. This suggests directionality to the development of UCS: that it arises by sarcomatous differentiation of already existing carcinoma of the uterus. We are pursuing additional assessments of point mutations and expression profiles to complement our existing analyses.

Citation Format: Aditya S. Deshpande, Zachary Weber, Raed Sulaiman, Natasha Flier, Cheryl Ageton, Mary Fagerness, Joseph Sulaiman, Gareth E. Davies, David Starks, Luis Rojas-Espaillat, Paul A. Scheet, Erik Ehli. Allelic imbalance analysis of uterine carcinosarcoma: An inquiry into the dual nature of the neoplasm. [abstract]. In: Proceedings of the 107th Annual Meeting of the American Association for Cancer Research; 2016 Apr 16-20; New Orleans, LA. Philadelphia (PA): AACR; Cancer Res 2016;76(14 Suppl):Abstract nr 131.

A genome-wide approach to children's aggressive behavior: The EAGLE consortium

Individual differences in aggressive behavior emerge in early childhood and predict persisting behavioral problems and disorders. Studies of antisocial and severe aggression in adulthood indicate substantial underlying biology. However, little attention has been given to genome-wide approaches of aggressive behavior in children. We analyzed data from nine population-based studies and assessed aggressive behavior using well-validated parent-reported questionnaires. This is the largest sample exploring children's aggressive behavior to date (N = 18,988), with measures in two developmental stages (N = 15,668 early childhood and N = 16,311 middle childhood/early adolescence). First, we estimated the additive genetic variance of children's aggressive behavior based on genome-wide SNP information, using genome-wide complex trait analysis (GCTA). Second, genetic associations within each study were assessed using a quasi-Poisson regression approach, capturing the highly right-skewed distribution of aggressive behavior. Third, we performed meta-analyses of genome-wide associations for both the total age-mixed sample and the two developmental stages. Finally, we performed a gene-based test using the summary statistics of the total sample. GCTA quantified variance tagged by common SNPs (10-54%). The meta-analysis of the total sample identified one region in chromosome 2 (2p12) at near genome-wide significance (top SNP rs11126630, P = 5.30 × 10(-8) ). The separate meta-analyses of the two developmental stages revealed suggestive evidence of association at the same locus. The gene-based analysis indicated association of variation within AVPR1A with aggressive behavior. We conclude that common variants at 2p12 show suggestive evidence for association with childhood aggression. Replication of these initial findings is needed, and further studies should clarify its biological meaning. © 2015 Wiley Periodicals, Inc.

Genomic Landscape of Colorectal Mucosa and Adenomas

The molecular basis of the adenoma-to-carcinoma transition has been deduced using comparative analysis of genetic alterations observed through the sequential steps of intestinal carcinogenesis. However, comprehensive genomic analyses of adenomas and at-risk mucosa are still lacking. Therefore, our aim was to characterize the genomic landscape of colonic at-risk mucosa and adenomas. We analyzed the mutation profile and copy number changes of 25 adenomas and adjacent mucosa from 12 familial adenomatous polyposis patients using whole-exome sequencing and validated allelic imbalances (AI) in 37 adenomas using SNP arrays. We assessed for evidence of clonality and performed estimations on the proportions of driver and passenger mutations using a systems biology approach. Adenomas had lower mutational rates than did colorectal cancers and showed recurrent alterations in known cancer driver genes (APC, KRAS, FBXW7, TCF7L2) and AIs in chromosomes 5, 7, and 13. Moreover, 80% of adenomas had somatic alterations in WNT pathway genes. Adenomas displayed evidence of multiclonality similar to stage I carcinomas. Strong correlations between mutational rate and patient age were observed in at-risk mucosa and adenomas. Our data indicate that at least 23% of somatic mutations are present in at-risk mucosa prior to adenoma initiation. The genomic profiles of at-risk mucosa and adenomas illustrate the evolution from normal tissue to carcinoma via greater resolution of molecular changes at the inflection point of premalignant lesions. Furthermore, substantial genomic variation exists in at-risk mucosa before adenoma formation, and deregulation of the WNT pathway is required to foster carcinogenesis. Cancer Prev Res; 9(6); 417-27. ©2016 AACR.

Pre-hospital pelvic girdle injury: Improving diagnostic accuracy in a physician-led trauma service

BACKGROUND

Examination of missed injuries in our physician-led pre-hospital trauma service indicated that the significant injuries missed were often pelvic fractures. We therefore conducted a study whose aim was to evaluate the pre-hospital diagnostic accuracy of pelvic girdle injuries, and how this would be affected by implementing the pelvic injury treatment guidelines recently published by the Faculty of Pre-Hospital Care.

STUDY DESIGN

All blunt trauma patients attended in a 5-month period were included in the study. The presence or absence of pelvic girdle injury on computed tomography (CT) or, if unavailable, pelvic X-ray was used as a primary outcome measure. A retrospective database and case note review was conducted to identify patients who had pelvic binder applied in the study period. For the purposes of the study, pelvic ring and acetabular fractures were grouped together as patients with suspected pelvic girdle injury that should be fitted with a pelvic binder in the pre-hospital setting. The sensitivity and specificity, relating to the presence of pelvic girdle injury in patients with pelvic binders, was calculated in order to determine pre-hospital diagnostic accuracy.

RESULTS

785 patients were attended during the study period. 170 met the study inclusion criteria. 26 (15.3%) sustained a pelvic girdle injury. 45 (26.5%) had a pelvic binder applied. There were eight missed fractures (31%), of which the majority (six) sustained less severe injuries that were managed non-operatively. Two patients required operative fixation. Radiological images and/or reports were available on 169 (99.4%) patients. As a test of the presence of pelvic fracture, pelvic binder application had a sensitivity of 0.69 (95% CI 0.50-0.85) and a specificity of 0.81 (95% CI 0.74-0.87).

CONCLUSIONS

Even with a careful clinical assessment and a low threshold for binder application, this study highlights the problems of distracting injury when trying to diagnose and manage pelvic fractures. By implementing the pelvic treatment guidelines published by the Faculty of Pre-hospital Care, the missed injury rate could be reduced from 31% to 8%.

Minimally invasive genomic and transcriptomic profiling of visceral cancers by next-generation sequencing of circulating exosomes

BACKGROUND

The ability to perform comprehensive profiling of cancers at high resolution is essential for precision medicine. Liquid biopsies using shed exosomes provide high-quality nucleic acids to obtain molecular characterization, which may be especially useful for visceral cancers that are not amenable to routine biopsies.

PATIENTS AND METHODS

We isolated shed exosomes in biofluids from three patients with pancreaticobiliary cancers (two pancreatic, one ampullary). We performed comprehensive profiling of exoDNA and exoRNA by whole genome, exome and transcriptome sequencing using the Illumina HiSeq 2500 sequencer. We assessed the feasibility of calling copy number events, detecting mutational signatures and identifying potentially actionable mutations in exoDNA sequencing data, as well as expressed point mutations and gene fusions in exoRNA sequencing data.

RESULTS

Whole-exome sequencing resulted in 95%-99% of the target regions covered at a mean depth of 133-490×. Genome-wide copy number profiles, and high estimates of tumor fractions (ranging from 56% to 82%), suggest robust representation of the tumor DNA within the shed exosomal compartment. Multiple actionable mutations, including alterations in NOTCH1 and BRCA2, were found in patient exoDNA samples. Further, RNA sequencing of shed exosomes identified the presence of expressed fusion genes, representing an avenue for elucidation of tumor neoantigens.

CONCLUSIONS

We have demonstrated high-resolution profiling of the genomic and transcriptomic landscapes of visceral cancers. A wide range of cancer-derived biomarkers could be detected within the nucleic acid cargo of shed exosomes, including copy number profiles, point mutations, insertions, deletions, gene fusions and mutational signatures. Liquid biopsies using shed exosomes has the potential to be used as a clinical tool for cancer diagnosis, therapeutic stratification and treatment monitoring, precluding the need for direct tumor sampling.

Synthesis of new thiazolo[4,5-d]pyrimidines as Corticotropin releasing factor modulators

Corticotropin-releasing factor (CRF) is a neurohormone that plays a crucial role in integrating the body's overall response to stress. It appears necessary and sufficient for the organism to mount functional, physiological and endocrine responses to stressors. CRF is released in response to various triggers such as chronic stress. The role of CRF and its involvement in these neurological disorders suggest that new drugs that can target the CRF function or bind to its receptors may represent a new development of neuropsychiatric medicines to treat various stress-related disorders including depression, anxiety and addictive disorders. Based on pharmacophore of the CRF1 receptor antagonists, a new series of thiazolo[4,5-d] pyrimidines were synthesized as Corticotropin-releasing factor (CRF) receptor modulators and the prepared compounds carry groups shown to produce optimum binding affinity to CRF receptors. Twenty two compounds were evaluated for their CRF1 receptor binding affinity in HEK 293 cell lines and two compounds 5o and 5s showed approximately 25% binding affinity to CRF1 receptors. Selected compounds (5c and 5f) were also evaluated for their effect on expression of genes associated with depression and anxiety disorders such as CRF1, CREB1, MAO-A, SERT, NPY, DatSLC6a3, and DBH and significant upregulation of CRF1 mRNA has been observed with compound 5c.

Genetic basis of a cognitive complexity metric

Relational complexity (RC) is a metric reflecting capacity limitation in relational processing. It plays a crucial role in higher cognitive processes and is an endophenotype for several disorders. However, the genetic underpinnings of complex relational processing have not been investigated. Using the classical twin model, we estimated the heritability of RC and genetic overlap with intelligence (IQ), reasoning, and working memory in a twin and sibling sample aged 15-29 years (N = 787). Further, in an exploratory search for genetic loci contributing to RC, we examined associated genetic markers and genes in our Discovery sample and selected loci for replication in four independent samples (ALSPAC, LBC1936, NTR, NCNG), followed by meta-analysis (N>6500) at the single marker level. Twin modelling showed RC is highly heritable (67%), has considerable genetic overlap with IQ (59%), and is a major component of genetic covariation between reasoning and working memory (72%). At the molecular level, we found preliminary support for four single-marker loci (one in the gene DGKB), and at a gene-based level for the NPS gene, having influence on cognition. These results indicate that genetic sources influencing relational processing are a key component of the genetic architecture of broader cognitive abilities. Further, they suggest a genetic cascade, whereby genetic factors influencing capacity limitation in relational processing have a flow-on effect to more complex cognitive traits, including reasoning and working memory, and ultimately, IQ.

Educational attainment influences levels of homozygosity through migration and assortative mating

Individuals with a higher education are more likely to migrate, increasing the chance of meeting a spouse with a different ancestral background. In this context, the presence of strong educational assortment can result in greater ancestry differences within more educated spouse pairs, while less educated individuals are more likely to mate with someone with whom they share more ancestry. We examined the association between educational attainment and F_roh (= the proportion of the genome consisting of runs of homozygosity [ROHs]) in ~2,000 subjects of Dutch ancestry. The subjects’ own educational attainment showed a nominally significant negative association with F_roh (p = .045), while the contribution of parental education to offspring F_roh was highly significant (father: p < 10^-5; mother: p = 9×10^-5), with more educated parents having offspring with fewer ROHs. This association was significantly and fully mediated by the physical distance between parental birthplaces (paternal education: p_mediation = 2.4 × 10^-4; maternal education: p_mediation = 2.3 × 10^-4), which itself was also significantly associated with F_roh (p = 9 × 10^-5). Ancestry-informative principal components from the offspring showed a significantly decreasing association with geography as parental education increased, consistent with the significantly higher migration rates among more educated parents. Parental education also showed a high spouse correlation (Spearman’s ρ = .66, p = 3 × 10^-262). We show that less educated parents are less likely to mate with the more mobile parents with a higher education, creating systematic differences in homozygosity due to ancestry differences not directly captured by ancestry-informative principal components (PCs). Understanding how behaviors influence the genomic structure of a population is highly valuable for studies on the genetic etiology of behavioral, cognitive, and social traits.

Intelligence: shared genetic basis between Mendelian disorders and a polygenic trait

Multiple inquiries into the genetic etiology of human traits indicated an overlap between genes underlying monogenic disorders (eg, skeletal growth defects) and those affecting continuous variability of related quantitative traits (eg, height). Extending the idea of a shared genetic basis between a Mendelian disorder and a classic polygenic trait, we performed an association study to examine the effect of 43 genes implicated in autosomal recessive cognitive disorders on intelligence in an unselected Dutch population (N=1316). Using both single-nucleotide polymorphism (SNP)- and gene-based association testing, we detected an association between intelligence and the genes of interest, with genes ELP2, TMEM135, PRMT10, and RGS7 showing the strongest associations. This is a demonstration of the relevance of genes implicated in monogenic disorders of intelligence to normal-range intelligence, and a corroboration of the utility of employing knowledge on monogenic disorders in identifying the genetic variability underlying complex traits.

Common genetic variants influence human subcortical brain structures

The highly complex structure of the human brain is strongly shaped by genetic influences. Subcortical brain regions form circuits with cortical areas to coordinate movement, learning, memory and motivation, and altered circuits can lead to abnormal behaviour and disease. To investigate how common genetic variants affect the structure of these brain regions, here we conduct genome-wide association studies of the volumes of seven subcortical regions and the intracranial volume derived from magnetic resonance images of 30,717 individuals from 50 cohorts. We identify five novel genetic variants influencing the volumes of the putamen and caudate nucleus. We also find stronger evidence for three loci with previously established influences on hippocampal volume and intracranial volume. These variants show specific volumetric effects on brain structures rather than global effects across structures. The strongest effects were found for the putamen, where a novel intergenic locus with replicable influence on volume (rs945270; P = 1.08 × 10(-33); 0.52% variance explained) showed evidence of altering the expression of the KTN1 gene in both brain and blood tissue. Variants influencing putamen volume clustered near developmental genes that regulate apoptosis, axon guidance and vesicle transport. Identification of these genetic variants provides insight into the causes of variability in human brain development, and may help to determine mechanisms of neuropsychiatric dysfunction.

Polygenic scores associated with educational attainment in adults predict educational achievement and ADHD symptoms in children

The American Psychiatric Association estimates that 3 to 7 per cent of all school aged children are diagnosed with attention deficit hyperactivity disorder (ADHD). Even after correcting for general cognitive ability, numerous studies report a negative association between ADHD and educational achievement. With polygenic scores we examined whether genetic variants that have a positive influence on educational attainment have a protective effect against ADHD. The effect sizes from a large GWA meta-analysis of educational attainment in adults were used to calculate polygenic scores in an independent sample of 12-year-old children from the Netherlands Twin Register. Linear mixed models showed that the polygenic scores significantly predicted educational achievement, school performance, ADHD symptoms and attention problems in children. These results confirm the genetic overlap between ADHD and educational achievement, indicating that one way to gain insight into genetic variants responsible for variation in ADHD is to include data on educational achievement, which are available at a larger scale.