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.

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.

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.

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.

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.

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.

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.

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.

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.

Genetic Differentiation of a New World Screwworm Fly Population from Uruguay Detected by SNPs, Mitochondrial DNA and Microsatellites in Two Consecutive Years.

The New World screwworm (NWS) fly, Cochliomyia hominivorax (Diptera: Calliphoridae), is an economically important ectoparasite currently distributed in South America and in the Caribbean basin. The successful eradication of this species in USA, Mexico and continental Central America was achieved by a control program based on the sterile insect technique (SIT). In order to implement a genetic control strategy over the NWS fly's current area of occurrence, first, it is necessary to understand the species dynamics and population structure. In order to address this objective, the spatial genetic structure of the NWS fly was previously reported in South America based on different genetic markers; however, to date, no study has investigated temporal changes in the genetic composition of its populations. In the current study, the temporal genetic structure of a NWS fly population from Uruguay was investigated through two consecutive samplings from the same locality over an interval of approximately 18 generations. The genetic structure was accessed with neutral and under selection SNPs obtained with genotyping-by-sequencing. The results gathered with these data were compared to estimates achieved with mitochondrial DNA sequences and eight microsatellite markers. Temporal changes in the genetic composition were revealed by all three molecular markers, which may be attributed to seasonal changes in the NWS fly's southern distribution. SNPs were employed for the first time for estimating the genetic structure in a NWS fly population; these results provide new clues and perspectives on its population genetic structure. This approach could have significant implications for the planning and implementation of management programs.

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.

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.

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.

Urban metagenomics uncover antibiotic resistance reservoirs in coastal beach and sewage waters.

BACKGROUND: Microbial communities present in environmental waters constitute a reservoir for antibiotic-resistant pathogens that impact human health. For this reason, a diverse variety of water environments are being analyzed using metagenomics to uncover public health threats. However, the composition of these communities along the coastal environment of a whole city, where sewage and beach waters are mixed, is poorly understood. RESULTS: We shotgun-sequenced 20 coastal areas from the city of Montevideo (capital of Uruguay) including beach and sewage water samples to characterize bacterial communities and their virulence and antibiotic resistance repertories. As expected, we found that sewage and beach environments present significantly different bacterial communities. This baseline allowed us to detect a higher prevalence and a more diverse repertory of virulence and antibiotic-resistant genes in sewage samples. Many of these genes come from well-known enterobacteria and represent carbapenemases and extended-spectrum betalactamases reported in hospital infections in Montevideo. Additionally, we were able to genotype the presence of both globally disseminated pathogenic clones and emerging antibiotic-resistant bacteria in sewage waters. CONCLUSIONS: Our study represents the first in using metagenomics to jointly analyze beaches and the sewage system from an entire city, allowing us to characterize antibiotic-resistant pathogens circulating in urban waters. The data generated in this initial study represent a baseline metagenomic exploration to guide future longitudinal (time-wise) studies, whose systematic implementation will provide useful epidemiological information to improve public health surveillance.

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.

High Genetic Diversity and No Population Structure of the New World Screwworm Fly Cochliomyia hominivorax (Diptera: Calliphoridae) on a Microgeographic Scale: Implications for Management Units.

The New World screwworm fly Cochliomyia hominivorax (Coquerel, 1858) (Diptera: Calliphoridae) is an important livestock pest endemic to the Americas that has been eradicated from North and continental Central America with a control program based on the Sterile Insect Technique (SIT). The establishment of target management units is a strategic step in the implementation of new control programs, which can be achieved using genetic studies of natural populations. Previous studies of New World screwworm fly populations were conducted on the continental scale and identified four main groups: two in South America and two in the Caribbean. However, studies within these groups are needed to determine which smaller geographic areas can be treated as management units. Here, we analyze the genetic variability distribution and the population demographic signals of the New World screwworm fly in a 6,000 km2 area located along the border of Brazil and Uruguay. This area was the subject of the first control pilot program conducted in South America. We studied eight microsatellite loci and sequences from two mitochondrial DNA regions in individuals sampled at 20-25 livestock breeding farms. We observed no population structure and found high genetic variability on the geographical scale sampled for both molecular markers. Our microsatellite data suggest that these populations are not in equilibrium, and demographic analyses based on mitochondrial data indicate population expansion. These results suggest that this geographic scale is not adequate for future New World screwworm fly management in South America.

Sequencing and Annotation of the Genome of Mycobacterium tuberculosis MYC004, a Strain Causing Meningitis in Mexico.

Mycobacterium tuberculosis strain MYC004 was isolated from a Mexican patient with tuberculous meningitis, the most aggressive form of tuberculosis. The draft genome sequence is the first of a meningeal strain of M. tuberculosis reported from Latin America and consists of 4,411,530 bp, including 4,251 protein-encoding genes.

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.