Relationships Between Wright's F ST and F IS Statistics in a Context of Wahlund Effect

Analysis of Endangered Andalusian Black Cattle (Negra Andaluza) Reveals Genetic Reservoir for Bovine Black Trunk

This comprehensive study on the Andalusian Black cattle breed reveals a substantial population decline, with the average herd size decreasing significantly from 305.54 to 88.28 animals per herd. This decline is primarily attributed to agricultural changes and the introduction of foreign meat-focused breeds. The male-to-female ratio shift is noteworthy, with more cows than bulls, impacting selection intensity for both genders. Inbreeding levels, though relatively low historically (5.94%) and currently (7.23%), raise concerns as 37.08% historically and 48.82% currently of the animals exhibit inbreeding. Positive assortative mating is evident, reflected by the increasing non-random mating coefficient (α). Key ancestors play a crucial role in shaping genetic diversity, with one ancestor significantly influencing the current genetic pool and the top 10 ancestors contributing substantially. Breed maintains a conservation index of 2.75, indicating relatively high genetic diversity. Recent conservation efforts have led to an increase in registered animals. The Cañadas Reales, historical transhumance routes, may have contributed to genetic connections among provinces. Challenges include the historical bottleneck, demographic changes, and potential impacts from reproductive practices. The Andalusian Black breed's conservation necessitates ongoing efforts in genealogical registration, targeted breeding programs, and collaborative initiatives to address the observed demographic shifts and ensure sustainable genetic diversity.

Mitochondrial DNA D-Loop Polymorphisms among the Galla Goats Reveals Multiple Maternal Origins with Implication on the Functional Diversity of the HSP70 Gene

Despite much attention given to the history of goat evolution in Kenya, information on the origin, demographic history, dispersal route, and genetic diversity of Galla goats remains unclear. Here, we examined the genetic background, diversity, demographic history, and population genetic variation of Galla goats using mtDNA D-loop and HSP70 single-nucleotide polymorphism markers. The results revealed 90 segregating sites and 68 haplotypes in a 600-bp mtDNA D-loop sequence. The overall mean mitochondrial haplotype diversity was 0.993. The haplotype diversities ranged between 0.8939 ± 0.0777 and 1.0000 ± 0.0221 in all populations supporting high genetic diversity. Mitochondrial phylogenetic analysis revealed three Galla goat haplogroups (A, G, and D), supporting multiple maternal ancestries, of which haplogroup A was the most predominant. Analysis of molecular variance (AMOVA) showed considerable variation within populations at 94.39%, evidence of high genetic diversity. Bimodal mismatch distribution patterns were observed while most populations recorded negative results for Tajima and Fu's Fs neutrality tests supporting population expansion. Genetic variation among populations was also confirmed using HSP70 gene fragment sequences, where six polymorphic sites which defined 21 haplotypes were discovered. Analysis of molecular variance revealed a significant FST index value of 0.134 and a high FIS index value of 0.746, an indication of inbreeding. This information will pave the way for conservation strategies and informed breeding to improve Galla or other goat breeds for climate-smart agriculture.

Commonly used Hardy-Weinberg equilibrium filtering schemes impact population structure inferences using RADseq data

Reduced representation sequencing (RRS) is a widely used method to assay the diversity of genetic loci across the genome of an organism. The dominant class of RRS approaches assay loci associated with restriction sites within the genome (restriction site associated DNA sequencing, or RADseq). RADseq is frequently applied to non‐model organisms since it enables population genetic studies without relying on well‐characterized reference genomes. However, RADseq requires the use of many bioinformatic filters to ensure the quality of genotyping calls. These filters can have direct impacts on population genetic inference, and therefore require careful consideration. One widely used filtering approach is the removal of loci that do not conform to expectations of Hardy–Weinberg equilibrium (HWE). Despite being widely used, we show that this filtering approach is rarely described in sufficient detail to enable replication. Furthermore, through analyses of in silico and empirical data sets we show that some of the most widely used HWE filtering approaches dramatically impact inference of population structure. In particular, the removal of loci exhibiting departures from HWE after pooling across samples significantly reduces the degree of inferred population structure within a data set (despite this approach being widely used). Based on these results, we provide recommendations for best practice regarding the implementation of HWE filtering for RADseq data sets.

High-throughput genotype-based population structure analysis of selected buffalo breeds

India is considered as the home tract of some of the best buffalo breeds. However, the genetic structure of the Indian river buffalo is poorly understood. Hence, there is a need to characterize the populations and understand the genetic structure of various buffalo breeds for selection and to design breeding strategies. In this study, we have analyzed genetic variability and population structure of seven buffalo breeds from their respective geographical regions using Axiom Buffalo Genotyping Array. Diversity, as measured by expected heterozygosity, ranged from 0.364 in Surti to 0.384 in Murrah breed, and pair-wise F_ST values revealed the lowest genetic distance between Murrah and Nili-Ravi (0.0022), while the highest between Surti and Pandharpuri (0.030). Principal component analysis and structure analysis unveiled the differentiation of Surti, Pandharpuri, and Jaffarabadi in first two principal components and at K = 4, respectively, while remaining breeds were grouped together as a separate single cluster and admixed. Murrah and Mehsana showed early linkage disequilibrium (LD) decay, while Surti breed showed late decay. In LD blocks to quantitative trait locis (QTLs) concordance analysis, 4.65% of concordance was observed with 873 LD blocks overlapped with 2,330 QTLs. Overall, total 4,090 markers were identified from all LD blocks for six types of traits. Results of this study indicated that these single-nucleotide polymorphism (SNP) markers could differentiate phenotypically distinct breeds like Surti, Pandharpuri, and Jaffarabadi but not others. So, there is a need to develop SNP chip based on SNP markers identified by sequence information of local breeds.

Clustering patterns mirror the geographical distribution and genetic history of Lemnos and Lesvos sheep populations

Elucidating the genetic variation and structure of Lemnos and Lesvos sheep is critical for maintaining local genetic diversity, ecosystem integrity and resilience of local food production of the two North Aegean islands. In the present study, we explored genetic diversity and differentiation as well as population structure of the Lemnos and Lesvos sheep. Furthermore, we sought to identify a small panel of markers with the highest discriminatory power to assign animals across islands. A total number of n = 424 (n = 307, Lemnos and n = 117, Lesvos) ewes, sampled from n = 24 herds dispersed at different geographic regions on the two islands, were genotyped with the 50K SNP array. Mean observed heterozygosity was higher (but not statistically significantly different) in Lesvos than in Lemnos population (0.384 vs. 0.377) while inbreeding levels were higher in Lemnos than Lesvos herds (0.065 vs. 0.031). Results of principal components along with that of admixture analysis and estimated genetic distances revealed genetic clusters corresponding to Lesvos and Lemnos origin and the existence of infrastructure within islands that were associated with geographical isolation and genetic history of the studied populations. In particular, genetic analyses highlighted three geographically isolated herds in Lemnos that are located at mountainous areas of the island and are characterized as representatives of the local sheep by historic data and reports. Admixture analysis also showed a shared genetic background between Lemnos and Lesvos sheep attributable to past gene flow. Little overall genetic differentiation was detected between the two island sheep populations, while 150 discriminatory SNPs could accurately assign animals to their origin. Present results are comparable with those reported in the worldwide sheep breeds, suggesting geography related genetic patterns across and within islands and the existence of the local Lemnos sheep.

The genomic ancestry of Jat Sikh population from Northwest India inferred from 15 autosomal STR markers using capillary electrophoresis

BACKGROUND

Autosomal STR typing using capillary electrophoresis is a reliable method for establishing parentage and for deciphering genomic ancestry.

AIM

This study was planned to show the genetic diversity of the Jat Sikh population, which is a widespread community of the Punjab region, and to assess its genetic relationship with existing Indian populations.

SUBJECTS AND METHODS

Blood samples of unrelated healthy individuals of the Jat Sikhs (n = 123) were used in this study. Fifteen autosomal STR markers along with the sex determination genetic marker Amelogenin were amplified using AmpFlSTR^®Identifiler^® Plus kit, and genetic analyser 3100 was used for genotyping.

RESULTS

A total of 246 alleles were observed with allele frequencies ranging from 0.004 to 0.447. The heterozygosity ranged from 0.659 to 0.886, and all studied loci were in Hardy-Weinberg Equilibrium (HWE). Fibrinogen A alpha (Aα) chain (FGA) was found to be the most polymorphic and also the most discriminating locus in the studied population. Neighbor-joining (NJ) tree, principal component analysis (PCA) plot, and Nei's Distance matrix revealed genetic affinity with the previously reported Jatt Sikh (Punjab) population and showed the outlier nature of this population compared with other Indian populations.

CONCLUSION

The data generated by this study enhance the database of Indian populations to be used in civil and forensic cases and also in other population-based genetic studies.

Assessing the genetic diversity in Argopecten nucleus (Bivalvia: Pectinidae), a functional hermaphrodite species with extremely low population density and self-fertilization: Effect of null alleles

Argopecten nucleus is a functional hermaphroditic pectinid species that exhibits self-fertilization, whose natural populations have usually very low densities. In the present study, the genetic diversity of a wild population from Neguanje Bay, Santa Marta (Colombia), was estimated using microsatellite markers, and the effect of the presence of null alleles on this estimation was assessed. A total of 8 microsatellite markers were developed, the first described for this species, and their amplification conditions were standardized. They were used to determine the genotype of 48 wild individuals from Naguanje Bay, and 1,010 individuals derived from the offspring of 38 directed crosses. For each locus, the frequencies of the identified alleles, including null alleles, were estimated using the statistical package Micro-Checker, and the parental genotypes were confirmed using segregation analysis. Three to 8 alleles per locus with frequencies from 0.001 to 0.632 were detected. The frequencies of null alleles ranged from 0.10 to 0.45, with Ho from 0.0 to 0.79, and He from 0.53 to 0.80. All loci were in H-W disequilibrium. The null allele frequencies values were high, with lower estimations using segregation analysis than estimated using Micro-Checker. The present results show high levels of population genetic diversity and indicate that null alleles were not the only cause of deviation from H-W equilibrium in all loci, suggesting that the wild population under study presents signs of inbreeding and Wahlund effect.

Null Alleles and FIS × FST Correlations

Three published papers in this journal have considered the proposition that, under a Wahlund effect caused by population mixture, a positive correlation is expected between single-locus values of FIS for a sample from the mixture and FST between the populations contributing to the mixture. Two of the papers assumed unbiased samples to estimate FST but did not consider possible effects of null alleles; the other paper focused on effects of nulls but used biased samples that also included Wahlund effects to estimate FST. The result is an information gap regarding scenarios that include null alleles but have unbiased estimates of FST. Simulations were used to fill this information gap, with the following results: 1) converting ~10% of alleles to nulls substantially reduced apparent heterozygosity and substantially increased FIS, with few exceptions; 2) adding null alleles also increased FST at most loci, although the effect was much more modest; 3) null alleles generally degraded correlations between FIS and FST, but the relationship remained relatively strong for FST ≥ 0.06; and 4) null alleles had only a small effect on correlations between r2, a measure of linkage disequilibrium between pairs of loci, and the product of FST values for those loci. These results argue for some caution in interpreting FIS × FST correlations under conditions where null alleles might be common and suggest that two-locus analyses might provide more robust assessments of Wahlund effects.

Revisiting FIS, FST, Wahlund Effects, and Null Alleles

Null alleles and Wahlund effects are well known causes of heterozygote deficits in empirical population genetics studies as compared to Hardy-Weinberg genotypic expectations. Some authors have theoretically studied the relationship of Wright's FIS computed from subsamples displaying a Wahlund effect and FST before the Wahlund effect, as can occasionally be obtained from populations of long-lived organisms. In the 2 subsample case, a positive relationship between these 2 parameters across loci would represent a signature of Wahlund effects. Nevertheless, for most organisms, getting 2 independent subsamples of the same cohort and population, one with a Wahlund effect and the other without, is almost never achieved and most of the time, empirical population geneticists only collect a single sample, with or without a Wahlund effect, or with or without null alleles. Another issue is that null allele increase FIS and FST altogether and thus may also create such correlation. In this article, I show that, for organisms collected in a single sample, which corresponds to the most common situation, Wahlund effects and null alleles affect the values of both FIS and FST though in the opposite direction. I also show that Wahlund effect produces no or weak positive correlation between the 2 F-statistics, while null alleles generate a strong positive correlation between them. Variation of these F-statistics is small and even minimized for FST under Wahlund effects as compared to null alleles. I finally propose a determination key to interpret data with heterozygote deficits.

Genetic analysis of African lions (Panthera leo) in Zambia support movement across anthropogenic and geographical barriers

The Luangwa Valley in eastern Zambia is a transverse offshoot of the Great Rift Valley system. This region appears to have an isolating effect as evidenced by suspected endemic subspecies, such as the Cookson's wildebeest and Thornicroft's giraffe. Recent mitochondrial DNA studies demonstrated that African lions in Zambia consist of two highly diverse eastern and western sub-populations. Herein, we report nuclear and mitochondrial DNA results from 409 lions that support this population substructure across Zambia but proposes only partial isolation of the Luangwa Valley with more movement between the populations than previously thought. Population assignment analysis identifies two populations with little evidence of admixture assigning lions to either the eastern or western sub-populations. A high occurrence of private alleles and clear evidence for a Wahlund effect further justify the presence of a highly structured population. But, while mitochondrial DNA analysis still shows little to no matrilineal gene flow (FST = 0.53) between sub-populations, microsatellite analysis suggests there is gene flow (FST = 0.04) with low but significant isolation-by-distance and an average of 6 migrants per generation. Evidence of isolation-by-distance is also found in factorial correspondence analysis with the Lower Zambezi National Park and eastern corridor clusters overlapping isolated clusters of the Luangwa Valley and western sub-population. From this evidence, the Luangwa Valley appears separated from the western sub-population with some dispersal through the southern regions of the eastern sub-population. Both the eastern and western sub-populations have high heterozygosity (0.68 and 0.69, respectively) and genetic diversity (0.47 and 0.50, respectively) values, indicative of genetically healthy populations.

Isolation by Distance, Source-Sink Population Dynamics and Dispersal Facilitation by Trade Routes: Impact on Population Genetic Structure of a Stored Grain Pest

Population genetic structure of agricultural pests can be impacted not only by geographic distance and the broader ecological and physical barriers but also by patterns related to where crops are produced and how they are moved after harvest. Stored-product pests, for instance, specialize in exploiting grains such as wheat and rice from on-farm storage through transportation to final processing at often geographically distant locations; therefore human-aided movement may impact their dispersal. Although stored product insects are associated with stored grain, they can also exploit resources in the surrounding environments so different ecological regions where the grain is grown and stored may also influence population structure. Here we used 1,156 SNP markers to investigate how geographic distance, ecological and agricultural variables can impact the genetic structure and gene flow of the stored food pest beetle Rhyzopertha dominica We found a substantial degree of admixture between weakly structured populations in the US. Ecological regions were more important in explaining R. dominica population structure than crop type, suggesting insect movement between wheat and rice grain distribution channels. We have also found a significant correlation between the genetic and geographical distance (i.e., isolation by distance). However, our modeling approach combining the ecological and management variables has highlighted the importance of the volume of grain received by a location in the dispersal dynamics of the pest. The first-generation migrant analysis offered additional supported to movement over great distances that are likely associated with grain movement. Our data suggest that a multitude of factors play small but significant parts in the movement dynamics of the pest. The beetles can take advantage of the source-sink dynamic of grain movement in the US, but also engage in a high rate of movement at the local scale. Understanding population structure for R. dominica will provide insights into the potential for local processes of adaptation and broader patterns of movement that will impact management programs and the potential for spread of resistance genes.

Genetic Differentiation of Glossina pallidipes Tsetse Flies in Southern Kenya

The tsetse fly Glossina pallidipes, the major vector of the parasite that causes animal African trypanosomiasis in Kenya, has been subject to intense control measures with only limited success. The G. pallidipes population dynamics and dispersal patterns that underlie limited success in vector control campaigns remain unresolved, and knowledge on genetic connectivity can provide insights, and thereby improve control and monitoring efforts. We therefore investigated the population structure and estimated migration and demographic parameters in G. pallidipes using genotypic data from 11 microsatellite loci scored in 250 tsetse flies collected from eight localities in Kenya. Clustering analysis identified two genetically distinct eastern and western clusters (mean between-cluster F ST = 0.202) separated by the Great Rift Valley. We also found evidence of admixture and migration between the eastern and western clusters, isolation by distance, and a widespread signal of inbreeding. We detected differences in population dynamics and dispersal patterns between the western and eastern clusters. These included lower genetic diversity (allelic richness; 7.48 versus 10.99), higher relatedness (percent related individuals; 21.4% versus 9.1%), and greater genetic differentiation (mean within-cluster F ST; 0.183 versus 0.018) in the western than the eastern cluster. Findings are consistent with the presence of smaller, less well-connected populations in Western relative to eastern Kenya. These data suggest that recent anthropogenic influences such as land use changes and vector control programs have influenced population dynamics in G. pallidipes in Kenya, and that vector control efforts should include some region-specific strategies to effectively control this disease vector.

Variants in the TNFA, IL6 and IFNG genes are associated with the dengue severity in a sample from Colombian population

Introducción. La composición genética del huésped determina, entre otros aspectos, el perfil clínico del dengue, lo cual se debería al efecto de variantes en los genes que codifican citocinas proinflamatorias.Objetivo. Evaluar la asociación entre las variantes de tres polimorfismos en los genes candidatos TNFA, IL6 e IFNG con la gravedad del dengue en una población colombiana.Materiales y métodos. Se evaluaron los polimorfismos rs1800750, rs2069843 y rs2069705 de los genes TNFA, IL6 e IFNG, respectivamente, en 226 pacientes con dengue. Los genotipos se tipificaron usando la reacción en cadena de la polimerasa (PCR) y los polimorfismos de la longitud de los fragmentos de restricción (Restriction Fragment Length Polymorphism, RFLP). Para determinar el riesgo de diferentes fenotipos del dengue, se compararon las frecuencias alélicas con la prueba de ji al cuadrado, y los genotipos y los haplotipos, con regresión logística. Por último, los análisis se ajustaron utilizando datos de autoidentificación o del componente genético ancestral.Resultados. El alelo A del rs2069843, ajustado por autoidentificación, se asoció con casos de dengue hemorrágico en afrocolombianos. En la muestra completa, dicho polimorfismo, ajustado por componente genético ancestral, fue reproducible. Además, hubo asociaciones significativas entre las combinaciones alélicas GGT y GAC de los rs1800750, rs2069843 y rs2069705 en pacientes con dengue hemorrágico, con ajuste por componente genético ancestral y sin él. Además, la combinación alélica AGC produjo 58,03 pg/ml más de interleucina 6 que la GGC, independientemente de los componentes genéticos europeo, amerindio y africano.Conclusión. Las variantes de los polimorfismos GGT y GAC de los rs1800750, rs2069843 y rs2069705 en los genes TNFA, IL6 e IFNG, respectivamente, se correlacionaron con la gravedad del dengue en esta muestra de población colombiana.