In vitro Silencing of Acetyl-CoA Carboxylase beta (ACACB) Gene Reduces Cholesterol Synthesis in Knockdown Chicken Myoblast Cells

The poultry industry provides cost-effective, healthy, and protein-enriched food for the growing population and achieving the nutritional security to the country. Excessive abdominal and subcutaneous fat deposition is one of the major setbacks to the poultry industry that reduces carcass yield and feed eciency. In chicken abdominal fat constitutes 20% of total body fat which make up 2–3% of live weight of the bird. In fatty acid metabolism, acetyl-CoA Carboxylase (ACC) is one of the key enzymes with two isoforms i.e. ACACA and ACACB each of which plays a different role. In chicken, ACACB is involved in the β-oxidation of fatty acids and thereby potentially regulating the quality of meat and egg. The RNAi strategy is widely used for silencing the target gene expression. In this study, we designed ve shRNA constructs and identied the most ecient shRNA molecule for silencing the ACACB gene under in vitro chicken embryo myoblast (CEM) primary cell culture system. After knocking down the ACACB gene, for understanding how fatty acid metabolism is regulated, we tracked the expression of key fatty acid metabolism genes like ACACA, FASN, SCD, ELOVL2, and CPT1. Also, checked the expression of immune response genes like IFNA, IFNB, and BLB1 in control as well as ACACB knockdown myoblast cells and observed no signicant difference. We observed the down-regulation of key fatty acid metabolism genes along with ACACB, which may leads to the less fat accumulation in CEM cells. We also estimated the cholesterol and triglycerides in control and ACACB knockdown myoblast cells and found a signicant difference between control and the knockdown cells. In vitro knockdown of the ACACB gene in a cell culture system by a short hairpin RNA (shRNA) expressing construct would help to produce a knockdown chicken with reduced fat deposition. lysate


Introduction
The major objective of poultry production is to provide cheaper, safe, healthy, and protein-enriched food for the growing population thereby achieving the nutritional security to the country. Due to the adoption of improved nutritional strategies and technological innovations creating better rearing conditions, the poultry industry has grown greatly in recent years by making the chicken attaining a nishing body weight of about 2 kg in a very short period. On the contrary, apart from resulting in excessive body fat deposition, this rapid growth rate leads to the high incidence of metabolic diseases and skeletal disorders with increased mortality. In broiler chicken abdominal fat constitutes 20% of total body fat and excessive body fat has been recognized as a major problem in the poultry industry, which make up 2-3% of the total live weight of the broiler chicken (Cahaner et al. 1986;Butterwith 1989; Crespo and Esteve-Garcia 2002; Haro 2005; Leenstra 1986). Fat has limited value both to the poultry producer and consumer. Hence, excessive abdominal and subcutaneous fat deposition is one of the major setbacks to the poultry industry that reduces carcass yield, feed e ciency and ultimately causing consumers non-acceptance of meat (Lippens 2003;Jennen 2004). These negative viewpoints are of signi cant concern for the rancher and processor as they can bring considerable economic losses (Leenstra 1986;Pym 1987; Gri n 1996; Zubair and Leeson 1996;Buys et al. 1998;Buyse 1999;Havenstein et al. 2003;Nikolova et al. 2007). Moreover, in recent years, preference for leaner meat has been increasing as the consumers are becoming more sensitive about the positive correlation between intake of fatty substances and onset of cardiovascular diseases (Leeson and Summers 1980;Pym 1987;Cable and Waldroup 1990).
In avian species, major site of fat deposition, the abdominal fat pad is the prime indicator of total body fat content (Becker et al. 1979; Thomas et al. 1983). The fat deposition in body tissues is the net result of absorption, de novo synthesis, and βoxidation of fatty acids (Saadoun and Leclercq 1986). As most of the traditional dietary approaches like restricted feeding techniques to combat excessive fattening in commercial broilers have resulted in decreased genetic potential for weight gain and failed in cost perspective, it would be more appropriate to emphasize the molecular regulation of fatty acid metabolism (Dunnington et al. 1986). In fatty acid metabolism, Acetyl-CoA Carboxylase (ACC) is the key enzyme with two isoforms, ACACA and ACACB each of which plays a different role concerning fat metabolism. The ACACA and ACACB catalyzes the carboxylation of acetyl-CoA to malonyl-CoA and malonyl-CoA generated via ACACA isoform mainly provides carbon units for fatty acid synthesis in lipogenic tissues like liver and adipose tissue whereas malonyl-CoA produced by ACACB isoform regulates fat oxidation by inhibiting carnitine palmitoyl transferase-1 (CPT-1), which controls the entry of long-chain fatty acids into the mitochondrial site of oxidation in nonlipogenic tissues like the heart and skeletal muscles (Abu-Elheiga 2000; Abu-Elheiga 2001). Hence, inhibition of ACACB may lead to reduced fat deposition by increasing fatty acid oxidation. Therefore, we can reduce the expression of ACACB for improving poultry meat in terms of quality (leaner meat) as well as quantity (weight gain). Feeding experiments on mice revealed that, apart from consuming more food, ACACB knockout mice had lower body fat than their wild-type counterparts and were protected from high fat-induced obesity ( Wang et al. 2010). Therefore, down-regulation of ACACB gene expression is a suitable approach to reduce the fat deposition by increased fatty acid oxidation. RNAi has been a standard method in both cultured cells and various model organisms for the controlled downregulation of gene expression (Tripathi et al. 2012).
In chicken, Acetyl CoA carboxylase beta (ACACB) gene is involved in fatty acid metabolism and thereby, potentially regulating the quality of meat and egg. It is considered that in vitro knockdown of the ACACB gene in cell culture system by developing a short hairpin RNA (shRNA) expressing construct would help, in devising suitable in vivo strategies for knocking down of the gene. This, in turn, might help to produce a knockdown chicken with reduced fat deposition. The present investigation was designed with the objective of silencing ACACB gene with potential shRNA molecule and its effect on other fat synthesis genes under in vitro cell culture system.

Animals
The experiment was conducted in control broiler chicken line maintained at the Institute farm, ICAR-Directorate of Poultry Research, Hyderabad, India. The fertile eggs were collected from this chicken line and incubated in the incubator for 9 days at 98-100 o C with 78-80% relative humidity and turning 6 times a day for embryo development from which embryonic muscle tissues were collected for chicken embryo myoblast cell culture. The entire study was approved by the Institute Animal Ethics committee (IAEC) and Institute Bio-safety Committee (IBSC) of ICAR-Directorate of Poultry Research, Hyderabad, India. All the bio-safety guidelines of IBSC were followed while conducting the experiments.
Designing and cloning shRNA molecules BLOCK-iT RNAi Designer (https://rnaidesigner.thermo sher.com/rnaiexpress/) is an online tool that was used to build ve separate shRNA sequences that target the ACACB gene's open reading frame (ORF) using Tuschl's motif pattern. The siRNA sequences were then transformed into shRNA with CGAA as the stem-loop sequence and sense loop antisense as the strand orientation. The nal shRNA has a 4-nucleotide 5′ overhang (CACC or AAAA) for directional cloning of the ds oligos encoding the shRNA of interest (Table 1). The IDT manufacturers (Integrated DNA Technologies, USA) synthesized the oligos, and scrambled shRNA oligos (lac z) were supplied by the Invitrogen manufacturers (Invitrogen, USA) used in this experiment. The integrity of the ds oligos was checked by loading 5µl annealed ds oligo (500 nM stock) into 4% agarose gel. Both annealed ds oligos (50 bp) and remaining unannealed ss oligos (25 bp) were observed in the gel electrophoresis (Fig. 1). The annealed oligos were cloned into pENTR™/U6 vector (Invitrogen, USA) according to the manual instructions. The transformed recombinant colonies on LB agar plates were screened by using forward (U6: 5′-GGACTATCATATGCTTACCG-3′) and reverse (M13: 5′-CAGGAA ACAGCTATGAC-3′) primers for checking the presence or absence of shRNA inserts (Fig. 1). The plasmid obtained from each pENTR TM /U6 entry construct was sequenced to con rm the sequence and orientation of the ds oligos insert.  (Fig. S1). Further, secondary structures of chicken ACACB mRNA were predicted using a web server Mfold 2.3 version (http://www.unafold.org/mfold/applications/rna-folding-form-v2.php) (Fig. S2).
The Oligowalk software in the RNA structure version was used to determine the thermodynamic properties regulating each shRNA molecule's binding a nity to its mRNA target region were envisaged (http://rna.urmc.rochester.edu/cgibin/server_exe/oligowalk/ oligowalk_form.cgi) ( Table 2). Hence, the following parameters were taken into account: ΔG overall (Overall Gibbs free energy change): The net energy (ΔG in kcal/mol) resulting due to binding of oligos to the target site when all energy contributions are considered which includes target structure breaking energy and oligo self-structure energy. The more negative value of ΔG indicates a more stable duplex. ΔG duplex : It measures the oligo target binding a nity. A more negative value of ΔG duplex indicates more stability of the duplex and vice versa. ΔG break−target (disruption energy): The energy cost for disrupting base pairs in mRNA target region so that the binding site becomes single-stranded and completely accessible. A more negative value denotes that the binding site is less accessible. ΔG intra−oligomer and ΔG inter−oligomer : The negative ΔG of stable structures was greater than that of unstable structures. The free energy changes/differences caused by unimolecular and bimolecular siRNA foldings. ΔΔG ends : It tests the free energy variations in base pairing between the two ends of the antisense strands in the siRNA duplex, i.e. the 3′ and 5′ ends of the antisense strands, also known as differential stability of siRNA duplex ends (DSSE). The main characteristic of an effective siRNA is the unstable 5′end (End-diff is more positive). Chicken embryo myoblasts (CEM) primary cell culture The 9-day old embryos were used for preparing the chicken embryo myoblast (CEM) primary cell culture (Sato et al. 2006). The collected eggs were cleaned with 70% ethanol and the broad end of the egg was cracked using sterile forceps and peeled off the white shell membrane to reveal the chorioallantoic membrane (CAM) below along with blood vessels. The sterile curved forceps were used for piercing the CAM and gently grasping the embryo under the head and lifted out, and transfer to the sterile 9-cm petri dish containing sterile phosphate buffer saline (PBS) and rinsed thoroughly. The head, limbs, and wings were removed by using scissors, and nally, the ventral side of the embryo was cut open to remove all internal organs and transferred in to a new second petri dish to dissect unwanted tissue like fat or necrotic material, debris and blood tinge was removed by washing with PBS several times. The sterile scissors were used to cut the tissue into ne 3-mm pieces. The minced tissue pieces were transferred to a sterile beaker containing a sterile magnetic bar and 10 ml of 2.5% trypsin, and the beaker was placed on the magnetic stirrer for stirring at about 100 rpm for less than 10 minutes at 37 0 C and it was further ltered through a sterile double-layered muslin cloth into a fresh beaker and the ltrate was centrifuged for 5 minutes at 1500 rpm. After discarding the supernatant the resulting pellet was resuspended in 5 ml of growth medium (DMEM, HiMedia) with fetal bovine serum (FBS) to stop the trypsin action. The cell suspension was diluted to 1x10 6 cells/ml in a growth medium with the help of a hemocytometer and seeded approximately 2x10 5 cells/cm 2 in each 25-cm 2 tissue culture ask. The tissue culture asks were incubated at 37 0 C with 5% CO 2 until a con uent monolayer was obtained.

Transfection of shRNA construct
The recombinant shRNA clone (pENTR TM /U6 Entry vector with respective shRNA insert) against ACACB gene and a control plasmid (plasmid containing scrambled shRNA) was transfected into the chicken embryo myoblast (CEM) primary cell culture using the Gene Pulser Xcell™ Electroporation system (Bio-Rad) to predict the activity of respective shRNAs in myoblasts. Two days before electroporation, the cells were transferred to a new 25-cm 2 tissue culture ask with fresh growth medium (DMEM) supplemented with 10% FBS and antibiotic antimycotic solution. The cells were grown up to the late-log phase, such that there will be 70-80% con uent on the day of the experiment. Adherent cells were rst trypsinized by adding 0.1ml/cm 2 trypsin, it was nally inactivated with complete medium, and the cells were harvested by centrifuging at 1500 rpm for 5 minutes at room temperature. The cell pellet was re-suspended in growth medium yielding an approximate cell concentration of 2.5 X 10 6 cells/ml medium. The transfection reactions of all the shRNAs along with scrambled shRNA were carried out in triplicates.
Approximately 0.4ml of the cell suspension was transferred into a 0.4 cm electroporation cuvette and the puri ed plasmid DNA was added to the cell suspension to a nal concentration of approximately 50µg/ml. The DNA and cell suspension were mixed in the cuvette and placed in the shock pod unit holder in the electroporation apparatus and a single square wave pulse was given at 110V with 25 milliseconds pulse length. After the pulse, cell suspension was transferred into a 25-cm 2 ask containing 5ml of growth medium. The asks were rocked gently to assure even distribution of the cells over the surface of the ask and incubated at 37 0 C in a CO 2 incubator. After 48 hours of transfection, the adherent cells were trypsinized and transferred into sterile 15ml conical tubes. The cells were harvested by centrifuging at 1000 rpm for 10 minutes and stored at -80 0 C until total RNA isolation.
Total RNA isolation and cDNA synthesis The total RNA was isolated from the CEM cell pellet using 1ml of Trizin (GCC Biotech, India), according to the manufacturer's protocol. After homogenization, the sample was incubated for 5 minutes at room temperature and then chloroform (200µl/sample) was added to the sample, shaken vigorously for 15 seconds, and incubated for 5 minutes at room temperature followed by centrifugation at 12000 rpm for 15 minutes at 4 0 C. The upper aqueous phase was transferred to a new 1.5ml RNase free sterile microcentrifuge tube and 500µl of isopropanol was added for precipitation. The tubes were incubated for 30 minutes at -20 0 C and then pelleted by centrifuge at 12000 rpm for 10 minutes at 4 0 C. The pellet was washed with 750µl of 75% ethanol, air-dried the pellet for 5 to 10 minutes, and then dissolved in 20µl of RNase-free water. The total RNA was treated with DNase I (HiMedia, India) to remove any trace amount of genomic DNA. The RNA quality and quantity were checked by 1.2% denatured agarose gel and Jenway™ Genova Nano Micro-volume Spectrophotometer (Fisher Scienti c, USA) respectively. For each sample, 2µg of total RNA was reverse transcribed using the Verso cDNA Synthesis Kit (Thermo Scienti c, USA) in a 20µl reaction using Oligo dT and random primers. The cDNAs were diluted at 1:3 with nuclease-free water before the qPCR analysis.
Primer designing and qRT-PCR A total of 9 genes were selected on the basis of their role in de novo fat synthesis and immune response viz. ACACA (Acetyl-CoA Carboxylase Alpha), ACACB (Acetyl-CoA Carboxylase Beta), FASN (Fatty acid synthase), SCD (Stearoyl-CoA Desaturase), ELOVL2 (ELOVL Fatty Acid Elongase 2), CPT1 (Carnitine palmitoyltransferase 1), IFNA (Interferon alpha), IFNB (Interferon Beta), and BLB1 (Major histocompatibility complex class II beta chain BLB1). The sequences were downloaded from NCBI and CDS region was identi ed by using the ExPASy translation tool (https://web.expasy.org/translate/) and the primers were designed based on CDS using IDT Primer Quest software (https://www.idtdna.com/Primerquest/Home/Index) ( Table 3). The expression levels of target genes (ACACA, ACACB, FASN, SCD, ELOVL2, CPT1, IFNA, IFNB, and BLB1) and reference gene Albumin (ALB) were quanti ed using thermal cycler Himedia Insta Q96™ with Bright Green 2X qPCR Master mix ROX (abm, Canada). The qPCR experiments were performed in a 10µl reaction volume [containing 1µl of diluted cDNA, 5µl of BrightGreen 2× qPCR MasterMix, 0.3µl of each primer] under the following program: an initial denaturation for 5 minutes at 95 0 C, followed by 40 cycles of ampli cation with the 30s of denaturation at 95 0 C and 60s of annealing/extension at 60 0 C. The dissociation curve was obtained by heating the amplicon from 55 to 95 0 C. All qPCR reactions were carried out in three biological replicates.
Non-template controls (NTC) were also included in each run for each gene.

Statistical analysis
The experiments were repeated twice, the relative expression of each gene was calculated by using 2 −ΔΔCt , and statistical analysis was carried out using the trial version of SPSS 25. A univariate general linear model with Tukey's HSD and DMRT as post hoc test was used to study the signi cant difference between different shRNA groups due to the knockdown effect of target genes. Data from representative experiments were presented as Mean ± SE for different samples with differences determined by least signi cant differences at 5% level (p ≤ 0.05).

Secondary structures and thermodynamic properties for shRNAs
The evaluation of the shRNA constructs 1 and 5 were devoid of any secondary structures in their antisense strands while shRNA constructs 2, 3, and 4 formed secondary structures (Fig. S1). Hence, the shRNAs 1 and 5 showed minimum free energy (MFE) of 0.00 kcal/mol whereas, the shRNA 2, 3, and 4 showed MFE of -1.40, 0.30, and 0.10 kcal/mol respectively, due to the secondary structure formation (Fig. S1). The stem loop structures were formed for all the ACACB mRNA target regions of anti-ACACB shRNA molecules (Fig. S2). The GC percent of shRNAs ranged from 43-52% where shRNA 3 had the lowest (43%) of all the shRNAs studied, while shRNA 2 had a higher percentage (52%). The predicted values (negative) about overall/net ΔG value, ΔG duplex, and ΔG break-target were found to be highest in shRNA 5, shRNA 1, and shRNA 2 respectively, while the lowest values were observed in shRNA 2, shRNA 2, and shRNA 5, respectively. The DSSE was found to be highest for shRNA 5 and lowest for shRNA 2 ( Table 2).

Cloning and con rmation of ACACB anti shRNA in pENTR TM /U6 Vector
The RNAi-Ready pENTR TM /U6 vector was used to ligate all the ve shRNA oligos of the ACACB gene and transformed into One Shot TOP10 chemically competent E. coli cells. The recombinant clones of all the shRNAs were con rmed by colony PCR (Fig. 1). Further, the plasmid DNA was isolated from the recombinant clones and con rmed by plasmid PCR (Fig. 1). Each recombinant pENTR TM /U6 construct was sequenced to con rm the sequence of the shRNA, which revealed the absence of mutations. The DNA was isolated from the transfected cell culture of each construct and used as a template in PCR to check for the presence of pENTR TM /U6 Entry vector contains shRNA and found a product of 293 bp length signifying the successful transfection ( Fig. 1).

Silencing e ciency of ACACB shRNAs
The qRT-PCR was performed with ACACB and ALB gene speci c primers with all ve shRNA treated CEM samples and scrambled shRNA treated samples was used as a control. The initiation phase, exponential phase, and plateau phase of ampli cation curves were all optimal, indicating that the product was ampli ed exponentially, i.e. uorescence emission was corresponding to the ampli ed template. Melting curve analysis was performed at the end of the qRT-PCR cycle to verify the speci city of ampli cation, revealed a single peak for all genes, suggesting that the PCR products were homogeneous. The knockdown performance of shRNA 1 and 5 against scrambled shRNA was found to be signi cantly different (P ≤ 0.05) in the qRT-PCR study. The fold change of ACACB gene in different shRNA constructs was 0.33, 0.75, 0.61, 0.62 and 0.31 in the cells with shRNA 1, shRNA 2, shRNA 3, shRNA 4 and shRNA 5 constructs respectively. In contrast to the scrambled shRNA, the knockdown (KD) percent of ACACB mRNA caused by different shRNA ranged from 69% (shRNA 5) to 25% (shRNA 2), respectively (Fig. 2).

Effect of shRNAs on immune response genes
In the knockdown cells, the relative expression of immune response genes such as IFNA, IFNB, and BLB1 was also monitored by qRT-PCR in target as well as control samples. However, when compared with scrambled shRNA, no signi cant difference was observed in all ve shRNA groups (Fig. 3).
Relative quanti cation of de novo fat synthetic genes Once the knockdown of ACACB gene was found, the relative expression of de novo fat synthetic genes such as, ACACA, FASN, SCD, ELOVL2, CPT1 have been studied. The ACACA, FASN, SCD and CPT1 genes were down-regulated and ELOVL2 gene was signi cantly up-regulated. The fold change of ACACA, FASN, SCD, CPT1 and ELOVL2 was 0.55, 0.08, 0.01, 0.05 and 2.87 respectively (Fig. 4).

Quanti cation of cholesterol and triglycerides
The cholesterol and triglycerides were estimated from the control and ACACB knockdown cell lysate and estimated by using Turbo Chem 100 automatic blood analyzer. We found signi cant (P < 0.01) reduction in cholesterol and triglyceride levels in ACACB knockdown cells at 47.71% and 34.91% respectively as compared to the control (Fig. 5).

Synthesis of potential shRNAs
The success of RNAi depends on the designing of the shRNA for speci c target recognition and minimization of off-target effects. In the present study, unique/speci c shRNA molecules were designed based on the Reynolds ranking criteria for the ACACB gene and also, the speci city of the shRNA sequences is important for the formation of RISC (Reynolds et al. 2004; Paddison et al. 2004). In addition to speci city, the G-C content plays an enormous role in the formation of duplex siRNA molecule. Low G-C content is known to decrease a nity for the target sequence, whereas higher G-C content interferes with RISC formation that eventually cleaves the mRNA molecules. The G-C content of the designed shRNA molecules in the present study was moderate for the duplex siRNA formation. In addition to the G-C content, the lack of internal repeats, an A/U-rich 5' end, Tuschl motifs, and other features were included which improves the silencing e ciency of siRNA (Fuchs et al. 2004).
Finally, a simple local alignment search tool (BLAST) was also employed to ensure that shRNA had no signi cant homologies with genes other than the target to avoid possible off-target effects. Further, silencing performance was also positively associated with the siRNA-mRNA duplex (G duplex ) stability (Pascut et al. 2015). Similarly, 1 and 5 shRNAs had higher G duplex values than shRNA 2, 3, and 4 (which had a low G duplex value), meaning that 1 and 5 shRNAs could bind the target site more e ciently. The knockdown e ciency of 1 and 5 shRNAs was high, supporting the above predictions. The duplex asymmetry (DSSE) and target site accessibility could improve knockdown e ciency by about 26% and 40%, respectively (Shao et al. 2007). For enhanced RNAi potency, the siRNAs suitable disruption energy and DSSE is < − 10 kcal/mol and > 0.0 kcal/mol, respectively (Shao et al. 2007). We noticed that the expected disruption energies and DSSE for all of the modeled shRNAs were less than − 10 kcal/mol and greater than 0.0 kcal/mol, respectively, suggesting that all shRNAs have the ability to silence genes. It has been discovered that the development of self-structures in shRNA (both unimolecular and bimolecular shRNA folding) reduces the equilibrium a nity for the target mRNA (Lu and Mathews, 2007). For all of the shRNAs the predicted selfstructure energies in this sample were low, indicating their e ciency in silencing.
Gene silencing caused by shRNAs is mostly due to sequence-speci c mRNA degradation by antisense/guide strand of shRNA (Martinez et al. 2002). The 'guide strand' is inserted into the active RNA-induced silencing complex (RISC) to locate the mRNA, which has a complementary sequence leading to the endonucleolytic cleavage of the target mRNA leading to gene silencing (Hannon, 2002). The degree of secondary structure in the antisense strand was of utmost importance in determining the highly However, shRNAs shorter than 30 bp can evade PKR activation (Robbins et al. 2006) and some experiments exhibited a signi cant increase in the expression of immune response genes including Scramble shRNA. It was also observed that activation of interferon response in goat myoblast cells due to exogenous administration of shRNA against the ACTRIIB gene (Patel et al. 2014). Interferon modulation in chicken embryonic myoblast cells varied between 46-112 folds of OAS1 and 2-7.2 folds for IFNβ compared to mock-transfected control due to anti myostatin shRNA (Tripathi et al. 2013). It has been reported that the shRNA-mediated myostatin knockdown in transgenic sheep showed increased MHCI expression (Hu et al. 2013).
There are several reports regarding shRNA-induced interferon responses suggested that a high level of shRNA expression might be due to the accumulation of unprocessed or aberrantly processed transcripts triggering interferon response (Stewart et al. 2008;Cao et al. 2005;Watanabe et al. 2006). The shRNA transfected porcine embryo cells, showed induction of the OAS1 and IFNβ genes by 1000 and 50 folds respectively (Stewart et al. 2008). Further, reported that the introduction of H1 and U6 promoter-based shRNA constructs by pronuclear microinjection led to induction of OAS1 gene and early embryo lethality (Bridge et al. 2003). In the present experiment, the expression of interferon genes IFN α, IFN β, and BLB1 were analyzed in control (scrambled shRNA treatment) and knockdown cells possessing different shRNA molecules to explore the impact of shRNA on immune function. It was observed that there was no signi cant difference of expression of IFN α, IFN β, BLB1 genes between knockdown and scrambled shRNA treated cells. Hence, shRNA molecules used in the present study have not been captured by the interferon mechanism in vitro as any foreign DNA fragment of a speci c length is normally detected and destroyed by the interferons. It may therefore be construed that these shRNA molecules have been very much effective to silence ACACB expression without interfering with the body's immune system. However, the introduction of short (< 30 nt) dsRNAs with 2-base 3 overhangs resembling dicer processing does not activate the interferon pathway and also shRNA expression from vectors in the nucleus resemble endogenous miRNA (Elbashir et al. 2001 Effect of ACACB gene silencing on de novo fat synthetic genes Now-a-days, researchers are more focused on ACACA and ACACB to understand the chemistry and biological activity because they are very important enzymes in fatty acid synthesis and oxidation. The ACACB gene is localized subcellularly on the mitochondrial membrane and is involved in the synthesis of malonyl-CoA, and this inhibits the CPT1, which plays an important role for controlling the two opposing pathways i.e. fatty acid synthesis and oxidation. In this study, we knocked down the ACACB gene in chicken myoblast cells and tracked the expression of key fatty acid metabolism genes such as ACACA, FASN, SCD, ELOVL2, and CPT1. We observed down-regulation of ACACA, FASN, SCD, and CPT1 gene and up-regulation of ELOVL2 in ACACB knockdown myoblast cells. The down-regulation of ACACA, FASN, SCD genes indicates suppression of the fatty acid synthesis, and up-regulation of ELOVL2 indicates the enhancement of long-chain fatty acids formation (Fig. 6). The CPT1 is a rate-limiting enzyme and down-regulation of this enzyme indicates blocking the β-oxidation for balancing the fatty acid synthesis and oxidation because of fewer fats accumulation in the tissues.

Effect of ACACB gene silencing on cholesterol and triglyceride synthesis
In chicken egg and meat the high cholesterol and triglycerides are most undesired components and which compels us to eat chicken products in less quantity. In conclusion, from the present study we have identi ed potential shRNA molecules against the ACACB gene where shRNA1 and shRNA5 showed more than 60% knockdown e ciency on the expression of the ACACB gene under in vitro myoblast cell culture system. The silencing of the ACACB gene showed to have a direct effect on the down-regulation of ACACA, FASN, SCD, and CPT1 genes, and up-regulation of ELOVL2 gene in myoblast cells. We suggest that these shRNA molecules may be used under in vivo system for the development of knockdown chicken having a potential of producing lean meat by silencing of expression of the ACACB gene. Compliance with ethical standards

Declarations
The entire study was approved by the Institute Animal Ethics committee (IAEC) and Institute Bio-safety Committee (IBSC) of ICAR-Directorate of Poultry Research, Hyderabad, India. All the bio-safety guidelines of IBSC were followed while conducting the experiments.
Con ict of interest The authors have no con ict of interest.    Diagrammatic representation of targeted gene expression analysis conducted in the present study. In chicken tissues, Acetylcoenzyme A carboxylase 1 (ACACA) and acetyl-coenzyme A carboxylase 2 (ACACB) play distinct roles in lipid metabolism. Diet fat, carbohydrate, and protein are digested, and the fatty acids (FA), glucose, and amino acids are transported to various tissues, including liver, adipose, and muscle. In the cytosol acetyl-CoA is carboxylated to malonyl-CoA by ACC1 and utilized through fatty acid synthase (FAS) and Stearoyl CoA Desaturase (SCD) reactions to generate palmitate and palmitoleate, which is utilized in the synthesis of triglycerides (TG) and VLDL. Also, acetyl-CoA is carboxylated by ACC2 at the mitochondrial membrane to form malonyl-CoA, which inhibits the CPT1 and reduces acyl-CoA transfer to mitochondria for β-oxidation. The down and up arrow (↓) indicates the down and up-regulation of genes in ACACB knockdown CEM primary cells compared to the control cells. Mean values were different at *P ≤ 0.05.

Supplementary Files
This is a list of supplementary les associated with this preprint. Click to download. supplementaryFigures.pptx