A genome‐wide association study suggests new candidate genes for milk production traits in Chinese Holstein cattle

A genome‐wide association study (GWAS) was conducted on 15 milk production traits in Chinese Holstein. The experimental population consisted of 445 cattle, each genotyped by the GGP (GeneSeek genomic profiling)‐BovineLD V3 SNP chip, which had 26 151 public SNPs in its manifest file. After data cleaning, 20 326 SNPs were retained for the GWAS. The phenotypes were estimated breeding values of traits, provided by a public dairy herd improvement program center that had been collected once a month for 3 years. Two statistical models, a fixed‐effect linear regression model and a mixed‐effect linear model, were used to estimate the association effects of SNPs on each of the phenotypes. Genome‐wide significant and suggestive thresholds were set at 2.46E‐06 and 4.95E‐05 respectively. The two statistical models concurrently identified two genome‐wide significant (P < 0.05) SNPs on milk production traits in this Chinese Holstein population. The positional candidate genes, which were the ones closest to these two identified SNPs, were EEF2K (eukaryotic elongation factor 2 kinase) and KLHL1 (kelch like family member 1). These two genes could serve as new candidate genes for milk yield and lactation persistence, yet their roles need to be verified in further function studies.     

Characterization of the porcine acyl-CoA synthetase long-chain 4 gene and its association with growth and meat quality traits

Summary Long-chain acyl-CoA synthetase (ACSL) catalyses the formation of long-chain acyl-CoA from fatty acid, ATP and CoA, activating fatty acids for subsequent reactions. Long-chain acyl-CoA synthetase thus plays an essential role in both lipid biosynthesis and fatty acid degradation. The ACSL4 gene was evaluated as a positional candidate gene for the quantitative trait loci (QTL) located between SW2456 and SW1943 on chromosome X. We have sequenced 4906 bp of the pig ACSL4 mRNA. Sequence analysis allowed us to identify 10 polymorphisms located in the 3'-UTR region and to elucidate two ACSL4 haplotypes. Furthermore, a QTL and an association study between polymorphisms of the ACSL4 gene and traits of interest were carried out in an Iberian x Landrace cross. We report QTL that have not been previously identified, and we describe an association of the ACSL4 polymorphisms with growth and percentage of oleic fatty acid. Finally, we have determined allelic frequencies in 140 pigs belonging to the Iberian, Landrace, Large White, Meishan, Pietrain, Duroc, Vietnamese, Peccary and Babirusa populations.     

Red and blond Mangalitza pigs display a signature of divergent directional selection in the SLC45A2 gene

The Mangalitza lard‐type pig breed is well known for its fat appearance and curly hair, and it is mainly distributed in Eastern Europe. Four main lines were created in the nineteenth century by artificial selection: Blond Mangalitza, Black Mangalitza, Swallow‐Belly Mangalitza and Red Mangalitza. The Swallow‐Belly line has a black coat combined with yellow‐blond throat and underbelly. In the current work, we aimed to investigate if the colourations of Mangalitza pigs are genetically determined by one or a few loci whose frequencies have been modified by artificial selection. The results of selection scans, with Hap FLK and BayeScan , and of a GWAS for coat colour highlighted the existence of one region on SSC16 (18–20 Mb) with potential effects on hair pigmentation (Red vs. Blond contrast). The analysis of the gene content of this region allowed us to detect the solute carrier family 45 member 2 (SLC45A2) locus as a candidate gene for this trait. The polymorphism of the SLC45A2 locus has been associated with reduced levels or the absence of melanin in several mammalian species. The genotyping of four missense polymorphisms evidenced that rs341599992:G > A and rs693695020:G > A SNPs are strongly but not fully associated with the red and blond coat colours of Mangalitza pigs, a result that was confirmed by performing a haplotype association test. The near fixation of alternative SLC45A2 genotypes in Red and Blond Mangalitza pigs provides a compelling example of the consequences of a divergent directional selection for coat colour in a domestic species.     

Genome‐wide association QTL mapping for teat number in a purebred population of Duroc pigs

Because of increasing litter size in Western pig breeds, additional teats are desirable to increase the capacity for nursing offspring. We applied genome‐wide SNP markers to detect QTL regions that affect teat number in a Duroc population. We phenotyped 1024 animals for total teat number. A total of 36 588 SNPs on autosomes were used in the analysis. The estimated heritability for teat number was 0.34 ± 0.05 on the basis of a genomic relationship matrix constructed from all SNP markers. Using a BayesC method, we identified a total of 18 QTL regions that affected teat number in Duroc pigs; 9 of the 18 regions were newly detected.     

Systematic analyses for candidate genes of milk production traits in water buffalo (Bubalus Bubalis)

Water buffalo (Bubalus bubalis) is of great economic importance as a provider of milk and meat in many countries. However, the milk yield of buffalo is much lower than that of Holstein cows. Selection of candidate genes related to milk production traits can be applied to improve buffalo milk performance. A systematic review of studies of these candidate genes will be greatly beneficial for researchers to timely and efficiently understand the research development of molecular markers for buffalo milk production traits. Here, we identified and classified the candidate genes associated with buffalo milk production traits. A total of 517 candidate genes have been identified as being associated with milk performance in different buffalo breeds. Nineteen candidate genes containing 47 mutation sites have been identified using the candidate gene approach. In addition, 499 candidate genes have been identified in six genome‐wide association studies (GWASes) including two studies performed with the bovine SNP chip and four studies with the buffalo SNP chip. Genes CTNND2 (catenin delta 2), APOB (apolipoprotein B), FHIT (fragile histidine triad) and ESRRG (estrogen related receptor gamma) were identified in at least two GWASes. These four genes, especially APOB, deserve further study to explore regulatory roles in buffalo milk production. With growth in the number of buffalo genomic studies, more candidate genes associated with buffalo milk production traits will be identified. Therefore, future studies, such as those investigating gene location and functional analyses, are necessary to facilitate the exploitation of genetic potential and the improvement of buffalo milk performance.     

Role of ACSL4 in the chemical-induced cell death in human proximal tubule epithelial HK-2 cells

Acyl-CoA synthetase long-chain family member 4 (ACSL4) activates polyunsaturated fatty acids (PUFAs) to produce PUFA-derived acyl-CoAs, which are utilised for the synthesis of various biological components, including phospholipids (PLs). Although the roles of ACSL4 in non-apoptotic programmed cell death ferroptosis are well-characterised, its role in the other types of cell death is not fully understood. In the present study, we investigated the effects of ACSL4 knockdown on the levels of acyl-CoA, PL, and ferroptosis in the human normal kidney proximal tubule epithelial (HK-2) cells. Liquid chromatography–tandem mass spectrometry (LC-MS/MS) analyses revealed that the knockdown of ACSL4 markedly reduced the levels of PUFA-derived acyl-CoA, but not those of other acyl-CoAs. In contrast with acyl-CoA levels, the docosahexaenoic acid (DHA)-containing PL levels were preferentially decreased in the ACSL4-knockdown cells compared with the control cells. Cell death induced by the ferroptosis inducers RSL3 and FIN56 was significantly suppressed by treatment with ferrostatin-1 or ACSL4 knockdown, and, unexpectedly, upon treating with a necroptosis inhibitor. In contrast, ACSL4 knockdown failed to suppress the other oxidative stress-induced cell deaths initiated by cadmium chloride and sodium arsenite. In conclusion, ACSL4 is involved in the biosynthesis of DHA-containing PLs in HK-2 cells and is specifically involved in the cell death induced by ferroptosis inducers.     

Genome‐wide association study reveals candidate genes associated with body measurement traits in Chinese Wagyu beef cattle

We performed a genome‐wide association study to identify candidate genes for body measurement traits in 463 Wagyu beef cattle typed with the Illumina Bovine HD 770K SNP array. At the genome‐wide level, we detected 18, five and one SNPs associated with hip height, body height and body length respectively. In total, these SNPs are within or near 11 genes, six of which (PENK, XKR4, IMPAD1, PLAG1, CCND2 and SNTG1) have been reported previously and five of which (CSMD3, LAP3, SYN3, FAM19A5 and TIMP3) are novel candidate genes that we found to be associated with body measurement traits. Further exploration of these candidate genes will facilitate genetic improvement in Chinese Wagyu beef cattle.     

Genome‐wide association studies for hematological traits in swine

Improving immune capacity may increase the profitability of animal production if it enables animals to better cope with infections. Hematological traits play pivotal roles in animal immune capacity and disease resistance. Thus far, few studies have been conducted using a high‐density swine SNP chip panel to unravel the genetic mechanism of the immune capability in domestic animals. In this study, using mixed model‐based single‐locus regression analyses, we carried out genome‐wide association studies, using the Porcine SNP60 BeadChip, for immune responses in piglets for 18 hematological traits (seven leukocyte traits, seven erythrocyte traits, and four platelet traits) after being immunized with classical swine fever vaccine. After adjusting for multiple testing based on permutations, 10, 24, and 77 chromosome‐wise significant SNPs were identified for the leukocyte traits, erythrocyte traits, and platelet traits respectively, of which 10 reached genome‐wise significance level. Among the 53 SNPs for mean platelet volume, 29 are located in a linkage disequilibrium block between 32.77 and 40.59 Mb on SSC6. Four genes of interest are located within the block, providing genetic evidence that this genomic segment may be considered a candidate region relevant to the platelet traits. Other candidate genes of interest for red blood cell, hemoglobin, and red blood cell volume distribution width also have been found near the significant SNPs. Our genome‐wide association study provides a list of significant SNPs and candidate genes that offer valuable information for future dissection of molecular mechanisms regulating hematological traits.     

Malic enzyme 1 genotype is associated with backfat thickness and meat quality traits in pigs

Malic enzyme 1 (ME1) is a part of the tricarboxylate shuttle that provides NADPH and acetyl-CoA required in fatty acid biosynthesis. The pig ME1 locus maps on the proximal end of chromosome 1, where a quantitative trait loci (QTL) affecting fat deposition has been previously described. We amplified fragments of 1457 and 1459 bp that corresponded to the complete coding region and the 3'-untranslated region (UTR), respectively, of the pig ME1 gene. The sequences of these two fragments in pigs from three breeds (Landrace, Large White and Piétrain) contained five single nucleotide polymorphisms (SNP) in the 3'-UTR: C1706T, G1762T, A1807C, C1857A and T1880A. Three haplotypes were found in two generations of a selected Landrace population: H1 (C1706 G1762 A1807 C1857 A1880), H2 (C1706 G1762 A1807 C1857 T1880) and H3 (T1706 T1762 C1807 A1857 T1880). Using Bayesian association analyses, significant associations (highest posterior density at 95%) between ME1 genotype and backfat (BF) thickness at 171 days and muscular pH were found in a Landrace population.     

Whole‐genome association study for fatty acid composition of oleic acid in Japanese Black cattle

Fatty acid composition, especially oleic acid (C18:1), plays an important role in the eating quality of meat in Japanese Black cattle. Therefore, the objective of this study was to identify loci associated with C18:1 in the intramuscular fat of the trapezius muscles in Japanese Black cattle using the Illumina BovineSNP50 BeadChip whole genome single nucleotide polymorphism (SNP) assay. We also evaluated the relationship between C18:1 and three fatty acid synthesis genes, fatty acid synthase (FASN), stearoyl‐CoA desaturase and sterol regulatory element‐binding protein‐1. In this experiment, we applied a mixed model and Genomic Control approach using selective genotyping to perform a genome‐wide association study. A total of 160 animals (80 animals with higher values and 80 animals with lower values), selected from 3356 animals based on corrected phenotype, were genotyped using the Illumina BovineSNP50 BeadChip and three fatty acid synthesis genes, and the quality of these SNPs was assessed. In this study, a total of 38 955 SNPs, which included SNPs in the three fatty acid synthesis genes, were used, and the estimated inflation factor was 1.06. In the studied population, a total of 32 SNPs, including the FASN gene, had significant effects, and in particular 30 SNPs of all significant SNPs were located between 49 and 55 Mbp on chromosome 19. This study is one of the first genome‐wide association studies for fatty acid composition in a cattle population using the recently released Illumina BovineSNP50 BeadChip.     

丙酸通路基因多态性与猪肉质及胴体性状的关联分析

为了挖掘新的猪肉品质及胴体性状的候选基因,揭示猪肉质及胴体性状的遗传机制,文章将丙酸代谢通路作为候选通路,将通路内基因与猪肉质及胴体性状进行关联分析。实验采用37头三元杂交商品猪作为研究对象,首次针对丙酸通路中7个基因的36个SNP位点利用SNaPshot方法进行基因分型,分别用最小二乘模型及MB-MDR模型与肉质及胴体性状进行关联分析。结果发现,基因PCCB、MUT、MCEE及ACSS2上的4个SNP位点分别与肌内脂肪含量、背膘厚等性状显著相关(P<0.05),ACSS2与猪脂肪含量显著相关;MCEE及MUT与猪的背膘厚显著相关;PCCB基因与脂重显著相关。通过MB-MDR方法检测到多个SNP位点具有互作效应,并与背膘厚、水分含量、脂肪含量显著相关(P<0.05)。另外,丙酸代谢通路中的基因间的互作效应对猪肉品质有显著影响。     

An Essential Role of the N-Terminal Region of ACSL1 in Linking Free Fatty Acids to Mitochondrial β-Oxidation in C2C12 Myotubes

Free fatty acids are converted to acyl-CoA by long-chain acyl-CoA synthetases (ACSLs) before entering into metabolic pathways for lipid biosynthesis or degradation. ACSL family members have highly conserved amino acid sequences except for their N-terminal regions. Several reports have shown that ACSL1, among the ACSLs, is located in mitochondria and mainly leads fatty acids to the β-oxidation pathway in various cell types. In this study, we investigated how ACSL1 was localized in mitochondria and whether ACSL1 overexpression affected fatty acid oxidation (FAO) rates in C2C12 myotubes. We generated an ACSL1 mutant in which the N-terminal 100 amino acids were deleted and compared its localization and function with those of the ACSL1 wild type. We found that ACSL1 adjoined the outer membrane of mitochondria through interaction of its N-terminal region with carnitine palmitoyltransferase-1b (CPT1b) in C2C12 myotubes. In addition, overexpressed ACSL1, but not the ACSL1 mutant, increased FAO, and ameliorated palmitate-induced insulin resistance in C2C12 myotubes. These results suggested that targeting of ACSL1 to mitochondria is essential in increasing FAO in myotubes, which can reduce insulin resistance in obesity and related metabolic disorders.     

Genome-wide associations for investigating time-dependent genetic effects for milk production traits in dairy cattle

Phenotypic variation in milk production traits has been described over the course of a lactation as well as between different parities. The objective of this study was to investigate whether variation in production is affected by different loci across lactations. A genome-wide association study (GWAS) using a 50-k SNP chip was conducted in 152 divergent German Holstein Friesian cows to test for association with milk production traits over different lactations. The first four lactations were analysed regarding milk yield, fat, protein, lactose, milk urea nitrogen yield and content as well as somatic cell score. Two approaches were used: (i) Wilmink curve parameters were used to assess the genetic effects over the course of a lactation and (ii) test-day yield deviations (YD) were used as a normative approach for a GWAS. The significant effects were largest for markers affecting curve parameters for which there was a statistical power <0.8 of detection even in this small design. While significant markers for YDs were detected in this study, the power to detect effects of a similar magnitude was only 0.11, suggesting that many loci may have been missed with this approach in the present design. Furthermore, all significant effects were specific for a single lactation, leading to the conclusion that the variance explained by a certain locus changes from lactation to lactation. We confirm the common evidence that most production traits vary in the degree of persistency after the peak as a result of genetic influence.     

Peroxisomal fission is induced during appressorium formation and is required for full virulence of the rice blast fungus

Peroxisomes are involved in various metabolic processes and are important for virulence in different pathogenic fungi. How peroxisomes rapidly emerge in the appressorium during fungal infection is poorly understood. Here, we describe a gene, PEF1, which can regulate peroxisome formation in the appressorium by controlling peroxisomal fission, and is required for plant infection in the rice blast fungus Magnaporthe oryzae. Targeted deletion of PEF1 resulted in a reduction in virulence and a delay in penetration and invasive growth in host cells. PEF1 was particularly expressed during appressorial development, and its encoding protein was co‐localized with peroxisomes during appressorial development. Compared with the massive vesicle‐shaped peroxisomes formed in the wild‐type appressorium, the Δpef1 mutant could only form stringy linked immature peroxisomes, suggesting that PEF1 was involved in peroxisomal fission during appressorium formation. We also found that the Δpef1 mutant could not utilize fatty acids efficiently, which can improve significantly the expression level of PEF1 and induce peroxisomal fission. As expected, the Δpef1 mutant showed reduced intracellular production of reactive oxygen species (ROS) during appressorium formation and induced ROS accumulation in host cells during infection. Taken together, PEF1‐mediated peroxisomal fission is important for fungal infection by controlling the number of peroxisomes in the appressorium.     

Genomic analysis of the differential response to experimental infection with porcine circovirus 2b

Porcine circovirus type 2 (PCV2) is the etiological agent of a group of associated diseases (PCVAD) that affect production efficiency and can lead to mortality. Using different crossbred lines of pigs, we analyzed host genetic variation of viral load, immune response and weight change following experimental infection with a PCV2b strain (n = 386). Pigs expressed variation in the magnitude and initiation of viremia and immune response recorded weekly until 28 days post‐infection. A higher viral load was correlated with weight gain (r = ?0.26, P < 0.0001) and presence of PCV2‐specific antibodies (IgM, r = 0.26–0.34, P < 0.0001; IgG, r = 0.17–0.20, P < 0.01). In genome‐wide association analyses of the responses at different time points, the proportions of phenotypic variation explained by combined effects of 56 433 SNPs were 34.8–59.4% for viremia, 10.1–59.5% for antibody response and 5.6–14.9% for weight change. Relationships between genomic prediction of overall viral load and weight gain during the first weeks of challenge were negative (?0.21 and ?0.24 respectively, P < 0.0001). Individuals that carried more favorable alleles across three SNPs on SSC9 (0.60 Mb) and SSC12 (6.8 and 18.2 Mb) partially explained this relationship, having lower viral load (P < 0.0001); lower viremia at day 14 (P < 0.0001), day 21 (P < 0.01) and day 28 (P < 0.05) and greater overall average daily gain during infection (ADG_i; P < 0.01), ADG_i at week 3 (P < 0.001) and week 4 (P < 0.01). These additive genetic relationships could lead to molecular solutions to improve animal health and reduce production costs.     

Genome‐wide association study of six quality traits reveals the association of the TaRPP13L1 gene with flour colour in Chinese bread wheat

Flour colour, kernel hardness, grain protein content and wet gluten content are important quality properties that determine end use in bread wheat. Here, a wheat 90K genotyping assay was used for a genome‐wide association study (GWAS) of the six quality‐related traits in Chinese wheat cultivars in eight environments over four years. A total of 846 significant single nucleotide polymorphisms (SNPs) were identified, explaining approximately 30% of the phenotypic variation on average, and 103 multienvironment‐significant SNPs were detected in more than four environments. Quantitative trait loci (QTL) mapping in the biparent population confirmed some important SNP loci. Moreover, it was determined that some important genes were associated with the six quality traits, including some known functional genes and annotated unknown functional genes. Of the annotated unknown functional genes, it was verified that TaRPP13L1 was associated with flour colour. Wheat cultivars or lines with TaRPP13L1‐B1a showed extremely significantly higher flour redness and lower yellowness than those with TaRPP13L1‐B1b in the Chinese wheat natural population and the doubled haploid (DH) population. Two tetraploid wheat lines with premature stop codons of the TaRPP13L1 gene mutagenized by ethyl methanesulfonate (EMS) showed extremely significantly higher flour redness and lower yellowness than wild type. Our data suggest that the TaRPP13L1 gene plays an important role in modulating wheat flour colour. This study provides useful information for further dissection of the genetic basis of flour colour and also provides valuable genes or genetic loci for marker‐assisted selection to improve the process of breeding quality wheat in China.     

Evaluation of tubulin β‐3 as a novel senescence‐associated gene in melanocytic malignant transformation

Malignant melanoma might develop from melanocytic nevi in which the growth‐arrested state has been broken. We analyzed the gene expression of young and senescent human melanocytes in culture and compared the gene expression data with a dataset from nevi and melanomas. A concordant altered gene expression was identified in 84 genes when comparing the growth‐arrested samples with proliferating samples. TUBB3, which encodes the microtubule protein tubulin β‐3, showed a decreased expression in senescent melanocytes and nevi and was selected for further studies. Depletion of tubulin β‐3 caused accumulation of cells in the G2/M phase and decreased proliferation and migration. Immunohistochemical assessment of tubulin β‐3 in benign lesions revealed strong staining in the superficial part of the intradermal components, which faded with depth. In contrast, primary melanomas exhibited staining without gradient in a disordered pattern and strong staining of the invasive front. Our results describe an approach to find clinically useful diagnostic biomarkers to more precisely identify cutaneous malignant melanoma and present tubulin β‐3 as a candidate marker.     

Genetic mapping of recurrent exertional rhabdomyolysis in a population of North American Thoroughbreds

Recurrent exertional rhabdomyolysis is a heritable disorder that results in painful skeletal muscle cramping with exercise in up to 10% of all Thoroughbred racehorses. Here, we report a genome‐wide association study with 48 282 SNPs analyzed among 48 case and 37 control Thoroughbreds. The most significant SNPs spanned approximately 13 Mb on ECA16, and the P‐value of the most significant SNP after correcting for population structure was 8.0 × 10^?6. This region on ECA16 was further evaluated by genotyping 247 SNPs in both the initial population and a second population of 34 case and 98 control Thoroughbreds. Several SNPs across the 13‐Mb region on ECA16 showed significance when each population was analyzed separately; however, the exact positions of the most significant SNPs within this region on ECA16 varied between populations. This variability in location may be attributed to lack of power owing to insufficient sample sizes within each population individually, or to the relative distribution of long, conserved haplotypes, characteristic of the Thoroughbred breed. Future genome‐wide association studies with additional horses would likely improve the power to resolve casual loci located on ECA16 and increase the likelihood of detecting any additional loci on other chromosomes contributing to disease susceptibility.