Effective Population Research Articles

FnR: R package for computing inbreeding and numerator relationship coefficients

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

Population genetic structure of the European pond turtle (Emys orbicularis) at the northern edge of the Pannonian Basin

Populations of plants and animals at the edge of a species’ range may exhibit lower genetic diversity due to a combination of historical (e.g., postglacial dispersal) and contemporary (e.g., low effective population size, strong genetic drift or bottlenecks) processes. In this study, we analysed the genetic structure of populations of the European pond turtle (Emys orbicularis) in the northern part of the Pannonian Basin (Slovakia, Hungary, Austria), which represents the edge of the species’ range. Using mitochondrial DNA (mtDNA) and microsatellite markers, we aimed to assign individuals and populations into phylogeographic lineages and to determine genetic variation and population structure. Our study revealed that most individuals from natural populations belong to mtDNA lineage II (widely distributed haplotype IIa, and likely endemic haplotypes IIl and IIo), which is typical for the whole Pannonian region. However, non-native haplotypes were detected in one Slovak population (haplotype Ib) and in the Austrian Donau Auen National Park (haplotype IVa). Microsatellite markers revealed high variability comparable to that in other parts of the range, suggesting that marginal populations may not necessarily have reduced genetic diversity. A genetically mixed population was found in the Donau Auen National Park, consisting of turtles belonging to mtDNA lineage II (taxonomically assessed as subspecies E. orbicularis orbicularis) and lineage IV (taxonomically assessed as E. orbicularis hellenica), which is native to the circum-Adriatic region. The admixture suggests intense hybridization and introgression between the two lineages (subspecies).

Microsatellite-based bottleneck analysis and migration events among four native Turkish goat breeds

Abstract. Molecular data analyzed by accurate statistical approaches not only have the potential to investigate demographic events faced in the past, e.g., migration, but they also offer significant hints such as bottleneck effects to take suitable measures for sustainable breeding in farm animals. In this study, the genetic bottleneck and migration events among four native Turkish goat populations, i.e., Hair, Honamlı, Kabakulak, and Norduz, were assessed using 367 multi-allelic data belonging to 28 microsatellite loci. The null hypothesis was not rejected for the Wilcoxon sign rank test in the infinite allele model, two-phase model, and stepwise mutation model, while a normal L-shaped distribution of allele frequencies was observed in terms of mode-shift indicators in four native Turkish goat populations. Both the Wilcoxon sign rank test and the mode-shift indicator demonstrated that Anatolian goat populations have had a lack of recent genetic bottlenecks and have maintained their effective population sizes over the generations. Moreover, the 95 % confidence interval confirmed that the effective population sizes of Hair, Honamlı, and Kabakulak may reach infinity, while the highest effective population size for Norduz was estimated at 794.5, when the lowest allele frequency was considered to be 0.01. Up to four migration events revealed a significant migration from Norduz to Hair and Kabakulak populations. In contrast, no migration from other populations to Norduz was observed, most probably due to its geographic isolation. The bottleneck results may serve as a guide for future management practices, whereas further studies, especially on a whole-genome basis, are needed to confirm migration events among Anatolian goat breeds.

Insights into the genomic homogeneity of Moroccan indigenous sheep breeds though the lens of runs of homozygosity

Numerous studies have indicated that Morocco’s indigenous sheep breeds are genetically homogenous, posing a risk to their survival in the challenging harsh climate conditions where they predominantly inhabit. To understand the genetic behind genetic homogeneity through the lens of runs of homozygosity (ROH), we analyzed the whole genome sequences of five indigenous sheep breeds (Beni Guil, Ouled Djellal, D’man, Sardi, Timahdite and Admixed).The results from principal component, admixture, Fst, and neighbour joining tree analyses consistently showed a homogenous genetic structure. This structure was characterized by an average length of 1.83 Mb for runs of homozygosity (ROH) segments, with a limited number of long ROH segments (24–48 Mb and > 48 Mb). The most common ROH segments were those ranging from 1–6 Mb. The most significant regions of homozygosity (ROH Islands) were mostly observed in two chromosomes, namely Chr1 and Chr5. Specifically, ROH Islands were exclusively discovered in the Ouled Djellal breed on Chr1, whereas Chr5 exhibited ROH Islands in all breeds. The analysis of ROH Island and iHS technique was employed to detect signatures of selection on Chr1 and Chr5. The results indicate that Chr5 had a high level of homogeneity, with the same genes being discovered across all breeds. In contrast, Chr1 displays some genetic variances between breeds. Genes identified on Chr5 included SLC39A1, IL23A, CAST, IL5, IL13, and IL4 which are responsible for immune response while genes identified on Chr1 include SOD1, SLAMF9, RTP4, CLDN1, and PRKAA2. ROH segment profile and effective population sizes patterns suggests that the genetic uniformity of studied breeds is the outcome of events that transpired between 250 and 300 generations ago. This research not only contributes to the understanding of ROH distribution across breeds but helps design and implement native sheep breeding and conservation strategies in Morocco. Future research, incorporating a broader sample size and utilizing the pangenome for reference, is recommended to further elucidate these breeds’ genomic landscapes and adaptive mechanisms.

Effective population size in field pea

BackgroundEffective population size (Ne) is a pivotal parameter in population genetics as it can provide information on the rate of inbreeding and the contemporary status of genetic diversity in breeding populations. The population with smaller Ne can lead to faster inbreeding, with little potential for genetic gain making selections ineffective. The importance of Ne has become increasingly recognized in plant breeding, which can help breeders monitor and enhance the genetic variability or redesign their selection protocols. Here, we present the first Ne estimates based on linkage disequilibrium (LD) in the pea genome.ResultsWe calculated and compared Ne using SNP markers from North Dakota State University (NDSU) modern breeding lines and United States Department of Agriculture (USDA) diversity panel. The extent of LD was highly variable not only between populations but also among different regions and chromosomes of the genome. Overall, NDSU had a higher and longer-range LD than the USDA that could extend up to 500 Kb, with a genome-wide average r2 of 0.57 (vs 0.34), likely due to its lower recombination rates and the selection background. The estimated Ne for the USDA was nearly three-fold higher (Ne = 174) than NDSU (Ne = 64), which can be confounded by a high degree of population structure due to the selfing nature of pea.ConclusionsOur results provided insights into the genetic diversity of the germplasm studied, which can guide plant breeders to actively monitor Ne in successive cycles of breeding to sustain viability of the breeding efforts in the long term.

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

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

Pedigree analysis of Korean native chickens: unraveling inbreeding and genetic diversity

This study assessed the trends in inbreeding, effective population size, and genetic diversity across six Korean native chicken lines using pedigree records from 54,383 chickens. Understanding these genetic parameters is significantly important for maintaining healthy and viable chicken populations. The primary objective was to analyze the pedigree data to assess the levels of inbreeding and genetic diversity and to evaluate the effective population size across the different lines. Pedigree analysis revealed that pedigree completeness peaked in the first generation and declined in subsequent generations for all lines. Line A exhibited a mean inbreeding coefficient of 0.0201, whereas the other lines displayed lower mean values ranging from 0.0009 to 0.0098, indicating that inbreeding levels were within an acceptable range and considered safe from extinction. Average relatedness consistently increased with time. Individual increases in inbreeding were the highest in Line A (0.62%), with smaller increases in the other lines ranging from 0.02 to 0.23%. Effective population sizes varied from 81 to 2500, with average coancestry within parental populations ranging from 0.0032 to 0.0290. The fe/fa ratio between 1.00 and 1.69 in the 6 lines suggested a moderate impact during bottleneck events, with subsequent populations recovering well. The genetic diversity loss due to genetic drift and unequal founder contributions ranged from 0.66-3.15%, indicating that considerable genetic variability remains within the populations. The results of this study have practical applications in the management and conservation of genetic resources in poultry breeding programs. By highlighting the importance of monitoring inbreeding and maintaining genetic diversity, the findings can help develop strategies to ensure the long-term sustainability of these chicken lines. This study provides valuable insights into the genetic management of Korean native chicken lines, emphasizing the need for strategic breeding practices to preserve genetic health and diversity.

Phylogeography and genetic structure of Papaver bracteatum populations in Iran based on genotyping-by-sequencing (GBS)

Papaver bracteatum, known for its high thebaine content and absence of morphine, has emerged as a promising alternative to opium poppy for codeine production. In this study, our objective was to create a diverse panel representing the natural variation of this species in Iran. To achieve this, we employed genotyping-by-sequencing to obtain genome-wide distributed single-nucleotide polymorphisms (SNPs) for phylogeographic analysis, population structure assessment, and evaluation of genetic diversity within P. bracteatum populations. A total of 244 P. bracteatum individuals from 13 distinct populations formed seven genetic groups, along with one highly admixed population. We observed a clear split between the populations inhabiting the Alborz Mts. in the east and Zagros Mts. in the west. In between these mountain ranges, the population of Kachal Mangan exhibited a high degree of genetic admixture between both genetic groups. At or after the end of the last glacial maximum, when climate conditions rapidly changed, all P. bracteatum populations experienced a strong demographic bottleneck reducing the already small effective population sizes further before they increased to their recent strengths. Our results suggest that the ongoing climate change together with human pressure on the species’ habitats and limited seed-dispersal ability are potential factors contributing today to rising genetic isolation of P. bracteatum populations. Our results provide genetic data that can be used for conservation measures to safeguard the species’ genetic diversity as a resource for future breeding approaches in this medicinally important species.

Population genomics and demographic modeling enhance our understanding of trophic level interactions in biological control

In this paper, we focus on the application of population genomics, including next-generation sequencing and demographic modeling, to enhance the understanding of historical and contemporary trophic level interactions that are a fundamental basis of biological control. We outline how the integration of ecological approaches with population genomics and demographic modeling can provide additional insights into existing trophic interactions and how natural enemies, their parasitoids, and pest species evolve through these interactions, potentially affecting future levels of biological control. We briefly review and summarize our recent population genomic studies of predatory ladybird beetles and discuss how our results provide guidelines for the use of population genomics techniques and population demographic modeling to enhance our understanding of trophic level interactions. Our intent is to demonstrate as a proof of concept, how these methods can supplement currently used techniques and add a complementary genomic aspect to ecological understandings of trophic level interactions, including estimates of effective population sizes of natural enemies and pest species, examination of correlations between these effective population sizes, and evolutionary changes in natural enemies and pest species which may influence levels of biological control. This in combination with experiments and simulations provide a well-rounded framework for using population genomics to inform multi-tier trophic level interactions.

Benchmarking and Optimization of Methods for the Detection of Identity-By-Descent in High-Recombining Plasmodium falciparum Genomes.

Genomic surveillance is crucial for identifying at-risk populations for targeted malaria control and elimination. Identity-by-descent (IBD) is increasingly being used in Plasmodium population genomics to estimate genetic relatedness, effective population size (N e ), population structure, and signals of positive selection. Despite its potential, a thorough evaluation of IBD segment detection tools for species with high recombination rates, such as P. falciparum, remains absent. Here, we perform comprehensive benchmarking of IBD callers - probabilistic (hmmIBD, isoRelate), identity-by-state-based (hap-IBD, phased IBD) and others (Refined IBD) - using population genetic simulations tailored for high recombination, and IBD quality metrics at both the IBD segment level and the IBD-based downstream inference level. Our results demonstrate that low marker density per genetic unit, related to high recombination relative to mutation, significantly compromises the accuracy of detected IBD segments. In genomes with high recombination rates resembling P. falciparum, most IBD callers exhibit high false negative rates for shorter IBD segments, which can be partially mitigated through optimization of IBD caller parameters, especially those related to marker density. Notably, IBD detected with optimized parameters allows for more accurate capture of selection signals and population structure; IBD-based N e inference is very sensitive to IBD detection errors, with IBD called from hmmIBD uniquely providing less biased estimates of N e in this context. Validation with empirical data from the MalariaGEN Pf 7 database, representing different transmission settings, corroborates these findings. We conclude that context-specific evaluation and parameter optimization are essential for accurate IBD detection in high-recombining species and recommend hmmIBD for quality-sensitive analysis, such as estimation of N e in these species. Our optimization and high-level benchmarking methods not only improve IBD segment detection in high-recombining genomes but also enhance overall genomic analysis, paving the way for more accurate genomic surveillance and targeted intervention strategies for malaria.

Evaluation of heritability partitioning approaches in livestock populations

BackgroundHeritability partitioning approaches estimate the contribution of different functional classes, such as coding or regulatory variants, to the genetic variance. This information allows a better understanding of the genetic architecture of complex traits, including complex diseases, but can also help improve the accuracy of genomic selection in livestock species. However, methods have mainly been tested on human genomic data, whereas livestock populations have specific characteristics, such as high levels of relatedness, small effective population size or long-range levels of linkage disequilibrium.ResultsHere, we used data from 14,762 cows, imputed at the whole-genome sequence level for 11,537,240 variants, to simulate traits in a typical livestock population and evaluate the accuracy of two state-of-the-art heritability partitioning methods, GREML and a Bayesian mixture model. In simulations where a single functional class had increased contribution to heritability, we observed that the estimators were unbiased but had low precision. When causal variants were enriched in variants with low (< 0.05) or high (> 0.20) minor allele frequency or low (below 1st quartile) or high (above 3rd quartile) linkage disequilibrium scores, it was necessary to partition the genetic variance into multiple classes defined on the basis of allele frequencies or LD scores to obtain unbiased results. When multiple functional classes had variable contributions to heritability, estimators showed higher levels of variation and confounding between certain categories was observed. In addition, estimators from small categories were particularly imprecise. However, the estimates and their ranking were still informative about the contribution of the classes. We also demonstrated that using methods that estimate the contribution of a single category at a time, a commonly used approach, results in an overestimation. Finally, we applied the methods to phenotypes for muscular development and height and estimated that, on average, variants in open chromatin regions had a higher contribution to the genetic variance (> 45%), while variants in coding regions had the strongest individual effects (> 25-fold enrichment on average). Conversely, variants in intergenic or intronic regions showed lower levels of enrichment (0.2 and 0.6-fold on average, respectively).ConclusionsHeritability partitioning approaches should be used cautiously in livestock populations, in particular for small categories. Two-component approaches that fit only one functional category at a time lead to biased estimators and should not be used.

Epidemiology of carbapenem-resistant Klebsiella pneumoniae in China and the evolving trends of predominant clone ST11: a multicentre, genome-based study.

Carbapenem-resistant Klebsiella pneumoniae (CRKP) is a major nosocomial infectious pathogen with rapidly increasing prevalence. The genomic epidemiological characteristics of CRKP nationwide, especially the evolving trends within the predominant clones, should be evaluated clearly. We collected 3415 K. pneumoniae strains from 28 hospitals across China. Antimicrobial susceptibility testing and WGS were performed. Subsequent genomic analyses, including sequence typing, K-locus (KL) identification, antimicrobial resistance gene screening, and virulence score assessment were performed. The phylogenetic relationship of clonal group 11 was determined based on core-genome analysis, and the presence of the pLVPK-like virulence plasmid in ST11 isolates was confirmed using plasmid core-gene analysis. Additionally, the trends of the ST11 lineage with different KL types on a global scale were investigated using Beast2. Of the K. pneumoniae strains, 708 were identified as CRKP isolates (20.7%), of which 97.7% were MDR. ST11 was the predominant clone, and KPC-2 was the prevalent carbapenemase in China, although the prevalence of specific clones and carbapenemases varied by geographic region. Among ST11 isolates, KL47 and KL64 were the predominant KL types, and KL64 gradually replaced KL47, with a higher percentage of KL64 isolates harbouring the pLVPK-like plasmid. Global genome data showed a significant increase in the effective population size of KL64 over the last 5 years. The prevalence of CRKP was very high in certain regions in China. The increasing convergence of virulence and resistance, particularly in ST11-KL64 isolates, should be given more attention and further investigation.

ONeSAMP 3.0: estimation of effective population size via single nucleotide polymorphism data from one population.

The genetic effective size (Ne) is arguably one of the most important characteristics of a population as it impacts the rate of loss of genetic diversity. Methods that estimate Ne are important in population and conservation genetic studies as they quantify the risk of a population being inbred or lacking genetic diversity. Yet there are very few methods that can estimate the Ne from data from a single population and without extensive information about the genetics of the population, such as a linkage map, or a reference genome of the species of interest. We present ONeSAMP 3.0, an algorithm for estimating Ne from single nucleotide polymorphism data collected from a single population sample using approximate Bayesian computation and local linear regression. We demonstrate the utility of this approach using simulated Wright-Fisher populations, and empirical data from five endangered Channel Island fox (Urocyon littoralis) populations to evaluate the performance of ONeSAMP 3.0 compared to a commonly used Ne estimator. Our results show that ONeSAMP 3.0 is broadly applicable to natural populations and is flexible enough that future versions could easily include summary statistics appropriate for a suite of biological and sampling conditions. ONeSAMP 3.0 is publicly available under the GNU General Public License at https://github.com/AaronHong1024/ONeSAMP_3.

Illuminating Genetic Diversity and Selection Signatures in Matou Goats through Whole-Genome Sequencing Analysis.

(1) Background: Matou goats, native to Hunan and Hubei provinces in China, are renowned for their exceptional meat and skin quality. However, a comprehensive whole-genome-based exploration of the genetic architecture of this breed is scant in the literature. (2) Methods: To address this substantial gap, we used whole-genome sequences of 20 Matou goats and compared them with published genomic data of 133 goats of different breeds across China. This comprehensive investigation sought to assess genetic diversity, population structure, and the presence of genomic selection signals. (3) Results: The whole genome of Matou goat populations yielded a substantial catalog of over 19 million single nucleotide polymorphisms (SNPs), primarily distributed within intergenic and intron regions. The phylogenetic tree analysis revealed distinct clades corresponding to each goat population within the dataset. Notably, this analysis positioned Matou goats in a closer genetic affinity with Guizhou White goats, compared to other recognized goat breeds. This observation was corroborated by principal component analysis (PCA) and admixture analysis. Remarkably, Matou goats exhibited diminished genetic diversity and a notable degree of inbreeding, signifying a reduced effective population size. Moreover, the study employed five selective sweep detection methods (including PI, CLR, PI-Ratio, Fst, and XP-EHH) to screen top signal genes associated with critical biological functions, encompassing cardiomyocytes, immunity, coat color, and meat quality. (4) Conclusions: In conclusion, this study significantly advances our understanding of the current genetic landscape and evolutionary dynamics of Matou goats. These findings underscore the importance of concerted efforts in resource conservation and genetic enhancement for this invaluable breed.

Mitonuclear compatibility is maintained despite relaxed selection on male mitochondrial DNA in bivalves with doubly uniparental inheritance.

Mitonuclear coevolution is common in eukaryotes, but bivalve lineages that have doubly uniparental inheritance (DUI) of mitochondria may be an interesting example. In this system, females transmit mtDNA (F mtDNA) to all offspring, while males transmit a different mtDNA (M mtDNA) solely to their sons. Molecular evolution and functional data suggest oxidative phosphorylation (OXPHOS) genes encoded in M mtDNA evolve under relaxed selection due to their function being limited to sperm only (vs. all other tissues for F mtDNA). This has led to the hypothesis that mitonuclear coevolution is less important for M mtDNA. Here, we use comparative phylogenetics, transcriptomics, and proteomics to understand mitonuclear interactions in DUI bivalves. We found nuclear OXPHOS proteins coevolve and maintain compatibility similarly with both F and M mtDNA OXPHOS proteins. Mitochondrial recombination did not influence mitonuclear compatibility and nuclear-encoded OXPHOS genes were not upregulated in tissues with M mtDNA to offset dysfunction. Our results support that selection maintains mitonuclear compatibility with F and M mtDNA despite relaxed selection on M mtDNA. Strict sperm transmission, lower effective population size, and higher mutation rates may explain the evolution of M mtDNA. Our study highlights that mitonuclear coevolution and compatibility may be broad features of eukaryotes.

Repeated horizontal acquisition of lagriamide-producing symbionts in Lagriinae beetles.

Microbial symbionts associate with multicellular organisms on a continuum from facultative associations to mutual codependency. In some of the oldest intracellular symbioses there is exclusive vertical symbiont transmission, and co-diversification of symbiotic partners over millions of years. Such symbionts often undergo genome reduction due to low effective population sizes, frequent population bottlenecks, and reduced purifying selection. Here, we describe multiple independent acquisition events of closely related defensive symbionts followed by genome erosion in a group of Lagriinae beetles. Previous work in Lagria villosa revealed the dominant genome-eroded symbiont of the genus Burkholderia produces the antifungal compound lagriamide and protects the beetle's eggs and larvae from antagonistic fungi. Here, we use metagenomics to assemble 11 additional genomes of lagriamide-producing symbionts from seven different host species within Lagriinae from five countries, to unravel the evolutionary history of this symbiotic relationship. In each host species, we detected one dominant genome-eroded Burkholderia symbiont encoding the lagriamide biosynthetic gene cluster (BGC). Surprisingly, however, we did not find evidence for host-symbiont co-diversification, or for a monophyly of the lagriamide-producing symbionts. Instead, our analyses support at least four independent acquisition events of lagriamide-encoding symbionts and subsequent genome erosion in each of these lineages. By contrast, a clade of plant-associated relatives retained large genomes but secondarily lost the lagriamide BGC. In conclusion, our results reveal a dynamic evolutionary history with multiple independent symbiont acquisitions characterized by high degree of specificity. They highlight the importance of the specialized metabolite lagriamide for the establishment and maintenance of this defensive symbiosis.

Molecular Assessment of Genetic Diversity and Bottleneck in Hair Goat Reared in Türkiye

Objective: Being the most preferred and geographically distributed in Türkiye, the Hair goat was screened at a molecular level to evaluate genetic diversity and population structure via microsatellite DNA markers. This paper also aimed to investigate the effects of genetic bottleneck to evaluate whether the Hair goat has maintained its effective population size in recent past. Material and Methods: A total of 411 Hair goats were sampled from farms participating in the "Hair Goat Breeding" project, initiated by the General Directorate of Agricultural Research and Policies in Aydın and Denizli provinces. Sampled animals were genotyped with 18 microsatellite loci to assess genetic diversity, population structure, and genetic bottleneck. Results: A total of 341 different alleles were observed across 18 microsatellite loci in which the highest number of alleles (26) and effective alleles (10.18) were detected in INRA005 and HSC loci, respectively. The average observed heterozygosity (0.73) was lower than the expected value (0.83), whereas all loci turned out to be highly informative (PIC&amp;gt;0.50). Factorial correspondence analysis separated animals into two groups, while a genetic admixture was detected between these groups. STRUCTURE analysis, on the other hand, confirmed that 411 animals were derived from three ancestral populations in which the third group is drawn due to admixed individuals. The Wilcoxon test and mode-shift indicator detected a lack of genetic bottleneck indicating that Hair goats reared in Türkiye have maintained their effective population size in recent past. Conclusion: This study validates that used microsatellite markers are highly polymorphic and could be utilized for revealing genetic diversity in different local goat breeds. The findings recovered in this study could be integrated into breeding and conservation programs, while further studies should adopt SNP array technologies and next-generation sequencing platforms to reveal deeper knowledge about the genetic diversity and population structure of Anatolian goat breeds.

Drivers of genomic diversity and phenotypic development in early phases of domestication in Hermetia illucens.

The black soldier fly (BSF), Hermetia illucens, has the ability to efficiently bioremediate organic waste into usable bio-compounds. Understanding the impact of domestication and mass rearing on fitness and production traits is therefore important for sustainable production. This study aimed to assess patterns of genomic diversity and its association to phenotypic development across early generations of mass rearing under two selection strategies: selection for greater larval mass (SEL lines) and no direct artificial selection (NS lines). Genome-wide single nucleotide polymorphism (SNP) data were generated using 2bRAD sequencing, while phenotypic traits relating to production and population fitness were measured. Declining patterns of genomic diversity were observed across three generations of captive breeding, with the lowest diversity recorded for the F3 generation of both selection lines, most likely due to founder effects. The SEL cohort displayed statistically significantly greater larval weight com the NS lines with pronounced genetic and phenotypic directional changes across generations. Furthermore, lower genetic and phenotypic diversity, particularly for fitness traits, were evident for SEL lines, illustrating the trade-off between selecting for mass and the resulting decline in population fitness. SNP-based heritability was significant for growth, but was low or non-significant for fitness traits. Genotype-phenotype correlations were observed for traits, but individual locus effect sizes where small and very few of these loci demonstrated a signature for selection. Pronounced genetic drift, due to small effective population sizes, is likely overshadowing the impacts of selection on genomic diversity and consequently phenotypic development. The results hold particular relevance for genetic management and selective breeding for BSF in future.

Ecological genetics of isolated loach populations indicate compromised adaptive potential

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

Plasmodium vivax populations in the western Greater Mekong Subregion evaluated using a genetic barcode.

An improved understanding of the Plasmodium vivax populations in the Great Mekong Subregion (GMS) is needed to monitor the progress of malaria elimination. This study aimed to use a P. vivax single nucleotide polymorphism (SNP) barcode to evaluate the population dynamics and explore the gene flow among P. vivax parasite populations in the western GMS (China, Myanmar and Thailand). A total of 315 P. vivax patient samples collected in 2011 and 2018 from four regions of the western GMS were genotyped for 42 SNPs using the high-throughput MassARRAY SNP genotyping technology. Population genetic analysis was conducted to estimate the genetic diversity, effective population size, and population structure among the P. vivax populations. Overall, 291 samples were successfully genotyped at 39 SNPs. A significant difference was observed in the proportion of polyclonal infections among the five P. vivax populations (P = 0.0012, Pearson Chi-square test, χ2 = 18.1), with western Myanmar having the highest proportion (96.2%, 50/52) in 2018. Likewise, the average complexity of infection was also highest in western Myanmar (1.31) and lowest in northeast Myanmar (1.01) in 2018. The older samples from western China in 2011 had the highest pairwise nucleotide diversity (π, 0.388 ± 0.046), expected heterozygosity (He, 0.363 ± 0.02), and the largest effective population size. In comparison, in the neighboring northeast Myanmar, the more recent samples in 2018 showed the lowest values (π, 0.224 ± 0.036; He, 0.220 ± 0.026). Furthermore, the 2018 northeast Myanmar parasites showed high and moderate genetic differentiation from other populations with FST values of 0.162-0.252, whereas genetic differentiation among other populations was relatively low (FST ≤ 0.059). Principal component analysis, phylogeny, and STRUCTURE analysis showed that the P. vivax population in northeast Myanmar in 2018 substantially diverged from other populations. Although the 42 SNP barcode is a valuable tool for tracking parasite origins of worldwide parasite populations, a more extended barcode with additional SNPs is needed to distinguish the more related parasite populations in the western GMS.