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

Comparative transcriptomics analysis identifies crucial genes and pathways during goose spleen development

As the largest peripheral lymphoid organ in poultry, the spleen plays an essential role in regulating the body's immune capacity. However, compared with chickens and ducks, information about the age- and breed-related changes in the goose spleen remains scarce. In this study, we systematically analyzed and compared the age-dependent changes in the morphological, histological, and transcriptomic characteristics between Landes goose (LG; Anser anser) and Sichuan White goose (SWG; Anser cygnoides). The results showed a gradual increase in the splenic weights for both LG and SWG until week 10, while their splenic organ indexes reached the peak at week 6. Meanwhile, the splenic histological indexes of both goose breeds continuously increased with age, reaching the highest levels at week 30. The red pulp (RP) area was significantly higher in SWG than in LG at week 0, while the splenic corpuscle (AL) diameter was significantly larger in LG than in SWG at week 30. At the transcriptomic level, a total of 1710 and 1266 differentially expressed genes (DEGs) between week 0 and week 30 were identified in spleens of LG and SWG, respectively. Meanwhile, a total of 911 and 808 DEGs in spleens between LG and SWG were identified at weeks 0 and 30, respectively. Both GO and KEGG enrichment analysis showed that the age-related DEGs of LG or SWG were dominantly enriched in the Cell cycle, TGF-beta signaling, and Wnt signaling pathways, while most of the breed-related DEGs were enriched in the Neuroactive ligand-receptor interaction, Cytokine-cytokine receptor interaction, ECM-receptor interaction, and metabolic pathways. Furthermore, through construction of protein-protein interaction networks using significant DEGs, it was inferred that three hub genes including BUB1, BUB1B, and TTK could play crucial roles in regulating age-dependent goose spleen development while GRIA2, GRIA4, and RYR2 could be crucial for the breed-specific goose spleen development. These data provide novel insights into the splenic developmental differences between Chinese and European domestic geese, and the identified crucial pathways and genes are helpful for a better understanding of the mechanisms regulating goose immune functions.

Characterization of circular RNA expression profiles in the age-related thymic involution of Magang goose

The thymus is a vital immune organ, but its function gradually declines with age. Circular RNAs (circRNAs) are related to the development of tissues and organs. In this study, bioinformatics analysis showed that 1329, 755, and 417 circRNAs were differentially expressed between the comparison groups of 6-month age (M6) and 20-embryo age (E20), 3-day post-hatch (P3), and 3-month age (M3) Magang geese, respectively. Among them, 167 circRNAs were differentially co-expressed between thymic development (E20, P3, and M3) and involution (M6). Functional analysis showed significant enrichment of phosphorylation and positive regulation of GTPase activity. Furthermore, pathway analysis has shown that glycerolipid metabolism and the Wnt signaling pathway are critical pathways in the thymic involution process. Finally, we constructed the competitive endogenous RNA (ceRNA) network. The results of this study suggest that circRNAs may be involved in the age-related thymic involution of the Magang goose.

Network Meta-Analysis of Chicken Microarray Data Following Avian Influenza Challenge—A Comparison of Highly and Lowly Pathogenic Strains

The current bioinformatics study was undertaken to analyze the transcriptome of chicken (Gallus gallus) after influenza A virus challenge. A meta-analysis was carried out to explore the host expression response after challenge with lowly pathogenic avian influenza (LPAI) (H1N1, H2N3, H5N2, H5N3 and H9N2) and with highly pathogenic avian influenza (HPAI) H5N1 strains. To do so, ten microarray datasets obtained from the Gene Expression Omnibus (GEO) database were normalized and meta-analyzed for the LPAI and HPAI host response individually. Different undirected networks were constructed and their metrics determined e.g., degree centrality, closeness centrality, harmonic centrality, subgraph centrality and eigenvector centrality. The results showed that, based on criteria of centrality, the CMTR1, EPSTI1, RNF213, HERC4L, IFIT5 and LY96 genes were the most significant during HPAI challenge, with PARD6G, HMG20A, PEX14, RNF151 and TLK1L having the lowest values. However, for LPAI challenge, ZDHHC9, IMMP2L, COX7C, RBM18, DCTN3, and NDUFB1 genes had the largest values for aforementioned criteria, with GTF3C5, DROSHA, ATRX, RFWD2, MED23 and SEC23B genes having the lowest values. The results of this study can be used as a basis for future development of treatments/preventions of the effects of avian influenza in chicken.

Safety and Transcriptome Analysis of Live Attenuated Brucella Vaccine Strain S2 on Non-pregnant Cynomolgus Monkeys Without Abortive Effect on Pregnant Cynomolgus Monkeys

Brucellosis, caused by Brucella spp., is an important zoonotic disease leading to enormous economic losses in livestock, posing a great threat to public health worldwide. The live attenuated Brucella suis (B. suis) strain S2, a safe and effective vaccine, is widely used in animals in China. However, S2 vaccination in animals may raise debates and concerns in terms of safety to primates, particularly humans. In this study, we used cynomolgus monkey as an animal model to evaluate the safety of the S2 vaccine strain on primates. In addition, we performed transcriptome analysis to determine gene expression profiling on cynomolgus monkeys immunized with the S2 vaccine. Our results suggested that the S2 vaccine was safe for cynomolgus monkeys. The transcriptome analysis identified 663 differentially expressed genes (DEGs), of which 348 were significantly upregulated and 315 were remarkably downregulated. The Gene Ontology (GO) classification and the Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis indicated that these DEGs were involved in various biological processes (BPs), including the chemokine signaling pathway, actin cytoskeleton regulation, the defense response, immune system processing, and the type-I interferon signaling pathway. The molecular functions of the DEGs were mainly comprised of 2'-5'-oligoadenylate synthetase activity, double-stranded RNA binding, and actin-binding. Moreover, the cellular components of these DEGs included integrin complex, myosin II complex, and blood microparticle. Our findings alleviate the concerns over the safety of the S2 vaccine on primates and provide a genetic basis for the response from a mammalian host following vaccination with the S2 vaccine.

Transcriptome analysis of Corynebacterium pseudotuberculosis-infected spleen of dairy goats

Caseous lymphadenitis is a chronic disease of goats caused by Corynebacterium pseudotuberculosis (C.pseudotuberculosis) which causes great harm to the dairy goats industry. In order to obtain detailed information about the pathogenesis and host immune response in C.pseudotuberculosis-infected goats, in this study, the gene expression difference of spleen tissue after infection with C.pseudotuberculosis was analyzed by high-throughput sequencing. Transcripts obtained over 412 700 462 clean reads after reassembly were 21 343 genes detected, of which 14 720 were known genes and 7623 new genes were predicted. There were 448 up-regulated and 519 down-regulated differentially expressed genes (DEGs). Gene Ontology (GO) analysis indicated that all of the DEGs were annotated into biological process, cellular component and molecular function. Most of these unigenes are annotated in cellular processes, the cell and binding. KEGG analysis of the DEGs showed that a total of 8733 DEGs unigenes were annotated into 459 pathways classified into 6 main categories. Most of these annotated unigenes were related to immune system response to the infectious diseases pathways. In addition, 14 DEGs were verified by quantitative real-time PCR. As the first, in vivo, RNAseq analysis of dairy goats and C.pseudotuberculosis infection, this study provides knowledge about the transcriptomics of spleen in C.pseudotuberculosis-infected goats, from which a complex molecular pathways and immune response mechanism are involved in C.pseudotuberculosis infection.

Midgut transcriptome analysis of Clostera anachoreta treated with lethal and sublethal Cry1Ac protoxin

Clostera anachoreta is one of the important Lepidoptera insect pests in forestry, especially in poplars woods in China, Europe, Japan, and India, and so forth, and also the target insect of Cry1Ac toxin and Bt plants. Six genes, HSC70, GNB2L/RACK1, PNLIP, BI1-like, arylphorin type 2, and PKM were found in this study, and they might be associated with the response to the Cry1Ac toxin, found by analyzing the transcriptome data. And the PI3K-Akt pathway was highly enriched in differentially expressed unigenes and linked to several crucial pathways, including the B-cell receptor signaling pathway, toll-like receptor pathway, and mitogen-activated protein kinase signaling pathway. They might be involved in the recovery stage of the damaged midgut during the response to sublethal doses of Cry1Ac toxin. This is the first study conducted to specifically investigate C. anachoreta response to Cry toxin stress using large-scale sequencing technologies, and the results highlighted some important genes and pathways that could be involved in Btcry1Ac resistance development or could serve as targets for biologically based control mechanisms of this insect pest.

Blood transcriptome resources of chinstrap (Pygoscelis antarcticus) and gentoo (Pygoscelis papua) penguins from the South Shetland Islands, Antarctica

The chinstrap (Pygoscelis antarcticus) and gentoo (P. papua) penguins are distributed throughout Antarctica and the sub-Antarctic islands. In this study, high-quality de novo assemblies of blood transcriptomes from these penguins were generated using the Illumina MiSeq platform. A total of 22.2 and 21.8 raw reads were obtained from chinstrap and gentoo penguins, respectively. These reads were assembled using the Oases assembly platform and resulted in 26,036 and 21,854 contigs with N50 values of 929 and 933 base pairs, respectively. Functional gene annotations through pathway analyses of the Gene Ontology, EuKaryotic Orthologous Groups, and Kyoto Encyclopedia of Genes and Genomes (KEGG) databases were performed for each blood transcriptome, resulting in a similar compositional order between the two transcriptomes. Ortholog comparisons with previously published transcriptomes from the Adélie (P. adeliae) and emperor (Aptenodytes forsteri) penguins revealed that a high proportion of the four penguins’ transcriptomes had significant sequence homology. Because blood and tissues of penguins have been used to monitor pollution in Antarctica, immune parameters in blood could be important indicators for understanding the health status of penguins and other Antarctic animals. In the blood transcriptomes, KEGG analyses detected many essential genes involved in the major innate immunity pathways, which are key metabolic pathways for maintaining homeostasis against exogenous infections or toxins. Blood transcriptome studies such as this may be useful for checking the immune and health status of penguins without sacrifice.

Transcriptome profiling reveals the immune response of goose T cells under selenium stimuli

The goose is an economically important poultry species and a principal natural host of avian viruses. This study aimed to determine the effects of selenium on the immune response of geese. Under selenium stimulation, gene expression profiling was investigated using transcriptome sequencing. The selenoproteins were promoted by selenium stimulation, while the heat shock proteins, interleukin and interferons were mainly down-regulated. After comparison, 2228 differentially expressed genes were primarily involved in immune and environmental response, and infectious disease and genetic information processing related pathways were identified. Specifically, the enzymes of the lysosomes which acted as a safeguard in preventing pathogens were mostly up-regulated and six randomly selected differentially expressed genes were validated by quantitative polymerase chain reaction. In addition, the most proportional increased transcription factor family basic helix-loop-helix (bHLH) located in the 5' flank of selenoprotein P-like protein for selenium metabolism was identified by response to the selenium stimulation in this study. These analyses show that selenium can promote immune function by activating selenoproteins, transcript factors and lysosome pathway related genes, while weakening cytokine content genes in geese.

C1q: A fresh look upon an old molecule

Originally discovered as part of C1, the initiation component of the classical complement pathway, it is now appreciated that C1q regulates a variety of cellular processes independent of complement activation. C1q is a complex glycoprotein assembled from 18 polypeptide chains, with a C-terminal globular head region that mediates recognition of diverse molecular structures, and an N-terminal collagen-like tail that mediates immune effector mechanisms. C1q mediates a variety of immunoregulatory functions considered important in the prevention of autoimmunity such as the enhancement of phagocytosis, regulation of cytokine production by antigen presenting cells, and subsequent alteration in T-lymphocyte maturation. Furthermore, recent advances indicate additional roles for C1q in diverse physiologic and pathologic processes including pregnancy, tissue repair, and cancer. Finally, C1q is emerging as a critical component of neuronal network refinement and homeostatic regulation within the central nervous system. This review summarizes the classical functions of C1q and reviews novel discoveries within the field.

OAZ1 knockdown enhances viability and inhibits ER and LHR transcriptions of granulosa cells in geese

An increasing number of studies suggest that ornithine decarboxylase antizyme 1 (OAZ1), which is regarded as a tumor suppressor gene, regulates follicular development, ovulation, and steroidogenesis. The granulosa cells in the ovary play a critical role in these ovarian functions. However, the action of OAZ1 mediating physiological functions of granulosa cells is obscure. OAZ1 knockdown in granulosa cells of geese was carried out in the current study. The effect of OAZ1 knockdown on polyamine metabolism, cell proliferation, apoptosis, and hormone receptor transcription of primary granulosa cells in geese was measured. The viability of granulosa cells transfected with the shRNA OAZ1 at 48 h was significantly higher than the control (p<0.05). The level of putrescine and spermidine in granulosa cells down-regulating OAZ1 was 7.04- and 2.11- fold higher compared with the control, respectively (p<0.05). The CCND1, SMAD1, and BCL-2 mRNA expression levels in granulosa cells down-regulating OAZ1 were each significantly higher than the control, respectively (p<0.05), whereas the PCNA and CASPASE 3 expression levels were significantly lower than the control (p<0.05). The estradiol concentration, ER and LHR mRNA expression levels were significantly lower in granulosa cells down-regulating OAZ1 compared with the control (p<0.05). Taken together, our results indicated that OAZ1 knockdown elevated the putrescine and spermidine contents and enhanced granulosa cell viability and inhibited ER and LHR transcriptions of granulosa cells in geese.

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.

The first whole genome and transcriptome of the cinereous vulture reveals adaptation in the gastric and immune defense systems and possible convergent evolution between the Old and New World vultures

BACKGROUND

The cinereous vulture, Aegypius monachus, is the largest bird of prey and plays a key role in the ecosystem by removing carcasses, thus preventing the spread of diseases. Its feeding habits force it to cope with constant exposure to pathogens, making this species an interesting target for discovering functionally selected genetic variants. Furthermore, the presence of two independently evolved vulture groups, Old World and New World vultures, provides a natural experiment in which to investigate convergent evolution due to obligate scavenging.

RESULTS

We sequenced the genome of a cinereous vulture, and mapped it to the bald eagle reference genome, a close relative with a divergence time of 18 million years. By comparing the cinereous vulture to other avian genomes, we find positively selected genetic variations in this species associated with respiration, likely linked to their ability of immune defense responses and gastric acid secretion, consistent with their ability to digest carcasses. Comparisons between the Old World and New World vulture groups suggest convergent gene evolution. We assemble the cinereous vulture blood transcriptome from a second individual, and annotate genes. Finally, we infer the demographic history of the cinereous vulture which shows marked fluctuations in effective population size during the late Pleistocene.

CONCLUSIONS

We present the first genome and transcriptome analyses of the cinereous vulture compared to other avian genomes and transcriptomes, revealing genetic signatures of dietary and environmental adaptations accompanied by possible convergent evolution between the Old World and New World vultures.

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.