title: LN-TCV+2L_1
design: the transcriptome sequencing was performed using Illumina Hi Seq TM2500 (Gene Denovo Biotechnology Co. China). RNA expression levels were determined using the fragments per kilobase of transcript per million mapped reads method. The fold-change method was used to identify RNAs that were differentially expressed after castration using the R package DEGseq (R-3.6.2). Genes with a log2 fold change >2 and an adjusted p <0.05 were deemed to be significantly differentially expressed.