Polyceraty (multi‐horns) in Damara sheep maps to ovine chromosome 2

Polyceraty (presence of multiple horns) is rare in modern day ungulates. Although not found in wild sheep, polyceraty does occur in a small number of domestic sheep breeds covering a wide geographical region. Damara are fat‐tailed hair sheep, from the south‐western region of Africa, which display polyceraty, with horn number ranging from zero to four. We conducted a genome‐wide association study for horn number with 43 Damara genotyped with 606 006 SNP markers. The analysis revealed a region with multiple significant SNPs on ovine chromosome 2, in a location different from the mutation for polled in sheep on chromosome 10. The causal mutation for polyceraty was not identified; however, the region associated with polyceraty spans nine HOXD genes, which are critical in embryonic development of appendages. Mutations in HOXD genes are implicated in polydactly phenotypes in mice and humans. There was no evidence for epistatic interactions contributing to polyceraty. This is the first report on the genetic mechanisms underlying polyceraty in the under‐studied Damara.     

Ovine pedomics: the first study of the ovine foot 16S rRNA-based microbiome

We report the first study of the bacterial microbiome of ovine interdigital skin based on 16S rRNA by pyrosequencing and conventional cloning with Sanger-sequencing. Three flocks were selected, one a flock with no signs of footrot or interdigital dermatitis, a second flock with interdigital dermatitis alone and a third flock with both interdigital dermatitis and footrot. The sheep were classified as having either healthy interdigital skin (H) and interdigital dermatitis (ID) or virulent footrot (VFR). The ovine interdigital skin bacterial community varied significantly by flock and clinical condition. The diversity and richness of operational taxonomic units was greater in tissue from sheep with ID than H or VFR-affected sheep. Actinobacteria, Bacteriodetes, Firmicutes and Proteobacteria were the most abundant phyla comprising 25 genera. Peptostreptococcus, Corynebacterium and Staphylococcus were associated with H, ID and VFR, respectively. Sequences of Dichelobacter nodosus, the causal agent of ovine footrot, were not amplified because of mismatches in the 16S rRNA universal forward primer (27F). A specific real-time PCR assay was used to demonstrate the presence of D. nodosus, which was detected in all samples including the flock with no signs of ID or VFR. Sheep with ID had significantly higher numbers of D. nodosus (10⁴–10⁹ cells per g tissue) than those with H or VFR feet.     

Genome-wide association study of footrot in Texel sheep

Background

This is the first study based on a genome-wide association approach that investigates the links between ovine footrot scores and molecular polymorphisms in Texel sheep using the ovine 50 K SNP array (42 883 SNPs (single nucleotide polymorphisms) after quality control). Our aim was to identify molecular predictors of footrot resistance.

Methods

This study used data from animals selected from a footrot-phenotyped Texel sheep population of 2229 sheep with an average of 1.60 scoring records per animal. From these, a subset of 336 animals with extreme trait values for footrot was selected for genotyping based on their phenotypic records. De-regressed estimated breeding values (EBV) for footrot were used as pseudo-phenotypes in the genome-wide association analysis.

Results

Seven SNPs were significant on a chromosome-wise level but the association analysis did not reveal any genome-wise significant SNPs associated with footrot. Based on the current state of knowledge of the ovine genome, it is difficult to clearly link the function of the genes that contain these significant SNPs with a potential role in resistance/susceptibility to footrot. Linkage disequilibrium (LD) was analysed as one of the factors that influence the power of detecting QTL (quantitative trait loci). A mean LD of 0.20 (r² at a distance of 50 kb between two SNPs) in the population analysed was estimated. LD declined from 0.15 to 0.07 and to 0.04 at distances between two SNPs of 100, 1000 and 2000 kb, respectively.

Conclusions

Based on a relatively small number of genotyped animals, this study is a first step to search for genomic regions that are involved in resistance to footrot using the ovine 50 K SNP array. Seven SNPs were found to be significant on a chromosome-wise level. No major genome-wise significant QTL were identified.     

Genome‐wide association reveals the locus responsible for four‐horned ruminant

Phenotypic variability in horn characteristics, such as their size, number and shape, offers the opportunity to elucidate the molecular basis of horn development. The objective of this study was to map the genetic determinant controlling the production of four horns in two breeds, Jacob sheep and Navajo‐Churro, and examine whether an eyelid abnormality occurring in the same populations is related. Genome‐wide association mapping was performed using 125 animals from the two breeds that contain two‐ and four‐horned individuals. A case–control design analysis of 570 712 SNPs genotyped with the ovine HD SNP Beadchip revealed a strong association signal on sheep chromosome 2. The 10 most strongly associated SNPs were all located in a region spanning Mb positions 131.9–132.6, indicating the genetic architecture underpinning the production of four horns is likely to involve a single gene. The closest genes to the most strongly associated marker (OAR2_132568092) were MTX2 and the HOXD cluster, located approximately 93 Kb and 251 Kb upstream respectively. The occurrence of an eyelid malformation across both breeds was restricted to polled animals and those carrying more than two horns. This suggests the eyelid abnormality may be associated with departures from the normal developmental production of two‐horned animals and that the two conditions are developmentally linked. This study demonstrated the presence of separate loci responsible for the polled and four‐horned phenotypes in sheep and advanced our understanding of the complexity that underpins horn morphology in ruminants.     

Genome‐wide screening for DNA variants associated with reading and language traits

Reading and language abilities are heritable traits that are likely to share some genetic influences with each other. To identify pleiotropic genetic variants affecting these traits, we first performed a genome‐wide association scan (GWAS) meta‐analysis using three richly characterized datasets comprising individuals with histories of reading or language problems, and their siblings. GWAS was performed in a total of 1862 participants using the first principal component computed from several quantitative measures of reading‐ and language‐related abilities, both before and after adjustment for performance IQ. We identified novel suggestive associations at the SNPs rs59197085 and rs5995177 (uncorrected P ≈ 10^–7 for each SNP), located respectively at the CCDC136/FLNC and RBFOX2 genes. Each of these SNPs then showed evidence for effects across multiple reading and language traits in univariate association testing against the individual traits. FLNC encodes a structural protein involved in cytoskeleton remodelling, while RBFOX2 is an important regulator of alternative splicing in neurons. The CCDC136/FLNC locus showed association with a comparable reading/language measure in an independent sample of 6434 participants from the general population, although involving distinct alleles of the associated SNP. Our datasets will form an important part of on‐going international efforts to identify genes contributing to reading and language skills.     

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.     

Validation of a cheap and simple nondestructive method for obtaining AFLPs and DNA sequences (mitochondrial and nuclear) in amphibians

The use of nondestructive methods for obtaining DNA from amphibians (e.g. buccal swabs) allows genetic studies to be performed without affecting the survival of the studied individuals. In this study, we compared two methods of nondestructive DNA sampling, buccal swabs and interdigital membrane or toe‐clipping, in several amphibian species of different size: Rhinella spinulosa, R. atacamensis, six species of the genus Telmatobius and Pleurodema thaul. We evaluated the integrity of the DNA extracted by sequencing fragments of mitochondrial and nuclear genes and by generating amplified fragment length polymorphisms markers (AFLPs). In all cases, we obtained an adequate amount of DNA (mean range 55–298 ng/μL). We obtained identical DNA sequences from buccal swab and interdigital membrane/toe‐clip for all individuals. The differences in the coding of AFLP markers between the tissues were similar to those reported for replicas of the same type of sample in similar analyses in other species of amphibians. In conclusion, the use of buccal swabs is a trustworthy and inexpensive method to obtain DNA for mitochondrial and nuclear sequencing and AFLP analyses. Given the types of markers evaluated, buccal swabs may be used for phylogenetic, phylogeographic and population genetic studies, even in small amphibians (<33 mm).     

Haplotype‐based genotyping‐by‐sequencing in oat genome research

In a de novo genotyping‐by‐sequencing (GBS) analysis of short, 64‐base tag‐level haplotypes in 4657 accessions of cultivated oat, we discovered 164741 tag‐level (TL) genetic variants containing 241224 SNPs. From this, the marker density of an oat consensus map was increased by the addition of more than 70000 loci. The mapped TL genotypes of a 635‐line diversity panel were used to infer chromosome‐level (CL) haplotype maps. These maps revealed differences in the number and size of haplotype blocks, as well as differences in haplotype diversity between chromosomes and subsets of the diversity panel. We then explored potential benefits of SNP vs. TL vs. CL GBS variants for mapping, high‐resolution genome analysis and genomic selection in oats. A combined genome‐wide association study (GWAS) of heading date from multiple locations using both TL haplotypes and individual SNP markers identified 184 significant associations. A comparative GWAS using TL haplotypes, CL haplotype blocks and their combinations demonstrated the superiority of using TL haplotype markers. Using a principal component‐based genome‐wide scan, genomic regions containing signatures of selection were identified. These regions may contain genes that are responsible for the local adaptation of oats to Northern American conditions. Genomic selection for heading date using TL haplotypes or SNP markers gave comparable and promising prediction accuracies of up to r = 0.74. Genomic selection carried out in an independent calibration and test population for heading date gave promising prediction accuracies that ranged between r = 0.42 and 0.67. In conclusion, TL haplotype GBS‐derived markers facilitate genome analysis and genomic selection in oat.     

Association between putative functional variants in the PSMB9 gene and risk of melanoma – re‐analysis of published melanoma genome‐wide association studies

To mine possibly hidden causal single‐nucleotide polymorphisms (SNPs) of melanoma, we investigated the association of SNPs in 76 M/G1 transition genes with melanoma risk using our published genome‐wide association study (GWAS) data set with 1804 melanoma cases and 1026 cancer‐free controls. We found multiple SNPs with P < 0.01 and performed validation studies for 18 putative functional SNPs in PSMB9 in two other GWAS data sets. Two SNPs (rs1351383 and rs2127675) were associated with melanoma risk in the GenoMEL data set (P = 0.013 and 0.004, respectively), but failed in validation using the Australian data set. Genotype–phenotype analysis revealed these two SNPs were significantly correlated with mRNA expression level of PSMB9. Further experiments revealed that SNP rs2071480, which is in high LD with rs1351383 and rs2127675, may have a weak effect on the promoter activity of PSMB9. Taken together, our data suggested that functional variants in PSMB9 may contribute to melanoma susceptibility.     

Genome‐wide DNA markers to support genetic management for domestication and commercial production in a large rodent,the Ghanaian grasscutter (Thryonomys swinderianus)

Domestication and commercial production of the grasscutter, Thryonomys swinderianus, a large rodent, represents an important opportunity to secure sustainable animal protein for local communities in West Africa. To support production, DNA markers are required for population diversity assessment, pedigree analysis and marker‐assisted selection. This study reports the application of double‐digest RAD sequencing to simultaneously discover and genotype SNP markers in 24 wild and recently domesticated grasscutters. An initial panel of 1209 SNP loci was characterised from a total of more than 21 000 candidate loci containing single SNPs. This genome‐wide resource represents the first application of its type to commercial production of a large rodent for food and advances the use of agricultural genomics in Ghana.     

Mutations in Ovis aries TMEM154 are associated with lower small ruminant lentivirus proviral concentration in one sheep flock

Small ruminant lentivirus (SRLV), also called ovine progressive pneumonia virus or maedi‐visna, is present in 24% of US sheep. Like human immunodeficiency virus, SRLV is a macrophage‐tropic lentivirus that causes lifelong infection. The production impacts from SRLV are due to a range of disease symptoms, including pneumonia, arthritis, mastitis, body condition wasting and encephalitis. There is no cure and no effective vaccine for preventing SRLV infection. However, breed differences in prevalence and proviral concentration indicate a genetic basis for susceptibility to SRLV. Animals with high blood proviral concentration show increased tissue lesion severity, so proviral concentration represents a live animal test for control post‐infection in terms of proviral replication and disease severity. Recently, it was found that sheep with two copies of TMEM154 haplotype 1 (encoding lysine at position 35) had lower odds of SRLV infection. In this study, we examined the relationship between SRLV control post‐infection and variants in two genes, TMEM154 and CCR5, in four flocks containing 1403 SRLV‐positive sheep. We found two copies of TMEM154 haplotype 1 were associated with lower SRLV proviral concentration in one flock (P < 0.02). This identified the same favorable diplotype for SRLV control post‐infection as for odds of infection. However, frequencies of haplotypes 2 and 3 were too low in the other three flocks to test. The CCR5 promoter deletion did not have consistent association with SRLV proviral concentration. Future work in flocks with more balanced allele frequencies is needed to confirm or refute TMEM154 association with control of SRLV post‐infection.     

The brown coat colour of Coppernecked goats is associated with a non‐synonymous variant at the TYRP1 locus on chromosome 8

The recent development of a goat SNP genotyping microarray enables genome‐wide association studies in this important livestock species. We investigated the genetic basis of the black and brown coat colour in Valais Blacknecked and Coppernecked goats. A genome‐wide association analysis using goat SNP50 BeadChip genotypes of 22 cases and 23 controls allowed us to map the locus for the brown coat colour to goat chromosome 8. The TYRP1 gene is located within the associated chromosomal region, and TYRP1 variants cause similar coat colour phenotypes in different species. We thus considered TYRP1 as a strong positional and functional candidate. We resequenced the caprine TYRP1 gene by Sanger and Illumina sequencing and identified two non‐synonymous variants, p.Ile478Thr and p.Gly496Asp, that might have a functional impact on the TYRP1 protein. However, based on the obtained pedigree and genotype data, the brown coat colour in these goats is not due to a single recessive loss‐of‐function allele. Surprisingly, the genotype distribution and the pedigree data suggest that the ⁴⁹⁶Asp allele might possibly act in a dominant manner. The ⁴⁹⁶Asp allele was present in 77 of 81 investigated Coppernecked goats and did not occur in black goats. This strongly suggests heterogeneity underlying the brown coat colour in Coppernecked goats. Functional experiments or targeted matings will be required to verify the unexpected preliminary findings.     

Immunochip analysis identifies association of the RAD50/IL13 region with human longevity

Human longevity is characterized by a remarkable lack of confirmed genetic associations. Here, we report on the identification of a novel locus for longevity in the RAD50/IL13 region on chromosome 5q31.1 using a combined European sample of 3208 long‐lived individuals (LLI) and 8919 younger controls. First, we performed a large‐scale association study on 1458 German LLI (mean age 99.0 years) and 6368 controls (mean age 57.2 years) by targeting known immune‐associated loci covered by the Immunochip. The analysis of 142 136 autosomal single nucleotide polymorphisms (SNPs) revealed an Immunochip‐wide significant signal (P_Immunochip = 7.01 × 10^–9) for the SNP rs2075650 in the TOMM40/APOE region, which has been previously described in the context of human longevity. To identify novel susceptibility loci, we selected 15 markers with P_Immunochip < 5 × 10^–4 for replication in two samples from France (1257 LLI, mean age 102.4 years; 1811 controls, mean age 49.1 years) and Denmark (493 LLI, mean age 96.2 years; 740 controls, mean age 63.1 years). The association at SNP rs2706372 replicated in the French study collection and showed a similar trend in the Danish participants and was also significant in a meta‐analysis of the combined French and Danish data after adjusting for multiple testing. In a meta‐analysis of all three samples, rs2706372 reached a P‐value of P_{Immunochip+Repl} = 5.42 × 10^?7 (OR = 1.20; 95% CI = 1.12–1.28). SNP rs2706372 is located in the extended RAD50/IL13 region. RAD50 seems a plausible longevity candidate due to its involvement in DNA repair and inflammation. Further studies are needed to identify the functional variant(s) that predispose(s) to a long and healthy life.     

New resources for genetic studies in Populus nigra: genome‐wide SNP discovery and development of a 12k Infinium array

Whole genome resequencing of 51 Populus nigra (L.) individuals from across Western Europe was performed using Illumina platforms. A total number of 1 878 727 SNPs distributed along the P. nigra reference sequence were identified. The SNP calling accuracy was validated with Sanger sequencing. SNPs were selected within 14 previously identified QTL regions, 2916 expressional candidate genes related to rust resistance, wood properties, water‐use efficiency and bud phenology and 1732 genes randomly spread across the genome. Over 10 000 SNPs were selected for the construction of a 12k Infinium Bead‐Chip array dedicated to association mapping. The SNP genotyping assay was performed with 888 P. nigra individuals. The genotyping success rate was 91%. Our high success rate was due to the discovery panel design and the stringent parameters applied for SNP calling and selection. In the same set of P. nigra genotypes, linkage disequilibrium throughout the genome decayed on average within 5–7 kb to half of its maximum value. As an application test, ADMIXTURE analysis was performed with a selection of 600 SNPs spread throughout the genome and 706 individuals collected along 12 river basins. The admixture pattern was consistent with genetic diversity revealed by neutral markers and the geographical distribution of the populations. These newly developed SNP resources and genotyping array provide a valuable tool for population genetic studies and identification of QTLs through natural‐population based genetic association studies in P. nigra.     

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.     

Copy number variations of CBF genes at the Fr‐A2 locus are essential components of winter hardiness in wheat

Winter hardiness is important for the adaptation of wheat to the harsh winter conditions in temperate regions and is thus also an important breeding goal. Here, we employed a panel of 407 European winter wheat cultivars to dissect the genetic architecture of winter hardiness. We show that copy number variation (CNV) of CBF (C‐repeat Binding Factor) genes at the Fr‐A2 locus is the essential component for winter survival, with CBF‐A14 CNV being the most likely causal polymorphism, accounting for 24.3% of the genotypic variance. Genome‐wide association mapping identified several markers in the Fr‐A2 chromosomal region, which even after accounting for the effects of CBF‐A14 copy number explained approximately 15% of the genotypic variance. This suggests that additional, as yet undiscovered, polymorphisms are present at the Fr‐A2 locus. Furthermore, CNV of Vrn‐A1 explained an additional 3.0% of the genotypic variance. The allele frequencies of all loci associated with winter hardiness were found to show geographic patterns consistent with their role in adaptation. Collectively, our results from the candidate gene analysis, association mapping and genome‐wide prediction show that winter hardiness in wheat is a quantitative trait, but with a major contribution of the Fr‐A2 locus.     

Genome‐wide detection of genetic loci and candidate genes for teat number and body conformation traits at birth in Chinese Sushan pigs

Body conformation at birth and teat number are economically important traits in the pig industry, as these traits are usually explored to evaluate the growth and reproductive potential of piglets. To detect genetic loci and candidate genes for these traits, we performed a GWAS on 269 pigs from a recently developed Chinese breed (Sushan) using 38  128 informative SNPs on the Affymetrix Porcine SNP 55K Array. In total, we detected one genome‐wide significant (P = 1.31e‐6) SNP for teat number on chromosome X and 15 chromosome‐wide significant SNPs for teat number, body weight, body length, chest circumference and cannon circumference at birth on chromosomes 1, 3, 4, 6, 7, 9, 10, 13, 14, 15, 17 and 18. The most significant SNP had an additive effect of 0.74 × total teat number, explaining 20% of phenotypic variance. Five significant SNPs resided in the previously reported quantitative trait loci for these traits and seven significant SNPs had a pleiotropic effect on multiple traits. Intriguingly, 12 of the genes nearest to the significant SNPs are functionally related to body conformation and teat number traits, including SPRED2, MKX, TMSB4X and ESR1. GO analysis revealed that candidate genes proximal to the significant SNPs were enriched in the G‐protein coupled receptor and steroid hormone‐mediated signaling pathway. Our findings shed light on the genetic basis of the measured traits and provide molecular markers especially for the genetic improvement of teat number in Sushan and related pigs.