Molecular evidence of the avocado defense response to Fusarium kuroshium infection: a deep transcriptome analysis using RNA-Seq

Fusarium kuroshium is a novel member of the Ambrosia Fusarium Clade (AFC) that has been recognized as one of the symbionts of the invasive Kuroshio shot hole borer, an Asian ambrosia beetle. This complex is considered the causal agent of Fusarium dieback, a disease that has severely threatened natural forests, landscape trees, and avocado orchards in the last 8 years. Despite the interest in this species, the molecular responses of both the host and F. kuroshium during the infection process and disease establishment remain unknown. In this work, we established an in vitro pathosystem using Hass avocado stems inoculated with F. kuroshium to investigate differential gene expression at 1, 4, 7 and 14 days post-inoculation. RNA-seq technology allowed us to obtain data from both the plant and the fungus, and the sequences obtained from both organisms were analyzed independently. The pathosystem established was able to mimic Fusarium dieback symptoms, such as carbohydrate exudation, necrosis, and vascular tissue discoloration. The results provide interesting evidence regarding the genes that may play roles in the avocado defense response to Fusarium dieback disease. The avocado data set comprised a coding sequence collection of 51,379 UniGenes, from which 2,403 (4.67%) were identified as differentially expressed. The global expression analysis showed that F. kuroshium responsive UniGenes can be clustered into six groups according to their expression profiles. The biologically relevant functional categories that were identified included photosynthesis as well as responses to stress, hormones, abscisic acid, and water deprivation. Additionally, processes such as oxidation-reduction, organization and biogenesis of the cell wall and polysaccharide metabolism were detected. Moreover, we identified orthologues of nucleotide-binding leucine-rich receptors, and their possible action mode was analyzed. In F. kuroshium, we identified 57 differentially expressed genes. Interestingly, the alcohol metabolic process biological category had the highest number of upregulated genes, and the enzyme group in this category may play an important role in the mechanisms of secondary metabolite detoxification. Hydrolytic enzymes, such as endoglucanases and a pectate lyase, were also identified, as well as some proteases. In conclusion, our research was conducted mainly to explain how the vascular tissue of a recognized host of the ambrosia complex responds during F. kuroshium infection since Fusarium dieback is an ambrosia beetle-vectored disease and many variables facilitate its establishment.

Download Full-text

HIV-1 Infection Transcriptomics: Meta-Analysis of CD4+ T Cells Gene Expression Profiles

Viruses ◽

10.3390/v13020244 ◽

2021 ◽

Vol 13 (2) ◽

pp. 244 ◽

Cited By ~ 1

Author(s):

Antonio Victor Campos Coelho ◽

Rossella Gratton ◽

João Paulo Britto de Melo ◽

José Leandro Andrade-Santos ◽

Rafael Lima Guimarães ◽

...

Keyword(s):

Gene Expression ◽

T Cells ◽

Differentially Expressed Genes ◽

Cd4 T Cells ◽

Expression Profiles ◽

Meta Analysis ◽

Differentially Expressed ◽

Rna Seq ◽

Infected Cells ◽

Hiv 1

HIV-1 infection elicits a complex dynamic of the expression various host genes. High throughput sequencing added an expressive amount of information regarding HIV-1 infections and pathogenesis. RNA sequencing (RNA-Seq) is currently the tool of choice to investigate gene expression in a several range of experimental setting. This study aims at performing a meta-analysis of RNA-Seq expression profiles in samples of HIV-1 infected CD4+ T cells compared to uninfected cells to assess consistently differentially expressed genes in the context of HIV-1 infection. We selected two studies (22 samples: 15 experimentally infected and 7 mock-infected). We found 208 differentially expressed genes in infected cells when compared to uninfected/mock-infected cells. This result had moderate overlap when compared to previous studies of HIV-1 infection transcriptomics, but we identified 64 genes already known to interact with HIV-1 according to the HIV-1 Human Interaction Database. A gene ontology (GO) analysis revealed enrichment of several pathways involved in immune response, cell adhesion, cell migration, inflammation, apoptosis, Wnt, Notch and ERK/MAPK signaling.

Download Full-text

Comparative transcriptomics analysis revealing flower trichome development during flower development in two Lonicera Japonica Thunb. cultivars using RNA-seq

10.21203/rs.3.rs-19879/v1 ◽

2020 ◽

Author(s):

Jianjun Li ◽

Chenglin Ye ◽

Cuifang Chang

Keyword(s):

Signaling Pathways ◽

Plant Hormones ◽

Flower Development ◽

Expression Profiles ◽

Differentially Expressed ◽

Molecular Networks ◽

Lonicera Japonica ◽

Rna Seq ◽

Trichome Development ◽

A Genome

Abstract Background: Trichomes comprise specialized multicellular structures that have the capacity to synthesize and secrete secondary metabolites and protect plants from biotic and abiotic stresses. However, little is known about the trichome formation mechanism during flower development in Lonicera Japonica Thunb.Results: Here, we present a genome-wide comparative transcriptome analysis between two L. japonica cultivars, toward the identification of biological processes and functional gene activities that occur during flowering stage trichome development. In this study, the density and average lengths of flower trichomes were at their highest during three green periods. Using the Illumina RNA-Seq method, we obtained 134,304 unigenes, 33,733 of which were differentially expressed. In an analysis of 40 differentially expressed unigenes (DEGs) involved in trichome development, 29 of these were transcription factors. The DEGs analysis of plant hormone signal transduction indicated that plant growth and development may be independent of GA and CTK signaling pathways, and plant stress may be independent of JA and ET signaling pathways. We successfully isolated key genes involved in the floral biosynthesis of odors, tastes, colors, and plant hormones, and proposed biosynthetic pathways for sesquiterpenoid, triterpenoid, monoterpenoid, flavonoid, and plant hormones. Furthermore, 82 DEGs were assigned to cell cycles and 2,616 were predicted as plant resistance genes (PRGs).Conclusions: This study provides a comprehensive characterization of the expression profiles of flower development during the seven developmental stages of L. japonica, thereby offering valuable insights into the molecular networks that underly flower development in L. japonica.

Download Full-text

Integrative transcriptome data mining for identification of core lncRNAs in breast cancer

PeerJ ◽

10.7717/peerj.7821 ◽

2019 ◽

Vol 7 ◽

pp. e7821 ◽

Cited By ~ 7

Author(s):

Xiaoming Zhang ◽

Jing Zhuang ◽

Lijuan Liu ◽

Zhengguo He ◽

Cun Liu ◽

...

Keyword(s):

Breast Cancer ◽

Gene Expression ◽

Prognostic Value ◽

Early Stage ◽

Expression Profiles ◽

Diagnostic Value ◽

The Cancer Genome Atlas ◽

Functional Enrichment ◽

Differentially Expressed ◽

Data Set

Background Cumulative evidence suggests that long non-coding RNAs (lncRNAs) play an important role in tumorigenesis. This study aims to identify lncRNAs that can serve as new biomarkers for breast cancer diagnosis or screening. Methods First, the linear fitting method was used to identify differentially expressed genes from the breast cancer RNA expression profiles in The Cancer Genome Atlas (TCGA). Next, the diagnostic value of all differentially expressed lncRNAs was evaluated using a receiver operating characteristic (ROC) curve. Then, the top ten lncRNAs with the highest diagnostic value were selected as core genes for clinical characteristics and prognosis analysis. Furthermore, core lncRNA-mRNA co-expression networks based on weighted gene co-expression network analysis (WGCNA) were constructed, and functional enrichment analysis was performed using the Database for Annotation, Visualization and Integrated Discovery (DAVID). The differential expression level and diagnostic value of core lncRNAs were further evaluated by using independent data set from Gene Expression Omnibus (GEO). Finally, the expression status and prognostic value of core lncRNAs in various tumors were analyzed based on Gene Expression Profiling Interactive Analysis (GEPIA). Results Seven core lncRNAs (LINC00478, PGM5-AS1, AL035610.1, MIR143HG, RP11-175K6.1, AC005550.4, and MIR497HG) have good single-factor diagnostic value for breast cancer. AC093850.2 has a prognostic value for breast cancer. AC005550.4 and MIR497HG can better distinguish breast cancer patients in early-stage from the advanced-stage. Low expression of MAGI2-AS3, LINC00478, AL035610.1, MIR143HG, and MIR145 may be associated with lymph node metastasis in breast cancer. Conclusion Our study provides candidate biomarkers for the diagnosis and prognosis of breast cancer, as well as a bioinformatics basis for the further elucidation of the molecular pathological mechanism of breast cancer.

Download Full-text

Comparative analysis of gene expression profiles in differentiated subcutaneous adipocytes between Jiaxing Black and Large White Pigs

10.21203/rs.3.rs-19837/v1 ◽

2020 ◽

Author(s):

Dawei Zhang ◽

Wenjing Wu ◽

Xin Huang ◽

Ke Xu ◽

Cheng Zheng ◽

...

Keyword(s):

Signaling Pathway ◽

Expression Profiles ◽

Mapk Signaling ◽

Fat Deposition ◽

Mapk Signaling Pathway ◽

Differentially Expressed ◽

Rna Seq ◽

Domestic Pig ◽

Large White ◽

Pig Breeds

Abstract Background: Chinese domestic pig breeds are reputed for pork quality, but their low ratio of lean-to-fat carcass weight decreases production efficiency. A better understanding of the genetic regulation network of SC fat tissue is necessary for the rational selection of Chinese domestic pig breeds. In the present study, SC adipocytes were isolated from Jiaxing Black pigs (a Chinese indigenous pig breed with redundant SC fat deposition) and Large White pigs (a lean-type pig breed with relatively low SC fat deposition) and the expression profiles of mRNAs and lncRNAs were compared by RNA-seq analysis to identify biomarkers correlated with the differences of SC fat deposition between the two breeds.Results: A total of 3,371 differentially expressed genes (DEGs) and 1,182 differentially expressed lncRNAs (DELs) were identified in SC adipocytes between Jiaxing Black (JX) and Large White (LW) pigs, which included 797 upregulated mRNAs, 2,574 downregulated mRNAs, 461 upregulated lncRNAs and 721 downregulated lncRNAs. Gene Ontology and KEGG pathway analyses revealed that the DEGs and DELs were mainly involved in the immune response, cell fate determination, PI3K-Akt signaling pathway and MAPK signaling pathway, which are known to be related to adipogenesis and lipid metabolism. The expression levels of DEGs and DELs according to the RNA-seq data were verified by quantitative PCR, which showed 81.8% consistency. The differences in MAPK pathway activity between JX and LW pigs was confirmed by western blot analysis, with <100-fold elevated p38 phosphorylation in JX pigs.Conclusions: This study offers a detailed characterization of mRNAs and lncRNAs in fat- and lean-type pig breeds. The activity of the MAPK signaling pathway was found to be associated with subcutaneous adipogenesis. These results greatly enhance our understanding of the molecular mechanisms regulating SC fat deposition in pigs.

Download Full-text

RNA-Seq Reveals the Expression Profiles of Long Non-Coding RNAs in Lactating Mammary Gland from Two Sheep Breeds with Divergent Milk Phenotype

Animals ◽

10.3390/ani10091565 ◽

2020 ◽

Vol 10 (9) ◽

pp. 1565

Author(s):

Zhiyun Hao ◽

Yuzhu Luo ◽

Jiqing Wang ◽

Jiang Hu ◽

Xiu Liu ◽

...

Keyword(s):

Mammary Gland ◽

Signaling Pathway ◽

Target Genes ◽

Mammary Epithelial Cells ◽

Expression Profiles ◽

Wnt Signaling Pathway ◽

Mammary Epithelial ◽

Differentially Expressed ◽

Rna Seq ◽

Non Coding Rnas

Long non-coding RNAs (lncRNAs) are a kind of non-coding RNA with >200 nucleotides in length. Some lncRNAs have been proven to have clear regulatory functions in many biological processes of mammals. However, there have been no reports on the roles of lncRNAs in ovine mammary gland tissues. In the study, the expression profiles of lncRNAs were studied using RNA-Seq in mammary gland tissues from lactating Small-Tailed Han (STH) ewes and Gansu Alpine Merino (GAM) ewes with different milk yield and ingredients. A total of 1894 lncRNAs were found to be expressed. Compared with the GAM ewes, the expression levels of 31 lncRNAs were significantly up-regulated in the mammary gland tissues of STH ewes, while 37 lncRNAs were remarkably down-regulated. Gene Ontogeny (GO) enrichment and Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis found that the target genes of differentially expressed lncRNAs were enriched in the development and proliferation of mammary epithelial cells, morphogenesis of mammary gland, ErbB signaling pathway, and Wnt signaling pathway. Some miRNA sponges of differentially expressed lncRNAs, reported to be associated with lactation and mammary gland morphogenesis, were found in a lncRNA-miRNA network. This study reveals comprehensive lncRNAs expression profiles in ovine mammary gland tissues, thereby providing a further understanding of the functions of lncRNAs in the lactation and mammary gland development of sheep.

Download Full-text

A-To-I RNA Editing By ADAR1 Is Essential For Hematopoiesis

Blood ◽

10.1182/blood.v122.21.1199.1199 ◽

2013 ◽

Vol 122 (21) ◽

pp. 1199-1199 ◽

Cited By ~ 1

Author(s):

Brian Liddicoat ◽

Robert Piskol ◽

Alistair Chalk ◽

Miyoko Higuchi ◽

Peter Seeburg ◽

...

Keyword(s):

Gene Expression ◽

Rna Editing ◽

Fetal Liver ◽

Expression Profiles ◽

In Utero ◽

Tamoxifen Treatment ◽

Rna Seq ◽

Data Set ◽

Hematopoietic Stem

Abstract The role of RNA and its regulation is becoming increasingly appreciated as a vital component of hematopoietic development. RNA editing by members of the Adenosine Deaminase Acting on RNA (ADAR) gene family is a form of post-transcriptional modification which converts genomically encoded adenosine to inosine (A-to-I) in double-stranded RNA. A-to-I editing by ADAR directly converts the sequence of the RNA substrate and can alter the structure, function, processing, and localization of the targeted RNA. ADAR1 is ubiquitously expressed and we have previously described essential roles in the development of hematopoietic and hepatic organs. Germline ablation of murine ADAR1 results in a significant upregulation of interferon (IFN) stimulated genes and embryonic death between E11.5 and E12.5 associated with fetal liver disintegration and failed hemopoiesis. To determine the biological importance of A-to-I editing by ADAR1, we generated an editing dead knock-in allele of ADAR1 (ADAR1E861A). Mice homozygous for the ADAR1E861A allele died in utero at ∼E13.5. The fetal liver (FL) was small and had significantly lower cellularity than in controls. Analysis of hemopoiesis demonstrated increased apoptosis and a loss of hematopoietic stem cells (HSC) and all mature lineages. Most notably erythropoiesis was severely impaired with ∼7-fold reduction across all erythrocyte progenitor populations compared to controls. These data are consistent with our previous findings that ADAR1 is essential for erythropoiesis (unpublished data) and suggest that the ADAR1E861A allele phenocopies the null allele in utero. To assess the requirement of A-to-I editing in adult hematopoiesis, we generated mice where we could somatically delete the wild-type ADAR1 allele and leave only ADAR1E861A expressed in HSCs (hScl-CreERAdar1fl/E861A). In comparison to hScl-CreERAdar1fl/+ controls, hScl-CreERAdar1fl/E861A mice were anemic and had severe leukopenia 20 days post tamoxifen treatment. Investigation of marrow hemopoiesis revealed a significant loss of all cells committed to the erythroid lineage in hScl-CreERAdar1fl/E861A mice, despite having elevated phenotypic HSCs. Upon withdrawal of tamoxifen diet, all blood parameters were restored to control levels within 12 weeks owing to strong selection against cells expressing only the ADAR1E861A allele. To understand the mechanism through which ADAR1 mediated A-to-I editing regulates hematopoiesis, RNA-seq was performed. Gene expression profiles showed that a loss of ADAR1 mediated A-to-I editing resulted in a significant upregulation of IFN signatures, consistent with the gene expression changes in ADAR1 null mice. To define substrates of ADAR1 we assessed A-to-I mismatches in the RNA-seq data sets. 3,560 previously known and 353 novel A-to-I editing sites were identified in our data set. However, no single editing substrate discovered could account for the IFN signature observed or the lethality of ADAR1E861A/E861A mice. These results demonstrate that ADAR1 mediated A-to-I editing is essential for the maintenance of both fetal and adult hemopoiesis in a cell-autonomous manner and a key suppressor of the IFN response in hematopoiesis. Furthermore the ADAR1E861A allele demonstrates the essential role of ADAR1 in vivo is A-to-I editing. Disclosures: Hartner: TaconicArtemis: Employment.

Download Full-text

PADGE: analysis of heterogeneous patterns of differential gene expression

Physiological Genomics ◽

10.1152/physiolgenomics.00259.2006 ◽

2007 ◽

Vol 32 (1) ◽

pp. 154-159 ◽

Cited By ~ 15

Author(s):

Li Li ◽

Amitabha Chaudhuri ◽

John Chant ◽

Zhijun Tang

Keyword(s):

Gene Expression ◽

Differential Gene Expression ◽

Web Application ◽

Statistical Power ◽

Profile Analysis ◽

Expression Profiles ◽

Differentially Expressed ◽

Cancer Genes ◽

Data Set ◽

Differential Gene

We have devised a novel analysis approach, percentile analysis for differential gene expression (PADGE), for identifying genes differentially expressed between two groups of heterogeneous samples. PADGE was designed to compare expression profiles of sample subgroups at a series of percentile cutoffs and to examine the trend of relative expression between sample groups as expression level increases. Simulation studies showed that PADGE has more statistical power than t-statistics, cancer outlier profile analysis (COPA) (Tomlins SA, Rhodes DR, Perner S, Dhanasekaran SM, Mehra R, Sun XW, Varambally S, Cao X, Tchinda J, Kuefer R, Lee C, Montie JE, Shah RB, Pienta KJ, Rubin MA, Chinnaiyan AM. Science 310: 644–648, 2005), and kurtosis (Teschendorff AE, Naderi A, Barbosa-Morais NL, Caldas C. Bioinformatics 22: 2269–2275, 2006). Application of PADGE to microarray data sets in tumor tissues demonstrated its utility in prioritizing cancer genes encoding potential therapeutic targets or diagnostic markers. A web application was developed for researchers to analyze a large gene expression data set from heterogeneous biological samples and identify differentially expressed genes between subsets of sample classes using PADGE and other available approaches. Availability: http://www.cgl.ucsf.edu/Research/genentech/padge/ .

Download Full-text

Circular RNA Expression Profiles and Features in NAFLD Mice: A Study Using RNA-Seq Data

10.21203/rs.3.rs-58810/v1 ◽

2020 ◽

Author(s):

Xinlu Yuan ◽

Jianjun Diao ◽

Anqing Du ◽

Song Wen ◽

Ligang Zhou ◽

...

Keyword(s):

Expression Profile ◽

Noncoding Rna ◽

Expression Profiles ◽

Circular Rna ◽

Differentially Expressed ◽

Rna Seq ◽

Mice Model ◽

Specific Expression ◽

Nonalcoholic Fatty Liver ◽

Sequencing Method

Abstract Background: Nonalcoholic fatty liver disease (NAFLD) is primarily characterized by the hepatic cholesterol accumulation. Circular RNA (circRNA), one of noncoding RNA, involves in many liver diseases progression. However, no recent studies on circRNA expression profiles in NAFLD have been reported previously.Methods: A NAFLD mouse model was constructed by providing high-fat diet (HFD) for 32 weeks. The circRNAs expression profile in normal mice and NAFLD mice were determined using high-output RNA sequencing method and bioinformatics methods, while the differentially expressed circRNAs were confirmed using Sanger sequencing and qRT-PCR. The circRNA-miRNA network was also predicted. The biological functions of circRNAs were annotated by Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG).Results: The results demonstrated the successful construction of NAFLD mice model by immunohistology and serology assay. In total, 93 dysregulated circRNAs were observed, including 57 upregulated circRNAs and 36 downregulated circRNAs, in the NAFLD group. The circRNA-miRNA network revealed the complex interaction between circRNAs and its potential miRNA targets in NAFLD. The characteristic of tissue-specific expression in circRNA was demonstrated. The differentially expressed circRNAs with important biological function were also annotated using GO and KEGG. Both DDAH1 and VAV3 genes were found to be associated with the NAFLD development.Conclusions: Taken together, this study demonstrated the circRNAs expression profile and features in NAFLD, which may provide potential biological markers for the pathogenesis of NAFLD.

Download Full-text

Puccinia triticina pathotypes THTT and THTS display complex transcript profiles on wheat cultivar Thatcher

10.21203/rs.2.12168/v5 ◽

2020 ◽

Author(s):

jie Wei ◽

Liping Cui ◽

na Zhang ◽

dongdong Du ◽

qingfang Meng ◽

...

Keyword(s):

Transcription Initiation ◽

Wheat Cultivar ◽

Early Stage ◽

Puccinia Triticina ◽

Differentially Expressed ◽

Post Inoculation ◽

Rna Seq ◽

Initiation Factors ◽

Transcript Profiles ◽

Important Disease

Abstract Background: Wheat leaf rust is an important disease worldwide. Understanding the pathogenic molecular mechanism of Puccinia triticina Eriks . ( Pt ) and the inconstant toxic region is critical for managing the disease. The present study aimed to analyze the pathogenic divergence between Pt isolates. Results: Total RNA was extracted from the wheat cultivar Thatcher infected by two Pt isolates, Tc361_1 (THTT) and Tc284_2 (THTS), at 144 hours post inoculation (hpi). The mRNA was then sequenced, and a total of 2,784 differentially expressed genes (DEGs) were detected. Forty-five genes were specifically expressed in THTT; these genes included transcription initiation factors and genes with transmembrane transporter activity and other genes. Twenty-six genes were specifically expressed in THTS, including genes with GTPase activity, ABC transporters and other genes. Fifty-four differentially expressed candidate effectors were screened from the two isolates. Two candidate effectors were chosen and validated on tobacco, and the results showed that they could inhibit necrosis induced by BAX. qRT-PCR of 12 significant DEGs was carried out to validate that the results are similar to those of RNA-seq at 144 hpi, to show the expression levels of these DEGs in the early stage and to elucidate the differences in expression between the two Pt pathotypes. Conclusion: The results obtained in this study showed that although the two pathotypes of THTT and THTS contribute similar virulence to wheat, there are a large number of genes participate in the interaction with the susceptible wheat cultivar Thatcher, and revealed the pathogenicity of rust is very complicated.

Download Full-text

RNA-seq reveals few differences in resistant and susceptible responses of barley to infection by the spot blotch pathogenBipolaris sorokiniana

10.1101/384529 ◽

2018 ◽

Cited By ~ 1

Author(s):

Matthew Haas ◽

Martin Mascher ◽

Claudia Castell-Miller ◽

Brian J. Steffenson

Keyword(s):

Expression Profiles ◽

Durable Resistance ◽

Defense Responses ◽

Spot Blotch ◽

Seedling Stage ◽

Differentially Expressed ◽

Adult Plant ◽

Rna Seq ◽

Adult Plant Stage ◽

Plant Stage

AbstractSpot blotch, caused byBipolaris sorokiniana(Sacc.) Shoem., is an economically important disease affecting barley (Hordeum vulgareL.). The disease has largely been controlled in the Upper Midwest region of the USA through a suite of quantitative trait loci (QTL) termed the Midwest Six-rowed Durable Resistance Haplotype (MSDRH). These QTL have been bred into all six-rowed Midwest barley cultivars, including the widely used cultivar Morex. We identified a gamma ray- induced Morex mutant (MUT) that exhibits spot blotch susceptibility at the seedling stage. This mutant also spontaneously develops extremely large necrotic lesions in the absence of the pathogen at the adult plant stage. Spot blotch susceptibility at the seedling stage and necrotic lesion formation at the adult plant stage are highly correlated. To start dissecting the molecular responses underlying the observed symptoms at the seedling stage, we conducted a time course RNA-seq experiment comparing the wild type (WT) and the mutant (MUT) Morex at 12, 24 and 36 h afterB. sorokinianainoculation. Mock-inoculated controls were also included. A total of 10,772 and 11,530 genes were differentially expressed between treatments for WT and MUT genotypes, respectively, while 277 and 195 genes were differentially expressed between fungal and mock-inoculated genotypes, respectively. The transcript expression profiles of WT and MUT Morex samples were similar for most treatments. Two genes whose expression was putatively knocked out in the MUT were identified: HORVU3Hr1G019920 (glycine-rich protein) and HORVU5Hr1G120850 (Long- chain-fatty-acid—CoA ligase 1). The latter appears to be genetically intact, but not expressed. Collectively, these data suggest that MUT susceptibility toB. sorokinianais a result of minor, rather than major, differences in the defense responses.

Download Full-text