Description
Purpose: With the advent of Next-generation sequencing (NGS), several novel genes/proteins and cellular pathways in wide variety of tissues has been discovered. The aim of this study are to perform uterine transcriptome profiling (RNA-seq) to determine differently expressed genes in laying and non-laying hens and to further validate the expression of candidate genes using real-time quantitative reverse transcription polymerase chain reaction (qRT–PCR) in laying, non-laying and molting hens. Methods: Uterine mRNA profiles of 35-60 weeks-old laying and non-laying hens, three each, were generated with NextSeq 500 sequencer in single-end mode with a read length of 1x76 bp. Raw sequencing reads were cleaned and trimmmed with Prinseq tool and good reads were aligned against the chicken reference gemone (Galgal 5.0) in Array Studio. Differential gene expression analysis was performed by the DESeq2 algorithm as implemented in Array Studio. The genes with at least two-fold change (FC) and Benjamini and Hochberg q-value < 0.05 were called differentially expressed. Results: Using an optimized data analysis workflow, we mapped about 32 million reads from layers and 28 million reads from non-layers to the chicken genome. A total of 19,152 gene transcripts were annotated from Ensembl alignment which represents 50.24% of the chicken genome assembly. Differential gene expression analysis showed 616 were differentially expressed between layer and non-layer hens. 229 DEGs were significantly up-regulated and 286 were significantly down-regulated in the laying hens when compared to the non-laying hens. Twelve candiate genes, linked to calcium remodeling, were identified by gene function analysis and validated using qPCR. MEPE, CALCB, OTOP2, STC2 and ATP2C2 were confirmed to be highly expressed in laying hens as compared to molting and non-laying hens. RNA-seq and qPCR data for relative gene expression were highly correlated (R2 =0.99). Conclusions: Our study reports the expression of four novel genes that are speculated to transport calcium ions across the uterine epithellium for eggshell mineralization. These genes can be used as quantitative basis of selecting hens with an improved eggshell quality. Overall design: Uterine mRNA profile of three 35-60 weeks old laying and non-laying hens was genereted using Illumina NextSeq 500