Quantitative evaluation of intraspecific genetic diversity in a natural fish population using environmental DNA analysis

Recent advances in environmental DNA (eDNA) analysis using high‐throughput sequencing (HTS) enable evaluation of intraspecific genetic diversity in a population. As the intraspecific genetic diversity provides invaluable information for wildlife conservation and management, there is an increasing demand to apply eDNA analysis to population genetics and the phylogeography by quantitative evaluation of intraspecific diversity. However, quantitative evaluations of intraspecific genetic diversity using eDNA is not straightforward because the number of eDNA sequence reads obtained by HTS may not be an index of the quantity of eDNA. In this study, to quantitatively evaluate genetic diversity using eDNA analysis, we applied a quantitative eDNA metabarcoding method using the internal standard DNAs. We targeted Ayu (Plecoglossus altivelis altivelis) and added internal standard DNAs with known copy numbers to each eDNA sample obtained from three rivers during the library preparation process. The sequence reads of each Ayu haplotype were successfully converted to DNA copy numbers based on the relationship between the copy numbers and sequence reads of the internal standard DNAs. In all rivers, the calculated copy number of each haplotype showed a significant positive correlation with the haplotype frequency estimated by a capture‐based survey. Furthermore, estimates of genetic indicators such as nucleotide diversity based on the eDNA copy numbers were comparable with those estimated based on a capture‐based study. Our results demonstrate that eDNA analysis with internal standard DNAs enables reasonable quantification of intraspecific genetic diversity, and this method could thus be a promising tool in the field of population genetics and phylogeography.     

Helping decision making for reliable and cost‐effective 2b‐RAD sequencing and genotyping analyses in non‐model species

High‐throughput sequencing has revolutionized population and conservation genetics. RAD sequencing methods, such as 2b‐RAD, can be used on species lacking a reference genome. However, transferring protocols across taxa can potentially lead to poor results. We tested two different IIB enzymes (AlfI and CspCI) on two species with different genome sizes (the loggerhead turtle Caretta caretta and the sharpsnout seabream Diplodus puntazzo) to build a set of guidelines to improve 2b‐RAD protocols on non‐model organisms while optimising costs. Good results were obtained even with degraded samples, showing the value of 2b‐RAD in studies with poor DNA quality. However, library quality was found to be a critical parameter on the number of reads and loci obtained for genotyping. Resampling analyses with different number of reads per individual showed a trade‐off between number of loci and number of reads per sample. The resulting accumulation curves can be used as a tool to calculate the number of sequences per individual needed to reach a mean depth ≥20 reads to acquire good genotyping results. Finally, we demonstrated that selective‐base ligation does not affect genomic differentiation between individuals, indicating that this technique can be used in species with large genome sizes to adjust the number of loci to the study scope, to reduce sequencing costs and to maintain suitable sequencing depth for a reliable genotyping without compromising the results. Here, we provide a set of guidelines to improve 2b‐RAD protocols on non‐model organisms with different genome sizes, helping decision‐making for a reliable and cost‐effective genotyping.     

Development and validation of an experimental and theoretical method for the chiral discrimination of dinotefuran

Dinotefuran is a low-cost agrochemical considered a highly toxic product. In this sense, there is a need for its constant environmental, biological, and food control, aiming to ensure its use to humans as well as to preserve biodiversity and ecosystems. In the present work, we developed an experimental and theoretical method for dinotefuran chiral discrimination. According to the main results, the dinotefuran enantioselective separation was efficiently optimized by high-performance liquid chromatography evaluating the influence of different percentage compositions in the mobile phase to improve the resolution of the peaks in the chromatogram. The novelty of this work was the proposition of a reduced molecular model for the chiral selector amylose-Tris-(3,5-dimethylphenylcarbamate) polysaccharide that was able to adequately describe at the molecular level its interaction with the dinotefuran enantiomers. Besides, the thermodynamic and structural parameters obtained via density functional theory calculations pointed out the chiral discrimination as well as the enantiomeric elution order of the analyte studied, confirming the experimental data, thus validating our proposed method. Finally, hydrogen bonds and repulsive interactions played a key role in the discrimination between the diastereomeric complexes, and consequently, for the dinotefuran enantioselective separation.     

Using density surface models to estimate spatio-temporal changes in population densities and trend

Precise measures of population abundance and trend are needed for species conservation; these are most difficult to obtain for rare and rapidly changing populations. We compare uncertainty in densities estimated from spatio–temporal models with that from standard design-based methods. Spatio–temporal models allow us to target priority areas where, and at times when, a population may most benefit. Generalised additive models were fitted to a 31-year time series of point-transect surveys of an endangered Hawaiian forest bird, the Hawai‘i ‘ākepa Loxops coccineus. This allowed us to estimate bird densities over space and time. We used two methods to quantify uncertainty in density estimates from the spatio–temporal model: the delta method (which assumes independence between detection and distribution parameters) and a variance propagation method. With the delta method we observed a 52% decrease in the width of the design-based 95% confidence interval (CI), while we observed a 37% decrease in CI width when propagating the variance. We mapped bird densities as they changed across space and time, allowing managers to evaluate management actions. Integrating detection function modelling with spatio–temporal modelling exploits survey data more efficiently by producing finer-grained abundance estimates than are possible with design-based methods as well as producing more precise abundance estimates. Model-based approaches require switching from making assumptions about the survey design to assumptions about bird distribution. Such a switch warrants consideration. In this case the model-based approach benefits conservation planning through improved management efficiency and reduced costs by taking into account both spatial shifts and temporal changes in population abundance and distribution.     

Cold and hungry: combined effects of low temperature and resource scarcity on an edge-of-range temperate primate,the golden snub-nose monkey

Both biotic and abiotic factors play important roles in influencing ecological distributions and niche limits. Where biotic and abiotic stressors co-occur in space and time, homeostatic systems face a scenario in which stressors can compound to impose a challenge that is greater than the sum of the separate factors. We studied the homeostatic strategies of the golden snub-nosed monkey Rhinopithecus roxellana, a species living in temperate deciduous forests at the edge of the global distribution range for folivorous primates, to cope with the co-occurrence of cold temperatures and resource scarcity during winter. We discovered that in winter the monkeys experience a dietary energy deficit of 101 kJ mbm⁻¹ d⁻¹ compared with calculated needs, despite increased feeding. This is partly offset by behavioral changes (reduced locomotion and increased resting) and reducing skin temperature by an average of 3.2°C through a cutaneous vasoconstriction to decrease heat loss. However, their major strategy is ingesting surplus energy and accumulating fat reserves when food was not limiting during summer and autumn. Their 14% of body mass lost over the winter represented an energy yield of 102 kJ mbm⁻¹ d⁻¹, which closely matched the calculated winter energy deficit of 101 kJ mbm⁻¹ d⁻¹. However, the latter value assumes that all the 75.41 kJ mbm⁻¹ d⁻¹ of protein ingested in winter was available for energy metabolism. This is almost certainly an over-estimate, suggesting that the study population was in negative energy balance over the study period. Our study therefore suggests that despite its suit of integrated homeostatic responses, the confluence of low temperatures and resource limitation during winter places this edge-of-range primate close the threshold of what is energetically viable. It also provides a framework for quantitative models predicting the vulnerability of temperate primates to global change.     

气候变化和经济发展对肾综合征出血热发生的影响

肾综合征出血热(hemorrhagic fever with renal syndrome, HFRS)是一种啮齿动物传播的自然疫源性疾病, 危害严重, 已成为全球重要的公共卫生问题。本研究采用数理统计模型及小波分析方法, 对陕西省西安市鄠邑区1984-2016年HFRS的发生与鼠类、气候和经济因素的关系进行分析, 探讨气候和经济因素对HFRS发生的影响。小波分析结果表明, 该地区的HFRS暴发史可能分为两个时期, 推测每个时期具有不同的主要宿主, 在2002年褐家鼠(Rattus norvegicus)可能取代黑线姬鼠(Apodemus agrarius)成为HFRS疫源地的主要宿主。广义可加模型模拟结果表明, HFRS的发生与1984-2001年黑线姬鼠密度间存在极显著非线性效应(F_2.06,9.02 = 102.415, P < 0.01), 两者间显现为正相关; 与2002-2016年的褐家鼠密度间呈正相关(F_1.67,9.02 = 73.929, P < 0.01); HFRS主要宿主的这种变化可能与当地气候变化和经济发展有关: HFRS的发生与年平均温度存在极显著的非线性效应(F_2.93,9.02 = 12.164, P < 0.01), 两者间呈负相关; 同样, HFRS的发生与上一年的国内生产总值(GDP)也存在显著非线性效应(F_1.70,9.02 = 2.917, P < 0.05), 两者间也呈负相关。结构方程模型通过直接和间接的影响途径证明了这种转移机制, 发现温度对HFRS发生有显著的直接负向影响以及通过褐家鼠的间接正向影响; GDP对HFRS发生有直接的负向影响。本研究表明HFRS的发生与气候变化和经济发展相关, 两者均能影响HFRS的暴发, 该结论有助于今后更好地对HFRS疾病进行预防和控制。