2.1. Samples
2.1.1. Cows
A total of 514 Chinese Holstein cows used for the association study were randomly collected from six farms in the North of China (Figure 1). The parities of the cows ranged from 1 to 4 and were milked three times a day, the lactation days were among 30-400 days, and SCS varied from 0 to 10. All the dairy cows were fed on the same lactation diet according to the recommendations of lactating Chinese Holsteins.
In addition, 123 healthy perinatal Holstein cows (from a single dairy farm in Hebei Province, China) were randomly selected for the folic acid supplementation experiment. The parities of the cows were 1 or 2, the estimated body weight of each cow was shown in S1 Table, and the expected date of calving of the cows were within 10 days. Accordingly, they were grouped into three experimental groups (groups A, B, and C) supplemented with different doses of folic acid (Group A: 0 mg/day; Group B: body weight (kg) × 0.24 mg/kg/day; Group C: body weight (kg) × 0.48 mg/kg/day). In this study, folic acid was coated by Oriental Kingherb Company (Beijing, China). The coated folic acid was fed to groups B and C for 14 days before calving and 7 days after calving.
Milk samples were aseptically collected and immediately sent to the official Dairy Center of China (Beijing, China) for the detection of somatic cell count (SCC). Peripheral blood samples were collected from the cows’ tail vein by caudal venipuncture and divided into two tubes. One tube containing coagulant and the other containing potassium ethylenediaminetetraacetic acid (EDTA K3E 15%, 0.12mL; Becton, Dickinson and Company) were used for DNA and RNA extraction respectively. Sera were isolated from the peripheral blood and sent to Beijing Huaying Biotechnology Research Institute for cytokine (NIBP, NF-κB, TNF-α, IL17, IL4, IFN- γ, IL6, and IL10) tests detected by radioimmune outfit, which was performed by R-911-automatic radioimmunoassay count instrument.
2.1.2. Bacterial strains
S. aureus (strain 90-1) was originally isolated from a subclinical case of bovine mastitis by our labs in Beijing, China. MRSA (strain W18), provided by Prof. Xin Wang (Northwest A&F University, Yangling, Shanxi, China), was isolated from a clinical case of bovine mastitis. S.aureus and MRSA (50 μL) were inoculated into 5 mL of tryptone soya broth (Beijing Land Bridged Technology Ltd.), and allowed to grow overnight for 24 h at 37 °C, 200 rpm/min. The concentrations of S. aureus and MRSA were determined by serial dilution and standard plate counting. In Brief, S. aureus and MRSA were diluted six gradients successively with DPBS. Then, 100 μL of the diluent was transferred into plate count agar (PCA, Beijing Land Bridged Technology Ltd) and spread with a glass spreader. Subsequently, the agar plates were incubated at 37 °C for 18–24 h. After the PCA plate culture of the S. aureus and MRSA diluents, the numbers of S. aureus and MRSA colonies were counted. Each diluent PCA was conducted in triplicate. Finally, the bacteria were diluted in DMEM to obtain 1×108 CFU/mL.
2.1.3. Cell culture
The bovine mammary epithelial cell line (Mac-T) was provided by Zhejiang University. Mac-T cells were re-suspended in the warm growth medium of DMEM with Glutamax (Gibco) supplemented with 10% fetal bovine serum and 100 U/mL penicillin and streptomycin (100 mg/mL, Gibco). Mac-T cells were seeded in a 25 cm2 tissue culture flask at 37 °C in a 5% CO2 humidified incubator, and the medium was changed once every 24 h. After 48 h, cells at 85% confluence were split by adding 1 mL of 0.25% trypsin/EDTA (Gibco) after washing with 2 mL of DPBS (Gibco). Mac-T cells were cultured up to a maximum of three passages to reduce the risk of aberrant expression caused by extended culturing. These cells were centrifuged in 5 mL of DMEM growth medium for 5 min at 1000 rpm/min, seeded at 5×105 cells in a 6-well cell culture plate (Corning), and then allowed to grow in a growth medium at 37 °C in 5% CO2 humidified incubator.
2.2. Total RNA extraction, reverse transcription and quantitative Real-Time PCR
Buffy coat of fresh blood and collected Mac-T cells were placed directly in 1 mL of Trizol reagent (Invitrogen), and total RNA was extracted following the manufacturer’s instructions. The mRNA was converted to cDNA with the PrimeScriptTM RT reagent kit (Takara). RT-qPCR was performed using the SYBR Green I Master kit (Roche Diagnostics GmbH) on the LightCycler® 480 II (Roche Diagnostics Ltd) following the manufacturer’s instructions. RT-qPCR reaction was performed using the following program: 95 °C for 10 min, followed by 45 cycles at 95 °C for 10 s, 60 °C for 10 s and 72 °C for 10 s. At the end of each run, a dissociation melt curve was determined. All melt curves showed a single peak and were consistent with the presence of a single amplicon. Bovine GAPDH was used as the house-keeping gene to normalize the target genes’ expression levels. Each cDNA sample was analyzed in duplicate. The relative mRNA expression levels of the target genes were calculated using the 2-△△Ct method. Primer sequences are as follows (obtained from Tsingke): TRAPPC9-forward: CTGCTCCGCTCGGTGAATGAC, TRAPPC9-Reverse: GCTTTACCGCCAGTTCCACCA; IL-1β-Forward: CTGTCGGACCCATATGAGC, IL-1β-Reverse: GCTCATGGAGAATATCACTTGTTG; IL-6-Forward: TCCTGAGAAACCTTGAGAAT, IL-6-Reverse: ATAAGTTGTGTGCCCAGTGG; IL-8-Forward: TGAAGCTGCAGTTCTGTCAAG, IL-8-Reverse: TTCTGCACCCACTTTTCCTTGG; GAPDH-Forward: GGTGCTGAGTATGTGGTGGA, GAPDH-Reverse: GGCATTGCTGACAATCTTGA.
2.3. DNA extraction and genotyping of the bovine samples
The genomic DNA was extracted from the whole blood of the dairy cows using the DP-318 Blood DNA kit (Tiangen Biotech Co) following the manufacturer’s instructions. DNA was assessed with the NanoDropTM ND-2000c Spectrophotometer (Thermo Scientific) and by 1% agarose gel electrophoresis.
The genotyping of the bovine samples consists of two steps. First, SNPs identification was conducted by DNA pooling and sequencing. Briefly, a total of 30 samples were randomly selected from 514 cows to construct a DNA pool with an equal DNA concentration of 50 ng/μL for each sample. PCRs were performed in a 25 μL volume containing 12.5 μL of PreMIX (Takara Biotechnology Co), 1 μL of genomic DNA (50 ng), 9.5 μL of ddH2O, and 1μL of each primer (10 pmol/μL). Cycling reaction conditions were as follows: initial denaturation at 95 °C for 5 min; followed by 34 cycles of 30 s at 95 °C, annealing from 52 °C to 62 °C for 35 s, and 72 °C for 1 min; and a final extension at 72 °C for 10 min. The PCR products were sequenced using the ABI3730XL DNA analyzer (Applied Biosystems). In total, four SNPs were identified in the new population.
Second, the identified SNPs were genotyped for 514 cows through the multiplex Snapshot method (Applied Biosystems). The primer sequences of the four identified SNPs were as follows (obtained from Tsingke Inc.): TRAPPC9-SNP1-Forward: AGGAAAGTGAAGTCGCTCAGT, TRAPPC9-SNP1-Reverse: CCTTTTATTGGAGGTCTATGTG; TRAPPC9-SNP2-Forward: TTACATTAACCTTCGGCAGGT, TRAPPC9-SNP2-Reverse: GACTCCAGGGTCTTGACCACA; TRAPPC9-SNP3-Forward: GCTTCAGTCTTCTGGACCCTCT, TRAPPC9-SNP3-Reverse: GATACCCCTGGTCCCCTTTAG and TRAPPC9-SNP4-Forward: TTTATCCTGACGATGTCTGCC, TRAPPC9-SNP4-Reverse: CTGCTGTGAGCCCAAAACTAT.
For SNP1, SNP2, SNP3 and SNP4, 492, 511, 492 and 501 samples were detected successfully (S2 Table). And Haploview (Version 4.1) was used for linkage disequilibrium analysis of identified SNPs.
2.4. DNA bisulfite treatment and hot-start PCR
The EZ DNA Methylation Golden kit was used for the sodium bisulfite conversion of genomic DNA following the manufacturer’s instructions (ZYMO Research). The PCR and sequencing primers of the TRAPPC9 gene used for DNA methylation detection were designed with Oligo 6.0 software and PSQ Assay Design software (Qiagen). Hot start PCR was carried out in a 40 μL system including 20 μL of hot-start PCR premix (ZYMO Research), 1 μL of forward primer (10 pmol/μL), 0.1 μL of reverse primer with the universal tail (10 pmol/μL, [19]), 1 μL of biotin-labeled universal primer (10 pmol/μL), and 4 μL of bisulfite-treated DNA. The PCR conditions were 95 °C for 15 min, followed 94 °C for 30 s, 56 °C for 30 s and 72 °C for 30 s for 45 cycles, and 72 °C for 10 min. The PCR products were detected using 2% agarose gel with ethidium bromide [19].
2.5. Quantitative DNA methylation analyses by pyrosequencing
The pyrosequencing assays were used to analyze the promoter methylation level of the TRAPPC9 gene quantitatively. Six CpG sites located in the promoter region of the bovine TRAPPC9 gene were tested to analyze the DNA methylation level. The Pyro Q-CpG system (Qiagen) was used to analyze DNA methylation following the manufacturer protocol. In brief, bisulfite PCR products were bound with Streptavidin Sepharose High Performance (GE Healthy care). The Sepharose beads grasping the biotin-labeled PCR products were purified in 70% ethanol for 5 s, denatured in a denaturing buffer for 5 s, and then washed in washing buffer for 15 s using the Pyrosequencing Vacuum Prep Tool (Qiagen). Subsequently, a 0.5 μM pyrosequencing primer was annealed to the purified single-stranded PCR product with annealing buffer (Qiagen). The levels of CpG methylation were expressed as the percentage mC/ (mC+C), and non-CpG cytosine residues were used as internal controls to verify bisulfite conversion. The primers measuring the methylation levels of bovine TRAPPC9 are as follows (obtained from Tsingke): R1 Region-Forward: GGGAAGAGTATAGATAATAGTTAGATAGT, R1 Region-Reverse: GGGACACCGCTGATCGTTTAAACAATCCACTCCAATTACTATTACC, and R1 Region-Sequencing: TTTTAAAGGAAAGGAAAGTG. R2 Region-Forward: GAGTTTGGAGTGGTTTTTTTAGG, R2 Region-Reverse: GGGACACCGCTGATCGTTTAAAATAAATTCCCCTTTTACTATC, and R2 Region-Sequencing: GTTTGGAGTGGTTTTTTTAGG.
2.6. Invasion assay of MRSA and S. aureus
The gentamicin protection assay was performed to analyze the invasion ability of MRSA and S. aureus into Mac-T cells as described previously [20]. In brief, Mac-T cells at 85% confluence in 6-well plates were cultured in modified growth medium at 37 °C in 5% CO2. The monolayers were stimulated with MRSA and S. aureus (1×108 CFU/mL) for 3, 6, 8 and 10 hr at a multiplicity of infections of 10 bacteria to 1 Mac-T cell (MOI = 10:1). Then, the Mac-T cells were washed three times with PBS and incubated in DMEM supplemented with 100 μg/mL gentamycin for 2 h in 5% CO2 at 37 °C without serum. Each experimental treatment was conducted in triplicate. These cells were washed three times in PBS and collected at different time points post-inoculation using 0.25% trypsin/EDTA for RNA extraction. The number of CFU/mL was determined by the qPCR technique.
2.7. Challenging Mac-T cells with S. aureus and MRSA strains
Mac-T cells at 85% confluence in 6-well plates were cultured in modified growth medium at 37 °C in 5% CO2. The monolayers were stimulated with S. aureus (90-1) or MRSA (W18) (1×108 CFU/mL) for 3, 6, 12 and 24 h at a multiplicity of infections of 10 bacteria to 1 Mac-T cell (MOI = 10:1). The Mac-T cells infected with bacteria alone and cultured in the same volume of DMEM growth media without any treatment served as controls. Each experimental treatment was conducted in triplicate. These cells were washed three times in PBS and collected at different time points post-inoculation for RNA extraction.
2.8. Cell viability assay
Cell viability was measured using the One Solution Cell Proliferation assay (Promega Corporation) to investigate the effect of folic acid on cell growth. In brief, the Mac-T cells were seeded at 5×104 cells/mL in 96-well plates at 37 °C in 5% CO2 for 24 h, and the cell medium included different concentrations of folic acid (0, 5, 10, 20, and 40 μg/mL). After 24 h, the cells were treated with 20 μL of MTS solutions/well at 37 °C with 5% CO2 for 4 hr. The Mac-T cells cultured in the medium alone and the medium without cells in sextuplicate served as zero adjustment controls. Optical density was measured at 490 nm on a microplate spectrophotometer (Tecan).
2.9. Mac-T cells treated with folic acid and followed bacterial infection
Mac-T cells at 85% confluence in 6-well plates were cultured in modified growth medium at 37 °C in 5% CO2. The monolayers were treated with folic acid at different doses (5, 10, 20, and 40 μg/mL) for 2 hr, and then stimulated with S. aureus (90-1) or MRSA (W18) (1×108 CFU/mL) for 6 hr at a multiplicity of infections of 10 bacteria to 1 Mac-T cell (MOI = 10:1). The Mac-T cells infected with S. aureus or MRSA alone and cultured in the same volume of DMEM growth media without any treatment served as controls. Each experimental treatment was conducted in triplicate. These cells were washed three times in PBS and collected at 6 hr post-inoculation for RNA extraction.
2.10. Cell transfection and RNA interference (RNAi)
For cell-based functional assays, cell transfections were conducted with Lipofectamine 2000 (Invitrogen) on the RNA interference method following the manufacturer’s manual. In brief, the Mac-T cells were seeded at 5×105 cells/mL in 6-well plates for 24 hr, and then treated with DMEM containing extra folic acid (5 μg/mL) (treated group) or only DMEM (control group) for 2 h. For RNAi, cells were co-transfected with synthetic small interfering RNA (siRNA, 100 pmol), Lipofectamine 2000 (5 μL) supplemented with warm DMEM containing 10% FBS in a humidified atmosphere with 5% CO2 at 37 °C at 2 h post-inoculation. The Mac-T cells were treated with folic acid and then RNAi (TRAPPC9), because the purpose of this study was to investigate the protective effect of folic acid on the expression of inflammatory cytokines (indicators of inflammation). The cells were stimulated with S. aureus or MRSA for 6 hr with MOI = 10:1 at 48 h post RNAi transfection. Each experimental treatment was performed in triplicate. The Mac-T cells were washed three times in PBS and collected at 6 hr post-inoculation for RNA extraction. The primers used for bovine RNAi are as follows (designed and synthesized by GenePharma): TRAPPC9 siRNA-Forward: GCAUGGAAGCGUCGGAAUUTT, and TRAPPC9 siRNA-Reverse: AAUUCCGACGCUUCCAUGCTT. Negative control siRNA-Forward: UUCUCCGAACGUGUCACGUTT, and Negative control siRNA-Reverse: ACGUGACACGUUCGGAGAATT.
2.11. Total protein extraction and Western Blot analysis
Mac-T cells were stimulated for 6 hr with S. aureus at 72 hr post-transfection and 2 hr pre-treated with folic acid (5 μg/mL). Finally, these cells were harvested with 500 μL of 0.25% trypsin/EDTA and washed twice in PBS for protein extraction. Total proteins from the cells were extracted by mammalian protein extraction reagent. The cellular proteins were separated by 10% sodium dodecyl sulfate -polyacrylamide gel electrophoresis (SDS-PAGE) and transferred onto a polyvinylidene fluoride membrane (Millipore, Bedford, USA). The membrane was blocked for 2 hr with 5% non-fat dry milk in Tris-buffered saline-Tween (TBST) at room temperature on a rotary shaker and then washed three times (10 min each) with TBST. The membrane was blotted with 1:500 specific primary anti-bovine antibodies followed by 1:5000 secondary antibody conjugated with horseradish peroxidase at room temperature for 1 hr. The membrane was again washed three times (10 min each), and blots were generated with the ECL Plus Western Blotting Detection System (Amersham Life Science). Rabbit anti-GAPDH polyclonal antibody (Good Here) was used as the loading control.
2.12. Statistical analysis
The number (n) of individual animals used per group is described in each figure panel or figure legend. The replicates of experiments are described in each figure legend. The data were analyzed using SAS 9.2. For a single comparison between two groups, t-test was used, and for multiple groups comparison, ANOVA followed by Tukey’s studentized range multiple tests (Tukey) was used. Graphs show mean ± SEM. The upper threshold for statistical significance for all experiments was set at P < 0.05. In summary, ANOVA followed by Tukey’s studentized range multiple tests (Tukey) were applied in the data of Figures 1A, 1B, 4E, 5B, 5C, 5E and S1C; and t-tests were applied in Figures 1C (left panel), 1F, 2C, 3B, 3D, 3E, 4A, 4C, 4D, S1A.
Association study of these four SNPs in TRAPPC9 with milk somatic cell score (SCS) and the concentration of NF-κB, IL-17, IL-6, IL-4, TNF-α and IFN-γ were analyzed using general linear model (SAS 9.0):
y = μ + hys + p + m + h + g + e
where y is the phenotype value, μ is the overall mean, hys is the effect of herd-year-season, p is the effect of parity, m is the effect of lactation stages, h is the effect of healthy status, g is the effect of SNP, and e is random residuals.
Healthy status was divided into three grades: Grade 1, health (SCC ≤ 20 × 104/mL), Grade 2, subclinical mastitis (20 × 104/mL < SCC < 50 × 104/mL), and Grade 3, clinical mastitis (SCC > 50 × 104/mL). Healthy status was not considered in the model when conducting association analyses between the SCS and SNPs of TRAPPC9.