Next Article in Journal
POLE-Mutant Colon Cancer Treated with PD-1 Blockade Showing Clearance of Circulating Tumor DNA and Prolonged Disease-Free Interval
Previous Article in Journal
Temperature-Related Short-Term Succession Events of Bacterial Phylotypes in Potter Cove, Antarctica
Previous Article in Special Issue
Genome-Wide Scan of Wool Production Traits in Akkaraman Sheep
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Genes
Volume 14
Issue 5
10.3390/genes14051053
genes-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Single-Locus and Multi-Locus Genome-Wide Association Studies Identify Genes Associated with Liver Cu Concentration in Merinoland Sheep

by
Olusegun O. Adeniyi
1,
Ivica Medugorac
2,
Ewa Grochowska
3,
Rolf-Alexander Düring
4 and
Gesine Lühken
1,*
1
Institute of Animal Breeding and Genetics, Justus Liebig University Giessen, Ludwigstrasse 21, 35390 Giessen, Germany
2
Population Genomics Group, Department of Veterinary Sciences, Ludwig Maximilian University Munich, Lena-Christ-Str. 48, 82152 Martinsried, Germany
3
Department of Animal Biotechnology and Genetics, Bydgoszcz University of Science and Technology, Mazowiecka 28 St., 85-084 Bydgoszcz, Poland
4
Institute of Soil Science and Soil Conservation, Interdisciplinary Research Center for Biosystems, Land Use and Nutrition (IFZ), Justus Liebig University Giessen, Heinrich-Buff-Ring 26-32, 35392 Giessen, Germany
*
Author to whom correspondence should be addressed.
Genes 2023, 14(5), 1053; https://doi.org/10.3390/genes14051053
Submission received: 30 March 2023 / Revised: 2 May 2023 / Accepted: 4 May 2023 / Published: 8 May 2023
(This article belongs to the Special Issue Genetics and Breeding of Small Ruminants)
Browse Figures
Versions Notes

Abstract

:
Economic losses due to copper intoxication or deficiency is a problem encountered by sheep farmers. The aim of this study was to investigate the ovine genome for genomic regions and candidate genes responsible for variability in liver copper concentration. Liver samples were collected from slaughtered lambs of the Merinoland breed from two farms, and used for measurement of copper concentration and genome-wide association study (GWAS). A total of 45,511 SNPs and 130 samples were finally used for analysis, in which single-locus and several multi-locus GWAS (SL-GWAS; ML-GWAS) methods were employed. Gene enrichment analysis was performed for identified candidate genes to detect gene ontology (GO) terms significantly associated with hepatic copper levels. The SL-GWAS and a minimum of two ML-GWAS identified two and thirteen significant SNPs, respectively. Within genomic regions surrounding identified SNPs, we observed nine promising candidate genes such as DYNC1I2, VPS35, SLC38A9 and CHMP1A. GO terms such as lysosomal membrane, mitochondrial inner membrane and sodium:proton antiporter activity were significantly enriched. Genes involved in these identified GO terms mediate multivesicular body (MVB) fusion with lysosome for degradation and control mitochondrial membrane permeability. This reveals the polygenic status of this trait and candidate genes for further studies on breeding for copper tolerance in sheep.

1. Introduction

Sheep farming is continuously challenged with the risks of copper intoxication and deficiency. Copper is an essential trace element for many species. It serves as a cofactor for the proper functioning of some copper-containing enzymes, so-called ”cuproenzymes” such as copper-–zinc superoxide dismutase, cytochrome c oxidase and hephaestin which function in the anti-oxidation of free radicals, energy metabolism and iron metabolism, respectively [1,2]. However, sheep are known to be impaired in their ability to excrete excess copper from the liver, which leads to copper intoxication [3]. On the other hand, copper deficiency is also a known problem in sheep husbandry [4].
Hepatic copper levels in sheep have been classified into diagnostic levels of deficient, marginal, adequate, high, and toxic [5,6]. Deficient and high diagnostic levels of hepatic copper have been associated with poor performance and health status of sheep [7,8]. Liver copper levels have been reported to vary between as well as within sheep breeds [9,10]. A report on Suffolk and Texel sheep fed a similar diet showed that the coefficient of variation for liver copper concentration within each breed was 40% and 34%, respectively [9]. Furthermore, results of selection for low and high copper status in Scottish Blackface lambs have been reported by Woolliams et al. [11]. According to these authors, plasma copper concentration was higher in high copper lines than in low copper lines. In another report by Knowles et al. [12], variation in copper status was observed for two flocks of the Romney breed suggesting the availability of lines with distinct copper metabolism. These findings suggest that copper status within sheep breeds may be influenced by genetic factors.
Association of genes such as copper transporter 1 (CTR1), chaperone for superoxide dismutase 1 (CCS), antioxidant 1 (ATOX1), metallothionein (MT), copper-transporting ATPase 1 (ATP7A) and copper-transporting ATPase 2 (ATP7B) with copper metabolism are well documented for humans [13,14,15,16,17]. CTR1, CCS and ATOX1 are proteins involved in copper transport from circulation, distribution to superoxide dismutase (SOD1) and transport to the secretory pathway and nucleus, respectively [14]. In addition, ATP7A/B supports the transfer of copper to the secretory pathways and mediates copper excretion by sequestration in vesicles/vacuoles, while MT functions as a copper scavenger in hepatocytes [14,15,17]. Variations in liver copper levels have been reported for humans with mutations in the ATP7A and/or ATP7B genes that result in Menkes and Wilson disease, respectively [18,19]. Likewise, mutation in the copper metabolism MURR1 domain-containing 1 (COMMD1) gene has been implicated in copper toxicity in dogs [2]. However, information on candidate genes and genetic markers associated with within-breed variation in liver copper concentration of sheep is unknown.
The Illumina Ovine 50K SNP chip is a powerful tool to analyze the sheep genome with respect to an association with diverse traits as confirmed in various single-locus genome-wide association studies (SL-GWAS) [20,21,22,23] and multi-locus GWAS (ML-GWAS) [24,25] in animals. SL-GWAS aim to identify associations of single nucleotide polymorphisms (SNPs) with varying phenotypes of individuals by assessing the differences in allelic frequencies of their genetic variants [26,27]. This is currently implemented in various mixed linear model methods such as MLMA [28], EMMAX [29] and GEMMA [30]. These methods are usually subjected to Bonferroni correction due to multiple testing, which may result in the exclusion of some important loci with small effects. In this context, ML-GWAS is a better alternative because it does not require the conservative Bonferroni correction, therefore leading to the identification of more trait-associated markers especially for polygenic traits [25,31]. Moreover, an earlier report indicated that a multi-locus GWAS approach is better suited for the analysis of complex traits [32]. As a solution, several ML-GWAS models have been developed, such as multi-locus random-SNP-effect mixed linear model (mrMLM) [33], the fast multi-locus random-SNP-effect mixed linear model (FASTmrMLM) [34], the fast multi-locus random-SNP-effect efficient mixed-model association (FASTmrEMMA) [35], polygenic-background-control-based least angle regression plus empirical Bayes (pLARmEB) [36], polygenic-background-control-based Kruskal–Wallis test plus empirical Bayes (pKWmEB) [37], and Iterative sure independence screening expectation maximization Bayesian least absolute shrinkage and selection operator (ISIS EM-BLASSO) [38]. These methods are implemented in two phases. The first phase involves the selection of potentially associated markers after analyses using various algorithms. For the second phase, all the effects in a model harbouring the selected markers are estimated by empirical Bayes, after which further identification by a likelihood ratio test for true quantitative trait nucleotides (QTNs) of all the non-zero effects is performed.
To the best of our knowledge, a genome-wide analysis regarding liver copper concentration in sheep has not been conducted until now. Therefore, the aim of this study was to identify candidate genes associated with variability in hepatic copper concentration in sheep using both SL-GWAS and ML-GWAS methods. This analysis was performed using liver samples from slaughtered lambs of the Merinoland sheep breed, a breed of high economic importance in Germany.

2. Materials and Methods

2.1. Sample and Data Collection

The animals were kept under standard farming regulations on two private sheep farms in Bavaria. For this study, only male lambs were sampled after slaughter for human consumption. Documented information on the date of birth and date of slaughter was obtained from the farmers. Regarding farm 1, 89 lambs from the flock with ages between 4 and 5.5 months were sampled after slaughter in batches over a period of 46 days. They were kept indoors with unlimited access to their dams and fed ad libitum with grass silage and compound feed consisting of lucerne pellets (5.53%), concentrate pellets (58.84%), barley (22.93%) and beet pulp (12.70%). Concentrate feeding of lambs started 14 days after birth. All dams had unhindered access to grass silage with occasional input of the above-mentioned compound feed. A total of 45 lambs were sampled from the flock located at farm 2. These lambs were kept on pasture and under uniform feeding conditions until slaughter. In addition, fresh water was made available ad libitum for all animals by the farmers. The nutrient composition of the feed, and the mineral content of the grass silage, components of compound feed, and pasture grass were determined by a nutrient analysis laboratory (Intertek Food Services GmbH, Linden, Germany) (Table S1a,b). Liver samples from all lambs were collected directly after slaughter by cutting approximately 20 g of sheep liver from the tip of the Lobus caudatus, packed separately in plastic tubes, labelled and stored at −20 °C prior to the determination of liver copper concentration and DNA extraction.

2.2. Determination of Liver Copper Concentration

Liver samples were freeze-dried, pulverized and stored in airtight tubes. Duplicate samples of approximately 0.5 g were digested in 5 mL deionized H2O, 5 mL HNO3 and 3 mL H2O2 with a StarT-1500 microwave oven (MLS GmbH Cooperation, Leutkirch, Germany). Conditions for the microwave digestion system are shown in Table S2. Digested samples were filtered into a 50 mL flask and filled up to a volume of 50 mL with deionized water. The copper concentration was determined with an inductively coupled plasma-optical emission spectrometer (ICP-OES; Agilent 720ES, Darmstadt, Germany) at a wavelength of 327.4 nm. Used operating parameters for the ICP-OES were published elsewhere [39]. For quality assurance reagent blanks were measured and a detection limit corresponding to 0.039 mg/kg was calculated (limit of detection = mean blank + 3 standard deviations of the mean blank value). Furthermore, certified reference materials (Bovine Liver ERM-BB185, European Commission, Joint Research Centre, Geel, Belgium) were used for testing the precision of measurement. The mean recovery percentage was 91.4% ± 2.6.

2.3. Genotyping and Quality Control

DNA was extracted from liver samples using the Macherey-Nagel NucleoSpin Tissue Kit (Düren, Germany) according to the manufacturer’s instructions. All 134 liver samples were genotyped with the Illumina Ovine 50k SNP BeadChip. Quality control (QC) was performed to exclude loci with minor allele frequency (<5%), SNPs with a low call rate (<95%), samples with more than 5% of missing genotypes, and SNPs for which Hardy–Weinberg equilibrium test p-values were lower than 1.0 × 10−6. After QC a total of 45,511 SNPs and 130 samples were left for analysis.

2.4. Single-Locus and Multi-Locus Genome-Wide Association Analyses

SL-GWAS and ML-GWAS were performed for all samples and data obtained after QC. All 130 samples were used to estimate the genome-wide association of markers with liver copper concentration. To perform a single-locus GWAS, genetic relatedness was estimated using “pcrelate” and “pcrelateToMatrix” functions in the GENESIS R package [40]. Next, a linear mixed model analysis for quantitative phenotypes was implemented using the same software. A genetic relationship matrix was included in the analysis as a random effect to account for cryptic relatedness. The model includes: y = X β + Z 1 a + Z 2 c + e where y is the vector of observed liver copper concentration (mg/kg dry matter, DM) at slaughter, ß is a vector of fixed effects including the intercept term, farm, age at slaughter, a is the vector of the additive effect (fixed) of the candidate SNP to be tested for association, c is the vector of polygenic effect (random) of all markers (as captured by the genetic relationship matrix calculated using all SNPs), and e is a vector of residuals. X, Z1 and Z2 are the incidence matrices relating observations to β, a and c, respectively. Due to the conservative nature of Bonferroni threshold [41] and possibility that associated SNP effects can be missed [42], a false discovery rate (FDR) p-value threshold was used [43]. FDR was set at 0.01, with the p-value level of significance defined as p = N/M × 0.01, where N is the number of SNP markers with p-values below the FDR threshold and M represents the total number of SNPs [44]. A Manhattan plot of -log10(p-value) was constructed with the package “ggplot” in R [45]. SNPs above the given p-value threshold were considered significant and selected for further analysis. The proportion of total variability explained by SNPs (VG) and error variance (VE) were computed using “varCompCI” function in GENESIS R package [40]. Estimation of heritability was performed from the variance components with the formula: h2 = VG / (VG + VE) according to Yang et al. [46].
Regarding ML-GWAS analysis, the mrMLM v4.0.2 [47] software implemented in R was used to determine SNPs associated with liver copper concentration in Merinoland sheep. Five ML-GWAS methods including mrMLM, FASTmrMLM, FASTmrEMMA, pLARmEB, and ISIS EM-BLASSO were performed according to methods described by Wang et al. [33], Tamba and Zhang [34], Wen et al. [35], Zhang et al. [36], and Tamba et al. [38], respectively. In the first stage and for the selection of potential SNPs, all SNPs were treated as random effects in the five methods. In the second stage, all selected SNPs in the first stage of each method were placed into one multi-locus model and the markers with effects above the logarithm of odds (LOD) threshold value were regarded as possible candidate SNPs associated with liver copper concentration. For all methods except pLARmEB, the critical p-values for the selection of SNP in the first stage were set at default. These values were employed in order to maximize power for SNP detection as suggested by the respective authors [33,34,35,38]. The default values were 0.01 for mrMLM, 0.01 for FASTmrMLM, 0.005 for FASTmrEMMA, and 0.01 for ISIS EM-BLASSO. With regards to pLARmEB analysis, a value of 50 was set as the number of potentially associated markers selected by least angle regression (LARS) according to Zhang et al. [36]. The critical LOD score was set to 3 in the second stage according to previous reports using ML-GWAS methods [24,25]. Covariates of farm and age at slaughter were included in the analysis, which was corrected for cryptic relatedness using K (genomic relatedness) matrix calculated using the R package mrMLM v4.0.2. To increase the accuracy of this result, SNPs identified by at least two ML-GWAS methods were selected for further analysis.

2.5. Gene Annotation and Enrichment

Genes located within 0.5 megabase (Mb) regions upstream and 0.5 Mb downstream of putative SNPs after SL-GWAS and ML-GWAS were identified according to the UCSC genome browser [48] with the selection of the Oar_v4.0 assembly in the National Center for Biotechnology Information (NCBI). Functional enrichment analysis was performed on all genes identified by both GWAS methods using the Database for Annotation, Visualization and Integrated Discovery (DAVID) software (https://david.ncifcrf.gov; accessed on 23 February 2023) [49,50,51]. The human gene annotation was selected for analysis due to a higher number of recognized genes in comparison to sheep gene annotations. The gene ontology (GO) annotation category including biological process (BP), cellular component (CC) and molecular function (MF) terms, with the parameter set to a minimum of two genes, was investigated in this study. Using the Fishers exact test, GO TERMS of enriched genes with a p-value threshold of ≤0.05 were considered significant [52]. A subset of genes with functional relationship to cellular transport and excretion using available gene annotations from GeneCards (http://www.genecards.org (accessed on 23 February 2023)) and Uniprot (http://www.uniprot.org (accessed on 23 February 2023)) databases were considered as functional candidate genes. Additionally, a second subset of all genes was considered as positional candidate genes when harbouring a significant SNP or located near such an SNP.

3. Results

3.1. Liver Copper Levels and Estimated Heritability

Descriptive statistics of hepatic copper concentration and age of slaughtered lambs whose samples were included in the SL-GWAS and ML-GWAS analyses are shown in Table 1. A heritability of 0.67 ± 0.29 was estimated for this trait.

3.2. Identification of Genomic Regions and Candidate Genes

Relating to the SL-GWAS analysis of markers associated with liver copper concentration in this study, two significant SNPs (OAR2_145591151.1 and s62875.1) located on chromosome 2 and below the given threshold (p-value = 8.77 × 10−5; −log10(p-value) = 4.06; λ = 0.98) were observed (Table 2; Figure 1). Likewise, 35 significant SNPs above the LOD score cutoff value of 3 were observed after the performance of ML-GWAS methods (Table S3). Among these, 13 significant SNPs distributed on chromosomes 1, 2, 3, 4, 6, 7, 14, 16 and 23 (Table 3; Figure 2a–e) were detected by a minimum of two ML-GWAS methods. The two SNPs observed using the SL-GWAS were also identified by ML-GWAS. Of these, only one SNP (OAR2_145591151.1) was observed by more than one ML-GWAS method. Concerning candidate genes identified with SL-GWAS, 4 positional candidate genes (Table 2) within putative SNP regions were detected. Furthermore, two of these candidate genes were identified as both positional and functional. As regards ML-GWAS and within possible SNP regions associated with liver copper concentration, we observed 21 potential candidate genes (Table 3). Among them were 9 functional candidate genes, with 4 genes identified as both positional and functional candidate genes.
Finally, our results revealed a total of 16 GO terms, for example integral component of membrane, lysosomal membrane, hydrogen ion transmembrane transport, mitochondrial inner membrane and sodium:proton antiporter activity were significantly enriched (p ≤ 0.05; Table 4).

4. Discussion

Copper homoeostasis is less efficient in sheep compared to other farm animals such as cattle and goats [3]. This impairment can more easily lead to copper intoxication or deficiency and is reported to vary within and between breeds [1,9,10]. The study of copper-related genes is well documented [14,15,16,17]. However, there is a scarcity of information on the gene or genes responsible for the control of hepatic copper levels in sheep. This study aimed at identifying possible SNP markers and candidate genes associated with liver copper variation in Merinoland sheep. Measured hepatic copper concentrations were within the marginal and normal levels expected for sheep [6]. The estimated, high heritability (h2 ± SE = 0.67 ± 0.29) confirmed that this trait is largely influenced by genetic factors. This value is in the same range as estimated in Merino sheep by Judson et al. [53] (h2 ± SE = 0.60 ± 0.32; n = 208). However, the high SE for h2 observed in our study indicates that the heritability value may be inflated due to cryptic relatedness since SNP-based heritability is reported to increase with a higher number of related samples [46]. Therefore, further work with a larger sample size of unrelated animals is needed to validate this finding.
Among the candidate genes identified in this study, three genes including dynein cytoplasmic 1 intermediate chain 2 (DYNC1I2), retromer complex component (VPS35) and solute carrier family 38 member 9 (SLC38A9) located on chromosomes 2, 14 and 16 have been reported to indirectly influence cellular transport and excretion [54,55,56]. The DYNC1I2 gene expresses a subunit of the cytoplasmic dynein motor protein which is important for microtubule-based transport towards the minus end [57]. Cytoplasmic dynein transports cargoes such as cytoskeleton filaments, endosomes and lysosomes and has been reported to bind dynactin, an important protein complex associated with dynein activity through its p150Glued subunit [58,59]. Likewise, dynactin is reported to interact with the copper P-type ATPase (ATP7B) via its subunit p62 [60]. Interestingly, lysosomal copper exocytosis has been directly associated with dynactin in a report by Polishchuk et al. [61]. Findings in this report suggest that the targeting of ATP7B to the canalicular surface in hepatocyte cells, which precedes lysosomal exocytosis, requires copper-dependent interaction with the p62 dynactin subunit. This suggests that changes in the activity of cytoplasmic dynein motor protein may determine hepatic copper levels by impacting lysosomal exocytosis. Besides VPS29 and VPS26, VPS35 is one of three core proteins of the retromer complex, involved in cargo sorting and endosomal trafficking including the homeostasis of transmembrane proteins located at the plasma membrane and within endosomes/ lysosomes [62]. The retromer complex is reported to bind sorting nexin 3 (SNX3) resulting in the cargo recognition and binding of divalent metal transporter 1 (also known as solute carrier family 11 member 2) [62,63]. DMT1 has been implicated in copper ion transport in the intestine with a reduced expression resulting in lower copper absorption [64]. In addition, the VPS35 associates with FAM21 (a subunit of the WASH complex) and aids the recruitment of the WASH complex to endosomes [65,66,67]. Earlier findings suggest that copper metabolism MURR1 domain–containing 1 (COMMD1) protein, which functions as a copper chaperone and binds to the copper-transporting P-type ATPases ATP7A/ATP7B [68,69], interacts with the WASH complex and modulates copper-dependent trafficking of ATP7A [68]. For perspective, ATP7A is an important copper transporter that regulates the absorption of copper in the intestinal tract [70]. Yeast cells are considered an ideal model for the genetic analysis of mammalian genomes [71,72]. Interestingly, a report by Sowada et al. [73] showed increased copper sensitivity of yeast cells lacking VPS35. Therefore, the interaction of VPS35 with the DMT1 and WASH complex may have an influence on liver copper concentration in Merinoland sheep. The SLC38A9 gene encodes a lysosomal amino acid transporter protein reported in the mechanistic target of rapamycin complex 1 (mTORC1) signalling [56,74,75]. The mTORC1 controls transcription factor EB (TFEB) which regulates the expression of genes involved in lysosomal degradation and biogenesis such as vacuolar H+-ATPase (V-ATPase) and metallothionein 1 (MT1) [76,77]. The MT1 protein has been identified as a chaperon involved in copper homoeostasis [78]. This indicates that the involvement of SLC38A9 in mTORC1 signalling and TFEB expression may impact liver copper levels in Merinoland sheep. Notably, TFEB expression has been associated with lysosomal copper exocytosis [61]. In this report, increased TFEB expression in hepatocyte cells resulted in increased lysosomal copper exocytosis suggesting a possible influence of TFEB expression on liver copper exocytosis.
Regarding GO analysis, our findings revealed that lysosomal membrane and midbody were enriched GO terms in this study. Genes observed in these GO terms including VPS35, SLC38A9 and charged multivesicular body protein 1 A (CHMP1A) are reportedly involved in cellular transport and excretion [62,75,79]. Generally, the endosomal sorting complex required for transport (ESCRT) is important for the sorting and delivery of ubiquitinated cargo such as membrane transporters to the lysosome for degradation [80]. The CHMP1A is the member of the ESCRT-III family required for the formation of multivesicular bodies (MVBs), which fuse with lysosomes for the degradation of its content [79,81]. Additionally, the ESCRT-III is involved in the abscission of budding vesicles which mature into MVBs [80]. These findings suggest a possible influence of ESCRT on the degradation of copper transporters and consequently liver copper concentration. Furthermore, some genes including matrix AAA peptidase subunit, paraplegin (SPG7) and ATP synthase membrane subunit f (ATP5MF) were observed as functional genes involved in the mitochondrial inner membrane, which is a significant GO term observed in this study. The SPG7 and ATP5MF genes are implicated in the opening of the mitochondrial permeability transition pore (PTP) [82,83,84]. The opening of the PTP results in loss of ATP production, mitochondrial swelling and cell death due to cytosolic calcium overload and oxidative stress [83]. Interestingly, another gene (solute carrier family 25 member 12 (SLC25A12)) identified in this GO term is also required for the calcium-dependent maintenance of ATP homoeostasis in the mitochondria [85,86]. The mitochondria are important cellular organelles implicated in copper homeostasis [87]. Earlier reports suggest that copper toxicity in cells results in swelling and rupture of the mitochondria, as well as induction of oxidation stress [3,88]. According to Bernardi et al. [89], mitochondrial swelling is associated with changes in the inner mitochondrial membrane permeability and consequent expansion of the inner membrane matrix. In the report by Haywood et al. [3], mitochondria swelling in hepatocytes from the North Ronaldsay sheep breed, which is considered highly susceptible to copper intoxication, was observed at a lower liver copper concentration when compared to a less susceptible sheep breed (Cambridge sheep). These observations suggest a possible association of hepatic copper levels with gene functions regulating mitochondrial membrane permeability. Finally, the GO terms hydrogen ion transmembrane transport and sodium:proton antiporter activity were also strongly enriched with involved genes such as solute carrier family 9 member B1 (SLC9B1) and solute carrier family 9 member B2 (SLC9B2). These genes encode Na+/H+ exchangers (NHEs) important in cytoplasmic and organelle pH [90]. The SLC9B1 gene is predominantly expressed in the testis while SLC9B2 is highly expressed in the liver and brain. In earlier studies on yeast NHE (Nhx1p), it was observed that Nhx1p localized to the trans-Golgi network compartments, late endosomes, and recycling endosomes, is critical for protein sorting and endocytic pathway [91,92]. In another report, Nhx1 was found to aid MVB fusion with lysosomes [93]. These findings suggest that NHEs may influence endosomal sorting and trafficking of important copper transporters and consequently copper concentration in hepatocytes.

5. Conclusions

The results revealed a number of identified SNPs potentially contributing to variability in hepatic copper levels in Merinoland sheep. The identified regions surrounding these SNPs harbour some promising functional candidate genes such as DYNC1I2, VPS35, SLC38A9 and CHMP1A associated with endosomal cargo sorting and trafficking, as well as lysosomal transport. Additionally, genes such as SPG7, ATP5MF and SLC25A12 were identified as functional genes involved in mitochondrial membrane permeability which has been associated with copper toxicity. Likewise, the genes SLC9B1 and SLC9B2 which are involved in luminal and intraluminal pH and MVB fusion to lysosome, were observed as potential candidate genes influencing hepatic copper levels in Merinoland sheep. In total, nine promising candidate genes were observed in this study. These genes need to be further investigated to ascertain their involvement in liver copper variation in Merinoland sheep in particular, and other sheep breeds in general. This study provides evidence for a polygenic trait with high heritability and delivers promising clues for further studies to identify potential causal variants that may be associated with variation in copper tolerance in sheep and used for practical breeding.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/genes14051053/s1, Table S1a: Nutrient composition of compound feed fed to Merinoland sheep on farm 1; Table S1b: Mineral content of grass, compound feed and pasture grass fed to Merinoland sheep on farm 1 and 2; Table S2: Operating conditions for microwave-assisted extraction of freeze-dried liver samples; Table S3: Significant SNPs associated with liver Cu concentration in Merinoland sheep detected by ML-GWAS analysis methods (LOD score > 3).

Author Contributions

Conceptualization, G.L. and O.O.A.; methodology, O.O.A. and I.M.; formal analysis, O.O.A.; resources, G.L. and R.-A.D.; writing—original draft preparation, O.O.A.; writing—review and editing, O.O.A., G.L., E.G. and I.M.; visualization, O.O.A.; supervision, G.L.; funding acquisition, G.L. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the H. Wilhelm Schaumann Stiftung and by the Deutsche Forschungsgemeinschaft (DFG, grant no. 62202147).

Data Availability Statement

Raw data are available on request from the authors.

Acknowledgments

The authors would like to thank the sheep farmers for providing animal data and samples after slaughter, and Carina Crispens and Stephanie Steitz for technical assistance.

Conflicts of Interest

The authors declare no conflict of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript; or in the decision to publish the results.

References

  1. Suttle, N.F. Mineral Nutrition of Livestock, 4th ed.; CABI: Cambridge, UK, 2010; ISBN 978-1-84593-472-9. [Google Scholar]
  2. De Bie, P.; Van De Sluis, B.; Klomp, L.; Wijmenga, C. The many faces of the copper metabolism protein MURR1/COMMD1. J. Hered. 2005, 96, 803–811. [Google Scholar] [CrossRef] [PubMed]
  3. Haywood, S.; Simpson, D.M.; Ross, G.; Beynon, R.J. The greater susceptibility of North Ronaldsay sheep compared with Cambridge sheep to copper-induced oxidative stress, mitochondrial damage and hepatic stellate cell activation. J. Comp. Pathol. 2005, 133, 114–127. [Google Scholar] [CrossRef] [PubMed]
  4. Suttle, N.F.; Jones, D.G. Micronutrient imbalance. In Diseases of Sheep, 4th ed.; Aitken, I.D., Ed.; Blackwell: Oxford, UK, 2007; pp. 377–394. [Google Scholar]
  5. Underwood, E.J.; Suttle, N.F. The Mineral Nutrition of Livestock, 3rd ed.; CABI: Cambridge, UK, 1999. [Google Scholar]
  6. Puls, R. Mineral Levels in Animal Health: Diagnostic Data; Sherpa International: Clearbrook, BC, Canada, 1994; ISBN 9780969342922. [Google Scholar]
  7. Kumaratilake, J.S. Chronic Copper Poisoning in Sheep: Liver Injury. JTEA 2014, 3, 1–22. [Google Scholar] [CrossRef]
  8. Sousa, I.K.F.d.; Hamad Minervino, A.H.; Sousa, R.D.S.; Chaves, D.F.; Soares, H.S.; Barros, I.d.O.; Araújo, C.A.S.C.d.; Júnior, R.A.B.; Ortolani, E.L. Copper deficiency in sheep with high liver iron accumulation. Vet. Med. Int. 2012, 2012, 207950. [Google Scholar] [CrossRef]
  9. Suttle, N.F.; Lewis, R.M.; Small, J.N.W. Effects of breed and family on rate of copper accretion in the liver of purebred Charollais, Suffolk and Texel lambs. Anim. Sci. 2002, 75, 295–302. [Google Scholar] [CrossRef]
  10. van der Berg, R.; Levels, F.H.; van der Schee, W. Breed differences in sheep with respect to the accumulation of copper in the liver. Vet. Q. 1983, 5, 26–31. [Google Scholar] [CrossRef]
  11. Woolliams, J.A.; Woolliams, C.; Suttle, N.F.; Jones, D.G.; Wiener, G. Studies on lambs from lines genetically selected for low and high copper status 2. Incidence of hypocuprosis on improved hill pasture. Anim. Sci. 1986, 43, 303–317. [Google Scholar] [CrossRef]
  12. Knowles, S.O.; Rounce, J.R.; Grace, N.D.; Lee, J.H. Variation in copper metabolism between two flocks of Romney sheep in response to increasing dietary copper. Proc. New Zealand Soc. Anim. Prod. 1998, 58, 195–198. [Google Scholar]
  13. Kaler, S.G. Inborn errors of copper metabolism. Handb. Clin. Neurol. 2013, 113, 1745–1754. [Google Scholar] [CrossRef] [PubMed]
  14. Lutsenko, S. Human copper homeostasis: A network of interconnected pathways. Curr. Opin. Chem. Biol. 2010, 14, 211–217. [Google Scholar] [CrossRef] [PubMed]
  15. Barry, A.N.; Shinde, U.; Lutsenko, S. Structural organization of human Cu-transporting ATPases: Learning from building blocks. J. Biol. Inorg. Chem. 2010, 15, 47–59. [Google Scholar] [CrossRef] [PubMed]
  16. Lutsenko, S.; Barnes, N.L.; Bartee, M.Y.; Dmitriev, O.Y. Function and regulation of human copper-transporting ATPases. Physiol. Rev. 2007, 87, 1011–1046. [Google Scholar] [CrossRef] [PubMed]
  17. Tao, T.Y.; Gitlin, J.D. Hepatic copper metabolism: Insights from genetic disease. Hepatology 2003, 37, 1241–1247. [Google Scholar] [CrossRef] [PubMed]
  18. Horn, N.; Wittung-Stafshede, P. ATP7A-Regulated Enzyme Metalation and Trafficking in the Menkes Disease Puzzle. Biomedicines 2021, 9, 391. [Google Scholar] [CrossRef] [PubMed]
  19. Fuentealba, C.; Aburto, E.M. Review: Animal models of copper-associated liver disease. Comp. Hepatol. 2003, 2, 5–16. [Google Scholar] [CrossRef] [PubMed]
  20. Yilmaz, O.; Kizilaslan, M.; Arzik, Y.; Behrem, S.; Ata, N.; Karaca, O.; Elmaci, C.; Cemal, I. Genome-wide association studies of preweaning growth and in vivo carcass composition traits in Esme sheep. J. Anim. Breed. Genet. 2022, 139, 26–39. [Google Scholar] [CrossRef]
  21. Lakhssassi, K.; Lahoz, B.; Sarto, P.; Iguácel, L.P.; Folch, J.; Alabart, J.L.; Serrano, M.; Calvo, J.H. Genome-Wide Association Study Demonstrates the Role Played by the CD226 Gene in Rasa Aragonesa Sheep Reproductive Seasonality. Animals 2021, 11, 1171. [Google Scholar] [CrossRef]
  22. Kijas, J.W.; Lenstra, J.A.; Hayes, B.; Boitard, S.; Porto Neto, L.R.; San Cristobal, M.; Servin, B.; McCulloch, R.; Whan, V.; Gietzen, K.; et al. Genome-wide analysis of the world’s sheep breeds reveals high levels of historic mixture and strong recent selection. PLoS Biol. 2012, 10, e1001258. [Google Scholar] [CrossRef]
  23. Kijas, J.W.; Townley, D.; Dalrymple, B.P.; Heaton, M.P.; Maddox, J.F.; McGrath, A.; Wilson, P.; Ingersoll, R.G.; McCulloch, R.; McWilliam, S.; et al. A genome wide survey of SNP variation reveals the genetic structure of sheep breeds. PLoS ONE 2009, 4, e4668. [Google Scholar] [CrossRef]
  24. Ding, R.; Yang, M.; Quan, J.; Li, S.; Zhuang, Z.; Zhou, S.; Zheng, E.; Hong, L.; Li, Z.; Cai, G.; et al. Single-Locus and Multi-Locus Genome-Wide Association Studies for Intramuscular Fat in Duroc Pigs. Front. Genet. 2019, 10, 619. [Google Scholar] [CrossRef]
  25. Zhao, Y.; Pu, Y.; Liang, B.; Bai, T.; Liu, Y.; Jiang, L.; Ma, Y. A study using single-locus and multi-locus genome-wide association study to identify genes associated with teat number in Hu sheep. Anim. Genet. 2022, 53, 203–211. [Google Scholar] [CrossRef] [PubMed]
  26. Uffelmann, E.; Huang, Q.Q.; Munung, N.S.; Vries, J.d.; Okada, Y.; Martin, A.R.; Martin, H.C.; Lappalainen, T.; Posthuma, D. Genome-wide association studies. Nat. Rev. Methods Prim. 2021, 1, 59. [Google Scholar] [CrossRef]
  27. Marees, A.T.; Kluiver, H.d.; Stringer, S.; Vorspan, F.; Curis, E.; Marie-Claire, C.; Derks, E.M. A tutorial on conducting genome-wide association studies: Quality control and statistical analysis. Int. J. Methods Psychiatr. Res. 2018, 27, e1608. [Google Scholar] [CrossRef]
  28. Yang, J.; Zaitlen, N.A.; Goddard, M.E.; Visscher, P.M.; Price, A.L. Advantages and pitfalls in the application of mixed-model association methods. Nat. Genet. 2014, 46, 100–106. [Google Scholar] [CrossRef] [PubMed]
  29. Kang, H.M.; Sul, J.H.; Service, S.K.; Zaitlen, N.A.; Kong, S.; Freimer, N.B.; Sabatti, C.; Eskin, E. Variance component model to account for sample structure in genome-wide association studies. Nat. Genet. 2010, 42, 348–354. [Google Scholar] [CrossRef]
  30. Zhou, X.; Stephens, M. Genome-wide efficient mixed-model analysis for association studies. Nat. Genet. 2012, 44, 821–824. [Google Scholar] [CrossRef]
  31. Li, C.; Fu, Y.; Sun, R.; Wang, Y.; Wang, Q. Single-Locus and Multi-Locus Genome-Wide Association Studies in the Genetic Dissection of Fiber Quality Traits in Upland Cotton (Gossypium hirsutum L.). Front. Plant Sci. 2018, 9, 1083. [Google Scholar] [CrossRef]
  32. Segura, V.; Vilhjálmsson, B.J.; Platt, A.; Korte, A.; Seren, Ü.; Long, Q.; Nordborg, M. An efficient multi-locus mixed-model approach for genome-wide association studies in structured populations. Nat. Genet. 2012, 44, 825–830. [Google Scholar] [CrossRef]
  33. Wang, S.-B.; Feng, J.-Y.; Ren, W.-L.; Huang, B.; Zhou, L.; Wen, Y.-J.; Zhang, J.; Dunwell, J.M.; Xu, S.; Zhang, Y.-M. Improving power and accuracy of genome-wide association studies via a multi-locus mixed linear model methodology. Sci. Rep. 2016, 6, 19444. [Google Scholar] [CrossRef]
  34. Tamba, C.L.; Zhang, Y.-M. A fast mrMLM algorithm for multi-locus genome-wide association studies. bioRxiv 2018, bioRxiv:341784. [Google Scholar] [CrossRef]
  35. Wen, Y.-J.; Zhang, H.; Ni, Y.-L.; Huang, B.; Zhang, J.; Feng, J.-Y.; Wang, S.-B.; Dunwell, J.M.; Zhang, Y.-M.; Wu, R. Methodological implementation of mixed linear models in multi-locus genome-wide association studies. Brief. Bioinform. 2017, 18, 906. [Google Scholar] [CrossRef] [PubMed]
  36. Zhang, J.; Feng, J.-Y.; Ni, Y.-L.; Wen, Y.-J.; Niu, Y.; Tamba, C.L.; Yue, C.; Song, Q.; Zhang, Y.-M. pLARmEB: Integration of least angle regression with empirical Bayes for multilocus genome-wide association studies. Heredity 2017, 118, 517–524. [Google Scholar] [CrossRef] [PubMed]
  37. Ren, W.-L.; Wen, Y.-J.; Dunwell, J.M.; Zhang, Y.-M. pKWmEB: Integration of Kruskal-Wallis test with empirical Bayes under polygenic background control for multi-locus genome-wide association study. Heredity 2018, 120, 208–218. [Google Scholar] [CrossRef] [PubMed]
  38. Tamba, C.L.; Ni, Y.-L.; Zhang, Y.-M. Iterative sure independence screening EM-Bayesian LASSO algorithm for multi-locus genome-wide association studies. PLoS Comput. Biol. 2017, 13, e1005357. [Google Scholar] [CrossRef] [PubMed]
  39. Öztan, S.; Düring, R.-A. Microwave assisted EDTA extraction-determination of pseudo total contents of distinct trace elements in solid environmental matrices. Talanta 2012, 99, 594–602. [Google Scholar] [CrossRef]
  40. Gogarten, S.M.; Sofer, T.; Chen, H.; Yu, C.; Brody, J.A.; Thornton, T.A.; Rice, K.M.; Conomos, M.P. Genetic association testing using the GENESIS R/Bioconductor package. Bioinformatics 2019, 35, 5346–5348. [Google Scholar] [CrossRef]
  41. Johnson, R.C.; Nelson, G.W.; Troyer, J.L.; Lautenberger, J.A.; Kessing, B.D.; Winkler, C.A.; O’Brien, S.J. Accounting for multiple comparisons in a genome-wide association study (GWAS). BMC Genom. 2010, 11, 724. [Google Scholar] [CrossRef]
  42. Nelson, C.P.; Goel, A.; Butterworth, A.S.; Kanoni, S.; Webb, T.R.; Marouli, E.; Zeng, L.; Ntalla, I.; Lai, F.Y.; Hopewell, J.C.; et al. Association analyses based on false discovery rate implicate new loci for coronary artery disease. Nat. Genet. 2017, 49, 1385–1391. [Google Scholar] [CrossRef]
  43. Benjamini, Y.; Hochberg, Y. Controlling the False Discovery Rate: A Practical and Powerful Approach to Multiple Testing. J. R. Stat. Soc. Ser. B 1995, 57, 289–300. [Google Scholar] [CrossRef]
  44. Wang, Y.; Ding, X.; Tan, Z.; Ning, C.; Xing, K.; Yang, T.; Pan, Y.; Sun, D.; Wang, C. Genome-Wide Association Study of Piglet Uniformity and Farrowing Interval. Front. Genet. 2017, 8, 194. [Google Scholar] [CrossRef]
  45. R Core Team. R: A Language and Environment for Statistical Computing; R Foundation for Statistical Computing: Vienna, Austria, 2023; Available online: https://www.R-project.org/ (accessed on 20 January 2023).
  46. Yang, J.; Zeng, J.; Goddard, M.E.; Wray, N.R.; Visscher, P.M. Concepts, estimation and interpretation of SNP-based heritability. Nat. Genet. 2017, 49, 1304–1310. [Google Scholar] [CrossRef] [PubMed]
  47. Zhang, Y.-W.; Tamba, C.L.; Wen, Y.-J.; Li, P.; Ren, W.-L.; Ni, Y.-L.; Gao, J.; Zhang, Y.-M. mrMLM v4.0.2: An R Platform for Multi-locus Genome-wide Association Studies. Genom. Proteom. Bioinform. 2020, 18, 481–487. [Google Scholar] [CrossRef] [PubMed]
  48. Karolchik, D.; Hinrichs, A.S.; Kent, W.J. The UCSC Genome Browser. Curr. Protoc. Hum. Genet. 2011, 71, 18.6.1–18.6.33. [Google Scholar] [CrossRef] [PubMed]
  49. Huang, W.D.; Sherman, B.T.; Lempicki, R.A. Bioinformatics enrichment tools: Paths toward the comprehensive functional analysis of large gene lists. Nucleic Acids Res. 2009, 37, 1–13. [Google Scholar] [CrossRef]
  50. Huang, W.D.; Sherman, B.T.; Lempicki, R.A. Systematic and integrative analysis of large gene lists using DAVID bioinformatics resources. Nat. Protoc. 2009, 4, 44–57. [Google Scholar] [CrossRef]
  51. Dennis, G.; Sherman, B.T.; Hosack, D.A.; Yang, J.; Gao, W.; Lane, H.C.; Lempicki, R.A. DAVID: Database for Annotation, Visualization, and Integrated Discovery. Genome Biol. 2003, 4, R60. [Google Scholar] [CrossRef]
  52. Rivals, I.; Personnaz, L.; Taing, L.; Potier, M.-C. Enrichment or depletion of a GO category within a class of genes: Which test? Bioinformatics 2007, 23, 401–407. [Google Scholar] [CrossRef]
  53. Judson, G.J.; Walkley, J.R.W.; James, P.J.; Kleemann, D.O.; Ponzoni, R.W. Genetic variation in trace element status of Merino sheep. Proc. Aust. Soc. Anim. Prod. 1994, 20, 438. [Google Scholar]
  54. Wang, J.; Fedoseienko, A.; Chen, B.; Burstein, E.; Jia, D.; Billadeau, D.D. Endosomal receptor trafficking: Retromer and beyond. Traffic 2018, 19, 578–590. [Google Scholar] [CrossRef]
  55. Liu, J.-J. Regulation of dynein-dynactin-driven vesicular transport. Traffic 2017, 18, 336–347. [Google Scholar] [CrossRef]
  56. Rebsamen, M.; Pochini, L.; Stasyk, T.; Araújo, M.E.G.d.; Galluccio, M.; Kandasamy, R.K.; Snijder, B.; Fauster, A.; Rudashevskaya, E.L.; Bruckner, M.; et al. SLC38A9 is a component of the lysosomal amino acid sensing machinery that controls mTORC1. Nature 2015, 519, 477–481. [Google Scholar] [CrossRef] [PubMed]
  57. Pfister, K.K. Distinct functional roles of cytoplasmic dynein defined by the intermediate chain isoforms. Exp. Cell Res. 2015, 334, 54–60. [Google Scholar] [CrossRef] [PubMed]
  58. McKenney, R.J.; Huynh, W.; Tanenbaum, M.E.; Bhabha, G.; Vale, R.D. Activation of cytoplasmic dynein motility by dynactin-cargo adapter complexes. Science 2014, 345, 337–341. [Google Scholar] [CrossRef] [PubMed]
  59. Morgan, J.L.; Song, Y.; Barbar, E. Structural dynamics and multiregion interactions in dynein-dynactin recognition. J. Biol. Chem. 2011, 286, 39349–39359. [Google Scholar] [CrossRef] [PubMed]
  60. Lim, C.M.; Cater, M.A.; Mercer, J.F.B.; La Fontaine, S. Copper-dependent interaction of dynactin subunit p62 with the N terminus of ATP7B but not ATP7A. J. Biol. Chem. 2006, 281, 14006–14014. [Google Scholar] [CrossRef]
  61. Polishchuk, E.V.; Concilli, M.; Iacobacci, S.; Chesi, G.; Pastore, N.; Piccolo, P.; Paladino, S.; Baldantoni, D.; van IJzendoorn, S.C.D.; Chan, J.; et al. Wilson disease protein ATP7B utilizes lysosomal exocytosis to maintain copper homeostasis. Dev. Cell 2014, 29, 686–700. [Google Scholar] [CrossRef]
  62. Chen, K.-E.; Healy, M.D.; Collins, B.M. Towards a molecular understanding of endosomal trafficking by Retromer and Retriever. Traffic 2019, 20, 465–478. [Google Scholar] [CrossRef]
  63. Lucas, M.; Gershlick, D.C.; Vidaurrazaga, A.; Rojas, A.L.; Bonifacino, J.S.; Hierro, A. Structural Mechanism for Cargo Recognition by the Retromer Complex. Cell 2016, 167, 1623–1635. [Google Scholar] [CrossRef]
  64. Arredondo, M.; Muñoz, P.; Mura, C.V.; Nùñez, M.T. DMT1, a physiologically relevant apical Cu1+ transporter of intestinal cells. Am. J. Physiol.-Cell Physiol. 2003, 284, C1525–C1530. [Google Scholar] [CrossRef]
  65. Helfer, E.; Harbour, M.E.; Henriot, V.; Lakisic, G.; Sousa-Blin, C.; Volceanov, L.; Seaman, M.N.J.; Gautreau, A. Endosomal recruitment of the WASH complex: Active sequences and mutations impairing interaction with the retromer. Biol. Cell 2013, 105, 191–207. [Google Scholar] [CrossRef]
  66. Harbour, M.E.; Breusegem, S.Y.; Seaman, M.N.J. Recruitment of the endosomal WASH complex is mediated by the extended ‘tail’ of Fam21 binding to the retromer protein Vps35. Biochem. J. 2012, 442, 209–220. [Google Scholar] [CrossRef] [PubMed]
  67. Jia, D.; Gomez, T.S.; Billadeau, D.D.; Rosen, M.K. Multiple repeat elements within the FAM21 tail link the WASH actin regulatory complex to the retromer. Mol. Biol. Cell 2012, 23, 2352–2361. [Google Scholar] [CrossRef] [PubMed]
  68. Phillips-Krawczak, C.A.; Singla, A.; Starokadomskyy, P.; Deng, Z.; Osborne, D.G.; Li, H.; Dick, C.J.; Gomez, T.S.; Koenecke, M.; Zhang, J.-S.; et al. COMMD1 is linked to the WASH complex and regulates endosomal trafficking of the copper transporter ATP7A. Mol. Biol. Cell 2015, 26, 91–103. [Google Scholar] [CrossRef]
  69. Materia, S.; Cater, M.A.; Klomp, L.W.J.; Mercer, J.F.B.; La Fontaine, S. Clusterin and COMMD1 independently regulate degradation of the mammalian copper ATPases ATP7A and ATP7B. J. Biol. Chem. 2012, 287, 2485–2499. [Google Scholar] [CrossRef]
  70. Cox, D.W.; Moore, S.D.P. Copper Transporting P-Type ATPases and Human Disease. J. Bioenergy 2002, 34, 333–338. [Google Scholar] [CrossRef]
  71. Akram, Z.; Ahmed, I.; Mack, H.; Kaur, R.; Silva, R.C.; Castilho, B.A.; Friant, S.; Sattlegger, E.; Munn, A.L. Yeast as a Model to Understand Actin-Mediated Cellular Functions in Mammals-Illustrated with Four Actin Cytoskeleton Proteins. Cells 2020, 9, 672. [Google Scholar] [CrossRef]
  72. Zhou, F.-L.; Li, S.C.; Zhu, Y.; Guo, W.-J.; Shao, L.-J.; Nelson, J.; Simpkins, S.; Yang, D.-H.; Liu, Q.; Yashiroda, Y.; et al. Integrating yeast chemical genomics and mammalian cell pathway analysis. Acta Pharmacol. Sin. 2019, 40, 1245–1255. [Google Scholar] [CrossRef] [PubMed]
  73. Sowada, N.; Stiller, B.; Kubisch, C. Increased copper toxicity in Saccharomyces cerevisiae lacking VPS35, a component of the retromer and monogenic Parkinson disease gene in humans. Biochem. Biophys. Res. Commun. 2016, 476, 528–533. [Google Scholar] [CrossRef]
  74. Lei, H.-T.; Mu, X.; Hattne, J.; Gonen, T. A conformational change in the N terminus of SLC38A9 signals mTORC1 activation. Structure 2021, 29, 426–432. [Google Scholar] [CrossRef]
  75. Wyant, G.A.; Abu-Remaileh, M.; Wolfson, R.L.; Chen, W.W.; Freinkman, E.; Danai, L.V.; Vander Heiden, M.G.; Sabatini, D.M. mTORC1 Activator SLC38A9 Is Required to Efflux Essential Amino Acids from Lysosomes and Use Protein as a Nutrient. Cell 2017, 171, 642–654. [Google Scholar] [CrossRef]
  76. Nnah, I.C.; Wang, B.; Saqcena, C.; Weber, G.F.; Bonder, E.M.; Bagley, D.; Cegli, R.d.; Napolitano, G.; Medina, D.L.; Ballabio, A.; et al. TFEB-driven endocytosis coordinates MTORC1 signaling and autophagy. Autophagy 2019, 15, 151–164. [Google Scholar] [CrossRef] [PubMed]
  77. Peña-Llopis, S.; Brugarolas, J. TFEB, a novel mTORC1 effector implicated in lysosome biogenesis, endocytosis and autophagy. Cell Cycle 2011, 10, 3987–3988. [Google Scholar] [CrossRef] [PubMed]
  78. Calvo, J.; Jung, H.; Meloni, G. Copper metallothioneins. IUBMB Life 2017, 69, 236–245. [Google Scholar] [CrossRef] [PubMed]
  79. Woodman, P. ESCRT-III on endosomes: New functions, new activation pathway. Biochem. J. 2016, 473, e5–e8. [Google Scholar] [CrossRef]
  80. Henne, W.M.; Buchkovich, N.J.; Emr, S.D. The ESCRT pathway. Dev. Cell 2011, 21, 77–91. [Google Scholar] [CrossRef]
  81. Tanikawa, S.; Mori, F.; Tanji, K.; Kakita, A.; Takahashi, H.; Wakabayashi, K. Endosomal sorting related protein CHMP2B is localized in Lewy bodies and glial cytoplasmic inclusions in α-synucleinopathy. Neurosci. Lett. 2012, 527, 16–21. [Google Scholar] [CrossRef]
  82. Bonora, M.; Giorgi, C.; Pinton, P. Molecular mechanisms and consequences of mitochondrial permeability transition. Nat. Rev. Mol. Cell Biol. 2022, 23, 266–285. [Google Scholar] [CrossRef]
  83. Shanmughapriya, S.; Rajan, S.; Hoffman, N.E.; Higgins, A.M.; Tomar, D.; Nemani, N.; Hines, K.J.; Smith, D.J.; Eguchi, A.; Vallem, S.; et al. SPG7 Is an Essential and Conserved Component of the Mitochondrial Permeability Transition Pore. Mol. Cell 2015, 60, 47–62. [Google Scholar] [CrossRef]
  84. Bonora, M.; Bononi, A.; Marchi, E.d.; Giorgi, C.; Lebiedzinska, M.; Marchi, S.; Patergnani, S.; Rimessi, A.; Suski, J.M.; Wojtala, A.; et al. Role of the c subunit of the FO ATP synthase in mitochondrial permeability transition. Cell Cycle 2013, 12, 674–683. [Google Scholar] [CrossRef]
  85. Rueda, C.B.; Llorente-Folch, I.; Gonzalez-Sanchez, P.; Contreras, L.; Juaristi, I.; Martinez-Valero, P.; Pardo, B.; Del Arco, A.; Satrustegui, J. Glutamate excitotoxicity and Ca2+-regulation of respiration: Role of the Ca2+ activated mitochondrial transporters (CaMCs). Biochim. Biophys. Acta (BBA)-Bioenerg. 2016, 1857, e17–e18. [Google Scholar] [CrossRef]
  86. Llorente-Folch, I.; Rueda, C.B.; Amigo, I.; Del Arco, A.; Saheki, T.; Pardo, B.; Satrústegui, J. Calcium-regulation of mitochondrial respiration maintains ATP homeostasis and requires ARALAR/AGC1-malate aspartate shuttle in intact cortical neurons. J. Neurosci. 2013, 33, 13957–13971. [Google Scholar] [CrossRef] [PubMed]
  87. Horn, D.; Barrientos, A. Mitochondrial copper metabolism and delivery to cytochrome c oxidase. IUBMB Life 2008, 60, 421–429. [Google Scholar] [CrossRef] [PubMed]
  88. Hosseini, M.-J.; Shaki, F.; Ghazi-Khansari, M.; Pourahmad, J. Toxicity of copper on isolated liver mitochondria: Impairment at complexes I, II, and IV leads to increased ROS production. Cell Biochem. Biophys. 2014, 70, 367–381. [Google Scholar] [CrossRef] [PubMed]
  89. Bernardi, P.; Krauskopf, A.; Basso, E.; Petronilli, V.; Blalchy-Dyson, E.; Di Lisa, F.; Forte, M.A. The mitochondrial permeability transition from in vitro artifact to disease target. FEBS J. 2006, 273, 2077–2099. [Google Scholar] [CrossRef]
  90. Anderegg, M.A.; Gyimesi, G.; Ho, T.M.; Hediger, M.A.; Fuster, D.G. The Less Well-Known Little Brothers: The SLC9B/NHA Sodium Proton Exchanger Subfamily-Structure, Function, Regulation and Potential Drug-Target Approaches. Front. Physiol. 2022, 13, 898508. [Google Scholar] [CrossRef] [PubMed]
  91. Kojima, A.; Toshima, J.Y.; Kanno, C.; Kawata, C.; Toshima, J. Localization and functional requirement of yeast Na+/H+ exchanger, Nhx1p, in the endocytic and protein recycling pathway. Biochim. Biophys. Acta 2012, 1823, 534–543. [Google Scholar] [CrossRef] [PubMed]
  92. Bowers, K.; Levi, B.P.; Patel, F.I.; Stevens, T.H. The sodium/proton exchanger Nhx1p is required for endosomal protein trafficking in the yeast Saccharomyces cerevisiae. Mol. Biol. Cell 2000, 11, 4277–4294. [Google Scholar] [CrossRef]
  93. Karim, M.A.; Brett, C.L. The Na+(K+)/H+ exchanger Nhx1 controls multivesicular body–vacuolar lysosome fusion. Mol. Biol. Cell 2018, 29, 317–325. [Google Scholar] [CrossRef]
Genes 14 01053 g001 550
Figure 1. Manhattan plot of single-locus genome-wide association study (SL-GWAS) of liver Cu concentration for Merinoland sheep. The dashed line represents the FDR threshold (−log10(p-value) = 4.06).
Figure 1. Manhattan plot of single-locus genome-wide association study (SL-GWAS) of liver Cu concentration for Merinoland sheep. The dashed line represents the FDR threshold (−log10(p-value) = 4.06).
Genes 14 01053 g001
Genes 14 01053 g002 550
Figure 2. Manhattan plots of the five multi-locus genome-wide association studies (ML-GWAS) (a) mrMLM, (b) FASTmrMLM, (c) FASTmrEMMA and (d) pLARmEB, (e) ISIS EM-BLASSO for liver Cu concentration in Merinoland sheep. Yellow lines with red dots represent significant SNPs above the LOD thresholds. Dashed lines represent the thresholds for significance (LOD > 3).
Figure 2. Manhattan plots of the five multi-locus genome-wide association studies (ML-GWAS) (a) mrMLM, (b) FASTmrMLM, (c) FASTmrEMMA and (d) pLARmEB, (e) ISIS EM-BLASSO for liver Cu concentration in Merinoland sheep. Yellow lines with red dots represent significant SNPs above the LOD thresholds. Dashed lines represent the thresholds for significance (LOD > 3).
Genes 14 01053 g002
Table
Table 1. Descriptive statistics and estimated heritability of liver copper concentration for Merinoland sheep lambs used in SL-GWAs and ML-GWAS analyses.
Table 1. Descriptive statistics and estimated heritability of liver copper concentration for Merinoland sheep lambs used in SL-GWAs and ML-GWAS analyses.
FarmN 1Mean (SE) Cu (mg/kg DM)MinMaxMean (SE) Age at Slaughter (Days)MinMaxh2(SE) 2
189159.54 ± 5.2267.67273.42150.76 ± 1.301191650.67(0.29)
24179.05 ± 6.8020.57188.28136.78 ± 3.0489172
1 N refers to number of samples; 2 h2 refers to heritability of liver copper concentration.
Table
Table 2. Significant SNPs and genes associated with liver Cu concentration for Merinoland sheep using SL-GWAS.
Table 2. Significant SNPs and genes associated with liver Cu concentration for Merinoland sheep using SL-GWAS.
CHRSNP 1Positionp-ValueGenes 2 (1 Mb)
2OAR2_145591151.11368971365.11 × 10−5SLC25A12, DYNC1I2
2s62875.11509051306.59 × 10−5ACVR1C, ACVR1
1 SNP positions are based on Oar_v4.0 assembly in NCBI. 2 Candidate genes: underlined = positional candidate genes; bold = functional candidate genes within of SNP within a 1 Mb region.
Table
Table 3. Significant SNPs and genes associated with liver Cu concentration for Merinoland sheep using ML-GWAS.
Table 3. Significant SNPs and genes associated with liver Cu concentration for Merinoland sheep using ML-GWAS.
CHRSNPPosition (BP) 1Genes (1 Mb) 2ML-GWAS Methods 3
1s66850.177025732GPR88, ATP5MF2, 3, 4, 5
1OAR1_285395930.1263664109COL6A1, COL6A22, 4, 5
2s05644.114172465FRRS1L1, 2, 4
2OAR2_145591151.1 *136897136SLC25A12, DYNC1I23, 4
3OAR3_132833292.1124516955KITLG1, 2, 4
4OAR4_77358490.173042599ZNF804B2, 3, 5
6s42668.121789865CENPE, SLC9B2, SLC9B22, 4
6OAR6_47263223.142315913GBA32, 3, 4, 5
7s25674.187187746NRXN32, 3, 4
7OAR7_101357352.193147667-1, 3, 5
14OAR14_14650208.114404415SPG7, CHMP1A, SHCBP1, VPS353, 5
16OAR16_25377664_X.123278766IL6ST, SLC38A92, 4
23OAR23_37101686.135080540GREB1L2, 4
1 SNP positions are based on Oar_v4.0 assembly in NCBI. 2 Candidate genes: underlined = positional candidate genes; bold = functional candidate genes of SNP within a 1 Mb region. 3 ML-GWAS Methods: 1 = mrMLM, 2 = FASTmrMLM, 3 = FASTmrEMMA, 4 = pLARmEB, 5 = ISIS EM-BLASSO; * SNP identified by both SL-GWAS and ML-GWAS.
Table
Table 4. Enriched GO terms determined by DAVID from genes identified in regions (1 Mb) of SNPs for liver copper variation in Merinoland sheep using SL-GWAS and ML-GWAS.
Table 4. Enriched GO terms determined by DAVID from genes identified in regions (1 Mb) of SNPs for liver copper variation in Merinoland sheep using SL-GWAS and ML-GWAS.
CategoryTermGenesp-Value 1
GOTERM_MF_DIRECTGO:0016361~activin receptor activity, type IACVR1, ACVR1C0.006
GOTERM_CC_DIRECTGO:0005765~lysosomal membraneCOL6A1, CHMP1A, VPS35, SLC38A90.007
GOTERM_CC_DIRECTGO:0048179~activin receptor complexACVR1, ACVR1C0.008
GOTERM_BP_DIRECTGO:1902600~hydrogen ion transmembrane transportSLC9B1, ATP5MF, SLC9B20.010
GOTERM_CC_DIRECTGO:0005743~mitochondrial inner membraneSPG7, SLC25A12, ATP5MF, SLC9B20.014
GOTERM_CC_DIRECTGO:0016021~integral component of membraneACVR1, KITLG, FRRS1L, GREB1L, NRXN3, SLC38A9, SPG7, IL6ST, SLC25A12, SLC9B1, ATP5MF, SLC9B20.014
GOTERM_MF_DIRECTGO:0015385~sodium:proton antiporter activitySLC9B1, SLC9B20.015
GOTERM_CC_DIRECTGO:0030496~midbodyCENPE, CHMP1A, SHCBP10.016
GOTERM_CC_DIRECTGO:0005828~kinetochore microtubuleCENPE, CHMP1A0.021
GOTERM_BP_DIRECTGO:0032924~activin receptor signalling pathwayACVR1, ACVR1C0.023
GOTERM_CC_DIRECTGO:0043235~receptor complexACVR1, ACVR1C, IL6ST0.023
GOTERM_CC_DIRECTGO:0097228~sperm principal pieceSLC9B1, SLC9B20.035
GOTERM_MF_DIRECTGO:0019838~growth factor bindingACVR1C, IL6ST0.040
GOTERM_MF_DIRECTGO:0030020~extracellular matrix structural constituent conferring tensile strengthCOL6A2, COL6A10.045
GOTERM_BP_DIRECTGO:0007080~mitotic metaphase plate congressionCENPE, CHMP1A0.047
GOTERM_BP_DIRECTGO:0001755~neural crest cell migrationACVR1, KITLG0.049
1 p = p-value significant at 0.05; Please note that GOTERM_MF, GOTERM_CC and GOTERM_BP refer to molecular function, cellular component and biological process, respectively.
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Adeniyi, O.O.; Medugorac, I.; Grochowska, E.; Düring, R.-A.; Lühken, G. Single-Locus and Multi-Locus Genome-Wide Association Studies Identify Genes Associated with Liver Cu Concentration in Merinoland Sheep. Genes 2023, 14, 1053. https://doi.org/10.3390/genes14051053

AMA Style

Adeniyi OO, Medugorac I, Grochowska E, Düring R-A, Lühken G. Single-Locus and Multi-Locus Genome-Wide Association Studies Identify Genes Associated with Liver Cu Concentration in Merinoland Sheep. Genes. 2023; 14(5):1053. https://doi.org/10.3390/genes14051053

Chicago/Turabian Style

Adeniyi, Olusegun O., Ivica Medugorac, Ewa Grochowska, Rolf-Alexander Düring, and Gesine Lühken. 2023. "Single-Locus and Multi-Locus Genome-Wide Association Studies Identify Genes Associated with Liver Cu Concentration in Merinoland Sheep" Genes 14, no. 5: 1053. https://doi.org/10.3390/genes14051053

APA Style

Adeniyi, O. O., Medugorac, I., Grochowska, E., Düring, R. -A., & Lühken, G. (2023). Single-Locus and Multi-Locus Genome-Wide Association Studies Identify Genes Associated with Liver Cu Concentration in Merinoland Sheep. Genes, 14(5), 1053. https://doi.org/10.3390/genes14051053

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop