Molecular typing of Campylobacter jejuni and C. coli from chickens and patients with gastritis or Guillain‐Barré syndrome based on multilocus sequence types and pulsed‐field gel electrophoresis patterns

Campylobacter jejuni has recently been noted as the most common cause of bacterial foodborne diseases in Japan. In the present study, we determined ST types of C. jejuni and Campylobacter coli isolated from chickens and patients with enteritis or GBS in Japan and Thailand. C. jejuni from chickens, enteritis, and GBS exhibited divergent ST types and included several novel types in addition to worldwide common types. C. coli from enteritis was also divergent. Novel ST types may represent unidentified native clones in each country. Pulsed‐field gel electrophoresis confirmed the above typing and demonstrated long‐term persistence and transmission.     

DNA identification and characterization of Campylobacter jejuni and Campylobacter coli isolated from caecal samples of chickens in Grenada

Aims: To speciate Campylobacter strains from the caeca of chickens in Grenada using PCR and to evaluate DNA‐based typing methods for the characterization of these isolates. Methods and Results: Isolates were speciated with two multiplex PCR assays and were typed with flaA‐RFLP, pulsed‐field gel electrophoresis (PFGE) and multilocus sequence typing (MLST). Results confirmed that Campylobacter coli strains were more predominant than Campylobacter jejuni strains. From 56 isolates, 18 were misidentified using biochemical tests. PFGE typing gave the highest discriminatory power among the methods used (Simpson’s index of diversity, D = 0·9061). However, the combination of flaA‐RFLP, PFGE and MLST results gave the highest discrimination for subtyping of these isolates (D = 0·9857). A band position tolerance of 4% in Bio Numerics was the most appropriate for the analysis of this database. MLST profiles were generally concordant with PFGE and/or flaA‐RFLP types. Several isolates exhibited new MLST sequence types (STs), and 43 of the 49 Camp. coli strains belonged to the ST‐828 clonal complex. Conclusions: Campylobacter coli was the most prevalent species isolated from broilers and layers in Grenada, and a combination of restriction and sequence methods was most appropriate for the typing of Camp. coli isolates. Campylobacter coli STs clustered with described poultry‐associated Camp. coli STs by phylogenetic analysis. Significance and Impact of the Study: Further studies to understand the predominance of Camp. coli within Campylobacter spp. from chickens in Grenada may help elucidate the epidemiology of these pathogens in chickens.     

Multi‐locus sequence types of Campylobacter carried by flies and slugs acquired from local ruminant faeces

Aims: To assess whether flies and slugs acquire strains of Campylobacter jejuni and Campylobacter coli present in local ruminant faeces. Methods and Results: Campylobacter was cultured from flies, slugs and ruminant faeces that were collected from a single farm in Scotland over a 19‐week period. The isolates were typed using multi‐locus sequence typing (MLST) and compared with isolates from cattle and sheep faeces. Campylobacter jejuni and Camp. coli were isolated from 5·8% (n = 155, average of 75 flies per pool) and 13·3% (n = 15, average of 8·5 slugs per pool) of pooled fly and slug samples, respectively. The most common sequence type (ST) in flies was Camp. coli ST‐962 (approx. 40%) regardless of the prevalence in local cattle (2·3%) or sheep (25·0%) faeces. Two positive slug pools generated the same ST that has not been reported elsewhere. Conclusions: Despite their low carriage rate, flies are able to acquire Campylobacter STs that are locally present, although the subset carried may be biased when compared to local source. Slugs were shown to carry a previously unreported Campylobacter ST. Significance and Impact of the Study: This study has demonstrated that flies carry viable Campylobacter and may contribute to the transfer of STs within and between groups of animals on farms. Further, they may therefore present a risk to human health via their contact with ready‐to‐eat foods or surfaces.     

Genome‐wide association of functional traits linked with Campylobacter jejuni survival from farm to fork

Campylobacter jejuni is a major cause of bacterial gastroenteritis worldwide, primarily associated with the consumption of contaminated poultry. C. jejuni lineages vary in host range and prevalence in human infection, suggesting differences in survival throughout the poultry processing chain. From 7343 MLST‐characterised isolates, we sequenced 600 C. jejuni and C. coli isolates from various stages of poultry processing and clinical cases. A genome‐wide association study (GWAS) in C. jejuni ST‐21 and ST‐45 complexes identified genetic elements over‐represented in clinical isolates that increased in frequency throughout the poultry processing chain. Disease‐associated SNPs were distinct in these complexes, sometimes organised in haplotype blocks. The function of genes containing associated elements was investigated, demonstrating roles for cj1377c in formate metabolism, nuoK in aerobic survival and oxidative respiration, and cj1368‐70 in nucleotide salvage. This work demonstrates the utility of GWAS for investigating transmission in natural zoonotic pathogen populations and provides evidence that major C. jejuni lineages have distinct genotypes associated with survival, within the host specific niche, from farm to fork.     

Comparison between the biofilm initiation of Campylobacter jejuni and Campylobacter coli strains to an inert surface using BioFilm Ring Test®

Aims: The adhesion to an inert surface (the first step of biofilm formation) of the two main pathogenic Campylobacter species, Campylobacter jejuni and Campylobacter coli, isolated from diverse origins, was compared. Methods and Results: Adhesion assays were conducted in 96‐well, polystyrene microtiter plates using the BioFilm Ring Test^® method. This new technique, based on magnetic bead entrapment, was shown to be suitable for analysing the adhesion of Campylobacter sp. strains by comparing the adhesion of four C. jejuni strains as revealed by the BioFilm Ring Test^® and immunodetection. Among the 46 strains tested, C. jejuni and C. coli displayed different adhesion capabilities ranging from no adhesion to strong adhesion. However, no strain of C. coli was strongly adherent, and statistically, C. coli adhered less to an inert surface than C. jejuni. In addition, strains isolated from animals or carcasses were less adherent than those isolated from food‐processing and clinical cases. Conclusions: These observations suggest that the food environment and the human body could have selected strains with greater adhesion. Significance and Impact of the Study: The adhesion capability of strains could partly explain the cross‐contamination or re‐contamination of food products by Campylobacter. This property could provide a mode of survival for Campylobacter in the food chain.     

Rapid host switching in generalist Campylobacter strains erodes the signal for tracing human infections

Campylobacter jejuni and Campylobacter coli are the biggest causes of bacterial gastroenteritis in the developed world, with human infections typically arising from zoonotic transmission associated with infected meat. Because Campylobacter is not thought to survive well outside the gut, host-associated populations are genetically isolated to varying degrees. Therefore, the likely origin of most strains can be determined by host-associated variation in the genome. This is instructive for characterizing the source of human infection. However, some common strains, notably isolates belonging to the ST-21, ST-45 and ST-828 clonal complexes, appear to have broad host ranges, hindering source attribution. Here whole-genome sequencing has the potential to reveal fine-scale genetic structure associated with host specificity. We found that rates of zoonotic transmission among animal host species in these clonal complexes were so high that the signal of host association is all but obliterated, estimating one zoonotic transmission event every 1.6, 1.8 and 12 years in the ST-21, ST-45 and ST828 complexes, respectively. We attributed 89% of clinical cases to a chicken source, 10% to cattle and 1% to pig. Our results reveal that common strains of C. jejuni and C. coli infectious to humans are adapted to a generalist lifestyle, permitting rapid transmission between different hosts. Furthermore, they show that the weak signal of host association within these complexes presents a challenge for pinpointing the source of clinical infections, underlining the view that whole-genome sequencing, powerful though it is, cannot substitute for intensive sampling of suspected transmission reservoirs.     

Ribosomal operon intergenic sequence region (ISR) heterogeneity in Campylobacter coli and Campylobacter jejuni

Aims: The intergenic sequence regions (ISR) between the 16S and 23S genes of Campylobacter jejuni and Campylobacter coli are markedly different for each species. However, in the genomic sequence for Camp. coli RM2228 , two rRNA operons have an ISR that is characteristic of Camp. coli, and the third operon is characteristic of Camp. jejuni. The aim of this study was to determine the prevalence of ISR heterogeneity in these organisms. Methods and Results: PCR primers were designed to yield a 327‐base pair (bp) product for Camp. coli and 166‐bp product for Camp. jejuni. A strain like Camp. coli RM2228 should yield products of both sizes. DNA from a panel of Camp. coli (n = 133) and Camp. jejuni (n = 134) isolates were tested. All of the isolates yielded products of the predicted size for the species. To verify the data for Camp. coli RM2228 , each ribosomal operon from the isolate was individually amplified by PCR and tested with the ISR primer pair. Products of both sizes were produced as predicted. Conclusions: The cross‐species heterogeneity of the ISR seen in Camp. coli RM2228 is uncommon. Significance and Impact of the Study: The heterogeneity must have been caused by horizontal gene transfer at a frequency lower than predicted from housekeeping gene data. Thus, it can be expected that species identification based on the ISR can be confused in rare isolates.     

Demonstration of persistent strains of Campylobacter jejuni within broiler farms over a 1‐year period in Lithuania

Aims: The aim of the study was to investigate the flock prevalence of Campylobacter jejuni and Campylobacter coli in broiler farms in Lithuania and to identify possible persistent strains of Camp. jejuni using amplified fragment length polymorphism (AFLP) typing method. Methods and Results: During 1 year, 42 broiler flocks from 9 broiler farms were examined to determine the prevalence of Campylobacter‐positive broiler flocks in Lithuania. Among 42 broiler flocks examined, 31 flocks (73·8%) were positive for Camp. jejuni and 17 flocks (40·48%) for Camp. coli. Campylobacter jejuni isolates were genotyped by AFLP method using BspDI and BglII restriction enzymes. Typing of 190 isolates generated 50 AFLP genotypes with the highest diversity of strains found in the summer season. Each farm showed one or more predominant AFLP types, and one AFLP type (A32) was found in five broiler farms over a 1‐year period. Conclusions: Campylobacter jejuni and Camp. coli are highly prevalent in broiler farms in Lithuania. Farm‐specific genotypes were identified in all farms examined. Type A32 was present and persisted in different broiler farms, and a common source of transmission of Camp. jejuni was suspected. Significance and Impact of the Study: For the first time, Camp. jejuni in broiler flocks has been genetically characterized in Lithuania. Persistent strains of Camp. jejuni were detected over one period at the beginning of broiler meat production chain and, therefore, the identification of contamination source of such strains and the mechanism of their particular ability to persist are crucial to establish effective control measures against Camp. jejuni infection in broiler farms.     

Niche segregation and genetic structure of Campylobacter jejuni populations from wild and agricultural host species

Bacterial populations can display high levels of genetic structuring but the forces that influence this are incompletely understood. Here, by combining modelling approaches with multilocus sequence data for the zoonotic pathogen Campylobacter, we investigated how ecological factors such as niche (host) separation relate to population structure. We analysed seven housekeeping genes from published C. jejuni and C. coli isolate collections from a range of food and wild animal sources as well as abiotic environments. By reconstructing genetic structure and the patterns of ancestry, we quantified C. jejuni host association, inferred ancestral populations, investigated genetic admixture in different hosts and determined the host origin of recombinant C. jejuni alleles found in hybrid C. coli lineages. Phylogenetically distinct C. jejuni lineages were associated with phylogenetically distinct wild birds. However, in the farm environment, phylogenetically distant host animals shared several C. jejuni lineages that could not be segregated according to host origin using these analyses. Furthermore, of the introgressed C. jejuni alleles found in C. coli lineages, 73% were attributed to genotypes associated with food animals. Our results are consistent with an evolutionary scenario where distinct Campylobacter lineages are associated with different host species but the ecological factors that maintain this are different in domestic animals such that phylogenetically distant animals can harbour closely related strains.     

A chromosome‐level genome assembly of the parasitoid wasp Pteromalus puparum

Parasitoid wasps represent a large proportion of hymenopteran species. They have complex evolutionary histories and are important biocontrol agents. To advance parasitoid research, a combination of Illumina short‐read, PacBio long‐read and Hi‐C scaffolding technologies was used to develop a high‐quality chromosome‐level genome assembly for Pteromalus puparum, which is an important pupal endoparasitoid of caterpillar pests. The chromosome‐level assembly has aided in studies of venom and detoxification genes. The assembled genome size is 338 Mb with a contig N50 of 38.7 kb and a scaffold N50 of 1.16 Mb. Hi‐C analysis assembled scaffolds onto five chromosomes and raised the scaffold N50 to 65.8 Mb, with more than 96% of assembled bases located on chromosomes. Gene annotation was assisted by RNA sequencing for the two sexes and four different life stages. Analysis detected 98% of the BUSCO (Benchmarking Universal Single‐Copy Orthologs) gene set, supporting a high‐quality assembly and annotation. In total, 40.1% (135.6 Mb) of the assembly is composed of repetitive sequences, and 14,946 protein‐coding genes were identified. Although venom genes play important roles in parasitoid biology, their spatial distribution on chromosomes was poorly understood. Mapping has revealed venom gene tandem arrays for serine proteases, pancreatic lipase‐related proteins and kynurenine–oxoglutarate transaminases, which have amplified in the P. puparum lineage after divergence from its common ancestor with Nasonia vitripennis. In addition, there is a large expansion of P450 genes in P. puparum. These examples illustrate how chromosome‐level genome assembly can provide a valuable resource for molecular, evolutionary and biocontrol studies of parasitoid wasps.     

A chromosome‐level genome assembly reveals the genetic basis of cold tolerance in a notorious rice insect pest,Chilo suppressalis

The rice stem borer, Chilo suppressalis, is one of the most damaging insect pests to rice production worldwide. Although C. suppressalis has been the focus of numerous studies examining cold tolerance and diapause, plant–insect interactions, pesticide targets and resistance, and the development of RNAi‐mediated pest management, the absence of a high‐quality genome has limited deeper insights. To address this limitation, we generated a draft C. suppressalis genome constructed from both Illumina and PacBio sequences. The assembled genome size was 824.35 Mb with a contig N₅₀ of 307 kb and a scaffold N₅₀ of 1.75 Mb. Hi‐C scaffolding assigned 99.2% of the bases to one of 29 chromosomes. Based on universal single‐copy orthologues (BUSCO), the draft genome assembly was estimated to be 97% complete and is predicted to encompass 15,653 protein‐coding genes. Cold tolerance is an extreme survival strategy found in animals. However, little is known regarding the genetic basis of the winter ecology of C. suppressalis. Here, we focused our orthologous analysis on those gene families associated with animal cold tolerance. Our finding provided the first genomic evidence revealing specific cold‐tolerant strategies in C. suppressalis, including those involved in glucose‐originated glycerol biosynthesis, triacylglycerol‐originated glycerol biosynthesis, fatty acid synthesis and trehalose transport‐intermediate cold tolerance. The high‐quality C. suppressalis genome provides a valuable resource for research into a broad range of areas in molecular ecology, and subsequently benefits developing modern pest control strategies.     

De novo genome assembly of the stress tolerant forest species Casuarina equisetifolia provides insight into secondary growth

Casuarina equisetifolia (C. equisetifolia), a conifer‐like angiosperm with resistance to typhoon and stress tolerance, is mainly cultivated in the coastal areas of Australasia. C. equisetifolia, making it a valuable model to study secondary growth associated genes and stress‐tolerance traits. However, the genome sequence is unavailable and therefore wood‐associated growth rate and stress resistance at the molecular level is largely unexplored. We therefore constructed a high‐quality draft genome sequence of C. equisetifolia by a combination of Illumina second‐generation sequencing reads and Pacific Biosciences single‐molecule real‐time (SMRT) long reads to advance the investigation of this species. Here, we report the genome assembly, which contains approximately 300 megabases (Mb) and scaffold size of N50 is 1.06 Mb. Additionally, gene annotation, assisted by a combination of prediction and RNA‐seq data, generated 29 827 annotated protein‐coding genes and 1983 non‐coding genes, respectively. Furthermore, we found that the total number of repetitive sequences account for one‐third of the genome assembly. Here we also construct the genome‐wide map of DNA modification, such as two novel forms N⁶‐adenine (6mA) and N4‐methylcytosine (4mC) at the level of single‐nucleotide resolution using single‐molecule real‐time (SMRT) sequencing. Interestingly, we found that 17% of 6mA modification genes and 15% of 4mC modification genes also included alternative splicing events. Finally, we investigated cellulose, hemicellulose, and lignin‐related genes, which were associated with secondary growth and contained different DNA modifications. The high‐quality genome sequence and annotation of C. equisetifolia in this study provide a valuable resource to strengthen our understanding of the diverse traits of trees.     

Frequency of indicator bacteria,Salmonella and diarrhoeagenic Escherichia coli pathotypes on ready‐to‐eat cooked vegetable salads from Mexican restaurants

The presence of coliform bacteria, faecal coliforms, Escherichia coli, diarrhoeagenic E. coli pathotypes (DEP) and Salmonella were determined in ready‐to‐eat cooked vegetable salads (RECS) from restaurants in Pachuca city, Mexico. The RECS were purchased from three types of restaurants: national chain restaurants (A), local restaurants (B) and small restaurants (C). Two restaurants for each A and B, and three for C, were included. Forty RECS samples were purchased at each A and B restaurant and 20 at each C restaurant. Of the overall total of 220 analysed samples, 100, 98·2, 72·3, 4·1 and 4·1% had coliform bacteria, faecal coliforms, E. coli, DEP and Salmonella, respectively. Identified DEP included enteropathogenic E. coli (EPEC), enterotoxigenic E. coli (ETEC) and Shiga toxin‐producing E. coli (STEC). The EPEC, ETEC and STEC were isolated each from 1·4% of samples. No E. coli O157:H7 were detected in any STEC‐positive samples. The analysis of Kruskal–Wallis anova and median test of microbiological data showed that the microbiological quality of RECS did not differ between the different restaurants (P > 0·05).

Significance and Impact of the Study

This is the first report regarding microbiological quality and Salmonella, enteropathogenic E. coli (EPEC), enterotoxigenic E. coli (ETEC) and Shiga toxin‐producing E. coli (STEC) isolation from ready‐to‐eat cooked vegetable salads from Mexican restaurants. Ready‐to‐eat cooked vegetable salads could be an important factor contributing to the endemicity of EPEC, ETEC and STEC, and Salmonella caused gastroenteritis in Mexico.     

Molecular monitoring of Escherichia coli O157: H7 sterilization rate using qPCR and propidium monoazide treatment

Propidium monoazide is a DNA‐intercalating dye. PMA‐qPCR has been reported as a novel method to detect live bacteria in complex samples. In this study, this method was used to monitor the sterilization effects of UHP, ultrasound and high PEF on Escherichia coli O157:H7. Our results showed that all three sterilization techniques are successful to kill viable E. coli O157:H7 cells under their appropriate conditions. PMA‐qPCR can effectively monitor the amount of DNA released from viable E. coli O157:H7 cells, and the results from PMA‐qPCR were highly consistent with those from plate counting after treatment with UHP, ultrasound and high PEF. The maximal ΔC_t between PMA‐qPCR and qPCR obtained in this study was 10·39 for UHP, 5·76 for ultrasound and 2·30 for high PEF. The maximal sterilization rates monitored by PMA‐qPCR were 99·92% for UHP, 99·99% for ultrasound and 100% for high PEF. Thus, PMA‐qPCR can be used to detect the sterilization effect on food and water supplies after treatment with UHP, ultrasound and high PEF.

Significance and Impact of the Study

The reliable detection of viable foodborne pathogenic bacteria in water and food is of great importance in our daily life. However, the traditional bacteria cultivation‐based methods are time‐consuming and difficult to monitor all viable bacteria because of the limitation of cultivation conditions. This study demonstrated that PMA‐qPCR technique is very effective to monitor viable E. coli O157:H7 after sterilization and will help to monitor the viable bacteria in food and water.     

Chromosome‐level genome assembly of Triplophysa tibetana,a fish adapted to the harsh high‐altitude environment of the Tibetan Plateau

Triplophysa is an endemic fish genus of the Tibetan Plateau in China. Triplophysa tibetana, which lives at a recorded altitude of ~4,000 m and plays an important role in the highland aquatic ecosystem, serves as an excellent model for investigating high‐altitude environmental adaptation. However, evolutionary and conservation studies of T. tibetana have been limited by scarce genomic resources for the genus Triplophysa. In the present study, we applied PacBio sequencing and the Hi‐C technique to assemble the T. tibetana genome. A 652‐Mb genome with 1,325 contigs with an N50 length of 3.1 Mb was obtained. The 1,137 contigs were further assembled into 25 chromosomes, representing 98.7% and 80.47% of all contigs at the base and sequence number level, respectively. Approximately 260 Mb of sequence, accounting for ~39.8% of the genome, was identified as repetitive elements. DNA transposons (16.3%), long interspersed nuclear elements (12.4%) and long terminal repeats (11.0%) were the most repetitive types. In total, 24,372 protein‐coding genes were predicted in the genome, and ~95% of the genes were functionally annotated via a search in public databases. Using whole genome sequence information, we found that T. tibetana diverged from its common ancestor with Danio rerio ~121.4 million years ago. The high‐quality genome assembled in this work not only provides a valuable genomic resource for future population and conservation studies of T. tibetana, but it also lays a solid foundation for further investigation into the mechanisms of environmental adaptation of endemic fishes in the Tibetan Plateau.     

The complete chloroplast genome of Myriophyllum spicatum reveals a 4‐kb inversion and new insights regarding plastome evolution in Haloragaceae

Myriophyllum, among the most species‐rich genera of aquatic angiosperms with ca. 68 species, is an extensively distributed hydrophyte lineage in the cosmopolitan family Haloragaceae. The chloroplast (cp) genome is useful in the study of genetic evolution, phylogenetic analysis, and molecular dating of controversial taxa. Here, we sequenced and assembled the whole chloroplast genome of Myriophyllum spicatum L. and compared it to other species in the order Saxifragales. The complete chloroplast genome sequence of M. spicatum is 158,858 bp long and displays a quadripartite structure with two inverted repeats (IR) separating the large single copy (LSC) region from the small single copy (SSC) region. Based on sequence identification and the phylogenetic analysis, a 4‐kb phylogenetically informative inversion between trnE‐trnC in Myriophyllum was determined, and we have placed this inversion on a lineage specific to Myriophyllum and its close relatives. The divergence time estimation suggested that the trnE‐trnC inversion possibly occurred between the upper Cretaceous (72.54 MYA) and middle Eocene (47.28 MYA) before the divergence of Myriophyllum from its most recent common ancestor. The unique 4‐kb inversion might be caused by an occurrence of nonrandom recombination associated with climate changes around the K‐Pg boundary, making it interesting for future evolutionary investigations.     

Rapid proliferation and nucleolar organizer targeting centromeric retrotransposons in cotton

Centromeric chromatin in most eukaryotes is composed of highly repetitive centromeric retrotransposons and satellite repeats that are highly variable even among closely related species. The evolutionary mechanisms that underlie the rapid evolution of centromeric repeats remain unknown. To obtain insight into the evolution of centromeric repeats following polyploidy, we studied a model diploid progenitor (Gossypium raimondii, D‐genome) of the allopolyploid (AD‐genome) cottons, G. hirsutum and G. barbadense. Sequence analysis of chromatin‐immunoprecipitated DNA showed that the G. raimondii centromeric repeats originated from retrotransposon‐related sequences. Comparative analysis showed that nine of the 10 analyzed centromeric repeats were absent from the centromeres in the A‐genome and related diploid species (B‐, F‐ and G‐genomes), indicating that they colonized the centromeres of D‐genome lineage after the divergence of the A‐ and D‐ ancestral species or that they were ancestrally retained prior to the origin of Gossypium. Notably, six of the nine repeats were present in both the A‐ and D‐subgenomes in tetraploid G. hirsutum, and increased in abundance in both subgenomes. This finding suggests that centromeric repeats may spread and proliferate between genomes subsequent to polyploidization. Two repeats, Gr334 and Gr359 occurred in both the centromeres and nucleolar organizer regions (NORs) in D‐ and AD‐genome species, yet localized to just the NORs in A‐, B‐, F‐, and G‐genome species. Contained within is a story of an established centromeric repeat that is eliminated and allopolyploidization provides an opportunity for reinvasion and reestablishment, which broadens our evolutionary understanding behind the cycles of centromeric repeat establishment and targeting.