Nutrient supply affects the mRNA expression profile of the porcine skeletal muscle

Integrated transcriptomics and proteomics analysis reveals muscle metabolism effects of dietary Ulva lactuca and ulvan lyase supplementation in weaned piglets

Seaweeds, including the green Ulva lactuca, can potentially reduce competition between feed, food, and fuel. They can also contribute to the improved development of weaned piglets. However, their indigestible polysaccharides of the cell wall pose a challenge. This can be addressed through carbohydrase supplementation, such as the recombinant ulvan lyase. The objective of our study was to assess the muscle metabolism of weaned piglets fed with 7% U. lactuca and 0.01% ulvan lyase supplementation, using an integrated transcriptomics (RNA-seq) and proteomics (LC-MS) approach. Feeding piglets with seaweed and enzyme supplementation resulted in reduced macronutrient availability, leading to protein degradation through the proteasome (PSMD2), with resulting amino acids being utilized as an energy source (GOT2, IDH3B). Moreover, mineral element accumulation may have contributed to increased oxidative stress, evident from elevated levels of antioxidant proteins like catalase, as a response to maintaining tissue homeostasis. The upregulation of the gene AQP7, associated with the osmotic stress response, further supports these findings. Consequently, an increase in chaperone activity, including HSP90, was required to repair damaged proteins. Our results suggest that enzymatic supplementation may exacerbate the effects observed from feeding U. lactuca alone, potentially due to side effects of cell wall degradation during digestion.

Transcriptome-based analysis of early post-mortem formation of pale, soft, and exudative (PSE) pork

Pale, soft, and exudative (PSE) meat can cause consumer dissatisfaction and economic losses. This study determined meat quality, glycolytic enzyme activity, and differential gene expression in the longissimus lumborum (LL) and semimembranosus (SM) of normal and PSE pork carcasses. The SM did not result in PSE meat. Hexokinase, lactate dehydrogenase, and pyruvate kinase activities were lower in the SM of PSE carcasses than in the normal carcasses. Functional enrichment analysis revealed that immune, inflammatory, and muscle fibre genes were significantly enriched in PSE pork. More specifically, PPP1R3G and MSS51 may be key genes regulating pork quality in the SM. Meanwhile, the differential expression of PLVAB, ADIPOQ, LEP, MYH4, MYH7, MYL3, MYL6B, FOS, ATF3, and HSPA6 may induce PSE formation in the LL. These results may provide insights into PSE pork formation mechanisms and reveal candidate genes for improving meat quality after validation.

Modeling microRNA-driven post-transcriptional regulation using exon-intron split analysis in pigs

The contribution of microRNAs (miRNAs) to mRNA post-transcriptional regulation has often been explored by the post hoc selection of downregulated genes and determining whether they harbor binding sites for miRNAs of interest. This approach, however, does not discriminate whether these mRNAs are also downregulated at the transcriptional level. Here, we have characterized the transcriptional and post-transcriptional changes in mRNA expression in two porcine tissues: gluteus medius muscle of fasted and fed Duroc gilts and adipose tissue of lean and obese Duroc-Göttingen minipigs. Exon-intron split analysis of RNA-seq data allowed us to identify downregulated mRNAs with high post-transcriptional signals in fed or obese states, and we assessed whether they harbor binding sites for upregulated miRNAs in any of these two physiological states. We found 26 downregulated mRNAs with high post-transcriptional signals in the muscle of fed gilts and 21 of these were predicted targets of miRNAs upregulated in fed pigs. For adipose tissue, 44 downregulated mRNAs in obese minipigs displayed high post-transcriptional signals, and 25 of these were predicted targets of miRNAs upregulated in the obese state. These results suggest that the contribution of miRNAs to mRNA repression is more prominent in the skeletal muscle system. Finally, we identified several genes that may play relevant roles in the energy homeostasis of the pig skeletal muscle (DKK2 and PDK4) and adipose (SESN3 and ESRRG) tissues. By differentiating transcriptional from post-transcriptional changes in mRNA expression, exon-intron split analysis provides a valuable view of the regulation of gene expression, complementary to canonical differential expression analyses.

Transcriptome shifts triggered by vitamin A and SCD genotype interaction in Duroc pigs

BACKGROUND

The composition of intramuscular fat depends on genetic and environmental factors, including the diet. In pigs, we identified a haplotype of three SNP mutations in the stearoyl-coA desaturase (SCD) gene promoter associated with higher content of monounsaturated fatty acids in intramuscular fat. The second of these three SNPs (rs80912566, C > T) affected a putative retinol response element in the SCD promoter. The effect of dietary vitamin A restriction over intramuscular fat content is controversial as it depends on the pig genetic line and the duration of the restriction. This study aims to investigate changes in the muscle transcriptome in SCD rs80912566 TT and CC pigs fed with and without a vitamin A supplement during the fattening period.

RESULTS

Vitamin A did not affect carcass traits or intramuscular fat content and fatty acid composition, but we observed an interaction between vitamin A and SCD genotype on the desaturation of fatty acids in muscle. As reported before, the SCD-TT pigs had more monounsaturated fat than the SCD-CC animals. The diet lacking the vitamin A supplement enlarged fatty acid compositional differences between SCD genotypes, partly because vitamin A had a bigger effect on fatty acid desaturation in SCD-CC pigs (positive) than in SCD-TT and SCD-TC animals (negative). The interaction between diet and genotype was also evident at the transcriptome level; the highest number of differentially expressed genes were detected between SCD-TT pigs fed with the two diets. The genes modulated by the diet with the vitamin A supplement belonged to metabolic and signalling pathways related to immunity and inflammation, transport through membrane-bounded vesicles, fat metabolism and transport, reflecting the impact of retinol on a wide range of metabolic processes.

CONCLUSIONS

Restricting dietary vitamin A during the fattening period did not improve intramuscular fat content despite relevant changes in muscle gene expression, both in coding and non-coding genes. Vitamin A activated general pathways of retinol response in a SCD genotype-dependant manner, which affected the monounsaturated fatty acid content, particularly in SCD-CC pigs.

The Circadian Physiology: Implications in Livestock Health

Circadian rhythms exist in almost all types of cells in mammals. Thousands of genes exhibit approximately 24 h oscillations in their expression levels, making the circadian clock a crucial regulator of their normal functioning. In this regard, environmental factors to which internal physiological processes are synchronized (e.g., nutrition, feeding/eating patterns, timing and light exposure), become critical to optimize animal physiology, both by managing energy use and by realigning the incompatible processes. Once the circadian clock is disrupted, animals will face the increased risks of diseases, especially metabolic phenotypes. However, little is known about the molecular components of these clocks in domestic species and by which they respond to external stimuli. Here we review evidence for rhythmic control of livestock production and summarize the associated physiological functions, and the molecular mechanisms of the circadian regulation in pig, sheep and cattle. Identification of environmental and physiological inputs that affect circadian gene expressions will help development of novel targets and the corresponding approaches to optimize production efficiency in farm animals.

Molecular Characterization of Paralichthys olivaceus MAF1 and Its Potential Role as an Anti-Viral Hemorrhagic Septicaemia Virus Factor in Hirame Natural Embryo Cells

MAF1 is a global suppressor of RNA polymerase III-dependent transcription, and is conserved from yeast to human. Growing evidence supports the involvement of MAF1 in the immune response of mammals, but its biological functions in fish are unknown. We isolated and characterized Maf1 from the olive flounder Paralichthys olivaceus (PoMaf1). The coding region of PoMaf1 comprised 738 bp encoding a 245-amino-acid protein. The deduced PoMAF1 amino acid sequence shared features with those of MAF1 orthologues from vertebrates. PoMaf1 mRNA was detected in all tissues examined, and the levels were highest in eye and muscle tissue. The PoMaf1 mRNA level increased during early development. In addition, the PoMaf1 transcript level decreased during viral hemorrhagic septicemia virus (VHSV) infection of flounder hirame natural embryo (HINAE) cells. To investigate the role of PoMaf1 in VHSV infection, single-cell-derived PoMaf1 knockout HINAE cells were generated using the clustered regularly interspaced short palindromic repeats/CRISPR-associated-9 (CRISPR/Cas9) system, and cell clones with complete disruption of PoMaf1 were selected. PoMaf1 disruption increased the VHSV glycoprotein (G) mRNA levels during VHSV infection of HINAE cells, implicating PoMAF1 in the immune response to VSHV infection. To our knowledge, this is the first study to characterize fish Maf1, which may play a role in the response to viral infection.

Genetic variations (eQTLs) in muscle transcriptome and mitochondrial genes, and trans-eQTL molecular pathways in feed efficiency from Danish breeding pigs

Feed efficiency (FE) is a key trait in pig production, as improvement in FE has positive economic and environmental impact. FE is a complex phenotype and testing animals for FE is costly. Therefore, there has been a desire to find functionally relevant single nucleotide polymorphisms (SNPs) as biomarkers, to improve our biological understanding of FE as well as accuracy of genomic prediction for FE. We have performed a cis- and trans- eQTL (expression quantitative trait loci) analysis, in a population of Danbred Durocs (N = 11) and Danbred Landrace (N = 27) using both a linear and ANOVA model based on muscle tissue RNA-seq. We analyzed a total of 1425x19179 or 2.7x10⁷ Gene-SNP combinations in eQTL detection models for FE. The 1425 genes were from RNA-Seq based differential gene expression analyses using 25880 genes related to FE and additionally combined with mitochondrial genes. The 19179 SNPs were from applying stringent quality control and linkage disequilibrium filtering on genotype data using a GGP Porcine HD 70k SNP array. We applied 1000 fold bootstrapping and enrichment analysis to further validate and analyze our detected eQTLs. We identified 13 eQTLs with FDR < 0.1, affecting several genes found in previous studies of commercial pig breeds. Examples include MYO19, CPT1B, ACSL1, IER5L, CPT1A, SUCLA2, CSRNP1, PARK7 and MFF. The bootstrapping results showed statistically significant enrichment (p-value<2.2x10^-16) of eQTLs with p-value < 0.01 in both cis and trans-eQTLs. Enrichment analysis of top trans-eQTLs revealed high enrichment for gene categories and gene ontologies associated with genomic context and expression regulation. This included transcription factors (p-value = 1.0x10^-13), DNA-binding (GO:0003677, p-value = 8.9x10^-14), DNA-binding transcription factor activity (GO:0003700,) nucleus gene (GO:0005634, p-value<2.2x10^-16), negative regulation of expression (GO:0010629, p-value<2.2x10^-16). These results would be useful for future genome assisted breeding of pigs to improve FE, and in the improved understanding of the functional mechanism of trans eQTLs.

Analysis of Transcriptome and miRNAome in the Muscle of Bamei Pigs at Different Developmental Stages

Simple Summary

The pigs is the most popular agricultural animal in the world. Muscle growth—which has the highest economic value in pigs—can be regulated by multiple genes and involves complex regulatory mechanisms. It is necessary to understand the dynamics of muscle transcriptome during development to understand the muscle development mechanism. However, the genes and miRNAs that play regulatory roles underlying differences in the meat quality of pigs remain unclear. In the current study, qRT-PCR, miRNA-Seq, and RNA-Seq were applied to analyze and verify muscle tissues of pigs from three different developmental stages and screened genes, miRNAs and pathways related to pig muscle development. This study focused on analyzing the mechanisms of muscle development and uncover the development differences in muscle from embryo to adult.

Abstract

The growth of skeletal muscle involves complex developmental processes that play an important part in the determinization of pork quality. The investigation of skeletal muscle mRNA or miRNA profiles is especially important for finding molecular approaches to improve meat quality in pig breeding. Therefore, we studied the transcriptome (mRNA and miRNA) profiles of skeletal muscle with RNA-Seq in three developmental stages of pigs: 65-day embryonic (E65), postnatal 0 days (natal) and 10 months (adult). We found 10,035, 9050 and 4841 differentially expressed (DE) genes for natal vs. E65, adult vs. E65 and adult vs. natal, 55, 101 and 85 DE miRNA for natal vs. E65, adult vs. E65 and adult vs. natal, respectively. In addition, the target genes of DE miRNA that was in a negative correlation with the corresponding miRNA in the same comparison group were selected for enrichment analysis. Gene Ontology terms were mainly classified into developmental processes. Pathway analysis revealed enrichment in the Rap1 signaling pathway, citrate cycle and oxidative phosphorylation and carbon. Finally, RT-PCR was employed for validating the level of expression of 11 DE miRNA and 14 DEGs. The transcriptome profiles of skeletal muscle from the different developmental stages of the Bamei pigs were obtained. From these data, hundreds of DE miRNA and mRNA, and the miRNA–mRNA regulatory network can provide valuable insights into further understanding of key molecular mechanisms and improving the meat quality in pig breeding.

Transcriptome Analysis of Skeletal Muscle in Pigs with Divergent Residual Feed Intake Phenotypes

Residual feed intake (RFI) is defined as the difference between the observed and expected feed intake for maintenance and growth requirements. In this study, the expression profiles of mRNAs and long noncoding RNAs (lncRNAs) from skeletal muscle in Duroc pigs with divergent RFI phenotypes were investigated by Illumina sequencing. Finally, a total of 2195 annotated lncRNAs and 1976 novel lncRNAs were obtained. About 210 mRNAs and 43 lncRNAs were differentially expressed among high and low RFI pigs. The differentially expressed mRNAs were potentially involved in the biological processes of lipid metabolism, extracellular matrix organization, cell proliferation, and cell adhesion. The lipolysis in skeletal muscle was increased in high RFI pigs, suggesting that high RFI pigs might need more energy than low RFI pigs. However, skeletal muscle development was increased in low RFI pigs. These results suggested that low RFI pigs might be more efficient in energy utilization during skeletal muscle growth. The function of lncRNA was also analyzed by target prediction. Nine lncRNAs might be candidate lncRNAs for the determination of RFI phenotype, by the regulation of the biological processes of lipid metabolism, cell proliferation, and cell adhesion. This study should facilitate a further understanding of the molecular mechanism for the determination of RFI phenotype in pigs.

Co-expression network analysis predicts a key role of microRNAs in the adaptation of the porcine skeletal muscle to nutrient supply

BACKGROUND

The role of non-coding RNAs in the porcine muscle metabolism is poorly understood, with few studies investigating their expression patterns in response to nutrient supply. Therefore, we aimed to investigate the changes in microRNAs (miRNAs), long intergenic non-coding RNAs (lincRNAs) and mRNAs muscle expression before and after food intake.

RESULTS

We measured the miRNA, lincRNA and mRNA expression levels in the gluteus medius muscle of 12 gilts in a fasting condition (AL-T0) and 24 gilts fed ad libitum during either 5 h. (AL-T1, N = 12) or 7 h. (AL-T2, N = 12) prior to slaughter. The small RNA fraction was extracted from muscle samples retrieved from the 36 gilts and sequenced, whereas lincRNA and mRNA expression data were already available. In terms of mean and variance, the expression profiles of miRNAs and lincRNAs in the porcine muscle were quite different than those of mRNAs. Food intake induced the differential expression of 149 (AL-T0/AL-T1) and 435 (AL-T0/AL-T2) mRNAs, 6 (AL-T0/AL-T1) and 28 (AL-T0/AL-T2) miRNAs and none lincRNAs, while the number of differentially dispersed genes was much lower. Among the set of differentially expressed miRNAs, we identified ssc-miR-148a-3p, ssc-miR-22-3p and ssc-miR-1, which play key roles in the regulation of glucose and lipid metabolism. Besides, co-expression network analyses revealed several miRNAs that putatively interact with mRNAs playing key metabolic roles and that also showed differential expression before and after feeding. One case example was represented by seven miRNAs (ssc-miR-148a-3p, ssc-miR-151-3p, ssc-miR-30a-3p, ssc-miR-30e-3p, ssc-miR-421-5p, ssc-miR-493-5p and ssc-miR-503) which putatively interact with the PDK4 mRNA, one of the master regulators of glucose utilization and fatty acid oxidation.

CONCLUSIONS

As a whole, our results evidence that microRNAs are likely to play an important role in the porcine skeletal muscle metabolic adaptation to nutrient availability.

Discovery and annotation of novel microRNAs in the porcine genome by using a semi-supervised transductive learning approach

Despite the broad variety of available microRNA (miRNA) prediction tools, their application to the discovery and annotation of novel miRNA genes in domestic species is still limited. In this study we designed a comprehensive pipeline (eMIRNA) for miRNA identification in the yet poorly annotated porcine genome and demonstrated the usefulness of implementing a motif search positional refinement strategy for the accurate determination of precursor miRNA boundaries. The small RNA fraction from gluteus medius skeletal muscle of 48 Duroc gilts was sequenced and used for the prediction of novel miRNA loci. Additionally, we selected the human miRNA annotation for a homology-based search of porcine miRNAs with orthologous genes in the human genome. A total of 20 novel expressed miRNAs were identified in the porcine muscle transcriptome and 27 additional novel porcine miRNAs were also detected by homology-based search using the human miRNA annotation. The existence of three selected novel miRNAs (ssc-miR-483, ssc-miR484 and ssc-miR-200a) was further confirmed by reverse transcription quantitative real-time PCR analyses in the muscle and liver tissues of Göttingen minipigs. In summary, the eMIRNA pipeline presented in the current work allowed us to expand the catalogue of porcine miRNAs and showed better performance than other commonly used miRNA prediction approaches. More importantly, the flexibility of our pipeline makes possible its application in other yet poorly annotated non-model species.

Effects of long-term feeding of rapeseed meal on skeletal muscle transcriptome, production efficiency and meat quality traits in Norwegian Landrace growing-finishing pigs

This study was performed to investigate the effects of dietary inclusion of 20% rapeseed meal (RSM) as an alternative to soybean meal (SBM) in a three-month feeding experiment with growing finishing pigs. Dietary alteration affected growth performance, several carcass traits and transcriptional responses in the skeletal muscle, but did not affect measured meat quality traits. In general, pigs fed the RSM test diet exhibited reduced growth performance compared to pigs on SBM control diet. Significant transcriptional changes in the skeletal muscle of growing pigs fed RSM diet were likely the consequence of an increased amount of fiber and higher polyunsaturated fatty acids, and presence of bioactive phytochemicals, such as glucosinolates. RNAseq pipeline using Tophat2-Cuffdiff identified 57 upregulated and 63 downregulated genes in RSM compared to SBM pigs. Significantly enriched among downregulated pathways was p53-mediated signalling involved in cellular proliferation, while activation of negative growth regulators (IER5, KLF10, BTG2, KLF11, RETREG1, PRUNE2) in RSM fed pigs provided further evidence for reduced proliferation and increased cellular death, in accordance with the observed reduction in performance traits. Upregulation of well-known metabolic controllers (PDK4, UCP3, ESRRG and ESRRB), involved in energy homeostasis (glucose and lipid metabolism, and mitochondrial function), suggested less available energy and nutrients in RSM pigs. Furthermore, several genes supported more pronounced proteolysis (ABTB1, OTUD1, PADI2, SPP1) and reduced protein synthesis (THBS1, HSF4, AP1S2) in RSM muscle tissue. In parallel, higher levels of NR4A3, PDK4 and FGF21, and a drop in adropin, ELOVL6 and CIDEC/FSP27 indicated increased lipolysis and fatty acid oxidation, reflective of lower dressing percentage. Finally, pigs exposed to RSM showed greater expression level of genes responsive to oxidative stress, indicated by upregulation of GPX1, GPX2, and TXNIP.

Polymorphisms of the cryptochrome 2 and mitoguardin 2 genes are associated with the variation of lipid-related traits in Duroc pigs

The genetic factors determining the phenotypic variation of porcine fatness phenotypes are still largely unknown. We investigated whether the polymorphism of eight genes (MIGA2, CRY2, NPAS2, CIART, ARNTL2, PER1, PER2 and PCK1), which display differential expression in the skeletal muscle of fasted and fed sows, is associated with the variation of lipid and mRNA expression phenotypes in Duroc pigs. The performance of an association analysis with the GEMMA software demonstrated that the rs330779504 SNP in the MIGA2 gene is associated with LDL concentration at 190 days (LDL₂, corrected P-value = 0.057). Moreover, the rs320439526 SNP of the CRY2 gene displayed a significant association with stearic acid content in the longissimus dorsi muscle (LD C18:0, corrected P-value = 0.015). Both SNPs were also associated with the mRNA levels of the corresponding genes in the gluteus medius skeletal muscle. From a biological perspective these results are meaningful because MIGA2 protein plays an essential role in mitochondrial fusion, a process tightly connected with the energy status of the cell, while CRY2 is a fundamental component of the circadian clock. However, inclusion of these two SNPs in chromosome-wide association analyses demonstrated that they are not located at the peaks of significance for the two traits under study (LDL₂ for rs330779504 and LD C18:0 for rs320439526), thus implying that these two SNPs do not have causal effects.

SUMOylation promotes survival and integration of neural stem cell grafts in ischemic stroke

Background

Neural stem cell (NSC)-based therapies hold great promise for treating diseases of the central nervous system (CNS). However, several fundamental problems still need to be overcome to fully exploit the clinical potential of NSC therapeutics. Chief among them is the limited survival of NSC grafts within hostile microenvironments.

Methods

Herein, we sought to engineer NSCs in an effort to increase graft survival within ischemic brain lesions via upregulation of global SUMOylation, a post-translational modification critically involved in mediating tolerance to ischemia/reperfusion.

Findings

NSCs overexpressing the SUMO E2-conjugase Ubc9 displayed resistance to oxygen-glucose-deprivation/restoration of oxygen/glucose (OGD/ROG) and enhanced neuronal differentiation in vitro, as well as increased survival and neuronal differentiation when transplanted in mice with transient middle cerebral artery occlusion in vivo.

Interpretation

Our work highlights a critical role for SUMOylation in NSC biology and identifies a biological pathway that can be targeted to increase the effectiveness of exogenous stem cell medicines in ischemic stroke.

Fund

Intramural Research Program of the NINDS/NIH, the Italian Multiple Sclerosis Foundation (FISM), the Bascule Charitable Trust, NIH-IRTA-OxCam and Wellcome Trust Research Training Fellowships.

•

Ubc9-overexpressing NSCs demonstrate enhanced neuronal differentiation.

•

Upregulating SUMOylation in NSCs increases resistance to ischemia/reperfusion in vitro.

•

Ubc9-overexpressing NSC grafts robustly integrate within the brain of mice post-stroke.

Short-term high-fat meal intake alters the expression of circadian clock-, inflammation-, and oxidative stress-related genes in human skeletal muscle

Dietary food, depending on timing, amount and composition can influence gene expression in various tissues. Here, we investigated the effect of high-fat meal diets of different compositions on the gene expression pattern of human skeletal muscle. Gene expression data of skeletal muscle samples from human volunteers prior and 4 h after the consumption of high lipid-containing meal consisting of either saturated-, monounsaturated- or polyunsaturated fatty acids were downloaded from the public repository. List of 843 differently expressed genes (DEGs) was generated. Functional analysis revealed that circadian rhythm-, inflammation- and oxidative stress-related genes are highly overrepresented among the DEGs. The magnitude of gene expression changes significantly increases with the saturation level of the dietary fatty acids and the majority of the DEGs are upregulated. We propose that, by altering circadian clock gene expression and inducing inflammation and oxidative stress, high lipid intake can contribute to muscle function decay in the long run.

Analysing the Expression of Eight Clock Genes in Five Tissues From Fasting and Fed Sows

In a previous study, we observed that circadian clock genes are differentially expressed in the skeletal muscle of fasting and fed sows. The goal of the current work was to investigate if these genes are also differentially expressed in tissues containing the central (hypothalamus) and peripheral (duodenum, dorsal fat, muscle, and liver) clocks. As animal material, we used 12 sows that fasted 12 h before slaughtering (T0) and 12 sows that were fed ad libitum 7 h prior slaughtering (T2). Tissue samples were collected immediately after slaughter and total RNA was subsequently extracted. The expression of the ARNTL, BHLHE40, CRY2, NPAS2, NR1D1, PER1, PER2, and SIK1 genes was measured by quantitative reverse transcription PCR. The numbers of clock genes showing differential expression before and after feeding varied depending on the tissue i.e., four in dorsal fat and duodenum, six in skeletal muscle, and seven in the liver. In contrast, none of the eight analysed genes displayed a significant differential expression in hypothalamus, the tissue where the central clock resides. This result supports that the differential expression of clock genes in the four tissues mentioned above is probably induced by nutrition and not by the central clock entrained by light. Moreover, we have observed that the NPAS2 and ARNTL genes display positive log2(FC) values in the five tissues under analysis, whilst the CRY2, PER1 (except dorsal fat) and PER2 (except hypothalamus) genes generally show negative log2(FC) values. Such result might be explained by the existence of a negative feedback loop between the ARNTL/NPAS2 and CRY/PER genes. Collectively, these results support that nutrition plays an important role in modulating the timing of porcine peripheral circadian clocks. Such regulation could be essential for coordinating the subsequent metabolic response to nutrient supply.