Primers for sequence characterization and polymorphism detection in the Atlantic salmon (Salmo salar) growth hormone 1 (GH1) gene

We report the identification of intraspecific sequence variation in the Atlantic salmon (Salmo salar) growth hormone 1 gene. Rapid and inexpensive assays for polymorphism detection were developed for 10 sites. Five of the assays detected single nucleotide polymorphisms (SNPs) using polymerase chain reaction–restriction fragment length polymorphism (PCR–RFLP) analyses, and five were indel polymorphisms, detected using fragment length analyses. The average within population frequency of the most common allele varied from 0.52 to 0.90, and the average within population heterozygosity varied from 0.02 to 0.37 in seven European salmon populations.     

The first report of polymorphisms and genetic characteristics of the prion protein gene (PRNP) in horses

ABSTRACT

Prion diseases have a wide host range, but prion-infected cases have never been reported in horses. Genetic polymorphisms that can directly impact the structural stability of horse prion protein have not been investigated thus far. In addition, we noticed that previous studies focusing on horse-specific amino acids and secondary structure predictions of prion protein were performed for limited parts of the protein. In this study, we found genetic polymorphisms in the horse prion protein gene (PRNP) in 201 Thoroughbred horses. The identified polymorphism was assessed to determine whether this polymorphism impedes stability of protein using PolyPhen-2, PROVEAN and PANTHER. In addition, we evaluated horse-specific amino acids in horse and mouse prion proteins using same methods. We found only one single nucleotide polymorphism (SNP) in the horse prion protein, and three annotation tools predicted that the SNP is benign. In addition, horse-specific amino acids showed different effects on horse and mouse prion proteins, respectively.

Abbreviations: PRNP: prion protein gene; SNP: single nucleotide polymorphism; CJD: Creutzfeldt-Jakob disease; CWD: chronic wasting disease; TME: transmissible mink encephalopathy; FSE: feline spongiform encephalopathy; MD: molecular dynamics; ER: endoplasmic reticulum; GPI: glycosylphosphatidylinositol; NMR: nuclear magnetic resonance; ORF: open reading frame; GWAS: genome-wide association study; NAPA: non-adaptive prion amplification; HMM: hidden Markov model; NCBI: National Center for Biotechnology Information     

Genome‐wide SNP analysis unveils genetic structure and phylogeographic history of snow sheep (Ovis nivicola) populations inhabiting the Verkhoyansk Mountains and Momsky Ridge (northeastern Siberia)

Insights into the genetic characteristics of a species provide important information for wildlife conservation programs. Here, we used the OvineSNP50 BeadChip developed for domestic sheep to examine population structure and evaluate genetic diversity of snow sheep (Ovis nivicola) inhabiting Verkhoyansk Range and Momsky Ridge. A total of 1,121 polymorphic SNPs were used to test 80 specimens representing five populations, including four populations of the Verkhoyansk Mountain chain: Kharaulakh Ridge–Tiksi Bay (TIK, n = 22), Orulgan Ridge (ORU, n = 22), the central part of Verkhoyansk Range (VER, n = 15), Suntar‐Khayata Ridge (SKH, n = 13), and Momsky Ridge (MOM, n = 8). We showed that the studied populations were genetically structured according to a geographic pattern. Pairwise F_ST values ranged from 0.044 to 0.205. Admixture analysis identified K = 2 as the most likely number of ancestral populations. A Neighbor‐Net tree showed that TIK was an isolated group related to the main network through ORU. TreeMix analysis revealed that TIK and MOM originated from two different ancestral populations and detected gene flow from MOM to ORU. This was supported by the f3 statistic, which showed that ORU is an admixed population with TIK and MOM/SKH heritage. Genetic diversity in the studied groups was increasing southward. Minimum values of observed (Ho) and expected (He) heterozygosity and allelic richness (Ar) were observed in the most northern population—TIK, and maximum values were observed in the most southern population—SKH. Thus, our results revealed clear genetic structure in the studied populations of snow sheep and showed that TIK has a different origin from MOM, SKH, and VER even though they are conventionally considered a single subspecies known as Yakut snow sheep (Ovis nivicola lydekkeri). Most likely, TIK was an isolated group during the Late Pleistocene glaciations of Verkhoyansk Range.     

A single‐nucleotide polymorphism‐based approach for rapid and cost‐effective genetic wolf monitoring in Europe based on noninvasively collected samples

Noninvasive genetics based on microsatellite markers has become an indispensable tool for wildlife monitoring and conservation research over the past decades. However, microsatellites have several drawbacks, such as the lack of standardisation between laboratories and high error rates. Here, we propose an alternative single‐nucleotide polymorphism (SNP)‐based marker system for noninvasively collected samples, which promises to solve these problems. Using nanofluidic SNP genotyping technology (Fluidigm), we genotyped 158 wolf samples (tissue, scats, hairs, urine) for 192 SNP loci selected from the Affymetrix v2 Canine SNP Array. We carefully selected an optimised final set of 96 SNPs (and discarded the worse half), based on assay performance and reliability. We found rates of missing data in this SNP set of <10% and genotyping error of ~1%, which improves genotyping accuracy by nearly an order of magnitude when compared to published data for other marker types. Our approach provides a tool for rapid and cost‐effective genotyping of noninvasively collected wildlife samples. The ability to standardise genotype scoring combined with low error rates promises to constitute a major technological advancement and could establish SNPs as a standard marker for future wildlife monitoring.     

The easy road to genome‐wide medium density SNP screening in a non‐model species: development and application of a 10 K SNP‐chip for the house sparrow (Passer domesticus)

With the advent of next generation sequencing, new avenues have opened to study genomics in wild populations of non‐model species. Here, we describe a successful approach to a genome‐wide medium density Single Nucleotide Polymorphism (SNP) panel in a non‐model species, the house sparrow (Passer domesticus), through the development of a 10 K Illumina iSelect HD BeadChip. Genomic DNA and cDNA derived from six individuals were sequenced on a 454 GS FLX system and generated a total of 1.2 million sequences, in which SNPs were detected. As no reference genome exists for the house sparrow, we used the zebra finch (Taeniopygia guttata) reference genome to determine the most likely position of each SNP. The 10 000 SNPs on the SNP‐chip were selected to be distributed evenly across 31 chromosomes, giving on average one SNP per 100 000 bp. The SNP‐chip was screened across 1968 individual house sparrows from four island populations. Of the original 10 000 SNPs, 7413 were found to be variable, and 99% of these SNPs were successfully called in at least 93% of all individuals. We used the SNP‐chip to demonstrate the ability of such genome‐wide marker data to detect population sub‐division, and compared these results to similar analyses using microsatellites. The SNP‐chip will be used to map Quantitative Trait Loci (QTL) for fitness‐related phenotypic traits in natural populations.     

Genome‐wide association study identifies variants in the CAPN9 gene associated with umbilical hernia in pigs

Pig umbilical hernia (UH) affects pig welfare and brings considerable economic loss to the pig industry. To date, the molecular mechanisms underlying pig UH are still poorly understood. To identify potential loci for susceptibility to this disease, we performed a genome‐wide association study in an Erhualian × Shaziling F₂ intercross population. A total of 45 animals were genotyped using Illumina Porcine SNP60 BeadChips. We observed a SNP (rs80993347) located in the calpain‐9 (CAPN9) gene on Sus scrofa chromosome 14 that was significantly associated with UH (P = 1.97 × 10^?10). Then, we identified a synonymous mutation rs321865883 (g.20164T>C) in exon 10 of the CAPN9 gene that distinguished two affected individuals (CC) from their normal full‐sibs (TC). Finally, quantitative polymerase chain reaction was explored to investigate the mRNA expression profile of the CAPN9 gene in 12 tissues in Yorkshire pigs at different developmental stages (3, 90 and 180 days). CAPN9 showed high expression levels in the gastrointestinal tract at these three growth stages. The results of this study indicate that the CAPN9 gene might be implicated in UH. Further studies are required to establish a role of CAPN9 in pig UH.     

Estimation of linkage disequilibrium and interspecific gene flow in Ficedula flycatchers by a newly developed 50k single‐nucleotide polymorphism array

With the access to draft genome sequence assemblies and whole‐genome resequencing data from population samples, molecular ecology studies will be able to take truly genome‐wide approaches. This now applies to an avian model system in ecological and evolutionary research: Old World flycatchers of the genus Ficedula, for which we recently obtained a 1.1 Gb collared flycatcher genome assembly and identified 13 million single‐nucleotide polymorphism (SNP)s in population resequencing of this species and its sister species, pied flycatcher. Here, we developed a custom 50K Illumina iSelect flycatcher SNP array with markers covering 30 autosomes and the Z chromosome. Using a number of selection criteria for inclusion in the array, both genotyping success rate and polymorphism information content (mean marker heterozygosity = 0.41) were high. We used the array to assess linkage disequilibrium (LD) and hybridization in flycatchers. Linkage disequilibrium declined quickly to the background level at an average distance of 17 kb, but the extent of LD varied markedly within the genome and was more than 10‐fold higher in ‘genomic islands’ of differentiation than in the rest of the genome. Genetic ancestry analysis identified 33 F₁ hybrids but no later‐generation hybrids from sympatric populations of collared flycatchers and pied flycatchers, contradicting earlier reports of backcrosses identified from much fewer number of markers. With an estimated divergence time as recently as <1 Ma, this suggests strong selection against F₁ hybrids and unusually rapid evolution of reproductive incompatibility in an avian system.     

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.     

A targeted genotyping approach enhances identification of variants in taste receptor and appetite/reward genes of potential functional importance for obesity‐related porcine traits

Taste receptors (TASRs) and appetite and reward (AR) mechanisms influence eating behaviour, which in turn affects food intake and risk of obesity. In a previous study, we used next generation sequencing to identify potentially functional mutations in TASR and AR genes and found indications for genetic associations between identified variants and growth and fat deposition in a subgroup of animals (n = 38) from the UNIK resource pig population. This population was created for studying obesity and obesity‐related diseases. In the present study we validated results from our previous study by investigating genetic associations between 24 selected single nucleotide variants in TASR and AR gene variants and 35 phenotypes describing obesity and metabolism in the entire UNIK population (n = 564). Fifteen variants showed significant association with specific obesity‐related phenotypes after Bonferroni correction. Six of the 15 genes, namely SIM1, FOS, TAS2R4, TAS2R9, MCHR2 and LEPR, showed good correlation between known biological function and associated phenotype. We verified a genetic association between potentially functional variants in TASR/AR genes and growth/obesity and conclude that the combination of identification of potentially functional variants by next generation sequencing followed by targeted genotyping and association studies is a powerful and cost‐effective approach for increasing the power of genetic association studies.     

Genome‐wide analyses suggest parallel selection for universal traits may eclipse local environmental selection in a highly mobile carnivore

Ecological and environmental heterogeneity can produce genetic differentiation in highly mobile species. Accordingly, local adaptation may be expected across comparatively short distances in the presence of marked environmental gradients. Within the European continent, wolves (Canis lupus) exhibit distinct north–south population differentiation. We investigated more than 67‐K single nucleotide polymorphism (SNP) loci for signatures of local adaptation in 59 unrelated wolves from four previously identified population clusters (northcentral Europe n = 32, Carpathian Mountains n = 7, Dinaric‐Balkan n = 9, Ukrainian Steppe n = 11). Our analyses combined identification of outlier loci with findings from genome‐wide association study of individual genomic profiles and 12 environmental variables. We identified 353 candidate SNP loci. We examined the SNP position and neighboring megabase (1 Mb, one million bases) regions in the dog (C. lupus familiaris) genome for genes potentially under selection, including homologue genes in other vertebrates. These regions included functional genes for, for example, temperature regulation that may indicate local adaptation and genes controlling for functions universally important for wolves, including olfaction, hearing, vision, and cognitive functions. We also observed strong outliers not associated with any of the investigated variables, which could suggest selective pressures associated with other unmeasured environmental variables and/or demographic factors. These patterns are further supported by the examination of spatial distributions of the SNPs associated with universally important traits, which typically show marked differences in allele frequencies among population clusters. Accordingly, parallel selection for features important to all wolves may eclipse local environmental selection and implies long‐term separation among population clusters.     

Effects of hydroxy‐delta‐5‐steroid dehydrogenase, 3 beta‐ and steroid delta‐isomerase 1 polymorphisms on fat androstenone level and gene expression in Duroc pigs

A high level of androstenone in porcine adipose tissue is a major factor contributing to boar taint. Porcine hydroxy‐delta‐5‐steroid dehydrogenase, 3 beta‐ and steroid delta‐isomerase 1 (3β‐HSD, also known as HSD3B1) plays a key role in the hepatic metabolism that catalyzes androstenone to β‐androstenol. Therefore, 3β‐HSD is a candidate gene for boar taint. This study aimed to investigate functional 3β‐HSD polymorphisms in Duroc pigs. We found eight single nucleotide polymorphisms (SNPs) in the full‐length porcine 3β‐HSD. Four of the SNPs had restriction enzyme sites, and we genotyped them in 147 uncastrated male Duroc pigs using a polymerase chain reaction–restriction fragment length polymorphism method. Pigs with the GG genotype at the g.165262G>A locus (SNP5) had significantly lower androstenone levels than did those with other genotypes (P = 0.030). SNP5 also was associated with differences in 3β‐HSD mRNA levels: pigs with the GG genotype had higher levels than those with other genotypes (P = 0.019). The SNP5 polymorphism could affect the hepatic catabolism of androstenone and consequently impact androstenone accumulation in the adipose tissue. Therefore, SNP5 in the 3β‐HSD of Duroc pigs could be a useful selective marker for decreasing boar taint.     

Replication factor C1 (RFC1) is required for double‐strand break repair during meiotic homologous recombination in Arabidopsis

Replication factor C1 (RFC1), which is conserved in eukaryotes, is involved in DNA replication and checkpoint control. However, a RFC1 product participating in DNA repair at meiosis has not been reported in Arabidopsis. Here, we report functional characterization of AtRFC1 through analysis of the rfc1–2 mutant. The rfc1–2 mutant displayed normal vegetative growth but showed silique sterility because the male gametophyte was arrested at the uninucleus microspore stage and the female at the functional megaspore stage. Expression of AtRFC1 was concentrated in the reproductive organ primordia, meiocytes and developing gametes. Chromosome spreads showed that pairing and synapsis were normal, and the chromosomes were broken when desynapsis began at late prophase I, and chromosome fragments remained in the subsequent stages. For this reason, homologous chromosomes and sister chromatids segregated unequally, leading to pollen sterility. Immunolocalization revealed that the AtRFC1 protein localized to the chromosomes during zygotene and pachytene in wild‐type but were absent in the spo11–1 mutant. The chromosome fragmentation of rfc1–2 was suppressed by spo11–1, indicating that AtRFC1 acted downstream of AtSPO11‐1. The similar chromosome behavior of rad51 rfc1–2 and rad51 suggests that AtRFC1 may act with AtRAD51 in the same pathway. In summary, AtRFC1 is required for DNA double‐strand break repair during meiotic homologous recombination of Arabidopsis.