A metagenomic comparison of endemic viruses from broiler chickens with runting stunting syndrome and from normal birds

Runting-stunting syndrome (RSS) in broiler chickens is an enteric disease that causes significant economic losses to poultry producers worldwide due to elevated feed conversion ratios, decreased body weight during growth, and excessive culling. Of specific interest are the viral agents associated with RSS which have been difficult to fully characterise to date. Past research into the aetiology of RSS has implicated a wide variety of RNA and DNA viruses however, to date, no individual virus has been identified as the main agent of RSS and the current opinion is that it may be caused by a community of viruses, collectively known as the virome. This paper attempts to characterise the viral pathogens associated with 2 – 3 week old RSS-affected and unaffected broiler chickens using next-generation sequencing and comparative metagenomics. Analysis of the viromes identified a total of 20 DNA & RNA viral families, along with 2 unidentified categories, comprised of 31 distinct viral genera and 7 unclassified genera. The most abundant viral families identified in this study were the Astroviridae, Caliciviridae, Picornaviridae, Parvoviridae, Coronaviridae, Siphoviridae, and Myoviridae. This study has identified historically significant viruses associated with the disease such as chicken astrovirus, avian nephritis virus, chicken parvovirus, and chicken calicivirus along with relatively novel viruses such as chicken megrivirus and sicinivirus 1 and will help expand the knowledge related to enteric disease in broiler chickens, provide insights into the viral constituents of a healthy avian gut, and identify a variety of enteric viruses and viral communities appropriate for further study.

Prior knowledge transfer across transcriptional data sets and technologies using compositional statistics yields new mislabelled ovarian cell line

Here, we describe gene expression compositional assignment (GECA), a powerful, yet simple method based on compositional statistics that can validate the transfer of prior knowledge, such as gene lists, into independent data sets, platforms and technologies. Transcriptional profiling has been used to derive gene lists that stratify patients into prognostic molecular subgroups and assess biomarker performance in the pre-clinical setting. Archived public data sets are an invaluable resource for subsequent in silico validation, though their use can lead to data integration issues. We show that GECA can be used without the need for normalising expression levels between data sets and can outperform rank-based correlation methods. To validate GECA, we demonstrate its success in the cross-platform transfer of gene lists in different domains including: bladder cancer staging, tumour site of origin and mislabelled cell lines. We also show its effectiveness in transferring an epithelial ovarian cancer prognostic gene signature across technologies, from a microarray to a next-generation sequencing setting. In a final case study, we predict the tumour site of origin and histopathology of epithelial ovarian cancer cell lines. In particular, we identify and validate the commonly-used cell line OVCAR-5 as non-ovarian, being gastrointestinal in origin. GECA is available as an open-source R package.

Efficient Genotyping of KRAS Mutant Non-Small Cell Lung Cancer Using a Multiplexed Droplet Digital PCR Approach.

Droplet digital PCR (ddPCR) can be used to detect low frequency mutations in oncogene-driven lung cancer. The range of KRAS point mutations observed in NSCLC necessitates a multiplex approach to efficient mutation detection in circulating DNA. Here we report the design and optimisation of three discriminatory ddPCR multiplex assays investigating nine different KRAS mutations using PrimePCR™ ddPCR™ Mutation Assays and the Bio-Rad QX100 system. Together these mutations account for 95% of the nucleotide changes found in KRAS in human cancer. Multiplex reactions were optimised on genomic DNA extracted from KRAS mutant cell lines and tested on DNA extracted from fixed tumour tissue from a cohort of lung cancer patients without prior knowledge of the specific KRAS genotype. The multiplex ddPCR assays had a limit of detection of better than 1 mutant KRAS molecule in 2,000 wild-type KRAS molecules, which compared favourably with a limit of detection of 1 in 50 for next generation sequencing and 1 in 10 for Sanger sequencing. Multiplex ddPCR assays thus provide a highly efficient methodology to identify KRAS mutations in lung adenocarcinoma.

Challenges in molecular testing in non-small-cell lung cancer patients with advanced disease

Lung cancer diagnostics have progressed greatly in the previous decade. Development of molecular testing to identify an increasing number of potentially clinically actionable genetic variants, using smaller samples obtained via minimally invasive techniques, is a huge challenge. Tumour heterogeneity and cancer evolution in response to therapy means that repeat biopsies or circulating biomarkers are likely to be increasingly useful to adapt treatment as resistance develops. We highlight some of the current challenges faced in clinical practice for molecular testing of EGFR, ALK, and new biomarkers such as PDL1. Implementation of next generation sequencing platforms for molecular diagnostics in non-small-cell lung cancer is increasingly common, allowing testing of multiple genetic variants from a single sample. The use of next generation sequencing to recruit for molecularly stratified clinical trials is discussed in the context of the UK Stratified Medicine Programme and The UK National Lung Matrix Trial.

Microbial phylogenetic and functional responses within acidified wastewater communities exhibiting enhanced phosphate uptake

Acid stimulated accumulation of insoluble phosphorus within microbial cells is highly beneficial to wastewater treatment but remains largely unexplored. Using single cell analyses and next generation sequencing, the response of active polyphosphate accumulating microbial communities under conditions of enhanced phosphorus uptake under both acidic and aerobic conditions was characterised. Phosphorus accumulation activities were highest under acidic conditions (pH 5.5 > 8.5), where a significant positive effect on bioaccumulation was observed at pH 5.5 when compared to pH 8.5. In contrast to the Betaproteobacteria and Actinobacteria dominated enhanced biological phosphorus removal process, the functionally active polyP accumulators at pH 5.5 belonged to the Gammaproteobacteria, with key accumulators identified as members of the families Aeromonadaceae and Enterobacteriaceae. This study demonstrated a significant enrichment of key polyphosphate kinase and exopolyphosphatase genes within the community metagenome after acidification, concomitant with an increase in P accumulation kinetics.

Exploration génétique de l'hypothyroïdie congénitale par dysgénésie thyroïdienne

L'hypothyroïdie congénitale par dysgénésie thyroïdienne (HCDT, ectopie dans plus de 80 %) a une prévalence de 1 cas sur 4000 naissances vivantes. L’HCDT est la conséquence d'une défaillance de la thyroïde embryonnaire à se différencier, à se maintenir ou à migrer vers sa localisation anatomique (partie antérieure du cou), qui aboutit à une absence totale de la thyroïde (athyréose) ou à une ectopie thyroïdienne (linguale ou sublinguale). Les HCDT sont principalement non-syndromiques (soit 98% des cas sont non-familiale), ont un taux de discordance de 92% chez les jumeaux monozygotes, et ont une prédominance féminine et ethnique (i.e., Caucasienne). La majorité des cas d’HCDT n’a pas de cause connue, mais est associée à un déficit sévère en hormones thyroïdiennes (hypothyroïdie). Des mutations germinales dans les facteurs de transcription liés à la thyroïde (NKX2.1, FOXE1, PAX8, NKX2.5) ont été identifiées dans seulement 3% des patients atteints d’HCDT sporadiques et l’analyse de liaisons exclue ces gènes dans les rares familles multiplex avec HCDT. Nous supposons que le manque de transmission familiale claire d’HCDT peut résulter de la nécessité d’au moins deux « hits » génétiques différents dans des gènes importants pour le développement thyroïdien. Pour répondre au mieux nos questions de recherche, nous avons utilisé deux approches différentes: 1) une approche gène candidat, FOXE1, seul gène impliqué dans l’ectopie dans le modèle murin et 2) une approche en utilisant les techniques de séquençage de nouvelle génération (NGS) afin de trouver des variants génétiques pouvant expliquer cette pathologie au sein d’une cohorte de patients avec HCDT. Pour la première approche, une étude cas-contrôles a été réalisée sur le promoteur de FOXE1. Il a récemment été découvert qu’une région du promoteur de FOXE1 est différentiellement méthylée au niveau de deux dinucléotides CpG consécutifs, définissant une zone cruciale de contrôle de l’expression de FOXE1. L’analyse d’association basée sur les haplotypes a révélé qu’un haplotype (Hap1: ACCCCCCdel1C) est associé avec le HCDT chez les Caucasiens (p = 5x10-03). Une réduction significative de l’activité luciférase est observée pour Hap1 (réduction de 68%, p<0.001) comparé au promoteur WT de FOXE1. Une réduction de 50% de l’expression de FOXE1 dans une lignée de cellules thyroïdienne humaine est suffisante pour réduire significativement la migration cellulaire (réduction de 55%, p<0.05). Un autre haplotype (Hap2: ACCCCCCC) est observé moins fréquemment chez les Afro-Américain comparés aux Caucasiens (p = 1.7x10-03) et Hap2 diminue l’activité luciférase (réduction de 26%, p<0.001). Deux haplotypes distincts sont trouvés fréquemment dans les contrôles Africains (Black-African descents). Le premier haplotype (Hap3: GTCCCAAC) est fréquent (30.2%) chez les contrôles Afro-Américains comparés aux contrôles Caucasiens (6.3%; p = 2.59 x 10-9) tandis que le second haplotype (Hap4: GTCCGCAC) est trouvé exclusivement chez les contrôles Afro-Américains (9.4%) et est absent chez les contrôles Caucasiens (P = 2.59 x 10-6). Pour la deuxième approche, le séquençage de l’exome de l’ADN leucocytaire entre les jumeaux MZ discordants n’a révélé aucune différence. D'où l'intérêt du projet de séquençage de l’ADN et l’ARN de thyroïdes ectopiques et orthotopiques dans lesquelles de l'expression monoallélique aléatoire dans a été observée, ce qui pourrait expliquer comment une mutation monoallélique peut avoir des conséquences pathogéniques. Finalement, le séquençage de l’exome d’une cohorte de 36 cas atteints d’HCDT a permis d’identifier de nouveaux variants probablement pathogéniques dans les gènes récurrents RYR3, SSPO, IKBKE et TNXB. Ces quatre gènes sont impliqués dans l’adhésion focale (jouant un rôle dans la migration cellulaire), suggérant un rôle direct dans les défauts de migration de la thyroïde. Les essais de migration montrent une forte diminution (au moins 60% à 5h) de la migration des cellules thyroïdiennes infectées par shRNA comparés au shCtrl dans 2 de ces gènes. Des zebrafish KO (-/- et +/-) pour ces nouveaux gènes seront réalisés afin d’évaluer leur impact sur l’embryologie de la thyroïde.

Genomics and virulence factors of Fusobacterium necrophorum

Fusobacterium necrophorum, a Gram negative, anaerobic bacterium, is a common cause of acute pharyngitis and tonsillitis and a rare cause of more severe infections of the head and neck. At the beginning of the project, there was no available genome sequence for F. necrophorum. The aim of this project was to sequence the F. necrophorum genome and identify and study its putative virulence factors contained using in silico and in vitro analysis. Type strains JCM 3718 and JCM 3724,F. necrophorum subspecies necrophorum (Fnn) and funduliforme (Fnf), respectively, and strain ARU 01 (Fnf), isolated from a patient with LS, were commercially sequenced by Roche 454 GS-FLX+ next generation sequencing and assembled into contigs using Roche GS Assembler. Sequence data was annotated semi-automatically, using the xBASE pipeline, BLASTp and Pfam. The F. necrophorum genome was determined to be approximately 2.1 – 2.3 Mb in size, with an estimated 1,950 ORFs and includes genes for a leukotoxin, ecotin, haemolysin, haemagglutinin, haemin receptor, adhesin and type Vb and Vc secretion systems. The prevalence of the leukotoxin gene was investigated in strains JCM 3718, JCM 3724 and ARU 01, as well as a clinical collection of 25 Fnf strains, identified using biochemical and molecular tests. The leukotoxin operon was found to be universal within the strain collection by PCR. HL-60 cells subjected to aliquots of concentrated high molecular weight culture supernatant, predicted to contain the secreted leukotoxins of strains JCM 3718, JCM 3724 and ARU 01, were killed in a dose-dependent manner. The cytotoxic effect of the leukotoxin against human donor white blood cells was also tested to validate the HL-60 assay. The differences in the results between the two assays were not statistically significant. Ecotin, a serine protease inhibitor, was found to be present in 100 % of the strain collection and had a highly conserved sequence with primary and secondary binding sites exposed on opposing sides of the protein. During enzyme inhibition studies, a purified recombinant F. necrophorum ecotin protein inhibited human neutrophil elastase, a protease that degrades bacteria at inflammation sites, and human plasma kallikrein, a component of the host clotting cascade. The recombinant ecotin also prolonged human plasma clotting times by up to 7-fold for the extrinsic pathway, and up to 40-fold for the intrinsic pathway. The genome sequence data provides important information about F. necrophorum type strains and enables comparative study between strains and subspecies. Results from the leukotoxin and ecotin assays can be used to build up an understanding of how the organism behaves during infection.

Role of fibroblast growth factors and their receptors in prostate cancer

Prostate cancer (PCa) is the most common non-cutaneous malignant disease among males in the developed countries. Radical prostatectomy (RP) is an effective therapy for most PCa patients with localized or locally invaded tumors but in some cases the cancer recurs after RP. PCa is a heterogeneous disease, which is regulated by many factors, such as androgen receptor (AR), estrogen receptors and  (ER and ER), fibroblast growth factors (FGFs) and their receptors (FGFRs). In this study, the role of ERβ, FGF8, FGF13 and FGFRL1 was investigated in PCa. Previous studies have suggested that ER is protective against PCa whereas FGF8 has been shown to induce PCa in transgenic mice. FGF13 and FGFRL1 are poorly understood members of the FGF and FGFR families, respectively. Transgenic mouse models were used to investigate the ability of inactivated ERβ to facilitate FGF8-induced prostate tumorigenesis. Human PCa tissue microarrays (TMAs) were used to study the expression pattern of FGF13 and FGFRL1 in PCa and the results were correlated to corresponding patient data. The targets and biological functions of FGF13 and FGFRL1 were characterized using experimental in vivo and in vitro models. The results show that deficiency of ERβ, which had been expected to have tumor suppressing capacity, seemed to influence epithelial differentiation but did not affect FGF8-induced prostate tumorigenesis. Analysis of the TMAs showed increased expression of FGF13 in PCa. The level of cytoplasmic FGF13 was associated with the PCa biochemical recurrence (BCR), demonstrated by increasing serum PSA value, and was able to act as an independent prognostic biomarker for PCa patients after RP. Expression of FGFRL1, the most recently identified FGFR, was also elevated in PCa. Cytoplasmic and nuclear FGFRL1 was associated with high Gleason score and Ki67 level whereas the opposite was true for the cell membrane FGFRL1. Silencing of FGFRL1 in PC-3M cells led to a strongly decreased growth rate of these cells as xenografts in nude mice and the experiments with PCa cell lines showed that FGFRL1 is able to modulate the FGF2- and FGF8-induced signaling pathways. The next generation sequencing (NGS) experiments with FGFRL1-silenced PC-3M cells revealed candidates for FGFRL1 target genes. In summary, these studies provide new data on the FGF/FGFR signaling pathways in normal and malignant prostate and suggest a potential role for FGF13 and FGFRL1 as novel prognostic markers for PCa patients. Keywords: FGF8, FGF13, FGFRL1, ERβ, prostate cancer, prognostic marker

High diversity of skin-associated bacterial communities of marine fishes is promoted by their high variability among body parts, individuals and species

Animal-associated microbiotas form complex communities, which are suspected to play crucial functions for their host fitness. However, the biodiversity of these communities, including their differences between host species and individuals, has been scarcely studied, especially in case of skin-associated communities. In addition, the intraindividual variability (i.e. between body parts) has never been assessed to date. The objective of this study was to characterize skin bacterial communities of two teleostean fish species, namely the European seabass (Dicentrarchus labrax) and gilthead seabream (Sparus aurata), using a high-throughput DNA sequencing method. In order to focus on intrinsic factors of host-associated bacterial community variability, individuals of the two species were raised in controlled conditions. Bacterial diversity was assessed using a set of four complementary indices, describing the taxonomic and phylogenetic facets of biodiversity and their respective composition (based on presence/absence data) and structure (based on species relative abundances) components. Variability of bacterial diversity was quantified at the interspecific, interindividual and intraindividual scales. We demonstrated that fish surfaces host highly diverse bacterial communities, whose composition was very different from that of surrounding bacterioplankton. This high total biodiversity of skin-associated communities was supported by the important variability, between host species, individuals and the different body parts (dorsal, anal, pectoral and caudal fins).

Bioinformatics analysis of the Pedobacter sp. NL19 genome

The last decades of the 20th century defined the genetic engineering advent, climaxing in the development of techniques, such as PCR and Sanger sequencing. This, permitted the appearance of new techniques to sequencing whole genomes, identified as next-generation sequencing. One of the many applications of these techniques is the in silico search for new secondary metabolites, synthesized by microorganisms exhibiting antimicrobial properties. The peptide antibiotics compounds can be classified in two classes, according to their biosynthesis, in ribosomal or nonribosomal peptides. Lanthipeptides are the most studied ribosomal peptides and are characterized by the presence of lanthionine and methylanthionine that result from posttranslational modifications. Lanthipeptides are divided in four classes, depending on their biosynthetic machinery. In class I, a LanB enzyme dehydrate serine and threonine residues in the C-terminus precursor peptide. Then, these residues undergo a cyclization step performed by a LanC enzyme, forming the lanthionine rings. The cleavage and the transport of the peptide is achieved by the LanP and LanT enzymes, respectively. Although, in class II only one enzyme, LanM, is responsible for the dehydration and cyclization steps and also only one enzyme performs the cleavage and transport, LanT. Pedobacter sp. NL19 is a Gram-negative bacterium, isolated from sludge of an abandon uranium mine, in Viseu (Portugal). Antibacterial activity in vitro was detected against several Gram-positive and Gram-negative bacteria. Sequencing and in silico analysis of NL19 genome revealed the presence of 21 biosynthetic clusters for secondary metabolites, including nonribosomal and ribosomal peptides biosynthetic clusters. Four lanthipeptides clusters were predicted, comprising the precursor peptides, the modifying enzymes (LanB and LanC), and also a bifunctional LanT. This result revealed the hybrid nature of the clusters, comprising characteristics from two distinct classes, which are poorly described in literature. The phylogenetic analysis of their enzymes showed that they clustered within the bacteroidetes clade. Furthermore, hybrid gene clusters were also found in other species of this phylum, revealing that it is a common characteristic in this group. Finally, the analysis of NL19 colonies by MALDI-TOF MS allowed the identification of a 3180 Da mass that corresponds to the predicted mass of a lanthipeptide encoded in one of the clusters. However, this result is not fully conclusive and further experiments are needed to understand the full potential of the compounds encoded in this type of clusters. In conclusion, it was determined that NL19 strain has the potential to produce diverse secondary metabolites, including lanthipeptides that were not functionally characterized so far.

Quasispecies dynamics and treatment outcome during early hepatitis C infection in a cohort of HIV-infected men

Hepatitis C virus (HCV) is emerging as one of the leading causes of morbidity and mortality in individuals infected with HIV and has overtaken AIDS-defining illnesses as a cause of death in HIV patient populations who have access to highly active antiretroviral therapy. For many years, the clonal analysis was the reference method for investigating viral diversity. In this thesis, a next generation sequencing (NGS) approach was developed using 454 pyrosequencing and Illumina-based technology. A sequencing pipeline was developed using two different NGS approaches, nested PCR, and metagenomics. The pipeline was used to study the viral populations in the sera of HCV-infected patients from a unique cohort of 160 HIV-positive patients with early HCV infection. These pipelines resulted in an improved understanding of HCV quasispecies dynamics, especially regarding studying response to treatment. Low viral diversity at baseline correlated with sustained virological response (SVR) while high viral diversity at baseline was associated with treatment failure. The emergence of new viral strains following treatment failure was most commonly associated with emerging dominance of pre-existing minority variants rather than re-infection. In the new era of direct-acting antivirals, next generation sequencing technologies are the most promising tool for identifying minority variants present in the HCV quasispecies populations at baseline. In this cohort, several mutations conferring resistance were detected in genotype 1a treatment-naïve patients. Further research into the impact of baseline HCV variants on SVR rates should be carried out in this population. A clearer understanding of the properties of viral quasispecies would enable clinicians to make improved treatment choices for their patients.

Ecology, physiology and performance in high-rate anaerobic digestion

The design demands on water and sanitation engineers are rapidly changing. The global population is set to rise from 7 billion to 10 billion by 2083. Urbanisation in developing regions is increasing at such a rate that a predicted 56% of the global population will live in an urban setting by 2025. Compounding these problems, the global water and energy crises are impacting the Global North and South alike. High-rate anaerobic digestion offers a low-cost, low-energy treatment alternative to the energy intensive aerobic technologies used today. Widespread implementation however is hindered by the lack of capacity to engineer high-rate anaerobic digestion for the treatment of complex wastes such as sewage. This thesis utilises the Expanded Granular Sludge Bed bioreactor (EGSB) as a model system in which to study the ecology, physiology and performance of high-rate anaerobic digestion of complex wastes. The impacts of a range of engineered parameters including reactor geometry, wastewater type, operating temperature and organic loading rate are systematically investigated using lab-scale EGSB bioreactors. Next generation sequencing of 16S amplicons is utilised as a means of monitoring microbial ecology. Microbial community physiology is monitored by means of specific methanogenic activity testing and a range of physical and chemical methods are applied to assess reactor performance. Finally, the limit state approach is trialled as a method for testing the EGSB and is proposed as a standard method for biotechnology testing enabling improved process control at full-scale. The arising data is assessed both qualitatively and quantitatively. Lab-scale reactor design is demonstrated to significantly influence the spatial distribution of the underlying ecology and community physiology in lab-scale reactors, a vital finding for both researchers and full-scale plant operators responsible for monitoring EGSB reactors. Recurrent trends in the data indicate that hydrogenotrophic methanogenesis dominates in high-rate anaerobic digestion at both full- and lab-scale when subject to engineered or operational stresses including low-temperature and variable feeding regimes. This is of relevance for those seeking to define new directions in fundamental understanding of syntrophic and competitive relations in methanogenic communities and also to design engineers in determining operating parameters for full-scale digesters. The adoption of the limit state approach enabled identification of biological indicators providing early warning of failure under high-solids loading, a vital insight for those currently working empirically towards the development of new biotechnologies at lab-scale.

'Abnormal vertical growth' : A disorder threatening the viability of the Australian macadamia industry

'Abnormal vertical growth' (AVG) was recognised in Australia as a dysfunction of macadamia (Macadamia spp.) in the mid-1990s. Affected trees displayed unusually erect branching, and poor flowering and yield. Since 2002, the commercial significance of AVG, its cause, and strategies to alleviate its affects, has been studied. The cause is still unknown, and AVG remains a serious threat to orchard viability. AVG affects both commercial and urban macadamia. It occurs predominantly in the warmer-drier production regions of Queensland and New South Wales. An estimated 100,000 orchard trees are affected, equating to an annual loss of $ 10.5 M. In orchards, AVG occurs as aggregations of affected trees, affected tree number can increase by 4.5% per year, and yield reduction can exceed 30%. The more upright cultivars 'HAES 344' and '741' are highly susceptible, while the more spreading cultivars 'A4', 'A16' and 'A268' show tolerance. Incidence is higher (p<0.05) in soils of high permeability and good drainage. No soil chemical anomaly has been found. Fine root dry weight of AVG trees (0-15 cm depth) was found lower (p<0.05) than non-AVG. Next generation sequencing has led to the discovery of a new Bacillus sp. and a bipartite Geminivirus, which may have a role in the disease. Trunk cinctures will increase (p<0.05) yield of moderately affected trees. Further research is needed to clarify whether a pathogen is the cause, the role of soil moisture in AVG, and develop a varietal solution.

Characterization, development and multiplexing of microsatellite markers in three commercially exploited reef fish and their application for stock identification

Thirty-four microsatellite loci were isolated from three reef fish species; golden snapper Lutjanus johnii, blackspotted croaker Protonibea diacanthus and grass emperor Lethrinus laticaudis using a next generation sequencing approach. Both IonTorrent single reads and Illumina MiSeq paired-end reads were used, with the latter demonstrating a higher quality of reads than the IonTorrent. From the 1–1.5 million raw reads per species, we successfully obtained 10–13 polymorphic loci for each species, which satisfied stringent design criteria. We developed multiplex panels for the amplification of the golden snapper and the blackspotted croaker loci, as well as post-amplification pooling panels for the grass emperor loci. The microsatellites characterized in this work were tested across three locations of northern Australia. The microsatellites we developed can detect population differentiation across northern Australia and may be used for genetic structure studies and stock identification.