Genome-wide analyses of individual differences in quantitatively assessed reading- and language-related skills in up to 34,000 people.

The use of spoken and written language is a fundamental human capacity. Individual differences in reading- and language-related skills are influenced by genetic variation, with twin-based heritability estimates of 30 to 80% depending on the trait. The genetic architecture is complex, heterogeneous, and multifactorial, but investigations of contributions of single-nucleotide polymorphisms (SNPs) were thus far underpowered. We present a multicohort genome-wide association study (GWAS) of five traits assessed individually using psychometric measures (word reading, nonword reading, spelling, phoneme awareness, and nonword repetition) in samples of 13,633 to 33,959 participants aged 5 to 26 y. We identified genome-wide significant association with word reading (rs11208009, P = 1.098 × 10-8) at a locus that has not been associated with intelligence or educational attainment. All five reading-/language-related traits showed robust SNP heritability, accounting for 13 to 26% of trait variability. Genomic structural equation modeling revealed a shared genetic factor explaining most of the variation in word/nonword reading, spelling, and phoneme awareness, which only partially overlapped with genetic variation contributing to nonword repetition, intelligence, and educational attainment. A multivariate GWAS of word/nonword reading, spelling, and phoneme awareness maximized power for follow-up investigation. Genetic correlation analysis with neuroimaging traits identified an association with the surface area of the banks of the left superior temporal sulcus, a brain region linked to the processing of spoken and written language. Heritability was enriched for genomic elements regulating gene expression in the fetal brain and in chromosomal regions that are depleted of Neanderthal variants. Together, these results provide avenues for deciphering the biological underpinnings of uniquely human traits.

Role of DOCK8 in Cytokine Storm Syndromes.

BACKGROUND: Cytokine storm syndromes (CSS), including hemophagocytic lymphohistiocytosis (HLH), are increasingly recognized as hyperinflammatory states leading to multi-organ failure and death. Familial HLH (FHL) in infancy results from homozygous genetic defects in perforin-mediated cytolysis by CD8 T-lymphocytes and natural killer (NK) cells. Later onset CSS are frequently associated with heterozygous defects in FHL genes, but genetic etiologies for most are unknown. We identified rare DOCK8 variants in CSS patients. OBJECTIVE: We explore the role of CSS patient-derived DOCK8 mutations on cytolytic activity in NK cells. We further study effects of DOCK8 deficiency in murine models of CSS. METHODS: DOCK8 cDNA from 2 unrelated CSS patients with different missense mutations were introduced into human NK-92 NK cells by foamy virus transduction. NK cell degranulation (CD107a), cytolytic activity against K562 target cells, and interferon-gamma (IFNγ) production were explored by flow cytometry. A third CSS patient DOCK8 mRNA splice acceptor site variant was explored by exon trapping. Dock8-/- mice were assessed for features of CSS (weight loss, splenomegaly, hepatic inflammation, cytopenias, and IFNγ levels) upon challenge with lymphocytic choriomeningitis virus (LCMV) and excess IL-18. RESULTS: Both patient DOCK8 missense mutations decreased cytolytic function in NK cells in a partial dominant-negative fashion in vitro. The patient DOCK8 splice variant disrupted mRNA splicing in vitro. LCMV infection promoted CSS in Dock8-/- mice and interacted with excess IL-18 limiting T-cell numbers while promoting CD8 T-cell hyperactivation. CONCLUSION: Mutations in DOCK8 may contribute to CSS-like hyperinflammatory states by altering cytolytic function in a threshold model of disease.

Asthmatic patients with high serum amyloid A have proinflammatory HDL: Implications for augmented systemic and airway inflammation.

RATIONALE: Serum amyloid A (SAA) is bound to high-density lipoproteins (HDL) in blood. Although SAA is increased in the blood of patients with asthma, it is not known whether this modifies asthma severity. OBJECTIVE: We sought to define the clinical characteristics of patients with asthma who have high SAA levels and assess whether HDL from SAA-high patients with asthma is proinflammatory. METHODS: SAA levels in serum from subjects with and without asthma were quantified by ELISA. HDLs isolated from subjects with asthma and high SAA levels were used to stimulate human monocytes and were intravenously administered to BALB/c mice. RESULTS: An SAA level greater than or equal to 108.8 µg/mL was defined as the threshold to identify 11% of an asthmatic cohort (n = 146) as being SAA-high. SAA-high patients with asthma were characterized by increased serum C-reactive protein, IL-6, and TNF-α; older age; and an increased prevalence of obesity and severe asthma. HDL isolated from SAA-high patients with asthma (SAA-high HDL) had an increased content of SAA as compared with HDL from SAA-low patients with asthma and induced the secretion of IL-6, IL-1ß, and TNF-α from human monocytes via a formyl peptide receptor 2/ATP/P2X purinoceptor 7 axis. Intravenous administration to mice of SAA-high HDL, but not normal HDL, induced systemic inflammation and amplified allergen-induced neutrophilic airway inflammation and goblet cell metaplasia. CONCLUSIONS: SAA-high patients with asthma are characterized by systemic inflammation, older age, and an increased prevalence of obesity and severe asthma. HDL from SAA-high patients with asthma is proinflammatory and, when intravenously administered to mice, induces systemic inflammation, and amplifies allergen-induced neutrophilic airway inflammation. This suggests that systemic inflammation induced by SAA-high HDL may augment disease severity in asthma.

Genome-wide association study meta-analysis of dizygotic twinning illuminates genetic regulation of female fecundity.

STUDY QUESTION: Which genetic factors regulate female propensity for giving birth to spontaneous dizygotic (DZ) twins? SUMMARY ANSWER: We identified four new loci, GNRH1, FSHR, ZFPM1, and IPO8, in addition to previously identified loci, FSHB and SMAD3. WHAT IS KNOWN ALREADY: The propensity to give birth to DZ twins runs in families. Earlier, we reported that FSHB and SMAD3 as associated with DZ twinning and female fertility measures. STUDY DESIGN, SIZE, DURATION: We conducted a genome-wide association meta-analysis (GWAMA) of mothers of spontaneous dizygotic (DZ) twins (8265 cases, 264 567 controls) and of independent DZ twin offspring (26 252 cases, 417 433 controls). PARTICIPANTS/MATERIALS, SETTING, METHODS: Over 700 000 mothers of DZ twins, twin individuals and singletons from large cohorts in Australia/New Zealand, Europe, and the USA were carefully screened to exclude twins born after use of ARTs. Genetic association analyses by cohort were followed by meta-analysis, phenome wide association studies (PheWAS), in silico and in vivo annotations, and Zebrafish functional validation. MAIN RESULTS AND THE ROLE OF CHANCE: This study enlarges the sample size considerably from previous efforts, finding four genome-wide significant loci, including two novel signals and a further two novel genes that are implicated by gene level enrichment analyses. The novel loci, GNRH1 and FSHR, have well-established roles in female reproduction whereas ZFPM1 and IPO8 have not previously been implicated in female fertility. We found significant genetic correlations with multiple aspects of female reproduction and body size as well as evidence for significant selection against DZ twinning during human evolution. The 26 top single nucleotide polymorphisms (SNPs) from our GWAMA in European-origin participants weakly predicted the crude twinning rates in 47 non-European populations (r = 0.23 between risk score and population prevalence, s.e. 0.11, 1-tail P = 0.058) indicating that genome-wide association studies (GWAS) are needed in African and Asian populations to explore the causes of their respectively high and low DZ twinning rates. In vivo functional tests in zebrafish for IPO8 validated its essential role in female, but not male, fertility. In most regions, risk SNPs linked to known expression quantitative trait loci (eQTLs). Top SNPs were associated with in vivo reproductive hormone levels with the top pathways including hormone ligand binding receptors and the ovulation cycle. LARGE SCALE DATA: The full DZT GWAS summary statistics will made available after publication through the GWAS catalog (https://www.ebi.ac.uk/gwas/). LIMITATIONS, REASONS FOR CAUTION: Our study only included European ancestry cohorts. Inclusion of data from Africa (with the highest twining rate) and Asia (with the lowest rate) would illuminate further the biology of twinning and female fertility. WIDER IMPLICATIONS OF THE FINDINGS: About one in 40 babies born in the world is a twin and there is much speculation on why twinning runs in families. We hope our results will inform investigations of ovarian response in new and existing ARTs and the causes of female infertility. STUDY FUNDING/COMPETING INTEREST(S): Support for the Netherlands Twin Register came from the Netherlands Organization for Scientific Research (NWO) and The Netherlands Organization for Health Research and Development (ZonMW) grants, 904-61-193, 480-04-004, 400-05-717, Addiction-31160008, 911-09-032, Biobanking and Biomolecular Resources Research Infrastructure (BBMRI.NL, 184.021.007), Royal Netherlands Academy of Science Professor Award (PAH/6635) to DIB, European Research Council (ERC-230374), Rutgers University Cell and DNA Repository (NIMH U24 MH068457-06), the Avera Institute, Sioux Falls, South Dakota (USA) and the National Institutes of Health (NIH R01 HD042157-01A1) and the Genetic Association Information Network (GAIN) of the Foundation for the National Institutes of Health and Grand Opportunity grants 1RC2 MH089951. The QIMR Berghofer Medical Research Institute (QIMR) study was supported by grants from the National Health and Medical Research Council (NHMRC) of Australia (241944, 339462, 389927, 389875, 389891, 389892, 389938, 443036, 442915, 442981, 496610, 496739, 552485, 552498, 1050208, 1075175). L.Y. is funded by Australian Research Council (Grant number DE200100425). The Minnesota Center for Twin and Family Research (MCTFR) was supported in part by USPHS Grants from the National Institute on Alcohol Abuse and Alcoholism (AA09367 and AA11886) and the National Institute on Drug Abuse (DA05147, DA13240, and DA024417). The Women's Genome Health Study (WGHS) was funded by the National Heart, Lung, and Blood Institute (HL043851 and HL080467) and the National Cancer Institute (CA047988 and UM1CA182913), with support for genotyping provided by Amgen. Data collection in the Finnish Twin Registry has been supported by the Wellcome Trust Sanger Institute, the Broad Institute, ENGAGE-European Network for Genetic and Genomic Epidemiology, FP7-HEALTH-F4-2007, grant agreement number 201413, National Institute of Alcohol Abuse and Alcoholism (grants AA-12502, AA-00145, AA-09203, AA15416, and K02AA018755) and the Academy of Finland (grants 100499, 205585, 118555, 141054, 264146, 308248, 312073 and 336823 to J. Kaprio). TwinsUK is funded by the Wellcome Trust, Medical Research Council, Versus Arthritis, European Union Horizon 2020, Chronic Disease Research Foundation (CDRF), Zoe Ltd and the National Institute for Health Research (NIHR) Clinical Research Network (CRN) and Biomedical Research Centre based at Guy's and St Thomas' NHS Foundation Trust in partnership with King's College London. For NESDA, funding was obtained from the Netherlands Organization for Scientific Research (Geestkracht program grant 10000-1002), the Center for Medical Systems Biology (CSMB, NVVO Genomics), Biobanking and Biomolecular Resources Research Infrastructure (BBMRI-NL), VU University's Institutes for Health and Care Research (EMGO+) and Neuroscience Campus Amsterdam, University Medical Center Groningen, Leiden University Medical Center, National Institutes of Health (NIH, ROI D0042157-01A, MH081802, Grand Opportunity grants 1 RC2 Ml-1089951 and IRC2 MH089995). Part of the genotyping and analyses were funded by the Genetic Association Information Network (GAIN) of the Foundation for the National Institutes of Health. Computing was supported by BiG Grid, the Dutch e-Science Grid, which is financially supported by NWO. Work in the Del Bene lab was supported by the Programme Investissements d'Avenir IHU FOReSIGHT (ANR-18-IAHU-01). C.R. was supported by an EU Horizon 2020 Marie Sklodowska-Curie Action fellowship (H2020-MSCA-IF-2014 #661527). H.S. and K.S. are employees of deCODE Genetics/Amgen. The other authors declare no competing financial interests. TRIAL REGISTRATION NUMBER: N/A.

Anorexia nervosa polygenic risk, beyond diagnoses: relationship with adolescent disordered eating and behaviors in an Australian female twin population.

BACKGROUND: It is well established that there is a substantial genetic component to eating disorders (EDs). Polygenic risk scores (PRSs) can be used to quantify cumulative genetic risk for a trait at an individual level. Recent studies suggest PRSs for anorexia nervosa (AN) may also predict risk for other disordered eating behaviors, but no study has examined if PRS for AN can predict disordered eating as a global continuous measure. This study aimed to investigate whether PRS for AN predicted overall levels of disordered eating, or specific lifetime disordered eating behaviors, in an Australian adolescent female population. METHODS: PRSs were calculated based on summary statistics from the largest Psychiatric Genomics Consortium AN genome-wide association study to date. Analyses were performed using genome-wide complex trait analysis to test the associations between AN PRS and disordered eating global scores, avoidance of eating, objective bulimic episodes, self-induced vomiting, and driven exercise in a sample of Australian adolescent female twins recruited from the Australian Twin Registry (N = 383). RESULTS: After applying the false-discovery rate correction, the AN PRS was significantly associated with all disordered eating outcomes. CONCLUSIONS: Findings suggest shared genetic etiology across disordered eating presentations and provide insight into the utility of AN PRS for predicting disordered eating behaviors in the general population. In the future, PRSs for EDs may have clinical utility in early disordered eating risk identification, prevention, and intervention.

Transmission, seroprevalence, and maternal-fetal impact of lymphocytic choriomeningitis virus.

Congenital infections can have devastating short- and long-term impacts on the developing fetus. Lymphocytic choriomeningitis virus (LCMV) is a zoonotic pathogen of concern that causes a severe congenital syndrome but is under-recognized and under-studied. Herein we review data on the natural animal reservoirs of LCMV, modes of transmission to humans, seroprevalence of LCMV worldwide in both pregnant and non-pregnant individuals, mechanisms of viral dissemination to placenta and fetus, and impact of climate change on viral transmission. We highlight opportunities to enhance awareness of congenital LCMV and provide recommendations for prevention and monitoring among at-risk pregnant people. IMPACT: Key message of the article: LCMV is a zoonotic virus that poses a major threat to maternal-fetal health. Adds to the existing literature: We comprehensively address transmission of LCMV from the natural reservoir to the pregnant individual, placenta, and fetus. Impact: Available data call for enhanced patient and provider awareness about congenital LCMV during pregnancy, as well as a need for efforts to better define the seroprevalence and impact of congenital LCMV worldwide.

The RAG1 Ubiquitin Ligase Domain Stimulates Recombination of TCRß and TCRα Genes and Influences Development of αß T Cell Lineages.

RAG1/RAG2 (RAG) endonuclease-mediated assembly of diverse lymphocyte Ag receptor genes by V(D)J recombination is critical for the development and immune function of T and B cells. The RAG1 protein contains a ubiquitin ligase domain that stabilizes RAG1 and stimulates RAG endonuclease activity in vitro. We report in this study that mice with a mutation that inactivates the Rag1 ubiquitin ligase in vitro exhibit decreased rearrangements and altered repertoires of TCRß and TCRα genes in thymocytes and impaired thymocyte developmental transitions that require the assembly and selection of functional TCRß and/or TCRα genes. These Rag1 mutant mice present diminished positive selection and superantigen-mediated negative selection of conventional αß T cells, decreased genesis of invariant NK T lineage αß T cells, and mature CD4+ αß T cells with elevated autoimmune potential. Our findings reveal that the Rag1 ubiquitin ligase domain functions in vivo to stimulate TCRß and TCRα gene recombination and influence differentiation of αß T lineage cells, thereby establishing replete diversity of αß TCRs and populations of αß T cells while restraining generation of potentially autoreactive conventional αß T cells.

Common genetic variants contribute to risk of rare severe neurodevelopmental disorders.

There are thousands of rare human disorders that are caused by single deleterious, protein-coding genetic variants1. However, patients with the same genetic defect can have different clinical presentations2-4, and some individuals who carry known disease-causing variants can appear unaffected5. Here, to understand what explains these differences, we study a cohort of 6,987 children assessed by clinical geneticists to have severe neurodevelopmental disorders such as global developmental delay and autism, often in combination with abnormalities of other organ systems. Although the genetic causes of these neurodevelopmental disorders are expected to be almost entirely monogenic, we show that 7.7% of variance in risk is attributable to inherited common genetic variation. We replicated this genome-wide common variant burden by showing, in an independent sample of 728 trios (comprising a child plus both parents) from the same cohort, that this burden is over-transmitted from parents to children with neurodevelopmental disorders. Our common-variant signal is significantly positively correlated with genetic predisposition to lower educational attainment, decreased intelligence and risk of schizophrenia. We found that common-variant risk was not significantly different between individuals with and without a known protein-coding diagnostic variant, which suggests that common-variant risk affects patients both with and without a monogenic diagnosis. In addition, previously published common-variant scores for autism, height, birth weight and intracranial volume were all correlated with these traits within our cohort, which suggests that phenotypic expression in individuals with monogenic disorders is affected by the same variants as in the general population. Our results demonstrate that common genetic variation affects both overall risk and clinical presentation in neurodevelopmental disorders that are typically considered to be monogenic.

Genetic meta-analysis of twin birth weight shows high genetic correlation with singleton birth weight.

Birth weight (BW) is an important predictor of newborn survival and health and has associations with many adult health outcomes, including cardiometabolic disorders, autoimmune diseases and mental health. On average, twins have a lower BW than singletons as a result of a different pattern of fetal growth and shorter gestational duration. Therefore, investigations into the genetics of BW often exclude data from twins, leading to a reduction in sample size and remaining ambiguities concerning the genetic contribution to BW in twins. In this study, we carried out a genome-wide association meta-analysis of BW in 42 212 twin individuals and found a positive correlation of beta values (Pearson's r = 0.66, 95% confidence interval [CI]: 0.47-0.77) with 150 previously reported genome-wide significant variants for singleton BW. We identified strong positive genetic correlations between BW in twins and numerous anthropometric traits, most notably with BW in singletons (genetic correlation [rg] = 0.92, 95% CI: 0.66-1.18). Genetic correlations of BW in twins with a series of health-related traits closely resembled those previously observed for BW in singletons. Polygenic scores constructed from a genome-wide association study on BW in the UK Biobank demonstrated strong predictive power in a target sample of Dutch twins and singletons. Together, our results indicate that a similar genetic architecture underlies BW in twins and singletons and that future genome-wide studies might benefit from including data from large twin registers.

Genome-wide Association Analysis in Humans Links Nucleotide Metabolism to Leukocyte Telomere Length.

Leukocyte telomere length (LTL) is a heritable biomarker of genomic aging. In this study, we perform a genome-wide meta-analysis of LTL by pooling densely genotyped and imputed association results across large-scale European-descent studies including up to 78,592 individuals. We identify 49 genomic regions at a false dicovery rate (FDR) < 0.05 threshold and prioritize genes at 31, with five highlighting nucleotide metabolism as an important regulator of LTL. We report six genome-wide significant loci in or near SENP7, MOB1B, CARMIL1, PRRC2A, TERF2, and RFWD3, and our results support recently identified PARP1, POT1, ATM, and MPHOSPH6 loci. Phenome-wide analyses in >350,000 UK Biobank participants suggest that genetically shorter telomere length increases the risk of hypothyroidism and decreases the risk of thyroid cancer, lymphoma, and a range of proliferative conditions. Our results replicate previously reported associations with increased risk of coronary artery disease and lower risk for multiple cancer types. Our findings substantially expand current knowledge on genes that regulate LTL and their impact on human health and disease.

GWAS of Dizygotic Twinning in an Enlarged Australian Sample of Mothers of DZ Twins.

Female fertility is a complex trait with age-specific changes in spontaneous dizygotic (DZ) twinning and fertility. To elucidate factors regulating female fertility and infertility, we conducted a genome-wide association study (GWAS) on mothers of spontaneous DZ twins (MoDZT) versus controls (3273 cases, 24,009 controls). This is a follow-up study to the Australia/New Zealand (ANZ) component of that previously reported (Mbarek et al., 2016), with a sample size almost twice that of the entire discovery sample meta-analysed in the previous article (and five times the ANZ contribution to that), resulting from newly available additional genotyping and representing a significant increase in power. We compare analyses with and without male controls and show unequivocally that it is better to include male controls who have been screened for recent family history, than to use only female controls. Results from the SNP based GWAS identified four genomewide significant signals, including one novel region, ZFPM1 (Zinc Finger Protein, FOG Family Member 1), on chromosome 16. Previous signals near FSHB (Follicle Stimulating Hormone beta subunit) and SMAD3 (SMAD Family Member 3) were also replicated (Mbarek et al., 2016). We also ran the GWAS with a dominance model that identified a further locus ADRB2 on chr 5. These results have been contributed to the International Twinning Genetics Consortium for inclusion in the next GWAS meta-analysis (Mbarek et al., in press).

Dennd5b-Deficient Mice are Resistant to PCSK9-Induced Hypercholesterolemia and Diet-Induced Hepatic Steatosis.

Dennd5b plays a pivotal role in intestinal absorption of dietary lipids in mice and is associated with body mass index in humans. This study examined the impact of whole-body Dennd5b deletion on plasma lipid concentrations, atherosclerosis, and hepatic lipid metabolism in mice. Hypercholesterolemia was induced in Dennd5b-/- mice by infection with an adeno-associated virus expressing the proprotein convertase subtilisin/kexin type 9 serine protease (PCSK9) gain-of-function mutation (PCSK9D377Y) and feeding a Western diet for 12 weeks. Body weight and plasma lipid concentrations were monitored over 12 weeks, and then aortic atherosclerosis and hepatic lipid content were quantified. Compared to Dennd5b+/+ mice, Dennd5b-/- mice were resistant to diet-induced weight gain and PCSK9-induced hypercholesterolemia. Atherosclerosis quantified by en face analysis and in aortic root sections, revealed significantly smaller lesions in Dennd5b-/- compared to Dennd5b+/+ mice. Additionally, Dennd5b-/- mice had significantly less hepatic lipid content (triglyceride and cholesterol) compared to Dennd5b+/+ mice. To gain insight into the basis for reduced hepatic lipids, quantitative PCR was used to measure mRNA abundance of genes involved in hepatic lipid metabolism. Key genes involved in hepatic lipid metabolism and lipid storage were differentially expressed in Dennd5b-/- liver including Pparg, Cd36, and Pnpla3. These findings demonstrate a significant impact of Dennd5b on plasma and hepatic lipid concentrations and resistance to PCSK9-induced hypercholesterolemia in the absence of Dennd5b.

Higher polygenic risk for melanoma is associated with improved survival in a high ultraviolet radiation setting.

BACKGROUND: The role of germline genetic factors in determining survival from cutaneous melanoma (CM) is not well understood. OBJECTIVE: To perform a genome-wide association study (GWAS) meta-analysis of melanoma-specific survival (MSS), and test whether a CM-susceptibility polygenic risk score (PRS) is associated with MSS. METHODS: We conducted two Cox proportional-hazard GWAS of MSS using data from the Melanoma Institute Australia, a high ultraviolet (UV) radiation setting (MIA; 5,762 patients with melanoma; 800 melanoma deaths) and UK Biobank (UKB: 5,220 patients with melanoma; 241 melanoma deaths), and combined them in a fixed-effects meta-analysis. Significant (P < 5 × 10-8) results were investigated in the Leeds Melanoma Cohort (LMC; 1,947 patients with melanoma; 370 melanoma deaths). We also developed a CM-susceptibility PRS using a large independent GWAS meta-analysis (23,913 cases, 342,870 controls). The PRS was tested for an association with MSS in the MIA and UKB cohorts. RESULTS: Two loci were significantly associated with MSS in the meta-analysis of MIA and UKB with lead SNPs rs41309643 (G allele frequency 1.6%, HR = 2.09, 95%CI = 1.61-2.71, P = 2.08 × 10-8) on chromosome 1, and rs75682113 (C allele frequency 1.8%, HR = 2.38, 95%CI = 1.77-3.21, P = 1.07 × 10-8) on chromosome 7. While neither SNP replicated in the LMC, rs75682113 was significantly associated in the combined discovery and replication sets. After adjusting for age at diagnosis, sex and the first ten principal components, a one standard deviation increase in the CM-susceptibility PRS was associated with improved MSS in the discovery meta-analysis (HR = 0.88, 95% CI = 0.83-0.94, P = 6.93 × 10-5; I2 = 88%). However, this was only driven by the high UV setting cohort (MIA HR = 0.84, 95% CI = 0.78-0.90). CONCLUSION: We found two loci potentially associated with MSS. Increased genetic susceptibility to develop CM is associated with improved MSS in a high UV setting.

Associations of keratinocyte cancers with snp variants in the sonic hedgehog pathway.

BACKGROUND: Sonic Hedgehog (SHH) pathway dysregulation is implicated in basal cell carcinoma (BCC) development. To evaluate the possible wider role of SHH gene variants in skin carcinogenesis, we assessed associations of genes in the SHH pathway with lifetime development of any keratinocyte cancer (KC), and with developing either BCCs or squamous cell carcinomas (SCCs) exclusively, in a 25-year prospective, population-based study of 1,621 Australians. METHODS: We genotyped 795 unrelated adults with available blood samples: 311 cases with any KC (186 developing BCCs-only, 55 SCCs-only, 70 BCCs and SCCs) and 484 controls. We compared allele frequencies of 158 independent SNPs across 43 SHH genes between cases and controls, and performed a gene-based analysis. RESULTS: We found associations between SNP rs4848627 (GLI2) (related to DNA synthesis in keratinocytes) and development of any KC (OR = 1.53; 95% CI = 1.06-2.13, P < 0.01) and SCCs exclusively (OR = 2.12; 95%CI = 1.39-3.23, P < 0.01). SNP rs3217882 located in CCND2 was associated with exclusive BCC development (OR = 1.43, CI = 1.12-1.82, P < 0.01). The gene-based analysis suggested an association of PRKACG (protein kinase cAMP-activated catalytic subunit gamma) with any KC (P = 0.013). CONCLUSION: We conclude that variants located in genes in the SHH pathway may are involved in SCC as well as BCC development.

Genome-wide association meta-analysis of nicotine metabolism and cigarette consumption measures in smokers of European descent.

Smoking behaviors, including amount smoked, smoking cessation, and tobacco-related diseases, are altered by the rate of nicotine clearance. Nicotine clearance can be estimated using the nicotine metabolite ratio (NMR) (ratio of 3'hydroxycotinine/cotinine), but only in current smokers. Advancing the genomics of this highly heritable biomarker of CYP2A6, the main metabolic enzyme for nicotine, will also enable investigation of never and former smokers. We performed the largest genome-wide association study (GWAS) to date of the NMR in European ancestry current smokers (n = 5185), found 1255 genome-wide significant variants, and replicated the chromosome 19 locus. Fine-mapping of chromosome 19 revealed 13 putatively causal variants, with nine of these being highly putatively causal and mapping to CYP2A6, MAP3K10, ADCK4, and CYP2B6. We also identified a putatively causal variant on chromosome 4 mapping to TMPRSS11E and demonstrated an association between TMPRSS11E variation and a UGT2B17 activity phenotype. Together the 14 putatively causal SNPs explained ~38% of NMR variation, a substantial increase from the ~20 to 30% previously explained. Our additional GWASs of nicotine intake biomarkers showed that cotinine and smoking intensity (cotinine/cigarettes per day (CPD)) shared chromosome 19 and chromosome 4 loci with the NMR, and that cotinine and a more accurate biomarker, cotinine + 3'hydroxycotinine, shared a chromosome 15 locus near CHRNA5 with CPD and Pack-Years (i.e., cumulative exposure). Understanding the genetic factors influencing smoking-related traits facilitates epidemiological studies of smoking and disease, as well as assists in optimizing smoking cessation support, which in turn will reduce the enormous personal and societal costs associated with smoking.

IFNs Drive Development of Novel IL-15-Responsive Macrophages.

Disruption in homeostasis of IL-15 is linked to poor maternal and fetal outcomes during pregnancy. The only cells described to respond to IL-15 at the early maternal-fetal interface have been NK cells. We now show a novel population of macrophages, evident in several organs but enriched in the uterus of mice and humans, expressing the ß-chain of the IL-15R complex (CD122) and responding to IL-15. CD122+ macrophages (CD122+Macs) are morphologic, phenotypic, and transcriptomic macrophages that can derive from bone marrow monocytes. CD122+Macs develop in the uterus and placenta with kinetics that mirror IFN activity at the maternal-fetal interface. M-CSF permits macrophages to express CD122, and IFNs are sufficient to drive expression of CD122 on macrophages. Neither type I nor type II IFNs are required to generate CD122+Macs, however. In response to IL-15, CD122+Macs activate the ERK signaling cascade and enhance production of proinflammatory cytokines after stimulation with the TLR9 agonist CpG. Finally, we provide evidence of human cells that phenocopy murine CD122+Macs in secretory phase endometrium during the implantation window and in first-trimester uterine decidua. Our data support a model wherein IFNs local to the maternal-fetal interface direct novel IL-15-responsive macrophages with the potential to mediate IL-15 signals critical for optimal outcomes of pregnancy.

A New Structural Model of Apolipoprotein B100 Based on Computational Modeling and Cross Linking.

ApoB-100 is a member of a large lipid transfer protein superfamily and is one of the main apolipoproteins found on low-density lipoprotein (LDL) and very low-density lipoprotein (VLDL) particles. Despite its clinical significance for the development of cardiovascular disease, there is limited information on apoB-100 structure. We have developed a novel method based on the "divide and conquer" algorithm, using PSIPRED software, by dividing apoB-100 into five subunits and 11 domains. Models of each domain were prepared using I-TASSER, DEMO, RoseTTAFold, Phyre2, and MODELLER. Subsequently, we used disuccinimidyl sulfoxide (DSSO), a new mass spectrometry cleavable cross-linker, and the known position of disulfide bonds to experimentally validate each model. We obtained 65 unique DSSO cross-links, of which 87.5% were within a 26 Å threshold in the final model. We also evaluated the positions of cysteine residues involved in the eight known disulfide bonds in apoB-100, and each pair was measured within the expected 5.6 Å constraint. Finally, multiple domains were combined by applying constraints based on detected long-range DSSO cross-links to generate five subunits, which were subsequently merged to achieve an uninterrupted architecture for apoB-100 around a lipoprotein particle. Moreover, the dynamics of apoB-100 during particle size transitions was examined by comparing VLDL and LDL computational models and using experimental cross-linking data. In addition, the proposed model of receptor ligand binding of apoB-100 provides new insights into some of its functions.

Combined analysis of keratinocyte cancers identifies novel genome-wide loci.

The keratinocyte cancers (KC), basal cell carcinoma (BCC) and squamous cell carcinoma (SCC) are the most common cancers in fair-skinned people. KC treatment represents the second highest cancer healthcare expenditure in Australia. Increasing our understanding of the genetic architecture of KC may provide new avenues for prevention and treatment. We first conducted a series of genome-wide association studies (GWAS) of KC across three European ancestry datasets from Australia, Europe and USA, and used linkage disequilibrium (LD) Score regression (LDSC) to estimate their pairwise genetic correlations. We employed a multiple-trait approach to map genes across the combined set of KC GWAS (total N = 47 742 cases, 634 413 controls). We also performed meta-analyses of BCC and SCC separately to identify trait specific loci. We found substantial genetic correlations (generally 0.5-1) between BCC and SCC suggesting overlapping genetic risk variants. The multiple trait combined KC GWAS identified 63 independent genome-wide significant loci, 29 of which were novel. Individual separate meta-analyses of BCC and SCC identified an additional 13 novel loci not found in the combined KC analysis. Three new loci were implicated using gene-based tests. New loci included common variants in BRCA2 (distinct to known rare high penetrance cancer risk variants), and in CTLA4, a target of immunotherapy in melanoma. We found shared and trait specific genetic contributions to BCC and SCC. Considering both, we identified a total of 79 independent risk loci, 45 of which are novel.

Ethnic Identity and Genome Wide Runs of Homozygosity.

It is long known that inbreeding increases the detrimental effects of recessive sequence variants in "Runs of Homozygosity" (ROHs). However, although the phenotypic association of ROH has been investigated for a variety of traits, the statistical power of the results often remains limited as a sufficiently high number of cases are available for only a restricted number of traits. In the present study, we aim to analyze the association of runs of homozygosity with the trait "in-group ethnic favoritism". This analysis assumes that if ethnic identity is important for an individual, that individual may tend to marry more frequently within their own group and therefore ROH are expected to increase. We hypothesize that an attitude preferring one's own ethnic group may be associated with a stronger tendency of inbreeding and, as a result, with more and longer ROHs. Accordingly, we investigated the association between the attitude to someone's own ethnicity and ROH, using the Wisconsin Longitudinal data (WLS, total N ~ 9000) as discovery data set and the Brisbane Twin data as replication data set (N ~ 8000). We find that both the number as well as the total length of homozygous segments are significantly positively associated with "in-group ethnic favoritism", independent of the method used for ROH calculation.

The transcription factors T-bet and Eomes control key checkpoints of natural killer cell maturation.

Natural killer (NK) cells play critical roles defending against tumors and pathogens. We show that mice lacking both transcription factors Eomesodermin (Eomes) and T-bet failed to develop NK cells. Developmental stability of immature NK cells constitutively expressing the death ligand TRAIL depended on T-bet. Conversely, maturation characterized by loss of constitutive TRAIL expression and induction of Ly49 receptor diversity and integrin CD49b (DX5(+)) required Eomes. Mature NK cells from which Eomes was deleted reverted to phenotypic immaturity if T-bet was present or downregulated NK lineage antigens if T-bet was absent, despite retaining expression of Ly49 receptors. Fetal and adult hepatic hematopoiesis restricted Eomes expression and limited NK development to the T-bet-dependent, immature stage, whereas medullary hematopoiesis permitted Eomes-dependent NK maturation in adult mice. These findings reveal two sequential, genetically separable checkpoints of NK cell maturation, the progression of which is metered largely by the anatomic localization of hematopoiesis.