Identification of copy number variations in three Chinese horse breeds using 70K single nucleotide polymorphism BeadChip array

Copy number variation (CNV), an essential form of genetic variation, has been increasingly recognized as one promising genetic marker in the analysis of animal genomes. Here, we used the Equine 70K single nucleotide polymorphism genotyping array for the genome‐wide detection of CNVs in 96 horses from three diverse Chinese breeds: Debao pony (DB), Mongolian horse (MG) and Yili horse (YL). A total of 287 CNVs were determined and merged into 122 CNV regions (CNVRs) ranging from 199 bp to 2344 kb in size and distributed in a heterogeneous manner on chromosomes. These CNVRs were integrated with seven existing reports to generate a composite genome‐wide dataset of 1558 equine CNVRs, revealing 69 (56.6%) novel CNVRs. The majority (69.7%) of the 122 CNVRs overlapped with 438 genes, whereas 30.3% were located in intergenic regions. Most of these genes were associated with common CNVRs, which were shared by divergent horse breeds. As many as 60, 42 and 91 genes overlapping with the breed‐specific ss were identified in DB, MG and YL respectively. Among these genes, FGF11, SPEM1, PPARG, CIDEB, HIVEP1 and GALR may have potential relevance to breed‐specific traits. These findings provide valuable information for understanding the equine genome and facilitating association studies of economically important traits with equine CNVRs in the future.     

拷贝数变异的全基因组关联分析

基因组拷贝数变异(copy number variations,CNVs)是指与基因组参考序列相比,基因组中≥1 kb的DNA片段插入、缺失和/或扩增,及其互相组合衍生出的复杂变异．由于其具有分布范围广、可遗传、相对稳定和高度异质性等特点,目前认为,CNVs是一种新的可以作为疾病易感标志的基因组DNA多态性,其变异引起的基因剂量改变可以导致表型改变．最近,一种基于CNVs的新的疾病易感基因鉴定策略——CNV全基因组关联分析开始出现,这一策略和传统的基于单核苷酸多态性的关联分析具有互补性,通过认识基因组结构变异可以认识复杂疾病的分子机制和遗传基础．     

Characterization of eight single nucleotide polymorphism markers in Aedes aegypti

Sequencing of part of seven genes from Aedes aegypti collected in 16 Brazilian cities revealed the existence of 53 single nucleotide polymorphisms (SNPs), representing one SNP every 52 base pairs. From these 53 SNPs, we selected eight that are independent and highly polymorphic. We describe the use of these markers for differentiation of Brazilian populations of A. aegypti. These are the first SNPs developed for delineating population structure in A. aegypti, and will be a useful complement to epidemiological studies.     

Characterization of 15 single nucleotide polymorphism markers for chimpanzees (Pan troglodytes)

We report the characterization of 15 new single nucleotide polymorphism markers for a threatened species, the chimpanzee (Pan troglodytes), developed using a targeted gene approach. These markers are derived from the Y chromosome and autosomal regions of the genome and show frequency differences between chimpanzee subspecies from central and western Africa. These single nucleotide polymorphism markers are the first to be designed for the genotyping of wild chimpanzee populations and will provide a useful addition to the genetic tools employed for the conservation management of this threatened species.     

Haplotype structure and copy number polymorphism of the beta‐defensin 7 genes in diverse chicken breeds

Beta‐defensins is a family of avian peptides related to the innate immune system. Copy number variation was recently reported for the avian beta‐defensin 7 gene (AvBD7) between the highly inbred Leghorn and Fayoumi lines. Here, we examined copy number variants in 35 different chicken breeds and found that 31 of them have at least the same representation of the duplicated AvBD7 allele. We also found haplotypes upstream of the AvBD6 regions that are strongly linked to the AvBD7 duplication. We observed a strong linkage disequilibrium spanning of the upstream region of the AvBD6 gene, with two SNPs being flanking markers to detect duplication of the AvBD7.     

A genome‐wide scan for copy number variations using high‐density single nucleotide polymorphism array in Simmental cattle

Copy number variations (CNVs) have recently been identified as promising sources of genetic variation, complementary to single nucleotide polymorphisms (SNPs). As a result, detection of CNVs has attracted a great deal of attention. In this study, we performed genome‐wide CNV detection using Illumina Bovine HD BeadChip (770k) data on 792 Simmental cattle. A total of 263 CNV regions (CNVRs) were identified, which included 137 losses, 102 gains and 24 regions classified as both loss and gain, covering 35.48 Mb (1.41%) of the bovine genome. The length of these CNVRs ranged from 10.18 kb to 1.76 Mb, with an average length of 134.78 kb and a median length of 61.95 kb. In 136 of these regions, a total of 313 genes were identified related to biological functions such as transmembrane activity and olfactory transduction activity. To validate the results, we performed quantitative PCR to detect nine randomly selected CNVRs and successfully confirmed seven (77.6%) of them. Our results present a map of cattle CNVs derived from high‐density SNP data, which expands the current CNV map of the cattle genome and provides useful information for investigation of genomic structural variation in cattle.     

Analysis of public single nucleotide polymorphisms in commercial pig populations

More than 5500 pig single nucleotide polymorphisms (SNPs) were recently identified and deposited in the public domain. To test the usefulness of these public SNPs, 109 SNPs were analysed for polymorphism within six commercial pig populations. A functional polymerase chain reaction (PCR) assay was obtained for 103 SNPs and it was possible to validate c. 59% by PCR-restriction fragment length polymorphism. Furthermore, polymorphism was found using a relatively limited number of genomic DNA samples, indicating that these polymorphisms are segregating at a useful frequency in these populations. The high percentage of validated markers demonstrates the utility of these public pig SNPs to identify loci responsible for economically important traits in commercial pig populations.     

Using blocks of linked single nucleotide polymorphisms as highly polymorphic genetic markers for parentage analysis

Single nucleotide polymorphisms (SNPs) are plentiful in most genomes and amenable to high throughput genotyping, but they are not yet popular for parentage or paternity analysis. The markers are bi-allelic, so individually they contain little information about parentage, and in nonmodel organisms the process of identifying large numbers of unlinked SNPs can be daunting. We explore the possibility of using blocks of between three and 26 linked SNPs as highly polymorphic molecular markers for reconstructing male genotypes in polyandrous organisms with moderate (five offspring) to large (25 offspring) clutches of offspring. Haplotypes are inferred for each block of linked SNPs using the programs Haplore and Phase 2.1. Each multi-SNP haplotype is then treated as a separate allele, producing a highly polymorphic, 'microsatellite-like' marker. A simulation study is performed using haplotype frequencies derived from empirical data sets from Drosophila melanogaster and Mus musculus populations. We find that the markers produced are competitive with microsatellite loci in terms of single parent exclusion probabilities, particularly when using six or more linked SNPs to form a haplotype. These markers contain only modest rates of missing data and genotyping or phasing errors and thus should be seriously considered as molecular markers for parentage analysis, particularly when the study is interested in the functional significance of polymorphisms across the genome.     

Evaluation of 16 loci to examine the cross‐species utility of single nucleotide polymorphism arrays

Large collections of single nucleotide polymorphisms (SNPs) have recently been identified from a number of livestock genomes. This raises the possibility that SNP arrays might be useful for analysis in related species for which few genetic markers are currently available. To address the likely success of such an approach, the aim of this study was to examine the threshold number and position of flanking mutations which act to prevent genotype calls being produced. Sequence diversity was measured across 16 loci containing SNPs known either to work successfully between species or fail between species. In pairwise comparisons between domestic and wild sheep, sequence divergence surrounding working SNP assays was significantly lower than that surrounding non‐functional assays. In addition, the location of flanking mismatches tended to be closer to the target SNP in loci that failed to generate genotype calls across species. The magnitude of sequence divergence observed for both working and non‐functional assays was compared with the divergence separating domestic sheep from European Mouflon, African Barbary, goat and cattle. The results suggest that the utility of SNP arrays for analysis of shared polymorphism will be restricted to closely related pairs of species. Analysis across more divergent species will, however, be successful for other objectives, such as the identification of the ancestral state of SNPs.     

Estimation of the number of SNP genetic markers required for parentage verification

Formulae were developed to compute exclusion probabilities for parentage confirmation for any number of diallelic markers under the assumption that the minor allele frequency (MAF) varied among markers, but has a uniform distribution. Three scenarios were analysed: a progeny with (1) a single putative parent; (2) two putative parents; and (3) one actual parent and one putative parent. Exclusion probabilities were computed for minimum values for the MAFs of 0.1, 0.2 and 0.3, and required either one or at least two conflicts for exclusion. The numbers of markers required to obtain 99% exclusion probabilities based on a single conflict for the three minimum MAFs were 54, 45 and 39 for scenario 1; 17, 16 and 15 for scenario 2; and 28, 25 and 24 for scenario 3. The requirement of at least two conflicts for exclusion increased the number of markers required by approximately 45% for all three scenarios and all three minimum MAFs. The results obtained by the analytical formulae were very close to results obtained by simulation and to values in the literature for specific marker sets.     

Genomic copy number variation in Mus musculus

Background

Copy number variation is an important dimension of genetic diversity and has implications in development and disease. As an important model organism, the mouse is a prime candidate for copy number variant (CNV) characterization, but this has yet to be completed for a large sample size. Here we report CNV analysis of publicly available, high-density microarray data files for 351 mouse tail samples, including 290 mice that had not been characterized for CNVs previously.

Results

We found 9634 putative autosomal CNVs across the samples affecting 6.87 % of the mouse reference genome. We find significant differences in the degree of CNV uniqueness (single sample occurrence) and the nature of CNV-gene overlap between wild-caught mice and classical laboratory strains. CNV-gene overlap was associated with lipid metabolism, pheromone response and olfaction compared to immunity, carbohydrate metabolism and amino-acid metabolism for wild-caught mice and classical laboratory strains, respectively. Using two subspecies of wild-caught Mus musculus, we identified putative CNVs unique to those subspecies and show this diversity is better captured by wild-derived laboratory strains than by the classical laboratory strains. A total of 9 genic copy number variable regions (CNVRs) were selected for experimental confirmation by droplet digital PCR (ddPCR).

Conclusion

The analysis we present is a comprehensive, genome-wide analysis of CNVs in Mus musculus, which increases the number of known variants in the species and will accelerate the identification of novel variants in future studies.

Electronic supplementary material

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

Delineation of the frequency and boundary of chromosomal copy number variations in paediatric neuroblastoma

Neuroblastoma, the most common solid tumour in early childhood, is characterized by very frequent chromosomal copy number variations (CNVs). While chromosome 2p amplification, 17q gain, 1p and 11q deletion in human neuroblastoma tissues are well-known, the exact frequencies and boundaries of the chromosomal CNVs have not been delineated. We analysed the publicly available single nucleotide polymorphism (SNP) array data which were originally generated by the Therapeutically Applicable Research to Generate Effective Treatments (TARGET) initiative, defined the frequencies and boundaries of chromosomes 2p11.2 – 2p25.3 amplification, 17q11.1-17q25.3 gain, 1p13.3-1p36.33 deletion and 11q13.3-11q25 deletion in neuroblastoma tissues, and identified chromosome 7q14.1 (Chr7:38254795-38346971) and chromosome 14q11.2 (Chr14:21637401-22024617) deletion in blood and bone marrow samples from neuroblastoma patients, but not in tumour tissues. Kaplan Meier analysis showed that double deletion of Chr7q14.1 and Chr14q11.2 correlated with poor prognosis in MYCN gene amplified neuroblastoma patients. In conclusion, the oncogenes amplified or gained and tumour suppressor genes deleted within the boundaries of chromosomal CNVs in tumour tissues should be studied for their roles in tumourigenesis and as therapeutic targets. Focal deletions of Chr7q14.1 and Chr14q11.2 together in blood and bone marrow samples from neuroblastoma patients can be used as a marker for poorer prognosis and more aggressive therapies.     

Six diagnostic single nucleotide polymorphism markers for detecting introgression between cutthroat and rainbow trouts

Ten primer pairs were screened to develop single nucleotide polymorphism (SNP) TaqMan assays that will distinguish California golden trout and some rainbow trouts (Oncorhynchus mykiss sspp., O. m. aguabonita) from the Paiute and Lahontan cutthroat trouts (Oncorhynchus clarkii seleniris, O. c. henshawi). From these 10 primer pairs, one mitochondrial and five nuclear fixed SNP differences were discovered and developed into TaqMan assays. These six assays will be useful for characterizing and monitoring hybridization between these groups. Additional Oncorhynchus clarkii sspp. and Oncorhynchus mykiss sspp. were assayed to determine if these assays are useful in closely related species.     

Single nucleotide polymorphism discovery in barley using autoSNPdb

Molecular markers are used to provide the link between genotype and phenotype, for the production of molecular genetic maps and to assess genetic diversity within and between related species. Single nucleotide polymorphisms (SNPs) are the most abundant molecular genetic marker. SNPs can be identified in silico , but care must be taken to ensure that the identified SNPs reflect true genetic variation and are not a result of errors associated with DNA sequencing. The SNP detection method autoSNP has been developed to identify SNPs from sequence data for any species. Confidence in the predicted SNPs is based on sequence redundancy, and haplotype co-segregation scores are calculated for a further independent measure of confidence. We have extended the autoSNP method to produce autoSNPdb, which integrates SNP and gene annotation information with a graphical viewer. We have applied this software to public barley expressed sequences, and the resulting database is available over the Internet. SNPs can be viewed and searched by sequence, functional annotation or predicted synteny with a reference genome, in this case rice. The correlation between SNPs and barley cultivar, expressed tissue type and development stage has been collated for ease of exploration. An average of one SNP per 240 bp was identified, with SNPs more prevalent in the 5' regions and simple sequence repeat (SSR) flanking sequences. Overall, autoSNPdb can provide a wealth of genetic polymorphism information for any species for which sequence data are available.     

Bioinformatics approaches and resources for single nucleotide polymorphism functional analysis

Since the initial sequencing of the human genome, many projects are underway to understand the effects of genetic variation between individuals. Predicting and understanding the downstream effects of genetic variation using computational methods are becoming increasingly important for single nucleotide polymorphism (SNP) selection in genetics studies and understanding the molecular basis of disease. According to the NIH, there are now more than four million validated SNPs in the human genome. The volume of known genetic variations lends itself well to an informatics approach. Bioinformaticians have become very good at functional inference methods derived from functional and structural genomics. This review will present a broad overview of the tools and resources available to collect and understand functional variation from the perspective of structure, expression, evolution and phenotype. Additionally, public resources available for SNP identification and characterisation are summarised.     

The development and characterization of a 57K single nucleotide polymorphism array for rainbow trout

In this study, we describe the development and characterization of the first high‐density single nucleotide polymorphism (SNP) genotyping array for rainbow trout. The SNP array is publically available from a commercial vendor (Affymetrix). The SNP genotyping quality was high, and validation rate was close to 90%. This is comparable to other farm animals and is much higher than previous smaller scale SNP validation studies in rainbow trout. High quality and integrity of the genotypes are evident from sample reproducibility and from nearly 100% agreement in genotyping results from other methods. The array is very useful for rainbow trout aquaculture populations with more than 40 900 polymorphic markers per population. For wild populations that were confounded by a smaller sample size, the number of polymorphic markers was between 10 577 and 24 330. Comparison between genotypes from individual populations suggests good potential for identifying candidate markers for populations' traceability. Linkage analysis and mapping of the SNPs to the reference genome assembly provide strong evidence for a wide distribution throughout the genome with good representation in all 29 chromosomes. A total of 68% of the genome scaffolds and contigs were anchored through linkage analysis using the SNP array genotypes, including ~20% of the genome assembly that has not been previously anchored to chromosomes.     

Development and application of genotyping technologies

With the completion of Human Genome Project,International HapMap Project and the publication of copy number variation in human genome,a great number of accurate,rapid,and cost-effective technologies for SNP analysis have been developed,promoting the research of the complex diseases.This article presents a review of widely used genotyping techniques,and the progress and prospect in the study of complex diseases in terms of the projects and achievements of Chinese National Human Genome Center at Shanghai(CHGC...     

Characterization of single nucleotide polymorphism markers for eelgrass (Zostera marina)

We characterized 37 single nucleotide polymorphism (SNP) makers for eelgrass Zostera marina. SNP markers were developed using existing EST (expressed sequence tag)-libraries to locate polymorphic loci and develop primers from the functional expressed genes that are deposited in The ZOSTERA database (V1.2.1). SNP loci were genotyped using a single-base-extension approach which facilitated high-throughput genotyping with minimal optimization time. These markers show a wide range of variability among 25 eelgrass populations and will be useful for population genetic studies including evaluation of population structure, historical demography, and phylogeography. Potential applications include haplotype inference of physically linked SNPs and identification of genes under selection for temperature and desiccation stress.     

Efficiency of selection, as measured by single nucleotide polymorphism variation, is dependent on inbreeding rate in Drosophila melanogaster

It is often hypothesized that slow inbreeding causes less inbreeding depression than fast inbreeding at the same absolute level of inbreeding. Possible explanations for this phenomenon include the more efficient purging of deleterious alleles and more efficient selection for heterozygote individuals during slow, when compared with fast, inbreeding. We studied the impact of inbreeding rate on the loss of heterozygosity and on morphological traits in Drosophila melanogaster. We analysed five noninbred control lines, 10 fast inbred lines and 10 slow inbred lines; the inbred lines all had an expected inbreeding coefficient of approximately 0.25. Forty single nucleotide polymorphisms in DNA coding regions were genotyped, and we measured the size and shape of wings and counted the number of sternopleural bristles on the genotyped individuals. We found a significantly higher level of genetic variation in the slow inbred lines than in the fast inbred lines. This higher genetic variation was resulting from a large contribution from a few loci and a smaller effect from several loci. We attributed the increased heterozygosity in the slow inbred lines to the favouring of heterozygous individuals over homozygous individuals by natural selection, either by associative over‐dominance or balancing selection, or a combination of both. Furthermore, we found a significant polynomial correlation between genetic variance and wing size and shape in the fast inbred lines. This was caused by a greater number of homozygous individuals among the fast inbred lines with small, narrow wings, which indicated inbreeding depression. Our results demonstrated that the same amount of inbreeding can have different effects on genetic variance depending on the inbreeding rate, with slow inbreeding leading to higher genetic variance than fast inbreeding. These results increase our understanding of the genetic basis of the common observation that slow inbred lines express less inbreeding depression than fast inbred lines. In addition, this has more general implications for the importance of selection in maintaining genetic variation.     

Genome‐wide copy number profiling using high‐density SNP array in chickens

Phenotypic diversity is a direct consequence resulting mainly from the impact of underlying genetic variation, and recent studies have shown that copy number variation (CNV) is emerging as an important contributor to both phenotypic variability and disease susceptibility. Herein, we performed a genome‐wide CNV scan in 96 chickens from 12 diversified breeds, benefiting from the high‐density Affymetrix 600 K SNP arrays. We identified a total of 231 autosomal CNV regions (CNVRs) encompassing 5.41 Mb of the chicken genome and corresponding to 0.59% of the autosomal sequence. The length of these CNVRs ranged from 2.6 to 586.2 kb with an average of 23.4 kb, including 130 gain, 93 loss and eight both gain and loss events. These CNVRs, especially deletions, had lower GC content and were located particularly in gene deserts. In particular, 102 CNVRs harbored 128 chicken genes, most of which were enriched in immune responses. We obtained 221 autosomal CNVRs after converting probe coordinates to Galgal3, and comparative analysis with previous studies illustrated that 153 of these CNVRs were regarded as novel events. Furthermore, qPCR assays were designed for 11 novel CNVRs, and eight (72.73%) were validated successfully. In this study, we demonstrated that the high‐density 600 K SNP array can capture CNVs with higher efficiency and accuracy and highlighted the necessity of integrating multiple technologies and algorithms. Our findings provide a pioneering exploration of chicken CNVs based on a high‐density SNP array, which contributes to a more comprehensive understanding of genetic variation in the chicken genome and is beneficial to unearthing potential CNVs underlying important traits of chickens.