Study location
The present study was conducted at the Mega Seed Project on Pig Unit, Indian Council of Agricultural Research (ICAR), Research Complex for North Eastern Hill Region, Nagaland Centre, Medziphema located at 25o45'24" N, 93o50'26" E of Eastern Himalayan subtropical agro-climatic zone. The maximum air temperature ranged from 19.7 °C to 37.5 °C and the minimum air temperature ranged from 6.1 °C to 27 °C, respectively. The area has a humid to rainy climate and receives heavy rainfall of approximately 1800-2500 mm per year.
Animals and management
The experiment was conducted on Tenyi-vo male pig (n=24) reared under an intensive management system (Supplementary Figure 1). The breeding stock of Tenyi-vo pig were procured from the native tract and maintained in an isomanagerial condition with a commercial concentrate ration. Male piglets were selected at the age of four weeks and maintained with the dam till weaning at the 45 days of age. At post weaning, group feeding of grower ration was practised with water ad libitum. The experiment was performed in compliance with the approved guidelines of the Institute Animal Ethics Committee, ICAR Research Complex for NEH region, Umiam.
Collection of testes
Tenyi-vo male pigs were divided into four groups (n=6/group) viz., G1 (30-45 days), G2 (60-65 days), G3 (80-100 days) and G4 (150-160 days). The body weight of Tenyi-vo boar at birth and weaning was 0.57±0.03 Kg and 3.81±0.21 kg, respectively. Testes were collected by open surgical castration under lignocaine local anesthesia and were transported to laboratory in a sterile thermal flask. After washing with phosphate buffer saline (PBS, pH 7.4), epididymis was dissected. The weight and volume of testis and weight of epididymis were recorded. A piece of testicular tissue was preserved in 10% neutral buffer formalin saline for histology; another piece of tissue weighing about 100 mg was collected in RNAlater (Qiagen, GmbH, Germany) and stored at -20°C for total RNA extraction. On the day of castration, body weight was recorded and blood was collected in heparinised vials to assay for steroid (T and E2) and thyroid (T3 and T4) and cortisol hormones.
Cauda epididymal spermiogram
The cauda epididymis of each testis was incised and the seminal fluid was aspirated into a micro-centrifuge tube. In the absence of free fluid in the cauda, the lumen was flushed with a few drops of PBS and the fluid was allowed to drain into the micro-centrifuge tube under gravitation. The sperm concentration was estimated by measuring the absorbance of diluted semen sample (2.9% sodium citrate buffer as diluent) at 650 nm in a spectrophotometer (Shimadzu, Japan) using the following the formula [concentration = (Dilution Factor) × (21.39 × (Absorbance) - 1.09)]. The viability and morphological abnormality were measured by eosin-nigrosin staining methods (Kondracki et al. 2012).
Histology of testis
Preserved tissue samples were trimmed and dehydrated by treatment with increasing alcohol concentration, followed by clearing with acetone and benzene. The tissues were embedded in paraffin and the serial sections of five µm thickness were cut through a microtome. In the end, tissue sections were stained with haematoxylin and eosin and examined under a light microscope.
Hormonal profile
Plasma T (Cat #AA-E1300), E2 (Cat # FR E-2000), T3 (Cat # TF E-2100), T4 (Cat # TF E-2200) and cortisol (Cat # MS E 5000) concentrations were assayed using commercial ELISA kits (LDN Diagnostica, Germany). The hormone concentration of the test samples was determined using four parameter logistic equations derived from the standard curve in the Multiskan ELISA reader (Thermo scientific Inc., USA).
Differential gene expression profile in testicular tissues
Total RNA was extracted using the RNA isolation kit (Cat # KT151A, RaflexTM Genei, India) and the concentration and purity of RNA was calculated using a micro-cuvet in a UV-VIS spectrophotometer (Shimadzu, Japan). RNA samples having A260/280 between1.8 to 2.0 were considered for cDNA synthesis. DNase I (1U/μg of RNA, Cat # EN0521, Thermo-Fischer scientific Inc. USA) treatment was performed to eliminate gDNA from the total RNA extracted. cDNA was synthesised using RevertAid M-MuLV reverse transcriptase (Cat # EP0442, Thermo Scientific, USA) and oligo(dT)18 primers (Cat # SO132, Thermo Scientific, USA). A reverse transcriptase control was run in parallel to validate the effectiveness of the DNAse I during cDNA sysnthesis. The primers for AR, CYP19A1, ERα and ERβ, IGF-1R, and endogenous control gene (GAPDH) were designed using the Integrated DNA Technology primer quest tool and custom synthesized (Europhin, India, Table 1). The PCR cyclic condition of the target genes was optimized using Taq PCR Master Mix (Cat # K0171, Thermo-Fischer scientific Inc, USA) in a gradient thermal cycler (Nexus gradient Master cycler, Eppendorf, Germany) and the amplicon size was verified by 1.5 per cent (w/v) agarose gel electrophoresis.
Each target gene was amplified using Maxima SYBR Green qPCR master mix (Cat# K0251, Fermentas, Thermo scientific, USA) in a real-time qPCR system (Pico Real 96, Thermo Scientific, USA). cDNA 10 ng (1 μL) templates were added to 0.2 μL forward and reverse primers (10 pmol each) and 5 μL 2X SYBR Green Master Mix at a final volume of 10 μL. The thermal cyclic conditions included hot start denaturation at 95°C for 15 min, 40 cycles of three segmented amplification and quantification programmes (denaturation at 95°C for 15 s, annealing at a specific temperature for 30 s, an extension at 72°C for 30 s and a melting stage by heating from 55 to 95°C at a rate 1°C/s with the acquisition of fluorescent data. No template control was used as a negative control in the qPCR reaction. GAPDH was used as an endogenous control to generate ∆Ct value of the target genes, because its expression was consistent across the age group. The intra-assay coefficient of variation of each target gene at testis was less than 3%. Gr. 1 pigs served as calibrator to generate ∆∆Ct and the fold change of each gene was determined by 2-∆∆Ct method (Livak and Schmittgen 2001).
Statistical analysis
The variables such as testicular weight, volume, epididymal weight, cauda epididymal spermiogram and hormone were analyzed by One-way ANOVA with Tukey’s post-hoc test. ΔCt values of each target gene in G1, G2, G3 and G4 was analyzed by non-parametric Kruskal Wallis test with Dunn’s post-hoc test. Significance was set at 95% and the results are presented as mean ±SEM. SPSS 16.0 was used for data analysis and graphpad prism 6.0 was used for constructing the bar diagrams.