Estimates Of Genetic Differentiation Research Articles

Comparative assessment of genotyping-by-sequencing and whole-exome sequencing for estimating genetic diversity and geographic structure in small sample sizes: insights from wild jaguar populations.

Biologists currently have an assortment of high-throughput sequencing techniques allowing the study of population dynamics in increasing detail. The utility of genetic estimates depends on their ability to recover meaningful approximations while filtering out noise produced by artifacts. In this study, we empirically compared the congruence of two reduced representation approaches (genotyping-by-sequencing, GBS, and whole-exome sequencing, WES) in estimating genetic diversity and population structure using SNP markers typed in a small number of wild jaguar (Panthera onca) samples from South America. Due to its targeted nature, WES allowed for a more straightforward reconstruction of loci compared to GBS, facilitating the identification of true polymorphisms across individuals. We therefore used WES-derived metrics as a benchmark against which GBS-derived indicators were compared, adjusting parameters for locus assembly and SNP filtering in the latter. We observed significant variation in SNP call rates across samples in GBS datasets, leading to a recurrent miscalling of heterozygous sites. This issue was further amplified by small sample sizes, ultimately impacting the consistency of summary statistics between genotyping methods. Recognizing that the genetic markers obtained from GBS and WES are intrinsically different due to varying evolutionary pressures, particularly selection, we consider that our empirical comparison offers valuable insights and highlights critical considerations for estimating population genetic attributes using reduced representation datasets. Our results emphasize the critical need for careful evaluation of missing data and stringent filtering to achieve reliable estimates of genetic diversity and differentiation in elusive wildlife species.

Population Genomics of the Blue Shark, Prionace glauca, Reveals Different Populations in the Mediterranean Sea and the Northeast Atlantic.

Populations of marine top predators have been sharply declining during the past decades, and one-third of chondrichthyans are currently threatened with extinction. Sustainable management measures and conservation plans of large pelagic sharks require knowledge on population genetic differentiation and demographic connectivity. Here, we present the case of the Mediterranean blue shark (Prionace glauca, L. 1758), commonly found as bycatch in longline fisheries and classified by the IUCN as critically endangered. The management of this species suffers from a scarcity of data about population structure and connectivity within the Mediterranean Sea and between this basin and the adjacent Northeast Atlantic. Here, we assessed the genetic diversity and spatial structure of blue shark from different areas of the Mediterranean Sea and the Northeast Atlantic through genome scan analyses. Pairwise genetic differentiation estimates (F ST) on 203 specimens genotyped at 14,713 ddRAD-derived SNPs revealed subtle, yet significant, genetic differences within the Mediterranean sampling locations, and between the Mediterranean Sea and the Northeast Atlantic Ocean. Genetic differentiation suggests some degree of demographic independence between the Western and Eastern Mediterranean blue shark populations. Furthermore, results show limited genetic connectivity between the Mediterranean and the Atlantic basins, supporting the hypothesis of two distinct populations of blue shark separated by the Strait of Gibraltar. Although reproductive interactions may be limited, the faint genetic signal of differentiation suggests a recent common history between these units. Therefore, Mediterranean blue sharks may function akin to a metapopulation relying upon local demographic processes and connectivity dynamics, whereby the limited contemporary gene flow replenishment from the Atlantic may interplay with currently poorly regulated commercial catches and large-scale ecosystem changes. Altogether, these results emphasise the need for revising management delineations applied to these critically endangered sharks.

Inter and intra population genetic variability in ducks under conservation programs

This study focuses on estimation of the inter and intra population genetic variability of 6 duck populations. Microsatellite loci were used to assess the genetic variation and population structure of 6 duck populations under a conservation program in Poland. DNA polymorphism was assessed using 24 microsatellite markers and 50 individuals from each population. Polymorphism information content (PIC), heterozygosity with 2 estimators of genetic differentiation (FST and GST), and Nei's standard genetic distance were calculated. The results showed that these 6 endangered duck populations showed high genetic polymorphism. Observed heterozygosity within populations ranged from 0.14 to 0.83, with the average value of 0.58. PIC within populations ranged from 0.038 (P-8 and P-9 lines) to 0.89 (LsA line). Average number of alleles in the studied populations ranged from 4.5 (KhO-1 line) to 7.3 (LsA line). Based on the results, the pairs of lines LsA: P-33 and P-8: P-9 were found to be the most related; and the most genetically distant group was KhO-1 line, which originated as a cross between Khaki Campbell with Orpington duck.

Environmental DNA reveals the genetic diversity and population structure of an invasive species in the Laurentian Great Lakes

Environmental DNA (eDNA) has been established as a noninvasive and efficient approach to sample genetic material from aquatic environments. Although most commonly used to determine species presence and measure biodiversity, eDNA approaches also hold great potential to obtain population-level genetic information from water samples. In this study, we sequenced a panel of multiallelic microsatellite markers from filtered water and fish tissue samples to uncover patterns of intraspecific diversity in the freshwater Round Goby (Neogobius melanostomus) across their invaded range in the Laurentian Great Lakes region. Although we found that the concentration of nuclear eDNA is lower than mitochondrial eDNA, we nonetheless detected over two-thirds of all nuclear alleles identified from genotyped tissues in our eDNA samples, with the greatest recovery of common alleles in the population. Estimates of allele frequencies and genetic variability within and between populations were detected from eDNA in patterns that were consistent with individual tissue-based estimates of genetic diversity and differentiation. The strongest genetic differentiation in both eDNA and tissues exists in an isolation by distance pattern. Our study demonstrates the potential for eDNA-based approaches to characterize key population parameters required to effectively monitor, manage, or sustain aquatic species.

Genetic Diversity and Core Germplasm Research of 144 Munake Grape Resources Using 22 Pairs of SSR Markers

The Munake grape is a local variety of grape that is widely distributed in Xinjiang, China. This study aims to clarify the genetic structure of the Munake grape population, characterize genetic differentiation and gene flow among populations, gather germplasm, and establish the core germplasm collection. In total, 144 samples were collected from eight geographic populations. Twenty-two SSR markers were used to characterize the genetic diversity as well as the genetic structure of Munake grape germplasm and to establish the core germplasm collection. At each site, the average number of effective alleles (Ne) was 5.019. Overall, genetic diversity was high in the various geographic populations of Munake grapes. Polymorphic information content (PIC) ranged from 0.501 to 0.908, with an average of 0.728. Estimates of genetic differentiation and gene flow indicated that the Artux population had significant genetic differences from the other populations. Screening results indicated that a sampling proportion of 95% of the sample was required to achieve 100% allelic coverage, or a sampling proportion of 65% for 95% allelic coverage. This analysis was based on conventional genetic diversity indicators, with a core germplasm diversity index of 95% coverage. Characterization of the genetic diversity of germplasm from 144 Munake grapes not only provides valuable resources for future genetic mapping and functional genome research, but also facilitates the utilization of core germplasm and molecular breeding of Munake grapes.

Phylogeography of 912 Cherry Accessions Insight into Independent Origins of Fruiting Cherries and Domestication Footprints of Cultivated Chinese Cherry (Prunus pseudocerasus Lindl.).

The subgenus Cerasus (Rosaceae) contain numerous fruit trees and ornamentals with high economic values. The origin and genetic divergence among various types of fruiting cherries always remain a perplexing issue. We employed three plastom fragments and ITS sequence matrices derived from 912 cherry accessions to elucidate the phylogeographic structure and genetic relationship among fruiting cherries, as well as the origin and domestication of cultivated Chinese cherry. The integration of haplotype genealogies, Approximate Bayesian computation (ABC) approach and estimation of genetic differentiation within and between different groups and lineages has facilitated the resolution of several previously unresolved questions. Firstly, distant phylogenetic relationships between Cerasus and Microcerasus accessions, as indicated by both nuclear and chloroplast data, suggested independent origins and evolution for these two taxa. Moreover, two distinct geographic origin centers (Europe and China) have been confirmed, with significant phylogeographic signals and high genetic differentiation observed between cherries from these regions. This may be attributed to long-term geographic isolation caused by Himalaya-Hengduan Mountains. Our phylogeographic analyses and ABC analysis suggested that cherries inhabiting in China may have undergone multiple hybridization events during the glacial refugia of the eastern edge and southern Himalaya-Hengduan Mountains, followed by rapid radiation throughout their current habitats during interglacial period. The discrepancy between nuclear and chloroplast data may be attributed to hybridization events and incomplete lineage sorting. Furthermore, we speculated that the domesticated Chinese cherries were derived from wild accessions in Longmenshan Fault Zones approximately 2600 years ago. We have also traced the domestication processes and dispersal routes of cultivated Chinese cherries.

Comparative genomics of rainbow trout (Oncorhynchus mykiss): Is the genetic architecture of migratory behavior conserved among populations?

Rainbow trout (Oncorhynchus mykiss) are a partially migratory species wherein some individuals undergo long-distance anadromous migrations, and others stay as residents in their native freshwater streams. The decision to migrate is known to be highly heritable, and yet, the underlying genes and alleles associated with migration are not fully characterized. Here we used a pooled approach of whole-genome sequence data from migratory and resident trout of two native populations-Sashin Creek, Alaska and Little Sheep Creek, Oregon-to obtain a genome-wide perspective of the genetic architecture of resident and migratory life history. We calculated estimates of genetic differentiation, genetic diversity, and selection between the two phenotypes to locate regions of interest and then compared these associations between populations. We identified numerous genes and alleles associated with life history development in the Sashin Creek population with a notable area on chromosome 8 that may play a critical role in the development of the migratory phenotype. However, very few alleles appeared to be associated with life history development in the Little Sheep Creek system, suggesting population-specific genetic effects are likely important in the development of anadromy. Our results indicate that a migratory life history is not controlled by a singular gene or region but supports the idea that there are many independent ways for a migratory phenotype to emerge in a population. Therefore, conserving and promoting genetic diversity in migratory individuals is paramount to conserving these populations. Ultimately, our data add to a growing body of literature that suggests that population-specific genetic effects, likely mediated through environmental variation, contribute to life history development in rainbow trout.

Comparing the Utility of Microsatellites and Single Nucleotide Polymorphisms in Conservation Genetics: Insights from a Study on Two Freshwater Fish Species in France

Biodiversity is facing an unprecedented crisis and substantial efforts are needed to conserve natural populations, especially in river ecosystems. The use of molecular tools to guide conservation practices in rivers has grown in popularity over the last decades, but the amount of precision and/or biological information that would be gained by switching from the traditional short tandem repeats (STRs) to the increasingly used single nucleotide polymorphisms (SNPs) is still debated. Here, we compared the usefulness of STRs and SNPs to study spatial patterns of genetic variability in two freshwater fish species (Phoxinus dragarum and Gobio occitaniae) in southern France. SNPs were obtained from a pool-seq procedure and mapped to new genome assemblies. They provided much more precise estimates of genetic diversity and genetic differentiation than STRs, but both markers allowed the detection of very similar genetic structures in each species, which could be useful for delineating conservation units. While both markers provided similar outcomes, there were two discrepancies in genetic structures that could, nonetheless, be explained by unrecorded stocking events. Overall, we demonstrated that SNPs are not unconditionally superior to STRs in the context of large-scale riverscape genetic conservation, and that the choice of marker should primarily be based on research questions and resources available.

Genetic diversity and population structure of a Peruvian cattle herd using SNP data.

New-generation sequencing technologies, among them SNP chips for massive genotyping, are useful for the effective management of genetic resources. To date, molecular studies in Peruvian cattle are still scarce. For the first time, the genetic diversity and population structure of a reproductive nucleus cattle herd of four commercial breeds from a Peruvian institution were determined. This nucleus comprises Brahman (N = 9), Braunvieh (N = 9), Gyr (N = 5), and Simmental (N = 15) breeds. Additionally, samples from a locally adapted creole cattle, the Arequipa Fighting Bull (AFB, N = 9), were incorporated. Female individuals were genotyped with the GGPBovine100K and males with the BovineHD. Quality control, and the proportion of polymorphic SNPs, minor allele frequency, expected heterozygosity, observed heterozygosity, and inbreeding coefficient were estimated for the five breeds. Admixture, principal component analysis (PCA), and discriminant analysis of principal components (DAPC) were performed. Also, a dendrogram was constructed using the Neighbor-Joining clustering algorithm. The genetic diversity indices in all breeds showed a high proportion of polymorphic SNPs, varying from 51.42% in Gyr to 97.58% in AFB. Also, AFB showed the highest expected heterozygosity estimate (0.41 ± 0.01), while Brahman the lowest (0.33 ± 0.01). Besides, Braunvieh possessed the highest observed heterozygosity (0.43 ± 0.01), while Brahman the lowest (0.37 ± 0.02), indicating that Brahman was less diverse. According to the molecular variance analysis, 75.71% of the variance occurs within individuals, whereas 24.29% occurs among populations. The pairwise genetic differentiation estimates (FST) between breeds showed values that ranged from 0.08 (Braunvieh vs. AFB) to 0.37 (Brahman vs. Braunvieh). Similarly, pairwise Reynold's distance ranged from 0.09 (Braunvieh vs. AFB) to 0.46 (Brahman vs. Braunvieh). The dendrogram, similar to the PCA, identified two groups, showing a clear separation between Bos indicus (Brahman and Gyr) and B. taurus breeds (Braunvieh, Simmental, and AFB). Simmental and Braunvieh grouped closely with the AFB cattle. Similar results were obtained for the population structure analysis with K = 2. The results from this study would contribute to the appropriate management, avoiding loss of genetic variability in these breeds and for future improvements in this nucleus. Additional work is needed to speed up the breeding process in the Peruvian cattle system.

A review on Q ST-F ST comparisons of seed plants: Insights for conservation.

Increased access to genome-wide data provides new opportunities for plant conservation. However, information on neutral genetic diversity in a small number of marker loci can still be valuable because genomic data are not available to most rare plant species. In the hope of bridging the gap between conservation science and practice, we outline how conservation practitioners can more efficiently employ population genetic information in plant conservation. We first review the current knowledge about neutral genetic variation (NGV) and adaptive genetic variation (AGV) in seed plants, regarding both within-population and among-population components. We then introduce the estimates of among-population genetic differentiation in quantitative traits (Q ST) and neutral markers (F ST) to plant biology and summarize conservation applications derived from Q ST-F ST comparisons, particularly on how to capture most AGV and NGV on both in-situ and ex-situ programs. Based on a review of published studies, we found that, on average, two and four populations would be needed for woody perennials (n=18) to capture 99% of NGV and AGV, respectively, whereas four populations would be needed in case of herbaceous perennials (n=14). On average, Q ST is about 3.6, 1.5, and 1.1 times greater than F ST in woody plants, annuals, and herbaceous perennials, respectively. Hence, conservation and management policies or suggestions based solely on inference on F ST could be misleading, particularly in woody species. To maximize the preservation of the maximum levels of both AGV and NGV, we suggest using maximum Q ST rather than average Q ST. We recommend conservation managers and practitioners consider this when formulating further conservation and restoration plans for plant species, particularly woody species.

Microsatellite and mitochondrial DNA markers unveil the genetic structure of Nellore Palla sheep of India

The present study is carried out to genetically characterise the Nellore Palla sheep and document the genetic divergence among the three varieties of the Nellore sheep breed. Thirty two unrelated individuals from the hometract of Palla, nine each of Brown and Jodipi were sampled for the molecular studies. Fluorscent microsatellite markers were employed to genetically characterise Palla population.The analysis revealed a total of 201 alleles across the 26 loci studied. The mean observed and effective number of alleles per locus was 7.73 and 5.82 respectively. The overall observed and expected heterozygosity across 26 loci was 0.52 and 0.69 respectively. The mean FIS value in the sampled population was 0.24 with a range varying between − 0.24 (MAF209) to 0.77 (INRA063). The mean PIC value and the Shannon’s index value across the 26 loci was 0.643 and 1.49 respectively. The study concluded a considerable amount of genetic diversity with in the Nellore Palla population. The genetic divergence among the three varieties was analysed using 728 bp mitochondrial control region. The analysis for polymorphism and divergence was done with DnaSP software and the results concluded a high polymorphism within the genetic groups. The genetic differentiation estimates revealed no genetic differentiation among the three varieties of Nellore sheep breed. The phylogenetic analysis indicated that all the haplotypes of the three varieties of Nellore were clustering with lineage A except for one haplotype of Palla variety which was clustering with lineage B.

Assessing population genetic structure of three New Zealand stream insects using mitochondrial and nuclear DNA markers.

Assessing genetic differentiation among natural populations can aid understanding of dispersal patterns and connectivity among habitats. Several molecular markers have become increasingly popular in determining population genetic structure for this purpose. Here, we compared the resolution of mitochondrial cytochrome c oxidase subunit I (COI) and nuclear single nucleotide polymorphism (SNP) markers for detecting population structure among stream insects at small spatial scales. Individuals of three endemic taxa-Coloburiscus humeralis (Ephemeroptera), Zelandobius confusus (Plecoptera), and Hydropsyche fimbriata (Trichoptera)-were collected from forested streams that flow across open pasture in the North Island of New Zealand. Both COI and SNP data indicated limited population structure across the study area, and small differences observed among these species were likely related to their putative dispersal abilities. For example, fine-scale genetic differentiation between and among neighbouring stream populations for H. fimbriata suggests that gene flow, and hence dispersal, may be more limited for this species relative to the others. Based on the generally similar results provided by both types of markers, we suggest that either COI or SNP markers can provide suitable initial estimates of fine-scale population genetic differentiation in stream insects.

Comparative Phylogeography of Birds Across the Tocantins–Araguaia Interfluve Reveals a New Biogeographic Suture in the Amazon Far East

Few phylogeographic studies have focused on understanding the role of the Tocantins–Araguaia Interfluve (TAI) in the diversification of Amazonian biodiversity. Herein we evaluate the phylogeographic relationships of 14 avian lineages present in the TAI and its two neighboring Amazonian Areas of Endemism: the Xingu (XAE) and Belém (BAE). Four alternative scenarios coupling degree of genetic differentiation and area relationships were tested: (1) populations distributed in TAI, BAE, and XAE are not genetically differentiated from each other (assumed as the null hypothesis); (2) populations from TAI are more closely related to those from BAE; (3) populations from TAI are more closely related to those from XAE; and (4) TAI populations represent independent or endemic lineages not present in either the BAE or XAE. Molecular analyses considered Bayesian Inference methods and Bayesian Information Criterion (BIC) evolutionary models, haplotype sharing and genetic differentiation estimates. We found three distinct phylogeographic patterns: (i) four lineages presented no population structure across XAE, TAI and BAE; (ii) six lineages were represented in the TAI with distinct phylogroups replacing each other between XAE and BAE, but with varying degrees of contact and evidence of gene-flow within the TAI; and (iii) for four lineages, the Tocantins River acted as a barrier isolating BAE phylogroups from those inhabiting both TAI and XAE. These different patterns demonstrate a heterogeneous response to the barrier effects posed by both the Tocantins and Araguaia rivers on the local fauna. Historical geomorphological and hydromorphological factors, such as the presence and absence of paleochannels and anastomosed channel stretches and variations in sedimentation rates support a dynamic history for the Araguaia-Tocantins floodplains, likely accounting for the observed heterogeneity in species’ specific responses. Finally, the scenario of phylogeographic breaks and population subdivision recovered herein along the Tocantins and Araguaia rivers, associated with the existence of contact zones and the occurrence of gene flow, define the TAI as hitherto unknown biogeographic suture zone, localized in southeasternmost Amazonia.

Population genomics reveal low differentiation and complex demographic histories in a highly fragmented and endangered freshwater mussel

Abstract Freshwater mussels are an important element of freshwater biodiversity and provide essential ecosystem services. However, mussels are among the most imperilled groups of organisms in the world. Although research has increased in recent years, information about range‐wide genetic diversity and historical demography of most species is lacking. One such species is Cumberlandian combshell Epioblasma brevidens, which is listed as endangered under the US Endangered Species Act. Genetic diversity of E. brevidens was analysed using a high‐resolution RADseq approach and included the previously overlooked Bear Creek population. Hypotheses were tested about population decline, comparative genetic diversity and population structure with model‐based approaches enabled by a genome‐scale dataset. Estimates of genetic differentiation among populations of E. brevidens were lower than past analyses, suggesting higher historical population connectivity than previously known. Demographic analyses indicate relatively recent splits among E. brevidens populations in the late Pleistocene to early Holocene, with clear founder effects in two populations. The Clinch River population has the highest genetic diversity and effective population size, despite demographic analyses revealing decline of this population since the Pleistocene. Analyses of both population structure and migration show evidence of past gene flow, but all populations are currently isolated by artificial barriers. Analyses indicate that populations began to decline before industrialization, but fragmentation and population extirpation has been exacerbated by modern habitat destruction. Relatively high genetic diversity in the Bear Creek population indicates that water quality improvements in the last 20 years have had a positive impact on population viability, offering promise for targeted management actions. In contrast, the Big South Fork population that has been presumed stable showed low genetic diversity and effective population size. Furthermore, genetic structure among sampled populations indicates that reintroduction efforts should use broodstock from as close to the reintroduction site as possible.

High gene flow in the silverlip pearl oyster Pinctada maxima between inshore and offshore sites near Eighty Mile Beach in Western Australia.

An understanding of stock recruitment dynamics in fisheries is fundamental to successful management. Pinctada maxima is a bivalve mollusc widely distributed throughout the Indo-Pacific and is the main species targeted for cultured pearl and pearl shell production in Australia. Pearl production in Australia relies heavily on wild-caught individuals, the majority of which come from the Eighty Mile Beach region near Broome in Western Australia. In this study, we used a genotyping by sequencing approach to explore fine-scale patterns of genetic connectivity among inshore shallow and offshore deep populations of P. maxima near Eighty Mile Beach. Our results revealed high-levels of gene flow among inshore and offshore sites and no differences in genetic diversity between depths. Global estimates of genetic differentiation were low (FST = 0.006) but significantly different from zero, and pairwise estimates of genetic differentiation among sites were significant in only 3% of comparisons. Moreover, Bayesian clustering detected no separation of inshore and offshore sample sites, and instead showed all samples to be admixed among sites, locations and depths. Despite an absence of any clear spatial clustering among sites, we identified a significant pattern of isolation by distance. In a dynamic environment like Eighty Mile Beach, genetic structure can change from year-to-year and successive dispersal and recruitment events over generations likely act to homogenize the population. Although we cannot rule out the null hypothesis of panmixia, our data indicate high levels of dispersal and connectivity among inshore and offshore fishing grounds.

Genetic diversity and population structure of four Nigerian indigenous cattle breeds.

A total of thirty-eight (38) Nigerian indigenous cattle were used to evaluate the genetic relatedness, diversity, and population structure of four indigenous cattle breeds. Blood samples were collected from the experimental animals into vacutainer tubes containing EDTA, and genomic DNA extracted, quantified, and evaluated for integrity. Fourteen (14) microsatellite primers were used for polymerase chain reaction (PCR) and PCR amplification performed under standard conditions followed by electrophoresis in 2.5% Metaphor Agarose gel. Genomic parameter estimates included allele number (Na), observed (Ho) and expected (He) heterozygosity, polymorphism information content (PIC), test of Hardy-Weinberg equilibrium, and genetic diversity; pairwise Nei's genetic distance, Wright's F-statistics (FIT, FST, and FIS), and gene flow (Nm); and breed relationship, population structure, and degree of admixture. A total of 112 alleles were detected and mean number of alleles was 4.02 ± 0.190, while mean fixation index was 0.461 ± 0.068. Mean Ho and He were 0.352 ± 0.05 and 0.605 ± 0.018, respectively. Pairwise estimates of genetic differentiation, FST, were significantly different (p < 0.001) implying distinct breeds. Estimates of Nm were less than 4 but greater than 1, indicating that the cattle breeds do not belong to one panmictic population. Estimates of pairwise genetic distance revealed that White Fulani and Sokoto Gudali were more closely related than Muturu and N'Dama. The results of STRUCTURE, principal coordinate, and phylogenetic analyses revealed four clusters which implies that the breeds were genetically distinct. It is recommended that the four cattle breeds can be used to develop composites with higher genetic potentials for beef production and resistance to endemic diseases and pests. Further efforts should be made to conserve and genetically improve these breeds to meet present and future production and breeding imperatives.

Genetic Diversity, Population Structure and Phylogeny of Indigenous Goats of Mongolia Revealed by SNP Genotyping.

Simple SummaryWe performed SNP genotyping of indigenous goats of Mongolia to explore their demographic history in the global context and to estimate their genetic risks. Recently, these risks have become a subject of concern due to recent climatic disasters and uncontrolled massive breeding. Various clustering methods demonstrated close genetic relations among Mongolian, Russian, Chinese, and West Asian breeds. Mongolian goats themselves exhibited low to moderate estimates of genetic differentiation. We identified genetic features highlighting the distinct origin and breeding history of Mongolian goat breeds, as well as traces of artificial selection and geographic isolation. However, none of them met formal criteria to be considered as endangered.Mongolian goats are of great interest for studying ancient migration routes and domestication, and also represent a good model of adaptability to harsh environments. Recent climatic disasters and uncontrolled massive breeding endangered the valuable genetic resources of Mongolian goats and raised the question of their conservation status. Meanwhile, Mongolian goats have never been studied on genomic scale. We used Illumina Goat SNP50 to estimate genetic risks in five Mongolian goat breeds (Buural, Ulgii Red, Gobi GS, Erchim, Dorgon) and explored phylogenic relationships among these populations and in the context of other breeds. Various clustering methods showed that Mongolian goats grouped with other Asian breeds and were especially close to some neighboring Russian and Chinese breeds. The Buural breed showed the lowest estimates of inbreeding and exhibited the shortest genetic distances within the other Mongolian breeds, especially, to Ulgii Red and Gobi GS. These three breeds formed a single core group, being weakly differentiated from each other. Among them, Gobi GS displayed obvious signs of inbreeding probably resulted from artificial selection pressure. Dorgon and especially Erchim goats stand apart from the other Mongolian breeds according to various types of analyses, and bear unique features pointing to different breeding histories or distinct origins of these breeds. All populations showed strong decline in effective population size. However, none of them met formal criteria to be considered as endangered breeds. The SNP data obtained in this study improved the knowledge of Mongolian goat breeds and could be used in future management decisions in order to preserve their genetic diversity.

Evaluation of Maternal Genetic Background of Two Hungarian Autochthonous Sheep Breeds Coming from Different Geographical Directions.

Simple SummaryBy the 19th century, adequately producing, independent domestic animal breeds had developed in many regions around the globe. However, from the middle of the 20th century they have largely been replaced by high-performing, specialized, single-purpose cosmopolitan breeds. Breed maintenance is an activity aimed at rescuing old breeds from the threat of extinction. This process includes recording the valuable traits of a rare breed, specific diversity conservation selection, and utilization in the original production environment. Additionally, it deals with the history of breeds and the study of their genetic makeup. The aim of this paper is to evaluate the maternal genetic background of two autochthonous sheep breeds in Hungary.The aim of our research was the evaluation of the maternal genetic background of two Hungarian autochthonous sheep breeds of different geographical origin. A major argument for the preservation of endangered animal breeds is their documented past and historical importance. These also include the registration of pedigree data. This is the first study to evaluate and compare Tsigai and Cikta sheep in Hungary. Our investigation is based on two complete sequences of mitochondrial DNA (cytochrome b gene and control region). Our research was performed on these two sheep breeds with markedly different breed histories and breed characteristics to determine a possible common maternal genetic background, as ultimately the origin of both breeds can be traced back to Asia Minor. Between 2015 and 2017, a total of 203 biological samples were taken using a newly introduced founder sampling method. We found that the prevailing haplogroup B accounted for over 80% of both breeds, strengthening the common ancestral root. However, the pairwise genetic differentiation estimates (KST) calculated using the sequence-based statistics for cytochrome b gene and control region were 0.034 and 0.021, respectively (both at level p < 0.05); thus, revealing genetic differentiation in both sequences between the Tsigai and Cikta. We note that the known different history of the breeds is clearly justified by the currently studied deviations in their maternal genetic background.

Phylogeography, Population Structure, and Historical Demography of Black Drum in North America

AbstractThe Black Drum Pogonias cromis is an abundant estuarine‐dependent fish that supports recreational and commercial fisheries. We examined the phylogeography and historical demography of the Black Drum along the Atlantic and Gulf of Mexico coasts of the United States using mitochondrial DNA sequences and nuclear DNA microsatellites. Black Drum exhibited high haplotype diversity and low nucleotide diversity across its range. We observed 249 haplotypes that resolved into two haplogroups: haplogroup 1 was most abundant along the Texas coast, and haplogroup 2 was more abundant in the Atlantic Ocean and the eastern Gulf of Mexico. Analysis of molecular variance, pairwise estimates of genetic differentiation, assessment of isolation by distance, Bayesian clustering, and multivariate ordination indicated that most population structure within Black Drum was associated with a western Gulf of Mexico versus Atlantic–eastern Gulf of Mexico split. However, analysis of molecular variance, haplotype diversity, and measures of genetic differentiation provided evidence of weak population structure on the Texas coast, perhaps due to a of combination of isolation by distance and limited emigration into and dispersal from the upper and lower Laguna Madre on the lower coast. Mismatch distributions, neutrality tests, and Bayesian skyline plots suggest that haplogroup 1 underwent rapid demographic and range expansion approximately 15,000 years ago, whereas the demographic and range expansion of haplogroup 2 was more gradual and began earlier (circa 25,000–30,000 years ago). Our results suggest that Black Drum may benefit from being managed as two stocks: (1) western Gulf of Mexico and (2) Atlantic Ocean plus eastern Gulf of Mexico. Additionally, weak but significant genetic divergence in the Laguna Madre coincides with known life history divergence in this region. Caution should be used regarding the management of Black Drum in Texas as a single demographic unit.

Assessment of genetic differentiation of populations by analysis of molecular variance (analytical review)

Different approaches to using the analysis of molecular variance (AMOVA) to assess the genetic differentiation of populations have been compared in the research. Data on 11 microsatellite loci of 84 bulls of seven breeds were used. The results were compared for three options of the AMOVA module of the GenAlEx 6.502 program: the allele distance matrix (calculated FST(W&amp;C)(=θ) statistics – variant AMOVA1); the genotype distance matrix (ΦPT– AMOVA2); and the allele size difference matrix (RST– AMOVA3). Similar summary estimates of the genetic differentiation of breeds were obtained: FST(W&amp;C)= 0.108, ΦPT= 0.115, RST= 0.110 (all with pperm≤ 0.001). Between the estimates of FST(W&amp;C)and ΦPTfor each locus, the correlation coefficient was 0.99 (pvalue&lt;0.0001); no statistically significant correlations with RSTwere found. A high correlation of FST(W&amp;C)and ΦPTwith the estimates of differentiation according to Nei’s (0.96) was found. Programs other than GenAlEx (Arlequin v.3.5, GenePop v.4.7.3, RST22) gave similar AMOVA estimates. The negative linear dependence of FST(W&amp;C)and ΦPTon the level of the average heterozygosity of the breed samples was established (R2= 0.6, rS= -0.75 for pvalue &lt; 0.02) and the absence of such dependence for RST(R2= 0.04, rS= -0.23 for pvalue= 0.47). The standardization of the FST(W&amp;C)and ΦPTestimates according to Hedrick’s eliminated this dependence and raised the initial estimates to 0.35 and 0.37, respectively. The latter were comparable to the estimates obtained by the Nei-Hedrick’s (0.364-0.375), Jost’s (0.292), and Morisit-Horn’s (0.308) methods. The Mantel correlations between the matrices of paired genetic distances (GD) calculated by different measures were &gt;0.9 in most cases. The projections of the GD matrices in the principal coordinate analysis (PCoA) on the 2D plane were generally similar. The PCoA identified a cluster of Holstein «ecotypes», a cluster of «Red» breeds, and a branch of the Jersey breed. In the two-factor AMOVA of data on clusters (as two «regions»), the interregional GD was 0.357; the differentiation of breeds within the «regions» did not exceed 0.027. Modeling the association of breeds with close to zero GD resulted in an increase in the number of alleles per locus in the «new» breeds by 29 %, and an increase in the combined estimate of genetic differentiation by 29-46 %. The results obtained can be used in the development of measures for the conservation of endangered breeds.