Transcriptomes of bovine ovarian follicular and luteal cells

Affymetrix Bovine GeneChip® Gene 1.0 ST Array RNA expression analysis was performed on four somatic ovarian cell types: the granulosa cells (GCs) and theca cells (TCs) of the dominant follicle and the large luteal cells (LLCs) and small luteal cells (SLCs) of the corpus luteum. The normalized linear microarray data was deposited to the NCBI GEO repository (GSE83524). Subsequent ANOVA determined genes that were enriched (≥2 fold more) or decreased (≤−2 fold less) in one cell type compared to all three other cell types, and these analyzed and filtered datasets are presented as tables. Genes that were shared in enriched expression in both follicular cell types (GCs and TCs) or in both luteal cells types (LLCs and SLCs) are also reported in tables. The standard deviation of the analyzed array data in relation to the log of the expression values is shown as a figure. These data have been further analyzed and interpreted in the companion article “Gene expression profiling of ovarian follicular and luteal cells provides insight into cellular identities and functions” (Romereim et al., 2017) [1].


a b s t r a c t
Affymetrix Bovine GeneChip s Gene 1.0 ST Array RNA expression analysis was performed on four somatic ovarian cell types: the granulosa cells (GCs) and theca cells (TCs) of the dominant follicle and the large luteal cells (LLCs) and small luteal cells (SLCs) of the corpus luteum. The normalized linear microarray data was deposited to the NCBI GEO repository (GSE83524). Subsequent ANOVA determined genes that were enriched (Z2 fold more) or decreased (r À 2 fold less) in one cell type compared to all three other cell types, and these analyzed and filtered datasets are presented as tables. Genes that were shared in enriched expression in both follicular cell types (GCs and TCs) or in both luteal cells types (LLCs and SLCs) are also reported in tables. The standard deviation of the analyzed array data in relation to the log of the expression values is shown as a figure. These data have been further analyzed and interpreted in the companion article "Gene expression profiling of

Value of the data
The follicular GC and TC transcriptome data can be compared to previously published bovine gene expression analyses for corroboration [2,3] and are valuable in metadata analyses investigating GC and TC transcriptomes at different stages or in different species.
There are no previously published transcriptomes available for large luteal cells (LLC) and small luteal cells (SLC). Therefore, the luteal cell gene expression data allow novel insight into these two cell types.
Lists of identified genes that are specifically enriched in each somatic ovarian cell type can inform future physiological research on the functions of the ovarian somatic cells.

Data
The normalized linear microarray data are available at the NCBI GEO repository: GSE83524 at http://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc ¼GSE83524 All tables listed below contain the Affymetrix probeset ID, gene symbol, description of the gene, fold changes, and relevant normalized linear microarray data. There are four microarray replicates for the GC cell type and three replicates for the TCs, LLCs, and SLCs. RNA targets detected at a lower intensity (suggesting low levels of expression) have a slightly higher standard deviation and a poor signal to noise ratio. The overall average standard deviation, however, is low.

Experimental design, materials and methods
Follicular granulosa (n ¼4) and theca cells (n ¼ 3) were isolated from estrogen active dominant follicles in ovaries of beef cows (75% Red Angus, 25% MARC III) from the physiology herd located at the University of Nebraska Agricultural Research and Development Center. The University of Nebraska-Lincoln Institutional Animal Care and Use Committee approved all procedures and facilities used in this experiment. Estrous cycles of cows were synchronized a modified Co-Synch protocol using gonadotropin releasing hormone (GnRH) and a controlled internal drug release device (CIDR; 1.38 g progesterone, Zoetis) for 7 days with a PGF2α (25 mg/mL; Lutalyse, Pfizer Animal Health) injection at CIDR removal [4]. Ovariectomy was performed approximately 36 h after CIDR removal. Upon ovariectomy, each dominant antral follicle was aspirated/dissected and the granulosa cells, theca cells, and follicular fluid were isolated as described previously [4]. Both granulosa and theca cells were homogenized in Tri-reagent (Sigma-Aldrich) for RNA isolation.
Luteal cells were isolated by elutriation from the corpora lutea of ovaries collected at a local abattoir (JBSSA, Omaha, NE) as described previously [5,6]. Each corpus luteum (n¼ 3 for both LLC and SLC) was digested with collagenase to dissociate the cells, which were then suspended in a solution of DMEM (calcium free, 3.0 g/L glucose, antibiotics, 25 mM HEPES, 3.8 g/L sodium bicarbonate; 0.1% BSA, 0.02 mg/mL deoxyribonuclease, pH 7.2) to a total volume of 30 mL. Elutriation was then performed with a Beckman Coulter Avanti-J20X centrifuge with a JE 5.0 rotor. The cells were applied to a Sanderson (Beckman) elutriation chamber and the eluates were collected with continuous flow with each fraction comprising 100 mL of eluate. The first cell fraction was collected at a flow rate of 16 mL/min with a centrifuge speed of 1800 rpm and contained predominantly red blood cells. The second collected fraction was collected at a flow rate of 16 mL/min and a rotor speed of 1400 rpm and contained mostly (SLCs). The third fraction was collected at 24 mL/min at a rotor speed of 1200 rpm and contained predominantly SLCs. The fourth fraction was collected at a flow rate of 30 mL/min and a speed of 680 rpm and contained a majority of large luteal cells (Z80%) along with some SLC and endothelial cells. SLCs and LLCs were concentrated with additional centrifugation of the relevant fractions, and the cells were homogenized in Tri-reagent (Sigma-Aldrich) for RNA isolation.
After RNA extraction, 200 ng RNA for each sample were submitted to the University of Nebraska Microarray Core facility where the Affymetrix Bovine GeneChip s Gene 1.0 ST Array RNA expression analysis was performed. The microarray data were normalized with Robust Multi-Array Averaging. Array analysis was then performed using the National Institute of Aging tool (http://lgsun.grc.nia.nih. gov/ANOVA/) for Analysis of Variance (ANOVA). Fellowship 2014-67011-22280 to HT. USDA is an equal opportunity provider and employer. Names are necessary to report factually on available data; however, the USDA neither guarantees nor warrants the standard of the product, and the use of names by the USDA implies no approval of the product to the exclusion of others that may also be suitable. The research was also supported by funds from the VA Research and Development service.