Transcriptome profiling identifies differentially expressed genes in Huoyan goose ovaries between the laying period and ceased period

Molecular and transcriptomic analysis of the ovary during laying and brooding stages in Zhedong white geese (Anser cygnoides domesticus)

1. This study investigates the molecular mechanisms affecting brooding in Zhedong white geese by examining differences in reproductive endocrine levels, ovarian histology and transcriptomics.2. Twenty 18-month-old Zhedong white geese were selected to examine their ovaries using histological, biochemical, molecular biological, and high-throughput sequencing techniques during the laying and brooding periods.3. The results showed that the number of atretic follicles and apoptotic cells in the ovaries increased significantly (p < 0.05), the levels of follicle-stimulating hormone, luteinising hormone, gonadotropin-releasing hormone and oestradiol decreased significantly (p < 0.05), and the level of prolactin increased significantly (p < 0.01) during the brooding stage.4. In broody geese, the expression of CASP3, CASP9, P53, BAX, and Cyt-c were considerably higher (p < 0.05), but BCL2 expression was significantly lower (p < 0.05).5. In ovarian tissues, 260 differentially expressed lncRNAs, 13 differentially expressed miRNA and 60 differentially expressed mRNA were all discovered using transcriptome sequencing analysis. Functional enrichment analysis revealed that the differentially expressed mRNA and non-coding RNA target genes were primarily involved in ECM-receptor interaction, cell adhesion, cardiac muscle contraction, mTOR signalling, and the calcium signalling pathway.6. In conclusion, follicular atrophy and apoptosis occurred in the ovaries and serum reproductive hormone levels were significantly changed during the brooding period of Zhedong white geese. COL3A1, COL1A2, GRIA1, RNF152, miR-192, and miR-194 may be important candidates for the regulation of brooding behaviour, with the mTOR signalling pathway playing a key role.

Integrative analysis of miRNA-mRNA expression in the brain during high temperature-induced masculinization of female Nile tilapia (Oreochromis niloticus)

Temperature is one of the most important non-genetic sex differentiation factors for fish. The technique of high temperature-induced sex reversal is commonly used in Nile tilapia (Oreochromis niloticus) culture, although the molecular regulatory mechanisms involved in this process remain unclear. The brain is an essential organ for the regulation of neural signals involved in germ cell differentiation and gonad development. To investigate the regulatory roles of miRNAs-mRNAs in the conversion of female to male Nile tilapia gender under high-temperature stress, we compared RNA-Seq data from brain tissues between a control group (28 °C) and a high temperature-treated group (36 °C). The result showed that a total of 123,432,984 miRNA valid reads, 288,202,524 mRNA clean reads, 1128 miRNAs, and 32,918 mRNAs were obtained. Among them, there were 222 significant differentially expressed miRNAs (DE miRNAs) and 810 differentially expressed mRNAs (DE mRNAs) between the two groups. Eight DE miRNAs and eight DE mRNAs were randomly selected, and their expression patterns were validated by qRT-PCR. The miRNA-mRNA co-expression network demonstrated that 40 DE miRNAs targeted 136 protein-coding genes. Functional enrichment analysis demonstrated that these genes were involved in several gonadal differentiation pathways, including the oocyte meiosis signaling pathway, progesterone-mediated oocyte maturation signaling pathway, cell cycle signaling pathway and GnRH signaling pathway. Then, an interaction network was constructed for 8 miRNAs (mir-137-5p, let-7d, mir-1388-5p, mir-124-4-5p, mir-1306, mir-99, mir-130b and mir-21) and 10 mRNAs (smc1al, itpr2, mapk1, ints8, cpeb1b, bub1, fbxo5, mmp14b, cdk1 and hrasb) involved in the oocyte meiosis signaling pathway. These findings provide novel information about the mechanisms underlying miRNA-mediated sex reversal in female Nile tilapia.

Transcriptome analysis of hypothalamus and pituitary tissues reveals genetic mechanisms associated with high egg production rates in Changshun green-shell laying hens

BACKGROUND

Changshun green-shell laying hens are unique to the Guizhou Province, China, and have high egg quality but relatively low yield. Egg production traits are regulated by the hypothalamus-pituitary-ovary axis. However, the underlying mechanism remains unclear. Thus, we conducted RNA sequencing of hypothalamic and pituitary tissues from low- and high-yielding Changshun green-shell laying hens to identify critical pathways and candidate genes involved in controlling the egg production rate.

RESULTS

More than 39 million clean reads per sample were obtained, and more than 82% were mapped to the Gallus gallus genome. Further analysis identified 1,817 and 1,171 differentially expressed genes (DEGs) in the hypothalamus and pituitary, respectively. Nineteen DEGs were upregulated in both the hypothalamus and pituitary of high-yielding chickens. The functions of these DEGs were mainly associated with ion transport or signal transduction. Gene set enrichment analysis revealed that the pathways enriched in the hypothalamus were mainly associated with gonadotropin-releasing hormone (GnRH) secretion, neurotransmitter release, and circadian rhythms. The pathways enriched in the pituitary were mainly associated with GnRH secretion, energy metabolism, and signal transduction. Five and four DEGs in the hypothalamus and pituitary, respectively, were selected randomly for qRT-PCR analysis. The expression trends determined via qRT-PCR were consistent with the RNA-seq results.

CONCLUSIONS

The current study identified 19 DEGs upregulated in both the hypothalamus and pituitary gland, which could provide an important reference for further studies on the molecular mechanisms underlying egg production in Changshun green-shell laying hens. In addition, enrichment analysis showed that GnRH secretion and signal transduction, especially neurotransmitter release, play crucial roles in the regulation of egg production.

DNA methylome and transcriptome identified Key genes and pathways involved in Speckled Eggshell formation in aged laying hens

BACKGROUND

The quality of poultry eggshells is closely related to the profitability of egg production. Eggshell speckles reflect an important quality trait that influences egg appearance and customer preference. However, the mechanism of speckle formation remains poorly understood. In this study, we systematically compared serum immune and antioxidant indices of hens laying speckled and normal eggs. Transcriptome and methylome analyses were used to elucidate the mechanism of eggshell speckle formation.

RESULTS

The results showed that seven differentially expressed genes (DEGs) were identified between the normal and speckle groups. Gene set enrichment analysis (GSEA) revealed that the expressed genes were mainly enriched in the calcium signaling pathway, focal adhesion, and MAPK signaling pathway. Additionally, 282 differentially methylated genes (DMGs) were detected, of which 15 genes were associated with aging, including ARNTL, CAV1, and GCLC. Pathway analysis showed that the DMGs were associated with T cell-mediated immunity, response to oxidative stress, and cellular response to DNA damage stimulus. Integrative analysis of transcriptome and DNA methylation data identified BFSP2 as the only overlapping gene, which was expressed at low levels and hypomethylated in the speckle group.

CONCLUSIONS

Overall, these results indicate that aging- and immune-related genes and pathways play a crucial role in the formation of speckled eggshells, providing useful information for improving eggshell quality.

Embryonic transcriptome unravels mechanisms and pathways underlying embryonic development with respect to muscle growth, egg production, and plumage formation in native and broiler chickens

Background: Muscle development, egg production, and plumage colors are different between native and broiler chickens. The study was designed to investigate why improved Aseel (PD4) is colorful, stronger, and grew slowly compared with the control broiler (CB).

Methods: A microarray was conducted using the 7th-day embryo (7EB) and 18th-day thigh muscle (18TM) of improved Aseel and broiler, respectively. Also, we have selected 24 Gallus gallus candidate reference genes from NCBI, and total RNA was isolated from the broiler, improved Aseel embryo tissues, and their expression profiles were studied by real-time quantitative PCR (qPCR). Furthermore, microarray data were validated with qPCR using improved Aseel and broiler embryo tissues.

Results: In the differential transcripts screening, all the transcripts obtained by microarray of slow and fast growth groups were screened by fold change ≥ 1 and false discovery rate (FDR) ≤ 0.05. In total, 8,069 transcripts were differentially expressed between the 7EB and 18TM of PD4 compared to the CB. A further analysis showed that a high number of transcripts are differentially regulated in the 7EB of PD4 (6,896) and fewer transcripts are differentially regulated (1,173) in the 18TM of PD4 compared to the CB. On the 7th^- and 18th-day PD4 embryos, 3,890, 3,006, 745, and 428 transcripts were up- and downregulated, respectively. The commonly up- and downregulated transcripts are 91 and 44 between the 7th- and 18th-day of embryos. In addition, the best housekeeping gene was identified. Furthermore, we validated the differentially expressed genes (DEGs) related to muscle growth, myostatin signaling and development, and fatty acid metabolism genes in PD4 and CB embryo tissues by qPCR, and the results correlated with microarray expression data.

Conclusion: Our study identified DEGs that regulate the myostatin signaling and differentiation pathway; glycolysis and gluconeogenesis; fatty acid metabolism; Jak-STAT, mTOR, and TGF-β signaling pathways; tryptophan metabolism; and PI3K-Akt signaling pathways in PD4. The results revealed that the gene expression architecture is present in the improved Aseel exhibiting embryo growth that will help improve muscle development, differentiation, egg production, protein synthesis, and plumage formation in PD4 native chickens. Our findings may be used as a model for improving the growth in Aseel as well as optimizing the growth in the broiler.

Transcriptomic analysis of the liver in aged laying hens with different eggshell strength

Eggshell is composed of a very ordered and mineralized structure and is important for egg quality. Eggshell strength is particularly important because of its direct association with economic outcomes and egg safety. Various factors related to laying hens and their environment affects eggshell strength. However, the molecular mechanisms of liver functions related to decreased eggshell strength of aged laying hens are largely unknown. Therefore, this study aimed to identify potential factors affecting eggshell strength in aged laying hens at the hepatic transcriptomic level. A total of five hundred 92-wk-old Hy-line Brown laying hens were screened to select those exhibiting the greatest variation in eggshell strength. Based on the final eggshell strength, 12 hens producing eggs with strong eggshell strength (SES) and weak eggshell strength (WES) were finally selected (n = 6) for liver tissue sampling. The RNA-sequencing was performed to identify differentially expressed genes (DEGs) between the 2 groups. We identified a total of 2,084 DEGs, of which 1,358 genes were upregulated and 726 genes were downregulated in the WES group compared with SES group. According to the Kyoto Encyclopedia of Genes and Genomes pathway enrichment analysis, the DEGs indicated the mammalian target of rapamycin signaling pathway, the Janus kinase-signal transducer and activator of transcription pathway, the mitogen‑activated protein kinase signaling pathway, and the insulin resistance pathways. Genes related to fatty liver disease were upregulated in WES group compared with SES group. In addition, expression of several genes associated with oxidative stress and bone resorption activity was altered in aged laying hens with different eggshell strength. Overall, these findings contribute to the identification of genes involved in different intensity of eggshell strength, enabling more understanding of the hepatic molecular mechanism underlying in decreased eggshell strength of aged laying hens.

MicroRNAs with non-additive expression in the ovary of hybrid hens target genes enriched in key reproductive pathways that may influence heterosis for egg laying traits

Heterosis has been extensively exploited in chicken breeding to improve laying traits in commercial hybrid stock. However, the molecular mechanisms underlying it remains elusive. This study characterizes the miRNAome in the pre-hierarchical follicles of purebred and hybrid laying hens, and investigate the functions of miRNAs with non-additive expression in the pre-hierarchical follicles as they modulate heterosis for egg number and clutch size. To achieve that aim, White Leghorn and Rhode Island Red chicken lines were reciprocally crossed to generate hybrids. The crossbreds demonstrated heterosis for egg number and clutch size, and pre-hierarchical follicles from 4 birds of each genotype were collected at 53 weeks of age. Mode of miRNA expression was characterized after miRNA sequencing. A total of 50 miRNAs including 30 novel ones, were found to exhibit non-additive expression. Dominance was the predominant mode of expression exhibited by majority of the miRNAs. Functional analysis of target genes of the known miRNAs with non-additive expression revealed Gene Ontology terms related to regulation of transcription, metabolic processes and gene expression. KEGG and REACTOME pathways including hedgehog, cellular senescence, wnt, TGF-β, progesterone-mediated oocyte maturation, oocyte meiosis, GnRH signaling, signal transduction and generic transcription, which can be linked to primordial follicle activation, growth and ovulation, were significantly enriched by target genes of miRNAs with non-additive expression. Majority of the genes enriched in these biological pathways were targeted by gga-miR-19a, gga-miR-19b, gga-miR-375, gga-miR-135a, and gga-miR-7 and 7b, thus, revealing their synergistic roles in enhancing processes that could influence heterosis for egg number and clutch size in hybrid hens.

Whole-Genome Sequencing Reveals the Genomic Characteristics and Selection Signatures of Hainan Black Goat

Goats have become one of the most adaptive and important livestock species distributed in developing countries in recent years. The Hainan Black goat is a native goat breed of the Hainan region that is generally well-liked by the local population and is thus raised in large numbers. However, the genomic diversity and selective signals of the Hainan Black goat have not been clearly elucidated yet. Therefore, in this study, we performed whole-genome resequencing of 16 Hainan Black goats and compared the results with those of 71 goats of 6 other breeds from different geographic regions. Principal component analysis (PCA) and phylogenetic analysis identified seven lineages for all goats. Hainan Black goats showed the most similarity with Leizhou goats and the least similarity with Boer goats. Selective sweep analysis identified candidate genes associated with various functions, including immune resistance to disease (TNFAIP2 (TNF alpha induced protein 2) and EXOC3L4 (exocyst complex component 3 like 4)), melanin biosynthetic process (CDH15 (cadherin 15), ASIP (agouti signaling protein), and PARD3 (par-3 family cell polarity regulator)), and light sensitivity (CNGB3 (cyclic nucleotide gated channel subunit beta 3) and CNBD1 (cyclic nucleotide binding domain containing 1)), underlying strong selection signatures in Hainan Black goats. The melanin biosynthetic process, circadian entrainment, regulation of cyclic adenosine 3,5-monophosphate (cAMP)-mediated signaling, and the Rap-1 signaling pathway were significantly enriched in Hainan Black and Alashan Cashmere goats. This result may be important for understanding each trait. Selection signature analysis revealed candidate single nucleotide polymorphisms (SNPs) and genes correlated with the traits of Hainan Black goats. Collectively, our results provide valuable insights into the genetic basis of specific traits correlated with the Hainan island climate, artificial selection in certain local goat breeds, and the importance of protecting breed resources.

Transcriptomics and Metabolomics Analysis of the Ovaries of High and Low Egg Production Chickens

Simple Summary

The ovarian tissues of different breeds of hens during egg production were investigated through transcriptomics and metabolomics to provide a more comprehensive understanding of the molecular mechanisms of the ovary during egg production. Four genes involved in egg production were predicted by the transcriptome, including P2RX1, INHBB, VIPR2, and FABP3, and several close metabolites associated with reproduction were identified in the metabolome, including 17α-hydroxyprogesterone, iloprost, spermidine and adenosine. Correlation analysis of specific differential genes and differential metabolites identified important gene-metabolite pairs VIPR2–Spermidine and P2RX1–Spermidine in the reproductive process.

Abstract

Egg production is a pivotal indicator for evaluating the fertility of poultry, and the ovary is an essential organ for egg production and plays an indispensable role in poultry production and reproduction. In order to investigate different aspects of egg production mechanisms in different poultry, in this study we performed a metabolomic analysis of the transcriptomic combination of the ovaries of two chicken breeds, the high-production Ninghai indigenous chickens and the low-production Wuliangshan black-boned chickens, to analyze the biosynthesis and potential key genes and metabolic pathways in the ovaries during egg production. We predicted four genes in the transcriptomic that are associated with egg production, namely P2RX1, INHBB, VIPR2, and FABP3, and identified three important pathways during egg production, “Calcium signaling pathway”, “Neuroactive ligand–receptor interaction” and “Cytokine–cytokine receptor interaction”, respectively. In the metabolomic 149 significantly differential metabolites were identified, 99 in the negative model and 50 in the positive model, of which 17α-hydroxyprogesterone, iloprost, spermidine, and adenosine are important metabolites involved in reproduction. By integrating transcriptomics and metabolomics, the correlation between specific differential genes and differential metabolites identified important gene-metabolite pairs “VIPR2-Spermidine” and “P2RX1-Spermidine” in egg production. In conclusion, these data provide a better understanding of the molecular differences between the ovaries of low- and high-production hens and provide a theoretical basis for further studies on the mechanics of poultry egg production.

Genetic Basis of Sexual Maturation Heterosis: Insights From Ovary lncRNA and mRNA Repertoire in Chicken

Sexual maturation is fundamental to the reproduction and production performance, heterosis of which has been widely used in animal crossbreeding. However, the underlying mechanism have long remained elusive, despite its profound biological and agricultural significance. In the current study, the reciprocal crossing between White Leghorns and Beijing You chickens were performed to measure the sexual maturation heterosis, and the ovary lncRNAs and mRNAs of purebreds and crossbreeds were profiled to illustrate molecular mechanism of heterosis. Heterosis larger than 20% was found for pubic space and oviduct length, whereas age at first egg showed negative heterosis in both crossbreeds. We identified 1170 known lncRNAs and 1994 putative lncRNAs in chicken ovary using a stringent pipeline. Gene expression pattern showed that nonadditivity was predominant, and the proportion of nonadditive lncRNAs and genes was similar between two crossbreeds, ranging from 44.24% to 49.15%. A total of 200 lncRNAs and 682 genes were shared by two crossbreeds, respectively. GO and KEGG analysis showed that the common genes were significantly enriched in the cell cycle, animal organ development, gonad development, ECM-receptor interaction, calcium signaling pathway and GnRH signaling pathway. Weighted gene co-expression network analysis (WGCNA) identified that 7 out of 20 co-expressed lncRNA-mRNA modules significantly correlated with oviduct length and pubic space. Interestingly, genes harbored in seven modules were also enriched in the similar biological process and pathways, in which nonadditive lncRNAs, such as MSTRG.17017.1 and MSTRG.6475.20, were strongly associated with nonadditive genes, such as CACNA1C and TGFB1 to affect gonad development and GnRH signaling pathway, respectively. Moreover, the results of real-time quantitative PCR (RT-qPCR) correlated well with the transcriptome data. Integrated with positive heterosis of serum GnRH and melatonin content detected in crossbreeds, we speculated that nonadditive genes involved in the GnRH signaling pathway elevated the gonad development, leading to the sexual maturation heterosis. We characterized a systematic landscape of ovary lncRNAs and mRNAs related to sexual maturation heterosis in chicken. The quantitative exploration of hybrid transcriptome changes lays foundation for genetic improvement of sexual maturation traits and provides insights into endocrine control of sexual maturation.

Transcriptome comparative analysis of ovarian follicles reveals the key genes and signaling pathways implicated in hen egg production

Background

Ovarian follicle development plays an important role in determination of poultry egg production. The follicles at the various developmental stages possess their own distinct molecular genetic characteristics and have different biological roles in chicken ovary development and function. In the each stage, several genes of follicle-specific expression and biological pathways are involved in the vary-sized follicular development and physiological events. Identification of the pivotal genes and signaling pathways that control the follicular development is helpful for understanding their exact regulatory functions and molecular mechanisms underlying egg-laying traits of laying hens.

Results

The comparative mRNA transcriptomic analysis of ovarian follicles at three key developmental stages including slow growing white follicles (GWF), small yellow follicles (SYF) of recruitment into the hierarchy, and differentiated large yellow follicles (LYF), was accomplished in the layers with lower and higher egg production. Totally, 137, 447, and 229 of up-regulated differentially expressed genes (DEGs), and 99, 97, and 157 of down-regulated DEGs in the GWF, SYF and LYF follicles, including VIPR1, VIPR2, ADRB2, and HSD17B1 were identified, respectively. Moreover, NDUFAB1 and GABRA1 genes, two most promising candidates potentially associated with egg-laying performance were screened out from the 13 co-expressed DEGs in the GWF, SYF and LYF samples. We further investigated the biological effects of NDUFAB1 and GABRA1 on ovarian follicular development and found that NDUFAB1 promotes follicle development by stimulating granulosa cell (GC) proliferation and decreasing cell apoptosis, increases the expression of CCND1 and BCL-2 but attenuates the expression of caspase-3, and facilitates steroidogenesis by enhancing the expression of STAR and CYP11A1. In contrast, GABRA1 inhibits GC proliferation and stimulates cell apoptosis, decreases the expression of CCND1, BCL-2, STAR, and CYP11A1 but elevates the expression of caspase-3. Furthermore, the three crucial signaling pathways such as PPAR signaling pathway, cAMP signaling pathway and neuroactive ligand-receptor interaction were significantly enriched, which may play essential roles in ovarian follicle growth, differentiation, follicle selection, and maturation.

Conclusions

The current study provided new molecular data for insight into the regulatory mechanism underlying ovarian follicle development associated with egg production in chicken.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-08213-w.

Polymorphism analysis and expression profile of the estrogen receptor 2 gene in Leizhou black duck

Our previous study on the ovarian transcriptomic analysis in Leizhou black duck revealed that the ESR2 gene was involved in hormone regulation in reproduction and the estrogen signaling pathway related to reproductive performance was enriched. This suggested that ESR2 may have a functional role in the reproductive performance of the Leizhou black duck. Thus, this study aimed at evaluating the polymorphism of the ESR2 gene and its association with egg-laying traits and the distribution pattern of ESR2 mRNA in laying and non-laying Leizhou black ducks. In this study, genomic DNA was extracted from blood samples of 101 Leizhou black ducks to identify single nucleotide polymorphisms (SNPs) of the ESR2 gene to elucidate molecular markers highly associated with egg-laying traits. Four each of laying and non-laying Leizhou black ducks were selected to collect different tissues to analyze the ESR2 gene expression. A total of 23 SNPs were identified and association analysis of the single SNP sites showed that SNPs g.56805646 T>C and exon 3-20G>A were significantly (P < 0.05) associated with egg weight. Ducks with CT and AG genotypes had significantly higher (P < 0.05) egg weights than their respective other genotypes. Haplotype association analysis of g.56805646 T>C and exon 3-20G>A showed that the haplotypes were significantly associated with egg weight. Higher egg weight was seen in individuals with H3H4 haplotypes. In the hypothalamus-pituitary-gonadal (HPG) axis, the results of qRT/PCR showed that ESR2 mRNA was significantly (P < 0.05) expressed in the ovaries of both duck groups than in the hypothalamus and pituitary. In the oviduct, ESR2 was significantly (P < 0.05) higher in the infundibulum and magnum of laying and non-laying ducks respectively. This study provides a molecular marker for selecting Leizhou black ducks for egg production. In addition, it offers theoretical knowledge for studying the related biological functions of the ESR2 gene at the cellular level.

Transcriptome analysis reveals key genes and pathways associated with egg production in Nandan-Yao domestic chicken

BACKGROUND

Egg production is a very important economic trait in chicken breeding, but its molecular mechanism is unclear until now. Nandan-Yao chicken (Gallus gallus domesticus) is a native breed in Guangxi province, China, which is famous for good meet quality, but with low egg production.

METHODS

To explore the molecular regulation related to egg production, high egg production (HEP) and low egg production (LEP) were divided according to the total egg number at 55 weeks, and the concentration of serum sex hormones was tested to evaluate the physiological function of ovary and uterus. RNA sequencing (RNA-Seq) was used to explore the transcriptome from the ovary and uterus of Nandan-Yao chicken.

RESULTS

The levels of serum sex hormone showed that concentrations of estradiol (E2), follicle-stimulating hormone (FSH), and luteotropic hormone (LH) were significantly higher in HEP than those in LEP (P < 0.01), while the concentration of testosterone (T) was significantly lower in HEP (P < 0.01). RNA-Seq analysis identified 901 and 2763 differentially expressed genes (DEGs) in ovary and uterus, respectively. Enrichment analysis showed that DEGs were significantly involved in the regulation of tight junction in the ovary (P < 0.05), while in uterus, DEGs were mainly enriched in the phagosome, ECM-receptor interaction, cell adhesion molecules (CAMs), focal adhesion, cardiac muscle contraction, cytokine-cytokine receptor interaction, and the regulation of MAPK signaling pathway (P < 0.05). Protein network interaction and function analyses revealed that FN1, FGF7, SOX2 identified from the ovary, and UQCRH, COX5A, FN1 from the uterus might be key candidate genes for egg production in Nandan-Yao chicken.

CONCLUSIONS

Our study provided key candidate genes and pathways involved in the egg-laying process of Nandan-Yao chicken and could help to further understand the molecular mechanisms of chicken reproduction.

Ovarian transcriptome profile from pre-laying period to broody period of Xupu goose

Xupu goose, a breed from Hunan province, produces high quality and quantity of meat and liver. However, its egg production rate is low, with poor reproductive traits but strong broody performance. These characteristics decrease the economic value of Xupu goose significantly. Here, RNA-seq was used to analyze the transcriptome changes of ovaries of Xupu goose at different stages to explore the molecular mechanism of reproduction from the pre-laying period to the broody period. A total of 258 genes were differentially expressed in the 3 stages. These genes are associated with inflammation, reproduction, mutual recognition and adhesion between cells, and cytoskeleton formation, and so on. In particular, we report, for the first time, the expression patterns of MRP126, serglycin, TXNIP, and FZD2 during the pre-laying, egg-laying, and broody periods of goose ovaries. Functional analysis by GO annotation revealed that GO terms were mainly involved in actin, cell signal transduction and regulation, and cellular components. Three pathways, including focal adhesion (gga04510), ECM-receptor interaction (gga04512), and N-Glycan biosynthesis (gga00510), were significantly enriched in the three groups. These findings provide a basis for further exploration of profiles of goose ovaries to improve egg production of Xupu goose.

Transcriptome Analysis of Ovarian Follicles Reveals Potential Pivotal Genes Associated With Increased and Decreased Rates of Chicken Egg Production

Egg production is an important economic trait in the commercial poultry industry. Ovarian follicle development plays a pivotal role in regulation of laying hen performance and reproductive physiology. However, the key genes and signaling pathways involved in the various-stages of laying hen follicular development remain poorly understood. In this study, transcriptomes of ovarian follicles at three developmental stages, the large white follicle (LWF), small yellow follicle (SYF), and large yellow follicle (LYF), were comparatively analyzed in hens with high (HR) and low (LR) egg-laying rates by RNA-sequencing. Eighteen cDNA libraries were constructed and a total of 236, 544, and 386 unigenes were significantly differentially expressed in the LWF, SYF, and LYF follicles of HR and LR hens, respectively. Among them, 47 co-transcribed differentially expressed genes (DEGs) in LWF and SYF, 68 co-expressed DEGs in SYF and LYF, and 54 co-expressed DEGs in LWF and LYF were mined. Thirteen co-expressed DEGs were found in LWF, SYF, and LYF follicles. Eighteen candidate genes, including P2RX1, CAB39L, BLK, CSMD3, GPR65, ADRB2, CSMD1, PLPP4, ATF3, PRLL, STMN3, RORB, PIK3R1, PERP1, ACSBG1, MRTO4, CDKN1A, and EDA2R were identified to be potentially related to egg production. Furthermore, Kyoto Encyclopedia of Genes and Genomes analysis indicated neuroactive ligand-receptor interaction, cell adhesion molecules, peroxisome proliferator-activated receptor pathway, and cAMP signaling pathway might elicit an important role in formation of egg-laying traits by influencing ovarian follicle development. This study represents the first transcriptome analysis of various-sized follicles between HR and LR hens. These results provide useful molecular evidence for elucidating the genetic mechanism underlying ovarian follicle development associated with egg production in chicken.

Comparative Transcriptome Profiling of Ovary Tissue between Black Muscovy Duck and White Muscovy Duck with High- and Low-Egg Production

The egg-laying rate is an important indicator for evaluating fertility of poultry. In order to better understand the laying mechanism of Muscovy ducks, gene expression profiles and pathways of ovarian tissues in high- and low-laying black (BH and BL) and white Muscovy ducks (WH and WL) during the peak production period were performed by using RNA-seq. The total number of reads produced for each ovarian sample ranged from 44,344,070 to 47,963,328. A total of 113, 619 and 87 differentially expressed genes (DEGs) were identified in BH-vs-WH, BL-vs-BH and BL-vs-WL, respectively. Among them, 54, 356 and 49 genes were up regulated and 59, 263 and 38 genes were down regulated. In addition, there were only 10 up-regulated genes in WL-vs-WH. In the comparison of DEGs in black and white Muscovy ducks, two co-expressed DEG genes were detected between BH-vs-WH and BL-vs-WL and seven DEGs were co-expressed between BL-vs-BH and WL-vs-WH. The RNA-Seq data were confirmed to be reliable by qPCR. Numerous DEGs known to be involved in ovarian development were identified, including TGFβ2, NGFR, CEBPD, CPEB2, POSTN, SMOC1, FGF18, EFNA5 and SDC4. Gene Ontology (GO) annotations indicated that DEGs related to ovarian development were mainly enriched in biological processes of "circadian sleep/wake cycle process," "negative regulation of transforming growth factor-β secretion," "positive regulation of calcium ion transport" in BH-vs-WH and "cell surface receptor signaling pathway," "Notch signaling pathway" and "calcium ion transport" in BL-vs-BH. Besides, "steroid biosynthetic process," "granulosa cell development" and "egg coat formation" were mainly enriched in BL-vs-WL and "reproduction," "MAPK cascade" and "mitotic cell cycle" were mainly enriched in WL-vs-WH. KEGG pathway analysis showed that the PI3K-Akt signaling pathway and ovarian steroidogenesis were the most enriched in Muscovy duck ovary transcriptome data. This work highlights potential genes and pathways that may affect ovarian development in Muscovy duck.

Ovarian transcriptomic analysis and follicular development of Leizhou black duck

This study investigated the factors that caused the differences in egg production during the development of ovarian follicles in Leizhou black ducks. Leizhou black ducks population was divided into 2 groups as high-yield group (HG) and low-yield group (LG). The number of eggs (NE), age at first egg (AFE), weight at first egg, and egg weight (EW) of both groups were recorded, and differences were analyzed using the t test. The logistic model was used to simulate the egg production curves to analyze the production rules. The ovarian follicles of both duck groups were collected to count the number of different grades sized follicles, weigh the ovaries, and observe follicular sections to analyze the developmental differences. Ovarian transcriptomic sequencing was performed to investigate differentially expressed genes and signal pathways in both duck groups. The results revealed a significant difference (P < 0.01) in the NE laid, AFE, and EW between both groups. Comparatively, HG had significantly more (P < 0.01) large yellow follicles (LYF) than LG. The density of medullary layer cells of the follicle section was greater in HG than LG ducks. Transcriptome sequencing revealed a total of 1,027 differentially expressed genes between the HG and LG ducks of which 495 genes were upregulated, and 532 genes were downregulated. Fifty genes were related to reproduction and reproductive processes. Kyoto Encyclopedia of Genes and Genomes–enriched signaling pathways revealed 274 signal pathways enriched in these differentially expressed genes of which the steroid biosynthesis pathway was significantly enriched. Analysis (Q < 0.05) showed that HSD3β → gonadotropin-releasing hormone (GnRH) and estrogen receptor (ESR) → LHβ/ERK1/2 were enriched in the steroid biosynthesis signal pathway. Follicle-stimulating hormone signal pathway mediated by HSD3β → GnRH and ESR → LHβ/ERK1/2 may be involved in ovarian follicle development to regulate LYF reserve process and affect its ovulation cycle, which in turn influence the egg production of Leizhou black ducks.

Transcriptome analysis reveals differentially expressed genes associated with high rates of egg production in chicken hypothalamic-pituitary-ovarian axis

The hypothalamic-pituitary-ovarian (HPO) axis regulates the breeding process cycle of laying hens. However, the key regulatory genes of the HPO axis and pathways that drive chicken egg laying performance remain elusive. A total of 856 Chinese Luhua chicken was raised and the highest two hundred and the lowest two hundred chicken egg production were considered as high egg production (HEP) and low egg production (LEP) according to the total egg number at 300 days of age, respectively. RNA-seq sequencing (RNA-Seq) was conducted to explore the chicken transcriptome from the hypothalamus, pituitary gland and ovary tissue of 6 Chinese Luhua chicken with 3 high and low-rate egg production. In total, 76.09 Gb RNA-seq sequences were generated from 15 libraries with an average of 5.07 Gb for each library. Further analysis showed that 414, 356 and 10 differentially expressed genes (DEGs) were identified in pituitary gland, ovary and hypothalamus between HEP and LEP chickens, respectively. In pituitary gland, DEGs were involve in regulation of cellular glucose homeostasis, Ras protein signal transduction, negative regulation of hormone secretion. In Ovary DEGs were mainly involved in embryonic organ development, regulation of canonical Wnt signaling, response to peptide hormone. Our study identified DEGs that regulate mTOR signaling pathway, Jak-STAT signaling pathway, Tryptophan metabolism and PI3K-Akt signaling pathways at HPO-axis in laying hens. These important data contribute to improve our understanding of reproductive biology of chicken and isolating effective molecular markers that can be used for genetic selection in Chinese domestic Luhua chicken.

Chromosome Missegregation in Single Human Oocytes Is Related to the Age and Gene Expression Profile

The growing trend for women to postpone childbearing has resulted in a dramatic increase in the incidence of aneuploid pregnancies. Despite the importance to human reproductive health, the events precipitating female age-related meiotic errors are poorly understood. To gain new insight into the molecular basis of age-related chromosome missegregation in human oocytes, we combined the transcriptome profiles of twenty single oocytes (derived from females divided into two groups according to age <35 and ≥35 years) with their chromosome status obtained by array comparative genomic hybridization (aCGH). Furthermore, we compared the transcription profile of the single oocyte with the surrounding cumulus cells (CCs). RNA-seq data showed differences in gene expression between young and old oocytes. Dysregulated genes play a role in important biological processes such as gene transcription regulation, cytoskeleton organization, pathways related to RNA maturation and translation. The comparison of the transcription profile of the oocyte and the corresponding CCs highlighted the differential expression of genes belonging to the G protein-coupled receptor superfamily. Finally, we detected the loss of a X chromosome in two oocytes derived from women belonging to the ≥35 years age group. These aneuploidies may be caused by the detriment of REEP4, an endoplasmic reticulum protein, in women aged ≥35 years. Here we gained new insight into the complex regulatory circuit between the oocyte and the surrounding CCs and uncovered a new putative molecular basis of age-related chromosome missegregation in human oocytes.

Comparative Genomics Reveals the Genetic Mechanisms of Musk Secretion and Adaptive Immunity in Chinese Forest Musk Deer

The Chinese forest musk deer (Moschus berezovskii; FMD) is an artiodactyl mammal and is both economically valuable and highly endangered. To investigate the genetic mechanisms of musk secretion and adaptive immunity in FMD, we compared its genome to nine other artiodactyl genomes. Comparative genomics demonstrated that eight positively selected genes (PSGs) in FMD were annotated in three KEGG pathways that were related to metabolic and synthetic activity of musk, similar to previous transcriptome studies. Functional enrichment analysis indicated that many PSGs were involved in the regulation of immune system processes, implying important reorganization of the immune system in FMD. FMD-specific missense mutations were found in two PSGs (MHC class II antigen DRA and ADA) that were classified as deleterious by PolyPhen-2, possibly contributing to immune adaptation to infectious diseases. Functional assessment showed that the FMD-specific mutation enhanced the ADA activity, which was likely to strengthen the immune defense against pathogenic invasion. Single nucleotide polymorphism-based inference showed the recent demographic trajectory for FMD. Our data and findings provide valuable genomic resources not only for studying the genetic mechanisms of musk secretion and adaptive immunity, but also for facilitating more effective management of the captive breeding programs for this endangered species.

Differential gene expression analysis of human cumulus cells

Objective

This study was performed to explore the possibility that each oocyte and its surrounding cumulus cells might have different genetic expression patterns that could affect human reproduction.

Methods

Differential gene expression analysis was performed for 10 clusters of cumulus cells obtained from 10 cumulus-oocyte complexes from 10 patients. Same procedures related to oocyte maturation, microinjection, and microarray analyses were performed for each group of cumulus cells. Two differential gene expression analyses were performed: one for the outcome of clinical pregnancy and one for the outcome of live birth.

Results

Significant genes resulting from these analyses were selected and the top 20 affected pathways in each group were analyzed. Circadian entrainment is determined to be the most affected pathway for clinical pregnancy, and proteoglycans in cancer pathway is the most affected pathway for live birth. Circadian entrainment is also amongst the 12 pathways that are found to be in top 20 affected pathways for both outcomes, and has both lowest p-value and highest number of times found count.

Conclusion

Although further confirmatory studies are necessary, findings of this study suggest that these pathways, especially circadian entrainment in cumulus cells, may be essential for embryo development and pregnancy.

Hypothalamic transcriptome analysis reveals the neuroendocrine mechanisms in controlling broodiness of Muscovy duck (Cairina moschata)

Broodiness, one of the maternal behaviors and instincts for natural breeding in birds, is an interesting topic in reproductive biology. Broodiness in poultry is characterized by persistent nesting, usually associated with cessation of egg laying. The study of avian broodiness is essential for bird conservation breeding and commercial poultry industry. In this study, we examined the hypothalamus transcriptome of Muscovy duck in three reproductive stages, including egg-laying anaphase (LA), brooding prophase (BP) and brooding metaphase (BM). Differences in gene expression during the transition from egg-laying to broodiness were examined, and 155, 379, 292 differently expressed genes (DEGs) were obtained by pairwise comparisons of LA-vs-BP, LA-vs-BM and BP-vs-BM, respectively (fold change≥1.5, P < 0.05). Gene Ontology Term (GO) enrichment analysis suggested a possible role of oxidative stress in the hypothalamus might invoke reproductive costs that potentially change genes expression. KEGG analysis revealed glutamatergic synapse, dopaminergic synapse, serotonergic synapse and GABAergic synapse pathway were significantly enriched, and regulator genes were identified. Eight gene expression patterns were illustrated by trend analysis and further clustered into three clusters. Additional six hub genes were identified through combining trend analysis and protein-protein interaction (PPI) analysis. Our results suggested that the cyclical mechanisms of reproductive function conversion include effects of oxidative stress, biosynthesis of neurotransmitters or their receptors, and interactions between glucocorticoids and thyroid hormones and regulatory genes. These candidate genes and biological pathways may be used as targets for artificial manipulation and marker-assisted breeding in the reproductive behavior.

Comparative proteomic analysis of hypothalamus tissue from Huoyan geese between pre-laying period and laying period using an iTRAQ-based approach

The hypothalamus plays a central role in controlling poultry endocrine and reproductive activities. So far there is limited information focused on the proteome profiles of the hypothalamus from geese during different stages of the egg-laying cycle. In order to identify proteins regulating the egg-laying process of Huoyan geese, we investigated the proteome profiles of the hypothalamus from Huoyan geese during the laying period and pre-laying period by applying an isobaric tags for relative and absolute quantitation (iTRAQ)-based proteomic technology. A total number of 3,337 were identified and quantified, of which 18 were significantly up-regulated and 16 were significantly down-regulated. These differentially expressed proteins were subjected to bioinformatics analyses based on the Gene Ontology annotation and Kyoto Encyclopedia of Genes and Genomes pathway. Some of these were revealed to be involved in hormone and neurotransmitter secretion, exocytosis, calcium ion transport and synaptic transmission. Subsequently, excitatory amino acid transporter 2, complexin-1 and inositol 1,4,5-trisphosphate receptor, type 3 were confirmed at the messenger RNA level using quantitative real-time RT-PCR. Then, the abundance change of these proteins was verified further using Western blotting analysis. These data may aid in elucidating the molecular mechanism of higher laying performance in Huoyan geese.

Transcriptome analysis revealed the possible regulatory pathways initiating female geese broodiness within the hypothalamic-pituitary-gonadal axis

Geese have the strongest tendency toward broodiness among all poultry. The mechanisms initiating broodiness within the goose hypothalamic-pituitary-gonadal axis (HPGA) are still unclear. Here, we reported the transcriptome differences between laying and initial nesting within the HPGA tissues of geese. We constructed a unigene database based on HPGA tissues and identified 128,148 unigenes, 100% of which have been annotated. By using Digital Gene Expression (DGE) sequencing, we screened 19, 110, 289, and 211 differentially expressed genes (DEGs) in the hypothalamus, pituitary gland, stroma ovarii, and follicles, respectively, between laying and nesting geese. Expression changes of hypocretin (HCRT) and pro-opiomelanocortin (POMC) in the hypothalamus of nesting geese may cause appetite reduction, which is possibly the first step and a prerequisite to initiate broodiness. In addition to prolactin (PRL), follicle-stimulating hormone (FSH) and luteinizing hormone (LH), genes including oxytocin-neurophysin (OXT), chordin-like protein 1 (CHRDL1) and growth hormone (GH), expressed in the pituitary gland, are new candidate molecules that may be involved in broodiness in geese. Heme oxygenase 1 (HMOX1) in the pituitary gland, the proto-oncogene c-Fos (FOS), heat shock protein 90-alpha (HSP90AA), and cyclin-dependent kinase 1 (CDK1) in the ovary that may consolidate and transduce signals regulating the HPGA during broodiness in geese.

Comparative proteomic analysis of pituitary glands from Huoyan geese between pre-laying and laying periods using an iTRAQ-based approach

In this study, we performed a comprehensive evaluation of the proteomic profile of the pituitary gland of the Huoyan goose during the laying period compared to the pre-laying period using an iTRAQ-based approach. Protein samples were prepared from pituitary gland tissues of nine pre-laying period and nine laying period geese. Then the protein samples from three randomly selected geese within each period were pooled in equal amounts to generate one biological sample pool. We identified 684 differentially expressed proteins, including 418 up-regulated and 266 down-regulated proteins. GO annotation and KEGG pathway analyses of these proteins were conducted. Some of these proteins were found to be associated with hormone and neurotransmitter secretion and transport, neuropeptide signalling and GnRH signalling pathways, among others. Subsequently, the modification of the abundance of three proteins (prolactin, chromogranin-A and ITPR3) was verified using western blotting. Our results will provide a new source for mining genes and gene products related to the egg-laying performance of Huoyan geese, and may provide important information for the conservation and utilization of local goose breeds.

Identification and analysis of differentially expressed long non-coding RNAs between multiparous and uniparous goat (Capra hircus) ovaries

Long non-coding RNAs (lncRNAs) play important roles in almost all biological processes. However, there is little information on the effects of lncRNAs on ovulation and lambing rates. In the present study, we used high-throughput RNA sequencing to identify differentially expressed lncRNAs between the ovaries of multiparous (Mul) and uniparous (Uni) Anhui White goats. Among the 107,255,422 clean reads, 183,754 lncRNAs were significantly differentially expressed between the Uni and Mul. Among them, 455 lncRNAs were co-expressed between the two samples, whereas, 157,523 lncRNAs were uniquely expressed in the Uni, and 25,776 uniquely lncRNAs were expressed in the Mul. Through Cis role analysis, 24 lncRNAs were predicted to overlap with cis-regulatory elements, which involved in Progesterone-mediated oocyte maturation, Steroid biosynthesis, Oocyte meiosis, and gonadotropin-releasing hormone (GnRH) signaling pathway. These 4 pathways were related to ovulation, and the KEGG pathway analysis on target genes of the differentially expressed lncRNAs confirmed this results. In addition, 10 lncRNAs harbored precursors of 40 miRNAs, such as TCONS_00320849 related to a mature miRNA sequence, miR-365a, which was reported to be related to proliferation, were annotated in the precursor analysis of miRNAs. The present expand the understanding of lncRNA biology and contribute to the annotation of the goat genome. The study will provide a resource for lncRNA studies of ovulation and lambing.

Ovarian transcriptomic analysis of Shan Ma ducks at peak and late stages of egg production

Objective

To assess the differences in ovarian transcriptomes in Shan Ma ducks between their peak and late stages of egg production, and to obtain new transcriptomic data of these egg-producing ducks.

Methods

The Illumina HiSeq 2000 system was used for high throughput sequencing of ovarian transcriptomes from Shan Ma ducks at their peak or late stages of egg production.

Results

Greater than 93% of the sequencing data had a base quality score (Q score) that was not less than 20 (Q20). From ducks at their peak stage of egg production, 42,782,676 reads were obtained, with 4,307,499,083 bp sequenced. From ducks at their late stage of egg production, 45,316,166 reads were obtained, with 4,562,063,363 bp sequenced. A comparison of the two datasets identified 2,002 differentially expressed genes, with 790 upregulated and 1,212 downregulated. Further analysis showed that 1,645 of the 2,002 differentially expressed genes were annotated in the non-redundant (NR) database, with 646 upregulated and 999 downregulated. Among the differentially expressed genes with annotations in the NR database, 696 genes were functionally annotated in the clusters of orthologous groups of proteins database, involving 25 functional categories. One thousand two hundred four of the differentially expressed genes with annotations in the NR database were functionally annotated in the gene ontology (GO) database, and could be divided into three domains and 56 categories. The three domains were cellular component, molecular function, and biological process. Among the genes identified in the GO database, 451 are involved in development and reproduction. Analysis of the differentially expressed genes with annotations in the NR database against the Kyoto encyclopedia of genes and genomes database revealed that 446 of the genes could be assigned to 175 metabolic pathways, of which the peroxisome proliferator-activated receptor signaling pathway, insulin signaling pathway, fructose and mannose metabolic pathways, gonadotropin releasing hormone signaling pathway and transforming growth factor beta signaling pathway were significantly enriched.

Conclusion

The differences in ovarian transcriptomes in Shan Ma ducks between their peak and late stages of egg production were elucidated, which greatly enriched the ovarian transcriptomic information of egg-producing ducks.

Expressed microRNA associated with high rate of egg production in chicken ovarian follicles

MicroRNA (miRNA) is a highly conserved class of small noncoding RNA about 19-24 nucleotides in length that function in a specific manner to post-transcriptionally regulate gene expression in organisms. Tissue miRNA expression studies have discovered a myriad of functions for miRNAs in various aspects, but a role for miRNAs in chicken ovarian tissue at 300 days of age has not hitherto been reported. In this study, we performed the first miRNA analysis of ovarian tissues in chickens with low and high rates of egg production using high-throughput sequencing. By comparing low rate of egg production chickens with high rate of egg production chickens, 17 significantly differentially expressed miRNAs were found (P < 0.05), including 11 known and six novel miRNAs. We found that all 11 known miRNAs were involved mainly in pathways of reproduction regulation, such as steroid hormone biosynthesis and dopaminergic synapse. Additionally, expression profiling of six randomly selected differentially regulated miRNAs were validated by quantitative real-time polymerase chain reaction (RT-qPCR). Some miRNAs, such as gga-miR-34b, gga-miR-34c and gga-miR-216b, were reported to regulate processes such as proliferation, cell cycle, apoptosis and metastasis and were expressed differentially in ovaries of chickens with high rates of egg production, suggesting that these miRNAs have an important role in ovary development and reproductive management of chicken. Furthermore, we uncovered that a significantly up-regulated miRNA-gga-miR-200a-3p-is ubiquitous in reproduction-regulation-related pathways. This miRNA may play a special central role in the reproductive management of chicken, and needs to be further studied for confirmation.

Molecular cloning of the SMAD4 gene and its mRNA expression analysis in ovarian follicles of the Yangzhou goose (Anser cygnoides)

Mothers against decapentaplegic homolog 4 (SMAD4) is an important protein in animal reproduction. It plays pivotal roles in cellular pathways, including apoptosis. The expression profile of the SMAD4 gene in goose ovarian follicles has not been reported. In this study, the SMAD4 coding sequence was cloned from the Yangzhou goose. A phylogenetic analysis was performed and mRNA expression was examined in various tissues using quantitative real-time PCR. An alternative splice form of SMAD4, SMAD4-b having 1656 bp, was identified. SMAD4-a mRNA was widely expressed in various healthy tissues, whereas SMAD4-b was very weakly expressed. SMAD4 mRNA in the ovary and oviduct was significantly higher than that in the pituitary and hypothalamus. SMAD4 mRNA expression analysis in hierarchical follicles showed that the level of SMAD4 mRNA was higher in large white follicles and post-ovulatory follicles than in the other follicles. The results indicate that SMAD4 might be involved in the recruitment of hierarchical follicles.

High-throughput sequencing reveals differential expression of miRNAs in prehierarchal follicles of laying and brooding geese

Broodiness is the primary factor influencing egg production in geese, in which several genes and miRNAs participate. Detailed spatiotemporal profiles of miRNAs encompassing follicle development levels, however, are lacking. In this study, we collected preovulatory follicles (classified as small white follicles, large white follicles, and small yellow follicles) from brooding and laying geese and aimed to analyze microRNA (miRNA or miR) during folliculogenesis. High-throughput sequencing and bioinformatics analysis were used to identify the miRNAs involved in follicle development. The let7 family, miR-10 family, and miR-143 family were abundant in these libraries, and they have been suggested to play a housekeeping role during folliculogenesis. Joint comparisons revealed 23 upregulated and 21 downregulated miRNAs (in at least two comparisons of follicles during brooding and laying, P < 0.1) in the laying stage. Unlike reproduction pathways reported for ovaries, GO and KEGG analysis suggested pathways for cell apoptosis and proliferation, such as the regulation of actin cytoskeleton, endocytosis, axon guidance, pathways in cancer, tight junctions, focal adhesion, the MAPK signaling pathway, cytokine-cytokine receptor interactions, and the Wnt signaling pathway in folliculogenesis. This study revealed the miRNAs that were directly involved in follicular atresia, and our results added to the understanding of the functional involvement of miRNAs during specific stages of follicle development.

Transcriptome Analysis and Identification of Differentially Expressed Transcripts of Immune-Related Genes in Spleen of Gosling and Adult Goose

The goose (Anser cygnoides), having high nutritional value, high-quality feathers and high economic benefit, is an economically important poultry species. However, the molecular mechanisms underlying the higher susceptibility to pathogens in goslings than in adult geese remains poorly understood. In this study, the histological sections of spleen tissue from a two-week-old gosling and an adult goose, respectively, were subjected to comparative analysis. The spleen of gosling was mainly composed of mesenchyma, accompanied by scattered lymphocytes, whereas the spleen parenchyma was well developed in the adult goose. To investigate goose immune-related genes, we performed deep transcriptome and gene expression analyses of the spleen samples using paired-end sequencing technology (Illumina). In total, 50,390 unigenes were assembled using Trinity software and TGICL software. Moreover, these assembled unigenes were annotated with gene descriptions and gene ontology (GO) analysis was performed. Through Kyoto encyclopedia of genes and genomes (KEGG) analysis, we investigated 558 important immune-relevant unigenes and 23 predicted cytokines. In addition, 22 immune-related genes with differential expression between gosling and adult goose were identified, among which the three genes showing largest differences in expression were immunoglobulin alpha heavy chain (IgH), mannan-binding lectin serine protease 1 isoform X1 (MASP1) and C-X-C chemokine receptor type 4 (CXCR4). Finally, of these 22 differentially expressed immune-related genes, seven genes, including tumor necrosis factor receptor superfamily member 13B (TNFRSF13B), C-C motif chemokine 4-like (CCL4), CXCR4, interleukin 2 receptor alpha (IL2RA), MHC class I heavy chain (MHCIα), transporter of antigen processing 2 (TAP2) IgH, were confirmed by quantitative real-time PCR (qRT-PCR). The expression levels of all the candidate unigenes were up-regulated in adult geese other than that of TNFRSF13B. The comparative analysis of the spleen transcriptomes of gosling and adult goose may promote better understanding of immune molecular development in goose.

De novo transcriptomic analysis of peripheral blood lymphocytes from the Chinese goose: gene discovery and immune system pathway description

Background

The Chinese goose is one of the most economically important poultry birds and is a natural reservoir for many avian viruses. However, the nature and regulation of the innate and adaptive immune systems of this waterfowl species are not completely understood due to limited information on the goose genome. Recently, transcriptome sequencing technology was applied in the genomic studies focused on novel gene discovery. Thus, this study described the transcriptome of the goose peripheral blood lymphocytes to identify immunity relevant genes.

Principal Findings

De novo transcriptome assembly of the goose peripheral blood lymphocytes was sequenced by Illumina-Solexa technology. In total, 211,198 unigenes were assembled from the 69.36 million cleaned reads. The average length, N50 size and the maximum length of the assembled unigenes were 687 bp, 1,298 bp and 18,992 bp, respectively. A total of 36,854 unigenes showed similarity by BLAST search against the NCBI non-redundant (Nr) protein database. For functional classification, 163,161 unigenes were comprised of three Gene Ontology (Go) categories and 67 subcategories. A total of 15,334 unigenes were annotated into 25 eukaryotic orthologous groups (KOGs) categories. Kyoto Encyclopedia of Genes and Genomes (KEGG) database annotated 39,585 unigenes into six biological functional groups and 308 pathways. Among the 2,757 unigenes that participated in the 15 immune system KEGG pathways, 125 of the most important immune relevant genes were summarized and analyzed by STRING analysis to identify gene interactions and relationships. Moreover, 10 genes were confirmed by PCR and analyzed. Of these 125 unigenes, 109 unigenes, approximately 87%, were not previously identified in the goose.

Conclusion

This de novo transcriptome analysis could provide important Chinese goose sequence information and highlights the value of new gene discovery, pathways investigation and immune system gene identification, and comparison with other avian species as useful tools to understand the goose immune system.