Novel tools for conservation genomics: comparing two high-throughput approaches for SNP discovery in the transcriptome of the European hake

The growing accessibility to genomic resources using next-generation sequencing (NGS) technologies has revolutionized the application of molecular genetic tools to ecology and evolutionary studies in non-model organisms. Here we present the case study of the European hake (Merluccius merluccius), one of the most important demersal resources of European fisheries. Two sequencing platforms, the Roche 454 FLX (454) and the Illumina Genome Analyzer (GAII), were used for Single Nucleotide Polymorphisms (SNPs) discovery in the hake muscle transcriptome. De novo transcriptome assembly into unique contigs, annotation, and in silico SNP detection were carried out in parallel for 454 and GAII sequence data. High-throughput genotyping using the Illumina GoldenGate assay was performed for validating 1,536 putative SNPs. Validation results were analysed to compare the performances of 454 and GAII methods and to evaluate the role of several variables (e.g. sequencing depth, intron-exon structure, sequence quality and annotation). Despite well-known differences in sequence length and throughput, the two approaches showed similar assay conversion rates (approximately 43%) and percentages of polymorphic loci (67.5% and 63.3% for GAII and 454, respectively). Both NGS platforms therefore demonstrated to be suitable for large scale identification of SNPs in transcribed regions of non-model species, although the lack of a reference genome profoundly affects the genotyping success rate. The overall efficiency, however, can be improved using strict quality and filtering criteria for SNP selection (sequence quality, intron-exon structure, target region score).

Application of SNPs for population genetics of nonmodel organisms:new opportunities and challenges

Recent improvements in the speed, cost and accuracy of next generation sequencing are revolutionizing the discovery of single nucleotide polymorphisms (SNPs). SNPs are increasingly being used as an addition to the molecular ecology toolkit in nonmodel organisms, but their efficient use remains challenging. Here, we discuss common issues when employing SNP markers, including the high numbers of markers typically employed, the effects of ascertainment bias and the inclusion of nonneutral loci in a marker panel. We provide a critique of considerations specifically associated with the application and population genetic analysis of SNPs in nonmodel taxa, focusing specifically on some of the most commonly applied methods.

Production of extended-spectrum β-lactamases and the potential indirect pathogenic role of Prevotella isolates from the cystic fibrosis respiratory microbiota

Extended-spectrum β-lactamase (ESBL) production and the prevalence of the β-lactamase-encoding gene blaTEM were determined in Prevotella isolates (n=50) cultured from the respiratory tract of adults and young people with cystic fibrosis (CF). Time-kill studies were used to investigate the concept of passive antibiotic resistance and to ascertain whether a β-lactamase-positive Prevotella isolate can protect a recognised CF pathogen from the action of ceftazidime in vitro. The results indicated that approximately three-quarters (38/50; 76%) of Prevotella isolates produced ESBLs. Isolates positive for ESBL production had higher minimum inhibitory concentrations (MICs) of β-lactam antibiotics compared with isolates negative for production of ESBLs (P<0.001). The blaTEM gene was detected more frequently in CF Prevotella isolates from paediatric patients compared with isolates from adults (P=0.002), with sequence analysis demonstrating that 21/22 (95%) partial blaTEM genes detected were identical to blaTEM-116. Furthermore, a β-lactamase-positive Prevotella isolate protected Pseudomonas aeruginosa from the antimicrobial effects of ceftazidime (P=0.03). Prevotella isolated from the CF respiratory microbiota produce ESBLs and may influence the pathogenesis of chronic lung infection via indirect methods, including shielding recognised pathogens from the action of ceftazidime.

Rapid emergence of multidrug resistant, H58-lineage Salmonella typhi in Blantyre, Malawi

INTRODUCTION: Between 1998 and 2010, S. Typhi was an uncommon cause of bloodstream infection (BSI) in Blantyre, Malawi and it was usually susceptible to first-line antimicrobial therapy. In 2011 an increase in a multidrug resistant (MDR) strain was detected through routine bacteriological surveillance conducted at Queen Elizabeth Central Hospital (QECH).

METHODS: Longitudinal trends in culture-confirmed Typhoid admissions at QECH were described between 1998-2014. A retrospective review of patient cases notes was conducted, focusing on clinical presentation, prevalence of HIV and case-fatality. Isolates of S. Typhi were sequenced and the phylogeny of Typhoid in Blantyre was reconstructed and placed in a global context.

RESULTS: Between 1998-2010, there were a mean of 14 microbiological diagnoses of Typhoid/year at QECH, of which 6.8% were MDR. This increased to 67 in 2011 and 782 in 2014 at which time 97% were MDR. The disease predominantly affected children and young adults (median age 11 [IQR 6-21] in 2014). The prevalence of HIV in adult patients was 16.7% [8/48], similar to that of the general population (17.8%). Overall, the case fatality rate was 2.5% (3/94). Complications included anaemia, myocarditis, pneumonia and intestinal perforation. 112 isolates were sequenced and the phylogeny demonstrated the introduction and clonal expansion of the H58 lineage of S. Typhi.

CONCLUSIONS: Since 2011, there has been a rapid increase in the incidence of multidrug resistant, H58-lineage Typhoid in Blantyre. This is one of a number of reports of the re-emergence of Typhoid in Southern and Eastern Africa. There is an urgent need to understand the reservoirs and transmission of disease and how to arrest this regional increase.

Characterisation of a secretory serine protease inhibitor (SjB6) from Schistosoma japonicum

BACKGROUND: Proteins belonging to the serine protease inhibitor (serpin) superfamily play essential physiological roles in many organisms. In pathogens, serpins are thought to have evolved specifically to limit host immune responses by interfering with the host immune-stimulatory signals. Serpins are less well characterised in parasitic helminths, although some are thought to be involved in mechanisms associated with host immune modulation. In this study, we cloned and partially characterised a secretory serpin from Schistosoma japonicum termed SjB6, these findings provide the basis for possible functional roles.

METHODS: SjB6 gene was identified through database mining of our previously published microarray data, cloned and detailed sequence and structural analysis and comparative modelling carried out using various bioinformatics and proteomics tools. Gene transcriptional profiling was determined by real-time PCR and the expression of native protein determined by immunoblotting. An immunological profile of the recombinant protein produced in insect cells was determined by ELISA.

RESULTS: SjB6 contains an open reading frame of 1160 base pairs that encodes a protein of 387 amino acid residues. Detailed sequence analysis, comparative modelling and structural-based alignment revealed that SjB6 contains the essential structural motifs and consensus secondary structures typical of inhibitory serpins. The presence of an N-terminal signal sequence indicated that SjB6 is a secretory protein. Real-time data indicated that SjB6 is expressed exclusively in the intra-mammalian stage of the parasite life cycle with its highest expression levels in the egg stage (p < 0.0001). The native protein is approximately 60 kDa in size and recombinant SjB6 (rSjB6) was recognised strongly by sera from rats experimentally infected with S. japonicum.

CONCLUSIONS: The significantly high expression of SjB6 in schistosome eggs, when compared to other life cycle stages, suggests a possible association with disease pathology, while the strong reactivity of sera from experimentally infected rats against rSjB6 suggests that native SjB6 is released into host tissue and induces an immune response. This study presents a comprehensive demonstration of sequence and structural-based analysis of a secretory serpin from a trematode and suggests SjB6 may be associated with important functional roles in S. japonicum, particularly in parasite modulation of the host microenvironment.

Functional characterization of SjB10, an intracellular serpin from Schistosoma japonicum

Serine protease inhibitors (serpin) play essential roles in many organisms. Mammalian serpins regulate the blood coagulation, fibrinolysis, inflammation and complement activation pathways. In parasitic helminths, serpins are less well characterized, but may also be involved in evasion of the host immune response. In this study, a Schistosoma japonicum serpin (SjB10), containing a 1212 bp open reading frame (ORF), was cloned, expressed and functionally characterized. Sequence analysis, comparative modelling and structural-based alignment revealed that SjB10 contains the essential structural motifs and consensus secondary structures of inhibitory serpins. Transcriptional profiling demonstrated that SjB10 is expressed in adult males, schistosomula and eggs but particularly in the cercariae, suggesting a possible role in cercarial penetration of mammalian host skin. Recombinant SjB10 (rSjB10) inhibited pancreatic elastase (PE) in a dose-dependent manner. rSjB10 was recognized strongly by experimentally infected rat sera indicating that native SjB10 is released into host tissue and induces an immune response. By immunochemistry, SjB10 localized in the S. japonicum adult foregut and extra-embryonic layer of the egg. This study provides a comprehensive demonstration of sequence and structural-based analysis of a functional S. japonicum serpin. Furthermore, our findings suggest that SjB10 may be associated with important functional roles in S. japonicum particularly in host-parasite interactions.

Transcriptional profiling of the oesophageal gland region of male worms of Schistosoma mansoni

The intestinal tract of schistosomes opens at the mouth and leads into the foregut or oesophageal region that is lined with syncytium continuous with the apical cytoplasm of the tegument. The oesophagus is surrounded by a specialised gland, the oesophageal gland. This gland releases materials into the lumen of the oesophagus and the region is thought to initiate the lysis of erythrocytes and neutralisation of immune effectors of the host. The oesophageal region is present in the early invasive schistosomulum, a stage potentially targetable by anti-schistosome vaccines. We used a 44k oligonucleotide microarray to identify highly up-regulated genes in microdissected frozen sections of the oesophageal gland of male worms of S. mansoni. We show that 122 genes were up-regulated 2-fold or higher in the oesophageal gland compared with a whole male worm tissue control. The enriched genes included several associated with lipid metabolism and transmembrane transport as well as some micro-exon genes. Since the oesophageal gland is important in the initiation of digestion and the fact that it develops early after invasion of the mammalian host, further study of selected highly up-regulated functionally important genes in this tissue may reveal new anti-schistosome intervention targets for schistosomiasis control.

Unexpected Activity of a Novel Kunitz-Type Inhibtior: Inhibition of Cysteine Proteases but not Serine Proteases

Kunitz-type (KT) protease inhibitors are low molecular weight proteins classically defined as serine protease inhibitors. We identified a novel secreted KT inhibitor associated with the gut and parenchymal tissues of the infective juvenile stage of Fasciola hepatica, a helminth parasite of medical and veterinary importance. Unexpectedly, recombinant KT inhibitor (rFhKT1) exhibited no inhibitory activity towards serine proteases but was a potent inhibitor of the major secreted cathepsin L cysteine proteases of F. hepatica, FhCL1 and FhCL2, and of human cathepsins L and K (Ki = 0.4 nM - 27 nM). FhKT1 prevented the auto-catalytic activation of FhCL1 and FhCL2 and formed stable complexes with the mature enzymes. Pull-down experiments from adult parasite culture medium showed that rFhKT1 interacts specifically with native secreted FhCL1, FhCL2 and FhCL5. Substitution of the unusual P1 Leu15 within the exposed reactive loop of FhKT1 for the more commonly found Arg (FhKT1Leu15/Arg15) had modest adverse effects on the cysteine protease inhibition but conferred potent activity against the serine protease trypsin (Ki = 1.5 nM). Computational docking and sequence analysis provided hypotheses for the exclusive binding of FhKT1 to cysteine proteases, the importance of the Leu15 in anchoring the inhibitor into the S2 active site pocket, and the inhibitor's selectivity towards FhCL1, FhCL2 and human cathepsins L and K. FhKT1 represents a novel evolutionary adaptation of KT protease inhibitors by F. hepatica, with its prime purpose likely in the regulation of the major parasite-secreted proteases and/or cathepsin L-like proteases of its host.

RAS mutations and cetuximab in locally advanced rectal cancer: results of the EXPERT-C trial.

Background: RAS mutations predict resistance to anti-epidermal growthfactor receptor (EGFR) monoclonal antibodies in metastatic colorectal cancer. We analysed RAS mutations in 30 non-metastatic rectal cancer patients treated with or without cetuximab within the 31 EXPERT-C trial.

Methods: Ninety of 149 patients with tumours available for analysis were KRAS/BRAF wild-type, and randomly assigned to capecitabine plus oxaliplatin (CAPOX) followed by chemoradiotherapy, surgery and adjuvant CAPOX or the same regimen plus cetuximab (CAPOX-C). Of these, four had a mutation of NRAS exon 3, and 84 were retrospectively analysed for additional KRAS (exon 4) and NRAS (exons 2/4) mutations by using bi-directional Sanger sequencing. The effect of cetuximab on study end-points in the RAS wild-type population was analysed.

Results: Eleven (13%) of 84 patients initially classified as KRAS/BRAF wild-type were found to have a mutation in KRAS exon 4 (11%) or NRAS exons 2/4 (2%). Overall, 78/149 (52%) assessable patients were RAS wild-type (CAPOX, n = 40; CAPOX-C, n = 38). In this population, after a median follow-up of 63.8 months, in line with the initial analysis, the addition of cetuximab was associated with numerically higher, but not statistically significant, rates of complete response (15.8% versus 7.5%, p = 0.31), 5-year progression-free survival (75.5% versus 67.5%, hazard ratio (HR) 0.61, p = 0.25) and 5-year overall survival (83.8% versus 70%, HR 0.54, p = 0.20).

Conclusions: RAS mutations beyond KRAS exon 2 and 3 were identified in 17% of locally advanced rectal cancer patients. Given the small sample size, no definitive conclusions on the effect of additional RAS mutations on cetuximab treatment in this setting can be drawn and further investigation of RAS in larger studies is warranted.

Frequency of mutations and polymorphisms in borderline ovarian tumors of known cancer genes.

Borderline ovarian tumors represent an understudied subset of ovarian tumors. Most studies investigating aberrations in borderline tumors have focused on KRAS/BRAF mutations. In this study, we conducted an extensive analysis of mutations and single-nucleotide polymorphisms (SNPs) in borderline ovarian tumors. Using the Sequenom MassArray platform, we investigated 160 mutations/polymorphisms in 33 genes involved in cell signaling, apoptosis, angiogenesis, cell cycle regulation and cellular senescence. Of 52 tumors analyzed, 33 were serous, 18 mucinous and 1 endometrioid. KRAS c.35G>A p.Gly12Asp mutations were detected in eight tumors (six serous and two mucinous), BRAF V600E mutations in two serous tumors, and PIK3CA H1047Y and PIK3CA E542K mutations in a serous and an endometrioid BOT, respectively. CTNNB1 mutation was detected in a serous tumor. Potentially functional polymorphisms were found in vascular endothelial growth factor (VEGF), ABCB1, FGFR2 and PHLPP2. VEGF polymorphisms were the most common and detected at four loci. PHLPP2 polymorphisms were more frequent in mucinous as compared with serous tumors (P=0.04), with allelic imbalance in one case. This study represents the largest and most comprehensive analysis of mutations and functional SNPs in borderline ovarian tumors to date. At least 25% of borderline ovarian tumors harbor somatic mutations associated with potential response to targeted therapeutics.