A new spike-in-based method for quantitative metabarcoding of soil fungi and bacteria.

Metabarcoding is a powerful tool to characterize biodiversity in biological samples. The interpretation of taxonomic profiles from metabarcoding data has been hindered by their compositional nature. Several strategies have been proposed to transform compositional data into quantitative data, but they have intrinsic limitations. Here, I propose a workflow based on bacterial and fungal cellular internal standards (spike-ins) for absolute quantification of the microbiota in soil samples. These standards were added to the samples before DNA extraction in amounts estimated after qPCRs, to target around 1-2% coverage in the sequencing run. In bacteria, proportions of spike-in reads in the sequencing run were very similar (< 2-fold change) to those predicted by the qPCR assessment, and for fungi they differed up to 40-fold. The low variation among replicates highlights the reproducibility of the method. Estimates based on multiple bacterial spike-ins were highly correlated (r = 0.99). Procrustes analysis evidenced significant biological effects on the community composition when normalizing compositional data. A protocol based on qPCR estimation of input amounts of cellular spikes is proposed as a cheap and reliable strategy for quantitative metabarcoding of biological samples.

Phenotypic Variability and Genetic Diversity of the Pathogenic Fungus Macrophomina phaseolina from Several Hosts and Host Specialization in Strawberry.

Macrophomina phaseolina, is a pathogenic soil-borne fungus that affects more than 500 plant species, causing various types of disease to several crops, among which is the crown and root rot disease in strawberry. Its wide variability has been characterized reiteratively in the literature, but little is known about its virulence mechanisms. Morphological, physiological, genetic and phytopathogenic parameters were evaluated among 32 isolates of Macrophomina from different hosts occurring in Argentina and Spain. Colony characteristics, average size of microsclerotia, chlorate phenotype and mycelial growth at different temperatures (5º-40 °C), and pH (3.0-8.0) were recorded. The morphological and physiological traits were heterogeneous and did not show any association with the genetic structure nor with their pathogenicity. Most of the isolates (71.9%) exhibited chlorate-sensitive phenotype. The optimal growth temperature range was between 25 °C and 35 °C, and the optimal pH varied between 4.0 and 6.0. The genetic structure analyzed with four DNA markers (EF-1α, ITS, CAL and TUB) showed little diversity among isolates of M. phaseolina, with no clear association with the country of origin, but a significant association with the host. Based on their phylogenetic affinity, one isolate was reclassified as M. pseudophaseolina and another one as M. tecta. It is the first report of M. pseudophaseolina causing charcoal rot on beans, in Argentina, and the first report of M. tecta outside Australia. Pathogenicity tests on strawberry plants revealed marked host specialization, being the isolates obtained from strawberry more virulent than those from other hosts.

Evaluating surrogates of genetic diversity for conservation planning.

Protected-area systems should conserve intraspecific genetic diversity. Because genetic data require resources to obtain, several approaches have been proposed for generating plans for protected-area systems (prioritizations) when genetic data are not available. Yet such surrogate-based approaches remain poorly tested. We evaluated the effectiveness of potential surrogate-based approaches based on microsatellite genetic data collected across the Iberian Peninsula for 7 amphibian and 3 reptilian species. Long-term environmental suitability did not effectively represent sites containing high genetic diversity (allelic richness). Prioritizations based on long-term environmental suitability had similar performance to random prioritizations. Geographic distances and resistance distances based on contemporary environmental suitability were not always effective surrogates for identification of combinations of sites that contain individuals with different genetic compositions. Our results demonstrate that population genetic data based on commonly used neutral markers can inform prioritizations, and we could not find an adequate substitute. Conservation planners need to weigh the potential benefits of genetic data against their acquisition costs.

Evaluación de los Sustitutos de la Diversidad Genética para la Planeación de la Conservación Resumen Los sistemas de áreas protegidas deberían conservar la diversidad genética intraespecífica. Ya que para obtener datos genéticos se requieren recursos, se han propuesto distintas estrategias para generar los planes para los sistemas de áreas protegidas (priorizaciones) cuando los datos genéticos no están disponibles. A pesar de lo anterior, dichas estrategias basadas en sustitutos han sido poco evaluadas. Evaluamos la efectividad del potencial de las estrategias basadas en sustitutos cuya base son los datos genéticos de microsatélites obtenidos en toda la Península Ibérica y correspondientes a siete especies de anfibios y a tres de reptiles. La idoneidad ambiental a largo plazo no representó efectivamente los sitios que contienen una diversidad genética alta (riqueza de alelos). Las priorizaciones basadas en la idoneidad ambiental a largo plazo tuvieron un desempeño similar a las priorizaciones aleatorias. Las distancias geográficas y las distancias de resistencia basadas en la idoneidad ambiental contemporánea no siempre fueron sustitutos efectivos para la identificación de las combinaciones de sitios que contienen individuos con composiciones genéticas diferentes. Nuestros resultados demuestran que los datos genéticos de una población basados en marcadores neutrales de uso común pueden informar a las priorizaciones y que no pudimos encontrar un sustituto adecuado. Los planificadores de la conservación necesitan sopesar los beneficios potenciales de los datos genéticos contra sus costos de adquisición.

Little genetic structure in a Bornean endemic small mammal across a steep ecological gradient.

Janzen's influential "mountain passes are higher in the tropics" hypothesis predicts restricted gene flow and genetic isolation among populations spanning elevational gradients in the tropics. Few studies have tested this prediction, and studies that focus on population genetic structure in Southeast Asia are particularly underrepresented in the literature. Here, we test the hypothesis that mountain treeshrews (Tupaia montana) exhibit limited dispersal across their broad elevational range which spans ~2,300 m on two peaks in Kinabalu National Park (KNP) in Borneo: Mt Tambuyukon (MT) and Mt Kinabalu (MK). We sampled 83 individuals across elevations on both peaks and performed population genomics analyses on mitogenomes and single nucleotide polymorphisms from 4,106 ultraconserved element loci. We detected weak genetic structure and infer gene flow both across elevations and between peaks. We found higher genetic differentiation on MT than MK despite its lower elevation and associated environmental variation. This implies that, contrary to our hypothesis, genetic structure in this system is not primarily shaped by elevation. We propose that this pattern may instead be the result of historical processes and limited upslope gene flow on MT. Importantly, our results serve as a foundational estimate of genetic diversity and population structure from which to track potential future effects of climate change on mountain treeshrews in KNP, an important conservation stronghold for the mountain treeshrew and other montane species.

Mitogenomes Reveal Multiple Colonization of Mountains by Rattus in Sundaland.

Tropical mountains are cradles of biodiversity and endemism. Sundaland, tropical Southeast Asia, hosts 3 species of Rattus endemic to elevations above 2000 m with an apparent convergence in external morphology: Rattus korinchi and R. hoogerwerfi from Sumatra, and R. baluensis from Borneo. A fourth one, R. tiomanicus, is restricted to lowland elevations across the whole region. The origins of these endemics are little known due to the absence of a robust phylogenetic framework. We use complete mitochondrial genomes from the 3 high altitude Rattus, and several related species to determine their relationships, date divergences, reconstruct their history of colonization, and test for selection on the mitochondrial DNA. We show that mountain colonization happened independently in Borneo (<390 Kya) and Sumatra (~1.38 Mya), likely from lowland lineages. The origin of the Bornean endemic R. baluensis is very recent and its genetic diversity is nested within the diversity of R. tiomanicus. We found weak evidence of positive selection in the high-elevation lineages and attributed the greater nonsynonymous mutations on these branches (specially R. baluensis) to lesser purifying selection having acted on the terminal branches in the phylogeny.

Soil Microbial Community Responses to Different Management Strategies in Almond Crop.

A comparative study of organic and conventional farming systems was conducted in almond orchards to determine the effect of management practices on their fungal and bacterial communities. Soils from two orchards under organic (OM) and conventional (CM), and nearby nonmanaged (NM) soil were analyzed and compared. Several biochemical and biological parameters were measured (soil pH, electrical conductivity, total nitrogen, organic material, total phosphorous, total DNA, and fungal and bacterial DNA copies). Massive parallel sequencing of regions from fungal ITS rRNA and bacterial 16 S genes was carried out to characterize their diversity in the soil. We report a larger abundance of bacteria and fungi in soils under OM, with a more balanced fungi:bacteria ratio, compared to bacteria-skewed proportions under CM and NM. The fungal phylum Ascomycota corresponded to around the 75% relative abundance in the soil, whereas for bacteria, the phyla Proteobacteria, Acidobacteriota and Bacteroidota integrated around 50% of their diversity. Alpha diversity was similar across practices, but beta diversity was highly clustered by soil management. Linear discriminant analysis effect size (LEfSE) identified bacterial and fungal taxa associated with each type of soil management. Analyses of fungal functional guilds revealed 3-4 times larger abundance of pathogenic fungi under CM compared to OM and NM treatments. Among them, the genus Cylindrocarpon was more abundant under CM, and Fusarium under OM.

Bacterial assemblage in Mediterranean salt marshes: Disentangling the relative importance of seasonality, zonation and halophytes.

Salt marshes gather a high diversity of prokaryotes across their environmental gradients. Most of this diversity and the factors determining their community assemblage are unknown. We massively sequenced a portion of the 16S gene to characterize the diversity of prokaryotes in soils from a salt marsh in Río Piedras, Southern Spain. We sampled in the four seasons, and in five plots dominated by a different halophyte (Spartina maritima, S. densiflora, Salicornia ramosissima, Arthrocaulon macrostachyum and Atriplex portulacoides) growing under different environmental conditions and representing different stages in the marsh ecological succession. Soil was sampled in their rhizosphere and adjacent bulk soil. We report the effects of different factors explaining prokaryotic beta diversity in the marsh: zonation (50 %), seasonality (14 %), and halophyte rhizosphere (7 %). Proteobacteria and Bacteroidota were the most abundant phyla. Firmicutes had a peak in winter and Desulfobacterota with other bacteria involved in sulfur cycling were abundant in the low marsh plots from S. maritima. Alpha diversity was highest in spring and decreased in winter. We detected a marked phylogenetic turnover between seasons and in rhizospheric soil respect to adjacent bulk soil for most pairwise comparisons. The effect of halophyte on its rhizosphere was species-specific, being S. maritima the species with more differentiated taxa between rhizosphere versus surrounding bulk soil. Our work highlights how the complex interaction between marsh zonation, seasonality and rhizosphere, onsets processes structuring bacterial community assemblage in salt marsh soils.

Challenging ecogeographical rules: Phenotypic variation in the Mountain Treeshrew (Tupaia montana) along tropical elevational gradients.

Bergmann's and Allen's rules were defined to describe macroecological patterns across latitudinal gradients. Bergmann observed a positive association between body size and latitude for endothermic species while Allen described shorter appendages as latitude increases. Almost two centuries later, there is still ongoing discussion about these patterns. Temperature, the common variable in these two rules, varies predictably across both latitude and elevation. Although these rules have been assessed extensively in mammals across latitude, particularly in regions with strong seasonality, studies on tropical montane mammals are scarce. We here test for these patterns and assess the variation of several other locomotory, diet-associated, body condition, and thermoregulatory traits across elevation in the Mountain Treeshrew (Tupaia montana) on tropical mountains in Borneo. Based on morphological measurements from both the field and scientific collections, we found a complex pattern: Bergmann's rule was not supported in our tropical mountain system, since skull length, body size, and weight decreased from the lowest elevations (<1000 m) to middle elevations (2000-2500 m), and then increased from middle elevations to highest elevations. Allen's rule was supported for relative tail length, which decreased with elevation, but not for ear and hindfoot length, with the former remaining constant and the latter increasing with elevation. This evidence together with changes in presumed diet-related traits (rostrum length, zygomatic breadth and upper tooth row length) along elevation suggest that selective pressures other than temperature, are playing a more important role shaping the morphological variation across the distribution of the Mountain Treeshrew. Diet, food acquisition, predation pressure, and/or intra- and inter-specific competition, are some of the potential factors driving the phenotypic variation of this study system. The lack of variation in body condition might suggest local adaptation of this species across its elevational range, perhaps due to generalist foraging strategies. Finally, a highly significant temporal effect was detected in several traits but not in others, representing the first phenotypic variation temporal trends described on treeshrews.

Markers for genetic change.

BACKGROUND AND AIMS: Wildlife conservation has focused primarily on species for the last decades. Recently, popular perception and laws have begun to recognize the central importance of genetic diversity in the conservation of biodiversity. How to incorporate genetic diversity in ongoing monitoring and management of wildlife is still an open question. METHODS: We tested a panel of multiplexed, high-throughput sequenced introns in the small mammal communities of two UNESCO World Heritage Sites on different continents to assess their viability for large-scale monitoring of genetic variability in a spectrum of diverse species. To enhance applicability across other systems, the bioinformatic pipeline for primer design was outlined. RESULTS: The number of loci amplified and amplification evenness decreased as phylogenetic distance increased from the reference taxa, yet several loci were still variable across multiple mammal orders. CONCLUSIONS: Genetic variability found is informative for population genetic analyses and for addressing phylogeographic and phylogenetic questions, illustrated by small mammal examples here.

Comparative assessment of range-wide patterns of genetic diversity and structure with SNPs and microsatellites: A case study with Iberian amphibians.

Reduced representation genome sequencing has popularized the application of single nucleotide polymorphisms (SNPs) to address evolutionary and conservation questions in nonmodel organisms. Patterns of genetic structure and diversity based on SNPs often diverge from those obtained with microsatellites to different degrees, but few studies have explicitly compared their performance under similar sampling regimes in a shared analytical framework. We compared range-wide patterns of genetic structure and diversity in two amphibians endemic to the Iberian Peninsula: Hyla molleri and Pelobates cultripes, based on microsatellite (18 and 14 loci) and SNP (15,412 and 33,140 loci) datasets of comparable sample size and spatial extent. Model-based clustering analyses with STRUCTURE revealed minor differences in genetic structure between marker types, but inconsistent values of the optimal number of populations (K) inferred. SNPs yielded more repeatable and less admixed ancestries with increasing K compared to microsatellites. Genetic diversity was weakly correlated between marker types, with SNPs providing a better representation of southern refugia and of gradients of genetic diversity congruent with the demographic history of both species. Our results suggest that the larger number of loci in a SNP dataset can provide more reliable inferences of patterns of genetic structure and diversity than a typical microsatellite dataset, at least at the spatial and temporal scales investigated.

Endemism and diversity of small mammals along two neighboring Bornean mountains.

Mountains offer replicated units with large biotic and abiotic gradients in a reduced spatial scale. This transforms them into well-suited scenarios to evaluate biogeographic theories. Mountain biogeography is a hot topic of research and many theories have been proposed to describe the changes in biodiversity with elevation. Geometric constraints, which predict the highest diversity to occur in mid-elevations, have been a focal part of this discussion. Despite this, there is no general theory to explain these patterns, probably because of the interaction among different predictors with the local effects of historical factors. We characterize the diversity of small non-volant mammals across the elevational gradient on Mount (Mt.) Kinabalu (4,095 m) and Mt. Tambuyukon (2,579 m), two neighboring mountains in Borneo, Malaysia. We documented a decrease in species richness with elevation which deviates from expectations of the geometric constraints and suggests that spatial factors (e.g., larger diversity in larger areas) are important. The lowland small mammal community was replaced in higher elevations (from above ~1,900 m) with montane communities consisting mainly of high elevation Borneo endemics. The positive correlation we find between elevation and endemism is concordant with a hypothesis that predicts higher endemism with topographical isolation. This supports lineage history and geographic history could be important drivers of species diversity in this region.

Effect of the enzyme and PCR conditions on the quality of high-throughput DNA sequencing results.

Library preparation protocols for high-throughput DNA sequencing (HTS) include amplification steps in which errors can build up. In order to have confidence in the sequencing data, it is important to understand the effects of different Taq polymerases and PCR amplification protocols on the DNA molecules sequenced. We compared thirteen enzymes in three different marker systems: simple, single copy nuclear gene and complex multi-gene family. We also tested a modified PCR protocol, which has been suggested to reduce errors associated with amplification steps. We find that enzyme choice has a large impact on the proportion of correct sequences recovered. The most complex marker systems yielded fewer correct reads, and the proportion of correct reads was greatly affected by the enzyme used. Modified cycling conditions did reduce the number of incorrect sequences obtained in some cases, but enzyme had a much greater impact on the number of correct reads. Thus, the coverage required for the safe identification of genotypes using one of the low quality enzymes could be seven times larger than with more efficient enzymes in a biallelic system with equal amplification of the two alleles. Consequently, enzyme selection for downstream HTS has important consequences, especially in complex genetic systems.

Preservation of RNA and DNA from mammal samples under field conditions.

Ecological and conservation genetics require sampling of organisms in the wild. Appropriate preservation of the collected samples, usually by cryostorage, is key to the quality of the genetic data obtained. Nevertheless, cryopreservation in the field to ensure RNA and DNA stability is not always possible. We compared several nucleic acid preservation solutions appropriate for field sampling and tested them on rat (Rattus rattus) blood, ear and tail tip, liver, brain and muscle. We compared the efficacy of a nucleic acid preservation (NAP) buffer for DNA preservation against 95% ethanol and Longmire buffer, and for RNA preservation against RNAlater (Qiagen) and Longmire buffer, under simulated field conditions. For DNA, the NAP buffer was slightly better than cryopreservation or 95% ethanol, but high molecular weight DNA was preserved in all conditions. The NAP buffer preserved RNA as well as RNAlater. Liver yielded the best RNA and DNA quantity and quality; thus, liver should be the tissue preferentially collected from euthanized animals. We also show that DNA persists in nonpreserved muscle tissue for at least 1 week at ambient temperature, although degradation is noticeable in a matter of hours. When cryopreservation is not possible, the NAP buffer is an economical alternative for RNA preservation at ambient temperature for at least 2 months and DNA preservation for at least 10 months.