Neutral Variation Research Articles

Immunogenetic-pathogen networks shrink in Tome’s spiny rat, a generalist rodent inhabiting disturbed landscapes

Anthropogenic disturbance may increase the emergence of zoonoses. Especially generalists that cope with disturbance and live in close contact with humans and livestock may become reservoirs of zoonotic pathogens. Yet, whether anthropogenic disturbance modifies host-pathogen co-evolutionary relationships in generalists is unknown. We assessed pathogen diversity, neutral genome-wide diversity (SNPs) and adaptive MHC class II diversity in a rodent generalist inhabiting three lowland rainforest landscapes with varying anthropogenic disturbance, and determined which MHC alleles co-occurred more frequently with 13 gastrointestinal nematodes, blood trypanosomes, and four viruses. Pathogen-specific selection pressures varied between landscapes. Genome-wide diversity declined with the degree of disturbance, while MHC diversity was only reduced in the most disturbed landscape. Furthermore, pristine forest landscapes had more functional important MHC–pathogen associations when compared to disturbed forests. We show co-evolutionary links between host and pathogens impoverished in human-disturbed landscapes. This underscores that parasite-mediated selection might change even in generalist species following human disturbance which in turn may facilitate host switching and the emergence of zoonoses.

Population genetic simulation: Benchmarking frameworks for non-standard models of natural selection.

Population genetic simulation has emerged as a common tool for investigating increasingly complex evolutionary and demographic models. Software capable of handling high-level model complexity has recently been developed, and the advancement of tree sequence recording now allows simulations to merge the efficiency and genealogical insight of coalescent simulations with the flexibility of forward simulations. However, frameworks utilizing these features have not yet been compared and benchmarked. Here, we evaluate various simulation workflows using the coalescent simulator msprime and the forward simulator SLiM, to assess resource efficiency and determine an optimal simulation framework. Three aspects were evaluated: (1) the burn-in, to establish an equilibrium level of neutral diversity in the population; (2) the forward simulation, in which temporally fluctuating selection is acting; and (3) the final computation of summary statistics. We provide typical memory and computation time requirements for each step. We find that the fastest framework, a combination of coalescent and forward simulation with tree sequence recording, increases simulation speed by over twenty times compared to classical forward simulations without tree sequence recording, although it does require six times more memory. Overall, using efficient simulation workflows can lead to a substantial improvement when modelling complex evolutionary scenarios-although the optimal framework ultimately depends on the available computational resources.

Local adaptation with gene flow in a highly dispersive shark.

Adaptive divergence in response to environmental clines are expected to be common in species occupying heterogeneous environments. Despite numerous advances in techniques appropriate for non-model species, gene-environment association studies in elasmobranchs are still scarce. The bronze whaler or copper shark (Carcharhinus brachyurus) is a large coastal shark with a wide distribution and one of the most exploited elasmobranchs in southern Africa. Here, we assessed the distribution of neutral and adaptive genomic diversity in C. brachyurus across a highly heterogeneous environment in southern Africa based on genome-wide SNPs obtained through a restriction site-associated DNA method (3RAD). A combination of differentiation-based genome-scan (outflank) and genotype-environment analyses (redundancy analysis, latent factor mixed models) identified a total of 234 differentiation-based outlier and candidate SNPs associated with bioclimatic variables. Analysis of 26,299 putatively neutral SNPs revealed moderate and evenly distributed levels of genomic diversity across sites from the east coast of South Africa to Angola. Multivariate and clustering analyses demonstrated a high degree of gene flow with no significant population structuring among or within ocean basins. In contrast, the putatively adaptive SNPs demonstrated the presence of two clusters and deep divergence between Angola and all other individuals from Namibia and South Africa. These results provide evidence for adaptive divergence in response to a heterogeneous seascape in a large, mobile shark despite high levels of gene flow. These results are expected to inform management strategies and policy at the national and regional level for conservation of C. brachyurus populations.

Power Quality Impacts of Grid-Tied PV Inverters on Low Voltage Distribution Networks a Smart OpenDSS Model to Find Power Quality Threshold Limits

As a renewable initiative, risen attention can be seen in solar photovoltaic (PV) installations in the world. Since tropical countries and even many other countries during summer experience higher irradiance levels, there is great potential for PV generation in solar panels. However, power quality can be impacted when the number of PV modules is increased in a network due to the emergence of harmonics, unbalanced voltages, voltage flicker, neutral voltage variations, etc. Therefore, it is essential to empirically analyse power quality parameters, i.e., total harmonic distortion (THD), voltage unbalance, etc., and find threshold margins of PV capacities in a network. Based on a case study conducted in Negombo, Sri Lanka, this work determines power quality impacts with reference to the number of PV interconnections in a distributed network. For the task, a low-voltage distribution network was chosen which was fed by a 250kVA, 11/0.4kV transformer with domestic loads and grid-tied PV inverters. A simulation model was developed in Open Distribution System Simulator (OpenDSS) with time-varying load patterns and PV generation. Consequently, a smart efficient method was introduced to model domestic loads with unique time-varying demand patterns. A determination criterion was established to derive snapshot load flow instants using time-varying load flow results to analyse voltage profiles along distribution feeders. Then, the model was enhanced to quantify the power quality parameters such as individual harmonic content, THD of voltage as well as current and neutral voltage variation. The results reveal that node voltage has improved with PV interception without violating the upper limit. The THD of voltage and current have slightly increased with the addition of PV inverters.

Weak Genetic Isolation and Putative Phenotypic Selection in the Wild Carnation Dianthus virgineus (Caryophyllaceae).

By relating genetic divergence at neutral loci, phenotypic variation, and geographic and environmental distances, it is possible to dissect micro-evolutionary scenarios involving natural selection and neutral evolution. In this work, we tested the patterns of intraspecific genetic and phenotypic variation along an elevational gradient, using Dianthus virgineus as study system. We genotyped genome-wide SNPs through ddRAD sequencing and quantified phenotypic variation through multivariate morphological variation. We assessed patterns of variation by testing the statistical association between genetic, phenotypic, geographic, and elevational distances and explored the role of genetic drift and selection by comparing the Fst and Pst of morphometric traits. We revealed a weak genetic structure related to geographic distance among populations, but we excluded the predominant role of genetic drift acting on phenotypic traits. A high degree of phenotypic differentiation with respect to genetic divergence at neutral loci allowed us to hypothesize the effect of selection, putatively fuelled by changing conditions at different sites, on morphological traits. Thus, natural selection acting despite low genetic divergence at neutral loci can be hypothesized as a putative driver explaining the observed patterns of variation.

Functional genomic diversity is correlated with neutral genomic diversity in populations of an endangered rattlesnake

Theory predicts that genetic erosion in small, isolated populations of endangered species can be assessed using estimates of neutral genetic variation, yet this widely used approach has recently been questioned in the genomics era. Here, we leverage a chromosome-level genome assembly of an endangered rattlesnake (Sistrurus catenatus) combined with whole genome resequencing data (N = 110 individuals) to evaluate the relationship between levels of genome-wide neutral and functional diversity over historical and future timescales. As predicted, we found positive correlations between genome-wide estimates of neutral genetic diversity (π) and inferred levels of adaptive variation and an estimate of inbreeding mutation load, and a negative relationship between neutral diversity and an estimate of drift mutation load. However, these correlations were half as strong for projected future levels of neutral diversity based on contemporary effective population sizes. Broadly, our results confirm that estimates of neutral genetic diversity provide an accurate measure of genetic erosion in populations of a threatened vertebrate. They also provide nuance to the neutral-functional diversity controversy by suggesting that while these correlations exist, anthropogenetic impacts may have weakened these associations in the recent past and into the future.

Demographic modeling of admixed Latin American populations from whole genomes

Demographic models of Latin American populations often fail to fully capture their complex evolutionary history, which has been shaped by both recent admixture and deeper-in-time demographic events. To address this gap, we used high-coverage whole-genome data from Indigenous American ancestries in present-day Mexico and existing genomes from across Latin America to infer multiple demographic models that capture the impact of different timescales on genetic diversity. Our approach, which combines analyses of allele frequencies and ancestry tract length distributions, represents a significant improvement over current models in predicting patterns of genetic variation in admixed Latin American populations. We jointly modeled the contribution of European, African, East Asian, and Indigenous American ancestries into present-day Latin American populations. We infer that the ancestors of Indigenous Americans and East Asians diverged ∼30 thousand years ago, and we characterize genetic contributions of recent migrations from East and Southeast Asia to Peru and Mexico. Our inferred demographic histories are consistent across different genomic regions and annotations, suggesting that our inferences are robust to the potential effects of linked selection. In conjunction with published distributions of fitness effects for new nonsynonymous mutations in humans, we show in large-scale simulations that our models recover important features of both neutral and deleterious variation. By providing a more realistic framework for understanding the evolutionary history of Latin American populations, our models can help address the historical under-representation of admixed groups in genomics research and can be a valuable resource for future studies of populations with complex admixture and demographic histories.

Lidar Observations of Predawn Thermosphere‐Ionosphere Na (TINa) Layers Over Boulder (40.13°N, 105.24°W): Annual Phase Variations and Correlation With Sunrise and Tidal Winds

AbstractWe have discovered that the peak phase time of predawn thermosphere‐ionosphere Na (TINa) layers (∼110–150 km altitude) undergoes clear annual variations with the earliest occurrence in summer and latest in winter over Boulder (40.13°N, 105.24°W), which are closely correlated to annual phase variations of sunrise and tidal winds. Such discoveries were enabled by the first characterization of 12 monthly composites of TINa layers from January through December using 7 years of lidar observations (2011–2017). Despite their tenuous densities, the predawn TINa layers have nearly 100% occurrence rate (160 out of 164 nights of observations). Monthly composites show downward‐phase‐progression TINa descending at similar phase speeds as Climatological Tidal Model of the Thermosphere tidal winds. These TINa layers occur in ion convergence but neutral divergence regions, modeled using tidal winds. These results support the formation mechanism (neutralization of converged TINa+ forming TINa) proposed previously and suggest that migrating tidal winds experience annual phase variations.

Divergence in the internal genital morphology of females and correlated divergence in male intromittent structures among populations of a millipede.

Our understanding of genital evolution comes largely from studies of male genitalia. Females have received far less attention because of the difficulties inherent in quantifying the shapes of their internal genital structures. Here we combine advances in micro-computed tomography with a new landmark free method of quantifying three-dimensional trait shape, to document patterns of divergence in female genital shape, and the correlated divergence of male genitalia among populations of the millipede Antichiropus variabilis. We used single-nucleotide polymorphisms to estimate levels of neutral genetic divergence among seven populations of millipede. Genetic divergence was high and correlated with geographic distance. Comparing phenotypic divergence in genital shape to neutral genetic divergence, we found that genital shape for both females and males has diverged more than would be expected from random drift, consistent with a pattern of directional selection. While there was significant covariation between female and male genital shape across populations, the magnitude of divergence in genital shape between the sexes was not correlated. Our results demonstrate the utility of using three-dimensional scanning technologies to examine female genital traits and add to a small but growing number of studies showing that like male genitalia, female genitalia can be under strong directional selection.

Signs of rapid evolution in an invasive forest species: Drivers of the incipient neutral, adaptive and phenotypic divergence

Evolutionary forces such as genetic drift and natural selection operate during the process of biological invasion. More specifically, nonnative populations suffer genetic and demographic bottlenecks, and are exposed to new environmental, climatic and biotic filters. Thus, neutral and/or adaptive drivers may contribute to genetic differentiation of the introduced germplasm. We hypothesize that the gene pool of fire-adapted species that become invasive in fire-prone ecosystems may differ from that of originally planted and native one because of novel selective forces and founder effects. Genetic impoverishment due to genetic drift could also be buffered by demographic processes as the production of a large aerial seed bank by serotiny. We conducted genomic analyses of the world’s most planted tree Pinus radiata using Single Nucleotide Polymorphisms, comparing introduced cohorts with and without fire exposure, and also with one population of the native range. Phenotypic traits of introduced cohorts were assessed for quantitative characters. In the introduced range we found lower expected heterozygosity and nucleotide diversity due to bottlenecks during introduction and commercial selection in orchards. We also detected decreased inbreeding, related to a serotinous seed bank that buffers bottlenecks. We found genomic regions under selection, both in the comparison between native and introduced range, and between introduced cohorts with and without exposure to fire. Within the introduced area, we observed possible fire-mediated diversifying selection for five quantitative traits. Our results show that both neutral and adaptive forces promote early divergence of the introduced species Pinus radiata. Also, fire acts as a rapid selection factor, prompting the incipient adaptive divergence of cohorts recruited after fires, facilitating invasion.

The impact of habitat loss and population fragmentation on genomic erosion

Habitat loss and population fragmentation pose severe threats to biodiversity and the survival of many species. Population isolation and the decline in effective population size lead to increased genetic drift and inbreeding. In turn, this reduces neutral diversity, and it also affects the genetic load of deleterious mutations. Here, we analyse the effect of such genomic erosion by designing a spatially explicit, individual based model in SLiM, simulating the effects of the recorded habitat loss in Mauritius over the past ~ 250 years. We show that the loss of neutral diversity (genome-wide heterozygosity) was barely noticeable during the first 100 years of habitat loss. Changes to the genetic load took even more time to register, and they only became apparent circa 200 years after the start of habitat decline. Although a considerable number of deleterious mutations were lost by drift, others increased in frequency. The masked load was thus converted into a realised load, which compromised individual fitness and population viability after much of the native habitat had been lost. Importantly, genomic erosion continued after the metapopulation had stabilised at low numbers. Our study shows that historic habitat loss can pose a sustained threat to populations also in future generations, even without further habitat loss. The UN’s Decade on Ecosystem Restoration needs to lead to transformative change to save species from future extinction, and this requires the urgent restoration of natural habitats.

Genetic drift versus natural selection affecting the evolution of spectral and functional traits of two key macrophytes: Phragmites australis and Nuphar lutea

Abstract Both genetic and phenotypic intraspecific diversity play a crucial role in the ecological and evolutionary dynamics of organisms. Several studies have compared phenotypic divergence (Pst) and differentiation of neutral loci (Fst) to infer the relative roles of genetic drift and natural selection in population differentiation (Pst–Fst comparison). For the first time, we have assessed and compared the genetic variation and differentiation at the leaf trait level in two key macrophytes, Phragmites australis and Nuphar lutea. To this aim, we quantified and described the genetic structure and phenotypic diversity of both species in five lake systems in north‐central Italy. We then investigated the relative roles of genetic drift and natural selection on leaf trait differentiation (Pst–Fst), assuming that Fst reflects divergence caused only by genetic drift while Pst also incorporates the effects of selective dynamics on the phenotype. In terms of genetic structure, the results for P. australis were in line with those observed for other Italian and European conspecific populations. Conversely, N. lutea showed a more complex genetic structure than expected at the site level, probably due to the combined effect of genetic isolation and its mixed mating system. Both species exhibited high variability in leaf functional traits within and among sites, highlighting a high degree of phenotypic plasticity. Pst–Fst comparisons showed a general tendency towards directional selection in P. australis and a more complex pattern in N. lutea. Indeed, the drivers of phenotypic differentiation in N. lutea showed a variable mix of stabilising and directional selection or neutral divergence at most sites. The prevalence of vegetative over generative reproduction leads P. australis populations to be dominated by a few clones that are well adapted to local conditions, including phenotypes that respond plastically to the environment. In contrast, in N. lutea, the interaction of a mixed mating system and geographical isolation among distant sites tends to reduce the effect of outbreeding depression and provides the genetic basis for adaptive capacity. The first joint analysis of the genetic structure of these two key macrophytes allowed a better understanding of the relative roles of genetic drift and natural selection in the diversification of phenotypic traits within habitats dominated by P. australis and N. lutea.

Inferring human neutral genetic variation from craniodental phenotypes.

There is a growing consensus that global patterns of modern human cranial and dental variation are shaped largely by neutral evolutionary processes, suggesting that craniodental features can be used as reliable proxies for inferring population structure and history in bioarchaeological, forensic, and paleoanthropological contexts. However, there is disagreement on whether certain types of data preserve a neutral signature to a greater degree than others. Here, we address this unresolved question and systematically test the relative neutrality of four standard metric and nonmetric craniodental data types employing an extensive computational genotype-phenotype comparison across modern populations from around the world. Our computation draws on the largest existing data sets currently available, while accounting for geographically structured environmental variation, population sampling uncertainty, disparate numbers of phenotypic variables, and stochastic variation inherent to a neutral model of evolution. Our results reveal that the four data types differentially capture neutral genomic variation, with highest signals preserved in dental nonmetric and cranial metric data, followed by cranial nonmetric and dental metric data. Importantly, we demonstrate that combining all four data types together maximizes the neutral genetic signal compared with using them separately, even with a limited number of phenotypic variables. We hypothesize that this reflects a lower level of genetic integration through pleiotropy between, compared to within, the four data types, effectively forming four different modules associated with relatively independent sets of loci. Therefore, we recommend that future craniodental investigations adopt holistic combined data approaches, allowing for more robust inferences about underlying neutral genetic variation.

Broad-scale variation in human genetic diversity levels is predicted by purifying selection on coding and non-coding elements.

Analyses of genetic variation in many taxa have established that neutral genetic diversity is shaped by natural selection at linked sites. Whether the mode of selection is primarily the fixation of strongly beneficial alleles (selective sweeps) or purifying selection on deleterious mutations (background selection) remains unknown, however. We address this question in humans by fitting a model of the joint effects of selective sweeps and background selection to autosomal polymorphism data from the 1000 Genomes Project. After controlling for variation in mutation rates along the genome, a model of background selection alone explains ~60% of the variance in diversity levels at the megabase scale. Adding the effects of selective sweeps driven by adaptive substitutions to the model does not improve the fit, and when both modes of selection are considered jointly, selective sweeps are estimated to have had little or no effect on linked neutral diversity. The regions under purifying selection are best predicted by phylogenetic conservation, with ~80% of the deleterious mutations affecting neutral diversity occurring in non-exonic regions. Thus, background selection is the dominant mode of linked selection in humans, with marked effects on diversity levels throughout autosomes.

Sex differences in deleterious genetic variants in a haplodiploid social insect.

Deleterious variants are selected against but can linger in populations at low frequencies for long periods of time, decreasing fitness and contributing to disease burden in humans and other species. Deleterious variants occur at low frequency but distinguishing deleterious variants from low-frequency neutral variation is challenging based on population genomics data alone. As a result, we have little sense of the number and identity of deleterious variants in wild populations. For haplodiploid species, it has been hypothesised that deleterious alleles will be directly exposed to selection in haploid males, but selection can be masked in diploid females when deleterious variants are recessive, resulting in more efficient purging of deleterious mutations in males. Therefore, comparisons of the differences between haploid and diploid genomes from the same population may be a useful method for inferring rare deleterious variants. This study provides the first formal test of this hypothesis. Using wild populations of Northern paper wasps (Polistes fuscatus), we find that males have fewer missense and nonsense variants per generation than females from the same population. Allele frequency differences are especially pronounced for rare missense and nonsense variants and these differences lead to a lower mutational load in males than females. Based on these data we infer that many highly deleterious mutations are segregating in the paper wasp population. Stronger selection against deleterious alleles in haploid males may have implications for adaptation in other haplodiploid insects and provides evidence that wild populations harbour abundant deleterious variants.

The influence of explicit local dynamics on range expansions driven by long-range dispersal.

Range expansions are common in natural populations. They can take such forms as an invasive species spreading into a new habitat or a virus spreading from host to host during a pandemic. When the expanding species is capable of dispersing offspring over long distances, population growth is driven by rare but consequential long-range dispersal events that seed satellite colonies far from the densely occupied core of the population. These satellites accelerate growth by accessing unoccupied territory, and also act as reservoirs for maintaining neutral genetic variation present in the originating population, which would ordinarily be lost to drift. Prior theoretical studies of dispersal-driven expansions have shown that the sequential establishment of satellites causes initial genetic diversity to be either lost or maintained to a level determined by the breadth of the distribution of dispersal distances. If the tail of the distribution falls off faster than a critical threshold, diversity is steadily eroded over time; by contrast, broader distributions with a slower falloff allow some initial diversity to be maintained for arbitrarily long times. However, these studies used lattice-based models and assumed an instantaneous saturation of the local carrying capacity after the arrival of a founder. Real-world populations expand in continuous space with complex local dynamics, which potentially allow multiple pioneers to arrive and establish within the same local region. Here, we evaluate the impact of local dynamics on the population growth and the evolution of neutral diversity using a computational model of range expansions with long-range dispersal in continuous space, with explicit local dynamics that can be controlled by altering the mix of local and long-range dispersal events. We found that many qualitative features of population growth and neutral genetic diversity observed in lattice-based models are preserved under more complex local dynamics, but quantitative aspects such as the rate of population growth, the level of maintained diversity, and the rate of decay of diversity all depend strongly on the local dynamics. Besides identifying situations in which modeling the explicit local population dynamics becomes necessary to understand the population structure of jump-driven range expansions, our results show that local dynamics affects different features of the population in distinct ways, and can be more or less consequential depending on the degree and form of long-range dispersal as well as the scale at which the population structure is measured.

Improving the use of plant genetic resources to sustain breeding programs’ efficiency

Genetic progress of crop plants is required to face human population growth and guarantee production stability in increasingly unstable environmental conditions. Breeding is accompanied by a loss in genetic diversity, which hinders sustainable genetic gain. Methodologies based on molecular marker information have been developed to manage diversity and proved effective in increasing long-term genetic gain. However, with realistic plant breeding population sizes, diversity depletion in closed programs appears ineluctable, calling for the introduction of relevant diversity donors. Although maintained with significant efforts, genetic resource collections remain underutilized, due to a large performance gap with elite germplasm. Bridging populations created by crossing genetic resources to elite lines prior to introduction into elite programs can manage this gap efficiently. To improve this strategy, we explored with simulations different genomic prediction and genetic diversity management options for a global program involving a bridging and an elite component. We analyzed the dynamics of quantitative trait loci fixation and followed the fate of allele donors after their introduction into the breeding program. Allocating 25% of total experimental resources to create a bridging component appears highly beneficial. We showed that potential diversity donors should be selected based on their phenotype rather than genomic predictions calibrated with the ongoing breeding program. We recommend incorporating improved donors into the elite program using a global calibration of the genomic prediction model and optimal cross selection maintaining a constant diversity. These approaches use efficiently genetic resources to sustain genetic gain and maintain neutral diversity, improving the flexibility to address future breeding objectives.

Genomic variation across Chinook salmon populations reveals effects of a duplication on migration alleles and supports fine scale structure.

The distribution of ecotypic variation in natural populations is influenced by neutral and adaptive evolutionary forces that are challenging to disentangle. This study provides a high-resolution portrait of genomic variation in Chinook salmon (Oncorhynchus tshawytscha) with emphasis on a region of major effect for ecotypic variation in migration timing. With a filtered data set of ~13million single nucleotide polymorphisms (SNPs) from low-coverage whole genome resequencing of 53 populations (3566 barcoded individuals), we contrasted patterns of genomic structure within and among major lineages and examined the extent of a selective sweep at a major effect region underlying migration timing (GREB1L/ROCK1). Neutral variation provided support for fine-scale structure of populations, while allele frequency variation in GREB1L/ROCK1 was highly correlated with mean return timing for early and late migrating populations within each of the lineages (r2 = .58-.95; p< .001). However, the extent of selection within the genomic region controlling migration timing was much narrower in one lineage (interior stream-type) compared to the other two major lineages, which corresponded to the breadth of phenotypic variation in migration timing observed among lineages. Evidence of a duplicated block within GREB1L/ROCK1 may be responsible for reduced recombination in this portion of the genome and contributes to phenotypic variation within and across lineages. Lastly, SNP positions across GREB1L/ROCK1 were assessed for their utility in discriminating migration timing among lineages, and we recommend multiple markers nearest the duplication to provide highest accuracy in conservation applications such as those that aim to protect early migrating Chinook salmon. These results highlight the need to investigate variation throughout the genome and the effects of structural variants on ecologically relevant phenotypic variation in natural species.

Where to draw the line? Expanding the delineation of conservation units to highly mobile taxa.

Conservation units (CUs) are an essential tool for maximizing evolutionary potential and prioritizing areas across a species' range for protection when implementing conservation and management measures. However, current workflows for identifying CUs on the basis of neutral and adaptive genomic variation largely ignore information contained in patterns of isolation by distance (IBD), frequently the primary signal of population structure in highly mobile taxa, such as birds, bats, and marine organisms with pelagic larval stages. While individuals located on either end of a species' distribution may exhibit clear genetic, phenotypic, and ecological differences, IBD produces subtle changes in allele frequencies across space, making it difficult to draw clear boundaries for conservation purposes in the absence of discrete population structure. Here, we highlight potential pitfalls that arise when applying common methods for delineating CUs to continuously distributed organisms and review existing methods for detecting subtle breakpoints in patterns of IBD that can indicate barriers to gene flow in highly mobile taxa. In addition, we propose a new framework for identifying CUs in all organisms, including those characterized by continuous genomic differentiation, and suggest several possible ways to harness the information contained in patterns of IBD to guide conservation and management decisions.

Characterization of Cow, Goat, and Water Buffalo Milk Fat Globule Lipids by High-Performance Thin Layer Chromatography

Ruminant milk is an essential part of the human diet and is widely accepted as a major nutrient source in developing countries. However, the polar and neutral lipid content variation in milk fat globules (MFG)among cow, goat, and water buffalo is poorly understood. This study used high-performance thin layer chromatography to identify and quantify five major polar (PL) and three neutral lipids (NL) from the MFG of cow, goat, and water buffalo. Optimal separation was achieved for PLs using chloroform: methanol: water (65:25:4), and hexane: diethyl ether: acetic acid (70:30:1) for NLs. The lower detectable (0.12 to 1.53 μg/mL) and quantification (0.12 to 1.53 μg/mL) limits indicated the high sensitivity of the method. Quantification at 540 nm showed the highest abundance of phosphatidylethanolamine and triglycerides. Fat globules were further characterized for size and microstructural properties, which revealed smaller globules in goats (0.99 ± 0.04 μm) than cows (1.85 ± 0.03 μm) and water buffaloes (2.91 ± 0.08 μm), indicating a negative correlation with PL but a positive correlation with NL. The variation in lipid quantity among different animal species suggests more research to support their selection as a suitable source for developing functional food to impact human health positively.