Pan-genome analysis of the emerging foodborne pathogen Cronobacter spp. suggests a species-level bidirectional divergence driven by niche adaptation

Background: Members of the genus Cronobacter are causes of rare but severe illness in neonates and preterm infants following the ingestion of contaminated infant formula. Seven species have been described and two of the species genomes were subsequently published. In this study, we performed comparative genomics on eight strains of Cronobacter, including six that we sequenced (representing six of the seven species) and two previously published, closed genomes.

Results: We identified and characterized the features associated with the core and pan genome of the genus Cronobacter in an attempt to understand the evolution of these bacteria and the genetic content of each species. We identified 84 genomic regions that are present in two or more Cronobacter genomes, along with 45 unique genomic regions. Many potentially horizontally transferred genes, such as lysogenic prophages, were also identified. Most notable among these were several type six secretion system gene clusters, transposons that carried tellurium, copper and/or silver resistance genes, and a novel integrative conjugative element.

Conclusions: Cronobacter have diverged into two clusters, one consisting of C. dublinensis and C. muytjensii (Cdub-Cmuy) and the other comprised of C. sakazakii, C. malonaticus, C. universalis, and C. turicensis, (Csak-Cmal-Cuni-Ctur) from the most recent common ancestral species. While several genetic determinants for plant-association and human virulence could be found in the core genome of Cronobacter, the four Cdub-Cmuy clade genomes contained several accessory genomic regions important for survival in a plant-associated environmental niche, while the Csak-Cmal-Cuni-Ctur clade genomes harbored numerous virulence-related genetic traits.

Genome-wide association analysis identifies 13 new risk loci for schizophrenia

Schizophrenia is an idiopathic mental disorder with a heritable component and a substantial public health impact. We conducted a multi-stage genome-wide association study (GWAS) for schizophrenia beginning with a Swedish national sample (5,001 cases and 6,243 controls) followed by meta-analysis with previous schizophrenia GWAS (8,832 cases and 12,067 controls) and finally by replication of SNPs in 168 genomic regions in independent samples (7,413 cases, 19,762 controls and 581 parent-offspring trios). We identified 22 loci associated at genome-wide significance; 13 of these are new, and 1 was previously implicated in bipolar disorder. Examination of candidate genes at these loci suggests the involvement of neuronal calcium signaling. We estimate that 8,300 independent, mostly common SNPs (95% credible interval of 6,300-10,200 SNPs) contribute to risk for schizophrenia and that these collectively account for at least 32% of the variance in liability. Common genetic variation has an important role in the etiology of schizophrenia, and larger studies will allow more detailed understanding of this disorder.

Genome-wide scans of three independent sets of 90 Irish multiplex schizophrenia families and follow-up of selected regions in all families provides evidence for multiple susceptibility genes

From our linkage study of Irish families with a high density of schizophrenia, we have previously reported evidence for susceptibility genes in regions 5q21-31, 6p24-21, 8p22-21, and 10p15-p11. In this report, we describe the cumulative results from independent genome scans of three a priori random subsets of 90 families each, and from multipoint analysis of all 270 families in ten regions. Of these ten regions, three (13q32, 18p11-q11, and 18q22-23) did not generate scores above the empirical baseline pairwise scan results, and one (6q13-26) generated a weak signal. Six other regions produced more positive pairwise and multipoint results. They showed the following maximum multipoint H-LOD (heterogeneity LOD) and NPL scores: 2p14-13: 0.89 (P = 0.06) and 2.08 (P = 0.02), 4q24-32: 1.84 (P = 0.007) and 1.67 (P = 0.03), 5q21-31: 2.88 (P= 0.0007), and 2.65 (P = 0.002), 6p25-24: 2.13 (P = 0.005) and 3.59 (P = 0.0005), 6p23: 2.42 (P = 0.001) and 3.07 (P = 0.001), 8p22-21: 1.57 (P = 0.01) and 2.56 (P = 0.005), 10p15-11: 2.04 (P = 0.005) and 1.78 (P = 0.03). The degree of 'internal replication' across subsets differed, with 5q, 6p, and 8p being most consistent and 2p and 10p being least consistent. On 6p, the data suggested the presence of two susceptibility genes, in 6p25-24 and 6p23-22. Very few families were positive on more than one region, and little correlation between regions was evident, suggesting substantial locus heterogeneity. The levels of statistical significance were modest, as expected from loci contributing to complex traits. However, our internal replications, when considered along with the positive results obtained in multiple other samples, suggests that most of these six regions are likely to contain genes that influence liability to schizophrenia.

Genome scan meta-analysis of schizophrenia and bipolar disorder, part II:Schizophrenia

Schizophrenia is a common disorder with high heritability and a 10-fold increase in risk to siblings of probands. Replication has been inconsistent for reports of significant genetic linkage. To assess evidence for linkage across studies, rank-based genome scan meta-analysis (GSMA) was applied to data from 20 schizophrenia genome scans. Each marker for each scan was assigned to 1 of 120 30-cM bins, with the bins ranked by linkage scores (1 = most significant) and the ranks averaged across studies (R(avg)) and then weighted for sample size (N(sqrt)[affected casess]). A permutation test was used to compute the probability of observing, by chance, each bin's average rank (P(AvgRnk)) or of observing it for a bin with the same place (first, second, etc.) in the order of average ranks in each permutation (P(ord)). The GSMA produced significant genomewide evidence for linkage on chromosome 2q (PAvgRnk

Extensive genome duplications in sturgeons:New evidence from microsatellite data

Summary: Genome duplications and polyploidization events are thought to have played relevant roles in the early stages of vertebrate evolution, in particular near the time of divergence of the lamprey lineage. Additional genome duplications, specifically in ray-finned fish, may have occurred before the divergence of the teleosts. The role of polyploidization in vertebrate genome evolution is a thriving area of research. Sturgeons (order Acipenseriformes) provide a unique model for the investigation of genome duplication, with existing species possessing 120, 250 or 360 chromosomes. In the present study, data from 240 sturgeon specimens representing 11 species were used for analysis of ploidy levels. Allele numbers were assessed at eleven microsatellite loci. The results provide further evidence for functional diploidy, tetraploidy and hexaploidy in species possessing 120, 250 and 360 chromosomes, respectively. The analysis also uncovered novel evidence for functional hexaploidy in the shortnose sturgeon (Acipenser brevirostrum). In conclusion, the process of functional genome reduction is demonstrated to be an on-going process in this fish lineage. © 2013 Blackwell Verlag GmbH.

Genome-wide association meta-analysis of human longevity identifies a novel locus conferring survival beyond 90 years of age

The genetic contribution to the variation in human lifespan is approximately 25%.  Despite the large number of identified disease-susceptibility loci, it is not known which loci influence population mortality.  We performed a genome-wide association meta-analysis of 7729 long-lived individuals of European descent (≥ 85 years) and 16121 younger controls (< 65 years) followed by replication in an additional set of 13060 long-lived individuals and 61156 controls. In addition, we performed a subset analysis in cases ≥ 90 years. We observed genome-wide significant association with longevity, as reflected by survival to ages beyond 90 years, at a novel locus, rs2149954, on chromosome 5q33.3 (OR = 1.10, P =1.74 x 10-8). We also confirmed association of rs4420638 on chromosome 19q13.32 (OR = 0.72, P = 3.40 x 10-36), representing the TOMM40/APOE/APOC1 locus. In a prospective meta-analysis (n = 34103) the minor allele of rs2149954 (T) on chromosome 5q33.3 associates with increased survival (HR = 0.95, P = 0.003). This allele has previously been reported to associate with low blood pressure in middle age. Interestingly, the minor allele (T) associates with decreased cardiovascular mortality risk, independent of blood pressure. We report on the first GWAS-identified longevity locus on chromosome 5q33.3 influencing survival in the general European population. The minor allele of this locus associates with low blood pressure in middle age, although the contribution of this allele to survival may be less dependent on blood pressure. Hence, the pleiotropic mechanisms by which this intragenic variation contributes to lifespan regulation have to be elucidated.

Genome-wide association study of urinary albumin excretion rate in patients with type 1 diabetes

Aims/hypothesis: An abnormal urinary albumin excretion rate (AER) is often the first clinically detectable manifestation of diabetic nephropathy. Our aim was to estimate the heritability and to detect genetic variation associated with elevated AER in patients with type 1 diabetes. 
 Methods: The discovery phase genome-wide association study (GWAS) included 1,925 patients with type 1 diabetes and with data on 24 h AER. AER was analysed as a continuous trait and the analysis was stratified by the use of antihypertensive medication. Signals with a p value <10^{−4 were followed up in 3,750 additional patients with type 1 diabetes from seven studies.} 
 Results: The narrow-sense heritability, captured with our genotyping platform, was estimated to explain 27.3% of the total AER variability, and 37.6% after adjustment for covariates. In the discovery stage, five single nucleotide polymorphisms in the GLRA3 gene were strongly associated with albuminuria (p < 5 × 10^{−8). In the replication group, a nominally significant association (p = 0.035) was observed between albuminuria and rs1564939 in GLRA3, but this was in the opposite direction. Sequencing of the surrounding genetic region in 48 Finnish and 48 UK individuals supported the possibility that population-specific rare variants contribute to the synthetic association observed at the common variants in GLRA3. The strongest replication (p = 0.026) was obtained for rs2410601 between the PSD3 and SH2D4A genes. Pathway analysis highlighted natural killer cell mediated immunity processes.} 
 Conclusions/interpretation: This study suggests novel pathways and molecular mechanisms for the pathogenesis of albuminuria in type 1 diabetes.

No evidence that runs of homozygosity are associated with schizophrenia in an Irish genome-wide association dataset

Runs of homozygosity (ROH), regions of the genome containing many consecutive homozygous SNPs, may represent two copies of a haplotype inherited from a common ancestor. A rare variant on this haplotype could thus be present in a homozygous and potentially recessive state. To detect rare risk variants for schizophrenia, we performed an ROH analysis in a homogeneous Irish genome wide association study (GWAS) dataset consisting of 1606 cases and 1794 controls. There was no genome-wide excess of ROH in cases compared to controls (p=0.7986). No consensus ROH at individual loci showed association with schizophrenia after genome-wide correction.

A genome-wide association study identifies new susceptibility loci for esophageal adenocarcinoma and Barrett's esophagus

Esophageal adenocarcinoma is a cancer with rising incidence and poor survival. Most such cancers arise in a specialized intestinal metaplastic epithelium, which is diagnostic of Barrett's esophagus. In a genome-wide association study, we compared esophageal adenocarcinoma cases (n = 2,390) and individuals with precancerous Barrett's esophagus (n = 3,175) with 10,120 controls in 2 phases. For the combined case group, we identified three new associations. The first is at 19p13 (rs10419226: P = 3.6 × 10(-10)) in CRTC1 (encoding CREB-regulated transcription coactivator), whose aberrant activation has been associated with oncogenic activity. A second is at 9q22 (rs11789015: P = 1.0 × 10(-9)) in BARX1, which encodes a transcription factor important in esophageal specification. A third is at 3p14 (rs2687201: P = 5.5 × 10(-9)) near the transcription factor FOXP1, which regulates esophageal development. We also refine a previously reported association with Barrett's esophagus near the putative tumor suppressor gene FOXF1 at 16q24 and extend our findings to now include esophageal adenocarcinoma.

Relationship between genome and epigenome - challenges and requirements for future research

Understanding the links between genetic, epigenetic and non-genetic factors throughout the lifespan and across generations and their role in disease susceptibility and disease progression offer entirely new avenues and solutions to major problems in our society. To overcome the numerous challenges, we have come up with nine major conclusions to set the vision for future policies and research agendas at the European level.

Genome-wide meta-analysis identifies 11 new loci for anthropometric traits and provides insights into genetic architecture

Approaches exploiting trait distribution extremes may be used to identify loci associated with common traits, but it is unknown whether these loci are generalizable to the broader population. In a genome-wide search for loci associated with the upper versus the lower 5th percentiles of body mass index, height and waist-to-hip ratio, as well as clinical classes of obesity, including up to 263,407 individuals of European ancestry, we identified 4 new loci (IGFBP4, H6PD, RSRC1 and PPP2R2A) influencing height detected in the distribution tails and 7 new loci (HNF4G, RPTOR, GNAT2, MRPS33P4, ADCY9, HS6ST3 and ZZZ3) for clinical classes of obesity. Further, we find a large overlap in genetic structure and the distribution of variants between traits based on extremes and the general population and little etiological heterogeneity between obesity subgroups.

Genome-wide association study identifies novel loci associated with resistance to bovine tuberculosis

Tuberculosis (TB) caused by Mycobacterium bovis is a re-emerging disease of livestock that is of major economic importance worldwide, as well as being a zoonotic risk there is significant heritability for host resistance to bovine TB (bTB) in dairy cattle. To identify resistance loci for bTB, we undertook a genome-wide association study in female Holstein-Friesian cattle with 592 cases and 559 age-matched controls from case herds. Cases and controls were categorised into distinct phenotypes: skin test and lesion positive vs skin test negative on multiple occasions, respectively these animals were genotyped with the Illumina BovineHD 700K BeadChip. Genome-wide rapid association using linear and logistic mixed models and regression (GRAMMAR), regional heritability mapping (RHM) and haplotype-sharing analysis identified two novel resistance loci that attained chromosome-wise significance, protein tyrosine phosphatase receptor T (PTPRT; P=4.8 × 10 -7) and myosin IIIB (MYO3B; P=5.4 × 10 -6). We estimated that 21% of the phenotypic variance in TB resistance could be explained by all of the informative single-nucleotide polymorphisms, of which the region encompassing the PTPRT gene accounted for 6.2% of the variance and a further 3.6% was associated with a putative copy number variant in MYO3B the results from this study add to our understanding of variation in host control of infection and suggest that genetic marker-based selection for resistance to bTB has the potential to make a significant contribution to bTB control.

Genome-wide significant associations in schizophrenia to ITIH3/4, CACNA1C and SDCCAG8, and extensive replication of associations reported by the Schizophrenia PGC

The Schizophrenia Psychiatric Genome-Wide Association Study Consortium (PGC) highlighted 81 single-nucleotide polymorphisms (SNPs) with moderate evidence for association to schizophrenia. After follow-up in independent samples, seven loci attained genome-wide significance (GWS), but multi-locus tests suggested some SNPs that did not do so represented true associations. We tested 78 of the 81 SNPs in 2640 individuals with a clinical diagnosis of schizophrenia attending a clozapine clinic (CLOZUK), 2504 cases with a research diagnosis of bipolar disorder, and 2878 controls. In CLOZUK, we obtained significant replication to the PGC-associated allele for no fewer than 37 (47%) of the SNPs, including many prior GWS major histocompatibility complex (MHC) SNPs as well as 3/6 non-MHC SNPs for which we had data that were reported as GWS by the PGC. After combining the new schizophrenia data with those of the PGC, variants at three loci (ITIH3/4, CACNA1C and SDCCAG8) that had not previously been GWS in schizophrenia attained that level of support. In bipolar disorder, we also obtained significant evidence for association for 21% of the alleles that had been associated with schizophrenia in the PGC. Our study independently confirms association to three loci previously reported to be GWS in schizophrenia, and identifies the first GWS evidence in schizophrenia for a further three loci. Given the number of independent replications and the power of our sample, we estimate 98% (confidence interval (CI) 78-100%) of the original set of 78 SNPs represent true associations. We also provide strong evidence for overlap in genetic risk between schizophrenia and bipolar disorder.

Whole-mitochondrial genome sequencing in primary open-angle glaucoma using massively parallel sequencing identifies novel and known pathogenic variants

Purpose:The aim of this study was to determine whether mutations in mitochondrial DNA play a role in high-pressure primary open-angle glaucoma (OMIM 137760) by analyzing new data from massively parallel sequencing of mitochondrial DNA.
Methods:Glaucoma patients with high-tension primary open-angle glaucoma and ethnically matched and age-matched control subjects without glaucoma were recruited. The entire human mitochondrial genome was amplified in two overlapping fragments by long-range polymerase chain reaction and used as a template for massively parallel sequencing on an Ion Torrent Personal Genome Machine. All variants were confirmed by conventional Sanger sequencing.
Results:Whole-mitochondrial genome sequencing was performed in 32 patients with primary open-angle glaucoma from India (n = 16) and Ireland (n = 16). In 16 of the 32 patients with primary open-angle glaucoma (50% of cases), there were 22 mitochondrial DNA mutations consisting of 7 novel mutations and 8 previously reported disease-associated sequence variants. Eight of 22 (36.4%) of the mitochondrial DNA mutations were in complex I mitochondrial genes.
Conclusion:Massively parallel sequencing using the Ion Torrent Personal Genome Machine with confirmation by Sanger sequencing detected a pathogenic mitochondrial DNA mutation in 50% of the primary open-angle glaucoma cohort. Our findings support the emerging concept that mitochondrial dysfunction results in the development of glaucoma and, more specifically, that complex I defects play a significant role in primary open-angle glaucoma pathogenesis.