Identification and functional analysis of Cochliomyia hominivorax U6 gene promoters.

The New World screwworm, Cochliomyia hominivorax, is an obligate parasite, which is a major pest of livestock. While the sterile insect technique was used very successfully to eradicate C. hominivorax from North and Central America, more cost-effective genetic methods will likely be needed in South America. The recent development of CRISPR/Cas9-based genetic approaches, such as homing gene drive, could provide a very efficient means for the suppression of C. hominivorax populations. One component of a drive system is the guide RNA(s) driven by a U6 gene promoter. Here, we have developed an in vivo assay to evaluate the activity of the promoters from seven C. hominivorax U6 genes. Embryos from the related blowfly Lucilia cuprina were injected with plasmid DNA containing a U6-promoter-guide RNA construct and a source of Cas9, either protein or plasmid DNA. Activity was assessed by the number of site-specific mutations in the targeted gene in hatched larvae. One promoter, Chom U6_b, showed the highest activity. These U6 gene promoters could be used to build CRISPR/Cas9-based genetic systems for the control of C. hominivorax.

Whole genomic comparative analysis of Streptococcus pneumoniae serotype 1 isolates causing invasive and non-invasive infections among children under 5 years in Casablanca, Morocco.

BACKGROUND: Streptococcus pneumoniae serotype 1 remains a leading cause of invasive pneumococcal diseases, even in countries with PCV-10/PCV-13 vaccine implementation. The main objective of this study, which is part of the Pneumococcal African Genome project (PAGe), was to determine the phylogenetic relationships of serotype 1 isolates recovered from children patients in Casablanca (Morocco), compared to these from other African countries; and to investigate the contribution of accessory genes and recombination events to the genetic diversity of this serotype. RESULTS: The genome average size of the six-pneumococcus serotype 1 from Casablanca was 2,227,119 bp, and the average content of coding sequences was 2113, ranging from 2041 to 2161. Pangenome analysis of the 80 genomes used in this study revealed 1685 core genes and 1805 accessory genes. The phylogenetic tree based on core genes and the hierarchical bayesian clustering analysis revealed five sublineages with a phylogeographic structure by country. The Moroccan strains cluster in two different lineages, the five invasive strains clusters altogether in a divergent clade distantly related to the non-invasive strain, that cluster with all the serotype 1 genomes from Africa. CONCLUSIONS: The whole genome sequencing provides increased resolution analysis of the highly virulent serotype 1 in Casablanca, Morocco. Our results are concordant with previous works, showing that the phylogeography of S. pneumoniae serotype 1 is structured by country, and despite the small size (six isolates) of the Moroccan sample, our analysis shows the genetic cohesion of the Moroccan invasive isolates.

simurg: simulate bacterial pangenomes in R.

MOTIVATION: The pangenome concept describes genetic variability as the union of genes shared in a set of genomes and constitutes the current paradigm for comparative analysis of bacterial populations. However, there is a lack of tools to simulate pangenome variability and structure using defined evolutionary models. RESULTS: We developed simurg, an R package that allows to simulate bacterial pangenomes using different combinations of evolutionary constraints such as gene gain, gene loss and mutation rates. Our tool allows the straightforward and reproducible simulation of bacterial pangenomes using real sequence data, providing a valuable tool for benchmarking of pangenome software or comparing evolutionary hypotheses. AVAILABILITY AND IMPLEMENTATION: The simurg package is released under the GPL-3 license, and is freely available for download from GitHub (https://github.com/iferres/simurg). SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Intraspecific genetic variation matters when predicting seagrass distribution under climate change.

Seagrasses play a vital role in structuring coastal marine ecosystems, but their distributional range and genetic diversity have declined rapidly in recent decades. To improve conservation of seagrass species, it is important to predict how climate change may impact their ranges. Such predictions are typically made with correlative species distribution models (SDMs), which can estimate a species' potential distribution under present and future climatic scenarios given species' presence data and climatic predictor variables. However, these models are typically constructed with species-level data, and thus ignore intraspecific genetic variability, which can give rise to populations with adaptations to heterogeneous climatic conditions. Here, we explore the link between intraspecific adaptation and niche differentiation in Thalassia hemprichii, a seagrass broadly distributed in the tropical Indo-Pacific Ocean and a crucial provider of habitat for numerous marine species. By retrieving and re-analysing microsatellite data from previous studies, we delimited two distinct phylogeographical lineages within the nominal species and found an intermediate level of differentiation in their multidimensional environmental niches, suggesting the possibility for local adaptation. We then compared projections of the species' habitat suitability under climate change scenarios using species-level and lineage-level SDMs. In the Central Tropical Indo-Pacific region, models for both levels predicted considerable range contraction in the future, but the lineage-level models predicted more severe habitat loss. Importantly, the two modelling approaches predicted opposite patterns of habitat change in the Western Tropical Indo-Pacific region. Our results highlight the necessity of conserving distinct populations and genetic pools to avoid regional extinction due to climate change and have important implications for guiding future management of seagrasses.

Ibrutinib therapy downregulates AID enzyme and proliferative fractions in chronic lymphocytic leukemia.

Activation-induced cytidine deaminase (AID) initiates somatic hypermutation and class switch recombination of the immunoglobulin genes. As a trade-off for its physiological function, AID also contributes to tumor development through its mutagenic activity. In chronic lymphocytic leukemia (CLL), AID is overexpressed in the proliferative fractions (PFs) of the malignant B lymphocytes, and its anomalous expression has been associated with a clinical poor outcome. Recent preclinical data suggested that ibrutinib and idelalisib, 2 clinically approved kinase inhibitors, increase AID expression and genomic instability in normal and neoplastic B cells. These results raise concerns about a potential mutagenic risk in patients receiving long-term therapy. To corroborate these findings in the clinical setting, we analyzed AID expression and PFs in a CLL cohort before and during ibrutinib treatment. We found that ibrutinib decreases the CLL PFs and, interestingly, also reduces AID expression, which correlates with dampened AKT and Janus Kinase 1 signaling. Moreover, although ibrutinib increases AID expression in a CLL cell line, it is unable to do so in primary CLL samples. Our results uncover a differential response to ibrutinib between cell lines and the CLL clone and imply that ibrutinib could differ from idelalisib in their potential to induce AID in treated patients. Possible reasons for the discrepancy between preclinical and clinical findings, and their effect on treatment safety, are discussed.

A relic of the past: current genetic patterns of the palaeoendemic tree Nothofagus macrocarpa were shaped by climatic oscillations in central Chile.

BACKGROUND AND AIMS: The Mediterranean-type forest of central Chile is considered a 'biodiversity hotspot' and a relic of a wider ancient distribution produced by past climatic oscillations. Nothofagus macrocarpa, commonly known as 'roble de Santiago', is a threatened palaeoendemic of this forest, poorly represented in the protected area system. This tree has been repeatedly misidentified as the sister species N. obliqua, which has affected its recognition and protection. Only a few populations of N. macrocarpa remain within a matrix of intensive land use that has been affected by recent forest fires. We tested the hypothesis that current populations of N. macrocarpa are a relic state of a previously widespread range, with the aim of contributing to its identification, its biogeographical history and the design of conservation measures using genetic information. METHODS: We analysed remnant N. macrocarpa forests using nuclear (nDNA) and chloroplast DNA (cpDNA) sequences, conducted phylogenetic and phylogeographical analyses to reconstruct its biogeographical history, and assessed microsatellites [simple sequence repeats (SSRs)] to determine contemporary patters of diversity within and among all remnant populations. We also examined the degree of past, current and potential future isolation of N. macrocarpa populations using ecological niche models (ENMs). KEY RESULTS: The species N. macrocarpa was confirmed by nDNA sequences, as previously suggested by chromosomal analysis. Small isolated populations of N. macrocarpa exhibited moderate to high genetic diversity according to SSRs. cpDNA analysis revealed a marked past latitudinal geographical structure, whereas analysis of SSRs did not find such current structure. ENM analyses revealed local expansion-contraction of the N. macrocarpa range during warmer periods, particularly in the northern and central ranges where basal-most cpDNA haplotypes were detected, and recent expansion to the south of the distribution. CONCLUSIONS: Genetic patterns confirm that N. macrocarpa is a distinct species and suggest a marked latitudinal relic structure in at least two evolutionarily significant units, despite contemporary among-population gene flow. This information must be considered when choosing individuals (seeds and/or propagules) for restoration purposes, to avoid the admixture of divergent genetic stocks.

Isolation and molecular characterization of four novel Neospora caninum strains.

Neospora caninum causes neosporosis, a leading cause of bovine abortion worldwide. Uruguay is a developing economy in South America that produces milk to feed seven times its population annually. Naturally, dairy production is paramount to the country's economy, and bovine reproductive failure impacts it profoundly. Recent studies demonstrated that the vast majority of infectious abortions in dairy cows are caused by N. caninum. To delve into the local situation and contextualize it within the international standing, we set out to characterize the Uruguayan N. caninum strains. For this, we isolated four distinct strains and determined by microsatellite typing that these represent three unique genetic lineages, distinct from those reported previously in the region or elsewhere. An unbiased analysis of the current worldwide genetic diversity of N. caninum strains known, whereby six typing clusters can be resolved, revealed that three of the four Uruguayan strains group closely with regional strains from Argentina and Brazil. The remaining strain groups in an unrelated genetic cluster, suggesting multiple origins of the local strains. Microsatellite typing of N. caninum DNA from fetuses opportunistically collected from local dairy farms correlated more often with one of the isolates. Overall, our results contribute to further understanding of genetic diversity among strains of N. caninum both regionally and worldwide.

Whole genome sequencing of the monomorphic pathogen Mycobacterium bovis reveals local differentiation of cattle clinical isolates.

BACKGROUND: Bovine tuberculosis (bTB) poses serious risks to animal welfare and economy, as well as to public health as a zoonosis. Its etiological agent, Mycobacterium bovis, belongs to the Mycobacterium tuberculosis complex (MTBC), a group of genetically monomorphic organisms featured by a remarkably high overall nucleotide identity (99.9%). Indeed, this characteristic is of major concern for correct typing and determination of strain-specific traits based on sequence diversity. Due to its historical economic dependence on cattle production, Uruguay is deeply affected by the prevailing incidence of Mycobacterium bovis. With the world's highest number of cattle per human, and its intensive cattle production, Uruguay represents a particularly suited setting to evaluate genomic variability among isolates, and the diversity traits associated to this pathogen. RESULTS: We compared 186 genomes from MTBC strains isolated worldwide, and found a highly structured population in M. bovis. The analysis of 23 new M. bovis genomes, belonging to strains isolated in Uruguay evidenced three groups present in the country. Despite presenting an expected highly conserved genomic structure and sequence, these strains segregate into a clustered manner within the worldwide phylogeny. Analysis of the non-pe/ppe differential areas against a reference genome defined four main sources of variability, namely: regions of difference (RD), variable genes, duplications and novel genes. RDs and variant analysis segregated the strains into clusters that are concordant with their spoligotype identities. Due to its high homoplasy rate, spoligotyping failed to reflect the true genomic diversity among worldwide representative strains, however, it remains a good indicator for closely related populations. CONCLUSIONS: This study introduces a comprehensive population structure analysis of worldwide M. bovis isolates. The incorporation and analysis of 23 novel Uruguayan M. bovis genomes, sheds light onto the genomic diversity of this pathogen, evidencing the existence of greater genetic variability among strains than previously contemplated.

Historical isolation and contemporary gene flow drive population diversity of the brown alga Sargassum thunbergii along the coast of China.

BACKGROUND: Long-term survival in isolated marginal seas of the China coast during the late Pleistocene ice ages is widely believed to be an important historical factor contributing to population genetic structure in coastal marine species. Whether or not contemporary factors (e.g. long-distance dispersal via coastal currents) continue to shape diversity gradients in marine organisms with high dispersal capability remains poorly understood. Our aim was to explore how historical and contemporary factors influenced the genetic diversity and distribution of the brown alga Sargassum thunbergii, which can drift on surface water, leading to long-distance dispersal. RESULTS: We used 11 microsatellites and the plastid RuBisCo spacer to evaluate the genetic diversity of 22 Sargassum thunbergii populations sampled along the China coast. Population structure and differentiation was inferred based on genotype clustering and pairwise F ST and allele-frequency analyses. Integrated genetic analyses revealed two genetic clusters in S. thunbergii that dominated in the Yellow-Bohai Sea (YBS) and East China Sea (ECS) respectively. Higher levels of genetic diversity and variation were detected among populations in the YBS than in the ECS. Bayesian coalescent theory was used to estimate contemporary and historical gene flow. High levels of contemporary gene flow were detected from the YBS (north) to the ECS (south), whereas low levels of historical gene flow occurred between the two regions. CONCLUSIONS: Our results suggest that the deep genetic divergence in S. thunbergii along the China coast may result from long-term geographic isolation during glacial periods. The dispersal of S. thunbergii driven by coastal currents may facilitate the admixture between southern and northern regimes. Our findings exemplify how both historical and contemporary forces are needed to understand phylogeographical patterns in coastal marine species with long-distance dispersal.

Genomic and clinical evidence uncovers the enterohepatic species Helicobacter valdiviensis as a potential human intestinal pathogen.

BACKGROUND: Helicobacter valdiviensis is a recently described enterohepatic species isolated from wild bird's fecal samples. Currently, its pathogenic potential and clinical significance are unknown mainly due to the lack of whole-genome sequences to compare with other helicobacters and the absence of specific screenings to determine its prevalence in humans. MATERIALS AND METHODS: The species type strain (WBE14T ) was whole-genome-sequenced, and comparative analyses were carried out including the genomes from other Helicobacter species to determine the exact phylogenetic position of H. valdiviensis and to study the presence and evolution of virulence determinants. In parallel, stools from diarrheic patients and healthy individuals were screened by PCR to assess the clinical incidence of H. valdiviensis. RESULTS: Helicobacter valdiviensis belongs to a monophyletic clade conformed by H. canadensis, H. pullorum, H. winghamensis, H. rodentium, and H. apodemus. Its predicted genome size is 2 176 246 bp., with 30% of G+C content and 2064 annotated protein-coding genes. The patterns of virulence factors in H. valdiviensis were similar to other enterohepatic species, but evidence of horizontal gene transfer from Campylobacter species was detected for key genes like those coding for the CDT subunits. Positive PCR results confirmed the presence of H. valdiviensis in 2 of 254 (0.78%) stools of patients with acute diarrhea while not a single sample was positive in healthy individuals. CONCLUSIONS: Horizontal gene transfer has contributed to shape the gene repertory of H. valdiviensis, which codes for virulence factors conserved in other pathogens that are well-known human pathogens. Additionally, the detection of H. valdiviensisDNA in diarrheic patients supports its role as a potential emergent intestinal pathogen. Further, sampling efforts are needed to uncover the clinical relevance of this species, which should be accomplished by the isolation of H. valdiviensis from ill humans and the obtention of whole genomes from clinical isolates.

Comparative transcriptome analysis of lufenuron-resistant and susceptible strains of Spodoptera frugiperda (Lepidoptera: Noctuidae).

BACKGROUND: The evolution of insecticide resistance in Spodoptera frugiperda (Lepidoptera: Noctuidae) has resulted in large economic losses and disturbances to the environment and agroecosystems. Resistance to lufenuron, a chitin biosynthesis inhibitor insecticide, was recently documented in Brazilian populations of S. frugiperda. Thus, we utilized large-scale cDNA sequencing (RNA-Seq analysis) to compare the pattern of gene expression between lufenuron-resistant (LUF-R) and susceptible (LUF-S) S. larvae in an attempt to identify the molecular basis behind the resistance mechanism(s) of S. frugiperda to this insecticide. RESULTS: A transcriptome was assembled using approximately 19.6 million 100 bp-long single-end reads, which generated 18,506 transcripts with a N50 of 996 bp. A search against the NCBI non-redundant database generated 51.1% (9,457) functionally annotated transcripts. A large portion of the alignments were homologous to insects, with the majority (45%) being similar to sequences of Bombyx mori (Lepidoptera: Bombycidae). Moreover, 10% of the alignments were similar to sequences of various species of Spodoptera (Lepidoptera: Noctuidae), with 3% of them being similar to sequences of S. frugiperda. A comparative analysis of the gene expression between LUF-R and LUF-S S. frugiperda larvae identified 940 differentially expressed transcripts (p ≤ 0.05, t-test; fold change ≥ 4). Six of them were associated with cuticle metabolism. Of those, four were overexpressed in LUF-R larvae. The machinery involved with the detoxification process was represented by 35 differentially expressed transcripts; 24 of them belonging to P450 monooxygenases, four to glutathione-S-transferases, six to carboxylases and one to sulfotransferases. RNA-Seq analysis was validated for a number of selected candidate transcripts by using quantitative real time PCR (qPCR). CONCLUSIONS: The gene expression profile of LUF-R larvae of S. frugiperda differs from LUF-S larvae. In general, gene expression is much higher in resistant larvae when compared to the susceptible ones, particularly for those genes involved with pathways for xenobiotic detoxification, mainly represented by P450 monooxygenases transcripts. Our data indicate that enzymes involved with the detoxification process, and mostly the P450, are one of the resistance mechanisms employed by the LUF-R S. frugiperda larvae against lufenuron.

Phylogeographic structure and deep lineage diversification of the red alga Chondrus ocellatus Holmes in the Northwest Pacific.

A major goal of phylogeographic analysis using molecular markers is to understand the ecological and historical variables that influence genetic diversity within a species. Here, we used sequences of the mitochondrial Cox1 gene and nuclear internal transcribed spacer to reconstruct its phylogeography and demographic history of the intertidal red seaweed Chondrus ocellatus over most of its geographical range in the Northwest Pacific. We found three deeply separated lineages A, B and C, which diverged from one another in the early Pliocene-late Miocene (c. 4.5-7.7 Ma). The remarkably deep divergences, both within and between lineages, appear to have resulted from ancient isolations, accelerated by random drift and limited genetic exchange between regions. The disjunct distributions of lineages A and C along the coasts of Japan may reflect divergence during isolation in scattered refugia. The distribution of lineage B, from the South China Sea to the Korean Peninsula, appears to reflect postglacial recolonizations of coastal habitats. These three lineages do not coincide with the three documented morphological formae in C. ocellatus, suggesting that additional cryptic species may exist in this taxon. Our study illustrates the interaction of environmental variability and demographic processes in producing lineage diversification in an intertidal seaweed and highlights the importance of phylogeographic approaches for discovering cryptic marine biodiversity.

Phenotypic, genomic and in planta characterization of Bacillus sensu lato for their phosphorus biofertilization and plant growth promotion features in soybean.

Bacillus sensu lato were screened for their capacity to mineralize organic phosphorus (P) and promote plant growth, improving nitrogen (N) and P nutrition of soybean. Isolates were identified through Type Strain Genome Server (TYGS) and Average Nucleotide Identity (ANI). ILBB95, ILBB510 and ILBB592 were identified as Priestia megaterium, ILBB139 as Bacillus wiedmannii, ILBB44 as a member of a sister clade of B. pumilus, ILBB15 as Peribacillus butanolivorans and ILBB64 as Lysinibacillus sp. These strains were evaluated for their capacity to mineralize sodium phytate as organic P and solubilize inorganic P in liquid medium. These assays ranked ILBB15 and ILBB64 with the highest orthophosphate production from phytate. Rhizocompetence and plant growth promotion traits were evaluated in vitro and in silico. Finally, plant bioassays were conducted to assess the effect of the co-inoculation with rhizobial inoculants on nodulation, N and P nutrition. These bioassays showed that B. pumilus, ILBB44 and P. megaterium ILBB95 increased P-uptake in plants on the poor substrate of sand:vermiculite and also on a more fertile mix. Priestia megaterium ILBB592 increased nodulation and N content in plants on the sand:vermiculite:peat mixture. Peribacillus butanolivorans ILBB15 reduced plant growth and nutrition on both substrates. Genomes of ILBB95 and ILBB592 were characterized by genes related with plant growth and biofertilization, whereas ILBB15 was differentiated by genes related to bioremediation. Priestia megaterium ILBB592 is considered as nodule-enhancing rhizobacteria and together with ILBB95, can be envisaged as prospective PGPR with the capacity to exert positive effects on N and P nutrition of soybean plants.

Datathons: fostering equitability in data reuse in ecology.

Approaches to rapidly collecting global biodiversity data are increasingly important, but biodiversity blind spots persist. We organized a three-day Datathon event to improve the openness of local biodiversity data and facilitate data reuse by local researchers. The first Datathon, organized among microbial ecologists in Uruguay and Argentina assembled the largest microbiome dataset in the region to date and formed collaborative consortia for microbiome data synthesis.

Genetic variability highlights the invasion route of the Lutzomyia longipalpis complex, the main vector of Visceral Leishmaniasis in Uruguay.

In the Americas, the sandfly Lutzomyia longipalpis is the main vector of the parasitic protozoa Leishmania infantum, the etiological agent of visceral leishmaniasis (VL). The Lu. longipalpis species complex is currently discontinuously distributed across the Neotropical region, from Mexico to the north of Argentina and Uruguay. During its continental spreading, it must have adapted to several biomes and temperature amplitudes, when founder events should have contributed to the high genetic divergence and geographical structure currently observed, reinforcing the speciation process. The first report of Lu. longipalpis in Uruguay was in 2010, calling the attention of Public Health authorities. Five years later, the parasite Le. infantum was recorded and in 2015 the first case of VL in canids was reported. Hitherto seven human cases by VL have been reported in Uruguay. Here, we publish the first DNA sequences from the mitochondrial genes ND4 and CYTB of Lu. longipalpis collected in Uruguay, and we used these molecular markers to investigate their genetic variability and population structure. We described four new ND4 haplotypes in a total of 98 (4/98) and one CYTB in a total of 77 (1/77). As expected, we were able to establish that the Lu. longipalpis collected in two localities (i.e. Salto and Bella Unión) from the north of Uruguay are closely related to the populations from neighbouring countries. We also propose that the possible route for the vector arrival to the region may have been through vegetation and forest corridors of the Uruguay River system, as well as it may have benefited from landscape modifications generated by commercial forestation. The ecological-scale processes shaping Lu. longipalpis populations, the identification of genetically homogeneous groups and the gene flow among them must be carefully investigated by using highly sensible molecular markers (i.e. genome wide SNPs) since it will help to the understanding of VL transmission and contribute to the planification of public policies on its control.

RNAi in Piezodorus guildinii (Hemiptera: Pentatomidae): Transcriptome Assembly for the Development of Pest Control Strategies.

Red-banded stink bug Piezodorus guildinii (P. guildinii) has been described as the most damaging stink bug regarding soybean crops, leading to seed injury, low germination percentages, and foliar retention, at low population densities. In recent years, RNA interference (RNAi), a conserved eukaryote silencing mechanism has been explored to develop species-selective pesticides. In this work, we evaluated RNAi in P. guildinii to develop new pest-control strategies. For this, we assembled and annotated a P. guildinii transcriptome from a pool of all developmental stages. Analysis of this transcriptome led to the identification of 56 genes related to the silencing process encompassing siRNA, miRNA, and piRNA pathways. To evaluate the functionality of RNAi machinery, P. guildinii adults were injected with 28 ng/mg of body weight of double stranded RNA (dsRNA) targeting vATPase A. A mortality of 35 and 51.6% was observed after 7 and 14 days, respectively, and a downregulation of vATPase A gene of 84% 72 h post-injection. In addition, Dicer-2 and Argonaute-2 genes, core RNAi proteins, were upregulated 1.8-fold 48 h after injection. These findings showed for the first time that RNAi is functional in P. guildinii and the silencing of essential genes has a significant effect in adult viability. Taken together, the work reported here shows that RNAi could be an interesting approach for the development of red-banded stink bug control strategies.

Phylogenetic and Multiple-Locus Variable number tandem repeat analysis of Mycobacterium avium subsp. paratuberculosis isolates from Argentina.

Paratuberculosis is a worldwide chronic enteric disease of ruminants, caused by Mycobacterium avium subsp. paratuberculosis (MAP). While MAP has been widely investigated all around the world, little is known about the different strains that circulate in each country. This study describes the genetic diversity of MAP isolates from different bovine and deer herds from Argentina, analyzed by Multiple-Locus Variable number tandem repeat Analysis (MLVA), as well as the phylogenetic relatedness between geographically distant isolates through Whole Genome Sequencing (WGS) and core-genome analysis. A total of 90 MAP isolates were analyzed. The results showed seven different MLVA genotypes, with almost 75% of them belonging to pattern INMV 1, described in all the herds studied. WGS results suggested the presence of a common INMV 1 strain circulating throughout the country. Our results allow confirming the coexistence of different strains in time and space and the mixed infections identified in some animals. These observations suggest the absence of animal monitoring prior to introduction to the herds and the need for a control program in the country. This study represents the first to report WGS of MAP strains in Argentina.

Genetic structure and demographic history of new world screwworm across its current geographic range.

The phylogeographical history of the pest fly screwworm, Cochliomyia hominivorax (Coquerel), was studied using partial mitochondrial DNA sequences of the control region, Cytochrome c oxidase (CO) subunit I and CO subunit II from 361 individuals collected across its current geographic range. Analyses showed marked genetic differentiation on a macrogeographic scale. The genetic diversity in the species is structured into four main "regional groups," corresponding to Cuba, the Dominican Republic, and the North and South Amazon region. Results indicated that the distribution of screwworm genetic diversity was mainly shaped by historical events, i.e., colonization of Caribbean islands, vicariance in the Amazon region and population expansion. Demographic history analyses revealed that the population expansion started approximately 20-25,000 yr ago and recently increased exponentially. We hypothesized that the initial period of expansion was probably associated with environmental amelioration in the late Pleistocene and the exponential increase with resource availability in recent times. The population expansion is probably responsible for the low divergence and the lack of genetic and geographic correlation in the South Amazon region but did not erase the genetic structure pattern on a continental scale. The screwworm is one of the most damaging livestock pests in South and Central America, and the pattern of genetic variability distribution reported here suggests that the Caribbean area and the North and South Amazon regions could be considered as independent units for future pest control programs.

Pangenome analysis reveals genetic isolation in Campylobacter hyointestinalis subspecies adapted to different mammalian hosts.

Campylobacter hyointestinalis is an emerging pathogen currently divided in two subspecies: C. hyointestinalis subsp. lawsonii which is predominantly recovered from pigs, and C. hyointestinalis subsp. hyointestinalis which can be found in a much wider range of mammalian hosts. Despite C. hyointestinalis being reported as an emerging pathogen, its evolutionary and host-associated diversification patterns are still vastly unexplored. For this reason, we generated whole-genome sequences of 13 C. hyointestinalis subsp. hyointestinalis strains and performed a comprehensive comparative analysis including publicly available C. hyointestinalis subsp. hyointestinalis and C. hyointestinalis subsp. lawsonii genomes, to gain insight into the genomic variation of these differentially-adapted subspecies. Both subspecies are distinct phylogenetic lineages which present an apparent barrier to homologous recombination, suggesting genetic isolation. This is further supported by accessory gene patterns that recapitulate the core genome phylogeny. Additionally, C. hyointestinalis subsp. hyointestinalis presents a bigger and more diverse accessory genome, which probably reflects its capacity to colonize different mammalian hosts unlike C. hyointestinalis subsp. lawsonii that is presumably host-restricted. This greater plasticity in the accessory genome of C. hyointestinalis subsp. hyointestinalis correlates to a higher incidence of genome-wide recombination events, that may be the underlying mechanism driving its diversification. Concordantly, both subspecies present distinct patterns of gene families involved in genome plasticity and DNA repair like CRISPR-associated proteins and restriction-modification systems. Together, our results provide an overview of the genetic mechanisms shaping the genomes of C. hyointestinalis subspecies, contributing to understand the biology of Campylobacter species that are increasingly recognized as emerging pathogens.

Reduced tillage, cover crops and organic amendments affect soil microbiota and improve soil health in Uruguayan vegetable farming systems.

Conventional tillage and mineral fertilization (CTMF) jeopardize soil health in conventional vegetable production systems. Using a field experiment established in Uruguay in 2012, we aimed to compare the soil restoration potential of organic fertilization (compost and poultry manure) combined with conventional tillage and cover crop incorporated into the soil (CTOF) or with reduced tillage and the use of cover crop as mulch (RTOF). In 2017, table beet was cultivated under CTMF, CTOF and RTOF, and yields, soil aggregate composition and nutrients, as well as soil and table beet rhizosphere microbiota (here: bacteria and archaea) were evaluated. Microbiota was studied by high-throughput sequencing of 16S rRNA gene fragments amplified from total community DNA. RTOF exhibited higher soil aggregation, soil organic C, nutrient availability and microbial alpha-diversity than CTMF, and became more similar to an adjacent natural undisturbed site. The soil microbiota was strongly shaped by the fertilization source which was conveyed to the rhizosphere and resulted in differentially abundant taxa. However, 229 amplicon sequencing variants were found to form the core table beet rhizosphere microbiota shared among managements. In conclusion, our study shows that after only 5 years of implementation, RTOF improves soil health under intensive vegetable farming systems.