Clonal and plasmidic dissemination of critical antimicrobial resistance genes through clinically relevant ExPEC and APEC-like lineages (ST) in the dairy cattle population of Québec, Canada.

Antimicrobial resistance can be effectively limited by improving the judicious use of antimicrobials in food production. However, its effect on the spread of AMR genes in animal populations is not well described. In the province of Québec, Canada, a new legislation implemented in 2019 has led to an unprecedented reduction in the use of critical antimicrobials in dairy production. We aimed to investigate the potential link between ESBL/AmpC E. coli isolated before and after legislation and to determine the presence of plasmids carrying genes responsible for critical AMR. We collected fecal samples from calves, cows, and manure pit from 87 Québec dairy farms approximately 2 years before and 2 years after the legislation came into effect. The whole genomes of 183 presumptive ESBL/AmpC E. coli isolated after cefotaxime enrichment were sequenced. Their phylogenetic characteristics (MLST, serogroup, cgMLST) and the presence of virulence and resistance genes and replicons were examined. A maximum likelihood phylogenetic tree was constructed based on single nucleotide polymorphism (SNPs). We identified 10 clonal lineages (same cgMLST) and 7 clones (SNPs ≤ 52). Isolates belonging to these clones could be found on different farms before and after the legislation, strongly suggesting a clonal spread of AMR genes in the population during this 4-year period. All isolates were multidrug resistant (MDR), with clone 2 being notable for the presence of macrolide, fluoroquinolone, and third-generation cephalosporin resistance genes. We also identified clinically relevant ExPEC (ST10) and APEC-like lineages (ST117, ST58, ST88) associated with the presence of ExPEC and APEC virulence genes, respectively. Our data also suggests the presence of one epidemic plasmid belonging to the IncY incompatibility group and carrying qnrs1 and blaCTX-M-15. We demonstrated that AMR genes spread through farms and can persist over a 4-year period in the dairy cattle population through both plasmids and E. coli clones, despite the restriction of critical antimicrobial use. MDR ExPEC and APEC-like STs are present in the normal microbiota of cattle (more frequently in calves). These data increase our knowledge on gene dissemination dynamics and highlight the fact that biosecurity measures should be enhanced in this industry to limit such dissemination.

Impact of Stagnation on the Diversity of Cyanobacteria in Drinking Water Treatment Plant Sludge.

Health-related concerns about cyanobacteria-laden sludge of drinking water treatment plants (DWTPs) have been raised in the past few years. Microscopic taxonomy, shotgun metagenomic sequencing, and microcystin (MC) measurement were applied to study the fate of cyanobacteria and cyanotoxins after controlled sludge storage (stagnation) in the dark in a full-scale drinking water treatment plant within 7 to 38 days. For four out of eight dates, cyanobacterial cell growth was observed by total taxonomic cell counts during sludge stagnation. The highest observed cell growth was 96% after 16 days of stagnation. Cell growth was dominated by potential MC producers such as Microcystis, Aphanocapsa, Chroococcus, and Dolichospermum. Shotgun metagenomic sequencing unveiled that stagnation stress shifts the cyanobacterial communities from the stress-sensitive Nostocales (e.g., Dolichospermum) order towards less compromised orders and potential MC producers such as Chroococcales (e.g., Microcystis) and Synechococcales (e.g., Synechococcus). The relative increase of cyanotoxin producers presents a health challenge when the supernatant of the stored sludge is recycled to the head of the DWTP or discharged into the source. These findings emphasize the importance of a strategy to manage cyanobacteria-laden sludge and suggest practical approaches should be adopted to control health/environmental impacts of cyanobacteria and cyanotoxins in sludge.

Shotgun Metagenomic Sequencing to Assess Cyanobacterial Community Composition following Coagulation of Cyanobacterial Blooms.

The excessive proliferation of cyanobacteria in surface waters is a widespread problem worldwide, leading to the contamination of drinking water sources. Short- and long-term solutions for managing cyanobacterial blooms are needed for drinking water supplies. The goal of this research was to investigate the cyanobacteria community composition using shotgun metagenomics in a short term, in situ mesocosm experiment of two lakes following their coagulation with ferric sulfate (Fe2(SO4)3) as an option for source water treatment. Among the nutrient paramenters, dissolved nitrogen was related to Microcystis in both Missisquoi Bay and Petit Lac St. François, while the presence of Synechococcus was related to total nitrogen, dissolved nitrogen, dissolved organic carbon, and dissolved phosphorus. Results from the shotgun metagenomic sequencing showed that Dolichospermum and Microcystis were the dominant genera in all of the mesocosms in the beginning of the sampling period in Missisquoi Bay and Petit Lac St. François, respectively. Potentially toxigenic genera such as Microcystis were correlated with intracellular microcystin concentrations. A principal component analysis showed that there was a change of the cyanobacterial composition at the genus level in the mesocosms after two days, which varied across the studied sites and sampling time. The cyanobacterial community richness and diversity did not change significantly after its coagulation by Fe2(SO4)3 in all of the mesocosms at either site. The use of Fe2(SO4)3 for an onsite source water treatment should consider its impact on cyanobacterial community structure and the reduction of toxin concentrations.

Soil Depth Significantly Shifted Microbial Community Structures and Functions in a Semiarid Prairie Agroecosystem.

The North American Great Plains cover a large area of the Nearctic ecozone, and an important part of this biome is semiarid. The sustainable intensification of agriculture that is necessary to produce food for an ever-increasing world population requires knowledge of the taxonomic and functional structure of the soil microbial community. In this study, we investigated the influence of soil depth on the composition and functions of the microbial communities hosted in agricultural soils of a semiarid agroecosystem, using metagenomic profiling, and compared them to changes in soil chemical and physical properties. Shotgun sequencing was used to determine the composition and functions of the soil microbial community of 45 soil samples from three soil depths (0-15 cm, 15-30 cm, and 30-60 cm) under different agricultural land use types (native prairie, seeded prairie, and cropland) in southwest Saskatchewan. Analysis of community composition revealed the declining abundance of phyla Verrucomicrobia, Bacteroidetes, Chlorophyta, Bacillariophyta, and Acidobacteria with soil depth, whereas the abundance of phyla Ascomycota, Nitrospirae, Planctomycetes, and Cyanobacteria increased with soil depth. Soil functional genes related to nucleosides and nucleotides, phosphorus (P) metabolism, cell division and cell cycle, amino acids and derivatives, membrane transport, and fatty acids were particularly abundant at 30-60 cm. In contrast, functional genes related to DNA and RNA metabolism, metabolism of nitrogen, sulfur and carbohydrates, and stress response were more abundant in the top soil depth. The RDA analysis of functional genes and soil physico-chemical properties revealed a positive correlation between phages and soil organic P concentrations. In the rooting zone of this semiarid agroecosystem, soil microbes express variable structural patterns of taxonomic and functional diversity at different soil depths. This study shows that the soil microbial community is structured by soil depth and physicochemical properties, with the middle soil depth being an intermediate transition zone with a higher taxonomic diversity. Our results suggest the co-existence of various microbial phyla adapted to upper and lower soil depths in an intermediate-depth transition zone.

A Glyphosate-Based Herbicide Cross-Selects for Antibiotic Resistance Genes in Bacterioplankton Communities.

Agrochemicals often contaminate freshwater bodies, affecting microbial communities that underlie aquatic food webs. For example, the herbicide glyphosate has the potential to indirectly select for antibiotic-resistant bacteria. Such cross-selection could occur if the same genes (encoding efflux pumps, for example) confer resistance to both glyphosate and antibiotics. To test for cross-resistance in natural aquatic bacterial communities, we added a glyphosate-based herbicide (GBH) to 1,000-liter mesocosms filled with water from a pristine lake. Over 57 days, we tracked changes in bacterial communities with shotgun metagenomic sequencing and annotated metagenome-assembled genomes (MAGs) for the presence of known antibiotic resistance genes (ARGs), plasmids, and resistance mutations in the enzyme targeted by glyphosate (enolpyruvyl-shikimate-3-phosphate synthase; EPSPS). We found that high doses of GBH significantly increased ARG frequency and selected for multidrug efflux pumps in particular. The relative abundance of MAGs after a high dose of GBH was predictable based on the number of ARGs in their genomes (17% of variation explained) and, to a lesser extent, by resistance mutations in EPSPS. Together, these results indicate that GBHs can cross-select for antibiotic resistance in natural freshwater bacteria. IMPORTANCE Glyphosate-based herbicides (GBHs) such as Roundup formulations may have the unintended consequence of selecting for antibiotic resistance genes (ARGs), as demonstrated in previous experiments. However, the effects of GBHs on ARGs remain unknown in natural aquatic communities, which are often contaminated with pesticides from agricultural runoff. Moreover, the resistance provided by ARGs compared to canonical mutations in the glyphosate target enzyme, EPSPS, remains unclear. Here, we performed a freshwater mesocosm experiment showing that a GBH strongly selects for ARGs, particularly multidrug efflux pumps. These selective effects were evident after just a few days, and the ability of bacteria to survive and thrive after GBH stress was predictable by the number of ARGs in their genomes and, to a lesser extent, by mutations in EPSPS. Intensive GBH application may therefore have the unintended consequence of selecting for ARGs in natural freshwater communities.

Single-colony sequencing reveals microbe-by-microbiome phylosymbiosis between the cyanobacterium Microcystis and its associated bacteria.

BACKGROUND: Cyanobacteria from the genus Microcystis can form large mucilaginous colonies with attached heterotrophic bacteria-their microbiome. However, the nature of the relationship between Microcystis and its microbiome remains unclear. Is it a long-term, evolutionarily stable association? Which partners benefit? Here we report the genomic diversity of 109 individual Microcystis colonies-including cyanobacteria and associated bacterial genomes-isolated in situ and without culture from Lake Champlain, Canada and Pampulha Reservoir, Brazil. RESULTS: We identified 14 distinct Microcystis genotypes from Canada, of which only two have been previously reported, and four genotypes specific to Brazil. Microcystis genetic diversity was much greater between than within colonies, consistent with colony growth by clonal expansion rather than aggregation of Microcystis cells. We also identified 72 bacterial species in the microbiome. Each Microcystis genotype had a distinct microbiome composition, and more closely related genotypes had more similar microbiomes. This pattern of phylosymbiosis could be explained by co-phylogeny in only two out of the nine most prevalent associated bacterial genera, Roseomonas and Rhodobacter. These phylogenetically associated genera could enrich the metabolic repertoire of Microcystis, for example by encoding the biosynthesis of complementary carotenoid molecules. In contrast, other colony-associated bacteria showed weaker signals of co-phylogeny, but stronger evidence of horizontal gene transfer with Microcystis. These observations suggest that acquired genes are more likely to be retained in both partners (Microcystis and members of its microbiome) when they are loosely associated, whereas one gene copy is sufficient when the association is physically tight and evolutionarily long-lasting. CONCLUSIONS: We have introduced a method for culture-free isolation of single colonies from nature followed by metagenomic sequencing, which could be applied to other types of microbes. Together, our results expand the known genetic diversity of both Microcystis and its microbiome in natural settings, and support their long-term, specific, and potentially beneficial associations. Video Abstract.

Can Cyanobacterial Diversity in the Source Predict the Diversity in Sludge and the Risk of Toxin Release in a Drinking Water Treatment Plant?

Conventional processes (coagulation, flocculation, sedimentation, and filtration) are widely used in drinking water treatment plants and are considered a good treatment strategy to eliminate cyanobacterial cells and cell-bound cyanotoxins. The diversity of cyanobacteria was investigated using taxonomic cell counts and shotgun metagenomics over two seasons in a drinking water treatment plant before, during, and after the bloom. Changes in the community structure over time at the phylum, genus, and species levels were monitored in samples retrieved from raw water (RW), sludge in the holding tank (ST), and sludge supernatant (SST). Aphanothece clathrata brevis, Microcystis aeruginosa, Dolichospermum spiroides , and Chroococcus minimus were predominant species detected in RW by taxonomic cell counts. Shotgun metagenomics revealed that Proteobacteria was the predominant phylum in RW before and after the cyanobacterial bloom. Taxonomic cell counts and shotgun metagenomic showed that the Dolichospermum bloom occurred inside the plant. Cyanobacteria and Bacteroidetes were the major bacterial phyla during the bloom. Shotgun metagenomics also showed that Synechococcus, Microcystis , and Dolichospermum were the predominant detected cyanobacterial genera in the samples. Conventional treatment removed more than 92% of cyanobacterial cells but led to cell accumulation in the sludge up to 31 times more than in the RW influx. Coagulation/sedimentation selectively removed more than 96% of Microcystis and Dolichospermum. Cyanobacterial community in the sludge varied from raw water to sludge during sludge storage (1-13 days). This variation was due to the selective removal of coagulation/sedimentation as well as the accumulation of captured cells over the period of storage time. However, the prediction of the cyanobacterial community composition in the SST remained a challenge. Among nutrient parameters, orthophosphate availability was related to community profile in RW samples, whereas communities in ST were influenced by total nitrogen, Kjeldahl nitrogen (N- Kjeldahl), total and particulate phosphorous, and total organic carbon (TOC). No trend was observed on the impact of nutrients on SST communities. This study profiled new health-related, environmental, and technical challenges for the production of drinking water due to the complex fate of cyanobacteria in cyanobacteria-laden sludge and supernatant.

Diversity Assessment of Toxic Cyanobacterial Blooms during Oxidation.

Fresh-water sources of drinking water are experiencing toxic cyanobacterial blooms more frequently. Chemical oxidation is a common approach to treat cyanobacteria and their toxins. This study systematically investigates the bacterial/cyanobacterial community following chemical oxidation (Cl2, KMnO4, O3, H2O2) using high throughput sequencing. Raw water results from high throughput sequencing show that Proteobacteria, Actinobacteria, Cyanobacteria and Bacteroidetes were the most abundant phyla. Dolichospermum, Synechococcus, Microcystis and Nostoc were the most dominant genera. In terms of species, Dolichospermum sp.90 and Microcystis aeruginosa were the most abundant species at the beginning and end of the sampling, respectively. A comparison between the results of high throughput sequencing and taxonomic cell counts highlighted the robustness of high throughput sequencing to thoroughly reveal a wide diversity of bacterial and cyanobacterial communities. Principal component analysis of the oxidation samples results showed a progressive shift in the composition of bacterial/cyanobacterial communities following soft-chlorination with increasing common exposure units (CTs) (0-3.8 mg·min/L). Close cyanobacterial community composition (Dolichospermum dominant genus) was observed following low chlorine and mid-KMnO4 (287.7 mg·min/L) exposure. Our results showed that some toxin producing species may persist after oxidation whether they were dominant species or not. Relative persistence of Dolichospermum sp.90 was observed following soft-chlorination (0.2-0.6 mg/L) and permanganate (5 mg/L) oxidation with increasing oxidant exposure. Pre-oxidation using H2O2 (10 mg/L and one day contact time) caused a clear decrease in the relative abundance of all the taxa and some species including the toxin producing taxa. These observations suggest selectivity of H2O2 to provide an efficient barrier against toxin producing cyanobacteria entering a water treatment plant.

Two cases of type-a Haemophilus influenzae meningitis within the same week in the same hospital are phylogenetically unrelated but recently exchanged capsule genes.

Haemophilus influenzae causes common and sometimes severe adult and pediatric disease including chronic obstructive respiratory disease, otitis media and infections of the central nervous system. Serotype b strains, with a b-type capsule, have been the historical cause of invasive disease, and the introduction of a serotype b-specific vaccine has led to their decline. However, unencapsulated or non-b-type H. influenzae infections are not prevented by the vaccine and appear to be increasing in frequency. Here we report two pediatric cases of severe central nervous system H. influenzae infection presenting to the same hospital in San Diego, California during the same week in January 2016. Due to good vaccine coverage in this part of the world, H. influenzae cases are normally rare and seeing two cases in the same week was unexpected. We thus suspected a recent transmission chain, and possible local outbreak. To test this hypothesis, we isolated and sequenced whole genomes from each patient and placed them in a phylogenetic tree spanning the known diversity of H. influenzae. Surprisingly, we found that the two isolates (SD2016_1 and SD2016_2) belonged to distantly related lineages, suggesting two independent transmission events and ruling out a local outbreak. Despite being distantly related, the two isolates belong to two different lineages that have exchanged capsule loci in the recent past. Therefore, as in other bacterial pathogens, capsule switching by horizontal gene transfer may be an important evolutionary mechanism of vaccine evasion in H. influenzae.

Coherence of Microcystis species revealed through population genomics.

Microcystis is a genus of freshwater cyanobacteria, which causes harmful blooms in ecosystems worldwide. Some Microcystis strains produce harmful toxins such as microcystin, impacting drinking water quality. Microcystis colony morphology, rather than genetic similarity, is often used to classify Microcystis into morphospecies. Yet colony morphology is a plastic trait, which can change depending on environmental and laboratory culture conditions, and is thus an inadequate criterion for species delineation. Furthermore, Microcystis populations are thought to disperse globally and constitute a homogeneous gene pool. However, this assertion is based on relatively incomplete characterization of Microcystis genomic diversity. To better understand these issues, we performed a population genomic analysis of 33 newly sequenced genomes mainly from Canada and Brazil. We identified 17 Microcystis clusters of genomic similarity, five of which correspond to monophyletic clades containing at least three newly sequenced genomes. Four out of these five clades match to named morphospecies. Notably, M. aeruginosa is paraphyletic, distributed across 12 genomic clusters, suggesting it is not a coherent species. A few clades of closely related isolates are specific to a unique geographic location, suggesting biogeographic structure over relatively short evolutionary time scales. Higher homologous recombination rates within than between clades further suggest that monophyletic groups might adhere to a Biological Species-like concept, in which barriers to gene flow maintain species distinctness. However, certain genes-including some involved in microcystin and micropeptin biosynthesis-are recombined between monophyletic groups in the same geographic location, suggesting local adaptation.

Characterising and predicting cyanobacterial blooms in an 8-year amplicon sequencing time course.

Cyanobacterial blooms occur in lakes worldwide, producing toxins that pose a serious public health threat. Eutrophication caused by human activities and warmer temperatures both contribute to blooms, but it is still difficult to predict precisely when and where blooms will occur. One reason that prediction is so difficult is that blooms can be caused by different species or genera of cyanobacteria, which may interact with other bacteria and respond to a variety of environmental cues. Here we used a deep 16S amplicon sequencing approach to profile the bacterial community in eutrophic Lake Champlain over time, to characterise the composition and repeatability of cyanobacterial blooms, and to determine the potential for blooms to be predicted based on time course sequence data. Our analysis, based on 135 samples between 2006 and 2013, spans multiple bloom events. We found that bloom events significantly alter the bacterial community without reducing overall diversity, suggesting that a distinct microbial community-including non-cyanobacteria-prospers during the bloom. We also observed that the community changes cyclically over the course of a year, with a repeatable pattern from year to year. This suggests that, in principle, bloom events are predictable. We used probabilistic assemblages of OTUs to characterise the bloom-associated community, and to classify samples into bloom or non-bloom categories, achieving up to 92% classification accuracy (86% after excluding cyanobacterial sequences). Finally, using symbolic regression, we were able to predict the start date of a bloom with 78-92% accuracy (depending on the data used for model training), and found that sequence data was a better predictor than environmental variables.

Mitochondrial Recombination and Introgression during Speciation by Hybridization.

Genome recombination is a major source of genotypic diversity and contributes to adaptation and speciation following interspecies hybridization. The contribution of recombination in these processes has been thought to be largely limited to the nuclear genome because organelles are mostly uniparentally inherited in animals and plants, which prevents recombination. Unicellular eukaryotes such as budding yeasts do, however, transmit mitochondria biparentally, suggesting that during hybridization, both parents could provide alleles that contribute to mitochondrial functions such as respiration and metabolism in hybrid populations or hybrid species. We examined the dynamics of mitochondrial genome transmission and evolution during speciation by hybridization in the natural budding yeast Saccharomyces paradoxus. Using population-scale mitochondrial genome sequencing in two endemic North American incipient species SpB and SpC and their hybrid species SpC*, we found that both parental species contributed to the hybrid mitochondrial genome through recombination. We support our findings by showing that mitochondrial recombination between parental types is frequent in experimental crosses that recreate the early step of this speciation event. In these artificial hybrids, we observed that mitochondrial genome recombination enhances phenotypic variation among diploid hybrids, suggesting that it could play a role in the phenotypic differentiation of hybrid species. Like the nuclear genome, the mitochondrial genome can, therefore, also play a role in hybrid speciation.

Plant Community and Nitrogen Deposition as Drivers of Alpha and Beta Diversities of Prokaryotes in Reconstructed Oil Sand Soils and Natural Boreal Forest Soils.

The Athabasca oil sand deposit is one of the largest single oil deposits in the world. Following surface mining, companies are required to restore soil-like profiles that can support the previous land capabilities. The objective of this study was to assess whether the soil prokaryotic alpha diversity (α-diversity) and ß-diversity in oil sand soils reconstructed 20 to 30 years previously and planted to one of three vegetation types (coniferous or deciduous trees and grassland) were similar to those found in natural boreal forest soils subject to wildfire disturbance. Prokaryotic α-diversity and ß-diversity were assessed using massively parallel sequencing of 16S rRNA genes. The ß-diversity, but not the α-diversity, differed between reconstructed and natural soils. Bacteria associated with an oligotrophic lifestyle were more abundant in natural forest soils, whereas bacteria associated with a copiotrophic lifestyle were more abundant in reconstructed soils. Ammonia-oxidizing archaea were most abundant in reconstructed soils planted with grasses. Plant species were the main factor influencing α-diversity in natural and in reconstructed soils. Nitrogen deposition, pH, and plant species were the main factors influencing the ß-diversity of the prokaryotic communities in natural and reconstructed soils. The results highlight the importance of nitrogen deposition and aboveground-belowground relationships in shaping soil microbial communities in natural and reconstructed soils.IMPORTANCE Covering over 800 km2, land disturbed by the exploitation of the oil sands in Canada has to be restored. Here, we take advantage of the proximity between these reconstructed ecosystems and the boreal forest surrounding the oil sand mining area to study soil microbial community structure and processes in both natural and nonnatural environments. By identifying key characteristics shaping the structure of soil microbial communities, this study improved our understanding of how vegetation, soil characteristics and microbial communities interact and drive soil functions.

Vibrio cholerae genomic diversity within and between patients.

Cholera is a severe, water-borne diarrhoeal disease caused by toxin-producing strains of the bacterium Vibrio cholerae. Comparative genomics has revealed 'waves' of cholera transmission and evolution, in which clones are successively replaced over decades and centuries. However, the extent of V. cholerae genetic diversity within an epidemic or even within an individual patient is poorly understood. Here, we characterized V. cholerae genomic diversity at a micro-epidemiological level within and between individual patients from Bangladesh and Haiti. To capture within-patient diversity, we isolated multiple (8 to 20) V. cholerae colonies from each of eight patients, sequenced their genomes and identified point mutations and gene gain/loss events. We found limited but detectable diversity at the level of point mutations within hosts (zero to three single nucleotide variants within each patient), and comparatively higher gene content variation within hosts (at least one gain/loss event per patient, and up to 103 events in one patient). Much of the gene content variation appeared to be due to gain and loss of phage and plasmids within the V. cholerae population, with occasional exchanges between V. cholerae and other members of the gut microbiota. We also show that certain intra-host variants have phenotypic consequences. For example, the acquisition of a Bacteroides plasmid and non-synonymous mutations in a sensor histidine kinase gene both reduced biofilm formation, an important trait for environmental survival. Together, our results show that V. cholerae is measurably evolving within patients, with possible implications for disease outcomes and transmission dynamics.

Mating type gene homologues and putative sex pheromone-sensing pathway in arbuscular mycorrhizal fungi, a presumably asexual plant root symbiont.

The fungal kingdom displays a fascinating diversity of sex-determination systems. Recent advances in genomics provide insights into the molecular mechanisms of sex, mating type determination, and evolution of sexual reproduction in many fungal species in both ancient and modern phylogenetic lineages. All major fungal groups have evolved sexual differentiation and recombination pathways. However, sexuality is unknown in arbuscular mycorrhizal fungi (AMF) of the phylum Glomeromycota, an ecologically vital group of obligate plant root symbionts. AMF are commonly considered an ancient asexual lineage dating back to the Ordovician, approximately 460 M years ago. In this study, we used genomic and transcriptomic surveys of several AMF species to demonstrate the presence of conserved putative sex pheromone-sensing mitogen-activated protein (MAP) kinases, comparable to those described in Ascomycota and Basidiomycota. We also find genes for high mobility group (HMG) transcription factors, homologous to SexM and SexP genes in the Mucorales. The SexM genes show a remarkable sequence diversity among multiple copies in the genome, while only a single SexP sequence was detected in some isolates of Rhizophagus irregularis. In the Mucorales and Microsporidia, the sexM gene is flanked by genes for a triosephosphate transporter (TPT) and a RNA helicase, but we find no evidence for synteny in the vicinity of the Sex locus in AMF. Nonetheless, our results, together with previous observations on meiotic machinery, suggest that AMF could undergo a complete sexual reproduction cycle.

Atractaspis aterrima toxins: the first insight into the molecular evolution of venom in side-stabbers.

Although snake venoms have been the subject of intense research, primarily because of their tremendous potential as a bioresource for design and development of therapeutic compounds, some specific groups of snakes, such as the genus Atractaspis, have been completely neglected. To date only limited number of toxins, such as sarafotoxins have been well characterized from this lineage. In order to investigate the molecular diversity of venom from Atractaspis aterrima-the slender burrowing asp, we utilized a high-throughput transcriptomic approach completed with an original bioinformatics analysis pipeline. Surprisingly, we found that Sarafotoxins do not constitute the major ingredient of the transcriptomic cocktail; rather a large number of previously well-characterized snake venom-components were identified. Notably, we recovered a large diversity of three-finger toxins (3FTxs), which were found to have evolved under the significant influence of positive selection. From the normalized and non-normalized transcriptome libraries, we were able to evaluate the relative abundance of the different toxin groups, uncover rare transcripts, and gain new insight into the transcriptomic machinery. In addition to previously characterized toxin families, we were able to detect numerous highly-transcribed compounds that possess all the key features of venom-components and may constitute new classes of toxins.

Mitochondrial genome rearrangements in glomus species triggered by homologous recombination between distinct mtDNA haplotypes.

Comparative mitochondrial genomics of arbuscular mycorrhizal fungi (AMF) provide new avenues to overcome long-lasting obstacles that have hampered studies aimed at understanding the community structure, diversity, and evolution of these multinucleated and genetically polymorphic organisms.AMF mitochondrial (mt) genomes are homogeneous within isolates, and their intergenic regions harbor numerous mobile elements that have rapidly diverged, including homing endonuclease genes, small inverted repeats, and plasmid-related DNA polymerase genes (dpo), making them suitable targets for the development of reliable strain-specific markers. However, these elements may also lead to genome rearrangements through homologous recombination, although this has never previously been reported in this group of obligate symbiotic fungi. To investigate whether such rearrangements are present and caused by mobile elements in AMF, the mitochondrial genomes from two Glomeraceae members (i.e., Glomus cerebriforme and Glomus sp.) with substantial mtDNA synteny divergence,were sequenced and compared with available glomeromycotan mitochondrial genomes. We used an extensive nucleotide/protein similarity network-based approach to investigated podiversity in AMF as well as in other organisms for which sequences are publicly available. We provide strong evidence of dpo-induced inter-haplotype recombination, leading to a reshuffled mitochondrial genome in Glomus sp. These findings raise questions as to whether AMF single spore cultivations artificially underestimate mtDNA genetic diversity.We assessed potential dpo dispersal mechanisms in AMF and inferred a robust phylogenetic relationship with plant mitochondrial plasmids. Along with other indirect evidence, our analyses indicate that members of the Glomeromycota phylum are potential donors of mitochondrial plasmids to plants.

Recruitment of glycosyl hydrolase proteins in a cone snail venomous arsenal: further insights into biomolecular features of Conus venoms.

Cone snail venoms are considered an untapped reservoir of extremely diverse peptides, named conopeptides, displaying a wide array of pharmacological activities. We report here for the first time, the presence of high molecular weight compounds that participate in the envenomation cocktail used by these marine snails. Using a combination of proteomic and transcriptomic approaches, we identified glycosyl hydrolase proteins, of the hyaluronidase type (Hyal), from the dissected and injectable venoms ("injectable venom" stands for the venom variety obtained by milking of the snails. This is in contrast to the "dissected venom", which was obtained from dissected snails by extraction of the venom glands) of a fish-hunting cone snail, Conus consors (Pionoconus clade). The major Hyal isoform, Conohyal-Cn1, is expressed as a mixture of numerous glycosylated proteins in the 50 kDa molecular mass range, as observed in 2D gel and mass spectrometry analyses. Further proteomic analysis and venom duct mRNA sequencing allowed full sequence determination. Additionally, unambiguous segment location of at least three glycosylation sites could be determined, with glycans corresponding to multiple hexose (Hex) and N-acetylhexosamine (HexNAc) moieties. With respect to other known Hyals, Conohyal-Cn1 clearly belongs to the hydrolase-type of Hyals, with strictly conserved consensus catalytic donor and positioning residues. Potent biological activity of the native Conohyals could be confirmed in degrading hyaluronic acid. A similar Hyal sequence was also found in the venom duct transcriptome of C. adamsonii (Textilia clade), implying a possible widespread recruitment of this enzyme family in fish-hunting cone snail venoms. These results provide the first detailed Hyal sequence characterized from a cone snail venom, and to a larger extent in the Mollusca phylum, thus extending our knowledge on this protein family and its evolutionary selection in marine snail venoms.

High-resolution picture of a venom gland transcriptome: case study with the marine snail Conus consors.

Although cone snail venoms have been intensively investigated in the past few decades, little is known about the whole conopeptide and protein content in venom ducts, especially at the transcriptomic level. If most of the previous studies focusing on a limited number of sequences have contributed to a better understanding of conopeptide superfamilies, they did not give access to a complete panorama of a whole venom duct. Additionally, rare transcripts were usually not identified due to sampling effect. This work presents the data and analysis of a large number of sequences obtained from high throughput 454 sequencing technology using venom ducts of Conus consors, an Indo-Pacific living piscivorous cone snail. A total of 213,561 Expressed Sequence Tags (ESTs) with an average read length of 218 base pairs (bp) have been obtained. These reads were assembled into 65,536 contiguous DNA sequences (contigs) then into 5039 clusters. The data revealed 11 conopeptide superfamilies representing a total of 53 new isoforms (full length or nearly full-length sequences). Considerable isoform diversity and major differences in transcription level could be noted between superfamilies. A, O and M superfamilies are the most diverse. The A family isoforms account for more than 70% of the conopeptide cocktail (considering all ESTs before clustering step). In addition to traditional superfamilies and families, minor transcripts including both cysteine free and cysteine-rich peptides could be detected, some of them figuring new clades of conopeptides. Finally, several sets of transcripts corresponding to proteins commonly recruited in venom function could be identified for the first time in cone snail venom duct. This work provides one of the first large-scale EST project for a cone snail venom duct using next-generation sequencing, allowing a detailed overview of the venom duct transcripts. This leads to an expanded definition of the overall cone snail venom duct transcriptomic activity, which goes beyond the cysteine-rich conopeptides. For instance, this study enabled to detect proteins involved in common post-translational maturation and folding, and to reveal compounds classically involved in hemolysis and mechanical penetration of the venom into the prey. Further comparison with proteomic and genomic data will lead to a better understanding of conopeptides diversity and the underlying mechanisms involved in conopeptide evolution.