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.     

A genome‐wide significant association on chromosome 2 for footrot resistance/susceptibility in Swiss White Alpine sheep

Footrot is one of the most important causes of lameness in global sheep populations and is characterized by a bacterial infection of the interdigital skin. As a multifactorial disease, its clinical representation depends not only on pathogen factors and environmental components but also on the individual resistance/susceptibility of the host. A genetic component has been shown in previous studies; however, so far no causative genetic variant influencing the risk of developing footrot has been identified. In this study, we genotyped 373 Swiss White Alpine sheep, using the ovine high‐density 600k SNP chip, in order to run a DNA‐based comparison of individuals with known clinical footrot status. We performed a case–control genome‐wide association study, which revealed a genome‐wide significant association for SNP rs418747104 on ovine chromosome 2 at 81.2 Mb. The three best associated SNP markers were located at the MPDZ gene, which codes for the multiple PDZ domain crumbs cell polarity complex component protein, also known as multi‐PDZ domain protein 1 (MUPP1). This protein is possibly involved in maintaining the barrier function and integrity of tight junctions. Therefore, we speculate that individuals carrying MPDZ variants may differ in their footrot resistance/susceptibility due to modified horn and interdigital skin integrity. In conclusion, our study reveals that MPDZ might represent a functional candidate gene, and further research is needed to explore its role in footrot affected sheep.     

Validation of genetic polymorphisms on BTA14 associated with carcass trait in a commercial Hanwoo population

The objective of this study was to validate the association of significant SNPs identified from a previous genome‐wide association study with carcass weight (CWT) in a commercial Hanwoo population. We genotyped 13 SNPs located on BTA14 in 867 steers from Korea Hanwoo feedlot bulls. Of these 13 SNPs, five SNPs, namely rs29021868, rs110061498, rs109546980, rs42404006 and rs42303720, were found to be significantly associated (P < 0.001) with CWT. These five significant markers spanned the 24.3 to 29.4 Mb region of BTA14. The most significant marker (rs29021868) for CWT in this study had a 13.07 kg allele substitution effect and accounted for 2.4% of the additive genetic variance in the commercial Hanwoo population. The SNP marker rs109546980 was found to be significantly associated with both CWT (P < 0.001) and eye muscle area (P < 0.001) and could potentially be exploited for marker‐assisted selection in Hanwoo cattle. We also genotyped the ss319607402 variation, which maps to intron2 of PLAG1 gene and which is already reported to be associated with height, to identify any significant association with carcass weight; however, no such association was observed in this Hanwoo commercial population.     

Effect of swine leukocyte antigen‐DQA gene variation on diarrhea in Large White,Landrace, and Duroc piglets

Piglet diarrhea is one of the most common factors that affects the benefits of the swine industry. Although recent studies have shown that exon 2 of SLA‐DQA is associated with piglet resistance to diarrhea, contributions of genetic variation in the additional exon coding regions of this gene remain unclear. Here, we investigated variation in exons 1, 3 and 4 of the SLA‐DQA gene and evaluated their effects on diarrheal infection in 425 suckling piglets. No variation was identified in exon 1. In exon 3, there were eight alleles detected, generated by 14 single nucleotide polymorphisms (SNPs) and three nucleotide deletions, eight SNPs being newly identified. Four allele sequences and three SNPs were identified in exon 4, only one SNP being newly identified. Statistical analysis showed that the genotypes of exon 3 are significantly associated with piglet diarrhea; indeed, genotypes DQA*wb01/wb02 and wb04/wb05 are clearly associated with resistance to piglet diarrhea, as they have the lowest probabilities of infection (P < 0.05). However, no significant association was found between the genotypes of exon 4 and diarrhea (P > 0.05). These results provide important new information concerning the level of genetic diversity at the SLA‐DQA locus and suggest that further genetic association studies of piglet diarrhea resistance should include analyses of both exons 2 and 3 of this locus.     

Genome‐wide association analysis and gene set enrichment analysis with SNP data identify genes associated with 305‐day milk yield in Holstein dairy cows

Milk production traits, such as 305‐day milk yield (305MY), have been under direct selection to improve production in dairy cows. Over the past 50 years, the average milk yield has nearly doubled, and over 56% of the increase is attributable to genetic improvement. As such, additional improvements in milk yield are still possible as new loci are identified. The objectives of this study were to detect SNPs and gene sets associated with 305MY in order to identify new candidate genes contributing to variation in milk production. A population of 781 primiparous Holstein cows from six central Washington dairies with records of 305MY and energy corrected milk were used to perform a genome‐wide association analysis (GWAA) using the Illumina BovineHD BeadChip (777 962 SNPs) to identify QTL associated with 305MY (P < 1.0 × 10^?5). A gene set enrichment analysis with SNP data (GSEA‐SNP) was performed to identify gene sets (normalized enrichment score > 3.0) and leading edge genes (LEGs) influencing 305MY. The GWAA identified three QTL comprising 34 SNPs and 30 positional candidate genes. In the GSEA‐SNP, five gene sets with 58 unique and 24 shared LEGs contributed to 305MY. Identification of QTL and LEGs associated with 305MY can provide additional targets for genomic selection to continue to improve 305MY in dairy cattle.     

Polymorphisms in calpastatin and mu‐calpain genes are associated with beef iron content

The objective of this study was to assess the association of markers in the calpastatin and mu‐calpain loci with iron in beef cattle muscle. The population consisted of 259 cross‐bred steers from Beefmaster, Brangus, Bonsmara, Romosinuano, Hereford and Angus sires. Total iron and heme iron concentrations were measured. Markers in the calpastatin (referred to as CAST) and mu‐calpain (referred to as CAPN4751) genes were used to assess their association with iron levels. The mean and standard error for iron and heme iron content in the population was 35.6 ± 1.3 μg and 27.1 ± 1.4 μg respectively. Significant associations (P < 0.01) of markers were observed for both iron and heme iron content. For CAST, animals with the CC genotype had higher levels of iron and heme iron in longissimus dorsi muscle. For CAPN4751, individuals with the TT genotype had higher concentrations of iron and heme iron than did animals with the CC and CT genotypes. Genotypes known to be associated with tougher meat were associated with higher levels of iron concentration.     

Genetic dissection of pre‐anthesis sub‐phase durations during the reproductive spike development of wheat

Flowering time is an important factor affecting grain yield in wheat. In this study, we divided reproductive spike development into eight sub‐phases. These sub‐phases have the potential to be delicately manipulated to increase grain yield. We measured 36 traits with regard to sub‐phase durations, determined three grain yield‐related traits in eight field environments and mapped 15 696 single nucleotide polymorphism (SNP, based on 90k Infinium chip and 35k Affymetrix chip) markers in 210 wheat genotypes. Phenotypic and genetic associations between grain yield traits and sub‐phase durations showed significant consistency (Mantel test; r = 0.5377, P < 0.001). The shared quantitative trait loci (QTLs) revealed by the genome‐wide association study suggested a close association between grain yield and sub‐phase duration, which may be attributed to effects on spikelet initiation/spikelet number (double ridge to terminal spikelet stage, DR‐TS) and assimilate accumulation (green anther to anthesis stage, GA‐AN). Moreover, we observed that the photoperiod‐sensitivity allele at the Ppd‐D1 locus on chromosome 2D markedly extended all sub‐phase durations, which may contribute to its positive effects on grain yield traits. The dwarfing allele at the Rht‐D1 (chromosome 4D) locus altered the sub‐phase duration and displayed positive effects on grain yield traits. Data for 30 selected genotypes (from among the original 210 genotypes) in the field displayed a close association with that from the greenhouse. Most importantly, this study demonstrated specific connections to grain yield in narrower time windows (i.e. the eight sub‐phases), rather than the entire stem elongation phase as a whole.     

Limited dispersal in an ectoparasitic mite,Laelaps giganteus,contributes to significant phylogeographic congruence with the rodent host,Rhabdomys

To explore how biogeography, parasite life history and host vagility influences evolutionary codivergences, we followed a comparative phylogeography approach using a host‐specific nonpermanent mite, Laelaps giganteus, that occurs on four rodent species within the genus Rhabdomys. A mtDNA COI haplotype network derived for 278 parasite specimens showed marked phylogeographic congruence with host distributions. Analysis of the less variable nuclear intron Tropomyosin was in part consistent with these results. Although distance‐based cophylogenetic analyses in axparafit failed to support significant mtDNA codivergences (P ≥ 0.02), event‐based analyses revealed significant cophylogeny between sampling localities of Rhabdomys and Laelaps using core‐pa (P = 0.046) and jane (P = 0.026; P = 0.00). These findings, in conjunction with the weak congruence previously reported among the permanent ectoparasitic lice Polyplax and Rhabdomys, suggest that host–parasite intimacy is not the most important driver of significant codivergence in our study system. Instead, the more restricted dispersal ability of L. giganteus, when compared to Polyplax, resulted in stronger spatial structuring and this could have resulted in significant codivergence. Host switching occurred predominantly on the edges of host distributions and was probably facilitated by climate‐induced range shifts. When host ranges shift, the phylogeographic structure of L. giganteus is not reflecting the host movements as most of the nest bound parasites do not disperse with the host (they miss the boat) and the genetic contribution of the few dispersing mite individuals is often overwhelmed by the large number of individuals already present in nests within the new environment (causing them to drown on arrival).     

Directed dispersal by rotational shepherding supports landscape genetic connectivity in a calcareous grassland plant

Directed dispersal by animal vectors has been found to have large effects on the structure and dynamics of plant populations adapted to frugivory. Yet, empirical data are lacking on the potential of directed dispersal by rotational grazing of domestic animals to mediate gene flow across the landscape. Here, we investigated the potential effect of large‐flock shepherding on landscape‐scale genetic structure in the calcareous grassland plant Dianthus carthusianorum, whose seeds lack morphological adaptations to dispersal to animals or wind. We found a significant pattern of genetic structure differentiating population within grazed patches of three nonoverlapping shepherding systems and populations of ungrazed patches. Among ungrazed patches, we found a strong and significant effect of isolation by distance (r = 0.56). In contrast, genetic distance between grazed patches within the same herding system was unrelated to geographical distance but significantly related to distance along shepherding routes (r = 0.44). This latter effect of connectivity along shepherding routes suggests that gene flow is spatially restricted occurring mostly between adjacent populations. While this study used nuclear markers that integrate gene flow by pollen and seed, the significant difference in the genetic structure between ungrazed patches and patches connected by large‐flock shepherding indicates the potential of directed seed dispersal by sheep across the landscape.     

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.     

Identifying genetic loci controlling neonatal passive transfer of immunity using a hybrid genotyping strategy

Colostrum intake is critical to a piglet's survival and can be measured by precipitating out the γ‐immunoglobulins from serum with ammonium sulfate (immunocrit). Genetic analysis of immunocrits on 5312 piglets indicated that the heritabilities (se) for direct and maternal effects were 0.13 (0.06) and 0.53 (0.08) respectively. To identify QTL for direct genetic effects, piglets with the highest and lowest immunocrits from 470 litters were selected. Six sets of DNA pools were created based on sire of the litter. These 12 DNA pools were applied to Illumina Porcine SNP60 BeadChips. Normalized X and Y values were analyzed. Three different SNP selection methods were used: deviation of the mean from high vs. low pools, the deviation adjusted for variance based on binomial theory and ANOVA. The 25 highest ranking SNPs were selected from each evaluation for further study along with 12 regions selected based on a five‐SNP window approach. Selected SNPs were individually genotyped in the 988 piglets included in pools as well as in 524 piglets that had intermediate immunocrits. Association analyses were conducted fitting an animal model using the estimated genetic parameters. Nineteen SNPs were nominally associated (P < 0.01) with immunocrit values, of which nine remained significant (P < 0.05) after Bonferroni correction, located in 16 genomic regions on 13 chromosomes. In conclusion, the pooling strategy reduced the cost to scan the genome by more than 80% and identified genomic regions associated with a piglet's ability to acquire γ‐immunoglobulin from colostrum. Each method to rank SNPs from the pooled analyses contributed unique validated markers, suggesting that multiple analyses will reveal more QTL than a single analysis.     

SNPs in SNCA,MCCC1, DLG2, GBF1 and MBNL2 are associated with Parkinson's disease in southern Chinese population

Numerous single nucleotide polymorphisms (SNPs), which have been identified as susceptibility factors for Parkinson's disease (PD) as per genome‐wide association studies, have not been fully characterized for PD patients in China. This study aimed to replicate the relationship between 12 novel SNPs of 12 genes and PD risk in southern Chinese population. Twelve SNPs of 12 genes were detected in 231 PD patients and 249 controls, using the SNaPshot technique. Meta‐analysis was used to assess heterogeneity of effect sizes between this study and published data. The impact of SNPs on gene expression was investigated by analysing the SNP‐gene association in the expression quantitative trait loci (eQTL) data sets. rs8180209 of SNCA (allele model: P = .047, OR = 0.77; additive model: P = .047, OR = 0.77), rs2270968 of MCCC1 (dominant model: P = .024, OR = 1.52), rs7479949 of DLG2 (recessive model; P = .019, OR = 1.52), rs10748818 of GBF1 (additive model: P < .001, OR = 0.37), and rs4771268 of MBNL2 (recessive model: P = .003, OR = 0.48) were replicated to be significantly associated with the increased risk of PD. Noteworthy, a meta‐analysis of previous studies suggested rs8180209, rs2270968, rs7479949 and rs4771268 were in line with those of our cohort. Our study replicated five novel functional SNPs in SNCA, MCCC1, DLG2, GBF1 and MBNL2 could be associated with increased risk of PD in southern Chinese population.     

Genetic markers that influence feed efficiency phenotypes also affect cattle temperament as measured by flight speed

Flight speed is a predictive indicator of cattle temperament and is associated with feed efficiency phenotypes. Genetic markers associated with both traits may assist with selection of calmer animals with improved economic value. A preliminary genome‐wide association study determined chromosomal regions on BTA9, and 17 were associated with flight speed. The genes quaking (QKI), glutamate receptor, ionotropic, AMPA 2 (GRIA2) and glycine receptor β (GLRB) were identified in these regions as potential functional candidates. Beef steers (n = 1057) were genotyped with SNPs located within and flanking these genes. One SNP located near QKI and one near GRIA2 were nominally associated with flight speed (P ≤ 0.05) although neither was significant after Bonferroni correction. Several studies have shown a correlation between flight speed and feed intake or gain; therefore, we also analyzed SNPs on BTA6:38–39 Mb known to be associated with average daily gain (ADG) and average daily feed intake (ADFI) for association with flight speed. Several SNPs on BTA6 were associated with flight speed (P ≤ 0.005), and three were significant after Bonferroni correction. These results suggest that the genes tested are unlikely to contribute to flight speed variation for our cattle population, but SNPs on BTA6 associated with ADG and ADFI may influence temperament. Use of these markers to select for economically important feed efficiency phenotypes may produce cattle with more desirable temperaments.     

Boron fertilizers in rape – a risk for honey bees?

Rape (Brassica napus L.) is foraged intensively by honey bees (Apis mellifera). Pesticide applications during bloom are sometimes combined with foliar boron fertilizer applications. Boron has insecticidal properties, and therefore, risk to honey bees cannot be excluded. This study was conducted to test whether foliar boron fertilizers could be hazardous for bees under real field conditions. Six colonies were transferred to a rape field in bloom which was treated with boron (1 kg/ha). Six control colonies were transferred to an untreated rape field approximately 7 km away. Performance parameters of the colonies were measured. Samples of honey and beebread were collected from all colonies before and after boron fertilizer application. The contents of boron and of Al, Cd, Cr, Fe, K, Mn, Ni, P, Pb, S and Zn were measured in honey by inductively coupled plasma mass spectroscopy (ICP MS) and by ICP–atomic emission spectroscopy (ICP‐OES). No significant differences were found in honey yield (P = 0.622), number of capped brood (P = 0.089), number of uncapped brood (P = 0.123) or number of bees (P = 0.87). Application of boron fertilizer did not affect the concentration of boron in honey (P = 0.656) or beebread (P = 0.665). The concentrations of other elements confirmed the suitability of rape nectar for bee nutrition. This study suggests that the application of foliar boron fertilizers in blooming rape is not hazardous for bee colonies.