Population Structure and Mating Dynamics in the Social Amoeba Dictyostelium discoideum

Successfully investigating the evolution and maintenance of sex and mating systems can often have as much to do with choosing the right study system as it has to do with asking the right questions. Dictyostelium discoideum has long been the focus of researchers interested in understanding a number of biological processes, such as motility, chemotaxis and development. More recently, attentions have shifted to include questions about the evolution of social and sexual interactions both within and between species. The D. discoideum life cycles, both asexual and sexual, are uniquely social, each requiring a costly sacrificial act. This offers an ideal system for exploring questions about kin recognition, conflict, and the evolution of multicellularity, as well as the evolution of differential sexual investment and mating types. This dissertation focused on understanding the phylogenetic and geographical relationships between clones in D. discoideum and identifying the social and selective pressures that shape its mating system. I introduce this mating system in Chapter 1. In Chapter 2, I investigated genetic variation and population structure in D. discoideum to identify possible factors that could affect interactions between clones. I used DNA sequence data and phylogenetic techniques to show that though D. discoideum clones form a monophyletic group, there is evidence of genetic differentiation among locations (FST = 0.242, P = 0.011), suggesting geographic or other barriers limit gene flow between populations. In chapter 3, I again looked for population structure, this time concentrating on gamete size and sex ratio, to understand selective pressures maintaining multiple mating types in D. discoideum. Evidence suggests that both balancing selection and drift are likely acting on the D. discoideum mating system. I found no differences in gamete size across the three mating types and also no genetic differentiation across three wild populations at the mating type locus. However, I found that mating type frequency varied across these populations, likely due to drift. Chapter 4 focused on understanding the social dynamics of mating in D. discoideum. During macrocyst formation, two cells of complementary mating types fuse to form a zygote. This zygote then consumes hundreds of surrounding amoebae, likely clones of the original two cells, for use as protection and food. I varied the frequencies at which two clones of differing mating types interacted to investigate the possibility that one mating type cheats another by differentially contributing to the cannibalized cells. Contrary to previous claims that mating type I induces mating type II, coercing it to contribute disproportionately more of these cannibalized cells during macrocyst production, I found that these cells are likely contributed relative to their frequency in the population, regardless of mating type. However, I did find evidence for differential contribution to macrocyst production between some pairs of clones, suggesting that cheating can happen between partners during sex, but is rare and clone-specific. Overall, these studies looked for evidence of underlying population structure in D. discoideum that could impact our understanding of social and sexual interactions in this species. I also applied questions about the maintenance of sex usually only asked in two-sex systems to the unique social sexual interactions within D. discoideum in order to expand the understanding of how mating systems evolve and are maintained in nature. I developed and used new tools and techniques for observing the processes important to understanding this unique system and identified genetic and social factors that could impact how individuals interact during both the asexual and sexual life cycles

Sex ratio and gamete size across eastern North America in Dictyostelium discoideum, a social amoeba with three sexes

Theory indicates that numbers of mating types should tend towards infinity or remain at two. The social amoeba, Dictyostelium discoideum, however, has three mating types. It is therefore a mystery how this species has broken the threshold of two mating types, but has not increased towards a much higher number. Frequency-dependent selection on rare types in combination with isogamy, a form of reproduction involving gametes similar in size, could explain the evolution of multiple mating types in this system. Other factors, such as drift, may be preventing the evolution of more than three. We first looked for evidence of isogamy by measuring gamete size associated with each type. We found no evidence of size dissimilarities between gametes. We then looked for evidence of balancing selection, by examining mating type distributions in natural populations and comparing genetic differentiation at the mating type locus to that at more neutral loci. We found that mating type frequency varied among the three populations we examined, with only one of the three showing an even sex ratio, which does not support balancing selection. However, we found more population structure at neutral loci than the mating type locus, suggesting that the three mating types are indeed maintained at intermediate frequencies by balancing selection. Overall, the data are consistent with balancing selection acting on D. discoideum mating types, but with a sufficiently weak rare sex advantage to allow for drift, a potential explanation for why these amoebae have only three mating types

A Common Variant at the 14q32 Endometrial Cancer Risk Locus Activates AKT1 through YY1 Binding.

A recent meta-analysis of multiple genome-wide association and follow-up endometrial cancer case-control datasets identified a novel genetic risk locus for this disease at chromosome 14q32.33. To prioritize the functional SNP(s) and target gene(s) at this locus, we employed an in silico fine-mapping approach using genotyped and imputed SNP data for 6,608 endometrial cancer cases and 37,925 controls of European ancestry. Association and functional analyses provide evidence that the best candidate causal SNP is rs2494737. Multiple experimental analyses show that SNP rs2494737 maps to a silencer element located within AKT1, a member of the PI3K/AKT/MTOR intracellular signaling pathway activated in endometrial tumors. The rs2494737 risk A allele creates a YY1 transcription factor-binding site and abrogates the silencer activity in luciferase assays, an effect mimicked by transfection of YY1 siRNA. Our findings suggest YY1 is a positive regulator of AKT1, mediating the stimulatory effects of rs2494737 increasing endometrial cancer risk. Identification of an endometrial cancer risk allele within a member of the PI3K/AKT signaling pathway, more commonly activated in tumors by somatic alterations, raises the possibility that well tolerated inhibitors targeting this pathway could be candidates for evaluation as chemopreventive agents in individuals at high risk of developing endometrial cancer.The QIMR Berghofer groups were supported by a Rio Tinto Ride to Conquer Cancer (RTCC)/Weekend to End Women's Cancers (WEWC) Grant and NHMRC project grants 1058415 to SLE and 1031333 to ABS. ABS is supported by an NHMRC Senior Research Fellowship (1061779). DFE is a Principal Research Fellow of CR-UK.This is the author accepted manuscript. It is currently under an indefinite embargo pending publication by Cell Press

Fine-mapping of the HNF1B multicancer locus identifies candidate variants that mediate endometrial cancer risk.

Common variants in the hepatocyte nuclear factor 1 homeobox B (HNF1B) gene are associated with the risk of Type II diabetes and multiple cancers. Evidence to date indicates that cancer risk may be mediated via genetic or epigenetic effects on HNF1B gene expression. We previously found single-nucleotide polymorphisms (SNPs) at the HNF1B locus to be associated with endometrial cancer, and now report extensive fine-mapping and in silico and laboratory analyses of this locus. Analysis of 1184 genotyped and imputed SNPs in 6608 Caucasian cases and 37 925 controls, and 895 Asian cases and 1968 controls, revealed the best signal of association for SNP rs11263763 (P = 8.4 × 10(-14), odds ratio = 0.86, 95% confidence interval = 0.82-0.89), located within HNF1B intron 1. Haplotype analysis and conditional analyses provide no evidence of further independent endometrial cancer risk variants at this locus. SNP rs11263763 genotype was associated with HNF1B mRNA expression but not with HNF1B methylation in endometrial tumor samples from The Cancer Genome Atlas. Genetic analyses prioritized rs11263763 and four other SNPs in high-to-moderate linkage disequilibrium as the most likely causal SNPs. Three of these SNPs map to the extended HNF1B promoter based on chromatin marks extending from the minimal promoter region. Reporter assays demonstrated that this extended region reduces activity in combination with the minimal HNF1B promoter, and that the minor alleles of rs11263763 or rs8064454 are associated with decreased HNF1B promoter activity. Our findings provide evidence for a single signal associated with endometrial cancer risk at the HNF1B locus, and that risk is likely mediated via altered HNF1B gene expression

Non-coding RNAs underlie genetic predisposition to breast cancer

Funder: National Breast Cancer Foundation; doi: http://dx.doi.org/10.13039/501100001026Abstract: Background: Genetic variants identified through genome-wide association studies (GWAS) are predominantly non-coding and typically attributed to altered regulatory elements such as enhancers and promoters. However, the contribution of non-coding RNAs to complex traits is not clear. Results: Using targeted RNA sequencing, we systematically annotated multi-exonic non-coding RNA (mencRNA) genes transcribed from 1.5-Mb intervals surrounding 139 breast cancer GWAS signals and assessed their contribution to breast cancer risk. We identify more than 4000 mencRNA genes and show their expression distinguishes normal breast tissue from tumors and different breast cancer subtypes. Importantly, breast cancer risk variants, identified through genetic fine-mapping, are significantly enriched in mencRNA exons, but not the promoters or introns. eQTL analyses identify mencRNAs whose expression is associated with risk variants. Furthermore, chromatin interaction data identify hundreds of mencRNA promoters that loop to regions that contain breast cancer risk variants. Conclusions: We have compiled the largest catalog of breast cancer-associated mencRNAs to date and provide evidence that modulation of mencRNAs by GWAS variants may provide an alternative mechanism underlying complex traits

Five endometrial cancer risk loci identified through genome-wide association analysis.

We conducted a meta-analysis of three endometrial cancer genome-wide association studies (GWAS) and two follow-up phases totaling 7,737 endometrial cancer cases and 37,144 controls of European ancestry. Genome-wide imputation and meta-analysis identified five new risk loci of genome-wide significance at likely regulatory regions on chromosomes 13q22.1 (rs11841589, near KLF5), 6q22.31 (rs13328298, in LOC643623 and near HEY2 and NCOA7), 8q24.21 (rs4733613, telomeric to MYC), 15q15.1 (rs937213, in EIF2AK4, near BMF) and 14q32.33 (rs2498796, in AKT1, near SIVA1). We also found a second independent 8q24.21 signal (rs17232730). Functional studies of the 13q22.1 locus showed that rs9600103 (pairwise r(2) = 0.98 with rs11841589) is located in a region of active chromatin that interacts with the KLF5 promoter region. The rs9600103[T] allele that is protective in endometrial cancer suppressed gene expression in vitro, suggesting that regulation of the expression of KLF5, a gene linked to uterine development, is implicated in tumorigenesis. These findings provide enhanced insight into the genetic and biological basis of endometrial cancer.I.T. is supported by Cancer Research UK and the Oxford Comprehensive Biomedical Research Centre. T.H.T.C. is supported by the Rhodes Trust and the Nuffield Department of Medicine. Funding for iCOGS infrastructure came from the European Community's Seventh Framework Programme under grant agreement 223175 (HEALTH-F2-2009-223175) (COGS), Cancer Research UK (C1287/A10118, C1287/A10710, C12292/A11174, C1281/A12014, C5047/A8384, C5047/A15007, C5047/A10692 and C8197/A16565), the US National Institutes of Health (R01 CA128978, U19 CA148537, U19 CA148065 and U19 CA148112), the US Department of Defense (W81XWH-10-1-0341), the Canadian Institutes of Health Research (CIHR) for the CIHR Team in Familial Risks of Breast Cancer, the Susan G. Komen Foundation for the Cure, the Breast Cancer Research Foundation and the Ovarian Cancer Research Fund. SEARCH recruitment was funded by a programme grant from Cancer Research UK (C490/A10124). Stage 1 and stage 2 case genotyping was supported by the NHMRC (552402 and 1031333). Control data were generated by the WTCCC, and a full list of the investigators who contributed to the generation of the data is available from the WTCCC website. We acknowledge use of DNA from the British 1958 Birth Cohort collection, funded by UK Medical Research Council grant G0000934 and Wellcome Trust grant 068545/Z/02; funding for this project was provided by the Wellcome Trust under award 085475. NSECG was supported by the European Union's Framework Programme 7 CHIBCHA grant and Wellcome Trust Centre for Human Genetics Core Grant 090532/Z/09Z, and CORGI was funded by Cancer Research UK. BCAC is funded by Cancer Research UK (C1287/A10118 and C1287/A12014). OCAC is supported by a grant from the Ovarian Cancer Research Fund thanks to donations by the family and friends of Kathryn Sladek Smith (PPD/RPCI.07) and the UK National Institute for Health Research Biomedical Research Centres at the University of Cambridge.This is the author accepted manuscript. The final version is available from Nature Publishing Group via http://dx.doi.org/10.1038/ng.356

Cerebral microbleeds and intracranial haemorrhage risk in patients anticoagulated for atrial fibrillation after acute ischaemic stroke or transient ischaemic attack (CROMIS-2):a multicentre observational cohort study

Background: Cerebral microbleeds are a potential neuroimaging biomarker of cerebral small vessel diseases that are prone to intracranial bleeding. We aimed to determine whether presence of cerebral microbleeds can identify patients at high risk of symptomatic intracranial haemorrhage when anticoagulated for atrial fibrillation after recent ischaemic stroke or transient ischaemic attack. Methods: Our observational, multicentre, prospective inception cohort study recruited adults aged 18 years or older from 79 hospitals in the UK and one in the Netherlands with atrial fibrillation and recent acute ischaemic stroke or transient ischaemic attack, treated with a vitamin K antagonist or direct oral anticoagulant, and followed up for 24 months using general practitioner and patient postal questionnaires, telephone interviews, hospital visits, and National Health Service digital data on hospital admissions or death. We excluded patients if they could not undergo MRI, had a definite contraindication to anticoagulation, or had previously received therapeutic anticoagulation. The primary outcome was symptomatic intracranial haemorrhage occurring at any time before the final follow-up at 24 months. The log-rank test was used to compare rates of intracranial haemorrhage between those with and without cerebral microbleeds. We developed two prediction models using Cox regression: first, including all predictors associated with intracranial haemorrhage at the 20% level in univariable analysis; and second, including cerebral microbleed presence and HAS-BLED score. We then compared these with the HAS-BLED score alone. This study is registered with ClinicalTrials.gov, number NCT02513316. Findings: Between Aug 4, 2011, and July 31, 2015, we recruited 1490 participants of whom follow-up data were available for 1447 (97%), over a mean period of 850 days (SD 373; 3366 patient-years). The symptomatic intracranial haemorrhage rate in patients with cerebral microbleeds was 9·8 per 1000 patient-years (95% CI 4·0–20·3) compared with 2·6 per 1000 patient-years (95% CI 1·1–5·4) in those without cerebral microbleeds (adjusted hazard ratio 3·67, 95% CI 1·27–10·60). Compared with the HAS-BLED score alone (C-index 0·41, 95% CI 0·29–0·53), models including cerebral microbleeds and HAS-BLED (0·66, 0·53–0·80) and cerebral microbleeds, diabetes, anticoagulant type, and HAS-BLED (0·74, 0·60–0·88) predicted symptomatic intracranial haemorrhage significantly better (difference in C-index 0·25, 95% CI 0·07–0·43, p=0·0065; and 0·33, 0·14–0·51, p=0·00059, respectively). Interpretation: In patients with atrial fibrillation anticoagulated after recent ischaemic stroke or transient ischaemic attack, cerebral microbleed presence is independently associated with symptomatic intracranial haemorrhage risk and could be used to inform anticoagulation decisions. Large-scale collaborative observational cohort analyses are needed to refine and validate intracranial haemorrhage risk scores incorporating cerebral microbleeds to identify patients at risk of net harm from oral anticoagulation. Funding: The Stroke Association and the British Heart Foundation

Genomic analyses identify hundreds of variants associated with age at menarche and support a role for puberty timing in cancer risk

The timing of puberty is a highly polygenic childhood trait that is epidemiologically associated with various adult diseases. Using 1000 Genomes Project-imputed genotype data in up to similar to 370,000 women, we identify 389 independent signals (P <5 x 10(-8)) for age at menarche, a milestone in female pubertal development. In Icelandic data, these signals explain similar to 7.4% of the population variance in age at menarche, corresponding to similar to 25% of the estimated heritability. We implicate similar to 250 genes via coding variation or associated expression, demonstrating significant enrichment in neural tissues. Rare variants near the imprinted genes MKRN3 and DLK1 were identified, exhibiting large effects when paternally inherited. Mendelian randomization analyses suggest causal inverse associations, independent of body mass index (BMI), between puberty timing and risks for breast and endometrial cancers in women and prostate cancer in men. In aggregate, our findings highlight the complexity of the genetic regulation of puberty timing and support causal links with cancer susceptibility

Multiple novel prostate cancer susceptibility signals identified by fine-mapping of known risk loci among Europeans

Genome-wide association studies (GWAS) have identified numerous common prostate cancer (PrCa) susceptibility loci. We have fine-mapped 64 GWAS regions known at the conclusion of the iCOGS study using large-scale genotyping and imputation in 25 723 PrCa cases and 26 274 controls of European ancestry. We detected evidence for multiple independent signals at 16 regions, 12 of which contained additional newly identified significant associations. A single signal comprising a spectrum of correlated variation was observed at 39 regions; 35 of which are now described by a novel more significantly associated lead SNP, while the originally reported variant remained as the lead SNP only in 4 regions. We also confirmed two association signals in Europeans that had been previously reported only in East-Asian GWAS. Based on statistical evidence and linkage disequilibrium (LD) structure, we have curated and narrowed down the list of the most likely candidate causal variants for each region. Functional annotation using data from ENCODE filtered for PrCa cell lines and eQTL analysis demonstrated significant enrichment for overlap with bio-features within this set. By incorporating the novel risk variants identified here alongside the refined data for existing association signals, we estimate that these loci now explain ∼38.9% of the familial relative risk of PrCa, an 8.9% improvement over the previously reported GWAS tag SNPs. This suggests that a significant fraction of the heritability of PrCa may have been hidden during the discovery phase of GWAS, in particular due to the presence of multiple independent signals within the same regio