CYP19A1 fine-mapping and Mendelian randomization: Estradiol is causal for endometrial cancer

Candidate gene studies have reported CYP19A1 variants to be associated with endometrial cancer and with estradiol (E2) concentrations. We analyzed 2937 single nucleotide polymorphisms (SNPs) in 6608 endometrial cancer cases and 37 925 controls and report the first genome wide-significant association between endometrial cancer and a CYP19A1 SNP (rs727479 in intron 2, P=4.8x10(-11)). SNP rs727479 was also among those most strongly associated with circulating E2 concentrations in 2767 post-menopausal controls (P=7.4x10(-8)). The observed endometrial cancer odds ratio per rs727479 A-allele (1.15, CI=1.11-1.21) is compatible with that predicted by the observed effect on E2 concentrations (1.09, CI=1.03-1.21), consistent with the hypothesis that endometrial cancer risk is driven by E2. From 28 candidate-causal SNPs, 12 co-located with three putative gene-regulatory elements and their risk alleles associated with higher CYP19A1 expression in bioinformatical analyses. For both phenotypes, the associations with rs727479 were stronger among women with a higher BMI (Pinteraction=0.034 and 0.066 respectively), suggesting a biologically plausible gene-environment interaction.

Shared genetic basis for migraine and ischemic stroke: A genome-wide analysis of common variants

OBJECTIVE To quantify genetic overlap between migraine and ischemic stroke (IS) with respect to common genetic variation. METHODS We applied 4 different approaches to large-scale meta-analyses of genome-wide data on migraine (23,285 cases and 95,425 controls) and IS (12,389 cases and 62,004 controls). First, we queried known genome-wide significant loci for both disorders, looking for potential overlap of signals. We then analyzed the overall shared genetic load using polygenic scores and estimated the genetic correlation between disease subtypes using data derived from these models. We further interrogated genomic regions of shared risk using analysis of covariance patterns between the 2 phenotypes using cross-phenotype spatial mapping. RESULTS We found substantial genetic overlap between migraine and IS using all 4 approaches. Migraine without aura (MO) showed much stronger overlap with IS and its subtypes than migraine with aura (MA). The strongest overlap existed between MO and large artery stroke (LAS; p = 6.4 x 10(-28) for the LAS polygenic score in MO) and between MO and cardioembolic stroke (CE; p = 2.7 x 10(-20) for the CE score in MO). CONCLUSIONS Our findings indicate shared genetic susceptibility to migraine and IS, with a particularly strong overlap between MO and both LAS and CE pointing towards shared mechanisms. Our observations on MA are consistent with a limited role of common genetic variants in this subtype.

Fine mapping of variants associated with endometriosis in the WNT4 region on chromosome 1p36

Genome-wide association studies show strong evidence of association with endometriosis for markers on chromosome 1p36 spanning the potential candidate genes WNT4, CDC42 and LINC00339. WNT4 is involved in development of the uterus, and the expression of CDC42 and LINC00339 are altered in women with endometriosis. We conducted fine mapping to examine the role of coding variants in WNT4 and CDC42 and determine the key SNPs with strongest evidence of association in this region. We identified rare coding variants in WNT4 and CDC42 present only in endometriosis cases. The frequencies were low and cannot account for the common signal associated with increased risk of endometriosis. Genotypes for five common SNPs in the region of chromosome 1p36 show stronger association signals when compared with rs7521902 reported in published genome scans. Of these, three SNPs rs12404660, rs3820282, and rs55938609 were located in DNA sequences with potential functional roles including overlap with transcription factor binding sites for FOXA1, FOXA2, ESR1, and ESR2. Functional studies will be required to identify the gene or genes implicated in endometriosis risk.

Genome-wide association and genetic functional studies identify autism susceptibility candidate 2 gene (AUTS2) in the regulation of alcohol consumption

Alcohol consumption is a moderately heritable trait, but the genetic basis in humans is largely unknown, despite its clinical and societal importance. We report a genome-wide association study meta-analysis of approximately 2.5 million directly genotyped or imputed SNPs with alcohol consumption (gram per day per kilogram body weight) among 12 population-based samples of European ancestry, comprising 26,316 individuals, with replication genotyping in an additional 21,185 individuals. SNP rs6943555 in autism susceptibility candidate 2 gene (AUTS2) was associated with alcohol consumption at genome-wide significance (P = 4 x 10(-8) to P = 4 x 10(-9)). We found a genotype-specific expression of AUTS2 in 96 human prefrontal cortex samples (P = 0.026) and significant (P < 0.017) differences in expression of AUTS2 in whole-brain extracts of mice selected for differences in voluntary alcohol consumption. Down-regulation of an AUTS2 homolog caused reduced alcohol sensitivity in Drosophila (P < 0.001). Our finding of a regulator of alcohol consumption adds knowledge to our understanding of genetic mechanisms influencing alcohol drinking behavior.

Association mapping

Association mapping seeks to identify marker alleles present at significantly different frequencies in cases carrying a particular disease or trait compared with controls. Genome-wide association studies are increasingly replacing candidate gene-based association studies for complex diseases, where a number of loci are likely to contribute to disease risk and the effect size of each particular risk allele is typically modest or low. Good study design is essential to the success of an association study, and factors such as the heritability of the disease under investigation, the choice of controls, statistical power, multiple testing and whether the association can be replicated need to be considered before beginning. Likewise, thorough quality control of the genotype data needs to be undertaken prior to running any association analyses. Finally, it should be kept in mind that a significant genetic association is not proof positive that a particular genetic locus causes a disease, but rather an important first step in discovering the genetic variants underlying a complex disease.

Genome-wide linkage analysis of multiple measures of neuroticism of 2 large cohorts from Australia and the Netherlands

CONTEXT People meeting diagnostic criteria for anxiety or depressive disorders tend to score high on the personality scale of neuroticism. Studying this personality dimension can give insights into the etiology of these important psychiatric disorders. OBJECTIVES To undertake a comprehensive genome-wide linkage study of neuroticism using large study samples that have been measured multiple times and to compare the results between countries for replication and across time within countries for consistency. DESIGN Genome-wide linkage scan. SETTING Twin individuals and their family members from Australia and the Netherlands. PARTICIPANTS Nineteen thousand six hundred thirty-five sibling pairs completed self-report questionnaires for neuroticism up to 5 times over a period of up to 22 years. Five thousand sixty-nine sibling pairs were genotyped with microsatellite markers. METHODS Nonparametric linkage analyses were conducted in MERLIN-REGRESS for the mean neuroticism scores averaged across time. Additional analyses were conducted for the time-specific measures of neuroticism from each country to investigate consistency of linkage results. RESULTS Three chromosomal regions exceeded empirically derived thresholds for suggestive linkage using mean neuroticism scores: 10p 5 Kosambi cM (cM) (Dutch study sample), 14q 103 cM (Dutch study sample), and 18q 117 cM (combined Australian and Dutch study sample), but only 14q retained significance after correction for multiple testing. These regions all showed evidence for linkage in individual time-specific measures of neuroticism and 1 (18q) showed some evidence for replication between countries. Linkage intervals for these regions all overlap with regions identified in other studies of neuroticism or related traits and/or in studies of anxiety in mice. CONCLUSIONS Our results demonstrate the value of the availability of multiple measures over time and add to the optimism reported in recent reviews for replication of linkage regions for neuroticism. These regions are likely to harbor causal variants for neuroticism and its related psychiatric disorders and can inform prioritization of results from genome-wide association studies.

Mapping the carriage of flaA-restriction fragment length polymorphism Campylobacter genotypes on poultry carcasses through the processing chain and comparison to clinical isolates

Poultry are considered a major source for campylobacteriosis in humans. A total of 1866 Campylobacter spp. isolates collected through the poultry processing chain were typed using flaA-restriction fragment length polymorphism to measure the impact of processing on the genotypes present. Temporally related human clinical isolates (n = 497) were also typed. Isolates were obtained from whole chicken carcass rinses of chickens collected before scalding, after scalding, before immersion chilling, after immersion chilling and after packaging as well as from individual caecal samples. A total of 32 genotypes comprising at least four isolates each were recognised. Simpson's Index of Diversity (D) was calculated for each sampling site within each flock, for each flock as a whole and for the clinical isolates. From caecal collection to after packaging samples the D value did not change in two flocks, decreased in one flock and increased in the fourth flock. Dominant genotypes occurred in each flock but their constitutive percentages changed through processing. There were 23 overlapping genotypes between clinical and chicken isolates. The diversity of Campylobacter is flock dependant and may alter through processing. This study confirms that poultry are a source of campylobacteriosis in the Australian population although other sources may contribute.

Genome-wide association analysis and differential expression analysis of resistance to Sclerotinia stem rot in Brassica napus

Brassica napus is one of the most important oil crops in the world, and stem rot caused by the fungus Sclerotinia sclerotiorum results in major losses in yield and quality. To elucidate resistance genes and pathogenesis-related genes, genome-wide association analysis of 347 accessions was performed using the Illumina 60K Brassica SNP (single nucleotide polymorphism) array. In addition, the detached stem inoculation assay was used to select five highly resistant (R) and susceptible (S) B. napus lines, 48 h postinoculation with S. sclerotiorum for transcriptome sequencing. We identified 17 significant associations for stem resistance on chromosomes A8 and C6, five of which were on A8 and 12 on C6. The SNPs identified on A8 were located in a 409-kb haplotype block, and those on C6 were consistent with previous QTL mapping efforts. Transcriptome analysis suggested that S. sclerotiorum infection activates the immune system, sulphur metabolism, especially glutathione (GSH) and glucosinolates in both R and S genotypes. Genes found to be specific to the R genotype related to the jasmonic acid pathway, lignin biosynthesis, defence response, signal transduction and encoding transcription factors. Twenty-four genes were identified in both the SNP-trait association and transcriptome sequencing analyses, including a tau class glutathione S-transferase (GSTU) gene cluster. This study provides useful insight into the molecular mechanisms underlying the plant's response to S. sclerotiorum.

Computational Methods for Reconstruction and Analysis of Genome-Scale Metabolic Networks

Metabolism is the cellular subsystem responsible for generation of energy from nutrients and production of building blocks for larger macromolecules. Computational and statistical modeling of metabolism is vital to many disciplines including bioengineering, the study of diseases, drug target identification, and understanding the evolution of metabolism. In this thesis, we propose efficient computational methods for metabolic modeling. The techniques presented are targeted particularly at the analysis of large metabolic models encompassing the whole metabolism of one or several organisms. We concentrate on three major themes of metabolic modeling: metabolic pathway analysis, metabolic reconstruction and the study of evolution of metabolism. In the first part of this thesis, we study metabolic pathway analysis. We propose a novel modeling framework called gapless modeling to study biochemically viable metabolic networks and pathways. In addition, we investigate the utilization of atom-level information on metabolism to improve the quality of pathway analyses. We describe efficient algorithms for discovering both gapless and atom-level metabolic pathways, and conduct experiments with large-scale metabolic networks. The presented gapless approach offers a compromise in terms of complexity and feasibility between the previous graph-theoretic and stoichiometric approaches to metabolic modeling. Gapless pathway analysis shows that microbial metabolic networks are not as robust to random damage as suggested by previous studies. Furthermore the amino acid biosynthesis pathways of the fungal species Trichoderma reesei discovered from atom-level data are shown to closely correspond to those of Saccharomyces cerevisiae. In the second part, we propose computational methods for metabolic reconstruction in the gapless modeling framework. We study the task of reconstructing a metabolic network that does not suffer from connectivity problems. Such problems often limit the usability of reconstructed models, and typically require a significant amount of manual postprocessing. We formulate gapless metabolic reconstruction as an optimization problem and propose an efficient divide-and-conquer strategy to solve it with real-world instances. We also describe computational techniques for solving problems stemming from ambiguities in metabolite naming. These techniques have been implemented in a web-based sofware ReMatch intended for reconstruction of models for 13C metabolic flux analysis. In the third part, we extend our scope from single to multiple metabolic networks and propose an algorithm for inferring gapless metabolic networks of ancestral species from phylogenetic data. Experimenting with 16 fungal species, we show that the method is able to generate results that are easily interpretable and that provide hypotheses about the evolution of metabolism.

Mapping the common antigenic determinants in avidin and streptavidin

An epitope scan analysis of the whole sequence of avidin and core streptavidin using polyclonal antibodies to these two antigens reveal the presence of multiple common epitopes in both the proteins. These antigenic determinants consist mostly of either identical or similar residues. The antibody recognition sites in both antigens are shown to be localized to homologous regions.

Whole-exome re-sequencing in a family quartet identifies POP1 mutations as the cause of a novel skeletal dysplasia

Recent advances in DNA sequencing have enabled mapping of genes for monogenic traits in families with small pedigrees and even in unrelated cases. We report the identification of disease-causing mutations in a rare, severe, skeletal dysplasia, studying a family of two healthy unrelated parents and two affected children using whole-exome sequencing. The two affected daughters have clinical and radiographic features suggestive of anauxetic dysplasia (OMIM 607095), a rare form of dwarfism caused by mutations of RMRP. However, mutations of RMRP were excluded in this family by direct sequencing. Our studies identified two novel compound heterozygous loss-of-function mutations in POP1, which encodes a core component of the RNase mitochondrial RNA processing (RNase MRP) complex that directly interacts with the RMRP RNA domains that are affected in anauxetic dysplasia. We demonstrate that these mutations impair the integrity and activity of this complex and that they impair cell proliferation, providing likely molecular and cellular mechanisms by which POP1 mutations cause this severe skeletal dysplasia. © 2011 Glazov et al.

Genetic Mechanisms Underlying Autism Spectrum Disorders

Autism is a childhood-onset developmental disorder characterized by deficits in reciprocal social interaction, verbal and non-verbal communication, and dependence on routines and rituals. It belongs to a spectrum of disorders (autism spectrum disorders, ASDs) which share core symptoms but show considerable variation in severity. The whole spectrum affects 0.6-0.7% of children worldwide, inducing a substantial public health burden and causing suffering to the affected families. Despite having a very high heritability, ASDs have shown exceptional genetic heterogeneity, which has complicated the identification of risk variants and left the etiology largely unknown. However, recent studies suggest that rare, family-specific factors contribute significantly to the genetic basis of ASDs. In this study, we investigated the role of DISC1 (Disrupted-in-schizophrenia-1) in ASDs, and identified association with markers and haplotypes previously associated with psychiatric phenotypes. We identified four polymorphic micro-RNA target sites in the 3 UTR of DISC1, and showed that hsa-miR-559 regulates DISC1 expression in vitro in an allele-specific manner. We also analyzed an extended autism pedigree with genealogical roots in Central Finland reaching back to the 17th century. To take advantage of the beneficial characteristics of population isolates to gene mapping and reduced genetic heterogeneity observed in distantly related individuals, we performed a microsatellite-based genome-wide screen for linkage and linkage disequilibrium in this pedigree. We identified a putative autism susceptibility locus on chromosome 19p13.3 and obtained further support for previously reported loci at 1q23 and 15q11-q13. To follow-up these findings, we extended our study sample from the same sub-isolate and initiated a genome-wide analysis of homozygosity and allelic sharing using high-density SNP markers. We identified a small number of haplotypes shared by different subsets of the genealogically connected cases, along with convergent biological pathways from SNP and gene expression data, which highlighted axon guidance molecules in the pathogenesis of ASDs. In conclusion, the results obtained in this thesis show that multiple distinct genetic variants are responsible for the ASD phenotype even within single pedigrees from an isolated population. We suggest that targeted resequencing of the shared haplotypes, linkage regions, and other susceptibility loci is essential to identify the causal variants. We also report a possible micro-RNA mediated regulatory mechanism, which might partially explain the wide-range neurobiological effects of the DISC1 gene.

Physical mapping of the genomic DNA of the Oryctes rhinoceros baculovirus, KI

A non-occluded baculovirus, OBV-KI has been isolated from the insect pest, Oryctes rhinoceros. The viral genome is estimated to be 123 kb, with a G + C content of 43 mol% and no detectible methylated bases. A restriction map of the OBV-KI genome for BamHI, EcoRI, HindIII, PstI, SalI and XbaI has been constructed.

Draft genome of a commonly misdiagnosed multidrug resistant pathogen Candida auris

Background: Candida auris is a multidrug resistant, emerging agent of fungemia in humans. Its actual global distribution remains obscure as the current commercial methods of clinical diagnosis misidentify it as C. haemulonii. Here we report the first draft genome of C. auris to explore the genomic basis of virulence and unique differences that could be employed for differential diagnosis. Results: More than 99.5 % of the C. auris genomic reads did not align to the current whole (or draft) genome sequences of Candida albicans, Candida lusitaniae, Candida glabrata and Saccharomyces cerevisiae; thereby indicating its divergence from the active Candida clade. The genome spans around 12.49 Mb with 8527 predicted genes. Functional annotation revealed that among the sequenced Candida species, it is closest to the hemiascomycete species Clavispora lusitaniae. Comparison with the well-studied species Candida albicans showed that it shares significant virulence attributes with other pathogenic Candida species such as oligopeptide transporters, mannosyl transfersases, secreted proteases and genes involved in biofilm formation. We also identified a plethora of transporters belonging to the ABC and major facilitator superfamily along with known MDR transcription factors which explained its high tolerance to antifungal drugs. Conclusions: Our study emphasizes an urgent need for accurate fungal screening methods such as PCR and electrophoretic karyotyping to ensure proper management of fungemia. Our work highlights the potential genetic mechanisms involved in virulence and pathogenicity of an important emerging human pathogen namely C. auris. Owing to its diversity at the genomic scale; we expect the genome sequence to be a useful resource to map species specific differences that will help develop accurate diagnostic markers and better drug targets.

Comparative molecular cytogenetic studies in the order Carnivora: mapping chromosomal rearrangements onto the phylogenetic tree

We have made a set of chromosome-specific painting probes for the American mink by degenerate oligonucleotide primed-PCR (DOP-PCR) amplification of flow-sorted chromosomes. The painting probes were used to delimit homologous chromosomal segments among human, red fox, dog, cat and eight species of the family Mustelidae, including the European mink, steppe and forest polecats, least weasel, mountain weasel, Japanese sable, striped polecat, and badger. Based on the results of chromosome painting and G-banding, comparative maps between these species have been established. The integrated map demonstrates a high level of karyotype conservation among mustelid species. Comparative analysis of the conserved chromosomal segments among mustelids and outgroup species revealed 18 putative ancestral autosomal segments that probably represent the ancestral chromosomes, or chromosome arms, in the karyotype of the most recent ancestor of the family Mustelidae. The proposed 2n = 38 ancestral Mustelidae karyotype appears to have been retained in some modern mustelids, e.g., Martes, Lutra, ktonyx, and Vormela. The derivation of the mustelid karyotypes from the putative ancestral state resulted from centric fusions, fissions, the addition of heterochromatic arms, and occasional pericentric inversions. Our results confirm many of the evolutionary conclusions suggested by other data and strengthen the topology of the carnivore phylogenetic tree through the inclusion of genome-wide chromosome rearrangements. Copyright (C) 2002 S. KargerAG, Basel.