Inbreeding Coefficient Research Articles

Contrasting population genetic structure of three semi‐terrestrial brachyuran crabs on the coast of the Japanese archipelago

AbstractAnthropogenic activities have reduced ecotones between the ocean and land, which is likely to threaten the population of brackish‐water brachyuran crabs. To assess the current status of these crabs, we examine the population genetic structures of three semi‐terrestrial brachyuran crabs widely distributed along the coast of the Japan and to clarify factors determining their genetic structures. We collected 184 Orisarma dehaani, 252 Chiromantes haematocheir, and 151 Helice tridens crabs from 36 localities of the Japanese archipelago. Genome‐wide SNP data from these crabs were analyzed using MIG‐seq. Bayesian clustering of STRUCTURE and DAPC analysis were used to identify genetically disturbed populations and to visualize genetic differentiation between local populations. Genetic population structure showed clear differentiation between populations on the Pacific coast of the Tohoku region and on other Japanese coasts in O. dehaani, but not in C. haematocheir or H. tridens. The inbreeding coefficient of O. dehaani was significantly higher on the Pacific coast of the Tohoku region compared to other Japanese coasts. C. haematocheir and H. tridens had homogeneous genetic structures along the Japanese coast, but showed genetic differentiation of a local population at their range limits. Thus, O. dehaani showed little gene flow and clear genetic differentiation between populations in the Tohoku Pacific region and those on other Japanese coasts due to ocean currents. Although such a regional differentiation was not found in C. haematocheir and H. tridens, one population of C. haematocheir was genetically isolated at the edge of its distribution range and likely vulnerable to environmental changes.

한우의 계획교배를 위한 근교계수가 도체형질에 미치는 영향 연구

한우의 계획교배를 위한 근교계수가 도체형질에 미치는 영향 연구

Genomic insights into the critically endangered King Island scrubtit.

Small, fragmented or isolated populations are at risk of population decline due to fitness costs associated with inbreeding and genetic drift. The King Island scrubtit Acanthornis magna greeniana is a critically endangered subspecies of the nominate Tasmanian scrubtit A. m. magna, with an estimated population of < 100 individuals persisting in three patches of swamp forest. The Tasmanian scrubtit is widespread in wet forests on mainland Tasmania. We sequenced the scrubtit genome using PacBio HiFi and undertook a population genomic study of the King Island and Tasmanian scrubtits using a double-digest restriction site-associated DNA (ddRAD) dataset of 5,239 SNP loci. The genome was 1.48 Gb long, comprising 1,518 contigs with an N50 of 7.715Mb. King Island scrubtits formed one of four overall genetic clusters, but separated into three distinct subpopulations when analysed independently of the Tasmanian scrubtit. Pairwise FST values were greater among the King Island scrubtit subpopulations than among most Tasmanian scrubtit subpopulations. Genetic diversity was lower and inbreeding coefficients were higher in the King Island scrubtit than all except one of the Tasmanian scrubtit subpopulations. We observed crown baldness in 8/15 King Island scrubtits, but 0/55 Tasmanian scrubtits. Six loci were significantly associated with baldness, including one within the DOCK11 gene which is linked to early feather development. Contemporary gene flow between King Island scrubtit subpopulations is unlikely, with further field monitoring required to quantify the fitness consequences of its small population size, low genetic diversity and high inbreeding. Evidence-based conservation actions can then be implemented before the taxon goes extinct.

Whole-genome selection signature differences between Chaohu and Ji'an red ducks.

Assessing the genetic structure of local varieties and understanding their genetic data are crucial for effective management and preservation. However, the genetic differences among local breeds require further explanation. To enhance our understanding of their population structure and genetic diversity, we conducted a genome-wide comparative study of Chaohu and Ji'an Red ducks using genome sequence and restriction site-associated DNA sequencing technology. Our analysis revealed a distinct genetic distinction between the two breeds, leading to divided groups. The phylogenetic tree for Chaohu duck displayed two branches, potentially indicating minimal impact from artificial selection. Additionally, our ROH (runs of homozygosity) analysis revealed that Chaohu ducks had a lower average inbreeding coefficient than Ji'an Red ducks. We identified several genomic regions with high genetic similarity in these indigenous duck breeds. By conducting a selective sweep analysis, we identified 574 candidate genes associated with muscle growth (BMP2, ITGA8, MYLK, and PTCH1), fat deposits (ELOVL1 and HACD2), and pigmentation (ASIP and LOC101797494). These results offer valuable insights for the further enhancement and conservation of Chinese indigenous duck breeds.

Functional role of genes in ROH sites in &lt;i&gt;Czech Golden&lt;/i&gt; chickens

Modern technologies make it possible to identify homozygous chromosome regions that have arisen as a result of animal selection. The article presents the results obtained in the genotyping of chickens using the Illumina Chicken 60KSNP iSelect Bead Chip. For the first time, a genome-wide analysis of the extended homozygous SNP sequences (ROH) was performed in the genome of Czech Golden hens. The average number of ROH segments in the chicken genome was 143 ± 8. ROH segments are usually randomly distributed in chicken chromosomes. It has been proposed to prohibit the use of heterozygous SNPs in ROH segments to prevent overestimation of ROH data. The average inbreeding coefficient in chickens calculated from ROH data was 0.34 ± 0.03. ROH islands were found on chromosomes 2, 3, 9 and 22.They contain genesassociatedwithimmunity, degradation of neurons, Fabricius bursa weight, obesity, feather pigmentation, and regulation of hit-shock genes. In the Czech Golden chicken breed, selection and the associated inbreeding have therefore influenced the genes involved in several biological processes.

Whole-genome resequencing reveals genomic variation and dynamics in Ethiopian indigenous goats

Ethiopia has about 52 million indigenous goats with marked phenotypic variability, which is the outcome of natural and artificial selection. Here, we obtained whole-genome sequence data of three Ethiopian indigenous goat populations (Arab, Fellata, and Oromo) from northwestern Ethiopia and analyzed their genome-wide genetic diversity, population structure, and signatures of selection. We included genotype data from four other Ethiopian goat populations (Abergelle, Keffa, Gumuz, and Woyto-Guji) and goats from Asia; Europe; and eastern, southern, western, and northern Africa to investigate the genetic predisposition of the three Ethiopian populations and performed comparative genomic analysis. Genetic diversity analysis showed that Fellata goats exhibited the lowest heterozygosity values (Ho = 0.288 ± 0.005 and He = 0.334 ± 0.0001). The highest values were observed in Arab goats (Ho = 0.310 ± 0.010 and He = 0.347 ± 4.35e−05). A higher inbreeding coefficient (FROH = 0.137 ± 0.016) was recorded for Fellata goats than the 0.105 ± 0.030 recorded for Arab and the 0.112 ± 0.034 recorded for Oromo goats. This indicates that the Fellata goat population should be prioritized in future conservation activities. The three goat populations showed the majority (∼63%) of runs of homozygosity in the shorter (100–150 Kb) length category, illustrating ancient inbreeding and/or small founder effects. Population relationship and structure analysis separated the Ethiopian indigenous goats into two distinct genetic clusters lacking phylogeographic structure. Arab, Fellata, Oromo, Abergelle, and Keffa represented one genetic cluster. Gumuz and Woyto-Guji formed a separate cluster and shared a common genetic background with the Kenyan Boran goat. Genome-wide selection signature analysis identified nine strongest regions spanning 163 genes influencing adaptation to arid and semi-arid environments (HOXC12, HOXC13, HOXC4, HOXC6, and HOXC9, MAPK8IP2), immune response (IL18, TYK2, ICAM3, ADGRG1, and ADGRG3), and production and reproduction (RARG and DNMT1). Our results provide insights into a thorough understanding of genetic architecture underlying selection signatures in Ethiopian indigenous goats in a semi-arid tropical environment and deliver valuable information for goat genetic improvement, conservation strategy, genome-wide association study, and marker-assisted breeding.

Analysis of genetic diversity and structure of endangered Dengchuan cattle population using a single-nucleotide polymorphism chip

This study aimed to evaluate the genetic diversity and structure within the Dengchuan cattle population and effectively protect and utilize their germplasm resources. Herein, the single-nucleotide polymorphisms (SNPs) of 100 Dengchuan cattle (46 bulls and 54 cows) were determined using the GGP Bovine 100K SNP Beadchip. The results showed that among the Dengchuan cattle, a total of 101,220 SNPs were detected, and there were 83,534 SNPs that passed quality control, of which 85.7% were polymorphic. The average genetic distance based on identity-by-state (IBS) within the conservation population of Dengchuan cattle was 0.26 ± 0.02. A total of 3,999 genome-length runs of homozygosity (ROHs) were detected in the Dengchuan cattle, with ROH lengths primarily concentrated in the range of 1–5 Mb, accounting for 87.02% of the total. The average inbreeding coefficient based on ROHs was 4.6%, within the conservation population of Dengchuan cattle, whereas it was 4.9% for bulls, and the Wright inbreeding coefficient (FIS) value was 2.4%, demonstrating a low level of inbreeding within the Dengchuan cattle population. Based on neighbor-joining tree analysis, the Dengchuan cattle could be divided into 16 families. In summary, the conservation population of Dengchuan cattle displays relatively abundant diversity and a moderate genetic relationship. Inbreeding was observed among a few individuals, but the overall inbreeding level of the population remained low. It is important to maintain this low level of inbreeding when introducing purebred bloodlines to expand the core group. This approach will ensure the long-term conservation of Dengchuan cattle germplasm resources and prevent loss of genetic diversity.

Genetic analysis of a Kaijiang duck conservation population through genome-wide scan

ABSTRACT 1. The Kaijiang duck is a native Chinese breed known for its excellent egg laying performance, killing-out percentage (88.57%), and disease resistance. The assessment of population genetic structure is the basis for understanding the genetics of indigenous breeds and for their protection and management. 2. In this study, whole-genome sequencing was performed on 60 Kaijiang ducks to identify genetic variations and investigate the population structure. Homozygosity (ROH) analysis was conducted to assess inbreeding levels in the population. 3. The study revealed a moderate level of inbreeding, indicated by an average inbreeding coefficient of 0.1043. This may impact the overall genetic diversity. 4. Genomic Regions of Interest identified included 168 genomic regions exhibiting high levels of autozygosity. These regions were associated with processes including muscle growth, pigmentation, neuromodulation, and growth and reproduction. 5. The significance of these pathways indicated their potential role in shaping the desirable traits of the Kaijiang duck. These findings provide insights into the genetic basis of the Kaijiang duck’s desirable traits and can inform future breeding and conservation efforts.

Genetic diversity and population structure of four Anatolian sheep revealed by genome-wide ddRADseq data

This is the first wide-scale study to reveal next-generation sequencing-based and comprehensive study to reveal genetic diversity, population structure, and phylogenetic relationships of some of Türkiye's indigenous fat-tailed sheep breeds. A total of 80 animals belonging to Akkaraman (AKR), Güney Karaman (GKR), Morkaraman (MKR), and Karakaş (KRK) sheep breeds were genotyped with genome-wide of 296.097 single nucleotide polymorphisms (SNPs) data recovered from double digest restriction site-associated DNA sequencing (ddRADseq) libraries processed by Illumina HiSeq X Ten instrument. Minor allele frequency (MAF) ranged from 0.311 (MKR) to 0.316 (AKR and GKR). Observed heterozygosity (HO) values were higher than expected values (HE) across all populations with a mean of 0.307. Negative inbreeding coefficient (FIS) values were observed in four sheep populations varying between −0.034 (AKR) to −0.060 (MKR). Cross-entropy criterion-based ADMIXTURE analysis demonstrated that four Anatolian sheep breeds were represented by two ancestral populations in which GKR and AKR were distinct, while KRK and MKR were mixed populations. Genetic distance-based tree, discriminant analysis of principal components (DAPC), and ADMIXTURE analyses confirmed that GKR, which is believed to be derived from the AKR breed, has become genetically different from other sheep populations. Moreover, no migration from and into the GKR breed was detected, while a migration edge was drawn from the KRK and MKR clade into the AKR breed by the TreeMix algorithm. This study confirmed that high-density SNP data obtained from advanced next-generation sequencing (NGS) platforms are efficient in genetically distinguishing animal breeds, even varieties. It is recommended that further NGS-based studies could be carried out to better characterize different indigenous livestock populations raised in Türkiye.

Genetic Diversity and Population Structure among Arabian Horse Genealogical Lineages in Bulgaria

The present research aimed to characterize the genetic diversity and relationships among extant Arabian horse sire lines in Bulgaria, using 15 equine microsatellite markers. The evaluation included 537 Arabian horses representing nine sire lines (SAKLAWI I, LATIF, SEANDERICH, IBRAHIM, SHABAB, DJEBEL MOUSA, KUHAILAN AFAS, BAIRACTAR, and SARHAN). The obtained results indicated that within these lines, the mean number of alleles ranged from 4.15 in SARHAN to 5.54 in SAKLAWI I and LATIF. The mean expected heterozygosity (He) ranged from 0.54 in the SEANDERICH line to 0.67 in SAKLAWI I. The inbreeding coefficient for the entire Arabian populations was rather low: FIS = −0.109, fluctuating from −0.204 in SHABAB to −0.041 in SAKLAWI I. The mean genetic differentiation, FST, was 0.096, demonstrating that nearly 90% of the total genetic variation was due to genetic differentiation within each population. STRUCTURE analysis indicated a genetic similarity between SHABAB and LATIF, between IBRAHIM and KUHAILAN AFAS, as well as between SAKLAWI I, SEANDERICH, and BAIRACTAR. This study of the genetic diversity of Arabian sire lines in Bulgaria can assist in developing a national strategy for the exclusion of non-purebred animals from breeding programs in order to preserve the genetic profile of the original Arabian lines.

Genetic polymorphism and evidence of signatures of selection in the Plasmodium falciparum circumsporozoite protein gene in Tanzanian regions with different malaria endemicity

BackgroundIn 2021 and 2023, the World Health Organization approved RTS,S/AS01 and R21/Matrix M malaria vaccines, respectively, for routine immunization of children in African countries with moderate to high transmission. These vaccines are made of Plasmodium falciparum circumsporozoite protein (PfCSP), but polymorphisms in the gene raise concerns regarding strain-specific responses and the long-term efficacy of these vaccines. This study assessed the Pfcsp genetic diversity, population structure and signatures of selection among parasites from areas of different malaria transmission intensities in Mainland Tanzania, to generate baseline data before the introduction of the malaria vaccines in the country.MethodsThe analysis involved 589 whole genome sequences generated by and as part of the MalariaGEN Community Project. The samples were collected between 2013 and January 2015 from five regions of Mainland Tanzania: Morogoro and Tanga (Muheza) (moderate transmission areas), and Kagera (Muleba), Lindi (Nachingwea), and Kigoma (Ujiji) (high transmission areas). Wright’s inbreeding coefficient (Fws), Wright’s fixation index (FST), principal component analysis, nucleotide diversity, and Tajima’s D were used to assess within-host parasite diversity, population structure and natural selection.ResultsBased on Fws (< 0.95), there was high polyclonality (ranging from 69.23% in Nachingwea to 56.9% in Muheza). No population structure was detected in the Pfcsp gene in the five regions (mean FST = 0.0068). The average nucleotide diversity (π), nucleotide differentiation (K) and haplotype diversity (Hd) in the five regions were 4.19, 0.973 and 0.0035, respectively. The C-terminal region of Pfcsp showed high nucleotide diversity at Th2R and Th3R regions. Positive values for the Tajima’s D were observed in the Th2R and Th3R regions consistent with balancing selection. The Pfcsp C-terminal sequences revealed 50 different haplotypes (H_1 to H_50), with only 2% of sequences matching the 3D7 strain haplotype (H_50). Conversely, with the NF54 strain, the Pfcsp C-terminal sequences revealed 49 different haplotypes (H_1 to H_49), with only 0.4% of the sequences matching the NF54 strain (Hap_49).ConclusionsThe findings demonstrate high diversity of the Pfcsp gene with limited population differentiation. The Pfcsp gene showed positive Tajima’s D values, consistent with balancing selection for variants within Th2R and Th3R regions. The study observed differences between the intended haplotypes incorporated into the design of RTS,S and R21 vaccines and those present in natural parasite populations. Therefore, additional research is warranted, incorporating other regions and more recent data to comprehensively assess trends in genetic diversity within this important gene. Such insights will inform the choice of alleles to be included in the future vaccines.

Runs of homozygosity analysis for selection signatures in the Yellow Korean native chicken.

Yellow Korean native chicken (KNC-Y) is one of the five pure Korean indigenous chicken breeds that were restored through a government project in 1992. KNC-Y is recognized for its superior egg production performance compared to other KNC lines. In this study, we performed runs of homozygosity (ROH) analysis to discover selection signatures associated with egg production traits in the KNC-Y population. A total of 675 DNA samples from KNC-Y were genotyped to generate single nucleotide polymorphism (SNP) data using custom 60K Affymetrix SNP chips. ROH analysis was performed using PLINK software, with predefined parameters set for the analysis. The threshold of ROH island was defined as the top 1% frequency of SNPs withing the ROH among the population. In the KNC-Y population, a total of 29,958 runs of homozygosity (ROH) fragments were identified. The average total length of ROH was 120.84 Mb, with each ROH fragment having an average length of 2.71 Mb. The calculated ROH-based inbreeding coefficient (FROH) was 0.13. Furthermore, we revealed the presence of ROH islands on chromosomes 1, 2, 4, 5, 7, 8, and 11. Within the identified regions, a total of 111 genes were annotated, and among them were genes related to economic traits, including PRMT3, ANO5, HDAC4, LSS, PLA2G4A, and PTGS2. Most of the overlapping quantitative trait locus regions with ROH islands were found to be associated with production traits. This study conducted a comprehensive analysis of ROH in the KNC-Y population. Notably, among the findings, the PTGS2 gene is believed to play a crucial role in influencing the laying performance of KNC-Y.

Genealogical structure and genetic variability: A step towards conservation of white Morada Nova sheep in semi-arid region

Morada Nova breed has a low effective herd, and its white variety is in risk of extinction in Brazilian semi-arid region. This study aimed to evaluate white Morada Nova sheep’s genetic variability and population structure in four flocks. These are the only flocks in the world. The assessed parameters of pedigree integrity and a genealogy record of 219 animals. 44.74 % did not have pedigree information on their mother or father, and 50.47 % did not have available information on their grandparents. 52.97 % had a pedigree in their first ancestry (obtained through the average of known parents), 23.74 % in the second, 13.58 % in the third, and 5.33 % in the fourth. The effective size was 10.7, considering the complete equivalent generation. The effective number of founding animals (fe) and ancestors (fa) was 10, indicating the effective participation of all founding animals in the population over the generations, reducing losses of genetic variability. Out of all the ancestors, only 3 explained 50 % of the population’s genetic variability, reflecting the low values of fa and fe. The population’s average inbreeding and kinship coefficients were 2.61 and 4.54. The F statistics of Wright, Fis, Fst, and Fit, were −0.012, 0.016, and 0.003. The genetic conservation index (GCI) had an average of 2.01 ± 1.29, with minimum and maximum values of 1.0 and 6.38, respectively. Approximately, 71 % of the population had GCI below 2.0, while only 5.9 % of sheep had an GCI bigger than 5.0. This indicates the absence of sub-structures inf the population. There is no population subdivision, considering the genealogical structure of the evaluated flocks. The inbreeding and kinship coefficient values are controlled and serve as a starting point for establishing a future management plan for the genetic variability of the flock. It is recommended to use strategies to increase the effective population size, to meet the minimum number recommended by FAO and to promote the conservation and genetic improvement of the breed.

225 Genomic diversity of Katahdin sheep and impacts on genomic prediction

Abstract In the late 1950s, Katahdin hair sheep were developed as a composite breed of medium size and moderate prolificacy, with potential to express resistance to gastrointestinal nematodes. With its increasing popularity and the recent adoption of genomic prediction in the Katahdin breeding program, evaluating the genetic diversity in this breed is necessary. Prediction accuracies rely on the strength of genetic relationships. That strength often increases with selection using genomically-enhanced estimated breeding values; however, it carries the risk of decreasing diversity. We aim to characterize the genetic diversity present in the breed and the subsequent impact on the accuracies of genomic predictions. Genotypes on Katahdin sheep from member flocks in the National Sheep Improvement Program were used. Data from two medium density 50k bead chips and one high density 600k bead chip were merged and limited to shared autosomal single nucleotide polymorphisms (SNP). After quality control, there were 9,705 animals and 31,984 SNP remaining. Only for analyses of inbreeding and heterozygosity, linkage disequilibrium (LD) pruning was performed; variants with an LD score exceeding 0.8 were removed with 29,904 SNP remaining. To determine the extent of population structure, principal component (PC) analysis was conducted with PC one and PC two explaining 4.9% and 3.8% of the variation in the population, respectively. When examining subpopulations, a single ancestral population explained 99.9% of the genetic variation among animals. When clustered by identity-by-state, all animals were placed in a single cluster. Inbreeding coefficients were estimated with four methods: method of moments, variance-standardized relationship minus one, excess homozygosity, and correlation between uniting gametes. Across methods, average inbreeding coefficients ranged from 1.14% to 1.20% with Pearson and Spearman correlations ranging from -0.51 to 0.91 (Table 1). Historical effective population sizes (Ne) were estimated for earlier generations, and a linear regression fit to predict more recent Ne. At the founding of the breed, historical Ne was estimated to be 310, and the current Ne was estimated to be 150. Effects of Ne on the accuracy of genomic predictions were examined for reference population sizes ranging from 1k to 25k for traits of low (0.1), moderate (0.3), and high (0.5) heritabilities. The corresponding number of effective chromosomal segments were 2NeL, where L was the length of the genome in morgans (26M). Accuracies of genomic predictions increased as the reference population size grew, with only marginal gains once the population exceeded 15k. Overall, there were no detectable genetic subpopulations in the Katahdin breed. Current levels of inbreeding were low. With a small Ne the breed can achieve high accuracy for genomic predictions by aiming for a reference population of 15k animals. Negative impacts on genetic diversity due to genomic selection are not of immediate concern.

Detecting inbreeding depression in structured populations

Measuring inbreeding and its consequences on fitness is central for many areas in biology including human genetics and the conservation of endangered species. However, there is no consensus on the best method, neither for quantification of inbreeding itself nor for the model to estimate its effect on specific traits. We simulated traits based on simulated genomes from a large pedigree and empirical whole-genome sequences of human data from populations with various sizes and structures (from the 1,000 Genomes project). We compare the ability of various inbreeding coefficients ([Formula: see text]) to quantify the strength of inbreeding depression: allele-sharing, two versions of the correlation of uniting gametes which differ in the weight they attribute to each locus and two identical-by-descent segments-based estimators. We also compare two models: the standard linear model and a linear mixed model (LMM) including a genetic relatedness matrix (GRM) as random effect to account for the nonindependence of observations. We find LMMs give better results in scenarios with population or family structure. Within the LMM, we compare three different GRMs and show that in homogeneous populations, there is little difference among the different [Formula: see text] and GRM for inbreeding depression quantification. However, as soon as a strong population or family structure is present, the strength of inbreeding depression can be most efficiently estimated only if i) the phenotypes are regressed on [Formula: see text] based on a weighted version of the correlation of uniting gametes, giving more weight to common alleles and ii) with the GRM obtained from an allele-sharing relatedness estimator.

Analysis of the Genetic Diversity and Mating System of the Endangered Plant Keteleeria davidiana var. calcarea

Keteleeria davidiana var. calcarea is an endangered plant endemic to China. It is mainly distributed in the karst areas of the Guangxi and Guizhou Provinces. It is characterized by small populations and intermittent distribution. This study aimed to explore the genetic diversity and mating system of wild populations of Keteleeria davidiana var. calcarea in fragmented habitats. To achieve this, we genotyped 46 maternal trees and 214 progenies from four fragmented populations of Keteleeria davidiana var. calcarea using nine pairs of microsatellite primers. The genetic diversity of Keteleeria davidiana var. calcarea (Ho = 0.68, He = 0.63) was lower than that of the species overall but higher than that of other Keteleeria plants. The incidence of unbiased expected heterozygosity (uHe) and allelic richness (Ar) was higher in the maternal generation than in the progeny. This suggests that the genetic diversity of the progeny was lower than that of the maternal generation. Keteleeria davidiana var. calcarea is divided into four populations, but there is significant genetic exchange between the populations according to STRUCTURE and gene flow analyses. The multilocus mating system analysis (MLTR) results indicate that the multilocus outcrossing rate (tm) was 0.902, the single-locus outcrossing rate (ts) was 0.606, the bimaternal inbreeding coefficient (tm-ts) was 0.295, and the coefficient of inbreeding depression(δ) was 0.904. These results suggest a certain degree of selfing and inbreeding in Keteleeria davidiana var. calcarea. To prevent problems associated with inbreeding and conserve the genetic diversity of Keteleeria davidiana var. calcarea, we recommend establishing seed gardens, using artificial pollination, and employing asexual propagation techniques for conservation intervention.

Analysis of the Genetic Relationship and Inbreeding Coefficient of the Hetian Qing Donkey through a Simplified Genome Sequencing Technology.

The Hetian Qing donkey is an excellent local donkey breed in Xinjiang. It is of great significance to accelerate breeding and the speed of breeding and rejuvenation, as well as to understand the genetic basis of the strategies and population. This study collected a total of 4 male donkeys and 28 female donkeys. It then obtained genotype data through Simplified Genomic Sequencing (GBS) technology for data analysis. The results detected a total of 55,399 SNP loci, and the genotype detection rate of individuals was ≥90%. A total of 45,557 SNP loci were identified through quality control, of which 95.5% were polymorphic. The average minimum allele frequency was 0.250. The average observed heterozygosity was 0.347. The average expected heterozygosity was 0.340. The average IBS (state homologous) genetic distance was 0.268. ROH: 49 (homozygous fragments), with 73.47% of the length between 1 and 5 Mb. The average per-strip ROH length was 1.75 Mb. The mean inbreeding coefficient was 0.003. The 32 Hetian green donkeys could be divided into six families. The number of individuals in each family is significant. To sum up, the Hetian Qing donkey population has low heterozygosity, few families, and large differences in the number of individuals in each family, which can easily cause a loss of genetic diversity. In the subsequent process of seed protection, seed selection should be conducted according to the divided pedigree to ensure the long-term protection of the genetic resources of Hetian green donkeys.

Benchmarking kinship estimation tools for ancient genomes using pedigree simulations.

There is growing interest in uncovering genetic kinship patterns in past societies using low-coverage palaeogenomes. Here, we benchmark four tools for kinship estimation with such data: lcMLkin, NgsRelate, KIN, and READ, which differ in their input, IBD estimation methods, and statistical approaches. We used pedigree and ancient genome sequence simulations to evaluate these tools when only a limited number (1 to 50 K, with minor allele frequency ≥0.01) of shared SNPs are available. The performance of all four tools was comparable using ≥20 K SNPs. We found that first-degree related pairs can be accurately classified even with 1 K SNPs, with 85% F1 scores using READ and 96% using NgsRelate or lcMLkin. Distinguishing third-degree relatives from unrelated pairs or second-degree relatives was also possible with high accuracy (F1 > 90%) with 5 K SNPs using NgsRelate and lcMLkin, while READ and KIN showed lower success (69 and 79% respectively). Meanwhile, noise in population allele frequencies and inbreeding (first-cousin mating) led to deviations in kinship coefficients, with different sensitivities across tools. We conclude that using multiple tools in parallel might be an effective approach to achieve robust estimates on ultra-low-coverage genomes.

Characterisation of runs of homozygosity and inbreeding coefficients in the red-brown Korean native chickens.

The analysis of runs of homozygosity (ROH) has been applied to assess the level of inbreeding and identify selection signatures in various livestock species. The objectives of this study were to characterize the ROH pattern, estimate the rate of inbreeding, and identify signatures of selection in the red-brown Korean native chickens. The Illumina 60K single nucleotide polymorphism chip data of 651 chickens was used in the analysis. Runs of homozygosity were analysed using the PLINK v1.9 software. Inbreeding coefficients were estimated using the GCTA software and their correlations were examined. Genomic regions with high levels of ROH were explored to identify selection signatures. A total of 32,176 ROH segments were detected in this study. The majority of the ROH segments were shorter than 4 Mb. The average ROH inbreeding coefficients (FROH) varied with the length of ROH segments. The means of inbreeding coefficients calculated from different methods were also variable. The correlations between different inbreeding coefficients were positive and highly variable (r = 0.18-1). Five ROH islands harbouring important quantitative trait loci were identified. This study assessed the level of inbreeding and patterns of homozygosity in Red-brown native Korean chickens. The results of this study suggest that the level of recent inbreeding is low which indicates substantial progress in the conservation of red-brown Korean native chickens. Additionally, Candidate genomic regions associated with important production traits were detected in homozygous regions.

Genome-wide single nucleotide polymorphisms reveal the genetic diversity and population structure of Creole goats from northern Peru

Goat farming constitutes a significant source of income for farmers in northern Peru. There is currently an absence of information about the genetics of Peruvian Creole goats that would enable us to understand their origins and genetic spread. The objective of this study was to estimate the genetic diversity of Creole goats from northern Peru using SNP markers. This study involved the collection of 192 male Creole goats from three key goat production geographical departments in northern Peru. These goat samples were genotyped using the GGPGoat70k SNP panel. To explore the genetic influence of other breeds on Peruvian Creole goats, our dataset was combined with previously published SNP genotypes. External data set includes multiple breeds genotypes sampled from Argentina, Brazil, Spain, and Alpine breed from Italy, France, and Switzerland. After quality control 52,832 autosomal SNPs were used to assess genetic diversity in the Peruvian goats. For the population structure analysis of the merged data 20,513 common SNPs were used. Estimations for expected heterozygosity (He), observed heterozygosity (Ho), and inbreeding coefficient (FIS) were computed for the Peruvian groups. AMOVA, principal component analysis and ADMIXTURE were conducted to evaluate the population structure in the two data sets, Peru and merged. The results revealed a considerable genetic diversity, with Ho values ranging from 0.40 to 0.41 for the Peruvian sampling groups, and inbreeding coefficient was notably low for Peruvian goat. The population structure analysis demonstrated a distinction (p < 0.05) from other breeds. These findings suggest a level of genetic differentiation of the Peruvian goat population among other breeds, although further research is needed considering samples from other Peruvian areas. We expect this study will contribute to define genetic management strategies to prevent the loss of genetic diversity in Peruvian goat populations and for upcoming advancements in this field.