Global patterns of human-wildlife spatial associations and implications for differentiating conservation strategies.

Understanding the global patterns of human and wildlife spatial associations is essential for pragmatic conservation implementation, yet analytical foundations and indicator-based assessments that would further this understanding are lacking. We integrated the global distributions of 30,664 terrestrial vertebrates and human pressures to map human-nature index (HNI) categories that indicate the extent and intensity of human-wildlife interactions. Along the 2 dimensions of biodiversity and human activity, the HNI allowed placement of terrestrial areas worldwide in one of 4 HNI categories: anthropic (human-dominated areas), wildlife-dominated (little human influence and rich in wildlife), co-occurring (substantial presence of humans and wildlife), and harsh-environment (limited presence of humans and wildlife) areas. The HNI varied considerably among taxonomic groups, and the leading driver of HNI was global climate patterns. Co-occurring regions were the most prevalent (35.9%), and wildlife-dominated and anthropic regions encompassed 26.45% and 6.50% of land area, respectively. Our results highlight the necessity for customizing conservation strategies to regions based on human-wildlife spatial associations and the distribution of existing protected area networks. Human activity and biodiversity should be integrated for complementary strategies to support conservation toward ambitious and pragmatic 30×30 goals.

Patrones globales de las asociaciones espaciales entre humanos y fauna y las implicaciones para la diferenciación de las estrategias de conservación Resumen Es esencial entender los patrones globales de asociaciones entre humanos y fauna para la implementación pragmática de la conservación. Aun así, son muy pocos los fundamentos analíticos y las evaluaciones basadas en indicadores que incrementarían este conocimiento. Integramos la distribución global de 30,664 vertebrados terrestres y presiones humanas para mapear las categorías del índice de naturaleza humana (INH) que indican la extensión e intensidad de las interacciones humano­fauna. El INH permitió la colocación de áreas terrestres en todo el mundo en las dos dimensiones de la biodiversidad y las actividades humanas dentro de una de las cuatro categorías del INH: áreas antrópicas (dominadas por humanos), dominadas por fauna (poca influencia humana y rica en fauna), co­ocurrentes (presencia sustancial de humanos y fauna) y de ambiente severo (presencia limitada de humanos y fauna). El INH varió considerablemente entre los taxones, y el factor principal fueron los patrones climáticos mundiales. Las regiones co­ocurrentes fueron las más frecuentes (35.9%) las regiones antrópicas y dominadas por fauna englobaron el 26.45% y 6.50% del área terrestre respectivamente. Nuestros resultados enfatizan la necesidad de personalizar las estrategias de conservación acorde a la región con base en las asociaciones espaciales entre humanos y fauna y la distribución de las redes existentes de áreas protegidas. La actividad humana y la biodiversidad deberían estar integradas para las estrategias complementarias para respaldar a la conservación hacia los objetivos ambiciosos y pragmáticos de 30 para el 30.

Genetics and Evolution of Bird Migration.

Bird migration has long been a subject of fascination for humankind and is a behavior that is both intricate and multifaceted. In recent years, advances in technology, particularly in the fields of genomics and animal tracking, have enabled significant progress in our understanding of this phenomenon. In this review, we provide an overview of the latest advancements in the genetics of bird migration, with a particular focus on genomics, and examine various factors that contribute to the evolution of this behavior, including climate change. Integration of research from the fields of genomics, ecology, and evolution can enhance our comprehension of the complex mechanisms involved in bird migration and inform conservation efforts in a rapidly changing world.

Whole-genome survey reveals extensive variation in genetic diversity and inbreeding levels among peregrine falcon subspecies.

In efforts to prevent extinction, resource managers are often tasked with increasing genetic diversity in a population of concern to prevent inbreeding depression or improve adaptive potential in a changing environment. The assumption that all small populations require measures to increase their genetic diversity may be unwarranted, and limited resources for conservation may be better utilized elsewhere. We test this assumption in a case study focused on the peregrine falcon (Falco peregrinus), a cosmopolitan circumpolar species with 19 named subspecies. We used whole-genome resequencing to generate over two million single nucleotide polymorphisms (SNPs) from multiple individuals of all peregrine falcon subspecies. Our analyses revealed extensive variation among subspecies, with many island-restricted and nonmigratory populations possessing lower overall genomic diversity, elevated inbreeding coefficients (F ROH)-among the highest reported, and extensive runs of homozygosity (ROH) compared to mainland and migratory populations. Similarly, the majority of subspecies that are either nonmigratory or restricted to islands show a much longer history of low effective population size (N e). While mutational load analyses indicated an increased proportion of homozygous-derived deleterious variants (i.e., drift load) among nonmigrant and island populations compared to those that are migrant or reside on the mainland, no significant differences in the proportion of heterozygous deleterious variants (i.e., inbreeding load) was observed. Our results provide evidence that high levels of inbreeding may not be an existential threat for some populations or taxa. Additional factors such as the timing and severity of population declines are important to consider in management decisions about extinction potential.

SARS-CoV-2 RNAs are processed into 22-nt vsRNAs in Vero cells.

The severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) has caused the global pandemic, resulting in great fatalities around the world. Although the antiviral roles of RNA interference (RNAi) have been well studied in plants, nematodes and insects, the antiviral roles of RNAi in mammalians are still debating as RNAi effect is suspected to be suppressed by interferon (IFN) signaling pathways in most cell types. To determine the role of RNAi in mammalian resistance to SARS-CoV-2, we studied the profiling of host small RNAs and SARS-CoV-2 virus-derived small RNAs (vsRNAs) in the early infection stages of Vero cells, an IFN-deficient cell line. We found that host microRNAs (miRNAs) were dysregulated upon SARS-CoV-2 infection, resulting in downregulation of microRNAs playing antiviral functions and upregulation of microRNAs facilitating viral proliferations. Moreover, vsRNA peaked at 22 nt at negative strand but not the positive strand of SARS-CoV-2 and formed successive Dicer-spliced pattern at both strands. Similar characteristics of vsRNAs were observed in IFN-deficient cell lines infected with Sindbis and Zika viruses. Together, these findings indicate that host cell may deploy RNAi pathway to combat SARS-CoV-2 infection in IFN-deficient cells, informing the alternative antiviral strategies to be developed for patients or tissues with IFN deficiency.

Arctic introgression and chromatin regulation facilitated rapid Qinghai-Tibet Plateau colonization by an avian predator.

The Qinghai-Tibet Plateau (QTP), possesses a climate as cold as that of the Arctic, and also presents uniquely low oxygen concentrations and intense ultraviolet (UV) radiation. QTP animals have adapted to these extreme conditions, but whether they obtained genetic variations from the Arctic during cold adaptation, and how genomic mutations in non-coding regions regulate gene expression under hypoxia and intense UV environment, remain largely unknown. Here, we assemble a high-quality saker falcon genome and resequence populations across Eurasia. We identify female-biased hybridization with Arctic gyrfalcons in the last glacial maximum, that endowed eastern sakers with alleles conveying larger body size and changes in fat metabolism, predisposing their QTP cold adaptation. We discover that QTP hypoxia and UV adaptations mainly involve independent changes in non-coding genomic variants. Our study highlights key roles of gene flow from Arctic relatives during QTP hypothermia adaptation, and cis-regulatory elements during hypoxic response and UV protection.

Evolutionary Conservation Genomics Reveals Recent Speciation and Local Adaptation in Threatened Takins.

Incorrect species delimitation will lead to inappropriate conservation decisions, especially for threatened species. The takin (Budorcas taxicolor) is a large artiodactyl endemic to the Himalayan-Hengduan-Qinling Mountains and is well known for its threatened status and peculiar appearance. However, the speciation, intraspecies taxonomy, evolutionary history, and adaptive evolution of this species still remain unclear, which greatly hampers its scientific conservation. Here, we de novo assembled a high-quality chromosome-level genome of takin and resequenced the genomes of 75 wild takins. Phylogenomics revealed that takin was positioned at the root of Caprinae. Population genomics based on the autosome, X chromosome, and Y chromosome SNPs and mitochondrial genomes consistently revealed the existence of two phylogenetic species and recent speciation in takins: the Himalayan takin (B. taxicolor) and the Chinese takin (B. tibetana), with the support of morphological evidence. Two genetically divergent subspecies were identified in both takin species, rejecting three previously proposed taxonomical viewpoints. Furthermore, their distribution boundaries were determined, suggesting that large rivers play important roles in shaping the genetic partition. Compared with the other subspecies, the Qinling subspecies presented the lowest genomic diversity, higher linkage disequilibrium, inbreeding, and genetic load, thus is in urgent need of genetic management and protection. Moreover, coat color gene (PMEL) variation may be responsible for the adaptive coat color difference between the two species following Gloger's rule. Our findings provide novel insights into the recent speciation, local adaptation, scientific conservation of takins, and biogeography of the Himalaya-Hengduan biodiversity hotspot.

RNA-combine: a toolkit for comprehensive analyses on transcriptome data from different sequencing platforms.

BACKGROUND: Understanding the transcriptome has become an essential step towards the full interpretation of the biological function of a cell, a tissue or even an organ. Many tools are available for either processing, analysing transcriptome data, or visualizing analysis results. However, most existing tools are limited to data from a single sequencing platform and only several of them could handle more than one analysis module, which are far from enough to meet the requirements of users, especially those without advanced programming skills. Hence, we still lack an open-source toolkit that enables both bioinformatician and non-bioinformatician users to process and analyze the large transcriptome data from different sequencing platforms and visualize the results. RESULTS: We present a Linux-based toolkit, RNA-combine, to automatically perform the quality assessment, downstream analysis of the transcriptome data generated from different sequencing platforms, including bulk RNA-seq (Illumina platform), single cell RNA-seq (10x Genomics) and Iso-Seq (PacBio) and visualization of the results. Besides, this toolkit is implemented with at least 10 analysis modules more than other toolkits examined in this study. Source codes of RNA-combine are available on GitHub: https://github.com/dongxuemin666/RNA-combine . CONCLUSION: Our results suggest that RNA-combine is a reliable tool for transcriptome data processing and result interpretation for both bioinformaticians and non-bioinformaticians.

Rapid Spread of Severe Fever with Thrombocytopenia Syndrome Virus by Parthenogenetic Asian Longhorned Ticks.

Severe fever with thrombocytopenia syndrome virus (SFTSV) is spreading rapidly in Asia. This virus is transmitted by the Asian longhorned tick (Haemaphysalis longicornis), which has parthenogenetically and sexually reproducing populations. Parthenogenetic populations were found in ≥15 provinces in China and strongly correlated with the distribution of severe fever with thrombocytopenia syndrome cases. However, distribution of these cases was poorly correlated with the distribution of populations of bisexual ticks. Phylogeographic analysis suggested that the parthenogenetic population spread much faster than bisexual population because colonization is independent of sexual reproduction. A higher proportion of parthenogenetic ticks was collected from migratory birds captured at an SFTSV-endemic area, implicating the contribution to the long-range movement of these ticks in China. The SFTSV susceptibility of parthenogenetic females was similar to that of bisexual females under laboratory conditions. These results suggest that parthenogenetic Asian longhorned ticks, probably transported by migratory birds, play a major role in the rapid spread of SFTSV.

Climate-driven flyway changes and memory-based long-distance migration.

Millions of migratory birds occupy seasonally favourable breeding grounds in the Arctic1, but we know little about the formation, maintenance and future of the migration routes of Arctic birds and the genetic determinants of migratory distance. Here we established a continental-scale migration system that used satellite tracking to follow 56 peregrine falcons (Falco peregrinus) from 6 populations that breed in the Eurasian Arctic, and resequenced 35 genomes from 4 of these populations. The breeding populations used five migration routes across Eurasia, which were probably formed by longitudinal and latitudinal shifts in their breeding grounds during the transition from the Last Glacial Maximum to the Holocene epoch. Contemporary environmental divergence between the routes appears to maintain their distinctiveness. We found that the gene ADCY8 is associated with population-level differences in migratory distance. We investigated the regulatory mechanism of this gene, and found that long-term memory was the most likely selective agent for divergence in ADCY8 among the peregrine populations. Global warming is predicted to influence migration strategies and diminish the breeding ranges of peregrine populations of the Eurasian Arctic. Harnessing ecological interactions and evolutionary processes to study climate-driven changes in migration can facilitate the conservation of migratory birds.

Spatial patterns and conservation of genetic and phylogenetic diversity of wildlife in China.

Genetic diversity and phylogenetic diversity reflect the evolutionary potential and history of species, respectively. However, the levels and spatial patterns of genetic and phylogenetic diversity of wildlife at the regional scale have largely remained unclear. Here, we performed meta-analyses of genetic diversity in Chinese terrestrial vertebrates based on three genetic markers and investigated their phylogenetic diversity based on a dated phylogenetic tree of 2461 species. We detected strong positive spatial correlations among mitochondrial DNA-based genetic diversity, phylogenetic diversity, and species richness. Moreover, the terrestrial vertebrates harbored higher genetic and phylogenetic diversity in South China and Southwest China than in other regions. Last, climatic factors (precipitation and temperature) had significant positive effects while altitude and human population density had significant negative impacts on levels of mitochondrial DNA-based genetic diversity in most cases. Our findings will help guide national-level genetic diversity conservation plans and a post-2020 biodiversity conservation framework.

Deeply digging the interaction effect in multiple linear regressions using a fractional-power interaction term.

In multiple regression Y ~ ß0 + ß1X1 + ß2X2 + ß3X1 X2 + É›., the interaction term is quantified as the product of X1 and X2. We developed fractional-power interaction regression (FPIR), using ßX1 M X2 N as the interaction term. The rationale of FPIR is that the slopes of Y-X1 regression along the X2 gradient are modeled using the nonlinear function (Slope = ß1 + ß3MX1 M-1 X2 N), instead of the linear function (Slope = ß1 + ß3X2) that regular regressions normally implement. The ranges of M and N are from -56 to 56 with 550 candidate values, respectively. We applied FPIR using a well-studied dataset, nest sites of the crested ibis (Nipponia nippon).We further tested FPIR by other 4692 regression models. FPIRs have lower AIC values (-302 ± 5003.5) than regular regressions (-168.4 ± 4561.6), and the effect size of AIC values between FPIR and regular regression is 0.07 (95% CI: 0.04-0.10). We also compared FPIR with complex models such as polynomial regression, generalized additive model, and random forest. FPIR is flexible and interpretable, using a minimum number of degrees of freedom to maximize variance explained. We have provided a new R package, interactionFPIR, to estimate the values of M and N, and suggest using FPIR whenever the interaction term is likely to be significant. • Introduced fractional-power interaction regression (FPIR) as Y ~ ß0 + ß1X1 + ß2X2 + ß3X1 M X2 N + É› to replace the current regression model Y ~ ß0 + ß1X1 + ß2X2 + ß3X1 X2 + É›; • Clarified the rationale of FPIR, and compared it with regular regression model, polynomial regression, generalized additive model, and random forest using regression models for 4692 species; • Provided an R package, interactionFPIR, to calculate the values of M and N, and other model parameters.

Genomic analysis of the domestication and post-Spanish conquest evolution of the llama and alpaca.

BACKGROUND: Despite their regional economic importance and being increasingly reared globally, the origins and evolution of the llama and alpaca remain poorly understood. Here we report reference genomes for the llama, and for the guanaco and vicuña (their putative wild progenitors), compare these with the published alpaca genome, and resequence seven individuals of all four species to better understand domestication and introgression between the llama and alpaca. RESULTS: Phylogenomic analysis confirms that the llama was domesticated from the guanaco and the alpaca from the vicuña. Introgression was much higher in the alpaca genome (36%) than the llama (5%) and could be dated close to the time of the Spanish conquest, approximately 500 years ago. Introgression patterns are at their most variable on the X-chromosome of the alpaca, featuring 53 genes known to have deleterious X-linked phenotypes in humans. Strong genome-wide introgression signatures include olfactory receptor complexes into both species, hypertension resistance into alpaca, and fleece/fiber traits into llama. Genomic signatures of domestication in the llama include male reproductive traits, while in alpaca feature fleece characteristics, olfaction-related and hypoxia adaptation traits. Expression analysis of the introgressed region that is syntenic to human HSA4q21, a gene cluster previously associated with hypertension in humans under hypoxic conditions, shows a previously undocumented role for PRDM8 downregulation as a potential transcriptional regulation mechanism, analogous to that previously reported at high altitude for hypoxia-inducible factor 1α. CONCLUSIONS: The unprecedented introgression signatures within both domestic camelid genomes may reflect post-conquest changes in agriculture and the breakdown of traditional management practices.

Convergent genomic signatures of flight loss in birds suggest a switch of main fuel.

Flight loss in birds is as characteristic of the class Aves as flight itself. Although morphological and physiological differences are recognized in flight-degenerate bird species, their contributions to recurrent flight degeneration events across modern birds and underlying genetic mechanisms remain unclear. Here, in an analysis of 295 million nucleotides from 48 bird genomes, we identify two convergent sites causing amino acid changes in ATGLSer321Gly and ACOT7Ala197Val in flight-degenerate birds, which to our knowledge have not previously been implicated in loss of flight. Functional assays suggest that Ser321Gly reduces lipid hydrolytic ability of ATGL, and Ala197Val enhances acyl-CoA hydrolytic activity of ACOT7. Modeling simulations suggest a switch of main energy sources from lipids to carbohydrates in flight-degenerate birds. Our results thus suggest that physiological convergence plays an important role in flight degeneration, and anatomical convergence often invoked may not.