(C) PLOS One
This story was originally published by PLOS One and is unaltered.
. . . . . . . . . .



Dissection of key factors correlating with H5N1 avian influenza virus driven inflammatory lung injury of chicken identified by single-cell analysis [1]

['Manman Dai', 'National', 'Regional Joint Engineering Laboratory For Medicament Of Zoonosis Prevention', 'Control', 'Guangdong Provincial Key Laboratory Of Zoonosis Prevention', 'College Of Veterinary Medicine', 'South China Agricultural University', 'Guangzhou', 'Sufang Zhu', 'Zhihao An']

Date: 2023-11

Abstract Chicken lung is an important target organ of avian influenza virus (AIV) infection, and different pathogenic virus strains lead to opposite prognosis. Using a single-cell RNA sequencing (scRNA-seq) assay, we systematically and sequentially analyzed the transcriptome of 16 cell types (19 clusters) in the lung tissue of chickens infected with H5N1 highly pathogenic avian influenza virus (HPAIV) and H9N2 low pathogenic avian influenza virus (LPAIV), respectively. Notably, we developed a valuable catalog of marker genes for these cell types. Compared to H9N2 AIV infection, H5N1 AIV infection induced extensive virus replication and the immune reaction across most cell types simultaneously. More importantly, we propose that infiltrating inflammatory macrophages (clusters 0, 1, and 14) with massive viral replication, pro-inflammatory cytokines (IFN-β, IL1β, IL6 and IL8), and emerging interaction of various cell populations through CCL4, CCL19 and CXCL13, potentially contributed to the H5N1 AIV driven inflammatory lung injury. Our data revealed complex but distinct immune response landscapes in the lung tissue of chickens after H5N1 and H9N2 AIV infection, and deciphered the potential mechanisms underlying AIV-driven inflammatory reactions in chicken. Furthermore, this article provides a rich database for the molecular basis of different cell-type responses to AIV infection.

Author summary The low-pathogenicity avian influenza virus (LPAIV), H9N2, and highly pathogenic avian influenza virus (HPAIV), H5N1 are the main epidemic subtypes, resulting in great economic losses on poultry and potentially threat to human. Thus, in-depth exploration of AIV pathogenesis in chickens is necessary for developing efficient control methods. In our study, using a single-cell RNA sequencing analysis of 16 cell types in the chicken lung, we revealed complicated and distinct immune response landscapes after H5N1 and H9N2 AIV infection, and identified key factors contributing to H5N1 AIV driven inflammatory lung injury. Therefore, our study potentially provides new targets and direction for AIV control.

Citation: Dai M, Zhu S, An Z, You B, Li Z, Yao Y, et al. (2023) Dissection of key factors correlating with H5N1 avian influenza virus driven inflammatory lung injury of chicken identified by single-cell analysis. PLoS Pathog 19(10): e1011685. https://doi.org/10.1371/journal.ppat.1011685 Editor: Kanta Subbarao, The Peter Doherty Institute and Melbourne University, AUSTRALIA Received: December 4, 2022; Accepted: September 13, 2023; Published: October 11, 2023 Copyright: © 2023 Dai et al. This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited. Data Availability: Sequencing data have been deposited in Genome Sequence Archive (GSA) of National Genomics Data Center under the accession number CRA006707 (https://ngdc.cncb.ac.cn/search/?dbId=gsa&q=CRA006707&page=1). Funding: This work was supported by grants from the National Natural Science Foundation of China (31830097 to ML); National Natural Science Foundation of China (32211530072 to MD); National Key R&D Program of China (2022YFD1801005 to MD); the National Natural Science Foundation of Guangdong Province (2022A1515012480 to MD); Science and Technology Plan Project of Guangzhou City (202201010512 to MD); Young Elite Scientists Sponsorship Program by CAST (2020QNRC001 to MD); Biotechnology and Biological Sciences Research Council (BBSRC) Newton Fund Joint Centre Awards on “UK-China Centre of Excellence for Research on Avian Diseases” (BBS/OS/NW/000007 to VN). Special fund for scientific innovation strategy-construction of high-level Academy of Agriculture Science-Distinguished Scholar (R2020PY-JC001 and 202110TD to ML). The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript. Competing interests: The authors have declared that no competing interests exist.

Introduction The low-pathogenic avian influenza virus (LPAIV), H9N2, and highly pathogenic avian influenza virus (HPAIV), H5N1 are the main epidemic subtypes in the ongoing virus circulation among Chinese poultry causing major economic losses in spite of the long-term vaccination programs [1–3]. Notably, the ongoing 2021–2022 wave of avian influenza H5N1 in Asia, Africa and Europe is unprecedented in its rapid spread and extremely high frequency of outbreaks in poultry [4]. More seriously, H5N1 virus can be transmitted from chickens to humans. Also, H9N2 virus can serve as the gene donor for H5N1, H7N9 and H10N8 viruses infecting humans [1–3,5–7]. Therefore, successful control of H5N1 HPAIV and H9N2 LPAIV in chickens is vital for the eradication of diseases and preventing infections in humans. Hence, it warrants deeper understanding of the determinant factors of avian influenza virus (AIV) infection and pathogenesis in chickens for developing efficient control methods. AIV infections of the chicken mainly occur via the respiratory route, and the lung is the important target organ. In lung tissue, AIV infection induces both antiviral and inflammatory factors which play crucial roles in host protection and immunopathogenesis [8,9]. Infection of chickens with the H9N2 LPAIV usually results in mild clinical signs whereas H5N1 HPAIV induces death within 36–48 hours, which is related to heightened inflammatory responses in the latter [10]. However, the specific immune cell types and inflammatory factors contributing to the immune injury of chicken lung are poorly understood. Besides, efficient replication of HPAIV compared to LPAIV has been correlated with tissue damage in the lung [11,12]. But the extent and nature of HPAIV and LPAIV infection in different cell types of chicken lung has not been elucidated. In this study, we want to explore the spectrum of cell infection types, and the immune reaction profiles of chicken lung at the peak of H5N1 HPAIV and H9N2 LPAIV infections. Single-cell RNA sequencing (scRNA-seq) has been widely applied to identify the involvement of different cell types and to investigate the immune responses under conditions of virus infection [13–15]. It is also a powerful tool for defining viral target cells via analyzing the viral mRNAs and host signature genes in a single cell [13,16–18]. In particular, scRNA-seq can precisely be used to examine the patterns of cytokine release in each immune cell and intercommunication with the ligand and receptor interaction map at the single cell level [13,19,20]. To our knowledge, this is the first study on the application of scRNA-seq technology to examine the responses of the major immune cell types in the lung tissue of chickens infected with H5N1 and H9N2 AIV. In order to explore more comprehensive and refined immune cell responses to H5N1 and H9N2 AIV infection, we sorted the major immune cell types, including MHC Class II antigen presenting cells and CD3 positive T cells from chicken lung mononuclear cells for performing scRNA-seq analysis. In this study, we analyzed 16,642 immune cells in the lung tissue isolated from chickens after AIV infection with different pathogenic strains at various time points. Our results revealed complex but distinct immune response landscapes in the lung tissue of chickens after H5N1 and H9N2 AIV infection. Our study has developed a valuable catalog of marker genes for identifying 16 cell types in the lung tissue of chicken via scRNA-seq. We also provided (for the first time) the key immune cell types and pro-inflammatory factors that contribute to H5N1 AIV-driven inflammatory injury in chicken lung tissue.

Materials and methods Ethics statement All animal experiments were carried out in ABSL-3 facilities in compliance with approved protocol (CNAS BL0011) by the biosafety committee of South China Agriculture University (Guangzhou, China). All animal procedures were performed according to the regulations and guidelines established by this committee and international standards for animal welfare. Virus and experimental animal infection Low pathogenic avian influenza H9N2 subtype HN strain (A/Chicken/Hunan/HN/2015), and high pathogenic avian influenza H5N1 subtype DK383 strain (A/Duck/Guangdong/383/2008) were isolated and identified by our research team [21,22]. Nine four-week-old, specific pathogen–free (SPF) White Leghorn chickens (Guangdong Da Hua Nong Animal Health Products Co., Ltd., Guangdong, China) were randomly assigned to three groups, namely, a H9N2-infected group, a H5N1-infected group and control group, each with 3 chickens. The H9N2-infected group was intranasally (i.n.) inoculated with HN strain (107 EID50 in 0.2 mL). The H5N1-infected group was intranasally inoculated with DK383 strain (106 EID50 in 0.2 mL). The control group was inoculated with 0.2 mL of phosphate buffered saline (PBS). Three chickens were dissected at 1-day post-inoculation (DPI) in the H5N1-infected group, 3 DPI in the H9N2-infected group, and 0 DPI in the control group, respectively. The virus titer in selected organs and virus shedding were detected as previously described [2], and the data are shown in S1 File. Immunohistochemistry (IHC) and Preparation of chicken lung mononuclear cell suspension Infected or mock-infected chickens were humanely sacrificed. Their lungs were then perfused with RPMI-1640 containing 2% fetal bovine serum (FBS) (Invitrogen, Waltham, MA, USA) via the right ventricle. Half of the entire pool of the lung was fixed in 4% paraformaldehyde for 24 hours and embedded in paraffin. Tissue slides were deparaffinized with xylene, hydrated in a graded series of alcohol solutions to distilled water, and blocked for endogenous peroxidase in 3% hydrogen peroxide. Sections were then stained with mouse anti-Influenza A virus NP primary antibody (GeneTex, Alton, CA, USA), followed by staining with anti-mouse-HRP secondary antibody (Zhongshan Goldenbridge, Beijing, China) and visualizing via optical microscopy (Olympus, Tokyo, Japan). Then, 0.1g from the other half of the lung tissue was kept for virus titer detection, and the remaining half of the lung pool was dissected and dissociated into single-cell suspensions, using the Lung Dissociation Kit (Miltenyi Biotec, Bergisch Gladbach, Germany) in combination with a gentleMACS dissociator (Miltenyi Biotec) and enzymatic dissociation [23]. Following enzymatic incubation, cells were forced through a 70-μm mesh cell strainer. After that, lung mononuclear cells were enriched from lung suspension with the tissue mononuclear cell kit (Haoyang, Tianjin, China), as described previously [2,24]. Cell viability and counting was performed using Trypan Blue and hemocytometer (Sigma-Aldrich). Sample and library preparation for 10x scRNAseq 3×107 lung cell pooled suspension was prepared from three chickens, with 1×107 lung cell suspension for each treatment group. MHC Class II and CD3 positive cells were sorted out from each pulmonary cell pool after incubation with the FITC-conjugated MHC Class II and APC-conjugated CD3 antibodies (SouthernBiotech, Birmingham, AL) using Fluorescence-activated cell sorting machine (FACS Aria II, Becton Dickinson, New Jersey, USA). The gating strategy for each population is shown in S1 Fig. Equal number of MHC Class II and CD3 positive cells in each pulmonary cell pool were mixed together, followed by passing through a 40 μm cell strainer (Biosharp, China), and the cell viability was above 80%. Subsequently, the cell density was adjusted to 1 x 106 cells/mL. High quality single cell suspension was subjected to encapsulation using a 10x Genomics v.3 kit (10x Genomics, USA). The library preparation and RNA-sequencing were completed by Gene Denovo (Guangzhou, China) as described previously [15]. An average of 32956 reads per cell in the H9N2 group, mean reads of 35911 per cell in the H5N1 group, and mean reads of 28419 per cell in the control group were obtained, respectively. FACS analysis of Macrophages in lung mononuclear cell suspension The percentage of macrophages from lung single cell suspension was detected via flow cytometry (CytoFLEX, Beckman Coulter, Brea, CA, USA) with the FITC-conjugated mouse anti-chicken MHC Class II and PE-conjugated mouse anti-chicken KUL01 monoclonal antibodies (Southern Biotech, Birmingham, USA) [25]. Data were analyzed by the software of FlowJo V10 (TreestarInc, Ashland, OR, USA). Absolute number of macrophages was calculated by multiplying the percentage of macrophages with the total cell number of single cell suspension isolated from lung. Macrophage sorting and generation of macrophage reference transcriptional profile via SMART-Seq2 based scRNA-seq Fluorescence-activated cell sorting machine (FACS Aria II, Becton Dickinson, New Jersey, USA) was used to sort a single cell into each well of a 96-well PCR plate containing 2.5μL of 10× Lysis Buffer (Vazyme# N711). For the isolation of macrophage, each pulmonary cell pool was stained with FITC-conjugated mouse anti-chicken MHC Class II and PE-conjugated mouse anti-chicken KUL01 monoclonal antibodies (Southern Biotech, Birmingham, USA). Macrophages were sorted following the manufacturer’s procedures. Herein, each sorted cell population was analyzed in four replicates, and 100 single cells were sorted in each replicate for subsequent SMART-Seq2 analysis. Four empty wells served as controls in each 96-well plate. Immediately after sorting, each plate was spun down to ensure immersion of cells into the lysis solution, snap-frozen on dry ice, and stored at −80°C until processing. Library construction and sequencing were completed by Gene Denovo (Guangzhou, China) as described previously [26]. Pearson’s correlation analysis was used to investigate the relationships between the SMART-Seq2 data of macrophages and lung clusters in scRNA-seq based on levels of gene expression. Quantitative Reverse Transcription PCR (qRT-PCR) Total RNA was extracted from lung cell suspensions or sorted macrophages for qRT-PCR. qRT-PCR was performed on an ABI7500 Real-Time PCR system (Applied Biosystems, Waltham, MA) using Universal SYBR R Green Supermix Kit reagents (Biorad, CA, USA) according to the manufacturer’s specifications. Primers used for qRT-PCR are listed in S2 File. Data analysis were performed using the 2-ΔΔCt method [27]. 10x scRNAseq sequencing data analysis As described previously [15], the raw scRNA-seq data were aligned, filtered, and normalized using Cell Ranger (10x Genomics) software (Cell Ranger 3.1.0), and the cDNA reads were mapped to the chicken genome of GRCg6a [28]. Only reads that were confidently mapped to the transcriptome were used for Unique Molecular Identifier (UMI) counting. Cells with unusually high numbers of UMIs (≥8000) or mitochondrial gene percentage (≥20%) were filtered out. Cells with <500 or >4000 gene counts were also excluded. Seurat (v4.0.4) is a popular R package that was used to implement the graph-based clustering approaches [29,30]. T-distributed Stochastic Neighbor Embedding (t-SNE) [31] or Uniform Manifold Approximation and Projection (UMAP) [32] in Seurat were used to visualize and explore these datasets. Other data analyses including standardization, difference of gene expression, and marker gene screening were also achieved by Seurat. Cluster marker analysis and cell type annotation Four strategies were used for cell-type annotation. Firstly, we identified cell types based on the top 5 expressed genes in each cluster and classical marker genes of chicken immune cells, e.g: CD4 T cells (CD3D, IL7R and CD4), B cells (BCL11A, Bu-1 and BLB2), CD8 T cells (CD3D, CD8A and GNLY), Cytotoxic T cells (CD3D, CD8A, GNLY, GranzymeA and IFNG), and Dendritic Cells (DCs, BLB2, CD80, CD86 and XCR1) [33,34]. Secondly, the identity of some cell clusters was manually annotated based on the specific expression of the marker genes that had orthologous genes as commonly known markers published in the CellMarker database (http://iocc.hrbmu.edu.cn/CellMarker/) [35], e.g: Th17 cells (CD3D, RORA, CCR6 and CCL20), Macrophages (IL1B, VCAN, and BLB2), and Smooth muscle cells (MYH11 and TAGLN). Next, identification of most of the remaining cell clusters was based on the specific expression of the marker genes reported in the literature, e.g: Epithelial cells (BMX and EHF) [36], type II alveolar epithelial cells (SFTPA1, SFTPA2 and SFTPC) [37–39], Fibroblasts (COL1A1, COL1A2 and COL6A1) [14], Vascular endothelial cells (ADGRL4, VWF, PODXL and SELECTIN) [40–42], M2 macrophages (BLB2, RNASE6 and VSIG4) [43,44], and Regulatory T cells (Tregs, CD3D, CD25, IL7R and TNFRSF1B) [45,46]. Last, we combined the marker genes in the CellMarker database with reported distinctive genes in the literature for annotating the remaining clusters, e.g: Macrophage like cells (BLB2, VCAN, IL1B, HSPH1 and DNAJA4) [47], and Th2 cells (CD3D, MAF, CCR4, DRD4, and KK34) [15]. The distribution of the characteristic gene expression in each cluster was then demonstrated using heat maps and bubble diagrams. The detailed marker gene information for 18 clusters and cell type annotation is shown in S3 and S4 Files. Viral Gene Analysis in Cell Clusters of H9N2-infected group, H5N1-infected group and control group To determine the state of H9N2 AIV or H5N1 AIV infection in each cell type, the expression of viral genes and viral load were analyzed in the different cell types in the H9N2-infected group, H5N1-infected group and control group. Viral gene expression was analyzed using Cell Ranger, based on the viral sequence of H9N2 subtype HN strain (A/Chicken/Hunan/HN/2015) or H5N1 subtype DK383 strain (A/Duck/Guangdong/383/2008) used in this study. The ‘viral genome load’ of a cell in scRNA-seq analysis was based on the number of UMIs that mapped to the AIV eight segmented mRNAs and expressed as a percentage of total UMI content of a given cell. To reduce the false positive rate of infected cells, only the cells with highly expressed AIV genes (at least one transcript per gene per cell) were defined as highly infected cells. Therefore, the cells in each cell type susceptible to AIV infection could be divided into highly infected cells (I, total UMI counts of viral transcripts ≥ 8), potential or lowly infected cells (P, total UMI counts of viral transcripts ≥ 1), and undetected cells (N, UMI counts of viral transcripts = 0) [13]. Differentially Expressed Genes (DEGs) Analysis in each cell type between H5N1 or H9N2-infected and uninfected groups To explore the response of each cell type to AIV, we further analyzed the DEGs between the H5N1 or H9N2-infected group and the control group using Seurat software. A hurdle model in MAST (Model-based Analysis of Single-cell Transcriptomics) [48] was used to find DEGs for a group in each cell type. DEGs between the H5N1or H9N2-infected group and control groups were identified by the following criteria: 1) |log2FC|≥0.36; 2) p value_adj ≤ 0.05; and 3) percentage of cells where the gene was detected in a specific cluster ≥ 10%. Identified DEGs were subsequently subjected to GO enrichment analysis. Ligand-receptor intercellular communication network analysis We used cellphone DB [49] to infer ligand-receptor pairs in chicken cells via orthologous gene analysis. A ligand/receptor with non-zero expression in more than 10% of cells in a particular cell population was deemed an “expressed” ligand/receptor. Then igraph in R was used to draw interaction maps as previously described [50]. In order to show the lost and/or new interaction after H5N1 or H9N2 AIV infection treatment, we calculated all the interactions for control and the treatment (H5N1 or H9N2 AIV) and then generated an interaction map from all interactions of control and the treatment. Then, we colored the edges as cyan (lost after treatment), purple (induced after treatment) and black (unaffected by treatment but may be differentially expressed). Statistical analysis Data were analyzed using GraphPad Prism 8.0 software (GraphPad Software Inc., San Diego, CA, USA). The results were presented as mean ± SEM. The paired t-test and one-way ANOVA were used for statistical comparison. Data were considered significant at * P < 0.05, ** P < 0.01, *** P < 0.001.

Discussion Multiple key questions related to the response of the immune cells to AIV infection in chicken lung have yet to be answered. In particular, the identity of the key immune cell types, antiviral and inflammatory factors contributing to the host protection and immunopathogenesis have not been systematically elucidated. Advances in scRNA-seq technology have allowed us to explore the atlas of the immune cells response to AIV infection in chicken lungs. In this study, we performed 10x scRNA-seq on the sorted immune cells including MHC Class II (antigen presenting cells) and CD3 (T cells) positive cells for characterizing the responses of major immune cell populations in the lung tissue isolated from chickens after H5N1 HPAIV(1DPI), H9N2 LPAIV (3DPI) and PBS treatment. Here, 19 distinct clusters in chicken lungs were identified following the analysis of scRNA-seq data. In addition to the 11 specialized immune cell types (14 clusters) that were identified and annotated, we also identified 5 clusters of non-immune cell types which may be accidentally sorted out during the experiment. In total, 16,642 immune cells were profiled and analyzed including various T cell subsets (CD8 T cell, Cytotoxic T cell, CD4 T cell, Th2, Th17 and Treg) and antigen presentation cells (B cell, DC, Macrophage, Macrophage like cell, and M2 macrophage) (S3 and S5 Files). Moreover, we verified that the SMART-Seq2 data of sorted macrophages was highly correlated to the marker-based annotation of the macrophage population (Clusters 0, 1 and 14) (Fig 5B), which further indicated the reliability of the annotation of scRNA-seq clusters in this study. Thus for the first time (to our knowledge), we provided a valuable catalog of marker genes for identifying 16 cell types in lung tissue of chickens via scRNA-seq (Fig 1 and S3 File). In future studies, more samples and cell numbers will be necessary for identifying the complete cell atlas of the chicken lung. Compared to LPAIV infection, the efficient virus replication of HPAIV can be correlated with tissue damage in the lung [11,12]. But the extent and nature of HPAIV and LPAIV infection in different cell types of chicken lung have not been elucidated. In this study, we found that H9N2 LPAIV seemed to be limited in infecting different cell types. Instead, H5N1 HPAIV widely infected all cell types identified by single-cell sequence (Fig 3). The difference in the tropism of cell types between LPAI and HPAI seems to align with the clinical outcome of chickens after LPAIV and HPAIV infection. In mice studies, it has been reported that influenza A Virus could infect all cell types including T cells [14,53], which is in agreement with our findings on chickens. Besides the widespread and efficient virus replication, H5N1 AIV infection induced a stronger antiviral immune response than H9N2 AIV infection in lung (Fig 2C and 2F, and S7 File). Moreover, the antiviral genes were induced and appeared in nearly all H5N1 infected cell types (Fig 2C), which indicated that the host response conferred the first line of defense in all key cell types. HE staining results revealed that H5N1 AIV, as against H9N2 AIV infection, induced pulmonary inflammatory damage with inflammatory cell infiltration (S3 Fig), which prompted us to investigate the key immune cells and genes contributing to the H5N1 AIV-driven pneumonia. By sequentially analyzing the transcriptome of various cell types, we identified that macrophage, M2 macrophage and macrophage-like populations generated abundant pro-inflammatory factors (Figs 4 and 5A). Meanwhile, the high levels of H5N1 virus infection were discovered in these cell populations simultaneously (Fig 3E and S8 File), which reminded that massive viral replication potentially induced an excessive immune response. Besides, we found that macrophages (Clusters 0, 1 and 14) were hugely increased in the lung from H9N2 or H5N1 AIV infected chickens based on the data from scRNA-seq (Fig 1F) and FACS analysis (Fig 1G and 1H) detected using KUL01 and MHC Class II antibodies [25,54]. Compared to other cell types, macrophages demonstrated the most DEGs, especially pro-inflammatory factors responding to H9N2 and H5N1 AIV infection (Figs 2 and 4). More importantly, we found that pro-inflammatory cytokines including IFN-β, IL1β, IL6 and IL8 [51,52] and chemokines containing CCL1 and CCL19, were significantly up-regulated in the macrophages of H5N1 group as against the H9N2 group as determined by qRT-PCR (Fig 5C and 5D). Therefore, we would consider that moderate expression levels of monocytes/macrophages and inflammatory factors favor lung damage repair after H9N2 AIV infection. However, uncontrollable viral replication and cytokine release would induce serious lung immune pathology after H5N1 AIV infection. Collectively, we especially identified infiltrating macrophages with massive viral replication contributing to the excessive cytokines released and immune injury after H5N1 infection. In mice studies, it has been reported that the communication networks between various cell types play an important role in influenza A virus (IAV)-induced cytokine storm and pneumonia [13]. Thus, we plotted the intercommunication predictions for cell types expressing the receptor or ligand of the pro-inflammatory factors. Intriguingly, we have identified the interaction categories based on the CCL4, CCL19, and CXCL13 and their receptors after H5N1 AIV challenge, as against H9N2 AIV infection (Fig 5E), which may contribute to the inflammatory lung injury. The mechanism by which the intercommunication across various cell clusters induces the release of inflammatory factors has yet to be explored. In summary, through scRNA-seq analysis, we demonstrated the key factors associated with the pathogenesis of AIV infection in chickens. Importantly, by sequentially analyzing the transcriptome of various cell types, we discovered that infiltrating macrophages with massive viral replication and emerging interaction of various cell populations through CCL4, CCL19 and CXCL13, may contribute to the H5N1 AIV driven inflammatory lung injury. Paralleled with the previous findings from the literature, our results and validation studies confirm the fidelity of our analyses and interpretations. Our data also provide extensive resources for future studies to address the function of identified cell types and in response to pathogenic infection in chickens.

Acknowledgments We are grateful to the South China Agricultural University’s high-level talent launch program and the "Fuji Peiyou" program of the College of Veterinary Medicine, South China Agricultural University. We are extremely appreciative of the help given by Gene Denovo Corp. during bioinformatics analysis.

[END]
---
[1] Url: https://journals.plos.org/plospathogens/article?id=10.1371/journal.ppat.1011685

Published and (C) by PLOS One
Content appears here under this condition or license: Creative Commons - Attribution BY 4.0.

via Magical.Fish Gopher News Feeds:
gopher://magical.fish/1/feeds/news/plosone/

You are viewing proxied material from magical.fish. The copyright of proxied material belongs to its original authors. Any comments or complaints in relation to proxied material should be directed to the original authors of the content concerned. Please see the disclaimer for more details.