Genomic analyses of Haemonchus contortus infection in sheep: abomasal fistulation and two Haemonchus strains do not substantially confound host gene expression in microarrays

From innate to adaptive immunity: Abomasal transcriptomic responses of merino sheep to Haemonchus contortus infection

Although many important mediators and critical pathways are found to be involved in host immune responses to Haemonchus contortus infection, the initial responses to infection in the naïve and in the previously exposed state have not been compared at the transcriptional level. To further understand the development of adaptive immunity to H. contortus infection, we compared the early abomasal gene expression patterns between a primary and a tertiary challenge for four lines of sheep to discover differentially expressed genes (DEGs). The sheep were from the resistant (R) and susceptible (S) lines of two flocks of sheep selected for divergent responses to gastro-intestinal parasites (HSF and TSF). The flocks have separate origins and were initiated using two different strains of Merino sheep. One of the DEGs, mast cell proteinase 1, had significantly lower expression in tertiary compared to primary infections for all four lines of sheep. This gene was not identified in previous studies where resistant and susceptible sheep samples were compared within infection time points. Comparing the differentially expressed genes (DEGs) for the two R lines reveals that responses differed very little between the primary and tertiary challenges for HSFR and only two genes were identified, in contrast to the TSFR where there were 134 genes identified including the two identified using the HSFR animals. Similarly, comparing the primary and tertiary challenges for HSFS identified 15 DEGs, whilst for TSFS there were 128 DEGs identified. It is surprising that so few genes respond similarly between the two challenge regimes across the four lines of sheep, and suggests significant differences in immune mechanisms between the two flocks (across the lines) and also between the lines within flocks. Our results offer a quantitative snapshot comparing the transcriptome in the ovine abomasum between primary and tertiary infections with H. contortus in both genetically resistant and susceptible sheep.

The potential for vaccines against scour worms of small ruminants

This review addresses the research landscape regarding vaccines against scour worms, particularly Trichostrongylus spp. and Teladorsagia circumcincta. The inability of past research to deliver scour-worm vaccines with reliable and reproducible efficacy has been due in part to gaps in knowledge concerning: (i) host-parasite interactions leading to development of type-2 immunity, (ii) definition of an optimal suite of parasite antigens, and (iii) rational formulation and administration to induce protective immunity against gastrointestinal nematodes (GIN) at the site of infestation. Recent 'omics' developments enable more systematic analyses. GIN genomes are reaching completion, facilitating "reverse vaccinology" approaches that have been used successfully for the Rhipicephalus australis vaccine for cattle tick, while methods for gene silencing and editing in GIN enable identification and validation of potential vaccine antigens. We envisage that any efficacious scour worm vaccine(s) would be adopted similarly to "Barbervax™" within integrated parasite management schemes. Vaccines would therefore effectively parallel the use of resistant animals, and reduce the frequency of drenching and pasture contamination. These aspects of integration, efficacy and operation require updated models and validation in the field. The conclusion of this review outlines an approach to facilitate an integrated research program.

Dynamic transcriptomic changes of goat abomasal mucosa in response to Haemonchus contortus infection

Gastrointestinal nematode (GIN) infections are one of the major constraints for grazing sheep and goat production worldwide. Genetic selection for resistant animals is a promising control strategy. Whole-transcriptome analysis via RNA-sequencing (RNA-seq) provides knowledge of the mechanisms responsible for complex traits such as resistance to GIN infections. In this study, we used RNA-seq to monitor the dynamics of the response of the abomasal mucosa of Creole goat kids infected with Haemonchus contortus by comparing resistant and susceptible genotypes. A total of 8 cannulated kids, 4 susceptible and 4 resistant to GIN, were infected twice with 10 000 L3 H. contortus. During the second infection, abomasal mucosal biopsies were collected at 0, 8, 15 and 35 days post-infection (dpi) from all kids for RNA-seq analysis. The resistant animals showed early activation of biological processes related to the immune response. The top 20 canonical pathways of differentially expressed genes for different comparison showed activation of the immune response through many relevant pathways including the Th1 response. Interestingly, our results showed a simultaneous time series activation of Th2 related genes in resistant compared to susceptible kids.

Transcriptome variation in response to gastrointestinal nematode infection in goats

Gastrointestinal nematodes (GIN) are a major constraint for small ruminant production. Due to the rise of anthelmintic resistance throughout the world, alternative control strategies are needed. The development of GIN resistance breeding programs is a promising strategy. However, a better understanding of the mechanisms underlying genetic resistance might lead to more effective breeding programmes. In this study, we compare transcriptome profiling of abomasal mucosa and lymph node tissues from non-infected, resistant and susceptible infected Creole goats using RNA-sequencing. A total of 24 kids, 12 susceptible and 12 GIN resistant based on the estimated breeding value, were infected twice with 10,000 L3 Haemonchus contortus. Physiological and parasitological parameters were monitored during infection. Seven weeks after the second infection, extreme kids (n = 6 resistant and 6 susceptible), chosen on the basis of the fecal egg counts (FEC), and 3 uninfected control animals were slaughtered. Susceptible kids had significantly higher FEC compared with resistant kids during the second infection with no differences in worm burden, male and female worm count or establishment rate. A higher number of differentially expressed genes (DEG) were identified in infected compared with non-infected animals in both abomasal mucosa (792 DEG) and lymph nodes (1726 DEG). There were fewer DEG in resistant versus susceptible groups (342 and 450 DEG, in abomasal mucosa and lymph nodes respectively). ‘Cell cycle’ and ‘cell death and survival’ were the main identified networks in mucosal tissue when comparing infected versus non-infected kids. Antigen processing and presentation of peptide antigen via major histocompatibility complex class I were in the top biological functions for the DEG identified in lymph nodes. The TGFβ1 gene was one of the top 5 upstream DEG in mucosal tissue. Our results are one of the fist investigating differences in the expression profile induced by GIN infection in goats.

Iron supplementation for lambs experimentally infected by Haemonchus contortus: response to anemia and iron store in the bone marrow

ABSTRACT: To evaluate the effect of supplementation of iron dextran on blood variables and iron metabolism in lambs experimentally infected by Haemonchus contortus, four experimental groups were used: uninfected and non-supplemented animals (GI); infected animals supplemented with iron (GII); uninfected animals supplemented with iron (GIII); and infected non-supplemented animals (GIV). Groups II and IV received 10,000 larvae (L3) of Haemonchus contortus, and groups II and III received three doses of iron dextran (20mg/kg) intramuscularly with seven days of interval. Blood and faeces samples were collected on days 10 (D10), 17 (D17), 24 (D24), and 31 (D31), in order to determine red blood cell counts, iron metabolism, and EPG. Infected animals developed anemia from D24 and anemia was more severe on D31. Animals from GII had higher hematocrit and hemoglobin concentration compared to animals of GIV on D31. Iron stores in the bone marrow were higher in GII and GIII compared to GI and GIV. The GIV showed lower seric levels of iron on D24 compared to the other groups. The iron supplementation reduces the severity of the anemia caused by infection with Haemonchus contortus in lambs, improving erythropoietic response after blood loss.

Immunity to Haemonchus contortus and Vaccine Development

Sheep are capable of developing protective immunity to Haemonchus contortus through repeated exposure to this parasite, although this immune protection is the result of a complex interaction among age, gender, physiological status, pregnancy, lactation, nutrition and innate and adaptive immunity in the host animal. There are multiple effectors of the protective immune response, which differ depending on the developmental stage of the parasite being targeted, and our understanding of the effector mechanisms has developed considerably in the 2000s. The rational design of vaccines based on 'natural' or 'exposed' antigens depends on an understanding of this exposure-induced immunity. However, the most effective current vaccines rely on protection via the induction of high circulating antibody levels to 'hidden' gut antigens of H. contortus. The success of this latter strategy has resulted in the launch of a vaccine, which is based on extracts of the parasite's gut, to aid in the control of Haemonchus in Australia. The development of recombinant subunit vaccines based on the components of the successful native vaccine has not yet been achieved and most of the recent successes with recombinant subunit vaccines have focussed on antigens unrelated to the gut antigens. The future integration of an understanding of the immunobiology of this parasite with advances in antigen identification, expression (or synthesis) and presentation is likely to be pivotal to the further development of these recombinant subunit vaccines. Recent progress in each of the components underpinning this integrated approach is summarized in this review.

Variations in T Cell Transcription Factor Sequence and Expression Associated with Resistance to the Sheep Nematode Teladorsagia circumcincta

This study used selected lambs that varied in their resistance to the gastrointestinal parasite Teladorsagia circumcincta. Infection over 12 weeks identified susceptible (high adult worm count, AWC; high fecal egg count, FEC; low body weight, BW; low IgA) and resistant sheep (no/low AWC and FEC, high BW and high IgA). Resistance is mediated largely by a Th2 response and IgA and IgE antibodies, and is a heritable characteristic. The polarization of T cells and the development of appropriate immune responses is controlled by the master regulators, T-bet (TBX21), GATA-3 (GATA3), RORγt (RORC2) and RORα (RORA); and several inflammatory diseases of humans and mice are associated with allelic or transcript variants of these transcription factors. This study tested the hypothesis that resistance of sheep to T. circumcincta is associated with variations in the structure, sequence or expression levels of individual master regulator transcripts. We have identified and sequenced one variant of sheep TBX21, two variants of GATA3 and RORC2 and five variants of RORA from lymph node mRNA. Relative RT-qPCR analysis showed that TBX21, GATA3 and RORC2 were not significantly differentially-expressed between the nine most resistant (AWC, 0; FEC, 0) and the nine most susceptible sheep (AWC, mean 6078; FEC, mean 350). Absolute RT-qPCR on all 45 animals identified RORVv5 as being significantly differentially-expressed (p = 0.038) between resistant, intermediate and susceptible groups; RORCv2 was not differentially-expressed (p = 0.77). Spearman’s rank analysis showed that RORAv5 transcript copy number was significantly negatively correlated with parameters of susceptibility, AWC and FEC; and was positively correlated with BW. RORCv2 was not correlated with AWC, FEC or BW but was significantly negatively correlated with IgA antibody levels. This study identifies the full length RORA variant (RORAv5) as important in controlling the protective immune response to T. circumcincta infection in sheep.

Variable exon usage of differentially-expressed genes associated with resistance of sheep to Teladorsagia circumcincta

The resistance and susceptibility of sheep to the common abomasal nematode parasite, Teladorsagia circumcincta is strongly associated with the differential polarization of the immune response. Resistant animals control larval colonization by the production of a protective antibody response regulated by Th2 T cells. Susceptible sheep respond to infection by developing an inflammatory Th1/Th17 response that fails to control infection. Previous microarray analysis identified genes associated with T cell polarization that were differentially expressed between the resistant and susceptible sheep. RT-qPCR confirmed the microarray data for ALOX15 and IL13. Both ALOX15 exon 9 and IL13 exon 4 were significantly increased in resistant animals and copy number RT-qPCR showed that expression levels of these exons were significantly negatively correlated with quantitative phenotypic traits, including abomasal worm counts and faecal egg counts. Sequencing of the intronic regions 5' to these genes failed to identify any potential genetic links to differential exon usage.

Oxidative stress by Haemonchus contortus in lambs: Influence of treatment with zinc edetate

The aim of the present study was to assess the effects of zinc edetate on the oxidative stress of lambs infected by Haemonchus contortus. Twenty-four lambs were allocated into four groups: Group I--uninfected animals; Group II--uninfected animals treated subcutaneously with zinc edetate; Group III--animals infected by H. contortus and Group IV--animals infected and treated. The oxidative stress index (OSI) and the eggs per gram of feces (EPG) were assessed after 10, 17, 24, 31 and 38 days post-infection. Based on the EPG and the quantity of adult H. contortus, the infection did not differ between groups III and IV. Zinc edetate reduced the OSI in Group IV in relation to Group I after 24 days post-infection, and in relation to group III after 31 days post-infection. Treatment with zinc edetate could help reduce the oxidative stress induced by H. contortus in lambs.

Identification of novel loci associated with gastrointestinal parasite resistance in a Red Maasai x Dorper backcross population

Gastrointestinal (GI) parasitic infection is the main health constraint for small ruminant production, causing loss of weight and/or death. Red Maasai sheep have adapted to a tropical environment where extreme parasite exposure is a constant, especially with highly pathogenic Haemonchus contortus. This breed has been reported to be resistant to gastrointestinal parasite infection, hence it is considered an invaluable resource to study associations between host genetics and resistance. The aim of this study was to identify polymorphisms strongly associated with host resistance in a double backcross population derived from Red Maasai and Dorper sheep using a SNP-based GWAS analysis. The animals that were genotyped represented the most resistant and susceptible individuals based on the tails of phenotypic distribution (10% each) for average faecal egg counts (AVFEC). AVFEC, packed cell volume (AVPCV), and live weight (AVLWT) were adjusted for fixed effects and co-variables, and an association analysis was run using EMMAX. Revised significance levels were calculated using 100,000 permutation tests. The top five significant SNP markers with - log10 p-values >3.794 were observed on five different chromosomes for AVFEC, and BLUPPf90/PostGSf90 results confirmed EMMAX significant regions for this trait. One of these regions included a cluster of significant SNP on chromosome (Chr) 6 not in linkage disequilibrium to each other. This genomic location contains annotated genes involved in cytokine signalling, haemostasis and mucus biosynthesis. Only one association detected on Chr 7 was significant for both AVPCV and AVLWT. The results generated here reveal candidate immune variants for genes involved in differential response to infection and provide additional SNP marker information that has potential to aid selection of resistance to gastrointestinal parasites in sheep of a similar genetic background to the double backcross population.

Exploring the abomasal lymph node transcriptome for genes associated with resistance to the sheep nematode Teladorsagia circumcincta

This study exploited Blackface lambs that varied in their resistance to the abomasal nematode parasite, Teladorsagia circumcincta. Infection of these lambs over 3 months identified susceptible (high adult worm count, high faecal egg count and low IgA antibody) and resistant animals that had excluded all parasites. Previous work had shown that susceptibility and resistance is dependent on the differential immune response to the parasite, which occurs within the abomasal (gastric) lymph node (ALN) that drains the site of infection. The Affymetrix ovine gene array was used to interrogate the transcriptome of the ALN to identify genes and physiological pathways associated with resistance. We used a bovine RT-qPCR array of 84 genes to validate the gene array, and also report digital gene expression analysis on the same tissues, reanalysed using the Oar v3.1 sheep genome assembly. These analyses identified Humoral Immune Response, Protein Synthesis, Inflammatory Response and Hematological System Development and Function as the two top-ranked networks associated with resistance. Central genes within these networks were IL4, IL5, IL13RA2 and in particular IL13, which confirmed that differential activation of Th2 polarized responses is critical to the resistance phenotype. Furthermore, in resistant sheep there was up-regulation of genes linked to control and suppression of inflammation. The identity of differentially-expressed chemokines and receptors in the resistant and susceptible sheep also begins to explain the cellular nature of the host response to infection. This work will greatly help in the identification of candidate genes as potential selectable markers of genetic resistance.

Use of the Affymetrix Human GeneChip array and genomic DNA hybridisation probe selection to study ovine transcriptomes

Affymetrix GeneChip microarrays are a powerful tool to study global gene expression profiles and have been used on many species. However, no sheep-specific Affymetrix GeneChip is currently available and the bovine array is fairly limited in its coverage and annotation. Previously, a probe-selection method based on hybridisation of genomic DNA (gDNA) was developed, which enables GeneChips to be used for species that they were not designed for. This approach can greatly increase the number of potential annotated genes that can be studied beyond that which is currently available, particularly for species that do not have comprehensively characterised genomes. In this study, we demonstrate that gDNA-based probe selection on the Affymetrix Human U133+2 GeneChip array can be used to study gene expression profiles in sheep tissues. More than 20 000 transcripts were detected in triplicate ovine skeletal muscle and liver samples, which is more than would be possible using the commercially available sheep-specific microarray. The majority of the top 15 differentially expressed genes for each tissue were known to either be expressed in a tissue-specific manner or relate to specific functions of that tissue. Gene ontology analysis of the differentially expressed genes revealed the expected differences in gene expression profiles between the two tissues. Besides demonstrating that the probe selection method can be used to study the ovine transcriptome, the benefits of this approach are that it can greatly increase the number of annotated and novel genes that can be studied beyond those currently possible using ovine- or bovine-specific microarrays. This same method also has the potential to allow the study of other species where species-specific microarrays are not available or whose genomes have not been comprehensively characterised.

Small ruminant resistance against gastrointestinal nematodes: a case of Haemonchus contortus

Gastrointestinal nematode (GIN) infections are a common constraint to small ruminant industry throughout the world, and among those, haemonchosis has its own significance. Control of GIN primarily relies on the use of anthelmintics, but this approach has become less reliable due to the development of resistance in GINs against commonly used anthelmintics and an increased consumer demand for environmentally friendly animal products. These issues have stimulated investigations to find alternative sustainable control strategies, which are less reliant on anthelmintic input. One of such strategies is breeding of small ruminants for their resistance to the GINs. The susceptibility and resistance of animals to GIN infections varies within and between breeds. Various parasitological, biochemical and immunological parameters are employed to evaluate natural resistance status of animals both in natural pasture and artificial infections. The immune mechanisms responsible for resistance are not completely understood, but it has a significant effect in inherited resistance. Relatively resistant or tolerant animals show better local and generalised immune response as compared to susceptible. Immune response against GINs is influenced by many physiological factors. Determination of specific genes linked with host resistance will provide a valuable approach to find out the molecular mechanism of host resistance to GINs. Resistance has been reported to reduce pasture contamination, which in turn reduces re-infection and thus the requirement of the frequent anthelmintic treatments. The efficiency of control can be increased through objective and accurate identification of genetically tolerant individuals by natural and artificial infections with GINs. Complete resistance is the ultimate solution, but this has generally been ignored as a commercial reality. This paper reviews the published reports on natural resistance in small ruminants and discusses the prospects of developing small ruminants, which could be resistant to GINs.

The gastrointestinal nematode Trichostrongylus colubriformis down-regulates immune gene expression in migratory cells in afferent lymph

Background

Gastrointestinal nematode (GIN) infections are the predominant cause of economic losses in sheep. Infections are controlled almost exclusively by the use of anthelmintics which has lead to the selection of drug resistant nematode strains. An alternative control approach would be the induction of protective immunity to these parasites. This study exploits an ovine microarray biased towards immune genes, an artificially induced immunity model and the use of pseudo-afferent lymphatic cannulation to sample immune cells draining from the intestine, to investigate possible mechanisms involved in the development of immunity.

Results

During the development of immunity to, and a subsequent challenge infection with Trichostrongylus colubriformis, the transcript levels of 2603 genes of cells trafficking in afferent intestinal lymph were significantly modulated (P < 0.05). Of these, 188 genes were modulated more than 1.3-fold and involved in immune function. Overall, there was a clear trend for down-regulation of many genes involved in immune functions including antigen presentation, caveolar-mediated endocytosis and protein ubiquitination. The transcript levels of TNF receptor associated factor 5 (TRAF5), hemopexin (HPX), cysteine dioxygenase (CDO1), the major histocompatability complex Class II protein (HLA-DMA), interleukin-18 binding protein (IL-18BP), ephrin A1 (EFNA1) and selenoprotein S (SELS) were modulated to the greatest degree.

Conclusions

This report describes gene expression profiles of afferent lymph cells in sheep developing immunity to nematode infection. Results presented show a global down-regulation of the expression of immune genes which may be reflective of the natural temporal response to nematode infections in livestock.

Genetic variation in resistance to mixed, predominantly Teladorsagia circumcincta nematode infections of sheep: from heritabilities to gene identification

In cool temperate areas, such as Scotland, sheep are infected by a variety of nematodes but the dominant nematode is Teladorsagia circumcincta. Resistant animals have one or more of the following features: fewer adult nematodes, more inhibited larvae, shorter adult nematodes and decreased production of nematode eggs. In lambs at the end of the first grazing season, the heritability of adult worm length is very strong, whereas the heritability of egg production is moderate. The heritability of worm number is low while there is no detectable genetic variation in the number of inhibited larvae. The major mechanisms underlying resistance to T. circumcincta appear to be the IgA mediated suppression of worm growth and the mast cell mediated regulation of worm number. Mast cell responses are slow to develop, possibly because they are responsible for protein loss and reduced growth of the host. Two genes have been repeatedly associated with resistance to T. Circumcincta: the MHC class II DRB1 locus on chromosome 20 and the interferon-gamma locus on chromosome 3. Although the causative mutations are still unknown both genes are plausible candidates.

Immunological profiles of Bos taurus and Bos indicus cattle infested with the cattle tick, Rhipicephalus (Boophilus) microplus

The cattle tick, Rhipicephalus (Boophilus) microplus, is a major threat to the improvement of cattle production in tropical and subtropical countries worldwide. Bos indicus cattle are naturally more resistant to infestation with the cattle tick than are Bos taurus breeds, although considerable variation in resistance occurs within and between breeds. It is not known which genes contribute to the resistant phenotype, nor have immune parameters involved in resistance to R. microplus been fully described for the bovine host. This study was undertaken to determine whether selected cellular and antibody parameters of the peripheral circulation differed between tick-resistant Bos indicus and tick-susceptible Bos taurus cattle following a period of tick infestations. This study demonstrated significant differences between the two breeds with respect to the percentage of cellular subsets comprising the peripheral blood mononuclear cell population, cytokine expression by peripheral blood leukocytes, and levels of tick-specific immunoglobulin G1 (IgG1) antibodies measured in the peripheral circulation. In addition to these parameters, the Affymetrix bovine genome microarray was used to analyze gene expression by peripheral blood leukocytes of these animals. The results demonstrate that the Bos indicus cattle developed a stabilized T-cell-mediated response to tick infestation evidenced by their cellular profile and leukocyte cytokine spectrum. The Bos taurus cattle demonstrated cellular and gene expression profiles consistent with a sustained innate, inflammatory response to infestation, although high tick-specific IgG1 titers suggest that these animals have also developed a T-cell response to infestation.

Microarrays and high-throughput transcriptomic analysis in species with incomplete availability of genomic sequences

Microarrays produce a measurement of gene expression based on the relative measures of dye intensities that correspond to the amount of target RNA. This technology is fast developing and its application is expanding from Homo sapiens to a wide number of species, where enough information on sequences and annotations exist. Anyway, the number of species for which a dedicated platform exists is not high. The use of heterologous array hybridization, screening for gene expression in one species using an array developed for another one, is still quite frequent, even though cross-species microarray hybridization has raised many arguments. Some methods which are high throughput and do not rely on knowledge of the DNA/RNA sequence exist, namely serial analysis of gene expression (SAGE), Massively Parallel Signature Sequencing (MPSS) and deep sequencing of full transcriptome. Although very powerful, particularly the latter, they are still quite costly and cumbersome methods. In some species where genome sequences are largely unknown, several anonymous sequences are deposited in gene banks as a result of Expressed Sequence Tags (ESTs) sequencing projects. The ESTs databases represent a valuable knowledge that can be exploited with some bioinformatic effort to build species-specific microarrays. We present here a method of high-density in situ synthesized microarrays starting from available EST sequences in, Ovis aries. Our data indicate that the method is very efficient and can be easily extended to other species of which genetic sequences are present in public databases, but neglected so far with advanced devices like microarrays. As a perspective, the approach can be applied also to species of which no sequences are available to date, thanks to high-throughput deep sequencing methods.

Sequential microarray to identify timing of molecular responses to Haemonchus contortus infection in sheep

Anthelmintics are currently the most common method of worm control. The emergence of worms with multiple-drug resistance and issues of residues in the food chain make alternative parasite control measures a priority. To develop improved and sustainable methods for controlling Haemonchus contortus such as genetic selection of resistant sheep, a better understanding of the host-parasite relationship is required. A trial was undertaken using sheep surgically implanted with abomasal fistulas to enable sequential biopsy of the abomasal mucosa during trickle infection with two strains of H. contortus. These were ivermectin-resistant CAVR and ivermectin-sensitive McMaster. From a gross parasitology perspective, this approach enabled the effect of developing immunity to be observed on both the establishment and maturation of two CAVR doses within and between groups. Since the only difference in parasite treatment between the groups was the staggering of the two CAVR doses, microarray results from biopsies taken on the same day in different groups were combined and compared between different biopsy dates to observe differential gene transcription over time. Differential gene transcription was detected by comparing transcription in our array data between different biopsy dates using a low P value screen (P<0.01) and by compiling a list of 82 immunoparasitology-related genes and examining transcription in this list with a higher P value screen (P<0.05). Our microarray data were validated in silico by comparison with intelectin 2, trefoil factor 3, calcium activated chloride channel and mucin 5 from other gene transcription studies and with phenotypic data such as the response by gammadelta T cells and immunoglobulins to H. contortus. The first four genes are involved in non-specific responses to infection and mucosal healing. These were upregulated at the early time points and intelectin 2 remained prominent throughout the trial. As the trial progressed, immunoglobulin genes became strongly upregulated. These included IgCgamma IgG2a heavy chain constant region, IGHE immunoglobulin heavy constant epsilon and IGHM immunoglobulin heavy constant mu.