Sex‐Linked Differentiation in Commercially Exploited Fishes: Rethinking Population Structure in Dynamic Marine Environments

Runs of homozygosity (ROHs) arise due the transmission from parents to offspring of segments that are either identical by decent (IBD) or identical by state (IBS). The former is due to consanguineous matings whereas the latter is due to demographic processes. ROHs reduce individual nucleotide diversity (θ) as a function of homozygosity, and thus ROH distributions and θ are expected to vary among species because inbreeding levels, recombination rates, and demographic histories vary widely. To help interpret genetic diversity within and among species, we utilized genome sequence data from 78 mammalian species to compare θ and ROH burden (i.e., number and length of ROHs in the genome) among groups of mammals to assess genomic signatures of inbreeding. We compared θ and ROHs: (i) among threatened and non-threatened mammals to determine the significance of contemporary conservation status; (ii) among carnivorous and non-carnivorous mammals to determine the relevance of trophic effects; (iii) relative to body size because mutation rates generally vary with body mass; and (iv) across mammals from different latitudes to test for gradients in genomic diversity (e.g., due to effects of historic climatic regimes). Our results illustrate the considerable variance in genomic diversity across mammals, and that trophic level, body mass, and latitude have significant effects on θ and ROH burden. However, conservation status was not a reliable indicator of genomic diversity. We argue that genetic or genomic diversity should be an explicit component of conservation status, as such diversity is critical to the long-term sustainability of populations, and anticipate that ROHs will become more commonly used to estimate inbreeding in wild animals.

BackgroundIn nucleus populations, regions of the genome that have a high frequency of runs of homozygosity (ROH) occur and are associated with a reduction in genetic diversity, as well as adverse effects on fitness. It is currently unclear whether, and to what extent, ROH stretches persist in the crossbred genome and how genomic management in the nucleus population might impact low diversity regions and its implications on the crossbred genome.MethodsWe calculated a ROH statistic based on lengths of 5 (ROH5) or 10 (ROH10) Mb across the genome for genotyped Landrace (LA), Large White (LW) and Duroc (DU) dams. We simulated crossbred dam (LA × LW) and market [DU × (LA × LW)] animal genotypes based on observed parental genotypes and the ROH frequency was tabulated. We conducted a simulation using observed genotypes to determine the impact of minimizing parental relationships on multiple diversity metrics within nucleus herds, i.e. pedigree-(A), SNP-by-SNP relationship matrix or ROH relationship matrix. Genome-wide metrics included, pedigree inbreeding, heterozygosity and proportion of the genome in ROH of at least 5 Mb. Lastly, the genome was split into bins of increasing ROH5 frequency and, within each bin, heterozygosity, ROH5 and length (Mb) of ROH were evaluated.ResultsWe detected regions showing high frequencies of either ROH5 and/or ROH10 across both LW and LA on SSC1, SSC4, and SSC14, and across all breeds on SSC9. Long haplotypes were shared across parental breeds and thus, regions of ROH persisted in crossbred animals. Averaged across replicates and breeds, progeny had higher levels of heterozygosity (0.0056 ± 0.002%) and lower proportion of the genome in a ROH of at least 5 Mb (−0.015 ± 0.003%) than their parental genomes when genomic relationships were constrained, while pedigree relationships resulted in negligible differences at the genomic level. Across all breeds, only genomic data was able to target low diversity regions.ConclusionsWe show that long stretches of ROH present in the parents persist in crossbred animals. Furthermore, compared to using pedigree relationships, using genomic information to constrain parental relationships resulted in maintaining more genetic diversity and more effectively targeted low diversity regions.Electronic supplementary materialThe online version of this article (doi:10.1186/s12711-016-0269-y) contains supplementary material, which is available to authorized users.

Genome-wide structural variations in Brazilian Senepol cattle, a tropically adapted taurine breed

The evolution of island populations in natural systems is driven by local adaptation and genetic drift. However, evolutionary pathways may be altered by humans in several ways. The wild boar (WB) (Sus scrofa) is an iconic game species occurring in several islands, where it has been strongly managed since prehistoric times. We examined genomic diversity at 49 803 single-nucleotide polymorphisms in 99 Sardinian WBs and compared them with 196 wild specimens from mainland Europe and 105 domestic pigs (DP; 11 breeds). High levels of genetic variation were observed in Sardinia (80.9% of the total number of polymorphisms), which can be only in part associated to recent genetic introgression. Both Principal Component Analysis and Bayesian clustering approach revealed that the Sardinian WB population is highly differentiated from the other European populations (FST=0.126–0.138), and from DP (FST=0.169). Such evidences were mostly unaffected by an uneven sample size, although clustering results in reference populations changed when the number of individuals was standardized. Runs of homozygosity (ROHs) pattern and distribution in Sardinian WB are consistent with a past expansion following a bottleneck (small ROHs) and recent population substructuring (highly homozygous individuals). The observed effect of a non-random selection of Sardinian individuals on diversity, FST and ROH estimates, stressed the importance of sampling design in the study of structured or introgressed populations. Our results support the heterogeneity and distinctiveness of the Sardinian population and prompt further investigations on its origins and conservation status.

When a population is isolated and composed of few individuals, genetic drift is the paramount evolutionary force and results in the loss of genetic diversity. Inbreeding might also occur, resulting in genomic regions that are identical by descent, manifesting as runs of homozygosity (ROHs) and the expression of recessive traits. Likewise, the genes underlying traits of interest can be revealed by comparing fixed SNPs and divergent haplotypes between affected and unaffected individuals. Populations of white-tailed deer (Odocoileus virginianus) on islands of Saint Pierre and Miquelon (SPM, France) have high incidences of leucism and malocclusions, both considered genetic defects; on the Florida Keys islands (USA) deer exhibit smaller body sizes, a polygenic trait. Here we aimed to reconstruct island demography and identify the genes associated with these traits in a pseudo case-control design. The two island populations showed reduced levels of genomic diversity and a build-up of deleterious mutations compared to mainland deer; there was also significant genome-wide divergence in Key deer. Key deer showed higher inbreeding levels, but not longer ROHs, consistent with long-term isolation. We identified multiple trait-related genes in ROHs including LAMTOR2 which has links to pigmentation changes, and NPVF which is linked to craniofacial abnormalities. Our mixed approach of linking ROHs, fixed SNPs and haplotypes matched a high number (~50) of a-priori body size candidate genes in Key deer. This suite of biomarkers and candidate genes should prove useful for population monitoring, noting all three phenotypes show patterns consistent with a complex trait and non-Mendelian inheritance.

Aim:Recessive genetic variation is thought to play a role in non-Hodgkin lymphoma (NHL) etiology. Runs of homozygosity (ROH), defined based on long, continuous segments of homozygous SNPs, can be used to estimate both measured and unmeasured recessive genetic variation. We sought to examine genome-wide homozygosity and NHL risk.Methods:We used data from eight genome-wide association studies of four common NHL subtypes: 3061 chronic lymphocytic leukemia (CLL), 3814 diffuse large B-cell lymphoma (DLBCL), 2784 follicular lymphoma (FL), and 808 marginal zone lymphoma (MZL) cases, as well as 9374 controls. We examined the effect of homozygous variation on risk by: (1) estimating the fraction of the autosome containing runs of homozygosity (FROH); (2) calculating an inbreeding coefficient derived from the correlation among uniting gametes (F3); and (3) examining specific autosomal regions containing ROH. For each, we calculated beta coefficients and standard errors using logistic regression and combined estimates across studies using random-effects meta-analysis.Results:We discovered positive associations between FROH and CLL (β = 21.1, SE = 4.41, P = 1.6 × 10−6) and FL (β = 11.4, SE = 5.82, P = 0.02) but not DLBCL (P = 1.0) or MZL (P = 0.91). For F3, we observed an association with CLL (β = 27.5, SE = 6.51, P = 2.4 × 10−5). We did not find evidence of associations with specific ROH, suggesting that the associations observed with FROH and F3 for CLL and FL risk were not driven by a single region of homozygosity.Conclusion:Our findings support the role of recessive genetic variation in the etiology of CLL and FL; additional research is needed to identify the specific loci associated with NHL risk.

The content of the 13th Mutation Detection meeting (Leiden, April 2015) is summarized in this report. Topics discussed at the meeting included current challenges of clinical NGS, advances in bioinformatics, data quality control, single cell analysis and RNA sequencing, among others. Social, ethical and regulatory challenges of genomic data handling and data sharing were the focus of an expert panel debate. The 14th International Symposium on Variants in the Genome will take place in Santiago de Compostela, June 5-8, 2017. http://isv.variome.org.

Background: Walleye (Sander vitreus), a valuable sportfish and an important ecological apex predator, exhibits genetic structuring across their range and localized structuring as a result of stocking. Methods: Walleye from 17 sampling locations across West Virginia were sequenced using a ddRAD protocol, generating various SNP datasets to assess population structuring and genomic diversity, with specific emphasis on the native Eastern Highlands strain. Different minor allele frequency filter thresholds were tested to assess impacts on genetic diversity and differentiation metrics. Results: High genetic differentiation was observed between the Eastern Highlands and Great Lakes strains, with further sub-structuring within the Eastern Highlands strain between the Ohio River populations and the other populations. Increasing MAF thresholds generally reduced the distinctiveness of clusters, but the overall inference of the number of clusters was minimally impacted. Genetic diversity metrics indicated some variability among Eastern Highlands walleye populations, with isolated populations, including the New River and Summersville Lake, showing higher inbreeding coefficients. MAF filters generally increased diversity metrics, but the trend of diversity metrics among populations remained relatively consistent. Several SNPs were found to be potentially undergoing selection, with the minor allele frequencies of these SNPs being found to be highest in Summersville Lake, highlighting potential adaptive divergence between the riverine populations and a large lentic system. Conclusions: The use of any MAF filter generated the same trends of population structuring and genomic diversity inferences regardless of the MAF threshold used. Further management of Eastern Highlands walleye in West Virginia needs to emphasize protecting the genetic integrity of the Kanawha River population and ongoing genomic screening of broodstock to conserve native genetic diversity.

Traditional Dutch chicken breeds are marginalised breeds of ornamental and cultural-historical importance. In the last decades, miniaturising of existing breeds (so called neo-bantam) has become popular and resulted in alternatives to original large breeds. However, while backcrossing is increasing the neo-bantams homozygosity, genetic exchange between breeders may increase their genetic diversity. We use the 60 K SNP array to characterise the genetic diversity, demographic history, and level of inbreeding of Dutch heritage breeds, and particularly of neo-bantams. Commercial white layers are used to contrast the impact of management strategy on genetic diversity and demography. A high proportion of alleles was found to be shared between large fowls and neo-bantams, suggesting gene flow during neo-bantams development. Population admixture analysis supports these findings, in addition to revealing introgression from neo-bantams of the same breed and of phenotypically similar breeds. The prevalence of long runs of homozygosity (ROH) confirms the importance of recent inbreeding. A high diversity in management, carried out in small breeding units explains the high heterogeneity in diversity and ROH profile displayed by traditional breeds compared to commercial lines. Population bottlenecks may explain the long ROHs in large fowls, while repetitive backcrossing for phenotype selection may account for them in neo-bantams. Our results highlight the importance of using markers to inform breeding programmes on potentially harmful homozygosity to prevent loss of genetic diversity. We conclude that bantamisation has generated unique and identifiable genetic diversity. However, this diversity can only be preserved in the near future through structured breeding programmes.

Structural variants (SVs) comprise an axis of genetic diversity with strong consequences for phenotype and fitness, making them a potentially important target for conservation genomics. Here we review how and why SVs can play a role in in conservation genomics; the different types of SVs and how they can affect phenotype; and how pangenomes and long-read sequencing are illuminating their evolution in populations, including small populations and those of conservation concern. SVs comprise multinucleotide mutations including insertions, deletions, transpositions, inversions, and other multinucleotide mutations, often overlapping genes and other functional genome regions. As a result, SVs often play important roles in phenotypic evolution and local adaptation and can contribute substantially to genetic load in inbred populations. However, our understanding of the factors influencing SV diversity in populations is still in its infancy and is complicated by the vast range of sizes, effects, and mechanisms of formation of these mutations. We argue that SVs are an important axis of genetic diversity that should be characterized alongside more traditional metrics of genetic diversity in conservation contexts. There are a number of analytical challenges to detecting and studying SVs, but analyses aimed at understanding the role of SVs in inbreeding load and population health are rapidly becoming realizable goals, accelerated by new technologies and analytical approaches. New tools, including population-scale long-read sequencing and pangenome approaches, are beginning to make SVs accessible in ways that can be readily applied in conservation settings.

BackgroundCharacterizing large genomic variants is essential to expanding the research and clinical applications of genome sequencing. While multiple data types and methods are available to detect these structural variants (SVs), they remain less characterized than smaller variants because of SV diversity, complexity, and size. These challenges are exacerbated by the experimental and computational demands of SV analysis. Here, we characterize the SV content of a personal genome with Parliament, a publicly available consensus SV-calling infrastructure that merges multiple data types and SV detection methods.ResultsWe demonstrate Parliament’s efficacy via integrated analyses of data from whole-genome array comparative genomic hybridization, short-read next-generation sequencing, long-read (Pacific BioSciences RSII), long-insert (Illumina Nextera), and whole-genome architecture (BioNano Irys) data from the personal genome of a single subject (HS1011). From this genome, Parliament identified 31,007 genomic loci between 100 bp and 1 Mbp that are inconsistent with the hg19 reference assembly. Of these loci, 9,777 are supported as putative SVs by hybrid local assembly, long-read PacBio data, or multi-source heuristics. These SVs span 59 Mbp of the reference genome (1.8%) and include 3,801 events identified only with long-read data. The HS1011 data and complete Parliament infrastructure, including a BAM-to-SV workflow, are available on the cloud-based service DNAnexus.ConclusionsHS1011 SV analysis reveals the limits and advantages of multiple sequencing technologies, specifically the impact of long-read SV discovery. With the full Parliament infrastructure, the HS1011 data constitute a public resource for novel SV discovery, software calibration, and personal genome structural variation analysis.Electronic supplementary materialThe online version of this article (doi:10.1186/s12864-015-1479-3) contains supplementary material, which is available to authorized users.

Hunan Province, located in Central-South China, has a hot and humid climate, which has shaped the unique characteristics of its cattle. In this study, we analyzed the genomic diversity of 110 indigenous Hunan cattle using whole-genome sequencing and found that they have a mixed ancestry of indicine and taurine. By grouping the cattle based on their collection regions (western, central, southeastern, and southern Hunan), we used an unsupervised three-component Gaussian model to classify the runs of homozygosity (ROH) and calculated the genomic inbreeding coefficient based on runs of homozygosity (FROH) to assess inbreeding levels. The results showed that western Hunan cattle had the highest level of hybridization, while southern Hunan cattle had the lowest. Through selective sweep analysis, we identified candidate genes and pathways related to environmental adaptation and homeostasis. Notably, the SLC5A2 gene showed strong selection signals across all four regions and exhibited a distinct haplotype compared to other referenced cattle breeds. Additionally, we detected introgression from wild species into Hunan domestic cattle and analyzed their Y-chromosome haplotypes.

ObjectiveYantai Black pig (YT), as a native population of the eastern China’s Jiaodong Peninsula of Shandong province, characterized by coarse feeding tolerance, strong disease resistance, early sexual maturity, high litter size, and superior meat quality. However, the genetic characteristics and variations underlying its crucial economic traits remain poorly understood.MethodsIn this study, we resequenced the whole genome of 17 YT individuals from distinctly different lineages breeding in three conservation farms to detect single nucleotide polymorphism (SNP) density, pairwise fixation index (FST), nucleotide diversity (π), runs of homozygosity (ROHs).ResultsOur findings revealed that YT has higher genomic diversity compared to Chinese partial indigenous pig populations and Western commercial pig populations, but lower diversity than Asian wild boars (AWB). Based on FST and values (top 1%), we identified 321 selected regions, encompassing 156 genes, between YT and AWB. Functional annotation analysis suggested that these genes are potentially responsible for growth, reproduction, and immune responses. The RBFOX3 and WDR27 genes were confirmed to be strong positively selected in YT’s breeding. Combining the results of selection sweeps and ROH islands of YT, three overlapping regions were detected. Furthermore, we found that the quantitative trait loci (QTLs) with the most overlapping regions were related to teat number, body weight, and mean corpuscular hemoglobin concentration.ConclusionWe characterized the genomic features and population structure and identified selection signals in genomic regions linked to important germplasm characteristics of YT. The insights gained from this study provide valuable references and a solid foundation for the preservation, breeding, and utilization of YT and its valuable genetic resources.

Genome-wide association studies have seen unprecedented success in identifying genetic loci that correlate with disease susceptibility and severity. Early phases of these studies have predominantly been performed in the Caucasian populations. The next phase in medical genetics is to extend the exploration across genetically diverse populations to leverage on larger sample sizes for locating smaller effects that may be present in most human populations. However, discoveries from these studies do not actually reveal the underlying functional changes to the human genome, but only point to broad regions stipulated by the extent of linkage disequilibrium (LD). Fine-mapping the functional variants can, however, be hampered by extensive LD, which can yield multiple perfect surrogates that are not distinguishable from the underlying causal variants, although several studies have illustrated the value of relying on multiple genetically diverse populations to narrow the candidate regions where the functional variants can be found in. Here, we explore the efficiency of trans-ethnic meta-analysis in discovering genetic association and in fine-mapping the causal variants by asking: are there any population diversity metrics that will be useful for: (i) identifying the populations or genomic regions where meta-analysis are likely to be more successful for discovering associations?; (ii) identifying the populations or loci to perform deep targeted sequencing for the purpose of fine-mapping causal variants? Our results indicate that simple metrics like the F(ST) or the population specificity of haplotypes are useful in trans-ethnic meta-analyses, while the degree of haplotype sharing and LD variation are informative of the efficiency in trans-ethnic fine-mapping.

A high-quality reference genome coupled with resequencing data is a promising strategy to address issues in conservation genomics. This has greatly enhanced the development of conservation plans for endangered species. Pangolins are fascinating animals with a variety of unique features. Unfortunately, they are the most trafficked wild animal in the world. In this study, we assembled a chromosome-scale genome with HiFi long reads and Hi-C short reads for the Chinese and Malayan pangolin and provided two new representative reference genomes for the pangolin species. We found a great improvement in the evaluation of genetic diversity and inbreeding based on these high-quality genomes and obtained different results for the detection of genome-wide extinction risks compared with genomes assembled using short reads. Moderate inbreeding and genetic diversity were reverified in these two pangolin species, except for one Malayan pangolin population with high inbreeding and low genetic diversity. Moreover, we identified a much higher inbreeding level (FROH = 0.54) in the Chinese pangolin individual from Taiwan Province compared with that from Mainland China, but more than 99.6% runs of homozygosity (ROH) fragments were restricted to less than 1 Mb, indicating that the high FROH in Taiwan Chinese pangolins may have accumulated from historical inbreeding events. Furthermore, our study is the first to detect relatively mild genetic purging in pangolin populations. These two high-quality reference genomes will provide valuable genetic resources for future studies and contribute to the protection and conservation of pangolins.