Diagnosis of vibriosis in the era of genomics: lessons from invertebrates

Global changes linked to increases in temperature and ocean acidification, but also to more direct anthropogenic influences such as aquaculture, have caused a worldwide increase in the reports of Vibrio-associated illnesses affecting humans and also animals such as shrimp and molluscs. Investigation of the emergence of Vibrio pathogenesis events requires the analysis of microbial evolution at the gene, genome and population levels, in order to identify genomic modifications linked to increased virulence, resistance and/or prevalence, or to recent host shift. From a more applied point of view, the elucidation of virulence mechanisms is a prerequisite to devising prophylactic methods to fight infectious agents. In comparison with human pathogens, fairly little is known about the requirements for virulence in vibrios pathogenic to animals. However, the advent of genome sequencing, especially next-generation technologies,the possibility of genetically manipulating most of the Vibrio strains, and the recent availability of standardised animals for experimental infections have now compensated for the considerable delay in advancement of the knowledge of non-model pathogens such as Vibrio and have led to new scientific questions.

Proteomic analysis of Rhizobium freirei PRF 81 responses to low pH.

Rhizobium freirei PRF 81 is employed in common bean commercial inoculants in Brazil, due to its outstanding efficiency in fixing nitrogen, competitiveness and tolerance to abiotic stresses. Among the environmental conditions faced by rhizobia in soils, acidity is perhaps the encountered most, especially in Brazil. So, we used proteomics based approaches to study the responses of PRF 81 to a low pH condition. R. freirei PRF 81 was grown in TY medium until exponential phase in two treatments: pH 6,8 and pH 4,8. Whole-cell proteins were extracted and separated by two-dimensional gel electrophoresis, using IPG-strips with pH range 4-7 and 12% polyacrilamide gels. The experiment was performed in triplicate. Protein spots were detected in the high-resolution digitized gel images and analyzed by Image Master 2D Platinum v 5.0 software. Relative volumes (%vol) of compared between the two conditions tested and were statistically evaluated (p ≤ 0.05). Even knowing that R. freirei PRF 81 can still grow in more acid conditions, pH 4.8 was chosen because didn´t affect significantly the bacterial growth kinetics, a factor that could compromise the analysis. Using a narrow pH range, the gel profiles displayed a better resolution and reprodutibility than using broader pH range. Spots were mostly concentrated between pH 5-7 and molecular masses between 17-95 kDa. From the six hundred well-defined spots analyzed, one hundred and sixty-three spots presented a significant change in % vol, indicating that the pH led to expressive changes in the proteome of R. freirei PRF 81. Of these, sixty-one were up-regulated and one hundred two was downregulated in pH 4.8 condition. Also, fourteen spots were only identified in the acid condition, while seven spots was exclusively detected in pH 6.8. Ninety-five differentially expressed spots and two exclusively detected in pH 4,8 were selected for Maldi-Tof identification. Together with the genome sequencing and the proteome analysis of heat stress, we will search for molecular determinants of PRF 81 related to capacity to adapt to stressful tropical conditions.

Coexistence of mcr-1 and blaNDM-1 in Escherichia coli from Venezuela

We studied the presence of the mobile colistin resistance gene mcr-1 in human, animal, and environmental Enterobacteriaceae samples from Cumana, Venezuela, that were collected in 2015. The mcr-1 gene was detected in 2/93 Escherichia coli isolates from swine (novel ST452) and human (ST19) samples that were resistant to colistin. Whole-genome sequencing and transformation experiments identified mcr-1 on an IncI2 plasmid. One of the isolates also bore the widely spread carbapenemase NDM-1. A One Health approach is necessary to further elucidate the flux of these high-risk genes.

From field to fermentation: characterisation and application of non-dairy cultures in dairy foods

This thesis investigates the phenotypic and genotypic diversity of non-dairy L. lactis strains and their application to dairy fermentations. A bank of non-dairy lactococci were isolated from grass, vegetables and the bovine rumen. Subsequent analysis of these L. lactis strains revealed seven strains to possess cremoris genotypes which did not correlate with their observed phenotypes. Multi-locus sequence typing (MLST) and average nucleotide identity (ANI) highlighted the genetic diversity of lactis and cremoris subspecies. The application of these non-dairy lactococci to cheese production was also assessed. In milk, non-dairy strains formed diverse volatile profiles and selected strains were used as adjuncts in a mini Gouda-type cheese system. Sensory analysis showed non-dairy strains to be strongly associated with the development of off-flavours and bitterness. However, microfluidisation appeared to reduce bitterness. A novel bacteriophage, ɸL47, was isolated using the grass isolate L. lactis ssp. cremoris DPC6860 as a host. The phage, a member of the Siphoviridae, possessed a long tail fiber, previously unseen in dairy lactococcal phages. Genome sequencing revealed ɸL47 to be the largest sequenced lactococcal phage to date and owing to the high % similarity with ɸ949, a second member of the 949 group. Finally, to identify and characterise specific genes which may be important in niche adaptation and for applications to dairy fermentations, comparative genome sequence analysis was performed on L. lactis from corn (DPC6853), the bovine rumen (DPC6853) and grass (DPC6860). This study highlights the contribution of niche specialisation to the intra-species diversity of L. lactis and the adaptation of this organism to different environments. In summary this thesis describes the genetic diversity of L. lactis strains from outside the dairy environment and their potential application in dairy fermentations.

Evidences for an opportunistic and endophytic lifestyle of the Bursaphelenchus xylophilus -associated bacteria Serratia marcescens PWN146 isolated from wilting Pinus pinaster

Pine wilt disease (PWD) results from the interaction of three elements: the pathogenic nematode, Bursaphelenchus xylophilus; the insect-vector, Monochamus sp.; and the host tree, mostly Pinus species. Bacteria isolated from B. xylophilus may be a fourth element in this complex disease. However, the precise role of bacteria in this interaction is unclear as both plant-beneficial and as plant-pathogenic bacteria may be associated with PWD. Using whole genome sequencing and phenotypic characterization, we were able to investigate in more detail the genetic repertoire of Serratia marcescens PWN146, a bacterium associated with B. xylophilus. We show clear evidence that S. marcescens PWN146 is able to withstand and colonize the plant environment, without having any deleterious effects towards a susceptible host (Pinus thunbergii), B. xylophilus nor to the nematode model C. elegans. This bacterium is able to tolerate growth in presence of xenobiotic/organic compounds, and use phenylacetic acid as carbon source. Furthermore, we present a detailed list of S. marcescens PWN146 potentials to interfere with plant metabolism via hormonal pathways and/or nutritional acquisition, and to be competitive against other bacteria and/or fungi in terms of resource acquisition or production of antimicrobial compounds. Further investigation is required to understand the role of bacteria in PWD. We have now reinforced the theory that B. xylophilus-associated bacteria may have a plant origin.

Phenotypic and genotypic characterisation of antimicrobial resistance in Escherichia coli indicator of animal, food and human origin.

This thesis presents AMR phenotypic evaluation and whole genome sequencing analysis of 288 Escherichia coli strains isolated from different sources (livestock, companion animal, wildlife, food and human) in Italy. Our data reflects general resistance trends in Europe, reporting tetracycline, ampicillin, sulfisoxazole and aminoglycosides resistance as the most common phenotypic AMR profile among livestock, pets, wildlife and humans. Identification of human and animal (livestock and companion animal) AMR profiles in niches with a rare (fishery, mollusc) or absent (vegetable, wild animal, wild boar) direct exposure to antimicrobials, suggests widespread environmental pollution with ARGs conferring resistance to these antimicrobials. Phenotypic resistance to highest priority critically important antimicrobials was mainly observed in food-producing animals and related food such as rabbit, poultry, beef and swine. Discrepancies between AMR phenotypic pattern and genetic profile were observed. In particular, phenotypic aminoglycoside, cephalosporin, meropenem, colistin resistance and ESBL profile did not have a genetic explanation in different cases. This data could suggest the diffusion of new genetic variants of ARGs, associated to these antimicrobial classes. Generally, our collection shows a virulence profile typical of extraintestinal pathogenic Escherichia coli (ExPEC) pathotype. Different pandemic and emerging ExPEC lineages were identified, in particular in poultry meat (ST10; ST23; ST69, ST117; ST131). Rabbit was suggested as a source of ST20-ST40 potential hybrid pathogens. Wildlife carried a high average number (10) of VAGs (mostly associated to ExPEC pathotype) and different predominant ExPEC lineages (ST23, ST117, ST648), suggesting its possible involvement in maintenance and diffusion of virulence determinants. In conclusion, our study provides important knowledge related to the phenotypic/genetic AMR and virulence profiles circulating in E. coli in Italy. The role of different niches in AMR dynamics has been discussed. In particular, food-producing animals are worthy of continued investigation as a source of potential zoonotic pathogens, meanwhile wildlife might contribute to VAGs spread.

Contribution of rare damaging variants in Autism Spectrum Disorder

Autism Spectrum Disorder (ASD) is a heterogeneous and highly heritable neurodevelopmental disorder with a complex genetic architecture, consisting of a combination of common low-risk and more penetrant rare variants. This PhD project aimed to explore the contribution of rare variants in ASD susceptibility through NGS approaches in a cohort of 106 ASD families including 125 ASD individuals. Firstly, I explored the contribution of inherited rare variants towards the ASD phenotype in a girl with a maternally inherited pathogenic NRXN1 deletion. Whole exome sequencing of the trio family identified an increased burden of deleterious variants in the proband that could modulate the CNV penetrance and determine the disease development. In the second part of the project, I investigated the role of rare variants emerging from whole genome sequencing in ASD aetiology. To properly manage and analyse sequencing data, a robust and efficient variant filtering and prioritization pipeline was developed, and by its application a stringent set of rare recessive-acting and ultra-rare variants was obtained. As a first follow-up, I performed a preliminary analysis on de novo variants, identifying the most likely deleterious variants and highlighting candidate genes for further analyses. In the third part of the project, considering the well-established involvement of calcium signalling in the molecular bases of ASD, I investigated the role of rare variants in voltage-gated calcium channels genes, that mainly regulate intracellular calcium concentration, and whose alterations have been correlated with enhanced ASD risk. Specifically, I functionally tested the effect of rare damaging variants identified in CACNA1H, showing that CACNA1H variation may be involved in ASD development by additively combining with other high risk variants. This project highlights the challenges in the analysis and interpretation of variants from NGS analysis in ASD, and underlines the importance of a comprehensive assessment of the genomic landscape of ASD individuals.

Liquid biopsy for prostate cancer molecular characterization: homologous recombination system and genomic profile

Pathogenic aberrations in homologous recombination DNA repair (HRR) genes occur in approximately 1 to 4 men with advanced prostate cancer (PCa). Treatment with PARP inhibitors (PARPi) has recently been introduced for metastatic castration-resistant PCa patients, increasing clinicians' interest in the molecular characterization of all PCa patients. The limitations of using old, low-quality tumor tissue for genetic analysis, which is very common for PCa, can be overcome by using liquid biopsy as an alternative biomarker source. In this study, we aimed to evaluate the detection of molecular alterations in HRR genes on liquid biopsy compared with tumor tissue from PCa patients. Secondarily, we explored the genomic instability score (GIS), and a broader range of gene alterations for in-depth characterization of the PCa cohort. Plasma samples were collected from 63 patients with PCa. Sophia Homologous Recombination Solution (targeting 16 HRR genes) and shallow whole genome sequencing (sWGS) were used for genomic analysis of tissue DNA and circulating tumor DNA (ct). A total of 33 alterations (mainly on TP53, ATM, CHEK2, CDK12, and BRCA1/2) were identified in 28,5% of PCa plasma patients. By integrating the mutational and sWGS data, the HRR status of PCa patients was determined and a concordance agreement of 85,7% was identified with tumor tissue. A median GIS of 15 was obtained, reaching a score of 63 in 2 samples with double alterations, BRCA1 and TP53. We explored the PCa mutation landscape, and the most significant enriched pathways identified were the sphingosine 1-phosphate (S1P) receptor signaling and the PI3K-AKT-mTOR pathway. HRR analysis on FFPE and liquid biopsy samples show high concordance, demonstrating that the noninvasive ctDNA-enriched plasma can be an optimal alternative source for molecular SNV and CNV analysis. In addition, the evaluation of GIS and pathway interaction should be considered for more comprehensive molecular characterization in PCa patients.

Genome-wide re-sequencing of multidrug-resistant Mycobacterium leprae Airaku-3.

Genotyping and molecular characterization of drug resistance mechanisms in Mycobacterium leprae enables disease transmission and drug resistance trends to be monitored. In the present study, we performed genome-wide analysis of Airaku-3, a multidrug-resistant strain with an unknown mechanism of resistance to rifampicin. We identified 12 unique non-synonymous single-nucleotide polymorphisms (SNPs) including two in the transporter-encoding ctpC and ctpI genes. In addition, two SNPs were found that improve the resolution of SNP-based genotyping, particularly for Venezuelan and South East Asian strains of M. leprae.

Sequencing and characterizing the genome of Estrella lausannensis as an undergraduate project: training students and biological insights.

With the widespread availability of high-throughput sequencing technologies, sequencing projects have become pervasive in the molecular life sciences. The huge bulk of data generated daily must be analyzed further by biologists with skills in bioinformatics and by "embedded bioinformaticians," i.e., bioinformaticians integrated in wet lab research groups. Thus, students interested in molecular life sciences must be trained in the main steps of genomics: sequencing, assembly, annotation and analysis. To reach that goal, a practical course has been set up for master students at the University of Lausanne: the "Sequence a genome" class. At the beginning of the academic year, a few bacterial species whose genome is unknown are provided to the students, who sequence and assemble the genome(s) and perform manual annotation. Here, we report the progress of the first class from September 2010 to June 2011 and the results obtained by seven master students who specifically assembled and annotated the genome of Estrella lausannensis, an obligate intracellular bacterium related to Chlamydia. The draft genome of Estrella is composed of 29 scaffolds encompassing 2,819,825 bp that encode for 2233 putative proteins. Estrella also possesses a 9136 bp plasmid that encodes for 14 genes, among which we found an integrase and a toxin/antitoxin module. Like all other members of the Chlamydiales order, Estrella possesses a highly conserved type III secretion system, considered as a key virulence factor. The annotation of the Estrella genome also allowed the characterization of the metabolic abilities of this strictly intracellular bacterium. Altogether, the students provided the scientific community with the Estrella genome sequence and a preliminary understanding of the biology of this recently-discovered bacterial genus, while learning to use cutting-edge technologies for sequencing and to perform bioinformatics analyses.

Development and validation of a Xanthomonas axonopodis pv. citri DNA microarray platform (XACarray) generated from the shotgun libraries previously used in the sequencing of this bacterial genome

Background. From shotgun libraries used for the genomic sequencing of the phytopathogenic bacterium Xanthomonas axonopodis pv. citri (XAC), clones that were representative of the largest possible number of coding sequences (CDSs) were selected to create a DNA microarray platform on glass slides (XACarray). The creation of the XACarray allowed for the establishment of a tool that is capable of providing data for the analysis of global genome expression in this organism. Findings. The inserts from the selected clones were amplified by PCR with the universal oligonucleotide primers M13R and M13F. The obtained products were purified and fixed in duplicate on glass slides specific for use in DNA microarrays. The number of spots on the microarray totaled 6,144 and included 768 positive controls and 624 negative controls per slide. Validation of the platform was performed through hybridization of total DNA probes from XAC labeled with different fluorophores, Cy3 and Cy5. In this validation assay, 86% of all PCR products fixed on the glass slides were confirmed to present a hybridization signal greater than twice the standard deviation of the deviation of the global median signal-to-noise ration. Conclusions. Our validation of the XACArray platform using DNA-DNA hybridization revealed that it can be used to evaluate the expression of 2,365 individual CDSs from all major functional categories, which corresponds to 52.7% of the annotated CDSs of the XAC genome. As a proof of concept, we used this platform in a previously work to verify the absence of genomic regions that could not be detected by sequencing in related strains of Xanthomonas. © 2010 Moreira et al; licensee BioMed Central Ltd.

The genome sequence of the plant pathogen Xylella fastidiosa: The Xylella fastidiosa consortium of the organization for nucleotide sequencing and analysis, Sao Paulo, Brazil

Xylella fastidiosa is a fastidious, xylem-limited bacterium that causes a range of economically important plant diseases. Here we report the complete genome sequence of X. fastidiosa clone 9a5c, which causes citrus variegated chlorosis - a serious disease of orange trees. The genome comprises a 52.7% GC-rich 2,679,305-base-pair (bp) circular chromosome and 'two plasmids of 51,158 bp and 1,285 bp. We can assign putative functions to47% of the 2,904 predicted coding regions. Efficient metabolic functions are predicted, with sugars as the principal energy and carbon source, supporting existence in the nutrient-poor xylem sap. The mechanisms associated with pathogenicity and virulence involve toxins, antibiotics and ion sequestration systems, as well as bacterium-bacterium and bacterium-host interactions mediated by a range of proteins. Orthologues of some of these proteins have only been identified in animal and human pathogens; their presence in X. fastidiosa indicates that the molecular basis for bacterial pathogenicity is both conserved and independent of host. At least 83 genes are bacteriophage-derived and include virulence-associated genes from other bacteria, providing direct evidence of phage-mediated horizontal gene transfer.

Development and validation of a Xanthomonas axonopodis pv. citri DNA microarray platform (XACarray) generated from the shotgun libraries previously used in the sequencing of this bacterial genome

Abstract Background From shotgun libraries used for the genomic sequencing of the phytopathogenic bacterium Xanthomonas axonopodis pv. citri (XAC), clones that were representative of the largest possible number of coding sequences (CDSs) were selected to create a DNA microarray platform on glass slides (XACarray). The creation of the XACarray allowed for the establishment of a tool that is capable of providing data for the analysis of global genome expression in this organism. Findings The inserts from the selected clones were amplified by PCR with the universal oligonucleotide primers M13R and M13F. The obtained products were purified and fixed in duplicate on glass slides specific for use in DNA microarrays. The number of spots on the microarray totaled 6,144 and included 768 positive controls and 624 negative controls per slide. Validation of the platform was performed through hybridization of total DNA probes from XAC labeled with different fluorophores, Cy3 and Cy5. In this validation assay, 86% of all PCR products fixed on the glass slides were confirmed to present a hybridization signal greater than twice the standard deviation of the deviation of the global median signal-to-noise ration. Conclusions Our validation of the XACArray platform using DNA-DNA hybridization revealed that it can be used to evaluate the expression of 2,365 individual CDSs from all major functional categories, which corresponds to 52.7% of the annotated CDSs of the XAC genome. As a proof of concept, we used this platform in a previously work to verify the absence of genomic regions that could not be detected by sequencing in related strains of Xanthomonas.

Genome and transcriptome sequencing identifies breeding targets in the orphan crop tef (Eragrostis tef)

Background: Tef (Eragrostis tef), an indigenous cereal critical to food security in the Horn of Africa, is rich in minerals and protein, resistant to many biotic and abiotic stresses and safe for diabetics as well as sufferers of immune reactions to wheat gluten. We present the genome of tef, the first species in the grass subfamily Chloridoideae and the first allotetraploid assembled de novo. We sequenced the tef genome for marker-assisted breeding, to shed light on the molecular mechanisms conferring tef's desirable nutritional and agronomic properties, and to make its genome publicly available as a community resource. Results: The draft genome contains 672 Mbp representing 87% of the genome size estimated from flow cytometry. We also sequenced two transcriptomes, one from a normalized RNA library and another from unnormalized RNASeq data. The normalized RNA library revealed around 38000 transcripts that were then annotated by the SwissProt group. The CoGe comparative genomics platform was used to compare the tef genome to other genomes, notably sorghum. Scaffolds comprising approximately half of the genome size were ordered by syntenic alignment to sorghum producing tef pseudo-chromosomes, which were sorted into A and B genomes as well as compared to the genetic map of tef. The draft genome was used to identify novel SSR markers, investigate target genes for abiotic stress resistance studies, and understand the evolution of the prolamin family of proteins that are responsible for the immune response to gluten. Conclusions: It is highly plausible that breeding targets previously identified in other cereal crops will also be valuable breeding targets in tef. The draft genome and transcriptome will be of great use for identifying these targets for genetic improvement of this orphan crop that is vital for feeding 50 million people in the Horn of Africa.