Functional profiling of long intergenic non-coding RNAs in fission yeast.

Eukaryotic genomes express numerous long intergenic non-coding RNAs (lincRNAs) that do not overlap any coding genes. Some lincRNAs function in various aspects of gene regulation, but it is not clear in general to what extent lincRNAs contribute to the information flow from genotype to phenotype. To explore this question, we systematically analysed cellular roles of lincRNAs in Schizosaccharomyces pombe. Using seamless CRISPR/Cas9-based genome editing, we deleted 141 lincRNA genes to broadly phenotype these mutants, together with 238 diverse coding-gene mutants for functional context. We applied high-throughput colony-based assays to determine mutant growth and viability in benign conditions and in response to 145 different nutrient, drug, and stress conditions. These analyses uncovered phenotypes for 47.5% of the lincRNAs and 96% of the protein-coding genes. For 110 lincRNA mutants, we also performed high-throughput microscopy and flow cytometry assays, linking 37% of these lincRNAs with cell-size and/or cell-cycle control. With all assays combined, we detected phenotypes for 84 (59.6%) of all lincRNA deletion mutants tested. For complementary functional inference, we analysed colony growth of strains ectopically overexpressing 113 lincRNA genes under 47 different conditions. Of these overexpression strains, 102 (90.3%) showed altered growth under certain conditions. Clustering analyses provided further functional clues and relationships for some of the lincRNAs. These rich phenomics datasets associate lincRNA mutants with hundreds of phenotypes, indicating that most of the lincRNAs analysed exert cellular functions in specific environmental or physiological contexts. This study provides groundwork to further dissect the roles of these lincRNAs in the relevant conditions.

Recombination and biased segregation of mitochondrial genomes during crossing and meiosis of different Schizosaccharomyces pombe strains.

During meiosis, tethering of parental mitochondria to opposite cell poles inhibits the mixing of mitochondria with different genomes and ensures uniparental inheritance in thestandard laboratory strain of fission yeast. We here investigate mitochondrial inheritance in crosses between natural isolates using tetrad dissection and next-generation sequencing. We find that colonies grown from single spores can sometimes carry a mix of mitochondrial genotypes, that mitochondrial genomes can recombine during meiosis, that in some cases tetrads do not follow the 2:2 segregation pattern, and that certain crosses may feature a weak bias towards one of the parents. Together, these findings paint a more nuanced picture of mitochondrial inheritance in the wild.

Longevity is determined by ETS transcription factors in multiple tissues and diverse species.

Ageing populations pose one of the main public health crises of our time. Reprogramming gene expression by altering the activities of sequence-specific transcription factors (TFs) can ameliorate deleterious effects of age. Here we explore how a circuit of TFs coordinates pro-longevity transcriptional outcomes, which reveals a multi-tissue and multi-species role for an entire protein family: the E-twenty-six (ETS) TFs. In Drosophila, reduced insulin/IGF signalling (IIS) extends lifespan by coordinating activation of Aop, an ETS transcriptional repressor, and Foxo, a Forkhead transcriptional activator. Aop and Foxo bind the same genomic loci, and we show that, individually, they effect similar transcriptional programmes in vivo. In combination, Aop can both moderate or synergise with Foxo, dependent on promoter context. Moreover, Foxo and Aop oppose the gene-regulatory activity of Pnt, an ETS transcriptional activator. Directly knocking down Pnt recapitulates aspects of the Aop/Foxo transcriptional programme and is sufficient to extend lifespan. The lifespan-limiting role of Pnt appears to be balanced by a requirement for metabolic regulation in young flies, in which the Aop-Pnt-Foxo circuit determines expression of metabolic genes, and Pnt regulates lipolysis and responses to nutrient stress. Molecular functions are often conserved amongst ETS TFs, prompting us to examine whether other Drosophila ETS-coding genes may also affect ageing. We show that five out of eight Drosophila ETS TFs play a role in fly ageing, acting from a range of organs and cells including the intestine, adipose and neurons. We expand the repertoire of lifespan-limiting ETS TFs in C. elegans, confirming their conserved function in ageing and revealing that the roles of ETS TFs in physiology and lifespan are conserved throughout the family, both within and between species.

Fitness Landscape of the Fission Yeast Genome.

The relationship between DNA sequence, biochemical function, and molecular evolution is relatively well-described for protein-coding regions of genomes, but far less clear in noncoding regions, particularly, in eukaryote genomes. In part, this is because we lack a complete description of the essential noncoding elements in a eukaryote genome. To contribute to this challenge, we used saturating transposon mutagenesis to interrogate the Schizosaccharomyces pombe genome. We generated 31 million transposon insertions, a theoretical coverage of 2.4 insertions per genomic site. We applied a five-state hidden Markov model (HMM) to distinguish insertion-depleted regions from insertion biases. Both raw insertion-density and HMM-defined fitness estimates showed significant quantitative relationships to gene knockout fitness, genetic diversity, divergence, and expected functional regions based on transcription and gene annotations. Through several analyses, we conclude that transposon insertions produced fitness effects in 66-90% of the genome, including substantial portions of the noncoding regions. Based on the HMM, we estimate that 10% of the insertion depleted sites in the genome showed no signal of conservation between species and were weakly transcribed, demonstrating limitations of comparative genomics and transcriptomics to detect functional units. In this species, 3'- and 5'-untranslated regions were the most prominent insertion-depleted regions that were not represented in measures of constraint from comparative genomics. We conclude that the combination of transposon mutagenesis, evolutionary, and biochemical data can provide new insights into the relationship between genome function and molecular evolution.

Correction: Signatures of inflammation and impending multiple organ dysfunction in the hyperacute phase of trauma: A prospective cohort study.

[This corrects the article DOI: 10.1371/journal.pmed.1002352.].

Vitamin D receptor genotype influences risk of upper respiratory infection.

SNP in the vitamin D receptor (VDR) gene is associated with risk of lower respiratory infections. The influence of genetic variation in the vitamin D pathway resulting in susceptibility to upper respiratory infections (URI) has not been investigated. We evaluated the influence of thirty-three SNP in eleven vitamin D pathway genes (DBP, DHCR7, RXRA, CYP2R1, CYP27B1, CYP24A1, CYP3A4, CYP27A1, LRP2, CUBN and VDR) resulting in URI risk in 725 adults in London, UK, using an additive model with adjustment for potential confounders and correction for multiple comparisons. Significant associations in this cohort were investigated in a validation cohort of 737 children in Manchester, UK. In all, three SNP in VDR (rs4334089, rs11568820 and rs7970314) and one SNP in CYP3A4 (rs2740574) were associated with risk of URI in the discovery cohort after adjusting for potential confounders and correcting for multiple comparisons (adjusted incidence rate ratio per additional minor allele ≥1·15, P for trend ≤0·030). This association was replicated for rs4334089 in the validation cohort (P for trend=0·048) but not for rs11568820, rs7970314 or rs2740574. Carriage of the minor allele of the rs4334089 SNP in VDR was associated with increased susceptibility to URI in children and adult cohorts in the United Kingdom.

Signatures of inflammation and impending multiple organ dysfunction in the hyperacute phase of trauma: A prospective cohort study.

BACKGROUND: Severe trauma induces a widespread response of the immune system. This "genomic storm" can lead to poor outcomes, including Multiple Organ Dysfunction Syndrome (MODS). MODS carries a high mortality and morbidity rate and adversely affects long-term health outcomes. Contemporary management of MODS is entirely supportive, and no specific therapeutics have been shown to be effective in reducing incidence or severity. The pathogenesis of MODS remains unclear, and several models are proposed, such as excessive inflammation, a second-hit insult, or an imbalance between pro- and anti-inflammatory pathways. We postulated that the hyperacute window after trauma may hold the key to understanding how the genomic storm is initiated and may lead to a new understanding of the pathogenesis of MODS. METHODS AND FINDINGS: We performed whole blood transcriptome and flow cytometry analyses on a total of 70 critically injured patients (Injury Severity Score [ISS] ≥ 25) at The Royal London Hospital in the hyperacute time period within 2 hours of injury. We compared transcriptome findings in 36 critically injured patients with those of 6 patients with minor injuries (ISS ≤ 4). We then performed flow cytometry analyses in 34 critically injured patients and compared findings with those of 9 healthy volunteers. Immediately after injury, only 1,239 gene transcripts (4%) were differentially expressed in critically injured patients. By 24 hours after injury, 6,294 transcripts (21%) were differentially expressed compared to the hyperacute window. Only 202 (16%) genes differentially expressed in the hyperacute window were still expressed in the same direction at 24 hours postinjury. Pathway analysis showed principally up-regulation of pattern recognition and innate inflammatory pathways, with down-regulation of adaptive responses. Immune deconvolution, flow cytometry, and modular analysis suggested a central role for neutrophils and Natural Killer (NK) cells, with underexpression of T- and B cell responses. In the transcriptome cohort, 20 critically injured patients later developed MODS. Compared with the 16 patients who did not develop MODS (NoMODS), maximal differential expression was seen within the hyperacute window. In MODS versus NoMODS, 363 genes were differentially expressed on admission, compared to only 33 at 24 hours postinjury. MODS transcripts differentially expressed in the hyperacute window showed enrichment among diseases and biological functions associated with cell survival and organismal death rather than inflammatory pathways. There was differential up-regulation of NK cell signalling pathways and markers in patients who would later develop MODS, with down-regulation of neutrophil deconvolution markers. This study is limited by its sample size, precluding more detailed analyses of drivers of the hyperacute response and different MODS phenotypes, and requires validation in other critically injured cohorts. CONCLUSIONS: In this study, we showed how the hyperacute postinjury time window contained a focused, specific signature of the response to critical injury that led to widespread genomic activation. A transcriptomic signature for later development of MODS was present in this hyperacute window; it showed a strong signal for cell death and survival pathways and implicated NK cells and neutrophil populations in this differential response.

Transient structural variations have strong effects on quantitative traits and reproductive isolation in fission yeast.

Large structural variations (SVs) within genomes are more challenging to identify than smaller genetic variants but may substantially contribute to phenotypic diversity and evolution. We analyse the effects of SVs on gene expression, quantitative traits and intrinsic reproductive isolation in the yeast Schizosaccharomyces pombe. We establish a high-quality curated catalogue of SVs in the genomes of a worldwide library of S. pombe strains, including duplications, deletions, inversions and translocations. We show that copy number variants (CNVs) show a variety of genetic signals consistent with rapid turnover. These transient CNVs produce stoichiometric effects on gene expression both within and outside the duplicated regions. CNVs make substantial contributions to quantitative traits, most notably intracellular amino acid concentrations, growth under stress and sugar utilization in winemaking, whereas rearrangements are strongly associated with reproductive isolation. Collectively, these findings have broad implications for evolution and for our understanding of quantitative traits including complex human diseases.

A CRISPR/Cas9-based method and primer design tool for seamless genome editing in fission yeast.

In the fission yeast Schizosaccharomyces pombe the prevailing approach for gene manipulations is based on homologous recombination of a PCR product that contains genomic target sequences and a selectable marker. The CRISPR/Cas9 system has recently been implemented in fission yeast, which allows for seamless genome editing without integration of a selection marker or leaving any other genomic 'scars'. The published method involves manual design of the single guide RNA (sgRNA), and digestion of a large plasmid with a problematic restriction enzyme to clone the sgRNA. To increase the efficiency of this approach, we have established and optimized a PCR-based system to clone the sgRNA without restriction enzymes into a plasmid with a dominant natMX6 (nourseothricin) selection marker. We also provide a web-tool, CRISPR4P, to support the design of the sgRNAs and the primers required for the entire process of seamless DNA deletion. Moreover, we report the preparation of G1-synchronized and cryopreserved S. pombe cells, which greatly increases the efficiency and speed for transformations, and may also facilitate standard gene manipulations. Applying this optimized CRISPR/Cas9-based approach, we have successfully deleted over 80 different non-coding RNA genes, which are generally lowly expressed, and have inserted 7 point mutations in 4 different genomic regions.

Double-blind randomised placebo-controlled trial of bolus-dose vitamin D3 supplementation in adults with asthma (ViDiAs).

RATIONALE: Asthma exacerbations are commonly precipitated by viral upper respiratory infections (URIs). Vitamin D insufficiency associates with susceptibility to URI in patients with asthma. Trials of vitamin D in adults with asthma with incidence of exacerbation and URI as primary outcome are lacking. OBJECTIVE: To conduct a randomised controlled trial of vitamin D3 supplementation for the prevention of asthma exacerbation and URI (coprimary outcomes). MEASUREMENTS AND METHODS: 250 adults with asthma in London, UK were allocated to receive six 2-monthly oral doses of 3 mg vitamin D3 (n=125) or placebo (n=125) over 1 year. Secondary outcomes included asthma control test and St George's Respiratory Questionnaire scores, fractional exhaled nitric oxide and concentrations of inflammatory markers in induced sputum. Subgroup analyses were performed to determine whether effects of supplementation were modified by baseline vitamin D status or genotype for 34 single nucleotide polymorphisms in 11 vitamin D pathway genes. MAIN RESULTS: 206/250 participants (82%) were vitamin D insufficient at baseline. Vitamin D3 did not influence time to first severe exacerbation (adjusted HR 1.02, 95% CI 0.69 to 1.53, p=0.91) or first URI (adjusted HR 0.87, 95% CI 0.64 to 1.16, p=0.34). No clinically important effect of vitamin D3 was seen on any of the secondary outcomes listed above. The influence of vitamin D3 on coprimary outcomes was not modified by baseline vitamin D status or genotype. CONCLUSIONS: Bolus-dose vitamin D3 supplementation did not influence time to exacerbation or URI in a population of adults with asthma with a high prevalence of baseline vitamin D insufficiency. TRIAL REGISTRATION NUMBER: NCT00978315 (ClinicalTrials.gov).

Polymorphisms in the WNK1 gene are associated with blood pressure variation and urinary potassium excretion.

WNK1--a serine/threonine kinase involved in electrolyte homeostasis and blood pressure (BP) control--is an excellent candidate gene for essential hypertension (EH). We and others have previously reported association between WNK1 and BP variation. Using tag SNPs (tSNPs) that capture 100% of common WNK1 variation in HapMap, we aimed to replicate our findings with BP and to test for association with phenotypes relating to WNK1 function in the British Genetics of Hypertension (BRIGHT) study case-control resource (1700 hypertensive cases and 1700 normotensive controls). We found multiple variants to be associated with systolic blood pressure, SBP (7/28 tSNPs min-p = 0.0005), diastolic blood pressure, DBP (7/28 tSNPs min-p = 0.002) and 24 hour urinary potassium excretion (10/28 tSNPs min-p = 0.0004). Associations with SBP and urine potassium remained significant after correction for multiple testing (p = 0.02 and p = 0.01 respectively). The major allele (A) of rs765250, located in intron 1, demonstrated the strongest evidence for association with SBP, effect size 3.14 mmHg (95%CI:1.23-4.9), DBP 1.9 mmHg (95%CI:0.7-3.2) and hypertension, odds ratio (OR: 1.3 [95%CI: 1.0-1.7]).We genotyped this variant in six independent populations (n = 14,451) and replicated the association between rs765250 and SBP in a meta-analysis (p = 7 x 10(-3), combined with BRIGHT data-set p = 2 x 10(-4), n = 17,851). The associations of WNK1 with DBP and EH were not confirmed. Haplotype analysis revealed striking associations with hypertension and BP variation (global permutation p<10(-7)). We identified several common haplotypes to be associated with increased BP and multiple low frequency haplotypes significantly associated with lower BP (>10 mmHg reduction) and risk for hypertension (OR<0.60). Our data indicates that multiple rare and common WNK1 variants contribute to BP variation and hypertension, and provide compelling evidence to initiate further genetic and functional studies to explore the role of WNK1 in BP regulation and EH.