Genetic diversity and signatures of selection in Icelandic horses and Exmoor ponies.

BACKGROUND: The Icelandic horse and Exmoor pony are ancient, native breeds, adapted to harsh environmental conditions and they have both undergone severe historic bottlenecks. However, in modern days, the selection pressures on these breeds differ substantially. The aim of this study was to assess genetic diversity in both breeds through expected (HE) and observed heterozygosity (HO) and effective population size (Ne). Furthermore, we aimed to identify runs of homozygosity (ROH) to estimate and compare genomic inbreeding and signatures of selection in the breeds. RESULTS: HO was estimated at 0.34 and 0.33 in the Icelandic horse and Exmoor pony, respectively, aligning closely with HE of 0.34 for both breeds. Based on genomic data, the Ne for the last generation was calculated to be 125 individuals for Icelandic horses and 42 for Exmoor ponies. Genomic inbreeding coefficient (FROH) ranged from 0.08 to 0.20 for the Icelandic horse and 0.12 to 0.27 for the Exmoor pony, with the majority of inbreeding attributed to short ROHs in both breeds. Several ROH islands associated with performance were identified in the Icelandic horse, featuring target genes such as DMRT3, DOCK8, EDNRB, SLAIN1, and NEURL1. Shared ROH islands between both breeds were linked to metabolic processes (FOXO1), body size, and the immune system (CYRIB), while private ROH islands in Exmoor ponies were associated with coat colours (ASIP, TBX3, OCA2), immune system (LYG1, LYG2), and fertility (TEX14, SPO11, ADAM20). CONCLUSIONS: Evaluations of genetic diversity and inbreeding reveal insights into the evolutionary trajectories of both breeds, highlighting the consequences of population bottlenecks. While the genetic diversity in the Icelandic horse is acceptable, a critically low genetic diversity was estimated for the Exmoor pony, which requires further validation. Identified signatures of selection highlight the differences in the use of the two breeds as well as their adaptive trait similarities. The results provide insight into genomic regions under selection pressure in a gaited performance horse breed and various adaptive traits in small-sized native horse breeds. This understanding contributes to preserving genetic diversity and population health in these equine populations.

Genetic divergence and truncation and simultaneous selection in inbred families (S1) of elephant grass for bioenergetic purposes via mixed models.

The State University of North Fluminense Darcy Ribeiro (UENF) has been developing for fifteen years a breeding program that aims at the development of new cultivars of elephant grass due to its high potential and the low availability of cultivars developed by genetic breeding programs that meet the needs of producers in the State of Rio de Janeiro. In this sense, inbred families were also obtained as a way of fixing potential alleles for traits related to production, as the inbreeding process apparently does not strongly affect elephant grass in aspects related to inbreeding depression. This study aimed to estimate genetic diversity, variance components and prediction of genotypic values in 11 (S1) elephant grass families, and perform the truncation and simultaneous selection of traits using the selection index, by mixed models. The experimental design consisted of randomized blocks with 11 (S1) families, three replications, and six plants per plot. For variables dry matter production, percentage of dry matter, plant height, stem diameter, number of tillers and leaf blade width, was performed the estimation of genetic parameters and selection of the best genotypes based selection index using mixed model. The descriptors were subjected to correlation analysis, distance matrices were generated by the Mahalanobis method, and individuals were grouped by the UPGMA method. In the selection via mixed models (REML/BLUP), families 6, 11, 8, 1, 3, 7, and 9 contributed most of the genotypes selected for the evaluated traits, indicating their high potential to generate superior genotype. The selection indices via mixed models indicated that the multiplicative index presented a greater selection gain.The phenotypic correlations showed the possibility of performing an indirect selection from six evaluated traits.The genotypes were separated into 18 groups by the Mahalanobis distance, allowing the observation of a wide genetic diversity. The most divergent and productive genotypes were self-fertilized to obtain the second generation (S2), continuing the development program.

Genome-wide detection of runs of homozygosity and heterozygosity in Tunchang pigs.

Tunchang pigs, mainly distributed throughout Hainan Province of China, are well-known for their superior meat quality, crude feed tolerance, and adaptability to high temperatures and humidity. Runs of homozygosity (ROH) can provide valuable information about the inbreeding coefficient in individuals and selection signals that may reveal candidate genes associated with key functional traits. Runs of heterozygosity (ROHet) are commonly associated with balance selection, which can help us understand the adaptive evolutionary history of domestic animals. In this study, we investigated ROHs and ROHets in 88 Tunchang pigs. We also compared the estimates of inbreeding coefficients in individuals calculated based on four methods. In summary, we detected a total of 16 ROH islands in our study, and 100 genes were found within ROH regions. These genes were correlated with economically important traits such as reproduction (e.g., SERPIND1, HIRA), meat quality (e.g., PI4KA, TBX1), immunity (e.g., ESS2, RANBP1), adaption to heat stress (TXNRD2 and DGCR8), and crude food tolerance (TRPM6). Moreover, we discovered 18 ROHet islands harbouring genes associated with reproduction (e.g., ARHGEF12, BMPR2), immune system (e.g., BRD4, DNMT3B). These findings may help us design effective breeding and conservation strategies for this unique breed.

Clinging on the brink: Whole genomes reveal human-induced population declines and severe inbreeding in the Critically Endangered Emirati Leaf-toed Gecko (Asaccus caudivolvulus).

Human-mediated habitat destruction has had a profound impact on increased species extinction rates and population declines worldwide. The coastal development in the United Arab Emirates (UAE) over the last two decades, serves as an example of how habitat transformation can alter the landscape of a country in just a few years. Here, we study the genomic implications of habitat transformation in the Critically Endangered Emirati Leaf-toed Gecko (Asaccus caudivolvulus), the only endemic vertebrate of the UAE. We generate a high-quality reference genome for this gecko, representing the first reference genome for the family Phyllodactylidae, and produce whole-genome resequencing data for 23 specimens from 10 different species of leaf-toed geckos. Our results show that A. caudivolvulus has consistently lower genetic diversity than any other Arabian species of Asaccus, suggesting a history of ancient population declines. However, high levels of recent inbreeding are recorded among populations in heavily developed areas, with a more than 50% increase in long runs of homozygosity within a 9-year period. Moreover, results suggest that this species does not effectively purge deleterious mutations, hence making it more vulnerable to future stochastic threats. Overall, results show that A. caudivolvulus is in urgent need of protection, and habitat preservation must be warranted to ensure the species' survival.

Ecological genetics of isolated loach populations indicate compromised adaptive potential.

Many endangered species live in fragmented and isolated populations with low genetic variability, signs of inbreeding, and small effective population sizes - all features elevating their extinction risk. The flat-headed loach (Oreonectes platycephalus), a small noemacheilid fish, is widely across southern China, but only in the headwaters of hillstreams; as a result, they are spatially isolated from conspecific populations. We surveyed single nucleotide polymorphisms in 16 Hong Kong populations of O. platycephalus to determine whether loach populations from different streams were genetically isolated from each other, showed low levels of genetic diversity, signs of inbreeding, and had small contemporary effective population sizes. Estimates of average observed heterozygosity (HO = 0.0473), average weighted nucleotide diversity (πw = 0.0546) and contemporary effective population sizes (Ne = 10.2 ~ 129.8) were very low, and several populations showed clear signs of inbreeding as judged from relatedness estimates. The degree of genetic differentiation among populations was very high (average FST = 0.668), even over short geographic distances (<1.5 km), with clear patterns of isolation by distance. These results suggest that Hong Kong populations of O. platycephalus have experienced strong genetic drift and loss of genetic variability because sea-level rise after the last glaciation reduced connectedness among paleodrainages, isolating populations in headwaters. All this, together with the fact that the levels of genetic diversity and contemporary effective population sizes within O. platycephalus populations are lower than most other freshwater fishes, suggests that they face high local extinction risk and have limited capacity for future adaptation.

Genetic diversity analysis of Inner Mongolia cashmere goats (Erlangshan subtype) based on whole genome re-sequencing.

BACKGROUND: Inner Mongolia cashmere goat (IMCG), renowned for its superior cashmere quality, is a Chinese indigenous goat breed that has been developed through natural and artificial selection over a long period. However, recently, the genetic resources of IMCGs have been significantly threatened by the introduction of cosmopolitan goat breeds and the absence of adequate breed protection systems. RESULTS: In order to assess the conservation effectiveness of IMCGs and efficiently preserve and utilize the purebred germplasm resources, this study analyzed the genetic diversity, kinship, family structure, and inbreeding of IMCGs utilizing resequencing data from 225 randomly selected individuals analyzed using the Plink (v.1.90), GCTA (v.1.94.1), and R (v.4.2.1) software. A total of 12,700,178 high-quality SNPs were selected through quality control from 34,248,064 SNP sites obtained from 225 individuals. The average minor allele frequency (MAF), polymorphic information content (PIC), and Shannon information index (SHI) were 0.253, 0.284, and 0.530, respectively. The average observed heterozygosity (Ho) and the average expected heterozygosity (He) were 0.355 and 0.351, respectively. The analysis of the identity by state distance matrix and genomic relationship matrix has shown that most individuals' genetic distance and genetic relationship are far away, and the inbreeding coefficient is low. The family structure analysis identified 10 families among the 23 rams. A total of 14,109 runs of homozygosity (ROH) were identified in the 225 individuals, with an average ROH length of 1014.547 kb. The average inbreeding coefficient, calculated from ROH, was 0.026 for the overall population and 0.027 specifically among the 23 rams, indicating a low level of inbreeding within the conserved population. CONCLUSIONS: The IMCGs exhibited moderate polymorphism and a low level of kinship with inbreeding occurring among a limited number of individuals. Simultaneously, it is necessary to prevent the loss of bloodline to guarantee the perpetuation of the IMCGs' germplasm resources.

FnR: R package for computing inbreeding and numerator relationship coefficients.

BACKGROUND: Inbreeding and relationship coefficients are essential for conservation and breeding programs. Whether dealing with a small conserved population or a large commercial population, monitoring the inbreeding rate and designing mating plans that minimize the inbreeding rate and maximize the effective population size is important. Free, open-source, and efficient software may greatly contribute to conservation and breeding programs and help students and researchers. Efficient methods exist for calculating inbreeding coefficients. Therefore, an efficient way of calculating the numerator relationship coefficients is via the inbreeding coefficients. i.e., the relationship coefficient between parents is twice the inbreeding coefficient of their progeny. A dummy progeny is introduced where no progeny exists for a pair of individuals. Calculating inbreeding coefficients is very fast, and finding whether a pair of individuals has a progeny and picking one from multiple progenies is computationally more demanding. Therefore, the R package introduces a dummy progeny for any pair of individuals whose relationship coefficient is of interest, whether they have a progeny or not. RESULTS: Runtime and peak memory usage were benchmarked for calculating relationship coefficients between two sets of 250 and 800 animals (200,000 dummy progenies) from a pedigree of 2,721,252 animals. The program performed efficiently (200,000 relationship coefficients, which involved calculating 2,721,252 + 200,000 inbreeding coefficients) within 3:45 (mm:ss). Providing the inbreeding coefficients (for real animals), the runtime was reduced to 1:08. Furthermore, providing the diagonal elements of D in A = TDT ' (d), the runtime was reduced to 54s. All the analyses were performed on a machine with a total memory size of 1 GB. CONCLUSIONS: The R package FnR is free and open-source software with implications in conservation and breeding programs. It proved to be time and memory efficient for large populations and many dummy progenies. Calculation of inbreeding coefficients can be resumed for new animals in the pedigree. Thus, saving the latest inbreeding coefficient estimates is recommended. Calculation of d coefficients (from scratch) was very fast, and there was limited value in storing those for future use.

Assessing the population structure and genetic variability of Kenyan native goats under extensive production system.

Indigenous goats are important to many livelihoods. Despite this, they are subjected to indiscriminate crossbreeding. This affects their genetic variability which is needed to survive in current regime of climate change. The study assessed population structure and genetic diversity of Galla and Small East African goats (SEA) using pedigree information. A total of 7384 animals, 5222 of the Galla and 2162 of the SEA breeds, born between the years 1983 and 2022, were utilized. Individuals with known parents were defined as reference population. From the results, the maximum generation traced for Galla and SEA populations was 14.6 and 14.5, respectively. However, only 6 and 5 generations for Galla and SEA were complete. Pedigree completeness increased with the increasing number of generations. The average generation interval (GI) for Galla and SEA was 3.84 ± 0.04 and 4.4 ± 0.13 years. The average increase in the rate of inbreeding per generation for Galla and SEA was 0.04 and 0.05, with the effective ancestors and founders (fa/fe) ratio being same (1.00) for both breeds. Fifty percent (50%) of the genetic variability in the populations was contributed by 3 and 1 ancestor for Galla SEA, respectively. The effective population size (Ne) was 5.19 and 4.77 for Galla and SEA. Therefore, the current breeding programs should be changed to avoid future genetic bottlenecks in this population. These findings offer an opportunity to enhance the current genetic status and management of Kenyan native goats and other regions with similar production systems.

Long-term health outcomes from inbreeding in a historical Swedish population: longevity, fertility, and impairments.

BACKGROUND: Inbreeding, arising from consanguinity between related parents, has been observed to impact the health of individuals, typically attributed to biological factors. Nevertheless, these effects may be influenced by the social and environmental conditions. The prevalence of consanguineous marriages increased in certain parts of Sweden after it became legal in 1844, which offers a unique opportunity to study and understand the effects of inbreeding on health. AIM: The objective of this study is to explore the potential impact of inbreeding on the longevity, fertility, and impairments of individuals born in the Skellefteå region, Sweden, between 1890 and 1905, with a follow-up period extending until 1950. SUBJECTS AND METHODS: The level of inbreeding is calculated using micro-level parish register data and related to longevity, fertility, and impairments using regression analysis. RESULTS: Inbreeding is shown to be associated with longevity, fertility, and impairments. It seems to affect the risk of stillbirth and impairments and male longevity and fertility. CONCLUSION: Inbreeding seems to have had a detrimental effect on some health outcomes in this historical population under study.

Evaluating genomic inbreeding of two Chinese yak (Bos grunniens) populations.

BACKGROUND: Yaks are a vital livestock in the Qinghai-Tibetan Plateau area for providing food products, maintaining sustainable ecosystems, and promoting cultural heritage. Because of uncontrolled mating, it is impossible to estimate inbreeding level of yak populations using the pedigree-based approaches. With the aims to accurately evaluate inbreeding level of two Chinese yak populations (Maiwa and Jiulong), we obtained genome-wide single nucleotide polymorphisms (SNPs) by DNA sequencing and calculated five SNP-by-SNP estimators ([Formula: see text], [Formula: see text], [Formula: see text], [Formula: see text], and [Formula: see text]), as well as two segment-based estimators of runs of homozygosity (ROH, [Formula: see text]) and homozygous-by-descent (HBD, [Formula: see text]). Functional implications were analyzed for the positional candidate genes located within the related genomic regions. RESULTS: A total of 151,675 and 190,955 high-quality SNPs were obtained from 71 Maiwa and 30 Jiulong yaks, respectively. Jiulong had greater genetic diversity than Maiwa in terms of allele frequency and nucleotide diversity. The two populations could be genetically distinguished by principal component analysis, with the mean differentiation index (Fst) of 0.0054. The greater genomic inbreeding levels of Maiwa yaks were consistently supported by all five SNP-by-SNP estimators. Based on simple proportion of homozygous SNPs ([Formula: see text]), a lower inbreeding level was indicated by three successfully sequenced old leather samples that may represent historical Maiwa yaks about five generations ago. There were 3304 ROH detected among all samples, with mean and median length of 1.97 Mb and 1.0 Mb, respectively. A total of 94 HBD segments were found among all samples, whereas 92 of them belonged to the shortest class with the mean length of 10.9 Kb. Based on the estimates of [Formula: see text] and [Formula: see text], however, there was no difference in inbreeding level between Maiwa and Jiulong yaks. Within the genomic regions with the significant Fst or enriched by ROH, we found several candidate genes and pathways that have been reported to be related to diverse production traits in farm animals. CONCLUSIONS: We successfully evaluated the genomic inbreeding level of two Chinese yak populations. Although different estimators resulted in inconsistent conclusions on their genomic inbreeding levels, our results may be helpful to implement the genetic conservation and utilization programs for the two yak populations.

Using high-density SNP data to unravel the origin of the Franches-Montagnes horse breed.

BACKGROUND: The Franches-Montagnes (FM) is the last native horse breed of Switzerland, established at the end of the 19th century by cross-breeding local mares with Anglo-Norman stallions. We collected high-density SNP genotype data (Axiom™ 670 K Equine genotyping array) from 522 FM horses, including 44 old-type horses (OF), 514 European Warmblood horses (WB) from Sweden and Switzerland (including a stallion used for cross-breeding in 1990), 136 purebred Arabians (AR), 32 Shagya Arabians (SA), and 64 Thoroughbred (TB) horses, as introgressed WB stallions showed TB origin in their pedigrees. The aim of the study was to ascertain fine-scale population structures of the FM breed, including estimation of individual admixture levels and genomic inbreeding (FROH) by means of Runs of Homozygosity. RESULTS: To assess fine-scale population structures within the FM breed, we applied a three-step approach, which combined admixture, genetic contribution, and FROH of individuals into a high-resolution network visualization. Based on this approach, we were able to demonstrate that population substructures, as detected by model-based clustering, can be either associated with a different genetic origin or with the progeny of most influential sires. Within the FM breed, admixed horses explained most of the genetic variance of the current breeding population, while OF horses only accounted for a small proportion of the variance. Furthermore, we illustrated that FM horses showed high TB admixture levels and we identified inconsistencies in the origin of FM horses descending from the Arabian stallion Doktryner. With the exception of WB, FM horses were less inbred compared to the other breeds. However, the relatively few but long ROH segments suggested diversity loss in both FM subpopulations. Genes located in FM- and OF-specific ROH islands had known functions involved in conformation and behaviour, two traits that are highly valued by breeders. CONCLUSIONS: The FM remains the last native Swiss breed, clearly distinguishable from other historically introgressed breeds, but it suffered bottlenecks due to intensive selection of stallions, restrictive mating choices based on arbitrary definitions of pure breeding, and selection of rare coat colours. To preserve the genetic diversity of FM horses, future conservation managements strategies should involve a well-balanced selection of stallions (e.g., by integrating OF stallions in the FM breeding population) and avoid selection for rare coat colours.

Assessing different metrics of pedigree and genomic inbreeding and inbreeding effect on growth, fertility, and feed efficiency traits in a closed-herd Nellore cattle population.

BACKGROUND: The selection of individuals based on their predicted breeding values and mating of related individuals can increase the proportion of identical-by-descent alleles. In this context, the objectives of this study were to estimate inbreeding coefficients based on alternative metrics and data sources such as pedigree (FPED), hybrid genomic relationship matrix H (FH), and ROH of different length (FROH); and calculate Pearson correlations between the different metrics in a closed Nellore cattle population selected for body weight adjusted to 378 days of age (W378). In addition to total FROH (all classes) coefficients were also estimated based on the size class of the ROH segments: FROH1 (1-2 Mb), FROH2 (2-4 Mb), FROH3 (4-8 Mb), FROH4 (8-16 Mb), and FROH5 (> 16 Mb), and for each chromosome (FROH_CHR). Furthermore, we assessed the effect of each inbreeding metric on birth weight (BW), body weights adjusted to 210 (W210) and W378, scrotal circumference (SC), and residual feed intake (RFI). We also evaluated the chromosome-specific effects of inbreeding on growth traits. RESULTS: The correlation between FPED and FROH was 0.60 while between FH and FROH and FH and FPED were 0.69 and 0.61, respectively. The annual rate of inbreeding was 0.16% for FPED, 0.02% for FH, and 0.16% for FROH. A 1% increase in FROH5 resulted in a reduction of up to -1.327 ± 0.495 kg in W210 and W378. Four inbreeding coefficients (FPED, FH, FROH2, and FROH5) had a significant effect on W378, with reductions of up to -3.810 ± 1.753 kg per 1% increase in FROH2. There was an unfavorable effect of FPED on RFI (0.01 ± 0.0002 kg dry matter/day) and of FROH on SC (-0.056 ± 0.022 cm). The FROH_CHR coefficients calculated for BTA3, BTA5, and BTA8 significantly affected the growth traits. CONCLUSIONS: Inbreeding depression was observed for all traits evaluated. However, these effects were greater for the criterion used for selection of the animals (i.e., W378). The increase in the genomic inbreeding was associated with a higher inbreeding depression on the traits evaluated when compared to pedigree-based inbreeding. Genomic information should be used as a tool during mating to optimize control of inbreeding and, consequently, minimize inbreeding depression in Nellore cattle.

Genetic diversity of United States Rambouillet, Katahdin and Dorper sheep.

BACKGROUND: Managing genetic diversity is critically important for maintaining species fitness. Excessive homozygosity caused by the loss of genetic diversity can have detrimental effects on the reproduction and production performance of a breed. Analysis of genetic diversity can facilitate the identification of signatures of selection which may contribute to the specific characteristics regarding the health, production and physical appearance of a breed or population. In this study, breeds with well-characterized traits such as fine wool production (Rambouillet, N = 745), parasite resistance (Katahdin, N = 581) and environmental hardiness (Dorper, N = 265) were evaluated for inbreeding, effective population size (Ne), runs of homozygosity (ROH) and Wright's fixation index (FST) outlier approach to identify differential signatures of selection at 36,113 autosomal single nucleotide polymorphisms (SNPs). RESULTS: Katahdin sheep had the largest current Ne at the most recent generation estimated with both the GONe and NeEstimator software. The most highly conserved ROH Island was identified in Rambouillet with a signature of selection on chromosome 6 containing 202 SNPs called in an ROH in 50 to 94% of the individuals. This region contained the DCAF16, LCORL and NCAPG genes that have been previously reported to be under selection and have biological roles related to milk production and growth traits. The outlier regions identified through the FST comparisons of Katahdin with Rambouillet and Dorper contained genes with known roles in milk production and mastitis resistance or susceptibility, and the FST comparisons of Rambouillet with Katahdin and Dorper identified genes related to wool growth, suggesting these traits have been under natural or artificial selection pressure in these populations. Genes involved in the cytokine-cytokine receptor interaction pathways were identified in all FST breed comparisons, which indicates the presence of allelic diversity between these breeds in genomic regions controlling cytokine signaling mechanisms. CONCLUSIONS: In this paper, we describe signatures of selection within diverse and economically important U.S. sheep breeds. The genes contained within these signatures are proposed for further study to understand their relevance to biological traits and improve understanding of breed diversity.

Multi-generational benefits of genetic rescue.

Genetic rescue-an increase in population fitness following the introduction of new alleles-has been proven to ameliorate inbreeding depression in small, isolated populations, yet is rarely applied as a conservation tool. A lingering question regarding genetic rescue in wildlife conservation is how long beneficial effects persist in admixed populations. Using data collected over 40 years from 1192 endangered Florida panthers (Puma concolor coryi) across nine generations, we show that the experimental genetic rescue implemented in 1995-via the release of eight female pumas from Texas-alleviated morphological, genetic, and demographic correlates of inbreeding depression, subsequently preventing extirpation of the population. We present unequivocal evidence, for the first time in any terrestrial vertebrate, that genetic and phenotypic benefits of genetic rescue remain in this population after five generations of admixture, which helped increase panther abundance (> fivefold) and genetic effective population size (> 20-fold). Additionally, even with extensive admixture, microsatellite allele frequencies in the population continue to support the distinctness of Florida panthers from other North American puma populations, including Texas. Although threats including habitat loss, human-wildlife conflict, and infectious diseases are challenges to many imperiled populations, our results suggest genetic rescue can serve as an effective, multi-generational tool for conservation of small, isolated populations facing extinction from inbreeding.

Genomic inbreeding analysis reveals resilience and genetic diversity in Indian yak populations.

The yak (Bos grunniens), renowned for its adaptability to extreme cold and hypoxic conditions, stands as a remarkable domestic animal crucial for sustaining livelihoods in harsh climates. We conducted a comprehensive analysis of the whole genome sequence data from three distinct Indian yak populations: Arunachali yak (n = 10), Himachali yak (n = 10), and Ladakhi yak (n = 10). The genomic data for Indian yaks were meticulously generated by our laboratory and compared with their Chinese counterpart, the Jinchuan yak (n = 8), for a more nuanced understanding. Our investigation revealed a total of 37,437 runs of homozygosity (ROH) segments in 34 animals representing four distinct yak populations. The Jinchuan yak population exhibited the highest proportion, constituting 80.8 % of total ROHs, predominantly as small segments (<0.1 Mb), accounting for 63 % of the overall ROHs. Further analysis uncovered a significantly higher degree of inbreeding in Chinese yaks compared to their Indian counterparts. The Indian yak populations, in contrast, demonstrated relatively lower and consistent levels of inbreeding. Moreover, we identified ROH hotspots that covered at least 60 % of individuals in our study, indicating their pivotal role in environmental adaptation. A total of five hotspot regions were detected, housing genes such as ENSBGRG00000015023 (WNT2), YIPF4, SPAST, TLN2, and DSG4. These genes are associated with traits including hair follicle initiation, nutrient stress response, microtubule assembly, development of cardiac muscle, hair follicle, and coat color. This observation strongly suggests that there is substantial selection acting on these genes, emphasizing their important role in environmental adaptation among yak populations.

Runs of homozygosity reveal contrasting histories of inbreeding across global lineages of the edible porcini mushroom, Boletus edulis.

Inbreeding, the mating of individuals that are related through common ancestry, is of central importance in evolutionary and conservation biology due to its impacts on individual fitness and population dynamics. However, while advanced genomic approaches have revolutionised the study of inbreeding in animals, genomic studies of inbreeding are rare in plants and lacking in fungi. We investigated global patterns of inbreeding in the prized edible porcini mushroom Boletus edulis using 225 whole genomes from seven lineages distributed across the northern hemisphere. Genomic inbreeding was quantified using runs of homozygosity (ROHs). We found appreciable variation both among and within lineages, with some individuals having over 20% of their genomes in ROHs. Much of this variation could be explained by a combination of elevation and latitude, and to a lesser extent by predicted habitat suitability during the last glacial maximum. In line with this, the majority of ROHs were short, reflecting ancient common ancestry dating back approximately 200-1700 generations ago, while longer ROHs indicative of recent common ancestry (less than approximately 50 generations ago) were infrequent. Our study reveals the inbreeding legacy of major climatic events in a widely distributed forest mutualist, aligning with prevailing theories and empirical studies of the impacts of historical glaciation events on the dominant forest tree species of the northern hemisphere.

Selfing Promotes Spread and Introgression of Segregation Distorters in Hermaphroditic Plants.

Segregation distorters (SDs) are genetic elements that distort the Mendelian segregation ratio to favor their own transmission and are able to spread even when they incur fitness costs on organisms carrying them. Depending on the biology of the host organisms and the genetic architecture of the SDs, the population dynamics of SDs can be highly variable. Inbreeding is considered an effective mechanism for inhibiting the spread of SDs in populations, and can evolve as a defense mechanism against SDs in some systems. However, we show that inbreeding in the form of selfing in fact promotes the spread of SDs acting as pollen killers in a toxin-antidote system in hermaphroditic plants by two mechanisms: (i) By reducing the effective recombination rate between killer and antidote loci in the two-locus system and (ii) by increasing the proportion of SD alleles in individual flowers, rather than in the general gene-pool. We also show that in rice (Oryza sativa L.), a typical hermaphroditic plant, all molecularly characterized SDs associated with pollen killing were involved in population hybridization and have introgressed across different species. Paradoxically, these loci, which are associated with hybrid incompatibility and can be thought of as Bateson-Dobzhansky-Muller incompatibility loci are expected to reduce gene-flow between species, in fact cross species boundaries more frequently than random loci, and may act as important drivers of introgression.

Genomic Consequences of Isolation and Inbreeding in an Island Dingo Population.

Dingoes come from an ancient canid lineage that originated in East Asia around 8,000 to 11,000 years BP. As Australia's largest terrestrial predator, dingoes play an important ecological role. A small, protected population exists on a world heritage listed offshore island, K'gari (formerly Fraser Island). Concern regarding the persistence of dingoes on K'gari has risen due to their low genetic diversity and elevated inbreeding levels. However, whole-genome sequence data is lacking from this population. Here, we include five new whole-genome sequences of K'gari dingoes. We analyze a total of 18 whole-genome sequences of dingoes sampled from mainland Australia and K'gari to assess the genomic consequences of their demographic histories. Long (>1 Mb) runs of homozygosity (ROHs)-indicators of inbreeding-are elevated in all sampled dingoes. However, K'gari dingoes showed significantly higher levels of very long ROH (>5 Mb), providing genomic evidence for small population size, isolation, inbreeding, and a strong founder effect. Our results suggest that, despite current levels of inbreeding, the K'gari population is purging strongly deleterious mutations, which, in the absence of further reductions in population size, may facilitate the persistence of small populations despite low genetic diversity and isolation. However, there may be little to no purging of mildly deleterious alleles, which may have important long-term consequences, and should be considered by conservation and management programs.

Using Runs of Homozygosity and Machine Learning to Disentangle Sources of Inbreeding and Infer Self-Fertilization Rates.

Runs of homozygosity (ROHs) are indicative of elevated homozygosity and inbreeding due to mating of closely related individuals. Self-fertilization can be a major source of inbreeding which elevates genome-wide homozygosity and thus should also create long ROHs. While ROHs are frequently used to understand inbreeding in the context of conservation and selective breeding, as well as for consanguinity of populations and their demographic history, it remains unclear how ROH characteristics are altered by selfing and if this confounds expected signatures of inbreeding due to demographic change. Using simulations, we study the impact of the mode of reproduction and demographic history on ROHs. We apply random forests to identify unique characteristics of ROHs, indicative of different sources of inbreeding. We pinpoint distinct features of ROHs that can be used to better characterize the type of inbreeding the population was subjected to and to predict outcrossing rates and complex demographic histories. Using additional simulations and four empirical datasets, two from highly selfing species and two from mixed-maters, we predict the selfing rate and validate our estimations. We find that self-fertilization rates are successfully identified even with complex demography. Population genetic summary statistics improve algorithm accuracy particularly in the presence of additional inbreeding, e.g. from population bottlenecks. Our findings highlight the importance of ROHs in disentangling confounding factors related to various sources of inbreeding and demonstrate situations where such sources cannot be differentiated. Additionally, our random forest models provide a novel tool to the community for inferring selfing rates using genomic data.

Measuring the Efficiency of Purging by non-random Mating in Human Populations.

Human populations harbor a high concentration of deleterious genetic variants. Here, we tested the hypothesis that non-random mating practices affect the distribution of these variants, through exposure in the homozygous state, leading to their purging from the population gene pool. To do so, we produced whole-genome sequencing data for two pairs of Asian populations exhibiting different alliance rules and rates of inbreeding, but with similar effective population sizes. The results show that populations with higher rates of inbred matings do not purge deleterious variants more efficiently. Purging therefore has a low efficiency in human populations, and different mating practices lead to a similar mutational load.