BackgroundEvery year the human population encounters epidemic outbreaks of influenza, and history reveals recurring pandemics that have had devastating consequences. The current work focuses on the development of a robust algorithm for detecting influenza strains that have a composite genomic architecture. These influenza subtypes can be generated through a reassortment process, whereby a virus can inherit gene segments from two different types of influenza particles during replication. Reassortant strains are often not immediately recognised by the adaptive immune system of the hosts and hence may be the source of pandemic outbreaks. Owing to their importance in public health and their infectious ability, it is essential to identify reassortant influenza strains in order to understand the evolution of this virus and describe reassortment pathways that may be biased towards particular viral segments. Phylogenetic methods have been used traditionally to identify reassortant viruses. In many studies up to now, the assumption has been that if two phylogenetic trees differ, it is because reassortment has caused them to be different. While phylogenetic incongruence may be caused by real differences in evolutionary history, it can also be the result of phylogenetic error. Therefore, we wish to develop a method for distinguishing between topological inconsistency that is due to confounding effects and topological inconsistency that is due to reassortment.ResultsThe current work describes the implementation of two approaches for robustly identifying reassortment events. The algorithms rest on the idea of significance of difference between phylogenetic trees or phylogenetic tree sets, and subtree pruning and regrafting operations, which mimic the effect of reassortment on tree topologies. The first method is based on a maximum likelihood (ML) framework (MLreassort) and the second implements a Bayesian approach (Breassort) for reassortment detection. We focus on reassortment events that are found by both methods. We test both methods on a simulated dataset and on a small collection of real viral data isolated in Hong Kong in 1999.ConclusionsThe nature of segmented viral genomes present many challenges with respect to disease. The algorithms developed here can effectively identify reassortment events in small viral datasets and can be applied not only to influenza but also to other segmented viruses. Owing to computational demands of comparing tree topologies, further development in this area is necessary to allow their application to larger datasets.

Colorectal cancer (CRC) is a leading cause of cancer-related death worldwide, and has a strong heritable basis. We report a genome-wide association analysis of 34,627 CRC cases and 71,379 controls of European ancestry that identifies SNPs at 31 new CRC risk loci. We also identify eight independent risk SNPs at the new and previously reported European CRC loci, and a further nine CRC SNPs at loci previously only identified in Asian populations. We use in situ promoter capture Hi-C (CHi-C), gene expression, and in silico annotation methods to identify likely target genes of CRC SNPs. Whilst these new SNP associations implicate target genes that are enriched for known CRC pathways such as Wnt and BMP, they also highlight novel pathways with no prior links to colorectal tumourigenesis. These findings provide further insight into CRC susceptibility and enhance the prospects of applying genetic risk scores to personalised screening and prevention.

We aimed to develop an efficient, flexible and scalable approach to diagnostic genome-wide sequence analysis of genetically heterogeneous clinical presentations. Here we present G2P ( www.ebi.ac.uk/gene2phenotype ) as an online system to establish, curate and distribute datasets for diagnostic variant filtering via association of allelic requirement and mutational consequence at a defined locus with phenotypic terms, confidence level and evidence links. An extension to Ensembl Variant Effect Predictor (VEP), VEP-G2P was used to filter both disease-associated and control whole exome sequence (WES) with Developmental Disorders G2P (G2P DD ; 2044 entries). VEP-G2P DD shows a sensitivity/precision of 97.3%/33% for de novo and 81.6%/22.7% for inherited pathogenic genotypes respectively. Many of the missing genotypes are likely false-positive pathogenic assignments. The expected number and discriminative features of background genotypes are defined using control WES. Using only human genetic data VEP-G2P performs well compared to other freely-available diagnostic systems and future phenotypic matching capabilities should further enhance performance.

MicroRNAs (miRNAs) repress the expression levels of genes by binding to mRNA transcripts, acting as master regulators of cellular processes. Differential expression of miRNAs has been linked to virus-associated diseases involving members of the Hepacivirus , Herpesvirus , and Retrovirus families. In contrast, limited biological and molecular information has been reported on the potential role of cellular miRNAs in the life cycle of influenza A viruses (infA). In this study, we hypothesize that elucidating the miRNA expression signatures induced by low-pathogenicity swine-origin infA (S-OIV) pandemic H1N1 (2009) and highly pathogenic avian-origin infA (A-OIV) H7N7 (2003) infections could reveal temporal and strain-specific miRNA fingerprints during the viral life cycle, shedding important insights into the potential role of cellular miRNAs in host-infA interactions. Using a microfluidic microarray platform, we profiled cellular miRNA expression in human A549 cells infected with S- and A-OIVs at multiple time points during the viral life cycle, including global gene expression profiling during S-OIV infection. Using target prediction and pathway enrichment analyses, we identified the key cellular pathways associated with the differentially expressed miRNAs and predicted mRNA targets during infA infection, including the immune system, cell proliferation, apoptosis, cell cycle, and DNA replication and repair. By identifying the specific and dynamic molecular phenotypic changes (microRNAome) triggered by S- and A-OIV infection in human cells, we provide experimental evidence demonstrating a series of temporal and strain-specific host molecular responses involving different combinatorial contributions of multiple cellular miRNAs. Our results also identify novel potential exosomal miRNA biomarkers associated with pandemic S-OIV and deadly A-OIV-host infection.

ObjectiveWe assessed the effect of surgical resection of colorectal cancer (CRC) on perioperative plasma vitamin D (25OHD) and C-reactive protein (CRP) level. We investigated the relationship between circulating vitamin D level and CRC survival.DesignWe sequentially sampled 92 patients undergoing CRC resection, and measured plasma 25OHD and CRP. For survival analyses, we assayed 25OHD and CRP in two temporally distinct CRC patient cohorts (n=2006, n=2100) and investigated the association between survival outcome, circulating vitamin D and systemic inflammatory response.ResultsSerial sampling revealed a postoperative fall (mean 17.3 nmol/L; p=3.6e-9) in plasma 25OHD (nadir days 1–2). CRP peaked 3–5 days postoperatively (143.1 mg/L; p=1.4e-12), yet the postoperative fall in 25OHD was independent of CRP. In cohort analyses, 25OHD was lower in the 12 months following operation (mean=48.8 nmol/L) than preoperatively (54.8 nmol/L; p=1.2e-5) recovering after 24 months (52.2 nmol/L; p=0.002). Survival analysis in American Joint Committee on Cancer stages I–III demonstrated associations between 25OHD tertile and CRC mortality (HR=0.69; 95% CI 0.46 to 0.91) and all-cause mortality (HR=0.68; 95% CI 0.50 to 0.85), and was independent of CRP. We observed interaction effects between plasma 25OHD and rs11568820 genotype (functional VDR polymorphism) with a strong protective effect of higher 25OHD only in patients with GG genotype (HR=0.51; 95% CI 0.21 to 0.81). We developed an online tool for predicted survival (https://apps.igmm.ed.ac.uk/mortalityCalculator/) that incorporates 25OHD with clinically useful predictive performance (area under the curve 0.77).ConclusionsCRC surgery induces a fall in circulating 25OHD. Plasma 25OHD level is a prognostic biomarker with low 25OHD associated with poorer survival, particularly in those with rs11568820 GG genotype. A randomised trial of vitamin D supplementation after CRC surgery has compelling rationale.

BackgroundWhilst observational studies establish that lower plasma 25-hydroxyvitamin D (25-OHD) levels are associated with higher risk of colorectal cancer (CRC), establishing causality has proven challenging. Since vitamin D is modifiable, these observations have substantial clinical and public health implications. Indeed, many health agencies already recommend supplemental vitamin D. Here, we explore causality in a large Mendelian randomisation (MR) study using an improved genetic instrument for circulating 25-OHD.MethodsWe developed a weighted genetic score for circulating 25-OHD using six genetic variants that we recently reported to be associated with circulating 25-OHD in a large genome-wide association study (GWAS) meta-analysis. Using this score as instrumental variable in MR analyses, we sought to determine whether circulating 25-OHD is causally linked with CRC risk. We conducted MR analysis using individual-level data from 10,725 CRC cases and 30,794 controls (Scotland, UK Biobank and Croatia). We then applied estimates from meta-analysis of 11 GWAS of CRC risk (18,967 cases; 48,168 controls) in a summary statistics MR approach.ResultsThe new genetic score for 25-OHD was strongly associated with measured plasma 25-OHD levels in 2821 healthy Scottish controls (P = 1.47 × 10− 11), improving upon previous genetic instruments (F-statistic 46.0 vs. 13.0). However, individual-level MR revealed no association between 25-OHD score and CRC risk (OR 1.03/unit log-transformed circulating 25-OHD, 95% CI 0.51–2.07, P = 0.93). Similarly, we found no evidence for a causal relationship between 25-OHD and CRC risk using summary statistics MR analysis (OR 0.91, 95% CI 0.69–1.19, P = 0.48).ConclusionsDespite the scale of this study and employing an improved score capturing more of the genetic contribution to circulating 25-OHD, we found no evidence for a causal relationship between circulating 25-OHD and CRC risk. Although the magnitude of effect for vitamin D suggested by observational studies can confidently be excluded, smaller effects sizes and non-linear relationships remain plausible. Circulating vitamin D may be a CRC biomarker, but a causal effect on CRC risk remains unproven.Electronic supplementary materialThe online version of this article (10.1186/s12916-018-1119-2) contains supplementary material, which is available to authorized users.

This is a pre-copyedited, author-produced version of an article accepted for publication in Nature Genetics following peer review. The version of record "Deciphering colorectal cancer genetics through multi-omic analysis of 100,204 cases and 154,587 controls of European and east Asian ancestries"

Site-specific variation in colorectal cancer (CRC) incidence, biology and prognosis are poorly understood. We sought to determine whether common genetic variants influencing CRC risk might exhibit topographical differences on CRC risk through regional differences in effects on gene expression in the large bowel mucosa. We conducted a site-specific genetic association study (10 630 cases, 31 331 controls) to identify whether established risk variants exert differential effects on risk of proximal, compared to distal CRC. We collected normal colorectal mucosa and blood from 481 subjects and assessed mucosal gene expression using Illumina HumanHT-12v4 arrays in relation to germline genotype. Expression quantitative trait loci (eQTLs) were explored by anatomical location of sampling. The rs3087967 genotype (chr11q23.1 risk variant) exhibited significant site-specific effectsrisk of distal CRC (odds ratio [OR] = 1.20, P = 8.20 Â 10 À20 ) with negligible effects on proximal CRC risk (OR = 1.05, P = .10). Expression of 1261 genes differed between proximal and distal colonic mucosa (top hit PRAC gene, fold-difference = 10, P = 3.48 Â 10 À57 ). In eQTL studies, rs3087967 genotype was associated with expression of 8 cis-and 21 trans-genes. Four of these (AKAP14, ADH5P4, ASGR2, RP11-342M1.7) showed differential effects by site, with strongest trans-eQTL signals in proximal colonic mucosa (eg, AKAP14, beta = 0.61, P = 5.02 Â 10 À5 ) and opposite signals in distal mucosa (AKAP14, beta = À0.17, P = .04). In summary, genetic variation at the chr11q23.1 risk locus imparts greater risk of distal rather than proximal CRC and exhibits site-specific differences in eQTL effects in normal mucosa. Topographical differences in genomic control over gene expression relevant to CRC risk may underlie site-specific variation in CRC.Results may inform individualised CRC screening programmes.