Genome-wide data from medieval German Jews show that the Ashkenazi founder event pre-dated the 14th century.

We report genome-wide data from 33 Ashkenazi Jews (AJ), dated to the 14th century, obtained following a salvage excavation at the medieval Jewish cemetery of Erfurt, Germany. The Erfurt individuals are genetically similar to modern AJ, but they show more variability in Eastern European-related ancestry than modern AJ. A third of the Erfurt individuals carried a mitochondrial lineage common in modern AJ and eight carried pathogenic variants known to affect AJ today. These observations, together with high levels of runs of homozygosity, suggest that the Erfurt community had already experienced the major reduction in size that affected modern AJ. The Erfurt bottleneck was more severe, implying substructure in medieval AJ. Overall, our results suggest that the AJ founder event and the acquisition of the main sources of ancestry pre-dated the 14th century and highlight late medieval genetic heterogeneity no longer present in modern AJ.

Declining autozygosity over time: An exploration in over 1 million individuals from three diverse cohorts.

Previous studies have hypothesized that autozygosity is decreasing over generational time. However, these studies were limited to relatively small samples (n < 11,000) lacking in diversity, which may limit the generalizability of their findings. We present data that partially support this hypothesis from three large cohorts of diverse ancestries, two from the US (All of Us, n = 82,474; the Million Veteran Program, n = 622,497) and one from the UK (UK Biobank, n = 380,899). Our results from a mixed-effect meta-analysis demonstrate an overall trend of decreasing autozygosity over generational time (meta-analyzed slope = -0.029, SE = 0.009, p = 6.03e-4). On the basis of our estimates, we would predict FROH to decline 0.29% for every 20-year increase in birth year. We determined that a model including an ancestry-by-country interaction term fit the data best, indicating that ancestry differences in this trend differ by country. We found further evidence to suggest a difference between the US and UK cohorts by meta-analyzing within country, observing a significant negative estimate in the US cohorts (meta-analyzed slope = -0.058, SE = 0.015, p = 1.50e-4) but a non-significant estimate in the UK (meta-analyzed slope = -0.001, SE = 0.008, p = 0.945). The association between autozygosity and birth year was substantially attenuated when accounting for educational attainment and income (meta-analyzed slope = -0.011, SE = 0.008, p = 0.167), suggesting they may partially account for decreasing autozygosity over time. Overall, we demonstrate decreasing autozygosity over time in a large, modern sample and speculate that this trend can be attributed to increases in urbanization and panmixia and differences in sociodemographic processes lead to country-specific differences in the rate of decline.

Increased homozygosity due to endogamy results in fitness consequences in a human population.

Recessive alleles have been shown to directly affect both human Mendelian disease phenotypes and complex traits. Pedigree studies also suggest that consanguinity results in increased childhood mortality and adverse health phenotypes, presumably through penetrance of recessive mutations. Here, we test whether the accumulation of homozygous, recessive alleles decreases reproductive success in a human population. We address this question among the Namibian Himba, an endogamous agro-pastoralist population, who until very recently practiced natural fertility. Using a sample of 681 individuals, we show that Himba exhibit elevated levels of "inbreeding," calculated as the fraction of the genome in runs of homozygosity (FROH). Many individuals contain multiple long segments of ROH in their genomes, indicating that their parents had high kinship coefficients. However, we do not find evidence that this is explained by first-cousin consanguinity, despite a reported social preference for cross-cousin marriages. Rather, we show that elevated haplotype sharing in the Himba is due to a bottleneck, likely in the past 60 generations. We test whether increased recessive mutation load results in observed fitness consequences by assessing the effect of FROH on completed fertility in a cohort of postreproductive women (n = 69). We find that higher FROH is significantly associated with lower fertility. Our data suggest a multilocus genetic effect on fitness driven by the expression of deleterious recessive alleles, especially those in long ROH. However, these effects are not the result of consanguinity but rather elevated background identity by descent.

Conservation management strategy impacts inbreeding and mutation load in scimitar-horned oryx.

In an age of habitat loss and overexploitation, small populations, both captive and wild, are increasingly facing the effects of isolation and inbreeding. Genetic management has therefore become a vital tool for ensuring population viability. However, little is known about how the type and intensity of intervention shape the genomic landscape of inbreeding and mutation load. We address this using whole-genome sequence data of the scimitar-horned oryx (Oryx dammah), an iconic antelope that has been subject to contrasting management strategies since it was declared extinct in the wild. We show that unmanaged populations are enriched for long runs of homozygosity (ROH) and have significantly higher inbreeding coefficients than managed populations. Additionally, despite the total number of deleterious alleles being similar across management strategies, the burden of homozygous deleterious genotypes was consistently higher in unmanaged groups. These findings emphasize the risks associated with deleterious mutations through multiple generations of inbreeding. As wildlife management strategies continue to diversify, our study reinforces the importance of maintaining genome-wide variation in vulnerable populations and has direct implications for one of the largest reintroduction attempts in the world.

Runs of homozygosity and selection signature analyses reveal putative genomic regions for artificial selection in layer breeding.

BACKGROUND: The breeding of layers emphasizes the continual selection of egg-related traits, such as egg production, egg quality and eggshell, which enhance their productivity and meet the demand of market. As the breeding process continued, the genomic homozygosity of layers gradually increased, resulting in the emergence of runs of homozygosity (ROH). Therefore, ROH analysis can be used in conjunction with other methods to detect selection signatures and identify candidate genes associated with various important traits in layer breeding. RESULTS: In this study, we generated whole-genome sequencing data from 686 hens in a Rhode Island Red population that had undergone fifteen consecutive generations of intensive artificial selection. We performed a genome-wide ROH analysis and utilized multiple methods to detect signatures of selection. A total of 141,720 ROH segments were discovered in whole population, and most of them (97.35%) were less than 3 Mb in length. Twenty-three ROH islands were identified, and they overlapped with some regions bearing selection signatures, which were detected by the De-correlated composite of multiple signals methods (DCMS). Sixty genes were discovered and functional annotation analysis revealed the possible roles of them in growth, development, immunity and signaling in layers. Additionally, two-tailed analyses including DCMS and ROH for 44 phenotypes of layers were conducted to find out the genomic differences between subgroups of top and bottom 10% phenotype of individuals. Combining the results of GWAS, we observed that regions significantly associated with traits also exhibited selection signatures between the high and low subgroups. We identified a region significantly associated with egg weight near the 25 Mb region of GGA 1, which exhibited selection signatures and has higher genomic homozygosity in the low egg weight subpopulation. This suggests that the region may be play a role in the decline in egg weight. CONCLUSIONS: In summary, through the combined analysis of ROH, selection signatures, and GWAS, we identified several genomic regions that associated with the production traits of layers, providing reference for the study of layer genome.

Genetic assessment and candidate genes identification for breed-specific characteristics of Qingyuan partridge chicken based on runs of homozygosity.

BACKGROUND: Several core breeding and supporting lines of the Qingyuan partridge chicken, a representative local chicken breed in China, have been developed over 20 years. Consequently, its economic traits related to growth and reproduction have been significantly improved by breeding selection and commercial utilization, but some characteristic traits, such as partridge feathers, high meat quality and sufficient flavor, have always been retained. However, effective methods for genetic assessment and functional gene exploration of similar trait groups are lacking. The presence of identical haplotype fragments transmitted from parent to offspring results in runs of homozygosity (ROH), which offer an efficient solution. In this study, genomes of 134 Qingyuan partridge chickens representing two breeding populations and one preserved population were re-sequenced to evaluate the genetic diversity and explore functional genes by analyzing the diversity, distribution, and frequency of ROH. RESULTS: The results showed a low level of genomic linkage and degree of inbreeding within both the bred and preserved populations, suggesting abundant genetic diversity and an adequate genetic potential of the Qingyuan partridge chicken. Throughout the long-term selection process, 21 genes, including GLI3, ANO5, BLVRA, EFNB2, SLC5A12, and SVIP, associated with breed-specific characteristics were accumulated within three ROH islands, whereas another 21 genes associated with growth traits including IRX1, IRX2, EGFR, TPK1, NOVA1, BDNF and so on were accumulated within five ROH islands. CONCLUSIONS: These findings provide new insights into the genetic assessment and identification of genes with breed-specific and selective characteristics, offering a solid genetic basis for breeding and protection of Qingyuan partridge chickens.

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.

Genome-wide detections for runs of homozygosity and selective signatures reveal novel candidate genes under domestication in chickens.

BACKGROUND: Indigenous chickens were developed through a combination of natural and artificial selection; essentially, changes in genomes led to the formation of these modern breeds via admixture events. However, their confusing genetic backgrounds include a genomic footprint regulating complex traits, which is not conducive to modern animal breeding. RESULTS: To better evaluate the candidate regions under domestication in indigenous chickens, we considered both runs of homozygosity (ROHs) and selective signatures in 13 indigenous chickens. The genomes of Silkie feather chickens presented the highest heterozygosity, whereas the highest inbreeding status and ROH number were found in Luhua chickens. Short ROH (< 1 Mb), were the principal type in all chickens. A total of 291 ROH islands were detected, and QTLdb mapping results indicated that body weight and carcass traits were the most important traits. An ROH on chromosome 2 covering VSTM2A gene was detected in 12 populations. Combined analysis with the Tajima's D index revealed that 18 genes (e.g., VSTM2A, BBOX1, and RYR2) were under selection and covered by ROH islands. Transcriptional analysis results showed that RYR2 and BBOX1 were specifically expressed in the heart and muscle tissue, respectively. CONCLUSION: Based on genome-wide scanning for ROH and selective signatures, we evaluated the genomic characteristics and detected significant candidate genes covered by ROH islands and selective signatures. The findings in this study facilitated the understanding of genetic diversity and provided valuable insights for chicken breeding and conservation strategies.

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.

Extensive Phylogenomic Discordance and the Complex Evolutionary History of the Neotropical Cat Genus Leopardus.

Even in the genomics era, the phylogeny of Neotropical small felids comprised in the genus Leopardus remains contentious. We used whole-genome resequencing data to construct a time-calibrated consensus phylogeny of this group, quantify phylogenomic discordance, test for interspecies introgression, and assess patterns of genetic diversity and demographic history. We infer that the Leopardus radiation started in the Early Pliocene as an initial speciation burst, followed by another in its subgenus Oncifelis during the Early Pleistocene. Our findings challenge the long-held notion that ocelot (Leopardus pardalis) and margay (L. wiedii) are sister species and instead indicate that margay is most closely related to the enigmatic Andean cat (L. jacobita), whose whole-genome data are reported here for the first time. In addition, we found that the newly sampled Andean tiger cat (L. tigrinus pardinoides) population from Colombia associates closely with Central American tiger cats (L. tigrinus oncilla). Genealogical discordance was largely attributable to incomplete lineage sorting, yet was augmented by strong gene flow between ocelot and the ancestral branch of Oncifelis, as well as between Geoffroy's cat (L. geoffroyi) and southern tiger cat (L. guttulus). Contrasting demographic trajectories have led to disparate levels of current genomic diversity, with a nearly tenfold difference in heterozygosity between Andean cat and ocelot, spanning the entire range of variability found in extant felids. Our analyses improved our understanding of the speciation history and diversity patterns in this felid radiation, and highlight the benefits to phylogenomic inference of embracing the many heterogeneous signals scattered across the genome.

High prevalence of multilocus pathogenic variation in neurodevelopmental disorders in the Turkish population.

Neurodevelopmental disorders (NDDs) are clinically and genetically heterogenous; many such disorders are secondary to perturbation in brain development and/or function. The prevalence of NDDs is > 3%, resulting in significant sociocultural and economic challenges to society. With recent advances in family-based genomics, rare-variant analyses, and further exploration of the Clan Genomics hypothesis, there has been a logarithmic explosion in neurogenetic "disease-associated genes" molecular etiology and biology of NDDs; however, the majority of NDDs remain molecularly undiagnosed. We applied genome-wide screening technologies, including exome sequencing (ES) and whole-genome sequencing (WGS), to identify the molecular etiology of 234 newly enrolled subjects and 20 previously unsolved Turkish NDD families. In 176 of the 234 studied families (75.2%), a plausible and genetically parsimonious molecular etiology was identified. Out of 176 solved families, deleterious variants were identified in 218 distinct genes, further documenting the enormous genetic heterogeneity and diverse perturbations in human biology underlying NDDs. We propose 86 candidate disease-trait-associated genes for an NDD phenotype. Importantly, on the basis of objective and internally established variant prioritization criteria, we identified 51 families (51/176 = 28.9%) with multilocus pathogenic variation (MPV), mostly driven by runs of homozygosity (ROHs) - reflecting genomic segments/haplotypes that are identical-by-descent. Furthermore, with the use of additional bioinformatic tools and expansion of ES to additional family members, we established a molecular diagnosis in 5 out of 20 families (25%) who remained undiagnosed in our previously studied NDD cohort emanating from Turkey.

Genetic architecture of inbreeding depression may explain its persistence in a population of wild red deer.

Inbreeding depression is of major concern in declining populations, but relatively little is known about its genetic architecture in wild populations, such as the degree to which it is composed of large or small effect loci and their distribution throughout the genome. Here, we combine fitness and genomic data from a wild population of red deer to investigate the genomic distribution of inbreeding effects. Based on the runs of homozygosity (ROH)-based inbreeding coefficient, FROH, we use chromosome-specific inbreeding coefficients (FROHChr) to explore whether the effect of inbreeding varies between chromosomes. Under the assumption that within an individual the probability of being identical-by-descent is equal across all chromosomes, we used a multi-membership model to estimate the deviation of FROHChr from the average inbreeding effect. This novel approach ensures effect sizes are not overestimated whilst maximising the power of our available dataset of >3000 individuals genotyped on >35,000 autosomal SNPs. We find that most chromosomes confer a minor reduction in fitness-related traits, which when these effects are summed, results in the observed inbreeding depression in birth weight, survival and lifetime breeding success. However, no chromosomes had a significant detrimental effect compared to the overall effect of inbreeding, indicating no major effect loci. We conclude that in this population, inbreeding depression is likely the result of multiple mildly or moderately deleterious mutations spread across all chromosomes, which are difficult to detect with statistical confidence. Such mutations will be inefficiently purged, which may explain the persistence of inbreeding depression in this population.

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.

Genomic implications of the repeated shift to self-fertilization across a species' geographic distribution.

The ability to self-fertilize often varies among closely related hermaphroditic plant species, though, variation can also exist within species. In the North American Arabidopsis lyrata, the shift from self-incompatibility (SI) to selfing established in multiple regions independently, mostly since recent postglacial range expansion. This has made the species an ideal model for the investigation of the genomic underpinnings of the breakdown of SI and its population genetic consequences. By comparing nearby selfing and outcrossing populations across the entire species' geographic distribution, we investigated variation at the self-incompatibility (S-)locus and across the genome. Furthermore, a diallel crossing experiment on one mixed-mating population was performed to gain insight in the genetics of mating system variation. We confirmed that the breakdown of SI had evolved in several S-locus backgrounds. The diallel suggested the involvement of binuclearly expressed parental genes with dominance relations. Though, the population-level genome-wide association study did not single out clear-cut candidate genes but several regions with one near the S-locus. On the implication side, selfing as compared to outcrossing populations had less than half of the genomic diversity, while the number of runs of homozygosity and their length scaled with the degree of inbreeding. The results highlight that mating system shifts to selfing, its genetic underpinning and the likely negative genomic consequences for evolutionary potential can be strongly interlinked with past range dynamics.

Population structure and genomic diversity of the Einsiedler horse.

The breeding history of the Einsiedler horse is closely connected with the Benedictine cloister Einsiedeln. In the mid-nineteenth century, it was decided to use European Warmblood stallions for cross-breeding and to abandon the selection of stallions. Since that time, it has only been possible to trace back the origin of Einsiedler horses using maternal ancestry information. Here, we collected high-density genotype data for European Warmblood horses (Selle Français, Swiss Warmblood and Einsiedler) and Franches-Montagnes horses, the last native Swiss horse breed, to unravel the current population structure of the Einsiedler horse. Using commonly applied methods to ascertain fine-scale population structures, it was not possible to clearly differentiate the Einsiedler from other European Warmblood horses. However, by means of runs of homozygosity (ROH) we were able to detect breed-specific ROH islands for the Einsiedler horse, including genes involved in domestication and adaptation to high altitude. Therefore, future breeding activities should involve the screening of these breed-specific ROH segments, the revival of cryopreserved sperm and the selection of Einsiedler stallions.

Characterizing Homozygous Variants in Bardet-Biedl Syndrome-Associated Genes Within Iranian Families: Unveiling a Founder Variant in BBS2, c.471G>A.

Bardet-Biedl syndrome (BBS) is a rare inherited ciliopathy disorder characterized by a broad spectrum of clinical symptoms such as retinal dystrophy, obesity, polydactyly, genitourinary and kidney anomalies, learning disability, and hypogonadism. The understanding of the variants involved in BBS-causing genes remains incomplete, highlighting the need for further research to develop a molecular diagnostic strategy for this syndrome. Singleton whole-exome sequencing (WES) was performed on sixteen patients. Our study revealed (1) nine patients carried eight homozygous pathogenic variants with four of them being novel (2) Specifically, a synonymous splicing variant (c.471G > A) in BBS2 gene in six patients with Baloch ethnicity. The identification of runs of homozygosity (ROH) calling was performed using the BCFtools/RoH software on WES data of patients harboring c.471G > A variant. The presence of shared homozygous regions containing the identified variant was confirmed in these patients. In-silico analysis predicted the effect of the c.471G > A variants on BBS2 mRNA splicing. This variant results in disrupted wild-type donor site and intron retention in the mature mRNA. (3) And a deletion of exons 14 to 17 in the BBS1 gene was identified in one patient by Copy-Number Variation (CNV) analysis using the ExomeDepth pipeline. Our results identified the founder variant c.471G > A in the BBS2 gene in the Baloch ethnicity of the Iranian population. This finding can guide the diagnostic approach of this syndrome in future studies.

Reproductive traits in Pura Raza Española mares manifest inbreeding depression from low levels of homozygosity.

Inbreeding depression is a genetic phenomenon associated with the loss of fitness and mean phenotypic performance due to mating between relatives. Historically, inbreeding coefficients have been estimated from pedigree information. However, the onset of genomic selection programs provides large datasets of individuals genotyped using SNP arrays, enabling more precise assessment of an individual's genomic-level inbreeding using genomic data. One of the traits most sensitive to issues stemming from increased inbreeding is reproduction. This is particularly important in equine, in which fertility is only moderate compared to other livestock species. To explore this further, we evaluated the effect of inbreeding on five reproductive traits (age at first foaling (AFF), average interval between foalings (AIF), total number of foalings (NF), productive life (PL) and reproductive efficiency (RE)) in Pura Raza Español mares using genomic data. Residual predicted phenotypes were obtained by purging these traits through the REML (wgResidual) and ssGREML (gResidual) approaches in reproductive data of 29,847 PRE mares using the BLUPF90+ program. Next, we used pedigree-based (Fped) and ROH-based genomic (FROH) inbreeding coefficients derived from 1018 animals genotyped with 61,271 SNPs to estimate the inbreeding depression (linear regression). Our results indicated significant levels of inbreeding depression for all reproductive traits, with the exception of the AIF trait when Fped was used. However, all traits were negatively affected by the increase in genomic inbreeding, and FROH was found to capture more inbreeding depression than Fped. Likewise, REML models (ssGREML) using genomic data for estimated predicted residual phenotypes resulted in higher variance explained by the model compared with the models not using genomics (REML). Finally, a segmented regression analysis was conducted to evaluate the effect of inbreeding depression, revealing that the levels of genealogical and genomic homozygosity do not manifest uniformly in reproductive traits. In contrast, the levels of inbreeding depression ranged from low to high as homozygosity increased. This analysis also showed that reproductive traits are very sensitive to inbreeding depression, even with relatively low levels of homozygosity.

Differential patterns in runs of homozygosity in two mice lines under divergent selection for environmental variability for birth weight.

Runs of homozygosity (ROH) are defined as long continuous homozygous stretches in the genome which are assumed to originate from a common ancestor. It has been demonstrated that divergent selection for variability in mice is possible and that low variability in birth weight is associated with robustness. To analyse ROH patterns and ROH-based genomic inbreeding, two mouse lines that were divergently selected for birth weight variability for 26 generations were used, with: 752 individuals for the high variability line (H-Line), 766 individuals for the low variability line (L-Line) and 74 individuals as a reference population. Individuals were genotyped using the high density Affymetrix Mouse Diversity Genotyping Array. ROH were identified using both the sliding windows (SW) and the consecutive runs (CR) methods. Inbreeding coefficients were calculated based on pedigree (FPED ) information, on ROH identified using the SW method (FROHSW ) and on ROH identified using the CR method (FROHCR ). Differences in genomic inbreeding were not consistent across generations and these parameters did not show clear differences between lines. Correlations between FPED and FROH were high, particularly for FROHSW . Moreover, correlations between FROHSW and FPED were even higher when ROH were identified with no restrictions in the number of heterozygotes per ROH. The comparison of FROH estimates between either of the selected lines were based on significant differences at the chromosome level, mainly in chromosomes 3, 4, 6, 8, 11, 15 and 19. ROH-based inbreeding estimates that were computed using longer homozygous segments had a higher relationship with FPED . Differences in robustness between lines were not attributable to a higher homozygosis in the L-Line, but maybe to the different distribution of ROH at the chromosome level between lines. The analysis identified a set of genomic regions for future research to establish the genomic basis of robustness.

Characterization of runs of homozygosity islands in American mink using whole-genome sequencing data.

The genome-wide analysis of runs of homozygosity (ROH) islands can be an effective strategy for identifying shared variants within a population and uncovering important genomic regions related to complex traits. The current study performed ROH analysis to characterize the genome-wide patterns of homozygosity, identify ROH islands and annotated genes within these candidate regions using whole-genome sequencing data from 100 American mink (Neogale vison). After sequence processing, variants were called using GATK and Samtools pipelines. Subsequent to quality control, 8,373,854 bi-allelic variants identified by both pipelines remained for further analysis. A total of 34,652 ROH segments were identified in all individuals, among which shorter segments (0.3-1 Mb) were abundant throughout the genome, approximately accounting for 84.39% of all ROH. Within these segments, we identified 63 ROH islands housing 156 annotated genes. The genes located in ROH islands were associated with fur quality (EDNRA, FGF2, FOXA2 and SLC24A4), body size/weight (MYLK4, PRIM2, FABP2, EYS and PHF3), immune capacity (IL2, IL21, PTP4A1, SEMA4C, JAK2, CCNA2 and TNIP3) and reproduction (ADAD1, KHDRBS2, INSL6, PGRMC2 and HSPA4L). Furthermore, Gene Ontology and KEGG pathway enrichment analyses revealed 56 and 9 significant terms (FDR-corrected p-value < 0.05), respectively, among which cGMP-PKG signalling pathway, regulation of actin cytoskeleton, and calcium signalling pathway were highlighted due to their functional roles in growth and fur characteristics. This is the first study to present ROH islands in American mink. The candidate genes from ROH islands and functional enrichment analysis suggest possible signatures of selection in response to the mink breeding targets, such as increased body length, reproductive performance and fur quality. These findings contribute to our understanding of genetic characteristics, and provide complementary information to assist with implementation of breeding strategies for genetic improvement in American mink.

Genetic diversity, population structure and origin of the native goats in Central Laos.

Maintaining genetic diversity and variation in livestock populations is critical for natural and artificial selection promoting genetic improvement while avoiding problems due to inbreeding. In Laos, there are concerns that there has been a decline in genetic diversity and a rise in inbreeding among native goats in their village-based smallholder system. In this study, we investigated the genetic diversity of Lao native goats in Phin, Songkhone and Sepon districts in Central Laos for the first time using Illumina's Goat SNP50 BeadChip. We also explored the genetic relationships between Lao goats with 163 global goat populations from 36 countries. Our results revealled a close genetic relationship between Lao native goats and Chinese, Mongolian and Pakistani goats, sharing ancestries with Guangfen, Jining Grey and Luoping Yellow breeds (China) and Teddi goats (Pakistan). The observed (Ho) and expected (He) heterozygosity were 0.292 and 0.303 (Laos), 0.288 and 0.288 (Sepon), 0.299 and 0.308 (Phin) and 0.289 and 0.305 (Songkhone), respectively. There was low to moderate genetic differentiation (FST: 0.011-0.043) and negligible inbreeding coefficients (FIS: -0.001 to 0.052) between goat districts. The runs of homozygosity (ROH) had an average length of 5.92-6.85 Mb, with short ROH segments (1-5 Mb length) being the most prevalent (66.34%). Longer ROH segments (20-40 and >40 Mb length categories) were less common, comprising only 4.81% and 1.01%, respectively. Lao goats exhibit moderate genetic diversity, low-inbreeding levels and adequate effective population size. Some genetic distinctions between Lao goats may be explained by geographic and cultural features.