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.

In Vivo Fast Nonlinear Microscopy Reveals Impairment of Fast Axonal Transport Induced by Molecular Motor Imbalances in the Brain of Zebrafish Larvae.

Cargo transport by molecular motors along microtubules is essential for the function of eukaryotic cells, in particular neurons in which axonal transport defects constitute the early pathological features of neurodegenerative diseases. Mainly studied in motor and sensory neurons, axonal transport is still difficult to characterize in neurons of the brain in absence of appropriate in vivo tools. Here, we measured fast axonal transport by tracing the second harmonic generation (SHG) signal of potassium titanyl phosphate (KTP) nanocrystals (nanoKTP) endocytosed by brain neurons of zebrafish (Zf) larvae. Thanks to the optical translucency of Zf larvae and to the perfect photostability of nanoKTP SHG, we achieved a high scanning speed of 20 frames (of ≈90 µm × 60 µm size) per second in Zf brain. We focused our study on endolysosomal vesicle transport in axons of known polarization, separately analyzing kinesin and dynein motor-driven displacements. To validate our assay, we used either loss-of-function mutations of dynein or kinesin 1 or the dynein inhibitor dynapyrazole and quantified several transport parameters. We successfully demonstrated that dynapyrazole reduces the nanoKTP mobile fraction and retrograde run length consistently, while the retrograde run length increased in kinesin 1 mutants. Taking advantage of nanoKTP SHG directional emission, we also quantified fluctuations of vesicle orientation. Thus, by combining endocytosis of nanocrystals having a nonlinear response, fast two-photon microscopy, and high-throughput analysis, we are able to finely monitor fast axonal transport in vivo in the brain of a vertebrate and reveal subtle axonal transport alterations. The high spatiotemporal resolution achieved in our model may be relevant to precisely investigate axonal transport impairment associated with disease models.

Inhibition of PP1 phosphatase activity by HBx: a mechanism for the activation of hepatitis B virus transcription.

The regulatory protein HBx is essential for hepatitis B virus (HBV) replication in vivo and for transcription of the episomal HBV genome. We previously reported that in infected cells HBx activates genes targeted by the transcription factor CREB [cyclic adenosine monophosphate (cAMP) response element-binding protein]. cAMP induces phosphorylation and activation of CREB, and CREB inactivation is promoted by protein phosphatase 1 (PP1), which binds to CREB through histone deacetylase 1 (HDAC1). We showed that CREB was recruited to HBV DNA. Phosphorylation induced by cAMP had a longer half-life when CREB was bound to the episomal HBV genome compared to when it was bound to the promoter of a host target gene not regulated by HBx, suggesting that the virus has developed a mechanism to favor its own transcription. This mechanism required HBx, which interacted with and inhibited PP1 to extend the half-life of CREB phosphorylation. Silencing of PP1 rescued replication of an HBx-deficient HBV genome, suggesting that HBx enhances viral transcription in part by neutralizing PP1 activity. Our results illustrate a previously unknown mechanism of HBV transcriptional activation by HBx in which HBx interferes with the inactivation of CREB by the PP1 and HDAC1 complex.

Stem cell-like micro-RNA signature driven by Myc in aggressive liver cancer.

Myc activation has been implicated in the pathogenesis of hepatoblastoma (HB), a rare embryonal neoplasm derived from liver progenitor cells. Here, microRNA (miR) expression profiling of 65 HBs evidenced differential patterns related to developmental stage and Myc activity. Undifferentiated aggressive HBs overexpressed the miR-371-3 cluster with concomitant down-regulation of the miR-100/let-7a-2/miR-125b-1 cluster, evoking an ES cell expression profile. ChIP and Myc inhibition assays in hepatoma cells demonstrated that both miR clusters are regulated by Myc in an opposite manner. We show that the two miR clusters exert antagonistic effects on cell proliferation and tumorigenicity. Moreover, their combined deregulation cooperated in modulating the hepatic tumor phenotype, implicating stem cell-like regulation of Myc-dependent miRs in poorly differentiated HBs. Importantly, a four-miR signature representative of these clusters efficiently stratified HB patients, and when applied to 241 hepatocellular carcinomas (HCCs), it identified invasive tumors with a poor prognosis. Our data argue that Myc-driven reprogramming of miR expression patterns contributes to the aggressive phenotype of liver tumors originating from hepatic progenitor cells.