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

Multi-omics data elucidate parasite-host-microbiota interactions and resistance to Haemonchus contortus in sheep

BACKGROUND

The integration of molecular data from hosts, parasites, and microbiota can enhance our understanding of the complex biological interactions underlying the resistance of hosts to parasites. Haemonchus contortus, the predominant sheep gastrointestinal parasite species in the tropics, causes significant production and economic losses, which are further compounded by the diminishing efficiency of chemical control owing to anthelmintic resistance. Knowledge of how the host responds to infection and how the parasite, in combination with microbiota, modulates host immunity can guide selection decisions to breed animals with improved parasite resistance. This understanding will help refine management practices and advance the development of new therapeutics for long-term helminth control.

METHODS

Eggs per gram (EPG) of feces were obtained from Morada Nova sheep subjected to two artificial infections with H. contortus and used as a proxy to select animals with high resistance or susceptibility for transcriptome sequencing (RNA-seq) of the abomasum and 50 K single-nucleotide genotyping. Additionally, RNA-seq data for H. contortus were generated, and amplicon sequence variants (ASV) were obtained using polymerase chain reaction amplification and sequencing of bacterial and archaeal 16S ribosomal RNA genes from sheep feces and rumen content.

RESULTS

The heritability estimate for EPG was 0.12. GAST, GNLY, IL13, MGRN1, FGF14, and RORC genes and transcripts were differentially expressed between resistant and susceptible animals. A genome-wide association study identified regions on chromosomes 2 and 11 that harbor candidate genes for resistance, immune response, body weight, and adaptation. Trans-expression quantitative trait loci were found between significant variants and differentially expressed transcripts. Functional co-expression modules based on sheep genes and ASVs correlated with resistance to H. contortus, showing enrichment in pathways of response to bacteria, immune and inflammatory responses, and hub features of the Christensenellaceae, Bacteroides, and Methanobrevibacter genera; Prevotellaceae family; and Verrucomicrobiota phylum. In H. contortus, some mitochondrial, collagen-, and cuticle-related genes were expressed only in parasites isolated from susceptible sheep.

CONCLUSIONS

The present study identified chromosome regions, genes, transcripts, and pathways involved in the elaborate interactions between the sheep host, its gastrointestinal microbiota, and the H. contortus parasite. These findings will assist in the development of animal selection strategies for parasite resistance and interdisciplinary approaches to control H. contortus infection in sheep.

Host gene expression in wildlife disease: making sense of species-level responses

Emerging infectious diseases are significant threats to wildlife conservation, yet the impacts of pathogen exposure and infection can vary widely among host species. As such, conservation biologists and disease ecologists have increasingly aimed to understand species-specific host susceptibility using molecular methods. In particular, comparative gene expression assays have been used to contrast the transcriptomic responses of disease-resistant and disease-susceptible hosts to pathogen exposure. This work usually assumes that the gene expression responses of disease-resistant species will reveal the activation of molecular pathways contributing to host defence. However, results often show that disease-resistant hosts undergo little gene expression change following pathogen challenge. Here, we discuss the mechanistic implications of these "null" findings and offer methodological suggestions for future molecular studies of wildlife disease. First, we highlight that muted transcriptomic responses with minimal immune system recruitment may indeed be protective for nonsusceptible hosts if they limit immunopathology and promote pathogen tolerance in systems where susceptible hosts suffer from genetic dysregulation. Second, we argue that overly narrow investigation of responses to pathogen exposure may overlook important, constitutively active molecular pathways that underlie species-specific defences. Finally, we outline alternative study designs and approaches that complement interspecific transcriptomic comparisons, including intraspecific gene expression studies and genomic methods to detect signatures of selection. Collectively, these insights will help ecologists extract maximal information from conservation-relevant transcriptomic data sets, leading to a deeper understanding of host defences and, ultimately, the implementation of successful conservation interventions.

Transcriptional profiling of the ovine abomasal lymph node reveals a role for timing of the immune response in gastrointestinal nematode resistance

Gastrointestinal nematodes are a serious cause of morbidity and mortality in grazing ruminants. The major ovine defence mechanism is acquired immunity, with protective immunity developing over time in response to infection. Nematode resistance varies both within and between breeds and is moderately heritable. A detailed understanding of the genes and mechanisms involved in protective immunity, and the factors that regulate this response, is required to aid both future breeding strategies and the development of effective and sustainable nematode control methods. The aim of this study was to compare the abomasal lymph node transcriptome of resistant and susceptible lambs in order to determine biological processes differentially expressed between resistant and susceptible individuals. Scottish Blackface lambs, with divergent phenotypes for resistance, were challenged with 30,000 Teladorsagia circumcincta larvae (L3), and abomasal lymph nodes recovered at 7 and 14days post-infection (dpi). High-throughput sequencing of cDNA from the abomasal lymph node was used to quantitatively sample the transcriptome with an average of 32 million reads per sample. A total of 194 and 144 genes were differentially expressed between resistant and susceptible lambs at 7 and 14 dpi respectively. Differentially expressed networks and biological processes were identified using Ingenuity Pathway Analysis. Genes involved in the inflammatory response, attraction of T lymphocytes and binding of leukocytes were more highly expressed in resistant animals at 7 dpi and in susceptible animals at 14 dpi indicating that resistant animals respond to infection earlier than susceptible animals. Twenty-four Single Nucleotide Polymorphisms (SNP) within 11 differentially expressed genes, were tested for association with gastrointestinal nematode resistance in the Scottish Blackface lambs. Four SNP, in 2 genes (SLC30A2 and ALB), were suggestively associated with faecal egg count. In conclusion, a large number of genes were differentially expressed in the abomasal lymph node of resistant and susceptible lambs responding to gastrointestinal nematode challenge. Resistant Scottish Blackface lambs appear to generate an earlier immune response to T. circumcincta. In susceptible lambs this response appears to be delayed. SNP in 2 differentially expressed genes were suggestively associated with faecal egg count indicating that differentially expressed genes may be considered candidate loci for mediating nematode resistance.

Relative maturity and the development of immunity to gastrointestinal nematodes in sheep: an overlooked paradigm?

Sheep display considerable variation in both the timing and magnitude of development of immunity to gastrointestinal nematodes (GIN). Onset of immunity is dependent on a number of factors, including antigenic stimulus, nutrition supply, age and size of the animals, the latter of which are confounded. Here, we review the factors associated with the development of immunity to GIN in sheep, particularly in the context of the role that relative maturity may have through applying the rules of genetic size scaling based on examples from published literature. Comparing animals based on their metabolic age, rather than chronological age, may provide an explanation for the timing of immune development and may reduce the variation in immune development that frequently is observed both between and within breeds. Further, this approach may help explain the phenotypic differences in animal performance between animals of varying immunological capacity to GIN through influences on mature body weight. As such, when considering factors influencing immune development to GIN, physiological age or relative maturity may be considered an overlooked paradigm. We propose it may be worthwhile to consider metabolic age when comparing the immune competence of animals to ensure the subjects are at an analogous stage of physiological development.

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.

Genome-wide association and genome partitioning reveal novel genomic regions underlying variation in gastrointestinal nematode burden in a wild bird

Identifying the genetic architecture underlying complex phenotypes is a notoriously difficult problem that often impedes progress in understanding adaptive eco-evolutionary processes in natural populations. Host-parasite interactions are fundamentally important drivers of evolutionary processes, but a lack of understanding of the genes involved in the host's response to chronic parasite insult makes it particularly difficult to understand the mechanisms of host life history trade-offs and the adaptive dynamics involved. Here, we examine the genetic basis of gastrointestinal nematode (Trichostrongylus tenuis) burden in 695 red grouse (Lagopus lagopus scotica) individuals genotyped at 384 genome-wide SNPs. We first use genome-wide association to identify individual SNPs associated with nematode burden. We then partition genome-wide heritability to identify chromosomes with greater heritability than expected from gene content, due to harbouring a multitude of additive SNPs with individually undetectable effects. We identified five SNPs on five chromosomes that accounted for differences of up to 556 worms per bird, but together explained at best 4.9% of the phenotypic variance. These SNPs were closely linked to genes representing a range of physiological processes including the immune system, protein degradation and energy metabolism. Genome partitioning indicated genome-wide heritability of up to 29% and three chromosomes with excess heritability of up to 4.3% (total 8.9%). These results implicate SNPs and novel genomic regions underlying nematode burden in this system and suggest that this phenotype is somewhere between being based on few large-effect genes (oligogenic) and based on a large number of genes with small individual but large combined effects (polygenic).

Variation in the Ovine Abomasal Lymph Node Transcriptome between Breeds Known to Differ in Resistance to the Gastrointestinal Nematode

Texel lambs are known to be more resistant to gastrointestinal nematode (GIN) infection than Suffolk lambs, with a greater ability to limit infection. The objectives of this study were to: 1) profile the whole transcriptome of abomasal lymph node tissue of GIN-free Texel and Suffolk lambs; 2) identify differentially expressed genes and characterize the immune-related biological pathways and networks associated with these genes. Abomasal lymph nodes were collected from Texel (n = 6) and Suffolk (n = 4) lambs aged 19 weeks that had been GIN-free since 6 weeks of age. Whole transcriptome profiling was performed using RNA-seq on the Illumina platform. At the time of conducting this study, a well annotated Ovine genome was not available and hence the sequence reads were aligned with the Bovine (UMD3.1) genome. Identification of differentially expressed genes was followed by pathway and network analysis. The Suffolk breed accounted for significantly more of the differentially expressed genes, (276 more highly expressed in Suffolk v 162 in Texel; P < 0.001). The four most significant differentially expressed pathways were all related to immunity and were classified as: Role of Pattern Recognition Receptors in Recognition of Bacteria and Viruses, Activation of IRF by Cytosolic Pattern Recognition Receptors, Role of RIG-I-like Receptors in Antiviral Innate Immunity, and Interferon Signaling. Of significance is the fact that all of these four pathways were more highly expressed in the Suffolk. These data suggest that in a GIN-free environment, Suffolk lambs have a more active immune profile relative to the Texel: this immune profile may contribute to the poorer efficiency of response to a GIN challenge in the Suffolk breed compared to the Texel breed.

Profile of whole blood gene expression following immune stimulation in a wild passerine

Background

Immunoecology aims to explain variation among hosts in the strength and efficacy of immunological defences in natural populations. This requires development of biomarkers of the activation of the immune system so that they can be collected non-lethally and sampled from small amounts of easily obtainable tissue. We used transcriptome profiling in wild greenfinches (Carduelis chloris) to detect whole blood transcripts that most profoundly indicate upregulation of antimicrobial defences during acute phase response. The more general aim of this study was to obtain a functional annotation of a substantial portion of the greenfinch transcriptome that would enable to gain access to more specific genomic tools in subsequent studies. The birds received either bacterial lipopolysaccharide or saline injections and RNA-seq transcriptional profiling was performed 12 h after treatment to provide initial functional annotation of the transcriptome and assess whole blood response to immune stimulation.

Results

A total of 66,084 transcripts were obtained from de novo Trinty assembly, out of which 23,153 could be functionally annotated. Only 1,911 of these were significantly upregulated or downregulated. The manipulation caused marked upregulation of several transcripts related to immune activation. These included avian-specific antimicrobial agents avidin and gallinacin, but also some more general host response genes, such as serum amyloid A protein, lymphocyte antigen 75 and copper-transporting ATPase 1. However, links with avian immunity for most differentially regulated transcripts remained rather hypothetical, as a large set of differentially expressed transcripts lacked functional annotation.

Conclusions

This appears to be the first large scale transcriptional profiling of immune function in passerine birds. The transcriptomic data obtained suggest novel markers for the assessment of the immunological state of wild passerines. Characterizing the function of those possible novel infection markers would assist future vertebrate genome annotation. The extensive sequence information collected enables to identify possible target and housekeeping genes needed to gain access to more specific genomic tools in future studies.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-533) contains supplementary material, which is available to authorized users.

Review: Genetics of helminth resistance in sheep

Karrow, N. A., Goliboski, K., Stonos, N., Schenkel, F. and Peregrine, A. 2014. Review: Genetics of helminth resistance in sheep. Can. J. Anim. Sci. 94: 1–9. Gastrointestinal helminth parasites are an important source of economic loss to sheep producers. A rapid increase in anthelmintic resistance has occurred around the globe; therefore, the industry is exploring alternative strategies such as genetic selection to control losses attributed to helminth infection. Since helminths have co-evolved with sheep for millions of years, natural selection for enhanced helminth resistance has occurred within certain breeds from various parts of the world. These breeds of sheep are being used to better understand the genetic aspects of helminth resistance. If the genetic variants that contribute to this phenotype can be identified, it may be possible to use selection strategies to introduce resistance alleles into other breeds or to increase their frequency within breeds. This review will provide an up-to-date overview of the pathology of helminth disease, the immune response to helminth infection, and the search for genes that confer helminth resistance.

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.

The immunology and genetics of resistance of sheep to Teladorsagia circumcincta

Teladorsagia circumcincta is one of the most economically important gastrointestinal nematode parasites of sheep in cool temperate regions, to which sheep show genetically-varying resistance to infection. This is a very common parasite and viable sheep production requires the extensive use of anthelmintic drugs. However, the emergence of drug-resistant parasites has stimulated the search for alternative control strategies to curb production losses. Lambs become infected soon after weaning and begin to control parasite burden within 8-10 weeks of continual infection. This control is an acquired characteristic mediated by the development of parasite-specific antibodies. This paper describes the immunology associated with resistance and susceptibility, focussing on differential T cell activation that regulates the production of specific effector mechanisms. It continues by summarizing the methods used to identify genes that could be exploited as molecular markers of selection for resistance. In particular it focusses on the link between understanding the molecular immunology of infection and the identification of candidate genes for selection.

Detecting genes for variation in parasite burden and immunological traits in a wild population: testing the candidate gene approach

Identifying the genes underlying phenotypic variation in natural populations can provide novel insight into the evolutionary process. The candidate gene approach has been applied to studies of a number of traits in various species, in an attempt to elucidate their genetic basis. Here, we test the application of the candidate gene approach to identify the loci involved in variation in gastrointestinal parasite burden, a complex trait likely to be controlled by many loci, in a wild population of Soay sheep. A comprehensive literature review, Gene Ontology databases, and comparative genomics resources between cattle and sheep were used to generate a list of candidate genes. In a pilot study, these candidates, along with 50 random genes, were then sequenced in two pools of Soay sheep; one with low gastrointestinal nematode burden and the other high, using a NimbleGen sequence capture experiment. Further candidates were identified from single nucleotide polymorphisms (SNPs) that were highly differentiated between high- and low-resistance sheep breeds. A panel of 192 candidate and control SNPs were then typed in 960 individual Soay sheep to examine whether they individually explained variation in parasite burden, as measured as faecal egg count, as well as two immune measures (Teladorsagia circumcincta-specific antibodies and antinuclear antibodies). The cumulative effect of the candidate and control SNPs were estimated by fitting genetic relationship matrices (GRMs) as random effects in animal models of the three traits. No more significant SNPs were identified in the pilot sequencing experiment and association study than expected by chance. Furthermore, no significant difference was found between the proportions of candidate or control SNPs that were found to be significantly associated with parasite burden/immune measures. No significant effect of the candidate or control gene GRMs was found. There is thus little support for the candidate gene approach to the identification of loci explaining variation in parasitological and immunological traits in this population. However, a number of SNPs explained significant variation in multiple traits and significant correlations were found between the proportions of variance explained by individual SNPs across multiple traits. The significant SNPs identified in this study may still, therefore, merit further investigation.

Exploring local immunological adaptation of two stickleback ecotypes by experimental infection and transcriptome-wide digital gene expression analysis

Understanding the extent of local adaptation in natural populations and the mechanisms that allow individuals to adapt to their native environment is a major avenue in molecular ecology research. Evidence for the frequent occurrence of diverging ecotypes in species that inhabit multiple ecological habitats is accumulating, but experimental approaches to understanding the biological pathways as well as the underlying genetic mechanisms are still rare. Parasites are invoked as one of the major selective forces driving evolution and are themselves dependent on the ecological conditions in a given habitat. Immunological adaptation to local parasite communities is therefore expected to be a key component of local adaptation in natural populations. Here, we use next-generation sequencing technology to compare the transcriptome-wide response of experimentally infected three-spined sticklebacks from a lake and a river population, which are known to evolve under selection by distinct parasite communities. By comparing overall gene expression levels as well as the activation of functional pathways in response to parasite exposure, we identified potential differences between the two stickleback populations at several levels. Our results suggest locally adapted patterns of gene regulation in response to parasite exposure, which may reflect different local optima in the trade-off between the benefits and the disadvantages of mounting an immune response because of quantitative differences of the local parasite communities.

Flowering time and transcriptome variation in Capsella bursa-pastoris (Brassicaceae)

• Flowering is a major developmental transition and its timing in relation to environmental conditions is of crucial importance to plant fitness. Understanding the genetic basis of flowering time variation is important to determining how plants adapt locally. • Here, we investigated flowering time variation of Capsella bursa-pastoris collected from different latitudes in China. We also used a digital gene expression (DGE) system to generate partial gene expression profiles for 12 selected samples. • We found that flowering time was highly variable and most strongly correlated with day length and winter temperature. Significant differences in gene expression between early- and late-flowering samples were detected for 72 candidate genes for flowering time. Genes related to circadian rhythms were significantly overrepresented among the differentially expressed genes. • Our data suggest that circadian rhythms and circadian clock genes play an important role in the evolution of flowering time, and C. bursa-pastoris plants exhibit expression differences for candidate genes likely to affect flowering time across the broad range of environments they face in China.

Relationship between susceptibility of Blackface sheep to Teladorsagia circumcincta infection and an inflammatory mucosal T cell response

Teladorsagia circumcincta is the most economically important gastrointestinal (abomasal) nematode parasite of sheep in cool temperate regions, to which sheep show genetically-varying resistance to infection. Lambs, from parents with genetic variation for resistance, were trickle infected with L3 larvae over 12 weeks. 45 lambs were identified with a range of susceptibilities as assessed by: adult worm count at post mortem, faecal egg count (FEC) and IgA antibody levels. This project investigated the correlation of T cell cytokine expression and resistance to infection at the mature stage of response, when the resistant lambs had excluded all parasites.

Histopathology showed only minor changes in resistant animals with a low level lymphocyte infiltration; but in susceptible lambs, major pathological changes were associated with extensive infiltration of lymphocytes, eosinophils and neutrophils.

Absolute quantitative RT-qPCR assays on the abomasal lymph node (ALN) revealed a significant positive correlation between IL6, IL21 and IL23A transcript levels with adult worm count and FEC. IL23A was also negatively correlated with IgA antibody levels. Significantly positive correlation of TGFB1 levels with adult worm count and FEC were also seen in the abomasal mucosa. These data are consistent with the hypothesis that the inability to control L3 larval colonization, adult worm infection and egg production is due to the activation of the inflammatory Th17 T cell subset.

Snapshot of spatio-temporal cytokine responses to single and co-infections with helminths and bacteria

Cytokines play a key role in maintaining communication between organs and in so doing modulate the interaction between concurrent infections. The extent of these effects depends on the properties of the organ infected and the intensity and type of infections. To determine systemic bystander effects among organs, IFN-γ, IL-4 and IL-10 gene expression was quantified at 7 days post-challenge in directly infected and uninfected organs during single and co-infections with the respiratory bacterium Bordetella bronchiseptica and the gastrointestinal helminths Graphidium strigosum and Trichostrongylus retortaeformis. Results showed that cytokine expression in a specific organ was influenced by the type of infection occurring in another organ, and this bystander effect was more apparent in some organs than others. Within the same organ the relative cytokine expression was consistent across infections, although some cytokines were more affected by bystander effects than others. For the infected gastrointestinal tract, a stronger cytokine response was observed in the tissue that harbored the majority of helminths (i.e. duodenum and fundus). Overall, co-infections altered the intensity but to a lesser extent the relative cytokine profile against the focal infection, indicating clear bystander effects and low organ compartmentalization. However, organs appear to actively modulate cytokine expression to avoid potential immuno-pathological consequences.

Identification of genes associated with nitrogen-use efficiency by genome-wide transcriptional analysis of two soybean genotypes

BACKGROUND

Soybean is a valuable crop that provides protein and oil. Soybean requires a large amount of nitrogen (N) to accumulate high levels of N in the seed. The yield and protein content of soybean seeds are directly affected by the N-use efficiency (NUE) of the plant, and improvements in NUE will improve yields and quality of soybean products. Genetic engineering is one of the approaches to improve NUE, but at present, it is hampered by the lack of information on genes associated with NUE. Solexa sequencing is a new method for estimating gene expression in the transcription level. Here, the expression profiles were analyzed between two soybean varieties in N-limited conditions to identify genes related to NUE.

RESULTS

Two soybean genotypes were grown under N-limited conditions; a low-N-tolerant variety (No.116) and a low-N-sensitive variety (No.84-70). The shoots and roots of soybeans were used for sequencing. Eight libraries were generated for analysis: 2 genotypes × 2 tissues (roots and shoots) × 2 time periods [short-term (0.5 to 12 h) and long-term (3 to 12 d) responses] and compared the transcriptomes by high-throughput tag-sequencing analysis. 5,739,999, 5,846,807, 5,731,901, 5,970,775, 5,476,878, 5,900,343, 5,930,716, and 5,862,642 clean tags were obtained for the eight libraries: L1, 116-shoot short-term; L2 84-70-shoot short-term; L3 116-shoot long-term; L4 84-70-shoot long-term; L5 116-root short-term; L6 84-70-root short-term; L7 116-root long-term;L8 84-70-root long-term; these corresponded to 224,154, 162,415, 191,994, 181,792, 204,639, 206,998, 233,839 and 257,077 distinct tags, respectively. The clean tags were mapped to the reference sequences for annotation of expressed genes. Many genes showed substantial differences in expression among the libraries. In total, 3,231 genes involved in twenty-two metabolic and signal transduction pathways were up- or down-regulated. Twenty-four genes were randomly selected and confirmed their expression patterns by quantitative RT-PCR; Twenty-one of the twenty-four genes showed expression patterns consistent with the Digital Gene Expression (DGE) data.

CONCLUSIONS

A number of soybean genes were differentially expressed between the low-N-tolerant and low-N-sensitive varieties under N-limited conditions. Some of these genes may be candidates for improving NUE. These findings will help to provide a detailed understanding of NUE mechanisms, and also provide a basis for breeding soybean varieties that are tolerant to low-N conditions.