Genomic organization of Tropomodulins 2 and 4 and unusual intergenic and intraexonic splicing of YL-1 and Tropomodulin 4

Background

Single nucleotide polymorphisms (SNPs) are the foundation of powerful complex trait and pharmacogenomic analyses. The availability of large SNP databases, however, has emphasized a need for inexpensive SNP genotyping methods of commensurate simpliCity, robustness, and scalability. We describe a solution-based, microtiter plate method for SNP genotyping of human genomic DNA. The method is based upon allele discrimination by ligation of open circle probes followed by rolling circle amplification of the signal using fluorescent primers. Only the probe with a 3' base complementary to the SNP is circularized by ligation.

Results

SNP scoring by ligation was optimized to a 100,000 fold discrimination against probe mismatched to the SNP. The assay was used to genotype 10 SNPs from a set of 192 genomic DNA samples in a high-throughput format. Assay directly from genomic DNA eliminates the need to preamplify the target as done for many other genotyping methods. The sensitivity of the assay was demonstrated by genotyping from 1 ng of genomic DNA. We demonstrate that the assay can detect a single molecule of the circularized probe.

Conclusions

Compatibility with homogeneous formats and the ability to assay small amounts of genomic DNA meets the exacting requirements of automated, high-throughput SNP scoring.     

GnRHR gene polymorphisms and their effects on reproductive performance in Chinese goats

In this study, polymorphisms in the goat GnRHR gene exon 1 were detected by PCR-SSCP and DNA sequencing methods in 786 individuals from two different goat breeds. Two haplotypes (A and B), two observed genotypes (AA and AB), and two single nucleotide polymorphisms (SNPs) were detected, which resulted in five amino acid substitutions. The frequencies of haplotypes A and B in the two goat breeds were 0.78–0.83 and 0.17–0.22, respectively. The SNP locus was in Hardy–Weinberg disequilibrium in the two goat breeds (P < 0.05). Polymorphisms of the GnRHR gene were shown to be associated with litter size in the two goat breeds. The SNPs in the goat GnRHR gene had significant effects on litter size (P < 0.05). Therefore, these results suggest that the GnRHR gene is a strong candidate gene that affects litter size in goat.     

SNP discovery in nonmodel organisms: strand bias and base‐substitution errors reduce conversion rates

Single nucleotide polymorphisms (SNPs) have become the marker of choice for genetic studies in organisms of conservation, commercial or biological interest. Most SNP discovery projects in nonmodel organisms apply a strategy for identifying putative SNPs based on filtering rules that account for random sequencing errors. Here, we analyse data used to develop 4723 novel SNPs for the commercially important deep‐sea fish, orange roughy (Hoplostethus atlanticus), to assess the impact of not accounting for systematic sequencing errors when filtering identified polymorphisms when discovering SNPs. We used SAMtools to identify polymorphisms in a velvet assembly of genomic DNA sequence data from seven individuals. The resulting set of polymorphisms were filtered to minimize ‘bycatch’—polymorphisms caused by sequencing or assembly error. An Illumina Infinium SNP chip was used to genotype a final set of 7714 polymorphisms across 1734 individuals. Five predictors were examined for their effect on the probability of obtaining an assayable SNP: depth of coverage, number of reads that support a variant, polymorphism type (e.g. A/C), strand‐bias and Illumina SNP probe design score. Our results indicate that filtering out systematic sequencing errors could substantially improve the efficiency of SNP discovery. We show that BLASTX can be used as an efficient tool to identify single‐copy genomic regions in the absence of a reference genome. The results have implications for research aiming to identify assayable SNPs and build SNP genotyping assays for nonmodel organisms.     

SNP discovery using Paired‐End RAD‐tag sequencing on pooled genomic DNA of Sisymbrium austriacum (Brassicaceae)

Single nucleotide polymorphisms SNPs are rapidly replacing anonymous markers in population genomic studies, but their use in non model organisms is hampered by the scarcity of cost‐effective approaches to uncover genome‐wide variation in a comprehensive subset of individuals. The screening of one or only a few individuals induces ascertainment bias. To discover SNPs for a population genomic study of the Pyrenean rocket (Sisymbrium austriacum subsp. chrysanthum), we undertook a pooled RAD‐PE (Restriction site Associated DNA Paired‐End sequencing) approach. RAD tags were generated from the PstI‐digested pooled genomic DNA of 12 individuals sampled across the species distribution range and paired‐end sequenced using Illumina technology to produce ~24.5 Mb of sequences, covering ~7% of the specie's genome. Sequences were assembled into ~76 000 contigs with a mean length of 323 bp (N₅₀ = 357 bp, sequencing depth = 24x). In all, >15 000 SNPs were called, of which 47% were annotated in putative genic regions based on homology with the Arabidopsis thaliana genome. Gene ontology (GO) slim categorization demonstrated that the identified SNPs covered extant genic variation well. The validation of 300 SNPs on a larger set of individuals using a KASPar assay underpinned the utility of pooled RAD‐PE as an inexpensive genome‐wide SNP discovery technique (success rate: 87%). In addition to SNPs, we discovered >600 putative SSR markers.     

Development and evaluation of single-nucleotide polymorphism markers in allotetraploid rapeseed (<Emphasis Type="Italic">Brassica napus</Emphasis> L.)

Single-nucleotide polymorphisms (SNPs) and insertion–deletions (INDELs) are currently the important classes of genetic markers for major crop species. In this study, methods for developing SNP markers in rapeseed (Brassica napus L.) and their in silico mapping and use for genotyping are demonstrated. For the development of SNP and INDEL markers, 181 fragments from 121 different gene sequences spanning 86 kb were examined. A combination of different selection methods (genome-specific amplification, hetero-duplex analysis and sequence analysis) allowed the detection of 18 singular fragments that showed a total of 87 SNPs and 6 INDELs between 6 different rapeseed varieties. The average frequency of sequence polymorphism was estimated to be one SNP every 247 bp and one INDEL every 3,583 bp. Most SNPs and INDELs were found in non-coding regions. Polymorphism information content values for SNP markers ranged between 0.02 and 0.50 in a set of 86 varieties. Using comparative genetics data for B. napus and Arabidopsis thaliana, an allocation of SNP markers to linkage groups in rapeseed was achieved: a unique location was determined for seven gene sequences; two and three possible locations were found for six and four sequences, respectively. The results demonstrate the usefulness of existing genomic resources for SNP discovery in rapeseed.     

Genomewide single nucleotide polymorphism discovery in Atlantic salmon (Salmo salar): validation in wild and farmed American and European populations

A considerable number of single nucleotide polymorphisms (SNPs) are required to elucidate genotype–phenotype associations and determine the molecular basis of important traits. In this work, we carried out de novo SNP discovery accounting for both genome duplication and genetic variation from American and European salmon populations. A total of 9 736 473 nonredundant SNPs were identified across a set of 20 fish by whole‐genome sequencing. After applying six bioinformatic filtering steps, 200 K SNPs were selected to develop an Affymetrix Axiom^® myDesign Custom Array. This array was used to genotype 480 fish representing wild and farmed salmon from Europe, North America and Chile. A total of 159 099 (79.6%) SNPs were validated as high quality based on clustering properties. A total of 151 509 validated SNPs showed a unique position in the genome. When comparing these SNPs against 238 572 markers currently available in two other Atlantic salmon arrays, only 4.6% of the SNP overlapped with the panel developed in this study. This novel high‐density SNP panel will be very useful for the dissection of economically and ecologically relevant traits, enhancing breeding programmes through genomic selection as well as supporting genetic studies in both wild and farmed populations of Atlantic salmon using high‐resolution genomewide information.     

Single‐nucleotide polymorphism discovery and diversity in the model legume Medicago truncatula

Extensive genomic resources are available in the model legume Medicago truncatula. Here, we present the discovery and design of the first array of single‐nucleotide polymorphism (SNP) markers in M. truncatula through large‐scale Sanger resequencing of genomic fragments spanning the genome, in a diverse panel of 16 M. truncatula accessions. Both anonymous fragments and fragments targeting candidate genes for flowering phenology and symbiosis were surveyed for nucleotide variation in almost 230 kb of unique genomic regions. A set of 384 SNP markers was designed for an Illumina's GoldenGate assay, genotyped on a collection of 192 inbred lines (CC192) representing the geographical range of the species and used to survey the diversity of two natural populations. Finally, 86% of the tested SNPs were of high quality and exhibited polymorphism in the CC192 collection. Even at the population level, we detected polymorphism for more than 50% of the selected SNPs. Analysis of the allele frequency spectrum in the CC192 showed a reduced ascertainment bias, mostly limited to very rare alleles (frequency <0.01). The substantial polymorphism detected at the species and population levels, the high marker quality and the potential to survey large samples of individuals make this set of SNP markers a valuable tool to improve our understanding of the effect of demographic and selective factors that shape the natural genetic diversity within the selfing species Medicago truncatula.     

Analysis of gene-derived SNP marker polymorphism in US wheat (Triticum aestivum L.) cultivars

In this study, we developed 359 detection primers for single nucleotide polymorphisms (SNPs) previously discovered within intron sequences of wheat genes and used them to evaluate SNP polymorphism in common wheat (Triticum aestivum L.). These SNPs showed an average polymorphism information content (PIC) of 0.18 among 20 US elite wheat cultivars, representing seven market classes. This value increased to 0.23 when SNPs were pre-selected for polymorphisms among a diverse set of 13 hexaploid wheat accessions (excluding synthetic wheats) used in the wheat SNP discovery project (). PIC values for SNP markers in the D genome were approximately half of those for the A and B genomes. D genome SNPs also showed a larger PIC reduction relative to the other genomes (P < 0.05) when US cultivars were compared with the more diverse set of 13 wheat accessions. Within those accessions, D genome SNPs show a higher proportion of alleles with low minor allele frequencies (<0.125) than found in the other two genomes. These data suggest that the reduction of PIC values in the D genome was caused by differential loss of low frequency alleles during the population size bottleneck that accompanied the development of modern commercial cultivars. Additional SNP discovery efforts targeted to the D genome in elite wheat germplasm will likely be required to offset the lower diversity of this genome. With increasing SNP discovery projects and the development of high-throughput SNP assay technologies, it is anticipated that SNP markers will play an increasingly important role in wheat genetics and breeding applications. Electronic supplementary material  The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

Polymorphisms in the promoter region of ESR2 gene and breast cancer susceptibility

Genetic variations like single nucleotide polymorphisms (SNPs) in genes involved in estrogen biosynthesis, metabolism and signal transduction have been suggested to affect breast cancer susceptibility. In this study we tested the hypothesis that polymorphisms in the promoter of ESR2 gene may be associated with increased risk for breast cancer. We analyzed three SNPs in the promoter region of human ESR2 gene by means of allele-specific tetra-primer PCR. A total of 318 sporadic breast cancer cases and 318 age-matched controls were included in the study. With regard to homozygous genotypes, women with sporadic breast cancer more frequently carried the CC genotype of ESR2 promoter SNP rs2987983 (OR 1.99, p = 0.005). Calculation of allele positivity demonstrated that presence of T allele of this SNP was more frequent in healthy women. Our data suggest that a SNP in the promoter region of ESR2 gene might be able to affect breast cancer risk. These results further support the emerging hypothesis that ERβ is an important factor in breast cancer development.     

AffyMAPSDetector: a software tool to characterize Affymetrix GeneChip™ expression arrays with respect to SNPs

Background  

Affymetrix gene expression arrays incorporate paired perfect match (PM) and mismatch (MM) probes to distinguish true signals from those arising from cross-hybridization events. A MM signal often shows greater intensity than a PM signal; we propose that one underlying cause is the presence of allelic variants arising from single nucleotide polymorphisms (SNPs). To annotate and characterize SNP contributions to anomalous probe binding behavior we have developed a software tool called AffyMAPSDetector.     

Evaluation of reduced subsets of single nucleotide polymorphisms for the prediction of age at puberty in sows

Genomic information could be used efficiently to improve traits that are expensive to measure, sex limited or expressed late in life. This study analyzed the phenotypic variation explained by major SNPs and windows for age at puberty in gilts, an indicator of reproductive longevity. A genome‐wide association study using 56 424 SNPs explained 25.2% of the phenotypic variation in age at puberty in a training set (n = 820). All SNPs from the top 10% of 1‐Mb windows explained 33.5% of the phenotypic variance compared to 47.1% explained by the most informative markers (n = 261). In an evaluation population, consisting of subsequent batches (n = 412), the predictive ability of all SNPs from the major 1‐Mb windows was higher compared to the variance captured by the most informative SNP from each of these windows. The phenotypic variance explained in the evaluation population varied from 12.3% to 36.8% when all SNPs from major windows were used compared to 6.5–23.7% explained by most informative SNPs. The correlation between phenotype and genomic prediction values based on SNP effects estimated in the training population was marginal compared to their effects retrained in the evaluation population for all (0.46–0.81) or most informative SNPs (0.30–0.65) from major windows. An increase in genetic gain of 20.5% could be obtained if genomic selection included both sexes compared to females alone. The pleiotropic role of major genes such as AVPR1A could be exploited in selection of both age at puberty and reproductive longevity.     

In Vitro vs In Silico Detected SNPs for the Development of a Genotyping Array: What Can We Learn from a Non-Model Species?

Background

There is considerable interest in the high-throughput discovery and genotyping of single nucleotide polymorphisms (SNPs) to accelerate genetic mapping and enable association studies. This study provides an assessment of EST-derived and resequencing-derived SNP quality in maritime pine (Pinus pinaster Ait.), a conifer characterized by a huge genome size (∼23.8 Gb/C).

Methodology/Principal Findings

A 384-SNPs GoldenGate genotyping array was built from i/ 184 SNPs originally detected in a set of 40 re-sequenced candidate genes (in vitro SNPs), chosen on the basis of functionality scores, presence of neighboring polymorphisms, minor allele frequencies and linkage disequilibrium and ii/ 200 SNPs screened from ESTs (in silico SNPs) selected based on the number of ESTs used for SNP detection, the SNP minor allele frequency and the quality of SNP flanking sequences. The global success rate of the assay was 66.9%, and a conversion rate (considering only polymorphic SNPs) of 51% was achieved. In vitro SNPs showed significantly higher genotyping-success and conversion rates than in silico SNPs (+11.5% and +18.5%, respectively). The reproducibility was 100%, and the genotyping error rate very low (0.54%, dropping down to 0.06% when removing four SNPs showing elevated error rates).

Conclusions/Significance

This study demonstrates that ESTs provide a resource for SNP identification in non-model species, which do not require any additional bench work and little bio-informatics analysis. However, the time and cost benefits of in silico SNPs are counterbalanced by a lower conversion rate than in vitro SNPs. This drawback is acceptable for population-based experiments, but could be dramatic in experiments involving samples from narrow genetic backgrounds. In addition, we showed that both the visual inspection of genotyping clusters and the estimation of a per SNP error rate should help identify markers that are not suitable to the GoldenGate technology in species characterized by a large and complex genome.     

Identification, characterisation and application of single nucleotide polymorphisms for diversity assessment in cassava (Manihot esculenta Crantz)

To monitor genetic diversity in the field it is important that it is measured accurately. Here, we elucidate the potential of single nucleotide polymorphisms (SNPs) for measuring genetic diversity in cassava. The nature and frequency of SNPs was characterised and their utility in genetic diversity assessment compared to that of simple sequence repeats (SSRs). This was achieved by direct sequencing of amplicons in diverse cassava varieties. A total of 26 SNPs were identified from quality sequences of nine genes, giving an estimated frequency of one SNP every 121 nucleotides. Nucleotide diversity ranged from 7.8 × 10⁻⁴ to 5.6 × 10⁻³. Average haplotype-based polymorphic information content (PIC = 0.414) was higher than for individual SNPs (PIC = 0.228). The Mantel test indicated interdependence (r = 0.219; P < 0.001) between SNP and SSR genotypic data. Individual SNPs had lower PIC values than SSRs. For this reason larger numbers of SNPs may be necessary to achieve the same level of discrimination among genotypes provided by SSRs.     

Spatial dynamics and mixing of bluefin tuna in the Atlantic Ocean and Mediterranean Sea revealed using next‐generation sequencing

The Atlantic bluefin tuna is a highly migratory species emblematic of the challenges associated with shared fisheries management. In an effort to resolve the species’ stock dynamics, a genomewide search for spatially informative single nucleotide polymorphisms (SNPs) was undertaken, by way of sequencing reduced representation libraries. An allele frequency approach to SNP discovery was used, combining the data of 555 larvae and young‐of‐the‐year (LYOY) into pools representing major geographical areas and mapping against a newly assembled genomic reference. From a set of 184,895 candidate loci, 384 were selected for validation using 167 LYOY. A highly discriminatory genotyping panel of 95 SNPs was ultimately developed by selecting loci with the most pronounced differences between western Atlantic and Mediterranean Sea LYOY. The panel was evaluated by genotyping a different set of LYOY (n = 326), and from these, 77.8% and 82.1% were correctly assigned to western Atlantic and Mediterranean Sea origins, respectively. The panel revealed temporally persistent differentiation among LYOY from the western Atlantic and Mediterranean Sea (F_ST = 0.008, p = .034). The composition of six mixed feeding aggregations in the Atlantic Ocean and Mediterranean Sea was characterized using genotypes from medium (n = 184) and large (n = 48) adults, applying population assignment and mixture analyses. The results provide evidence of persistent population structuring across broad geographic areas and extensive mixing in the Atlantic Ocean, particularly in the mid‐Atlantic Bight and Gulf of St. Lawrence. The genomic reference and genotyping tools presented here constitute novel resources useful for future research and conservation efforts.     

Effect of polymorphisms in the CAMKMT gene on growth traits in Ujumqin sheep

Our previous genome‐wide association study in sheep revealed that OAR3‐84073899.1 (SNP31) in intron 8 of the CAMKMT gene was significantly associated with post‐weaning gain at the genomic level. Herein, we performed a replication study to investigate single nucleotide polymorphisms (SNPs) within the CAMKMT gene exons, and 1000 bp of the 5′‐ and 3′‐intranslated regions (UTRs) and their associations with growth traits in Ujumqin sheep. Five SNPs were identified through DNA pool sequencing technology: SNP26 in the 5′‐UTR, SNP06 in exon 5, SNP07 in exon 8 and SNP27 and SNP28 in the 3′‐UTR. Six SNPs, including SNP31 in intron 8, were genotyped in the validation group of 343 Ujumqin sheep, and each SNP was classified into three genotypes. The chi‐square test suggested that all the variations were in Hardy–Weinberg equilibrium (P > 0.05) except for SNP28 and SNP31. Linkage disequilibrium analysis showed that SNP07 and SNP31 were strongly linked. An association analysis suggested that SNP06 was significantly associated with chest girth at 6 months of age (P < 0.05). SNP07 exhibited significant correlation with body weight and chest girth at 4 months of age and with body weight, chest girth and chest width at 6 months of age (P < 0.05). SNP27 was highly associated with body weight and chest girth at 4 months of age (P < 0.05), and SNP28 was extremely significantly associated with body weight and chest girth at 4 months of age and with chest girth at 6 months of age (P < 0.01). SNP31 was significantly associated with body weight and shin circumference at 4 months of age and with post‐weaning gain (P < 0.05). Association analysis of the combined effect of SNP07 and SNP31 showed significant correlation with body weight and chest girth at four of months of age (P < 0.05) and with body weight and chest girth at 6 months of age (P < 0.05). These results indicate that the SNPs could be used as meritorious and available genetic markers in growth traits breeding and that the CAMKMT gene may be one of the key candidate genes that affect Ujumqin economic traits.     

Temporal stability and assignment power of adaptively divergent genomic regions between herring (Clupea harengus) seasonal spawning aggregations

Atlantic herring (Clupea harengus), a vital ecosystem component and target of the largest Northwest Atlantic pelagic fishery, undergo seasonal spawning migrations that result in elusive sympatric population structure. Herring spawn mostly in fall or spring, and genomic differentiation was recently detected between these groups. Here we used a subset of this differentiation, 66 single nucleotide polymorphisms (SNPs) to analyze the temporal dynamics of this local adaptation and the applicability of SNP subsets in stock assessment. We showed remarkable temporal stability of genomic differentiation corresponding to spawning season, between samples taken a decade apart (2005 N = 90 vs. 2014 N = 71) in the Gulf of St. Lawrence, and new evidence of limited interbreeding between spawning components. We also examined an understudied and overexploited herring population in Bras d'Or lake (N = 97); using highly reduced SNP panels (N_SNPs > 6), we verified little‐known sympatric spawning populations within this unique inland sea. These results describe consistent local adaptation, arising from asynchronous reproduction in a migratory and dynamic marine species. Our research demonstrates the efficiency and precision of SNP‐based assessments of sympatric subpopulations; and indeed, this temporally stable local adaptation underlines the importance of such fine‐scale management practices.     

Development and preliminary evaluation of a genomewide single nucleotide polymorphisms resource generated by RAD‐seq for the small yellow croaker (Larimichthys polyactis)

Recent advances in high‐throughput sequencing technologies have offered the possibility to generate genomewide sequence data to delineate previously unidentified genetic structure, obtain more accurate estimates of demographic parameters and to evaluate potential adaptive divergence. Here, we identified 27 556 single nucleotide polymorphisms for the small yellow croaker (Larimichthys polyactis) using restriction‐site‐associated DNA (RAD) sequencing of 24 individuals from two populations. Significant sources of genetic variation were identified, with an average nucleotide diversity (π) of 0.00105 ± 0.000425 across individuals, and long‐term effective population size was thus estimated to range between 26 172 and 261 716. According to the results, no differentiation between the two populations was detected based on the SNP data set of top quality score per contig or neutral loci. However, the two analysed populations were highly differentiated based on SNP data set of both top F_ST value per contig and the outlier SNPs. Moreover, local adaptation was highlighted by an F_ST‐based outlier tests implemented in LOSITAN and a total of 538 potentially locally selected SNPs were identified. blast2go annotation of contigs containing the outlier SNPs yielded hits for 37 (66%) of 56 significant blastx matches. Candidate genes for local adaptation constituted a wide array of biological functions, including cellular response to oxidative stress, actin filament binding, ion transmembrane transport and synapse assembly. The generated SNP resources in this study provided a valuable tool for future population genetics and genomics studies of L. polyactis.