Analysis of Genetic Diversity in Romanian Carpatina Goats Using SNP Genotyping Data.

Animal husbandry is one of man's oldest occupations. It began with the domestication of animals and developed continuously, in parallel with the evolution of human society. The selection and improvement of goats in Romania was not a clearly defined objective until around 1980. In recent years, with the increasing economic value given to goats, breeding programs are becoming established. In Romania, a few goat genetic studies using microsatellites and mtDNA have been carried out; however, a systematic characterization of the country's goat genomic resources remains missing. In this study, we analyzed the genetic variability of Carpatina goats from four distinct geographical areas (northern, north-eastern, eastern and southern Romania), using the Illumina OvineSNP60 (RefSeq ARS1) high-density chip for 67 goats. Heterozygosity values, inbreeding coefficients and effective population size across all autosomes were calculated for those populations that inhabit high- and low-altitude and high- and low-temperature environments. Diversity, as measured by expected heterozygosity (HE), ranged from 0.413 in the group from a low-temperature environment to 0.420 in the group from a high-temperature environment. Within studied groups, the HT (high temperature) goats were the only group with a positive but low average inbreeding coefficient value, which was 0.009. After quality control (QC) analysis, 46,965 SNPs remained for analysis (MAF < 0.01). LD was calculated for each chromosome separately. The Ne has been declining since the time of domestication, having recently reached 123, 125, 185 and 92 for the HA (high altitude), LA (low altitude), HT (high temperature) and LT (low temperature) group, respectively. Our study revealed a low impact of inbreeding in the Carpatina population, and the Ne trend also indicated a steep decline in the last hundred years. These results will contribute to the genetic improvement of the Carpatina breed.

Multiple genomic solutions for local adaptation in two closely related species (sheep and goats) facing the same climatic constraints.

The question of how local adaptation takes place remains a fundamental question in evolutionary biology. The variation of allele frequencies in genes under selection over environmental gradients remains mainly theoretical and its empirical assessment would help understanding how adaptation happens over environmental clines. To bring new insights to this issue we set up a broad framework which aimed to compare the adaptive trajectories over environmental clines in two domesticated mammal species co-distributed in diversified landscapes. We sequenced the genomes of 160 sheep and 161 goats extensively managed along environmental gradients, including temperature, rainfall, seasonality and altitude, to identify genes and biological processes shaping local adaptation. Allele frequencies at putatively adaptive loci were rarely found to vary gradually along environmental gradients, but rather displayed a discontinuous shift at the extremities of environmental clines. Of the 430 candidate adaptive genes identified, only 6 were orthologous between sheep and goats and those responded differently to environmental pressures, suggesting different putative mechanisms involved in local adaptation in these two closely related species. Interestingly, the genomes of the 2 species were impacted differently by the environment, genes related to signatures of selection were most related to altitude, slope and rainfall seasonality for sheep, and summer temperature and spring rainfall for goats. The diversity of candidate adaptive pathways may result from a high number of biological functions involved in the adaptations to multiple eco-climatic gradients, and a differential role of climatic drivers on the two species, despite their co-distribution along the same environmental gradients. This study describes empirical examples of clinal variation in putatively adaptive alleles with different patterns in allele frequency distributions over continuous environmental gradients, thus showing the diversity of genetic responses in adaptive landscapes and opening new horizons for understanding genomics of adaptation in mammalian species and beyond.

Comparative phylogeography reveals dissimilar genetic differentiation patterns in two sympatric amphibian species.

Global climate change is expected to have a profound effect on species distribution. Due to the temperature constraints, some narrow niche species could shift their narrow range to higher altitudes or latitudes. In this study, we explored the correlation between species traits, genetic structure, and geographical range size. More specifically, we analyzed how these variables are affected by differences in fundamental niche breadth or dispersal ability in the members of two sympatrically distributed stream-dwelling amphibian species (frog, Quasipaa yei; salamander, Pachyhynobius shangchengensis), in Dabie Mountains, East China. Both species showed relatively high genetic diversity in most geographical populations and similar genetic diversity patterns (JTX, low; BYM, high) correlation with habitat changes and population demography. Multiple clustering analyses were used to disclose differentiation among the geographical populations of these two amphibian species. Q. yei disclosed the relatively shallow genetic differentiation, while P. shangchengensis showed an opposite pattern. Under different historical climatic conditions, all ecological niche modeling disclosed a larger suitable habitat area for Q. yei than for P. shangchengensis; these results indicated a wider environment tolerance or wider niche width of Q. yei than P. shangchengensis. Our findings suggest that the synergistic effects of environmental niche variation and dispersal ability may help shape genetic structure across geographical topology, particularly for species with extremely narrow distribution.

The fast-food effect: costs of being a generalist in a human-dominated landscape.

Agricultural expansion in Southeast Asia has converted most natural landscapes into mosaics of forest interspersed with plantations, dominated by the presence of generalist species that benefit from resource predictability. Dietary shifts, however, can result in metabolic alterations and the exposure of new parasites that can impact animal fitness and population survival. Our study focuses on the Asian water monitor lizard (Varanus salvator), one of the largest predators in the Asian wetlands, as a model species to understand the health consequences of living in a human-dominated landscape in Sabah, Malaysian Borneo. We evaluated the effects of dietary diversity on the metabolism of monitor lizards and the impact on the composition of their parasite communities in an oil palm-dominated landscape. Our results showed that (1) rodent-dominated diets were associated with high levels of lipids, proteins and electrolytes, akin to a fast-food-based diet of little representativeness of the full nutritional requirements, but highly available, and (2) lizards feeding on diverse diets hosted more diverse parasite communities, however, at overall lower parasite prevalence. Furthermore, we observed that the effect of diet on lipid concentration differed depending on the size of individual home ranges, suggesting that sedentarism plays an important role in the accumulation of cholesterol and triglycerides. Parasite communities were also affected by a homogeneous dietary behaviour, as well as by habitat type. Dietary diversity had a negative effect on both parasite richness and prevalence in plantations, but not in forested areas. Our study indicates that human-dominated landscapes can pose a negative effect on generalist species and hints to the unforeseen health consequences for more vulnerable taxa using the same landscapes. Thus, it highlights the potential role of such a widely distributed generalist as model species to monitor physiological effects in the ecosystem in an oil palm-dominated landscape.

Herbivorous dietary selection shown by hawfinch (Coccothraustes coccothraustes) within mixed woodland habitats.

Knowledge of diet and dietary selectivity is vital, especially for the conservation of declining species. Accurately obtaining this information, however, is difficult, especially if the study species feeds on a wide range of food items within heterogeneous and inaccessible environments, such as the tree canopy. Hawfinches (Coccothraustes coccothraustes), like many woodland birds, are declining for reasons that are unclear. We investigated the possible role that dietary selection may have in these declines in the UK. Here, we used a combination of high-throughput sequencing of 261 hawfinch faecal samples assessed against tree occurrence data from quadrats sampled in three hawfinch population strongholds in the UK to test for evidence of selective foraging. This revealed that hawfinches show selective feeding and consume certain tree genera disproportionally to availability. Positive selection was shown for beech (Fagus), cherry (Prunus), hornbeam (Carpinus), maples (Acer) and oak (Quercus), while Hawfinch avoided ash (Fraxinus), birch (Betula), chestnut (Castanea), fir (Abies), hazel (Corylus), rowan (Sorbus) and lime (Tilia). This approach provided detailed information on hawfinch dietary choice and may be used to predict the effects of changing food resources on other declining passerines populations in the future.

Multi-marker DNA metabarcoding reveals spatial and sexual variation in the diet of a scarce woodland bird.

Avian diet can be affected by site-specific variables, such as habitat, as well as intrinsic factors such as sex. This can lead to dietary niche separation, which reduces competition between individuals, as well as impacting how well avian species can adapt to environmental variation. Estimating dietary niche separation is challenging, due largely to difficulties in accurately identifying food taxa consumed. Consequently, there is limited knowledge of the diets of woodland bird species, many of which are undergoing serious population declines. Here, we show the effectiveness of multi-marker fecal metabarcoding to provide in-depth dietary analysis of a declining passerine in the UK, the Hawfinch (Coccothraustes coccothraustes). We collected fecal samples from (n = 262) UK Hawfinches prior to, and during, the breeding seasons in 2016-2019. We detected 49 and 90 plant and invertebrate taxa, respectively. We found Hawfinch diet varied spatially, as well as between sexes, indicating broad dietary plasticity and the ability of Hawfinches to utilize multiple resources within their foraging environments.

Reconstruction of the Major Maternal and Paternal Lineages in the Feral New Zealand Kaimanawa Horses.

New Zealand has the fourth largest feral horse population in the world. The Kaimanawas (KHs) are feral horses descended from various domestic horse breeds released into the Kaimanawa ranges in the 19th and 20th centuries. Over time, the population size has fluctuated dramatically due to hunting, large-scale farming and forestry. Currently, the herd is managed by an annual round-up, limiting the number to 300 individuals to protect the native ecosystem. Here, we genotyped 96 KHs for uniparental markers (mitochondrial DNA, Y-chromosome) and assessed their genetic similarity with respect to other domestic horses. We show that at least six maternal and six paternal lineages contributed unequally to the KH gene pool, and today's KH population possibly represents two sub-populations. Our results indicate that three horse breeds, namely Welsh ponies, Thoroughbreds and Arabian horses had a major influence in the genetic-makeup of the extant KH population. We show that mitochondrial genetic diversity in KHs (π = 0.00687 ± 0.00355) is closer to that of the Sable Island horses (π = 0.0034 ± 0.00301), and less than other feral horse populations around the world. Our current findings, combined with ongoing genomic research, will provide insight into the population-specific genetic variation and inbreeding among KHs. This will largely advance equine research and improve the management of future breeding programs of these treasured New Zealand horse.

Molecular insights into high-altitude adaption and acclimatisation of Aporrectodea caliginosa.

Here, we explore the high-altitude adaptions and acclimatisation of Aporrectodea caliginosa Population diversity is assessed through mitochondrial barcoding, identifying closely related populations across the island of Pico (Azores). We present the first megabase N50 assembly size (1.2 Mbp) genome for A. caliginosa High- and low-altitude populations were exposed experimentally to a range of oxygen and temperature conditions, simulating altitudinal conditions, and the transcriptomic responses explored. SNP densities are assessed to identify signatures of selective pressure and their link to differentially expressed genes. The high-altitude A. caliginosa population had lower differential expression and fewer co-expressed genes between conditions, indicating a more condition-refined epigenetic response. Genes identified as under adaptive pressure through Fst and nucleotide diversity in the high-altitude population clustered around the differentially expressed an upstream environmental response control gene, HMGB1. The high-altitude population of A. caliginosa indicated adaption and acclimatisation to high-altitude conditions and suggested resilience to extreme weather events. This mechanistic understanding could help offer a strategy in further identifying other species capable of maintaining soil fertility in extreme environments.

A population genetic analysis of the Critically Endangered Madagascar big-headed turtle, Erymnochelys madagascariensis across captive and wild populations.

Erymnochelys madagascariensis is a Critically Endangered turtle endemic to Madagascar. Anthropogenic activity has depleted the wild population by 70% in the last century, and effective conservation management is essential to ensuring its persistence. Captive breeding was implemented to augment depleted populations in the southern part of Ankarafantsika National Park (ANP), when no genetic data were available for E. madagascariensis. It is unknown how much of the natural population's diversity is encapsulated in captivity. We used eight microsatellite loci and fragments of two mitochondrial genes to identify the genetic structure of E. madagascariensis in the wild. Captive bred turtles were compared with wild populations in order to assess the representativeness of this ex situ conservation strategy for ANP. Six microsatellite clusters, ten cytochrome b, and nine COI haplotypes were identified across wild populations, with high genetic divergence found between populations in two groups of watersheds. Captive bred individuals represent three out of six sampled microsatellite clusters found in the wild and just one mitochondrial haplotype, possibly due to genetic drift. To improve genetic representation, the strategy of frequent interchange between captive and wild breeders within ANP should be revitalised and, as originally planned, hatchlings or juveniles should not be released beyond ANP.

The effect of oil palm-dominated landscapes on the home range and distribution of a generalist species, the Asian water monitor.

The Asian water monitor lizard, Varanus salvator, is one of the largest predators in Southeast Asia which persists in human-dominated landscapes and, as such, is a suitable model to understand the behavioral plasticity of generalists in anthropogenic landscapes. We used Local Convex Hull with adaptive algorithm to estimate the home range size of 14 GPS-tagged individuals, followed by a MAXENT approach and community prey composition to understand the habitat preferences within the landscape. We estimated larger home ranges in forest than in oil palm plantations, as well as a larger diversity and abundance of mammals. Core home ranges were always linked to water bodies. However, the use of underproductive oil palm, freshwater swamp forest, and degraded forest by monitor lizards were higher than other kind of vegetation. This suitable habitat is proportionally larger in forest (73.7%) than in oil palm plantations (39.6%). Generalized estimation equation models showed that, while full home range size was negatively associated with the abundance of mammals, core areas depicted a positive association with mammal abundance, as well as with the proportion of suitable habitat within the home range. Besides having smaller home ranges in oil palm plantations, our findings suggest that limited suitable habitat availability forces the Asian water monitor lizard's population to establish only one or very few core areas. Contrastingly, under the protection of forest, they have more core areas, widely dispersed within larger home ranges. We conclude that regardless the plasticity of the species, human-dominated landscapes are altering natural patterns of home range establishment in the monitor lizard's population, creating a potential ecological trap where conditions may not remain favorable for them in the long run. A deeper understanding of the ecological implications on the species and the prey community is advisable.

The complete mitochondrial genome of rare and Critically Endangered Anilany helenae (Microhylidae) of Madagascar.

Anilany helenae is a Critically Endangered frog native to the central highlands of Madagascar. Due to ongoing habitat loss of its known range, this species' population is considered declining, while little is known about its ecology, behavior, and taxonomy. Within the context of developing tools that can aid the conservation of Madagascar's amphibian fauna, and add to the continued understanding of their taxonomy, we assembled its complete mitochondrial genome (Genbank Accession number MZ751042). This contributes the first complete mitochondrial genome of a microhylid from Madagascar, despite there being over 100 species in the Cophylinae subfamily alone. Anilany helenae's circular mitochondrial genome is 17,519 bp long, contains 37 genes, and exhibits differences in gene arrangement compared with other microhylids, including the placement of protein coding genes nad1 and nad2. A phylogeny of the 13 protein coding genes of the few Madagascan anuran mitogenomes available, along with species from Africa and East Asia, places A. helenae along with the New Guinean Mantophryne lateralis in a basal position with respect to the other microhylids in the tree.

An unexpected new red-bellied Stumpffia (Microhylidae) from forest fragments in central Madagascar highlights remaining cryptic diversity.

The Madagascan endemic subfamily Cophylinae in the family Microhylidae, is an example of a taxonomic group for which much is still to be discovered. Indeed, the cophyline frogs present a large portion of Madagascar's cryptic and microendemic amphibian diversity, yet they remain understudied. A new red-bellied species of the microhylid frog genus Stumpffia is described from the central plateau of Madagascar. Visual encounter surveys in Ambohitantely and Anjozorobe in 2019 and 2020 identified this previously unknown Stumpffia species, which closely resembles Stumpffiakibomena known from Andasibe in the east. Stumpffialynnae sp. nov. adds another species to the red-bellied species complex, differing from S.kibomena by genetic differentiation in the mitochondrial 16S rRNA gene (3.6-3.9%) and distinct nuclear RAG1 haplotypes, as well as strongly by its advertisement call. The new species is known from across Ambohitantely Special Reserve and Anjozorobe Angavo protected area, but is known only from one complete specimen and eight individual tissue samples. Based on the rarity of the species, the small number of locations in which it has been found, and its disappearing forest habitat, its IUCN Red List classification is suggested as "Endangered". This species is the first Stumpffia described from Madagascar's central plateau, highlighting the importance of conserving the remnant forest fragments in this area and the ongoing need to survey and protect this threatened habitat type.

MEDI: Macronutrient Extraction and Determination from invertebrates, a rapid, cheap and streamlined protocol.

Macronutrients, comprising carbohydrates, proteins and lipids, underpin many ecological processes, but their quantification in ecological studies is often inaccurate and laborious, requiring large investments of time and bulk samples, which make individual-level studies impossible. This study presents Macronutrient Extraction and Determination from Invertebrates (MEDI), a protocol for the direct, rapid and relatively low-cost determination of macronutrient content from single small macroinvertebrates.Macronutrients were extracted by a sequential process of soaking in 1:12 chloroform:methanol solution to remove lipid and then solubilising tissue in 0.1 M NaOH. Proteins, carbohydrates and lipids were determined by colorimetric assays from the same individual specimens.The limits of detection of MEDI with the equipment and conditions used were 0.067, 0.065 and 0.006 mg/ml for proteins, carbohydrates and lipids respectively. Adjusting the volume of reagents used for extraction and determination can broaden the range of concentrations that can be detected. MEDI successfully identified taxonomic differences in macronutrient content between five insect species.Macronutrient Extraction and Determination from Invertebrates can directly and rapidly determine macronutrient content in tiny (dry mass ~3 mg) and much larger individual invertebrates. Using MEDI, the total macronutrient content of over 50 macroinvertebrates can be determined within around 3 days of collection at a cost of ~$1.35 per sample.

Spatial dynamics of Chinese Muntjac related to past and future climate fluctuations.

Climate fluctuations in the past and in the future are likely to result in population expansions, shifts, or the contraction of the ecological niche of many species, and potentially leading to the changes in their geographical distributions. Prediction of suitable habitats has been developed as a useful tool for the assessment of habitat suitability and resource conservation to protect wildlife. Here, we model the ancestral demographic history of the extant modern Chinese Muntjac Muntiacus reevesi populations using approximate Bayesian computation (ABC) and used the maximum entropy model to simulate the past and predict the future spatial dynamics of the species under climate oscillations. Our results indicated that the suitable habitats for the M. reevesi shifted to the Southeast and contracted during the Last Glacial Maximum, whereas they covered a broader and more northern position in the Middle Holocene. The ABC analyses revealed that the modern M. reevesi populations diverged in the Middle Holocene coinciding with the significant contraction of the highly suitable habitat areas. Furthermore, our predictions suggest that the potentially suitable environment distribution for the species will expand under all future climate scenarios. These results indicated that the M. reevesi diverged in the recent time after the glacial period and simultaneously as its habitat's expanded in the Middle Holocene. Furthermore, the past and future climate fluctuation triggered the change of Chinese muntjac spatial distribution, which has great influence on the Chinese muntjac's population demographic history.

The critical role of natural forest as refugium for generalist species in oil palm-dominated landscapes.

In Borneo, oil palm plantations have replaced much of natural resources, where generalist species tend to be the principal beneficiaries, due to the abundant food provided by oil palm plantations. Here, we analyse the distribution of the Asian water monitor lizard (Varanus salvator) population within an oil palm-dominated landscape in the Kinabatangan floodplain, Malaysian Borneo. By using mark-recapture methods we estimated its population size, survival, and growth in forest and plantation habitats. We compared body measurements (i.e. body weight and body length) of individuals living in forest and oil palm habitats as proxy for the population's health status, and used general least squares estimation models to evaluate its response to highly fragmented landscapes in the absence of intensive hunting pressures. Contrary to previous studies, the abundance of lizards was higher in the forest than in oil palm plantations. Recruitment rates were also higher in the forest, suggesting that these areas may function as a source of new individuals into the landscape. While there were no morphometric differences among plantation sites, we found significant differences among forested areas, where larger lizards were found inhabiting forest adjacent to oil palm plantations. Although abundant in food resources, the limited availability of refugia in oil palm plantations may intensify intra-specific encounters and competition, altering the body size distribution in plantation populations, contrary to what happens in the forest. We conclude that large patches of forest, around and within oil palm plantations, are essential for the dynamics of the monitor lizard population in the Kinabatangan floodplain, as well as a potential source of individuals to the landscape. We recommend assessing this effect in other generalist species, as well as the impact on the prey communities, especially to reinforce the establishment of buffer zones and corridors as a conservation strategy within plantations.

Local Ancestry to Identify Selection in Response to Trypanosome Infection in Baoulé x Zebu Crossbred Cattle in Burkina Faso.

The genomes of crossbred (admixed) individuals are a mosaic of ancestral haplotypes formed by recombination in each generation. The proportion of these ancestral haplotypes in certain genomic regions can be responsible for either susceptibility or tolerance against pathogens, and for performances in production traits. Using a medium-density genomic marker panel from the Illumina Bovine SNP50 BeadChip, we estimated individual admixture proportions for Baoulé x Zebu crossbred cattle in Burkina Faso, which were tested for trypanosome infection by direct ELISA from blood samples. Furthermore, we calculated local ancestry deviation from average for each SNP across 29 autosomes to identify potential regions under selection in the trypanotolerant Baoulé cattle and their crossbreds. We identified significant deviation from the local average ancestry (above 5 and 10% genome-wide thresholds) on chromosomes 8 and 19 in the positive animals, while the negative ones showed higher deviation on chromosomes 6, 19, 21, and 22. Some candidate genes on chromosome 6 (PDGFRA) and chromosome 19 (CDC6) have been found associated to trypanotolerance in West African taurines. Screening for F ST outliers in trypanosome positive/negative animals we detected seven variants putatively under selection. Finally, we identified a minimum set of highly ancestry informative markers for routine admixture testing. The results of this study contribute to a better understanding of the genetic basis of trypanotolerance in Baoulé cattle and their crossbreeds. Furthermore, we provide a small informative marker set to monitor admixture in this valuable indigenous breed. As such, our results are important for conserving the genetic uniqueness and trypanotolerance of Baoulé cattle, as well as for the improvement of Baoulé and Zebu crossbreds in specific community-based breeding programs.

Detection of selection signatures in the genome of a farmed population of anadromous rainbow trout (Oncorhynchus mykiss).

Domestication processes and artificial selection are likely to leave signatures that can be detected at a molecular level in farmed rainbow trout (Oncorhynchus mykiss). These signatures of selection are genomic regions that contain functional genetic variants conferring a higher fitness to their bearers. We genotyped 749 rainbow trout from a commercial population using a rainbow trout Axiom 57 K SNP array panel and identified putative genomic regions under selection using the pcadapt, Composite Likelihood Ratio (CLR) and Integrated Haplotype Score (iHS) methods. After applying quality-control pipelines and statistical analyses, we detected 12, 96 and 16 SNPs putatively under selection, associated with 96, 781 and 115 candidate genes, respectively. Several of these candidate genes were associated with growth, early development, reproduction, behavior and immune system traits. In addition, some of the SNPs were found in interesting regions located in autosomal inversions on Omy05 and Omy20. These findings could represent a genome-wide map of selection signatures in farmed rainbow trout and could be important in explaining domestication and selection for genetic traits of commercial interest.

Innate and Adaptive Immune Genes Associated with MERS-CoV Infection in Dromedaries.

The recent SARS-CoV-2 pandemic has refocused attention to the betacoronaviruses, only eight years after the emergence of another zoonotic betacoronavirus, the Middle East respiratory syndrome coronavirus (MERS-CoV). While the wild source of SARS-CoV-2 may be disputed, for MERS-CoV, dromedaries are considered as source of zoonotic human infections. Testing 100 immune-response genes in 121 dromedaries from United Arab Emirates (UAE) for potential association with present MERS-CoV infection, we identified candidate genes with important functions in the adaptive, MHC-class I (HLA-A-24-like) and II (HLA-DPB1-like), and innate immune response (PTPN4, MAGOHB), and in cilia coating the respiratory tract (DNAH7). Some of these genes previously have been associated with viral replication in SARS-CoV-1/-2 in humans, others have an important role in the movement of bronchial cilia. These results suggest similar host genetic pathways associated with these betacoronaviruses, although further work is required to better understand the MERS-CoV disease dynamics in both dromedaries and humans.

The complete mitogenome of the Mountain chicken frog, Leptodactylus fallax.

The mountain chicken frog (Leptodactylus fallax) is a critically endangered frog native to the Caribbean islands of Dominica and Montserrat. Over the past 25 years their populations have declined by over 85%, largely due to a chytridiomycosis outbreak that nearly wiped out the Montserratian population. Within the context of developing tools that can aid in the conservation of the mountain chicken frog, we assembled its complete mitochondrial genome, contributing the first complete mitogenome of the genus Leptodactylus (Genbank Accession number MW260634). The circular genome is 18,669 bp long and contains 37 genes. A phylogenetic analysis reveals that L. fallax forms a clade with Leptodactylus melanonotus, highlighting the close relationship of Leptodactylus spp. relative to other species from the superfamily Hyloidea included in the analysis.