Expression of genes in gastrointestinal and lymphatic tissues during parasite infection in sheep genetically resistant or susceptible to Trichostrongylus colubriformis and Haemonchus contortus

Resistance to an acute gastrointestinal nematode (GIN) infection is dependent on the ability of the host to recognise the parasite and mount a protective Th2 response. It is hypothesised that lambs which are genetically susceptible to GIN will differentially up-regulate Th1-type genes and therefore remain susceptible to chronic parasitism compared with genetically resistant lambs which will differentially up-regulate Th2-type genes and clear the parasite infection. Two selection flocks, in which lines of Merino sheep produced lambs genetically resistant or susceptible to GIN, were acutely challenged once or thrice with either Haemonchus contortus or Trichostrongylus colubriformis. Faecal-egg counts (FECs), and plasma and tissue anti-parasite (H. contortus or T. colubriformis) antibody isotype responses showed that resistant animals challenged three times with T. colubriformis established a protective Th2 response (negligible FEC, IgG1 and IgE) whereas susceptible animals required multiple challenges to establish a significant IgG1 response despite FECs remaining high. Trichostrongylus colubriformis elicited a more pronounced host response than H. contortus. RNA extracted from tissues at the site of each parasite infection and associated lymph nodes were interrogated by microarray and quantitative PCR analyses to correlate host gene expression to FECs and antibody responses. The IFN-γ inducible gene cxcl10 was up-regulated in the susceptible line of the Trichostrongylus selection flock sheep after a tertiary challenge with the parasites H. contortus and T. colubriformis. However, a uniform pattern of genes was not up-regulated in resistant animals from both selection flocks during both parasite infections, suggesting that the mode of host resistance to these parasites is different, although some similarities in host susceptibility were apparent.     

Proteomic analysis of the abomasal mucosal response following infection by the nematode, Haemonchus contortus, in genetically resistant and susceptible sheep

Sheep have a variable ability to resist gastrointestinal nematode infection, but the key factors mediating this response are poorly defined. Here we report the first large-scale application of quantitative proteomic technologies to define proteins that are differentially abundant between sheep selectively bred to have an enhanced (resistant) or reduced (susceptible) ability to eliminate nematodes. Samples were collected from the abomasal mucosa three days after experimental challenge with the nematode, Haemonchus contortus. This timing reflects the initial interaction of host and parasite, and the tissue represents the immediate interface. We identified and quantified more than 4400 unique proteins, of which 158 proteins showed >1.5 fold difference between the resistant and susceptible sheep. Trefoil factor 2, a member of RAS oncogene family (RAP1A) and ring finger protein 126 were amongst the proteins found to be highly abundant in the abomasal surface of resistant sheep, whereas adenosine deaminase and the gastrokine-3 like precursor were found at higher levels in susceptible sheep. Construction of gut proteome interaction networks identified mitochondrial function and energetic partitioning as important components of an effective nematode eliminating response. The differentially abundant proteins may be useful targets for phenotypic tests that aim to identify sheep with an enhanced ability to resist nematode infection.     

Identification of potential functional variants underlying ovine resistance to gastrointestinal nematode infection by using RNA-Seq

In dairy sheep flocks from Mediterranean countries, replacement and adult ewes are the animals most affected by gastrointestinal nematode (GIN) infections. In this study, we have exploited the information derived from an RNA-Seq experiment with the aim of identifying potential causal mutations related to GIN resistance in sheep. Considering the RNA-Seq samples from 12 ewes previously classified as six resistant and six susceptible animals to experimental infection by Teladorsagia circumcincta, we performed a variant calling analysis pipeline using two different types of software, gatk version 3.7 and Samtools version 1.4. The variants commonly identified by the two packages (high-quality variants) within two types of target regions – (i) QTL regions previously reported in sheep for parasite resistance based on SNP-chip or sequencing technology studies and (ii) functional candidate genes selected from gene expression studies related to GIN resistance in sheep – were further characterised to identify mutations with a potential functional impact. Among the genes harbouring these potential functional variants (930 and 553 respectively for the two types of regions), we identified 111 immune-related genes in the QTL regions and 132 immune-related genes from the initially selected candidate genes. For these immune-related genes harbouring potential functional variants, the enrichment analyses performed highlighted significant GO terms related to apoptosis, adhesion and inflammatory response, in relation to the QTL related variants, and significant disease-related terms such as inflammation, adhesion and necrosis, in relation to the initial candidate gene list. Overall, the study provides a valuable list of potential causal mutations that could be considered as candidate causal mutations in relation to GIN resistance in sheep. Future studies should assess the role of these suggested mutations with the aim of identifying genetic markers that could be directly implemented in sheep breeding programmes considering not only production traits, but also functional traits such as resistance to GIN infections.     

Gene copy number variations as signatures of adaptive evolution in the parthenogenetic,plant‐parasitic nematode Meloidogyne incognita

Adaptation to changing environmental conditions represents a challenge to parthenogenetic organisms, and until now, how phenotypic variants are generated in clones in response to the selection pressure of their environment remains poorly known. The obligatory parthenogenetic root‐knot nematode species Meloidogyne incognita has a worldwide distribution and is the most devastating plant‐parasitic nematode. Despite its asexual reproduction, this species exhibits an unexpected capacity of adaptation to environmental constraints, for example, resistant hosts. Here, we used a genomewide comparative hybridization strategy to evaluate variations in gene copy numbers between genotypes of M. incognita resulting from two parallel experimental evolution assays on a susceptible vs. resistant host plant. We detected gene copy number variations (CNVs) associated with the ability of the nematodes to overcome resistance of the host plant, and this genetic variation may reflect an adaptive response to host resistance in this parthenogenetic species. The CNV distribution throughout the nematode genome is not random and suggests the occurrence of genomic regions more prone to undergo duplications and losses in response to the selection pressure of the host resistance. Furthermore, our analysis revealed an outstanding level of gene loss events in nematode genotypes that have overcome the resistance. Overall, our results support the view that gene loss could be a common class of adaptive genetic mechanism in response to a challenging new biotic environment in clonal animals.     

Digital gene expression analysis of gastrointestinal helminth resistance in Scottish blackface lambs

Digital gene expression (DGE) analysis offers a route to gene discovery which by-passes the need to develop bespoke arrays for nonmodel species, and is therefore a potentially valuable tool for molecular ecologists. Scottish blackface sheep, which vary in resistance to the common abomasal parasitic nematode Teladorsagia circumcincta, were trickle-infected with L3 larvae over 3 months to mimic the natural progression of infection. DGE was performed on abomasal lymph node tissue after the resolution of infection in resistant animals. Susceptible (low resistance) animals showed a large number of differentially expressed genes associated with inflammation and cell activation, but generally few differentially regulated genes in either the susceptible or the resistant group were directly involved in the adaptive immune function. Our results are consistent with the hypothesis that both resistance and susceptibility are active responses to infection and that susceptibility is associated with dysfunction in T cell differentiation and regulation.     

Heterologous expression of the Mi-1.2 gene from tomato confers resistance against nematodes but not aphids in eggplant

The Mi-1.2 gene in tomato (Solanum lycopersicum) is a member of the nucleotide-binding leucine-rich repeat (NBLRR) class of plant resistance genes, and confers resistance against root-knot nematodes (Meloidogyne spp.), the potato aphid (Macrosiphum euphorbiae), and the sweet potato whitefly (Bemisia tabaci). Mi-1.2 mediates a rapid local defensive response at the site of infection, although the signaling and defensive pathways required for resistance are largely unknown. In this study, eggplant (S. melongena) was transformed with Mi-1.2 to determine whether this gene can function in a genetic background other than tomato. Eggplants that carried Mi-1.2 displayed resistance to the root-knot nematode Meloidogyne javanica but were fully susceptible to the potato aphid, whereas a susceptible tomato line transformed with the same transgene was resistant to nematodes and aphids. This study shows that Mi-1.2 can confer nematode resistance in another Solanaceous species. It also indicates that the requirements for Mi-mediated aphid and nematode resistance differ. Potentially, aphid resistance requires additional genes that are not conserved between tomato and eggplant.     

Evaluation of reference genes for real-time PCR studies of Brazilian Somalis sheep infected by gastrointestinal nematodes

Precise normalization with reference genes is necessary, in order to obtain reliable relative expression data in response to gastrointestinal nematode infection. By using sheep from temperate regions as models, three reference genes, viz., ribosomal protein LO (RPLO), glyceraldehyde 3-phosphate dehydrogenase (GAPDH) and succinate dehydrogenase complex subunit A (SDHA), were investigated in the abomasum, abomasal lymph nodes and small intestine of Brazilian Somalis sheep, either resistant or susceptible to gastrointestinal nematodes infections. Real time PCR was carried out by using SYBR Green I dye, and gene stability was tested by geNorm. RPLO was an ideal reference gene, since its expression was constant across treatments, presented lower variation, and was ranked as the most stable in abomasum and lymph node tissues. On the other hand, SDHA was the most stable in the small intestine followed by RPLO and GAPDH. These findings demonstrate the importance of correctly choosing reference genes prior to relative quantification. In addition, we determined that reference genes used in sheep from temperate regions, when properly tested, can be applied in animals from tropical regions such as the Brazilian Somalis sheep.     

Possibilities to breed for resistance to nematode parasite infections in small ruminants in tropical production systems

Gastrointestinal nematode parasitism is the most important disease affecting livestock production systems in developing countries, particularly small ruminant production systems. Of particular importance are infections with the strongyle Haemonchus contortus. Integrated disease control strategies are required, including improved management, nutrition and wise use of anthelmintic chemicals. Increasingly, selection of sheep or goats for improved nematode resistance is viewed as a valuable option to complement other control measures. Breeding for resistance is possible because of the existence of extensive genetic variation in resistance, both within and between breeds of sheep and goats. Such breeding schemes are most likely to be based on choice of appropriate breeds adapted to the local environmental conditions, followed by phenotypic selection for resistance. Goal and selection objective traits are likely to include performance (e.g. growth rate) under conditions of parasite challenge, faecal egg count (FEC) and measures of anaemia. With current technologies, genetic markers are likely to be too expensive and logistically difficult to incorporate into breeding schemes in tropical or developing countries. Genotype by environment interactions may be expected, particularly when comparing animals in environments that differ in the extent of parasite challenge or differ in the quality of available nutrition. However, there is no reason to expect antagonistic genetic relationships between performance and resistance, and selection indices should be readily constructed that improve both performance and resistance. If FEC is decreased, then pasture contamination should also decrease, leading to additional benefits for all sheep grazing the same pasture. Finally, breeding for nematode resistance should lead to lasting and sustained improvements in resistance or tolerance. There is no empirical evidence to suggest that nematodes will evolve rapidly in response to resistant hosts, and mathematical models based on genetic and biological principles also suggest that resistance should be sustainable.     

Nematospiroides dubius: two H-2-linked genes influence levels of resistance to infection in mice

Strains of mice sharing common H-2 haplotypes but different genetic backgrounds, and H-2 congenic strains of mice differing only at H-2 genes were studied to assess the role of H-2 and non-H-2 genes in immunity to challenge infections with the nematode parasite Nematospiroides dubius. Strains of mice sharing the H-2k haplotype were uniformly more susceptible to challenge than strains expressing H-2q alleles, regardless of genetic background. However, in some cases strains of mice sharing the k or q haplotypes differed significantly in levels of resistance. Therefore, non-H-2 genes must influence the response observed. H-2 cogenic strains of mice differed markedly in their ability to resist challenge infections. Mice sharing the C57BL/10 background but expressing k alleles were very susceptible to challenge, while the H-2q, H-2f, and H-2s, haplotypes were associated with resistance. Studies of H-2 congenic recombinant strains of mice suggested that two H-2 genes influence the antiparasite response. One of these genes maps to the left of E alpha and the other to the D-end of the H-2 complex. It is concluded also that the unique configuration of H-2 genes in F1 hybrids contributes to increased resistance to challange.     

Catalase gene is associated with facial eczema disease resistance in sheep

Facial eczema (FE) is a hepatogenous photosensitization disease of ruminant animals, particularly in sheep which vary widely in their susceptibility to the disease. The liver damage is caused by the mycotoxin, sporidesmin. There is evidence that the toxicity of sporidesmin is due to its ability to generate 'active oxygen' species. We evaluated the catalase gene, which encodes an enzyme with antioxidant functions, as a candidate for determining the susceptibility of sheep to the disease. Two microsatellite markers, OarSHP3 and OarSHP4, which flank the sheep catalase gene, were isolated from a Yeast Artificial Chromosome (YAC) clone. These markers mapped the catalase locus by linkage to ovine chromosome 15. Eleven informative markers spaced throughout chromosome 15, inclusive of the catalase marker OarSHP4, gave no significant linkage with the disease traits when analysed in four outcross resource pedigrees. However, OarSHP3 and OarSHP4 allele frequencies showed significant differences between FE resistant and susceptible selection-lines. Comparison of sequences of catalase cDNAs from sheep of resistant and susceptible lines showed only two silent mutations. A single nucleotide polymorphisms (KP1) in exon 6 of the catalase gene also showed significant differences in allele frequencies between the selection lines. The lack of evidence for linkage in outcross pedigrees, but the significant association in the genetic lines, implies that catalase is involved in determining the susceptibility of sheep to facial eczema, and that the candidate gene's effect is probably recessive or minor.     

A strategy to identify positional candidate genes conferring Marek''s disease resistance by integrating DNA microarrays and genetic mapping

Marker-assisted selection (MAS) to enhance genetic resistance to Marek's disease (MD), a herpesvirus-induced T cell cancer in chicken, is an attractive alternative to augment control with vaccines. Our earlier studies indicate that there are many quantitative trait loci (QTL) containing one or more genes that confer genetic resistance to MD. Unfortunately, it is difficult to sufficiently resolve these QTL to identify the causative gene and generate tightly linked markers. One possible solution is to identify positional candidate genes by virtue of gene expression differences between MD resistant and susceptible chicken using deoxyribonucleic acid (DNA) microarrays followed by genetic mapping of the differentially-expressed genes. In this preliminary study, we show that DNA microarrays containing approximately 1200 genes or expressed sequence tags (ESTs) are able to reproducibly detect differences in gene expression between the inbred ADOL lines 63 (MD resistant) and 72 (MD susceptible) of uninfected and Marek's disease virus (MDV)-infected peripheral blood lymphocytes. Microarray data were validated by quantitative polymerase chain reaction (PCR) and found to be consistent with previous literature on gene induction or immune response. Integration of the microarrays with genetic mapping data was achieved with a sample of 15 genes. Twelve of these genes had mapped human orthologues. Seven genes were located on the chicken linkage map as predicted by the human-chicken comparative map, while two other genes defined a new conserved syntenic group. More importantly, one of the genes with differential expression is known to confer genetic resistance to MD while another gene is a prime positional candidate for a QTL.     

Identification and characterization of salt responsive miRNA-SSR markers in rice (Oryza sativa)

Salinity is an important abiotic stress that affects agricultural production and productivity. It is a complex trait that is regulated by different molecular mechanisms. miRNAs are non-coding RNAs which are highly conserved and regulate gene expression. Simple sequence repeats (SSRs) are robust molecular markers for studying genetic diversity. Although several SSR markers are available now, challenge remains to identify the trait-specific SSRs which can be used for marker assisted breeding. In order to understand the genetic diversity of salt responsive-miRNA genes in rice, SSR markers were mined from 130 members of salt-responsive miRNA genes of rice and validated among the contrasting panels of tolerant as well as susceptible rice genotypes, each with 12 genotypes. Although 12 miR-SSRs were found to be polymorphic, only miR172b-SSR was able to differentiate the tolerant and susceptible genotypes in 2 different groups. It had also been found that miRNA genes were more diverse in susceptible genotypes than the tolerant one (as indicated by polymorphic index content) which might interfere to form the stem-loop structure of premature miRNA and their subsequent synthesis in susceptible genotypes. Thus, we concluded that length variations of the repeats in salt responsive miRNA genes may be responsible for a possible sensitivity to salinity adaptation. This is the first report of characterization of trait specific miRNA derived SSRs in plants.     

Salmonella carrier state in chicken: comparison of expression of immune response genes between susceptible and resistant animals

Asymptomatic Salmonella enterica serovar Enteritidis carrier state in poultry has serious consequences on food safety and public health due to the risks of food poisoning following consumption of contaminated products. An understanding the mechanisms of persistence of Salmonella in the digestive tract of chicken can be achieved by a better knowledge of the defects in the control of infection in susceptible versus resistant animals. The gene expression of innate immune response factors including anti-microbial molecules, inflammatory and anti-infectious cytokines was studied in the caecal lymphoid tissue associated with the carrier state. Expression levels of these genes were assessed by real-time PCR and were compared in two inbred lines of chickens differing in resistance to the carrier state following oral inoculation of S. enterica serovar Enteritidis at 1 week of age. No correlation was observed between resistance/susceptibility to caecal carrier state and level of interleukin (IL)-1beta, IL-8, IL-18, inducible NO synthase (iNOS) and natural resistance associated macrophage protein 1 (NRAMP1). A high baseline level of defensin gene expression was recorded in young animals from the susceptible line. In contrast, a significantly low expression of interferon-gamma (IFN-gamma) gene was observed in these susceptible infected animals in comparison to resistant ones and healthy counterparts. IFN-gamma expression level represents a valuable indication of immunodeficiency associated with persistence of Salmonella in the chicken digestive tract, and IFN-gamma thus represents a factor to consider in the development of prophylactic measures for the reduction of Salmonella carrier state.     

A genomics-informed,SNP association study reveals FBLN1 and FABP4 as contributing to resistance to fleece rot in Australian Merino sheep

Background  

Fleece rot (FR) and body-strike of Merino sheep by the sheep blowfly Lucilia cuprina are major problems for the Australian wool industry, causing significant losses as a result of increased management costs coupled with reduced wool productivity and quality. In addition to direct effects on fleece quality, fleece rot is a major predisposing factor to blowfly strike on the body of sheep. In order to investigate the genetic drivers of resistance to fleece rot, we constructed a combined ovine-bovine cDNA microarray of almost 12,000 probes including 6,125 skin expressed sequence tags and 5,760 anonymous clones obtained from skin subtracted libraries derived from fleece rot resistant and susceptible animals. This microarray platform was used to profile the gene expression changes between skin samples of six resistant and six susceptible animals taken immediately before, during and after FR induction. Mixed-model equations were employed to normalize the data and 155 genes were found to be differentially expressed (DE). Ten DE genes were selected for validation using real-time PCR on independent skin samples. The genomic regions of a further 5 DE genes were surveyed to identify single nucleotide polymorphisms (SNP) that were genotyped across three populations for their associations with fleece rot resistance.     

A Study of Gene Expression in the Nematode Resistant Wild Peanut Relative, Arachis stenosperma, in Response to Challenge with Meloidogyne arenaria

Peanut (Arachis hypogaea) is amongst the most important legume crops in the world. One of its main yield constraints is the root-knot nematode Meloidogyne arenaria. A number of wild Arachis species, including A. stenosperma, are resistant to nematodes, and are a potential source of new resistance alleles for cultivated peanut. Using in silico subtraction of ESTs and macroarray analysis, we identified genes differentially expressed in A. stenosperma roots during its resistance response to M. arenaria. The three most differentially expressed genes [Auxin Repressed Protein (AsARP), Cytokinin Oxidase (AsCKX) and Metallothionein Type 2 (AsMET2)] were further analyzed using northern-blot and showed distinct expression profiles in the resistant A. stenosperma and susceptible A. hypogaea, both after, and sometimes even before, challenge with nematodes. Of the three most differentially expressed genes, AsARP and AsCKX are potentially involved in plant hormonal balance, and AsMET2 may be related to the reactive oxygen reaction triggered by the hypersensitive response (HR).