Dataset Information

Assessing the recovery of Y chromosome microsatellites with population genomic data using Papio and Theropithecus genomes.

ABSTRACT: Y chromosome markers can shed light on male-specific population dynamics but for many species no such markers have been discovered and are available yet, despite the potential for recovering Y-linked loci from available genome sequences. Here, we investigated how effective available bioinformatic tools are in recovering informative Y chromosome microsatellites from whole genome sequence data. In order to do so, we initially explored a large dataset of whole genome sequences comprising individuals at various coverages belonging to different species of baboons (genus: Papio) using Y chromosome references belonging to the same genus and more distantly related species (Macaca mulatta). We then further tested this approach by recovering Y-STRs from available Theropithecus gelada genomes using Papio and Macaca Y chromosome as reference sequences. Identified loci were validated in silico by a) comparing within-species relationships of Y chromosome lineages and b) genotyping male individuals in available pedigrees. Each STR was selected not to extend in its variable region beyond 100 base pairs, so that loci can be developed for PCR-based genotyping of non-invasive DNA samples. In addition to assembling a first set of Papio and Theropithecus Y-specific microsatellite markers, we released TYpeSTeR, an easy-to-use script to identify and genotype Y chromosome STRs using population genomic data which can be modulated according to available male reference genomes and genomic data, making it widely applicable across taxa.

SUBMITTER: Mutti G

PROVIDER: S-EPMC10449864 | biostudies-literature | 2023 Aug

REPOSITORIES: biostudies-literature

ACCESS DATA

Publications

Assessing the recovery of Y chromosome microsatellites with population genomic data using Papio and Theropithecus genomes.

Mutti Giacomo G Oteo-Garcia Gonzalo G Caldon Matteo M da Silva Maria Joana Ferreira MJF Minhós Tânia T Cowlishaw Guy G Gottelli Dada D Huchard Elise E Carter Alecia A Martinez Felipe I FI Raveane Alessandro A Capelli Cristian C

Scientific reports 20230824 1

Y chromosome markers can shed light on male-specific population dynamics but for many species no such markers have been discovered and are available yet, despite the potential for recovering Y-linked loci from available genome sequences. Here, we investigated how effective available bioinformatic tools are in recovering informative Y chromosome microsatellites from whole genome sequence data. In order to do so, we initially explored a large dataset of whole genome sequences comprising individual ...[more]

PMID: 37620368

Similar Datasets

Project description:Understanding the dynamics of population recovery in threatened species requires robust longitudinal monitoring datasets. However, evidence-based decision-making is often impeded by variable data collection approaches, necessitating critical evaluation of restricted available baselines. The Hainan gibbon, the world's rarest primate, had possibly declined to only seven or eight individuals in 1978 at Bawangling National Nature Reserve but has experienced subsequent population growth. Past population estimates lack detailed reporting of survey effort, and multiple conflicting estimates are available, hindering assessment of gibbon recovery. We investigated all reported estimates of Bawangling gibbon population size from 1978 to 2022, to evaluate the biological signal of population trends and the extent to which noise associated with varying survey effort, reporting and estimation may mask or misrepresent any underlying signal. This longitudinal dataset demonstrates that the Bawangling population experienced a series of bottlenecks and recoveries, with three successive periods of growth interspersed by population crashes (1978-1989, 1989-2000 and 2000-2022). The rate of gibbon population recovery was progressively slower over time in each successive period of growth, and this potential decline in recovery rate following serial bottlenecks suggests that additional management strategies may be required alongside "nature-based solutions" for this species. However, population viability analysis suggests the 1978 founder population is unlikely to have been as low as seven individuals, raising concerns for interpreting reported historical population counts and understanding the dynamics of the species' recovery. We caution against overinterpreting potential signals within "messy" conservation datasets, and we emphasise the crucial importance of standardised replicable survey methods and transparent reporting of data and effort in all future surveys of Hainan gibbons and other highly threatened species.

Project description:In the human genome, it has been estimated that considerably more sequence is under natural selection in non-coding regions [such as transcription-factor binding sites (TF-binding sites) and non-coding RNAs (ncRNAs)] compared to protein-coding ones. However, less attention has been paid to them. To study selective pressure on non-coding elements, we use next-generation sequencing data from the recently completed pilot phase of the 1000 Genomes Project, which, compared to traditional methods, allows for the characterization of a full spectrum of genomic variations, including single-nucleotide polymorphisms (SNPs), short insertions and deletions (indels) and structural variations (SVs). We develop a framework for combining these variation data with non-coding elements, calculating various population-based metrics to compare classes and subclasses of elements, and developing element-aware aggregation procedures to probe the internal structure of an element. Overall, we find that TF-binding sites and ncRNAs are less selectively constrained for SNPs than coding sequences (CDSs), but more constrained than a neutral reference. We also determine that the relative amounts of constraint for the three types of variations are, in general, correlated, but there are some differences: counter-intuitively, TF-binding sites and ncRNAs are more selectively constrained for indels than for SNPs, compared to CDSs. After inspecting the overall properties of a class of elements, we analyze selective pressure on subclasses within an element class, and show that the extent of selection is associated with the genomic properties of each subclass. We find, for instance, that ncRNAs with higher expression levels tend to be under stronger purifying selection, and the actual regions of TF-binding motifs are under stronger selective pressure than the corresponding peak regions. Further, we develop element-aware aggregation plots to analyze selective pressure across the linear structure of an element, with the confidence intervals evaluated using both simple bootstrapping and block bootstrapping techniques. We find, for example, that both micro-RNAs (particularly the seed regions) and their binding targets are under stronger selective pressure for SNPs than their immediate genomic surroundings. In addition, we demonstrate that substitutions in TF-binding motifs inversely correlate with site conservation, and SNPs unfavorable for motifs are under more selective constraints than favorable SNPs. Finally, to further investigate intra-element differences, we show that SVs have the tendency to use distinctive modes and mechanisms when they interact with genomic elements, such as enveloping whole gene(s) rather than disrupting them partially, as well as duplicating TF motifs in tandem.

Dataset Information

Assessing the recovery of Y chromosome microsatellites with population genomic data using Papio and Theropithecus genomes.

Publications

Assessing the recovery of Y chromosome microsatellites with population genomic data using Papio and Theropithecus genomes.

Similar Datasets

OmicsDI is part of the ELIXIR infrastructure

Tweets