Depauperate genetic variability detected in the American and European bison using genomic techniques

A reduced SNP panel optimised for non-invasive genetic assessment of a genetically impoverished conservation icon, the European bison

The European bison was saved from the brink of extinction due to considerable conservation efforts since the early twentieth century. The current global population of > 9500 individuals is the result of successful ex situ breeding based on a stock of only 12 founders, resulting in an extremely low level of genetic variability. Due to the low allelic diversity, traditional molecular tools, such as microsatellites, fail to provide sufficient resolution for accurate genetic assessments in European bison, let alone from non-invasive samples. Here, we present a SNP panel for accurate high-resolution genotyping of European bison, which is suitable for a wide variety of sample types. The panel accommodates 96 markers allowing for individual and parental assignment, sex determination, breeding line discrimination, and cross-species detection. Two applications were shown to be utilisable in further Bos species with potential conservation significance. The new SNP panel will allow to tackle crucial tasks in European bison conservation, including the genetic monitoring of reintroduced populations, and a molecular assessment of pedigree data documented in the world's first studbook of a threatened species.

Ancient Faunal History Revealed by Interdisciplinary Biomolecular Approaches

Starting four decades ago, studies have examined the ecology and evolutionary dynamics of populations and species using short mitochondrial DNA fragments and stable isotopes. Through technological and analytical advances, the methods and biomolecules at our disposal have increased significantly to now include lipids, whole genomes, proteomes, and even epigenomes. At an unprecedented resolution, the study of ancient biomolecules has made it possible for us to disentangle the complex processes that shaped the ancient faunal diversity across millennia, with the potential to aid in implicating probable causes of species extinction and how humans impacted the genetics and ecology of wild and domestic species. However, even now, few studies explore interdisciplinary biomolecular approaches to reveal ancient faunal diversity dynamics in relation to environmental and anthropogenic impact. This review will approach how biomolecules have been implemented in a broad variety of topics and species, from the extinct Pleistocene megafauna to ancient wild and domestic stocks, as well as how their future use has the potential to offer an enhanced understanding of drivers of past faunal diversity on Earth.

In Search of Species-Specific SNPs in a Non-Model Animal (European Bison (Bison bonasus))-Comparison of De Novo and Reference-Based Integrated Pipeline of STACKS Using Genotyping-by-Sequencing (GBS) Data

The European bison is a non-model organism; thus, most of its genetic and genomic analyses have been performed using cattle-specific resources, such as BovineSNP50 BeadChip or Illumina Bovine 800 K HD Bead Chip. The problem with non-specific tools is the potential loss of evolutionary diversified information (ascertainment bias) and species-specific markers. Here, we have used a genotyping-by-sequencing (GBS) approach for genotyping 256 samples from the European bison population in Bialowieza Forest (Poland) and performed an analysis using two integrated pipelines of the STACKS software: one is de novo (without reference genome) and the other is a reference pipeline (with reference genome). Moreover, we used a reference pipeline with two different genomes, i.e., Bos taurus and European bison. Genotyping by sequencing (GBS) is a useful tool for SNP genotyping in non-model organisms due to its cost effectiveness. Our results support GBS with a reference pipeline without PCR duplicates as a powerful approach for studying the population structure and genotyping data of non-model organisms. We found more polymorphic markers in the reference pipeline in comparison to the de novo pipeline. The decreased number of SNPs from the de novo pipeline could be due to the extremely low level of heterozygosity in European bison. It has been confirmed that all the de novo/Bos taurus and Bos taurus reference pipeline obtained SNPs were unique and not included in 800 K BovineHD BeadChip.

New data on the genetic diversity of European bison Bison bonasus (Linnaeus, 1758) in Belarus

The paper presents data on the assessment of the genetic diversity of five subpopulations of the European bison (Bison bonasus) in Belarus — from the National Park “Belovezhskaya Pushcha”, the National Park “Pripyatsky”, the Osipovichi district (Mogilev area), SEI “Berezinsky Biosphere Reserve” and the Grodno region. In general, the work includes 30 samples of muscle tissue from the collection of Gene bank of wild fauna in SSPA “SPC NAS of Belarus on Bioresources” (Minsk, Belarus). Microsatellites were used as markers to assess genetic diversity, structure, and search for signs of a sharp decline in the size of bison subpopulations in the past. A total of 11 microsatellite markers were used, recommended by the Food and Agriculture Organization of the United Nations for cattle research. The analysis of B. bonasus subpopulation from the NP “Pripyatsky” showed signs of passing through the genetic bottleneck. All studied subpopulations are characterized by a similarly low genetic diversity level in all analyzed indicators (mean number of alleles, allelic diversity, observed and expected heterozygosity). The expected heterozygosity (He) for the three subpopulations from the NP “Belovezhskaya Pushcha”, the NP “Pripyatsky” and from the Osipovichi district ranged from 0.37 to 0.39. For the studied subpopulations, the values of the fixation index were negative. The assessment of the presence of genetic structuring between the subpopulations of bison from the NP “Belovezhskaya Pushcha”, the NP “Pripyatsky” and from the Osipovichi district based on the values of such indexes as Fst and DJost which showed no signs of genetic differentiation, which is also confirmed by principal coordinates analysis (PCoA). The European bison conservation in Belarus has required tremendous efforts in the past. So far, even though the impressively large population size reached in Belarus, B. bonasus status still should not be considered as stable, which is closely linked to aspects of its overall low genetic diversity. Our research confirmed the low genetic variability of Belarusian subpopulations. Therefore, the more extensive research concentrated on identifying genetic diversity is necessary to ensure the beneficial control of gene flow and register a potential correlation of unfavorable gene variants with possible inbreeding depression. These attempts are required to lay the groundwork for the management and protection of the European bison in Belarus.

Prevalence of Ineffective Haplotypes at the Rice Blast Resistance (R) Gene Loci in Chinese Elite Hybrid Rice Varieties Revealed by Sequence-Based Molecular Diagnosis

Multiple haplotypes at the same rice blast R-gene locus share extremely high sequence similarity, which makes the gene diagnostic method using molecular markers less effective in differentiation from one another. The composition and distribution pattern of deployed R genes/haplotypes in elite rice varieties has not been extensively analyzed. In this study, we employed PCR amplification and sequencing approach for the diagnosis of R-gene haplotypes in 54 Chinese elite rice varieties. A varied number of functional and nonfunctional haplotypes of 4 target major R-gene loci, i.e., Pi2/9, Pi5, Pik, and Pib, were deduced by referring to the reference sequences of known R genes. Functional haplotypes accounted for relatively low frequencies for the Pi2/9 (15%) and Pik (9%) loci but for relatively high frequencies for the Pi5 (50%) and Pib (54%) loci. Intriguingly, significant frequencies of 33%, 39%, 46% of non-functional haplotypes at the Pi2/9, Pik, and Pib loci, respectively, with traceable original donors were identified, suggesting that they were most likely unintentionally spread by using undesirable donors in various breeding programs. In the case of Pi5 locus, only a single haplotype, i.e., Pi5 was identified. The reactions of 54 rice varieties to the differential isolates were evaluated, which showed a good correlation to the frequency of cognate avirulence (Avr) genes or haplotypes in the differential isolates. Four R genes, i.e., Pi2, Piz-t, Pi50, and Pikm were found to contribute significantly to the resistance of the elite rice varieties. Other two genes, Pi9 and Pikh, which were not utilized in rice varieties, showed promising values in breeding durable resistance due to their high resistance frequencies to the contemporary rice blast population. The sequence-based molecular diagnosis provided a promising approach for the identification and verification of haplotypes in different R-gene loci and effective R genes valuable for breeding durable rice resistance to rice blast.

Complex Admixture Preceded and Followed the Extinction of Wisent in the Wild

Retracing complex population processes that precede extreme bottlenecks may be impossible using data from living individuals. The wisent (Bison bonasus), Europe’s largest terrestrial mammal, exemplifies such a population history, having gone extinct in the wild but subsequently restored by captive breeding efforts. Using low coverage genomic data from modern and historical individuals, we investigate population processes occurring before and after this extinction. Analysis of aligned genomes supports the division of wisent into two previously recognized subspecies, but almost half of the genomic alignment contradicts this population history as a result of incomplete lineage sorting and admixture. Admixture between subspecies populations occurred prior to extinction and subsequently during the captive breeding program. Admixture with the Bos cattle lineage is also widespread but results from ancient events rather than recent hybridization with domestics. Our study demonstrates the huge potential of historical genomes for both studying evolutionary histories and for guiding conservation strategies.

A comprehensive benchmarking study of protocols and sequencing platforms for 16S rRNA community profiling

Background

In the last 5 years, the rapid pace of innovations and improvements in sequencing technologies has completely changed the landscape of metagenomic and metagenetic experiments. Therefore, it is critical to benchmark the various methodologies for interrogating the composition of microbial communities, so that we can assess their strengths and limitations. The most common phylogenetic marker for microbial community diversity studies is the 16S ribosomal RNA gene and in the last 10 years the field has moved from sequencing a small number of amplicons and samples to more complex studies where thousands of samples and multiple different gene regions are interrogated.

Results

We assembled 2 synthetic communities with an even (EM) and uneven (UM) distribution of archaeal and bacterial strains and species, as metagenomic control material, to assess performance of different experimental strategies. The 2 synthetic communities were used in this study, to highlight the limitations and the advantages of the leading sequencing platforms: MiSeq (Illumina), The Pacific Biosciences RSII, 454 GS-FLX/+ (Roche), and IonTorrent (Life Technologies). We describe an extensive survey based on synthetic communities using 3 experimental designs (fusion primers, universal tailed tag, ligated adaptors) across the 9 hypervariable 16S rDNA regions. We demonstrate that library preparation methodology can affect data interpretation due to different error and chimera rates generated during the procedure. The observed community composition was always biased, to a degree that depended on the platform, sequenced region and primer choice. However, crucially, our analysis suggests that 16S rRNA sequencing is still quantitative, in that relative changes in abundance of taxa between samples can be recovered, despite these biases.

Conclusion

We have assessed a range of experimental conditions across several next generation sequencing platforms using the most up-to-date configurations. We propose that the choice of sequencing platform and experimental design needs to be taken into consideration in the early stage of a project by running a small trial consisting of several hypervariable regions to quantify the discriminatory power of each region. We also suggest that the use of a synthetic community as a positive control would be beneficial to identify the potential biases and procedural drawbacks that may lead to data misinterpretation. The results of this study will serve as a guideline for making decisions on which experimental condition and sequencing platform to consider to achieve the best microbial profiling.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-2194-9) contains supplementary material, which is available to authorized users.

Genome-wide association study for posthitis in the free-living population of European bison (Bison bonasus)

Background

About 5–6% of the European bison (Bison bonasus) males are affected by posthitis (necrotic inflammation of the prepuce) and die in the wild forest. Despite many years of study, pathogenesis of this disease has not yet been determined. The main aim of the study was to find SNP markers significantly associated with the incidence of posthitis and mine the genome for candidate genes potentially involved in the development of the disease.

Results

It was shown that relatively small number of SNPs effects reached genome-wide significance after false discovery rate (FDR) correction. Among 25 significant markers, the highest effects were found for two SNPs (rs110456748 and rs136792896) located at the distance of 23846 bp and 37742 bp, respectively, from OR10A3 gene (olfactory receptor genes), known to be involved in atopic dermatitis in humans. It was also observed that five other significant SNP markers were located in the proximity of candidate genes involved in severe diseases of skin tissue and cancer/tumour development of epithelial or testicular germ cells, which suggest their potential participation in the posthitis. The 25 investigated SNPs showed marked differences in allelic and genotypic frequencies between the healthy and affected bison groups.

Conclusions

The 2 Mb region of the BTA15 chromosome is involved in genetic background of posthitis and should be closer examined to find causal mutations helpful in better understanding of the disease ethology and to control its incidence in the future.

Reviewers

This article was reviewed by Prof. Lev Klebanov and Dr. Fyodor Kondrashov.

Targeted amplicon sequencing (TAS): a scalable next-gen approach to multilocus, multitaxa phylogenetics

Next-gen sequencing technologies have revolutionized data collection in genetic studies and advanced genome biology to novel frontiers. However, to date, next-gen technologies have been used principally for whole genome sequencing and transcriptome sequencing. Yet many questions in population genetics and systematics rely on sequencing specific genes of known function or diversity levels. Here, we describe a targeted amplicon sequencing (TAS) approach capitalizing on next-gen capacity to sequence large numbers of targeted gene regions from a large number of samples. Our TAS approach is easily scalable, simple in execution, neither time-nor labor-intensive, relatively inexpensive, and can be applied to a broad diversity of organisms and/or genes. Our TAS approach includes a bioinformatic application, BarcodeCrucher, to take raw next-gen sequence reads and perform quality control checks and convert the data into FASTA format organized by gene and sample, ready for phylogenetic analyses. We demonstrate our approach by sequencing targeted genes of known phylogenetic utility to estimate a phylogeny for the Pancrustacea. We generated data from 44 taxa using 68 different 10-bp multiplexing identifiers. The overall quality of data produced was robust and was informative for phylogeny estimation. The potential for this method to produce copious amounts of data from a single 454 plate (e.g., 325 taxa for 24 loci) significantly reduces sequencing expenses incurred from traditional Sanger sequencing. We further discuss the advantages and disadvantages of this method, while offering suggestions to enhance the approach.

A global view of 54,001 single nucleotide polymorphisms (SNPs) on the Illumina BovineSNP50 BeadChip and their transferability to water buffalo

The Illumina BovineSNP50 BeadChip features 54,001 informative single nucleotide polymorphisms (SNPs) that uniformly span the entire bovine genome. Among them, 52,255 SNPs have locations assigned in the current genome assembly (Btau_4.0), including 19,294 (37%) intragenic SNPs (i.e., located within genes) and 32,961 (63%) intergenic SNPs (i.e., located between genes). While the SNPs represented on the Illumina Bovine50K BeadChip are evenly distributed along each bovine chromosome, there are over 14,000 genes that have no SNPs placed on the current BeadChip. Kernel density estimation, a non-parametric method, was used in the present study to identify SNP-poor and SNP-rich regions on each bovine chromosome. With bandwidth = 0.05 Mb, we observed that most regions have SNP densities within 2 standard deviations of the chromosome SNP density mean. The SNP density on chromosome X was the most dynamic, with more than 30 SNP-rich regions and at least 20 regions with no SNPs. Genotyping ten water buffalo using the Illumina BovineSNP50 BeadChip revealed that 41,870 of the 54,001 SNPs are fully scored on all ten water buffalo, but 6,771 SNPs are partially scored on one to nine animals. Both fully scored and partially/no scored SNPs are clearly clustered with various sizes on each chromosome. However, among 43,687 bovine SNPs that were successfully genotyped on nine and ten water buffalo, only 1,159 were polymorphic in the species. These results indicate that the SNPs sites, but not the polymorphisms, are conserved between two species. Overall, our present study provides a solid foundation to further characterize the SNP evolutionary process, thus improving understanding of within- and between-species biodiversity, phylogenetics and adaption to environmental changes.