Population Genetic Research Articles

Validation of the Amharic version of Cognitive Assessment Interview (CAI-A) in people with schizophrenia in Ethiopia

Assessing cognition with interview-based measures could be a low-resource alternative to traditional cognitive tests. We previously adapted the Cognitive Assessment Interview (CAI) into Amharic (CAI-A) for use with people with schizophrenia (PWS) in Ethiopia. This study examined the convergent and structural validity of the CAI-A in a group of 350 PWS sub-sampled from the Neuropsychiatric Genetics of African Populations – Psychosis (NeuroGAP-Psychosis) study, who fulfilled the inclusion criteria. Data were analyzed using confirmatory factor analysis (CFA), Spearman's correlation coefficient (ρ), multiple regression, and Item Response Theory (IRT). A one-factor solution best fits the items in the tool (factor loadings between 0.58 and 0.79), suggesting structural validity. The total score of the CAI-A moderately correlated with functioning (ρ = 0.44, p < 0.001) and symptom dimensions (ρ between 0.38 and 0.46, p < 0.001), suggesting convergent validity. Multiple regression showed that age (β = −0.06, 95 % CI (−0.12, 0.00), p = 0.044), the duration of illness (β = 0.08, 95 % CI (0.01, 0.14), p = 0.033), and medication side effects (β = 0.35, 95 % CI (0.21, 0.50), p < 0.001) were positively and significantly associated with the CAI-A total score. The IRT analysis suggested that the tool best functions among participants with moderate to severe impairment (difficulty coefficient between 0.05 and 2.73). We found that the CAI-A is a valid tool for use in Ethiopia. The moderate correlation with symptom and functional measures suggested that self-reported cognitive symptoms parallel other symptom dimensions and functional disability. The CAI-A can be used in clinical practice and research activities in PWS in Ethiopia when subjective assessment of cognition is desired.

Genomic insights into the diversity of mango cultivars in Saudi Arabia: Fluorescent SSR marker analysis

ObjectivesThis study presents a comprehensive genetic analysis of twelve mango cultivars from Saudi Arabia, investigating their genetic diversity, population structure, and potential loci under selection. MethodsUtilizing ten SSR loci, we conducted a thorough assessment of allele variation, revealing diverse patterns across loci and defined genetic diversity among cultivars. ResultsNotably, certain loci exhibited extensive allele diversity, with Loci M02 and M06 displaying five and seven alleles, respectively, while others showed more limited profiles, such as Locus M01 with only two alleles. Surprisingly, clustering patterns based on genetic relatedness did not consistently align with anticipated geographic groupings, suggesting additional factors influencing genetic structure beyond geography. Further analysis of allele sharing provided insights into the genetic relatedness among cultivars, highlighting both close genetic affinities and significant divergence. Subsequent Allele Sharing Matrix analysis unveiled sub-clustering within clusters, revealing finer-scale genetic structures within the population. Assessment of genetic differentiation using Analysis of Molecular Variance demonstrated significant variations among clusters and individuals, with variations among clusters ranging from 3% to 53%. Additionally, testing for loci under selection identified potential markers, such as locus M11, with observed heterozygosity greater than expected (1.024) and a relatively low observed FST value of 0.2695. ConclusionsOverall, our findings contribute to the characterization of mango genetic resources in Saudi Arabia, providing crucial information for breeding programs, germplasm conservation efforts, and agricultural practices. By leveraging advanced genomic technologies and interdisciplinary collaborations, future research can delve deeper into the genetic basis of adaptive traits and agronomic characteristics, paving the way for innovative strategies to enhance mango production and sustainability in Saudi Arabia.

Forward and backward modeling of cultural evolutionary processes

Modern cultural evolution theory adopts a variety of concepts and methods developed in mathematical biology, in particular population genetics theory. In addition to forward-looking approaches such as two-locus models, backward-looking approaches such as coalescent theory, which describe ancestral states of the current population, have played an important role in population genetics. Here, we show how forward and backward approaches can be applied to two examples in cultural evolution. The first example deals with the number of cultural traits, which can be analyzed by a backward approach. Cultural coalescent theory illustrates many unique aspects of cultural processes such as information transfer (i.e., social learning) from many donors (“cultural parents”), including or excluding the biological parents. Theory predicts that many cultural traits of intermediate (or higher) frequency (popularity) can exist that are surprisingly old. Many unsolved issues remain, however, such as how to incorporate social structure as well as natural and/or cultural selection, which we believe to be an urgent agenda. The second example is the punishment of sibling incest, a problem to which cultural coalescent theory cannot currently be usefully applied. By analyzing forward recursions, we show that punishment is ineffective in suppressing incest, unless the incestuous inclination is diminished by the presence of punishers. Based on these examples, we discuss the merits and demerits of forward and backward approaches.

Models of Fluctuating Selection Between Generations: A Solution for the Theoretical Inconsistency.

The theory of selection fluctuation between generations has been a topic with much activities in population genetics and molecular evolution in 1970's. Most studies suggested that, as the result of fluctuating selection between generations, the frequency of an (on average) neutral mutation may fluctuate around 0.5 during the long-term evolution before it was ultimately fixed or lost. However, this pattern can only be derived from a specific type Wright-Fisher additive model, coined by the Nei-Yokoyama puzzle. In this commentary, I revisited this issue and figured out a theoretical assumption that has never been claimed explicitly, the notion of reference phenotype. Consider one locus with two-alleles: A is the wildtype allele and A' is the mutation. The fluctuating selection model actually requires a constraint that one of three genotypes (AA, AA', or A'A') must maintain a constant fitness without fluctuating between generations. It appears that the balancing selection at a frequency of 0.5 emerges only when the heterozygote (AA') is the reference genotype. Because it is difficult to determine which genotype could be the reference genotype in a real population, a desirable population genetics model should take all three possibilities into account. To this end, I propose a mixture model, where each genotype has a certain chance to be the reference genotype. My analysis showed that the emergence of balancing selection depends on the relative proportions of three different reference genotypes.

Analysis of genetic diversity and population structure of Iranian Artemisia annua L. accessions using IRAP-based molecular markers and revealing retroelement insertional polymorphism

Analysis of genetic diversity and population structure of Iranian Artemisia annua L. accessions using IRAP-based molecular markers and revealing retroelement insertional polymorphism

Genomic Single Nucleotide Polymorphism (SNP) markers and mitochondrial haplotypes illuminate the origins of Crown-of-Thorns Starfish (Acanthaster solaris) outbreaks in the South China Sea.

Outbreaks of the coral predator Crown-of-Thorns Starfish (CoTS) pose a severe threat to coral reefs in the Indo-Pacific Ocean. In 2018, the South China Sea (SCS) experienced significant CoTS outbreaks, leading to extensive coral mortality across the Xisha, Zhongsha, Dongsha, and Nansha Islands, severely impacting the coral reef ecosystem. To explore the origins of these outbreaks, we conducted a comprehensive genomic analysis using data from genomic single nucleotide polymorphism sites (SNPs) and mitochondrial haplotypes. Our analysis reveals that CoTS populations in the SCS, which exhibit moderate genetic diversity and may have undergone positive selection or population expansion. There was limited genetic differentiation among CoTS populations from XS, ZS, and NS groups. Especially between the XS and ZS groups, there was almost no genetic differentiation. The populations from XS, ZS, and NS groups have strong genetic connections with populations in Vietnam and the Philippines. There was high gene flow from Vietnam to the Xisha Islands and from the Philippines to the Nansha Islands, suggesting that the CoTS populations in these regions primarily originate from these neighboring countries. The comprehensive analyses of SNP and mitochondrial genomes have provided valuable insights into the population genetics of CoTS. This research has generated significant genomic resources and facilitated important studies on the genetics of the CoTS species. By identifying potential source populations and understanding the genetic basis of their spread, managers can develop more effective conservation strategies to protect vulnerable coral reef ecosystems in the SCS.

Modeling the velocity of evolving lineages and predicting dispersal patterns

Accurate estimation of the dispersal velocity or speed of evolving organisms is no mean feat. In fact, existing probabilistic models in phylogeography or spatial population genetics generally do not provide an adequate framework to define velocity in a relevant manner. For instance, the very concept of instantaneous speed simply does not exist under one of the most popular approaches that models the evolution of spatial coordinates as Brownian trajectories running along a phylogeny. Here, we introduce a family of models-the so-called Phylogenetic Integrated Velocity (PIV) models-that use Gaussian processes to explicitly model the velocity of evolving lineages instead of focusing on the fluctuation of spatial coordinates over time. We describe the properties of these models and show an increased accuracy of velocity estimates compared to previous approaches. Analyses of West Nile virus data in the United States indicate that PIV models provide sensible predictions of the dispersal of evolving pathogens at a one-year time horizon. These results demonstrate the feasibility and relevance of predictive phylogeography in monitoring epidemics in time and space.

Extent and patterns of morphological and molecular genetic diversity and population structure of Nigerian Taro cultivars.

Genetic diversity is crucial for conservation efforts as well as breeding programs targeted at the development of improved varieties. Taro, a climate-resilient crop, plays a vital role in the nutritional and economic livelihoods of many households in Nigeria, but its yield is very low due to inadequate genetic improvement efforts. A diversity assessment of Nigerian taro is therefore required to create a premise for its improvement in yield, quality and disease tolerance. In this study, the genetic diversity and population structure of 490 taro cultivars comprising two main gene pools: Dasheen (215) and Eddoe (275), collected from farmers and marketers across seven states in Nigeria was assessed using 3047 Diversity Array Technology single nucleotide polymorphism (DArT-SNP) markers. A subset of 114 taro cultivars, comprising 30 Dasheens and 84 Eddoes were further phenotyped using 24 agro-morphological descriptors. Both phenotypic and molecular characterization revealed higher genetic diversity among the Eddoes than Dasheens. Estimates of gene flow (Nm = 0.353) revealed intermixing of cultivars among the States of collection, with the highest gene flow occurring between cultivars from Anambra and Ondo states and the lowest between Anambra and Kwara states. Population structure and Ward's minimum variance hierarchical cluster based on DArT-SNPs identified four groups, one comprising Dasheen and three comprising Eddoe cultivars. Hierarchical clustering based on phenotypic traits delineated three clusters. Variation between gene pools (49%) was higher than within gene pools (32%). Variation among States of collection was high (41%), while variation among individuals within gene pools (18%) and States of collection (19%) was relatively low. Correlation between phenotypic and genotypic diversity assessments was low (r = 0.01), indicating that both approaches were necessary for assessing genetic diversity in taro. However, genotypic assessment provided better information about genetic diversity of the taro cultivars. This is the first study that represented germplasm collection across the major taro growing regions of Nigeria. The findings from this study based on agro-morphological characterization and DArT-SNP genotyping are critical for genetic characterization, conservation and breeding of taro in Nigeria, mainly initiating hybridization between the two genepools after careful assessment of ploidy levels of the accessions collected in this study. This will facilitate in developing improved taro varieties with desirable traits, such as higher yield, better disease resistance, and improved nutritional quality.

A Comparative Metabolomics Study of the Potential Marker Compounds in Feces from Different Hybrid Offspring of Huainan Pigs

As a notable native Chinese genetic population, the Huainan pig has an exceptional meat quality but a low percentage of lean meat and subpar genetic performance. To better exploit the superior genetic traits of the Huainan pig and address knowledge gaps regarding the optimization of its hybrid offspring, this study used Huainan pigs as the maternal line and bred them with Yorkshire, Landrace, and Berkshire sires. This approach produced three hybrid combinations: Yorkshire × Huainan (YH), Landrace × Huainan (LH), and Berkshire × Huainan (BH). The body size, fat ratio, and average backfat thickness of these hybrid progeny were evaluated under the same feeding management and nutritional circumstances. The results revealed that the average backfat thickness of YH was significantly lower than that of LH and BH. In order to better understand the causes of these variations, fecal samples were taken from three pigs in each group for metabolomic analysis. A total of 2291 metabolites were identified, including benzene derivatives (16.6%), amino acids and their metabolites (14.5%), and organic acids (13.4%), with pyruvaldehyde and norethindrone acetate elevated in YH compared to LH and BH. In addition, the three hybrid pig groups commonly exhibited differences in the “glycerophospholipid metabolism” pathway. This variation may also contribute to differences in their fat ratio and backfat thickness. Our findings provide a novel perspective on the role of hybrid vigor in advancing the genetic population of Huainan pigs, while also revealing the unique metabolic characteristics of the YH with regard to fat deposition. This study is expected to enhance the conservation and effective utilization of genetic resources within the Huainan pig population.

Whole-genome sequencing reveals genetic structure and adaptive genes in Nepalese buffalo breeds.

Indigenous buffaloes, as the important livestock species contributing to economy of the country, are the lifeline of livelihood in Nepal. They are distributed across diverse geographical regions of the country and have adapted to various feeding, breeding, and management conditions. The larger group of these native buffalo breeds are present in narrow and stiff hilly terrains. Their dispersal indicates a possible environmental adaptation mechanism, which is crucial for the conservation of these breeds. We utilized whole-genome sequencing (WGS) to investigate the genetic diversity, population structure, and selection signatures of Nepalese indigenous buffaloes. We compared 66 whole-genome sequences with 118 publicly available sequences from six river and five swamp buffalo breeds. Genomic diversity parameters indicated genetic variability level in the Nepalese buffaloes comparable to those of Indian breeds, and population genetic structure revealed distinct geography-mediated genetic differentiation among these breeds. We used locus-specific branch length analysis (LSBL) for genome-wide scan, which revealed a list of potentially selected genes in Lime and Parkote breeds that inhabit the hilly region. A gene ontology (GO) analysis discovered that many GO terms were associated with cardiac function regulation. Furthermore, complementary analyses of local selection signatures, tissue expression profiles, and haplotype differences identified candidate genes, including KCNE1, CSF1R, and PDGFRB, related to the regulation of cardiac and pulmonary functions. This study is a comprehensive WGS-based genetic analysis of the native Nepalese buffalo breeds. Our study suggested that the Nepalese "hilly" buffaloes, especially the Lime and Parkote breeds, have undergone some characteristic genetic changes and evolved increased cardiac and pulmonary function for their adaptation to the steep hilly terrains of the country.

Identification and Gene Fine Mapping of the Bisexual Sterility Mutant Meiosis Abnormal Bisexual Sterility 1 in Rice

Exploring the genes regulating rice fertility is of great value for studying the molecular mechanisms of rice reproductive development and production practices. In this study, we identified a sterile mutant from the mutant library induced by ethyl methanesulfonate (EMS), designated as meiosis abnormal bisexual sterility 1 (mabs1). The mabs1 mutant exhibits no phenotypic differences from the wild-type during the vegetative growth phase but shows complete sterility during the reproductive growth phase. Phenotypic observations revealed that both pollen and embryo sac fertility are lost in mabs1. Notably, in mabs1, the development of the anther inner and outer walls, tapetum degeneration, and callose synthesis and degradation all proceed normally, yet meiosis fails to form normal tetrads. Genetic analysis indicated that this mutant trait is controlled by a single recessive nuclear gene. By constructing a genetic segregation population, we successfully mapped the MABS1 gene to a 49 kb region between primer markers Y7 and Y9 on chromosome 1. Resequencing revealed a single-nucleotide substitution in the exon of the LOC_Os01g66170 gene, which resulted in a change from Valine to Isoleucine. Subsequent sequencing of this locus in both wild-type and mabs1 mutants confirmed this mutation. Therefore, we have identified the gene at LOC_Os01g66170 as a candidate for MABS1, a previously unreported novel gene involved in rice meiosis. Through RT-qPCR, we found that the expression levels of multiple meiosis-related genes were significantly changed in the mabs1 mutant. Therefore, we believe that MABS1 is also involved in the process of rice meiosis. This study lays the groundwork for a functional study of MABS1.

Genetic Diversity and Population Structure Analysis in Guar

Guar [Cyamopsis tetragonoloba (L.) Taub] was domesticated in India and Pakistan. It is mainly self-pollinated, bushy, and deeply tap-rooted. Guar seed endosperm contains galactomannan gum, which is used in many food products, pharmaceuticals, cosmetics, explosives, meat products, and pet foods, and in the textile industry, yet its genetic diversity remains largely underexplored. Using 7000 high-quality single nucleotide polymorphism (SNP) markers acquired from genotyping by sequencing (GBS), we analyzed the genetic diversity and population structure in 225 guar accessions from India, Pakistan, and the United States. Structure Harvester revealed that K = 3 had the best delta K, whereas K = 2 had the second-highest delta K. Three major genetic clusters (K = 3) were identified using population structure analysis, utilizing an admixture model: 156 accessions (69.3%) were classified into Q1, 23 accessions (10.2%) in Q2, and 16 accessions (7.1%) in Q3. The remaining 30 accessions (13.3%) were included in the admixture. In all three of the subpopulations at K = 3, most of the guar accessions came from India. We also found that these clusters mostly correlated with geographic origins. Results showed that the Q2 and Q3 subpopulations included only guar accessions from India. Genetic resources from Q2 and Q3 may represent an untapped reservoir for introducing beneficial variety into the U.S. breeding populations. This genetic diversity and population structure analysis of the guar gene pool will be of interest to conduct allele mining and donor parent selection for the development of new and better guar germplasm for desired traits.

What are microsatellites and how to choose the best tool: a user-friendly review of SSR and 74 SSR mining tools

Microsatellites, also known as SSR or STR, are essential molecular markers in genomic research, playing crucial roles in genetic mapping, population genetics, and evolutionary studies. Their applications range from plant breeding to forensics, highlighting their diverse utility across disciplines. Despite their widespread use, traditional methods for SSR analysis are often laborious and time-consuming, requiring significant resources and expertise. To address these challenges, a variety of computational tools for SSR analysis have been developed, offering faster and more efficient alternatives to traditional methods. However, selecting the most appropriate tool can be daunting due to rapid technological advancements and the sheer number of options available. This study presents a comprehensive review and analysis of 74 SSR tools, aiming to provide researchers with a valuable resource for SSR analysis tool selection. The methodology employed includes thorough literature reviews, detailed tool comparisons, and in-depth analyses of tool functionality. By compiling and analyzing these tools, this study not only advances the field of genomic research but also contributes to the broader scientific community by facilitating informed decision-making in the selection of SSR analysis tools. Researchers seeking to understand SSRs and select the most appropriate tools for their projects will benefit from this comprehensive guide. Overall, this study enhances our understanding of SSR analysis tools, paving the way for more efficient and effective SSR research in various fields of study.

Genetic Variation and Population Structure of Clonorchis sinensis: An In Silico Analysis

Clonorchis sinensis is an important zoonotic parasite that is mainly prevalent in China, Korea, Vietnam and the Russian Far East. To explore the genetic variation and population structure of C. sinensis, an in silico analysis was conducted based on mitochondrial cytochrome c oxidase subunit 1 (COX1), ribosomal internal transcribed spacer 1 (ITS1) and ribosomal internal transcribed spacer 2 (ITS2) sequences. The sequences obtained from NCBI were truncated for further analyses, including haplotype network, phylogenetic, gene flow, diversity and neutrality analyses. The results showed that there were 20, 11 and 4 haplotypes for COX1, ITS1 and ITS2, respectively. The results of both the haplotype network and phylogenetic analyses indicated that the haplotypes for each type of sequence from the same country were not all clustered together. Haplotype diversity values were all lower than 0.5. Values of nucleotide diversity were higher than 0.005, except for ITS2. Tajima’s D and Fu’s Fs values were all negative, and p-values showed significant differences, indicating that the population of C. sinensis is growing. Fst values were all lower than 0.05. In conclusion, this study found that there are specific variations of C. sinensis in different countries, and the population of this parasite is growing with less genetic variation. The findings provide a crucial foundation for understanding the molecular epidemiology and population dynamics of C. sinensis.

Two fitness inference schemes compared using allele frequencies from 1,068,391 sequences sampled in the UK during the COVID-19 pandemic.

Throughout the course of the SARS-CoV-2 pandemic, genetic variation has contributed to the spread and persistence of the virus. For example, various mutations have allowed SARS-CoV-2 to escape antibody neutralization or to bind more strongly to the receptors that it uses to enter human cells. Here, we compared two methods that estimate the fitness effects of viral mutations using the abundant sequence data gathered over the course of the pandemic. Both approaches are grounded in population genetics theory but with different assumptions. One approach, tQLE, features an epistatic fitness landscape and assumes that alleles are nearly in linkage equilibrium. Another approach, MPL, assumes a simple, additive fitness landscape, but allows for any level of correlation between alleles. We characterized differences in the distributions of fitness values inferred by each approach and in the ranks of fitness values that they assign to sequences across time. We find that in a large fraction of weeks the two methods are in good agreement as to their top-ranked sequences, \textit{i.e.} as to which sequences observed that week are most fit. We also find that agreement between the ranking of sequences varies with genetic unimodality in the population in a given week.&#xD.

Multisample motif discovery and visualization for tandem repeats.

Tandem Repeats (TR) occupy a significant portion of the human genome and are the source of polymorphism due to variations in sizes and motif compositions. Some of these variations have been associated with various neuropathological disorders, highlighting the clinical importance of assessing the motif structure of TRs. Moreover, assessing the TR motif variation can offer valuable insights into evolutionary dynamics and population structure. Previously, characterizations of TRs have been limited by short-read sequencing technology, which lacks the ability to accurately capture the full TR sequences. As long-read sequencing becomes more accessible and can capture the full complexity of TRs, there is now also a need for tools to characterize and analyze TRs using long-read data across multiple samples. In this study, we present MotifScope, a novel algorithm for characterization and visualization of TRs based on a de novo k-mer approach for motif discovery. Comparative analysis against established tools reveals that MotifScope can identify a greater number of motifs and more accurately represent the underlying repeat sequence. Moreover, MotifScope has been specifically designed to enable motif composition comparisons across assemblies of different individuals, as well as across long-read sequencing reads within an individual, through combined motif discovery and sequence alignment. We showcase potential applications of MotifScope in diverse fields, including population genetics, clinical settings, and forensic analyses.

The evolution of dbSNP: 25 years of impact in genomic research.

The Single Nucleotide Polymorphism Database (dbSNP), established in 1998 by the National Center for Biotechnology Information (NCBI), has been a critical resource in genomics for cataloging small genetic variations. Originally focused on single nucleotide polymorphisms (SNPs), dbSNP has since expanded to include a variety of genetic variants, playing a key role in genome-wide association studies (GWAS), population genetics, pharmacogenomics, and cancer research. Over 25 years, dbSNP has grown to include more than 4.4 billion submitted SNPs and 1.1 billion unique reference SNPs, providing essential data for identifying disease-related genetic variants and studying human diversity. Integrating large-scale projects like 1000 Genomes, gnomAD, TOPMed, and ALFA has expanded dbSNP's catalog of human genetic variation, increasing its usefulness for research and clinical applications. Keeping up with advancements such as next-generation sequencing and cloud-based infrastructure, dbSNP remains a cornerstone of genetic research supporting continued discoveries in precision medicine and population genomics. DATABASE URL: https://www.ncbi.nlm.nih.gov/snp.

Genetic diversity and population structure of local avocado (Persea americana mill.) from northern Tanzania assessed using SSR markers

Genetic diversity and population structure of local avocado (Persea americana mill.) from northern Tanzania assessed using SSR markers

Genome-wide assessment of genetic variation and population structure in cultivated vanilla from Madagascar

Genome-wide assessment of genetic variation and population structure in cultivated vanilla from Madagascar

Inference of Genetic Diversity, Population Structure, and Selection Signatures in Xiangxi White Buffalo of China Through Whole-Genome Resequencing

(1) Background: Buffaloes are crucial livestock species for food and service in tropical and subtropical regions. Buffalo genetics, particularly in indigenous Chinese breeds such as the Xiangxi white buffalo (XWB), remains an intriguing area of study due to its unique traits and regional significance. (2) Methods: This investigation utilized the whole-genome sequences of twenty XWBs (newly sequenced), along with eighty published whole-genome sequences of other buffalo breeds (including Guizhou white buffalo, river buffalo, and Chinese buffalo in the Yangtze River). Using whole-genome sequencing analysis technology, the population structure, genomic diversity, and selection signatures of XWB were determined. (3) Results: This study revealed that the XWB, being phylogenetically positioned in the middle and lower reaches of the Yangtze River, exhibited substantial genomic diversity. Employing four selection sweep detection methods (CLR, iHS, π-ratio, and FST), several genes were positively identified for adaptive traits in the XWB, including coat color phenotypes (ASIP, KIT), the nervous system (GRIK2), reproduction (KCNIP4), growth and development (IFNAR1, BMP6, HDAC9, MGAT4C, and SLC30A9), the body (LINGO2, LYN, and FLI1), immunity (IRAK3 and MZB1), and lactation (TP63, LPIN1, SAE1). (4) Conclusions: In conclusion, this study enhances our understanding of the genetic distinctiveness and adaptive traits of XWB, highlighting selection signatures crucial for future breeding and conservation and ensuring sustainable use of this vital livestock resource.