Identification of Signatures of Positive Selection That Have Shaped the Genomic Landscape of South African Pig Populations

South Africa boasts a diverse range of pig populations, encompassing intensively raised commercial breeds, as well as indigenous and village pigs reared under low-input production systems. The aim of this study was to investigate how natural and artificial selection have shaped the genomic landscape of South African pig populations sampled from different genetic backgrounds and production systems. For this purpose, the integrated haplotype score (iHS), as well as cross population extended haplotype homozygosity (XP-EHH) and Lewontin and Krakauer's extension of the Fst statistic based on haplotype information (HapFLK) were utilised. Our results revealed several population-specific signatures of selection associated with the different production systems. The importance of natural selection in village populations was highlighted, as the majority of genomic regions under selection were identified in these populations. Regions under natural and artificial selection causing the distinct genetic footprints of these populations also allow for the identification of genes and pathways that may influence production and adaptation. In the context of intensively raised commercial pig breeds (Large White, Kolbroek, and Windsnyer), the identified regions included quantitative loci (QTLs) associated with economically important traits. For example, meat and carcass QTLs were prevalent in all the populations, showing the potential of village and indigenous populations' ability to be managed and improved for such traits. Results of this study therefore increase our understanding of the intricate interplay between selection pressures, genomic adaptations, and desirable traits within South African pig populations.

Resistome, mobilome, virulome analysis and phylogenomics of Enterococcus faecalis isolated from raw muscle foods of beef origin in Gauteng, South Africa

Enterococcus faecalis is a ubiquitous bacterium found in various environments, including processed beef meat, and is known for its importance in both food safety and public health. This pivotal significance stems not solely from its virulence but also from its adeptness in eliciting multidrug-resistant infections in humans. The aim of this study was to investigate the population structure, resistome, mobilome, and virulome of E. faecalis obtained from processed beef meat sources in South Africa. A total of eight genomes sequenced in this study were examined, alongside 78 publicly available, high-quality genomes of E. faecalis, with a comprehensive analysis conducted to identify antimicrobial resistance (AMR) determinants, virulence factors, and mobile genetic elements (MGE). Six distinct sequence types (STs) (ST79, ST860, ST40, ST238, ST21, and ST700) and 41 core virulence factors were found across all the genomes. The virulence factors included genes encoding adherence (ace, asa1, Ef0485, ebpA, ebpB, ebpC, srtC); exoenzyme (Ef3023, Ef0818, gelE, sprE); immunomodulation (cpsA, cpsB, cpsC, cpsD, cpsE, cpsF, cpsG, cpsH, cpsI, cpsK), and biofilm formation (bopD, fsrA, fsrB, fsrC). In addition, AMR genes were identified across all genomes, which include aminoglycoside resistance (ant(6)-Ia), trimethoprim resistance (dfrA), drug and biocide resistance (efrA and efrB), multidrug efflux pump (emeA), clindamycin quinupristin-dalfopristin, dalfopristin resistance (lsaA), and tetracycline resistance (tetM). The genomes of E. faecalis sequenced here contained a variety of MGEs, including Insertion Sequences (ISs), transposons, prophages, and plasmids, which may have facilitated genetic exchange within and between these species. The results highlight that beef meat products act as a reservoir for virulent E. faecalis strains possessing antibiotic-resistance traits. This study provides insight into the genomic characteristics, antimicrobial resistance genes, virulence factors, and genetic mobile elements associated with eight E. faecalis isolates from processed beef meat in the Gauteng province of South Africa.

Genome sequence of an indole-3-acetic acid and siderophores producing Pseudomonas monsensis SARCC-3054 with the potential as a microbial inoculant for enhanced crop production

The genome of Pseudomonas monsensis strain SARCC-3054 was sequenced after being confirmed as a potential plant growth-promoting rhizobacteria in both in vitro and in vivo assays. The 6.3 MB genome has a GC content of 60.2% and is divided into 59 contigs that contain several plant beneficial genes and proteins.

Genus-wide genomic characterization of Macrococcus: insights into evolution, population structure, and functional potential

Introduction

Macrococcus species have been isolated from a range of mammals and mammal-derived food products. While they are largely considered to be animal commensals, Macrococcus spp. can be opportunistic pathogens in both veterinary and human clinical settings. This study aimed to provide insight into the evolution, population structure, and functional potential of the Macrococcus genus, with an emphasis on antimicrobial resistance (AMR) and virulence potential.

Methods

All high-quality, publicly available Macrococcus genomes (n = 104, accessed 27 August 2022), plus six South African genomes sequenced here (two strains from bovine clinical mastitis cases and four strains from beef products), underwent taxonomic assignment (using four different approaches), AMR determinant detection (via AMRFinderPlus), and virulence factor detection (using DIAMOND and the core Virulence Factor Database).

Results

Overall, the 110 Macrococcus genomes were of animal commensal, veterinary clinical, food-associated (including food spoilage), and environmental origins; five genomes (4.5%) originated from human clinical cases. Notably, none of the taxonomic assignment methods produced identical results, highlighting the potential for Macrococcus species misidentifications. The most common predicted antimicrobial classes associated with AMR determinants identified across Macrococcus included macrolides, beta-lactams, and aminoglycosides (n = 81, 61, and 44 of 110 genomes; 73.6, 55.5, and 40.0%, respectively). Genes showing homology to Staphylococcus aureus exoenzyme aureolysin were detected across multiple species (using 90% coverage, n = 40 and 77 genomes harboring aureolysin-like genes at 60 and 40% amino acid [AA] identity, respectively). S. aureus Panton-Valentine leucocidin toxin-associated lukF-PV and lukS-PV homologs were identified in eight M. canis genomes (≥40% AA identity, >85% coverage). Using a method that delineates populations using recent gene flow (PopCOGenT), two species (M. caseolyticus and M. armenti) were composed of multiple within-species populations. Notably, M. armenti was partitioned into two populations, which differed in functional potential (e.g., one harbored beta-lactamase family, type II toxin-antitoxin system, and stress response proteins, while the other possessed a Type VII secretion system; PopCOGenT p < 0.05).

Discussion

Overall, this study leverages all publicly available Macrococcus genomes in addition to newly sequenced genomes from South Africa to identify genomic elements associated with AMR or virulence potential, which can be queried in future experiments.

Draft Genome Sequence of Pectobacterium brasiliense Isolated from Potato in South Africa

This study reports a draft genome of a phytopathogenic bacterium, Pectobacterium brasiliense, isolated from potato in South Africa. The total reported length of the genome is 4,897,858 bp, contained in 172 contigs with 4,378 genes. The GC content of the genome is 51.6%.

Multiple-Locus Variable-Number Tandem Repeat Analysis Genotypes of Listeria monocytogenes Isolated from Farms, Abattoirs, and Retail in Gauteng Province, South Africa

ABSTRACT

The use of multiple-locus variable-number analysis (MLVA) of tandem repeats (TRs) for subtyping Listeria monocytogenes has proven to be reliable and fast. This study determined the MLVA genotypes of 60 isolates of L. monocytogenes recovered from cattle farms, abattoirs, and retail outlets in Gauteng province, South Africa. The distribution of the 60 L. monocytogenes isolates analyzed by type of sample was as follows: raw beef (28, 46.7%), ready-to-eat beef products (9, 15.0%), beef carcass swabs (9, 15.0%), cattle environment (6, 10.0%), and cattle feces (8, 13.3%). The serogroups of the isolates were determined using PCR and the MLVA genotypes based on six selected loci. The frequency of the 60 serogroups detected was as follows: 1/2a-3a (IIa) (27, 45.0%); 4b-4d-4e (1Vb) (24, 40.0%); 1/2c-3c (IIc) (8, 13.3%); and 1/2b-3b (IIb) (1, 1.7%). MLVA successfully clustered genetically related isolates and differentiated nonrelated isolates, irrespective of their sources, sample types, and serogroups, as demonstrated by 16 MLVA pattern types detected. For serogroup 4b-4d-4e (IVb), there was no variation in TRs LM-TR2, LM-TR4, and LM-TR6, which each contained only one allele (02, 00, and 93, respectively). However, across the sources and sample types of isolates, there was variation in serogroup 4b-4d-4e (IVb): LM-TR1 contained 00, 03, and 05; LM-TR3 contained 14, 20, and 22; and LM-TR5 contained 14, 21, and 25. Similar patterns of variation in the TRs were detected in the other serogroups (1/2a-3a, 1/2b-3b, and 1/2c-3c). BioNumeric data analysis identified at least five types in Gauteng province. MLVA epidemiologically clustered the related isolates and differentiated unrelated isolates.

HIGHLIGHTS

Genomic Sequencing of Bacillus cereus Sensu Lato Strains Isolated from Meat and Poultry Products in South Africa Enables Inter- and Intranational Surveillance and Source Tracking

Members of the Bacillus cereus sensu lato species complex, also known as the B. cereus group, vary in their ability to cause illness but are frequently isolated from foods, including meat products; however, food safety surveillance efforts that use whole-genome sequencing (WGS) often neglect these potential pathogens. Here, we evaluate the surveillance and source tracking potential of WGS as applied to B. cereus sensu lato by (i) using WGS to characterize B. cereus sensu lato strains isolated during routine surveillance of meat products across South Africa (n = 25) and (ii) comparing the genomes sequenced here to all publicly available, high-quality B. cereus sensu lato genomes (n = 2,887 total genomes). Strains sequenced here were collected from meat products obtained from (i) retail outlets, processing plants, and butcheries across six South African provinces (n = 23) and (ii) imports held at port of entry (n = 2). The 25 strains sequenced here were partitioned into 15 lineages via in silico seven-gene multilocus sequence typing (MLST). While none of the South African B. cereus sensu lato strains sequenced here were identical to publicly available genomes, six MLST lineages contained multiple strains sequenced in this study, which were identical or nearly identical at the whole-genome scale (≤3 core single nucleotide polymorphisms). Five MLST lineages contained (nearly) identical genomes collected from two or three South African provinces; one MLST lineage contained nearly identical genomes from two countries (South Africa and the Netherlands), indicating that B. cereus sensu lato can spread intra- and internationally via foodstuffs. IMPORTANCE Nationwide foodborne pathogen surveillance programs that use high-resolution genomic methods have been shown to provide vast public health and economic benefits. However, Bacillus cereus sensu lato is often overlooked during large-scale routine WGS efforts. Thus, to our knowledge, no studies to date have evaluated the potential utility of WGS for B. cereus sensu lato surveillance and source tracking in foodstuffs. In this preliminary proof-of-concept study, we applied WGS to B. cereus sensu lato strains collected via South Africa's national surveillance program of domestic and imported meat products, and we provide strong evidence that B. cereus sensu lato can be disseminated intra- and internationally via the agro-food supply chain. Our results showcase that WGS has the potential to be used for source tracking of B. cereus sensu lato in foods, although future WGS and metadata collection efforts are needed to ensure that B. cereus sensu lato surveillance initiatives are on par with those of other foodborne pathogens.

Metatranscriptomic Analyses Reveal the Functional Role of Botrytis cinerea in Biochemical and Textural Changes during Noble Rot of Grapevines

Botrytis cinerea, can lead to the formation of noble rot (NR) of grape berries under certain environmental conditions, resulting in favored metabolic and physical changes necessary for producing highly regarded botrytized wines. The functional genes involved in the textural and biochemical processes are still poorly characterized. We generated and analyzed metatranscriptomic data from healthy (H) berries and from berries representing the four stages of NR from the Tokaj wine region in Hungary over three months. A weighted gene co-expression network analysis (WGCNA) was conducted to link B. cinerea functional genes to grape berry physical parameters berry hardness (BH), berry skin break force (F_sk), berry skin elasticity (E_sk), and the skin break energy (W_sk). Clustered modules showed that genes involved in carbohydrate and protein metabolism were significantly enriched in NR, highlighting their importance in the grape berry structural integrity. Carbohydrate active enzymes were particularly up-regulated at the onset of NR (during the transition from phase I to II) suggesting that the major structural changes occur early in the NR process. In addition, we identified genes expressed throughout the NR process belonging to enriched pathways that allow B. cinerea to dominate and proliferate during this state, including sulphate metabolizing genes and genes involved in the synthesis of antimicrobials.

A 16S Next Generation Sequencing Based Molecular and Bioinformatics Pipeline to Identify Processed Meat Products Contamination and Mislabelling

Simple Summary

Meat adulteration and fraud encompasses the deliberate fraudulent addition or substitution of proteins of animal or plant origin in edible products primarily for economic gain. The mitochondrial 16S ribosomal (rRNA) gene was used to identify species that are present in pure and processed meat samples. The meat samples were sequenced using an Illumina sequencing platform, and bioinformatics analysis was carried out for species identification. The results indicated that pork was the major contaminant in most of the meat samples. The bioinformatics pipeline demonstrated its specificity through identification of species specific and quantification of the contamination levels across all samples. Food business operators and regulatory sectors can validate this method for food fraud checks and manage any form of mislabeling in the animal or plant protein food ecosystem.

Abstract

Processed meat is a target in meat adulteration for economic gain. This study demonstrates a molecular and bioinformatics diagnostic pipeline, utilizing the mitochondrial 16S ribosomal RNA (rRNA) gene, to determine processed meat product mislabeling through Next-Generation Sequencing. Nine pure meat samples were collected and artificially mixed at different ratios to verify the specificity and sensitivity of the pipeline. Processed meat products (n = 155), namely, minced meat, biltong, burger patties, and sausages, were collected across South Africa. Sequencing was performed using the Illumina MiSeq sequencing platform. Each sample had paired-end reads with a length of ±300 bp. Quality control and filtering was performed using BBDuk (version 37.90a). Each sample had an average of 134,000 reads aligned to the mitochondrial genomes using BBMap v37.90. All species in the artificial DNA mixtures were detected. Processed meat samples had reads that mapped to the Bos (90% and above) genus, with traces of reads mapping to Sus and Ovis (2–5%) genus. Sausage samples showed the highest level of contamination with 46% of the samples having mixtures of beef, pork, or mutton in one sample. This method can be used to authenticate meat products, investigate, and manage any form of mislabeling.

Genomic Characterization of Endemic and Ecdemic Non-typhoidal Salmonella enterica Lineages Circulating Among Animals and Animal Products in South Africa

In Africa, the burden of illness caused by non-typhoidal Salmonella enterica is disproportionally high; however, whole-genome sequencing (WGS) efforts are overwhelmingly concentrated in world regions with lower burdens. While WGS is being increasingly employed in South Africa to characterize Salmonella enterica, the bulk of these efforts have centered on characterizing human clinical strains. Thus, very little is known about lineages circulating among animals in the country on a genomic scale. Here, we used WGS to characterize 63 Salmonella enterica strains isolated from livestock, companion animals, wildlife, and animal products in South Africa over a 60-year period. Genomes were assigned to serotypes Dublin, Hadar, Enteritidis, and Typhimurium (n = 18, 8, 13, and 24 strains, respectively) and sequence types (STs) ST10 (all S. Dublin), ST33 (all S. Hadar), ST11/ST366 (n = 12 and 1 S. Enteritidis, respectively), and ST19/ST34 (n = 23 and 1 S. Typhimurium, respectively; via seven-gene multi-locus sequence typing). Within-ST phylogenies were constructed using genomes sequenced in this study, plus publicly available genomes representative of each ST's (i) global (n = 2,802 and 1,569 S. Dublin and Hadar genomes, respectively) and (ii) African (n = 716 and 343 S. Enteritidis and Typhimurium genomes, respectively) population. For S. Dublin ST10, a largely antimicrobial-susceptible, endemic lineage circulating among humans, animals, and food in South Africa was identified, as well as a lineage that was likely recently introduced from the United States. For S. Hadar ST33, multiple South African lineages harboring streptomycin and tetracycline resistance-conferring genes were identified. African S. Enteritidis ST11 could be primarily partitioned into one largely antimicrobial-susceptible and one largely multidrug-resistant (MDR) clade, with South African isolates confined to the largely antimicrobial-susceptible clade. S. Typhimurium ST19/ST34 strains sequenced here were distributed across the African S. Typhimurium ST19/ST34 phylogeny, representing a diverse range of lineages, including numerous MDR lineages. Overall, this study provides critical insights into endemic and ecdemic non-typhoidal Salmonella enterica lineages circulating among animals, foods, and humans in South Africa and showcases the utility of WGS in characterizing animal-associated strains from a world region with a high salmonellosis burden.

Assessing single-nucleotide polymorphism selection methods for the development of a low-density panel optimized for imputation in South African Drakensberger beef cattle

A major obstacle in applying genomic selection (GS) to uniquely adapted local breeds in less-developed countries has been the cost of genotyping at high densities of single-nucleotide polymorphisms (SNP). Cost reduction can be achieved by imputing genotypes from lower to higher densities. Locally adapted breeds tend to be admixed and exhibit a high degree of genomic heterogeneity thus necessitating the optimization of SNP selection for downstream imputation. The aim of this study was to quantify the achievable imputation accuracy for a sample of 1,135 South African (SA) Drakensberger cattle using several custom-derived lower-density panels varying in both SNP density and how the SNP were selected. From a pool of 120,608 genotyped SNP, subsets of SNP were chosen (1) at random, (2) with even genomic dispersion, (3) by maximizing the mean minor allele frequency (MAF), (4) using a combined score of MAF and linkage disequilibrium (LD), (5) using a partitioning-around-medoids (PAM) algorithm, and finally (6) using a hierarchical LD-based clustering algorithm. Imputation accuracy to higher density improved as SNP density increased; animal-wise imputation accuracy defined as the within-animal correlation between the imputed and actual alleles ranged from 0.625 to 0.990 when 2,500 randomly selected SNP were chosen vs. a range of 0.918 to 0.999 when 50,000 randomly selected SNP were used. At a panel density of 10,000 SNP, the mean (standard deviation) animal-wise allele concordance rate was 0.976 (0.018) vs. 0.982 (0.014) when the worst (i.e., random) as opposed to the best (i.e., combination of MAF and LD) SNP selection strategy was employed. A difference of 0.071 units was observed between the mean correlation-based accuracy of imputed SNP categorized as low (0.01 < MAF ≤ 0.1) vs. high MAF (0.4 < MAF ≤ 0.5). Greater mean imputation accuracy was achieved for SNP located on autosomal extremes when these regions were populated with more SNP. The presented results suggested that genotype imputation can be a practical cost-saving strategy for indigenous breeds such as the SA Drakensberger. Based on the results, a genotyping panel consisting of ~10,000 SNP selected based on a combination of MAF and LD would suffice in achieving a <3% imputation error rate for a breed characterized by genomic admixture on the condition that these SNP are selected based on breed-specific selection criteria.

South African Buffalo-Derived Theileria parva Is Distinct From Other Buffalo and Cattle-Derived T. parva

Theileria parva is a protozoan parasite transmitted by the brown-eared ticks, Rhipicephalus appendiculatus and Rhipicephalus zambeziensis. Buffaloes are the parasite's ancestral host, with cattle being the most recent host. The parasite has two transmission modes namely, cattle-cattle and buffalo-cattle transmission. Cattle-cattle T. parva transmission causes East Coast fever (ECF) and January disease syndromes. Buffalo to cattle transmission causes Corridor disease. Knowledge on the genetic diversity of South African T. parva populations will assist in determining its origin, evolution and identify any cattle-cattle transmitted strains. To achieve this, genomic DNA of blood and in vitro culture material infected with South African isolates (8160, 8301, 8200, 9620, 9656, 9679, Johnston, KNP2, HL3, KNP102, 9574, and 9581) were extracted and paired-end whole genome sequencing using Illumina HiSeq 2500 was performed. East and southern African sample data (Chitongo Z2, Katete B2, Kiambu Z464/C12, Mandali Z22H10, Entebbe, Nyakizu, Katumba, Buffalo LAWR, and Buffalo Z5E5) was also added for comparative purposes. Data was analyzed using BWA and SAMtools variant calling with the T. parva Muguga genome sequence used as a reference. Buffalo-derived strains had higher genetic diversity, with twice the number of variants compared to cattle-derived strains, confirming that buffaloes are ancestral reservoir hosts of T. parva. Host specific SNPs, however, could not be identified among the selected 74 gene sequences. Phylogenetically, strains tended to cluster by host with South African buffalo-derived strains clustering with buffalo-derived strains. Among the buffalo-derived strains, South African strains were genetically divergent from other buffalo-derived strains indicating possible geographic sub-structuring. Geographic sub- structuring was also observed within South Africa strains. The knowledge generated from this study indicates that to date, ECF is not circulating in buffalo from South Africa. It also shows that T. parva has historically been present in buffalo from South Africa before the introduction of ECF and was not introduced into buffalo during the ECF epidemic.

Microbial Communities of Meat and Meat Products: An Exploratory Analysis of the Product Quality and Safety at Selected Enterprises in South Africa

Consumption of food that is contaminated by microorganisms, chemicals, and toxins may lead to significant morbidity and mortality, which has negative socioeconomic and public health implications. Monitoring and surveillance of microbial diversity along the food value chain is a key component for hazard identification and evaluation of potential pathogen risks from farm to the consumer. The aim of this study was to determine the microbial diversity in meat and meat products from different enterprises and meat types in South Africa. Samples (n = 2017) were analyzed for Yersinia enterocolitica, Salmonella species, Listeria monocytogenes, Campylobacter jejuni, Campylobacter coli, Staphylococcus aureus, Clostridium perfringens, Bacillus cereus, and Clostridium botulinum using culture-based methods. PCR was used for confirmation of selected pathogens. Of the 2017 samples analyzed, microbial ecology was assessed for selected subsamples where next generation sequencing had been conducted, followed by the application of computational methods to reconstruct individual genomes from the respective sample (metagenomics). With the exception of Clostridium botulinum, selective culture-dependent methods revealed that samples were contaminated with at least one of the tested foodborne pathogens. The data from metagenomics analysis revealed the presence of diverse bacteria, viruses, and fungi. The analyses provide evidence of diverse and highly variable microbial communities in products of animal origin, which is important for food safety, food labeling, biosecurity, and shelf life limiting spoilage by microorganisms.

Q Fever: Seroprevalence, Risk Factors in Slaughter Livestock and Genotypes of Coxiella burnetii in South Africa

Q fever is a neglected zoonosis in South Africa, causing significant losses in livestock and game animals through reproductive disorders. However, there are limited studies on the extent of Coxiella burnetii infections in livestock in South Africa. Further, there is also lack of knowledge about the types of C. burnetii strains that are currently circulating in the country. Therefore, a cross-sectional, abattoir-based study was conducted to determine the seroprevalence of C. burnetii and associated risk factors, and to characterize C. burnetii strains from slaughter livestock at red meat abattoirs in Gauteng, South Africa. Of the 507 animals tested, 6.9% (95% CI: 4.9-9.5%) were positive for antibodies against C. burnetii. The seroprevalence was 9.4% (31/331) in cattle, 4.3% (3/69) in sheep, and 0.9% (1/107) in pigs. Out of the 63 tissue samples from 35 seropositive animals including material from two sheep aborted fetuses from Mangaung district (Free State province), 12.7% (8/63) tested positive by IS1111 PCR. Genotyping of the eight PCR-positive tissues from eight animals by MLVA revealed two novel genotypes, not available in Coxiella MLVA databases. It is concluded that slaughter animals pose a risk of exposing abattoir and farm workers to C. burnetii in South Africa.

Runs of homozygosity analysis of South African sheep breeds from various production systems investigated using OvineSNP50k data

Background

Population history, production system and within-breed selection pressure impacts the genome architecture resulting in reduced genetic diversity and increased frequency of runs of homozygosity islands. This study tested the hypothesis that production systems geared towards specific traits of importance or natural or artificial selection pressures influenced the occurrence and distribution of runs of homozygosity (ROH) in the South African sheep population. The Illumina OvineSNP50 BeadChip was used to genotype 400 sheep belonging to 13 breeds from South Africa representing mutton, pelt and mutton and wool dual-purpose breeds, including indigenous non-descript breeds that are reared by smallholder farmers. To get more insight into the autozygosity and distribution of ROH islands of South African breeds relative to global populations, 623 genotypes of sheep from worldwide populations were included in the analysis. Runs of homozygosity were computed at cut-offs of 1–6 Mb, 6–12 Mb, 12–24 Mb, 24–48 Mb and > 48 Mb, using the R package detectRUNS. The Golden Helix SVS program was used to investigate the ROH islands.

Results

A total of 121,399 ROH with mean number of ROH per animal per breed ranging from 800 (African White Dorper) to 15,097 (Australian Poll Dorset) were obtained. Analysis of the distribution of ROH according to their size showed that, for all breeds, the majority of the detected ROH were in the short (1–6 Mb) category (88.2%). Most animals had no ROH > 48 Mb. Of the South African breeds, the Nguni and the Blackhead Persian displayed high ROH based inbreeding (F_ROH) of 0.31 ± 0.05 and 0.31 ± 0.04, respectively. Highest incidence of common runs per SNP across breeds was observed on chromosome 10 with over 250 incidences of common ROHs. Mean proportion of SNPs per breed per ROH island ranged from 0.02 ± 0.15 (island ROH224 on chromosome 23) to 0.13 ± 0.29 (island ROH175 on chromosome 15). Seventeen (17) of the islands had SNPs observed in single populations (unique ROH islands). The MacArthur Merino (MCM) population had five unique ROH islands followed by Blackhead Persian and Nguni with three each whilst the South African Mutton Merino, SA Merino, White Vital Swakara, Karakul, Dorset Horn and Chinese Merino each had one unique ROH island. Genes within ROH islands were associated with predominantly metabolic and immune response traits and predomestic selection for traits such as presence or absence of horns.

Conclusions

Overall, the frequency and patterns of distribution of ROH observed in this study corresponds to the breed history and implied selection pressures exposed to the sheep populations under study.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-020-07314-2.

Microbial and functional diversity of Cyclopia intermedia rhizosphere microbiome revealed by analysis of shotgun metagenomics sequence data

Cyclopia spp., commonly referred to as honeybush due to the honey scented flowers, are indigenous legumes mainly growing in the Cape Floristic Region of the Western Cape, South Africa. Dozens of species, including Cyclopia intermedia, C. subternata, C. plicata, C. genistoides are used to make the well-known, popular and widely enjoyed beverage called 'honeybush tea'. In the past, most rhizosphere microbial studies associated with Cyclopia spp. focused mainly on the taxonomy and diversity of the root nodule associated symbiotic nitrogen fixing rhizobia. The work presented here is the first report on the microbial and functional diversity of rhizosphere microbiome associated with Cyclopia intermedia. Metagenomic shotgun sequencing was performed on the rhizosphere soil sample collected from this Cyclopia sp. using illumina Hiseq 2500 platform which resulted in an α- diversity of 312 species. Analysis of the metagenome sequence using the Metagenomic analysis server (MG-RAST) indicated that bacteria constitute the dominant domain followed by Eukaryota, Archaea and other sequences derived from fungi and viruses. Functional diversity of the metagenome based on analysis using the Cluster Orthologous Group (COG) method showed metabolism as the most important function in the community. The raw sequence data is uploaded in FASTQ format on MG-RAST server with ID mgm4855911.3 which can be accessed at http://www.mg-rast.org/linkin.cgi?project=mgp90368. The data on the microbial and functional diversity of the rhizosphere community of Cyclopia intermedia generates a baseline information about the microbial ecology of this indigenous legume. The microbial profile data can also be used as indicators of soil health characteristic of the rhizosphere of this important legume.

Population Structure of Non-ST6 Listeria monocytogenes Isolated in the Red Meat and Poultry Value Chain in South Africa

Meat products have been implicated in many listeriosis outbreaks globally, however there is a dearth of information on the diversity of L. monocytogenes isolates circulating in food products in South Africa. The aim of this study was to investigate the population structure of L. monocytogenes isolated in the meat value chain within the South African market. Based on whole-genome sequence analysis, a total of 217 isolates were classified into two main lineage groupings namely lineages I (n = 97; 44.7%) and II (n = 120; 55.3%). The lineage groups were further differentiated into IIa (n = 95, 43.8%), IVb (n = 69, 31.8%), IIb (n = 28, 12.9%), and IIc (n = 25, 11.5%) sero-groups. The most abundant sequence types (STs) were ST204 (n = 32, 14.7%), ST2 (n = 30, 13.8%), ST1 (n = 25, 11.5%), ST9 (n = 24, 11.1%), and ST321 (n = 21, 9.7%). In addition, 14 clonal complex (CCs) were identified with over-representation of CC1, CC3, and CC121 in "Processed Meat-Beef", "RTE-Poultry", and "Raw-Lamb" meat categories, respectively. Listeria pathogenic islands were present in 7.4% (LIPI-1), 21.7% (LIPI-3), and 1.8% (LIPI-4) of the isolates. Mutation leading to premature stop codons was detected in inlA virulence genes across isolates identified as ST121 and ST321. The findings of this study demonstrated a high-level of genomic diversity among L. monocytogenes isolates recovered across the meat value chain control points in South Africa.

Phages Actively Challenge Niche Communities in Antarctic Soils

By modulating the structure, diversity, and trophic outputs of microbial communities, phages play crucial roles in many biomes. In oligotrophic polar deserts, the effects of katabatic winds, constrained nutrients, and low water availability are known to limit microbial activity. Although phages may substantially govern trophic interactions in cold deserts, relatively little is known regarding the precise ecological mechanisms. Here, we provide the first evidence of widespread antiphage innate immunity in Antarctic environments using metagenomic sequence data from hypolith communities as model systems. In particular, immunity systems such as DISARM and BREX are shown to be dominant systems in these communities. Additionally, we show a direct correlation between the CRISPR-Cas adaptive immunity and the metavirome of hypolith communities, suggesting the existence of dynamic host-phage interactions. In addition to providing the first exploration of immune systems in cold deserts, our results suggest that phages actively challenge niche communities in Antarctic polar deserts. We provide evidence suggesting that the regulatory role played by phages in this system is an important determinant of bacterial host interactions in this environment.IMPORTANCE In Antarctic environments, the combination of both abiotic and biotic stressors results in simple trophic levels dominated by microbiomes. Although the past two decades have revealed substantial insights regarding the diversity and structure of microbiomes, we lack mechanistic insights regarding community interactions and how phages may affect these. By providing the first evidence of widespread antiphage innate immunity, we shed light on phage-host dynamics in Antarctic niche communities. Our analyses reveal several antiphage defense systems, including DISARM and BREX, which appear to dominate in cold desert niche communities. In contrast, our analyses revealed that genes which encode antiphage adaptive immunity were underrepresented in these communities, suggesting lower infection frequencies in cold edaphic environments. We propose that by actively challenging niche communities, phages play crucial roles in the diversification of Antarctic communities.

De novo Assembly of Transcriptomes From a B73 Maize Line Introgressed With a QTL for Resistance to Gray Leaf Spot Disease Reveals a Candidate Allele of a Lectin Receptor-Like Kinase

Gray leaf spot (GLS) disease in maize, caused by the fungus Cercospora zeina, is a threat to maize production globally. Understanding the molecular basis for quantitative resistance to GLS is therefore important for food security. We developed a de novo assembly pipeline to identify candidate maize resistance genes. Near-isogenic maize lines with and without a QTL for GLS resistance on chromosome 10 from inbred CML444 were produced in the inbred B73 background. The B73-QTL line showed a 20% reduction in GLS disease symptoms compared to B73 in the field (p = 0.01). B73-QTL leaf samples from this field experiment conducted under GLS disease pressure were RNA sequenced. The reads that did not map to the B73 or C. zeina genomes were expected to contain novel defense genes and were de novo assembled. A total of 141 protein-coding sequences with B73-like or plant annotations were identified from the B73-QTL plants exposed to C. zeina. To determine whether candidate gene expression was induced by C. zeina, the RNAseq reads from C. zeina-challenged and control leaves were mapped to a master assembly of all of the B73-QTL reads, and differential gene expression analysis was conducted. Combining results from both bioinformatics approaches led to the identification of a likely candidate gene, which was a novel allele of a lectin receptor-like kinase named L-RLK-CML that (i) was induced by C. zeina, (ii) was positioned in the QTL region, and (iii) had functional domains for pathogen perception and defense signal transduction. The 817AA L-RLK-CML protein had 53 amino acid differences from its 818AA counterpart in B73. A second "B73-like" allele of L-RLK was expressed at a low level in B73-QTL. Gene copy-specific RT-qPCR confirmed that the l-rlk-cml transcript was the major product induced four-fold by C. zeina. Several other expressed defense-related candidates were identified, including a wall-associated kinase, two glutathione s-transferases, a chitinase, a glucan beta-glucosidase, a plasmodesmata callose-binding protein, several other receptor-like kinases, and components of calcium signaling, vesicular trafficking, and ethylene biosynthesis. This work presents a bioinformatics protocol for gene discovery from de novo assembled transcriptomes and identifies candidate quantitative resistance genes.

Selection of marker genes for genetic barcoding of microorganisms and binning of metagenomic reads by Barcoder software tools

Background

Metagenomic approaches have revealed the complexity of environmental microbiomes with the advancement in whole genome sequencing displaying a significant level of genetic heterogeneity on the species level. It has become apparent that patterns of superior bioactivity of bacteria applicable in biotechnology as well as the enhanced virulence of pathogens often requires distinguishing between closely related species or sub-species. Current methods for binning of metagenomic reads usually do not allow for identification below the genus level and generally stops at the family level.

Results

In this work, an attempt was made to improve metagenomic binning resolution by creating genome specific barcodes based on the core and accessory genomes. This protocol was implemented in novel software tools available for use and download from http://bargene.bi.up.ac.za/. The most abundant barcode genes from the core genomes were found to encode for ribosomal proteins, certain central metabolic genes and ABC transporters. Performance of metabarcode sequences created by this package was evaluated using artificially generated and publically available metagenomic datasets. Furthermore, a program (Barcoding 2.0) was developed to align reads against barcode sequences and thereafter calculate various parameters to score the alignments and the individual barcodes. Taxonomic units were identified in metagenomic samples by comparison of the calculated barcode scores to set cut-off values. In this study, it was found that varying sample sizes, i.e. number of reads in a metagenome and metabarcode lengths, had no significant effect on the sensitivity and specificity of the algorithm. Receiver operating characteristics (ROC) curves were calculated for different taxonomic groups based on the results of identification of the corresponding genomes in artificial metagenomic datasets. The reliability of distinguishing between species of the same genus or family by the program was nearly perfect.

Conclusions

The results showed that the novel online tool BarcodeGenerator (http://bargene.bi.up.ac.za/) is an efficient approach for generating barcode sequences from a set of complete genomes provided by users. Another program, Barcoder 2.0 is available from the same resource to enable an efficient and practical use of metabarcodes for visualization of the distribution of organisms of interest in environmental and clinical samples.

A reservoir of 'historical' antibiotic resistance genes in remote pristine Antarctic soils

BACKGROUND

Soil bacteria naturally produce antibiotics as a competitive mechanism, with a concomitant evolution, and exchange by horizontal gene transfer, of a range of antibiotic resistance mechanisms. Surveys of bacterial resistance elements in edaphic systems have originated primarily from human-impacted environments, with relatively little information from remote and pristine environments, where the resistome may comprise the ancestral gene diversity.

METHODS

We used shotgun metagenomics to assess antibiotic resistance gene (ARG) distribution in 17 pristine and remote Antarctic surface soils within the undisturbed Mackay Glacier region. We also interrogated the phylogenetic placement of ARGs compared to environmental ARG sequences and tested for the presence of horizontal gene transfer elements flanking ARGs.

RESULTS

In total, 177 naturally occurring ARGs were identified, most of which encoded single or multi-drug efflux pumps. Resistance mechanisms for the inactivation of aminoglycosides, chloramphenicol and β-lactam antibiotics were also common. Gram-negative bacteria harboured most ARGs (71%), with fewer genes from Gram-positive Actinobacteria and Bacilli (Firmicutes) (9%), reflecting the taxonomic composition of the soils. Strikingly, the abundance of ARGs per sample had a strong, negative correlation with species richness (r = - 0.49, P < 0.05). This result, coupled with a lack of mobile genetic elements flanking ARGs, suggests that these genes are ancient acquisitions of horizontal transfer events.

CONCLUSIONS

ARGs in these remote and uncontaminated soils most likely represent functional efficient historical genes that have since been vertically inherited over generations. The historical ARGs in these pristine environments carry a strong phylogenetic signal and form a monophyletic group relative to ARGs from other similar environments.

Draft genome sequence of Thermoactinomyces sp. strain AS95 isolated from a Sebkha in Thamelaht, Algeria

The members of the genus Thermoactinomyces are known for their protein degradative capacities. Thermoactinomyces sp. strain AS95 is a Gram-positive filamentous bacterium, isolated from moderately saline water in the Thamelaht region of Algeria. This isolate is a thermophilic aerobic bacterium with the capacity to produce extracellular proteolytic enzymes. This strain exhibits up to 99 % similarity with members of the genus Thermoactinomyces, based on 16S rRNA gene sequence similarity. Here we report on the phenotypic features of Thermoactinomyces sp. strain AS95 together with the draft genome sequence and its annotation. The genome of this strain is 2,558,690 bp in length (one chromosome, but no plasmid) with an average G + C content of 47.95 %, and contains 2550 protein-coding and 60 RNA genes together with 64 ORFs annotated as proteases.

The Geobacillus Pan-Genome: Implications for the Evolution of the Genus

The genus Geobacillus is comprised of a diverse group of spore-forming Gram-positive thermophilic bacterial species and is well known for both its ecological diversity and as a source of novel thermostable enzymes. Although the mechanisms underlying the thermophilicity of the organism and the thermostability of its macromolecules are reasonably well understood, relatively little is known of the evolutionary mechanisms, which underlie the structural and functional properties of members of this genus. In this study, we have compared 29 Geobacillus genomes, with a specific focus on the elements, which comprise the conserved core and flexible genomes. Based on comparisons of conserved core and flexible genomes, we present evidence of habitat delineation with specific Geobacillus genomes linked to specific niches. Our analysis revealed that Geobacillus and Anoxybacillus share a high proportion of genes. Moreover, the results strongly suggest that horizontal gene transfer is a major factor deriving the evolution of Geobacillus from Bacillus, with genetic contributions from other phylogenetically distant taxa.

Pre_GI: a global map of ontological links between horizontally transferred genomic islands in bacterial and archaeal genomes

The Predicted Genomic Islands database (Pre_GI) is a comprehensive repository of prokaryotic genomic islands (islands, GIs) freely accessible at http://pregi.bi.up.ac.za/index.php. Pre_GI, Version 2015, catalogues 26 744 islands identified in 2407 bacterial/archaeal chromosomes and plasmids. It provides an easy-to-use interface which allows users the ability to query against the database with a variety of fields, parameters and associations. Pre_GI is constructed to be a web-resource for the analysis of ontological roads between islands and cartographic analysis of the global fluxes of mobile genetic elements through bacterial and archaeal taxonomic borders. Comparison of newly identified islands against Pre_GI presents an alternative avenue to identify their ontology, origin and relative time of acquisition. Pre_GI aims to aid research on horizontal transfer events and materials through providing data and tools for holistic investigation of migration of genes through ecological niches and taxonomic boundaries.

Mainstreams of horizontal gene exchange in enterobacteria: consideration of the outbreak of enterohemorrhagic E. coli O104:H4 in Germany in 2011

BACKGROUND

Escherichia coli O104:H4 caused a severe outbreak in Europe in 2011. The strain TY-2482 sequenced from this outbreak allowed the discovery of its closest relatives but failed to resolve ways in which it originated and evolved. On account of the previous statement, may we expect similar upcoming outbreaks to occur recurrently or spontaneously in the future? The inability to answer these questions shows limitations of the current comparative and evolutionary genomics methods.

PRINCIPAL FINDINGS

The study revealed oscillations of gene exchange in enterobacteria, which originated from marine γ-Proteobacteria. These mobile genetic elements have become recombination hotspots and effective 'vehicles' ensuring a wide distribution of successful combinations of fitness and virulence genes among enterobacteria. Two remarkable peculiarities of the strain TY-2482 and its relatives were observed: i) retaining the genetic primitiveness by these strains as they somehow avoided the main fluxes of horizontal gene transfer which effectively penetrated other enetrobacteria; ii) acquisition of antibiotic resistance genes in a plasmid genomic island of β-Proteobacteria origin which ontologically is unrelated to the predominant genomic islands of enterobacteria.

CONCLUSIONS

Oscillations of horizontal gene exchange activity were reported which result from a counterbalance between the acquired resistance of bacteria towards existing mobile vectors and the generation of new vectors in the environmental microflora. We hypothesized that TY-2482 may originate from a genetically primitive lineage of E. coli that has evolved in confined geographical areas and brought by human migration or cattle trade onto an intersection of several independent streams of horizontal gene exchange. Development of a system for monitoring the new and most active gene exchange events was proposed.