Assessing rates of parasite coinfection and spatiotemporal strain variation via metabarcoding: Insights for the conservation of European turtle doves Streptopelia turtur.

Understanding the frequency, spatiotemporal dynamics and impacts of parasite coinfections is fundamental to developing control measures and predicting disease impacts. The European turtle dove (Streptopelia turtur) is one of Europe's most threatened bird species. High prevalence of infection by the protozoan parasite Trichomonas gallinae has previously been identified, but the role of this and other coinfecting parasites in turtle dove declines remains unclear. Using a high-throughput sequencing approach, we identified seven strains of T. gallinae, including two novel strains, from ITS1/5.8S/ITS2 ribosomal sequences in turtle doves on breeding and wintering grounds, with further intrastrain variation and four novel subtypes revealed by the iron-hydrogenase gene. High spatiotemporal turnover was observed in T. gallinae strain composition, and infection was prevalent in all populations (89%-100%). Coinfection by multiple Trichomonas strains was rarer than expected (1% observed compared to 38.6% expected), suggesting either within-host competition, or high mortality of coinfected individuals. In contrast, coinfection by multiple haemosporidians was common (43%), as was coinfection by haemosporidians and T. gallinae (90%), with positive associations between strains of T. gallinae and Leucocytozoon suggesting a mechanism such as parasite-induced immune modulation. We found no evidence for negative associations between coinfections and host body condition. We suggest that longitudinal studies involving the recapture and investigation of infection status of individuals over their lifespan are crucial to understand the epidemiology of coinfections in natural populations.

Development of a multiplex microsatellite marker set for the study of the solitary red mason bee, Osmia bicornis (Megachilidae).

BACKGROUND: Solitary bees, such as the red mason bee (Osmia bicornis), provide important ecosystem services including pollination. In the face of global declines of pollinator abundance, such haplodiploid Hymenopterans have a compounded extinction risk due to the potential for limited genetic diversity. In order to assess the genetic diversity of Osmia bicornis populations, we developed microsatellite markers and characterised them in two populations. METHODS AND RESULTS: Microsatellite sequences were mined from the recently published Osmia bicornis genome, which was assembled from DNA extracted from a single male bee originating from the United Kingdom. Sequences were identified that contained dinucleotide, trinucleotide, and tetranucleotide repeat regions. Seventeen polymorphic microsatellite markers were designed and tested, sixteen of which were developed into four multiplex PCR sets to facilitate cheap, fast and efficient genotyping and were characterised in unrelated females from Germany (n = 19) and England (n = 14). CONCLUSIONS: The microsatellite markers are highly informative, with a combined exclusion probability of 0.997 (first parent), which will enable studies of genetic structure and diversity to inform conservation efforts in this bee.

Subspecies hybridization as a potential conservation tool in species reintroductions.

Reintroductions are a powerful tool for the recovery of endangered species. However, their long-term success is strongly influenced by the genetic diversity of the reintroduced population. The chances of population persistence can be improved by enhancing the population's adaptive ability through the mixing of individuals from different sources. However, where source populations are too diverse the reintroduced population could also suffer from outbreeding depression or unsuccessful admixture due to behavioural or genetic barriers. For the reintroduction of Asiatic wild ass Equus hemionus ssp. in Israel, a breeding core was created from individuals of two different subspecies (E. h. onager & E. h. kulan). Today the population comprises approximately 300 individuals and displays no signs of outbreeding depression. The aim of this study was a population genomic evaluation of this conservation reintroduction protocol. We used maximum likelihood methods and genetic clustering analyses to investigate subspecies admixture and test for spatial autocorrelation based on subspecies ancestry. Further, we analysed heterozygosity and effective population sizes in the breeding core prior to release and the current wild population. We discovered high levels of subspecies admixture in the breeding core and wild population, consistent with a significant heterozygote excess in the breeding core. Furthermore, we found no signs of spatial autocorrelation associated with subspecies ancestry in the wild population. Inbreeding and variance effective population size estimates were low. Our results indicate no genetic or behavioural barriers to admixture between the subspecies and suggest that their hybridization has led to greater genetic diversity in the reintroduced population. The study provides rare empirical evidence of the successful application of subspecies hybridization in a reintroduction. It supports use of intraspecific hybridization as a tool to increase genetic diversity in conservation translocations.

Microsatellite characterisation and sex-typing in two invasive parakeet species, the monk parakeet Myiopsitta monachus and ring-necked parakeet Psittacula krameri.

Invasive species can have wide-ranging negative impacts, and an understanding of the process and success of invasions can be vital to determine management strategies, mitigate impacts and predict range expansions of such species. Monk parakeets (Myiopsitta monachus) and ring-necked parakeets (Psittacula krameri) are both widespread invasive species, but there has been little research into the genetic and social structure of these two species despite the potential links with invasion success. The aim of this study was to isolate novel microsatellite loci from the monk parakeet and characterise them in both monk and ring-necked parakeets in order to facilitate future investigations into their behaviour and population ecology. Sex-typing markers were also tested in both species. Of the 20 microsatellite loci assessed in 24 unrelated monk parakeets, 16 successfully amplified and were polymorphic displaying between 2 and 14 alleles (mean = 8.06). Expected heterozygosity ranged from 0.43 to 0.93 and observed heterozygosity ranged from 0.23 to 0.96. Nine of the 20 loci also successfully amplified and were polymorphic in the ring-necked parakeet, displaying between 2 and 10 alleles. Suitable markers to sex both species and a Z-linked microsatellite locus were identified. A multiplex marker set was validated for monk parakeets. These novel microsatellite loci will facilitate fine and broad-scale population genetic analyses of these two widespread invasive species.

Genomic associations with bill length and disease reveal drift and selection across island bird populations.

Island species provide excellent models for investigating how selection and drift operate in wild populations, and for determining how these processes act to influence local adaptation and speciation. Here, we examine the role of selection and drift in shaping genomic and phenotypic variation across recently separated populations of Berthelot's pipit (Anthus berthelotii), a passerine bird endemic to three archipelagos in the Atlantic. We first characterized genetic diversity and population structuring that supported previous inferences of a history of recent colonizations and bottlenecks. We then tested for regions of the genome associated with the ecologically important traits of bill length and malaria infection, both of which vary substantially across populations in this species. We identified a SNP associated with variation in bill length among individuals, islands, and archipelagos; patterns of variation at this SNP suggest that both phenotypic and genotypic variation in bill length is largely shaped by founder effects. Malaria was associated with SNPs near/within genes involved in the immune response, but this relationship was not consistent among archipelagos, supporting the view that disease resistance is complex and rapidly evolving. Although we found little evidence for divergent selection at candidate loci for bill length and malaria resistance, genome scan analyses pointed to several genes related to immunity and metabolism as having important roles in divergence and adaptation. Our findings highlight the utility and challenges involved with combining association mapping and population genetic analysis in nonequilibrium populations, to disentangle the effects of drift and selection on shaping genotypes and phenotypes.

Correction to: A multiplex marker set for microsatellite typing and sexing of sooty terns Onychoprion fuscatus.

Following publication of the original article [1], one of the authors reported that his name was listed incorrectly, and that he would like his name to appear as S. James Reynolds instead of Silas James Reynolds. The latter format would confuse citations as all his previous publications are in the former format.

Similar yet different: co-analysis of the genetic diversity and structure of an invasive nematode parasite and its invasive mammalian host.

Animal parasitic nematodes can cause serious diseases and their emergence in new areas can be an issue of major concern for biodiversity conservation and human health. Their ability to adapt to new environments and hosts is likely to be affected by their degree of genetic diversity, with gene flow between distinct populations counteracting genetic drift and increasing effective population size. The raccoon roundworm (Baylisascaris procyonis), a gastrointestinal parasite of the raccoon (Procyon lotor), has increased its global geographic range after being translocated with its host. The raccoon has been introduced multiple times to Germany, but not all its populations are infected with the parasite. While fewer introduced individuals may have led to reduced diversity in the parasite, admixture between different founder populations may have counteracted genetic drift and bottlenecks. Here, we analyse the population genetic structure of the roundworm and its raccoon host at the intersection of distinct raccoon populations infected with B. procyonis. We found evidence for two parasite clusters resulting from independent introductions. Both clusters exhibited an extremely low genetic diversity, suggesting small founding populations subjected to inbreeding and genetic drift with no, or very limited, genetic influx from population admixture. Comparison of the population genetic structures of both host and parasite suggested that the parasite spread to an uninfected raccoon founder population. On the other hand, an almost perfect match between cluster boundaries also suggested that the population genetic structure of B. procyonis has remained stable since its introduction, mirroring that of its raccoon host.

Identification of 24 new microsatellite loci in the sweat bee Lasioglossum malachurum (Hymenoptera: Halictidae).

OBJECTIVE: The objective here is to identify highly polymorphic microsatellite loci for the Palaearctic sweat bee Lasioglossum malachurum. Sweat bees (Hymenoptera: Halictidae) are widespread pollinators that exhibit an unusually large range of social behaviours from non-social, where each female nests alone, to eusocial, where a single queen reproduces while the other members of the colony help to rear her offspring. They thus represent excellent models for understanding social evolution. RESULTS: 24 new microsatellite loci were successfully optimized. When amplified across 23-40 unrelated females, the number of alleles per locus ranged from 3 to 17 and the observed heterozygosities 0.45 to 0.95. Only one locus showed evidence of significant deviation from Hardy-Weinberg equilibrium. No evidence of linkage disequilibrium was found. These 24 loci will enable researchers to gain greater understanding of colony relationships within this species, an important model for the study of eusociality. Furthermore, 22 of the same loci were also successfully amplified in L. calceatum, suggesting that these loci may be useful for investigating the ecology and evolution of sweat bees in general.

A multiplex marker set for microsatellite typing and sexing of sooty terns Onychoprion fuscatus.

OBJECTIVES: Seabirds have suffered dramatic population declines in recent decades with one such species being the sooty tern Onychoprion fuscatus. An urgent call to re-assess their conservation status has been made given that some populations, such as the one on Ascension Island, South Atlantic, have declined by over 80% in three generations. Little is known about their population genetics, which would aid conservation management through understanding ecological processes and vulnerability to environmental change. We developed a multiplex microsatellite marker set for sooty terns including sex-typing markers to assist population genetics studies. RESULTS: Fifty microsatellite loci were isolated and tested in 23 individuals from Ascension Island. Thirty-one were polymorphic and displayed between 4 and 20 alleles. Three loci were Z-linked and two autosomal loci deviated from Hardy-Weinberg equilibrium. The remaining 26 autosomal loci together with three sex-typing makers were optimised in seven polymerase chain reaction plexes. These 26 highly polymorphic markers will be useful for understanding genetic structure of the Ascension Island population and the species as a whole. Combining these with recently developed microsatellite markers isolated from Indian Ocean birds will allow for assessment of global population structure and genetic diversity.

Widespread gene flow between oceans in a pelagic seabird species complex.

Global-scale gene flow is an important concern in conservation biology as it has the potential to either increase or decrease genetic diversity in species and populations. Although many studies focus on the gene flow between different populations of a single species, the potential for gene flow and introgression between species is understudied, particularly in seabirds. The only well-studied example of a mixed-species, hybridizing population of petrels exists on Round Island, in the Indian Ocean. Previous research assumed that Round Island represents a point of secondary contact between Atlantic (Pterodroma arminjoniana) and Pacific species (Pterodroma neglecta and Pterodroma heraldica). This study uses microsatellite genotyping and tracking data to address the possibility of between-species hybridization occurring outside the Indian Ocean. Dispersal and gene flow spanning three oceans were demonstrated between the species in this complex. Analysis of migration rates estimated using bayesass revealed unidirectional movement of petrels from the Atlantic and Pacific into the Indian Ocean. Conversely, structure analysis revealed gene flow between species of the Atlantic and Pacific oceans, with potential three-way hybrids occurring outside the Indian Ocean. Additionally, geolocation tracking of Round Island petrels revealed two individuals travelling to the Atlantic and Pacific. These results suggest that interspecific hybrids in Pterodroma petrels are more common than was previously assumed. This study is the first of its kind to investigate gene flow between populations of closely related Procellariiform species on a global scale, demonstrating the need for consideration of widespread migration and hybridization in the conservation of threatened seabirds.

The characterisation of microsatellite markers reveals tetraploidy in the Greater Water Parsnip, Sium latifolium (Apiaceae).

BACKGROUND: The Greater Water Parsnip, Sium latifolium (Apiaceae), is a marginal aquatic perennial currently endangered in England and consequently the focus of a number of conservation translocation projects. Microsatellite markers were developed for S. latifolium to facilitate comparison of genetic diversity and composition between natural and introduced populations. RESULTS: We selected 65 S. latifolium microsatellite (MiSeq) sequences and designed primer pairs for these. Primer sets were tested in 32 individuals. We found 15 polymorphic loci that amplified consistently. For the selected 15 loci, the number of alleles per locus ranged from 8 to 17. For all loci, S. latifolium individuals displayed up to four alleles indicating polyploidy in this species. CONCLUSIONS: These are the first microsatellite loci developed for S. latifolium and each individual displayed 1-4 alleles per locus, suggesting polyploidy in this species. These markers provide a valuable resource in evaluating the population genetic composition of this endangered species and thus will be useful for guiding conservation and future translocations of the species.

Blood transcriptomes and de novo identification of candidate loci for mating success in lekking great snipe (Gallinago media).

We assembled the great snipe blood transcriptome using data from fourteen lekking males, in order to de novo identify candidate genes related to sexual selection, and determined the expression profiles in relation to mating success. The three most highly transcribed genes were encoding different haemoglobin subunits. All tended to be overexpressed in males with high mating success. We also called single nucleotide polymorphisms (SNPs) from the transcriptome data and found considerable genetic variation for many genes expressed during lekking. Among these, we identified 14 polymorphic candidate SNPs that had a significant genotypic association with mating success (number of females mated with) and/or mating status (mated or not). Four of the candidate SNPs were found in HBAA (encoding the haemoglobin α-chain). Heterozygotes for one of these and one SNP in the gene PABPC1 appeared to enjoy higher mating success compared to males homozygous for either of the alleles. In a larger data set of individuals, we genotyped 38 of the identified SNPs but found low support for consistent selection as only one of the zygosities of previously identified candidate SNPs and none of their genotypes were associated with mating status. However, candidate SNPs generally showed lower levels of spatial genetic structure compared to noncandidate markers. We also scored the prevalence of avian malaria in a subsample of birds. Males infected with avian malaria parasites had lower mating success in the year of sampling than noninfected males. Parasite infection and its interaction with specific genes may thus affect performance on the lek.

A new marker based on the avian spindlin gene that is able to sex most birds, including species problematic to sex with CHD markers.

We have developed a new marker (Z43B) that can be successfully used to identify the sex of most birds (69%), including species difficult or impossible to sex with other markers. We utilized the zebra finch Taeniopygia guttata EST microsatellite sequence (CK309496) which displays sequence homology to the 5' untranslated region (UTR) of the avian spindlin gene. This gene is known to be present on the Z and W chromosomes. To maximize cross-species utility, the primer set was designed from a consensus sequence created from homologs of CK309496 that were isolated from multiple distantly related species. Both the forward and reverse primer sequences were 100% identical to 14 avian species, including the Z chromosome of eight species and the chicken Gallus gallus W chromosome, as well as the saltwater crocodile Crocodylus porosus. The Z43B primer set was assessed by genotyping individuals of known sex belonging to 61 non-ratite species and a single ratite. The Z and W amplicons differed in size making it possible to distinguish between males (ZZ) and females (ZW) for the majority (69%) of non-ratite species tested, comprising 10 orders of birds. We predict that this marker will be useful for obtaining sex-typing data for ca 6,869 species of birds (69% of non-ratites but not galliforms). A wide range of species could be sex-typed including passerines, shorebirds, eagles, falcons, bee-eaters, cranes, shags, parrots, penguins, ducks, and a ratite species, the brown kiwi, Apteryx australis. Those species sexed include species impossible or problematic to sex-type with other markers (magpie, albatross, petrel, eagle, falcon, crane, and penguin species). Zoo Biol. 35:533-545, 2016. © 2016 The Authors. Zoo Biology published by Wiley Periodicals, Inc.

Quantitative Trait Locus Analysis of Mating Behavior and Male Sex Pheromones in Nasonia Wasps.

A major focus in speciation genetics is to identify the chromosomal regions and genes that reduce hybridization and gene flow. We investigated the genetic architecture of mating behavior in the parasitoid wasp species pair Nasonia giraulti and Nasonia oneida that exhibit strong prezygotic isolation. Behavioral analysis showed that N. oneida females had consistently higher latency times, and broke off the mating sequence more often in the mounting stage when confronted with N. giraulti males compared with males of their own species. N. oneida males produce a lower quantity of the long-range male sex pheromone (4R,5S)-5-hydroxy-4-decanolide (RS-HDL). Crosses between the two species yielded hybrid males with various pheromone quantities, and these males were used in mating trials with females of either species to measure female mate discrimination rates. A quantitative trait locus (QTL) analysis involving 475 recombinant hybrid males (F2), 2148 reciprocally backcrossed females (F3), and a linkage map of 52 equally spaced neutral single nucleotide polymorphism (SNP) markers plus SNPs in 40 candidate mating behavior genes revealed four QTL for male pheromone amount, depending on partner species. Our results demonstrate that the RS-HDL pheromone plays a role in the mating system of N. giraulti and N. oneida, but also that additional communication cues are involved in mate choice. No QTL were found for female mate discrimination, which points at a polygenic architecture of female choice with strong environmental influences.

A multiplex microsatellite set for non-invasive genotyping and sexing of the osprey (Pandion haliaetus).

During the 1950s and 1970s the osprey (Pandion haliaetus) experienced a dramatic population crash and remains of conservation concern in several parts of the world. We isolated 37 microsatellite loci and assessed these in ospreys sampled in the UK and Norway (using mouth swabs/feathers). From 26 loci variable in four ospreys, we selected 13, combined these into two multiplex-PCR sets and included a sex-typing marker. Additional markers confirmed sexes. In 17 ospreys, feather-sampled in central Norway, we found 3-10 alleles per locus. The 13 loci are autosomal (heterozygotes were present in both sexes) and observed heterozygosities ranged from 0.24 to 0.94. The combined probability of identity for the 13 loci was 8.0 × 10-12. These microsatellite loci will be useful for genetic monitoring, parentage analysis and population genetic studies of the osprey.

Oceanography and life history predict contrasting genetic population structure in two Antarctic fish species.

Understanding the key drivers of population connectivity in the marine environment is essential for the effective management of natural resources. Although several different approaches to evaluating connectivity have been used, they are rarely integrated quantitatively. Here, we use a 'seascape genetics' approach, by combining oceanographic modelling and microsatellite analyses, to understand the dominant influences on the population genetic structure of two Antarctic fishes with contrasting life histories, Champsocephalus gunnari and Notothenia rossii. The close accord between the model projections and empirical genetic structure demonstrated that passive dispersal during the planktonic early life stages is the dominant influence on patterns and extent of genetic structuring in both species. The shorter planktonic phase of C. gunnari restricts direct transport of larvae between distant populations, leading to stronger regional differentiation. By contrast, geographic distance did not affect differentiation in N. rossii, whose longer larval period promotes long-distance dispersal. Interannual variability in oceanographic flows strongly influenced the projected genetic structure, suggesting that shifts in circulation patterns due to climate change are likely to impact future genetic connectivity and opportunities for local adaptation, resilience and recovery from perturbations. Further development of realistic climate models is required to fully assess such potential impacts.

Characterization of the house sparrow (Passer domesticus) transcriptome: a resource for molecular ecology and immunogenetics.

The house sparrow (Passer domesticus) is an important model species in ecology and evolution. However, until recently, genomic resources for molecular ecological projects have been lacking in this species. Here, we present transcriptome sequencing data (RNA-Seq) from three different house sparrow tissues (spleen, blood and bursa). These tissues were specifically chosen to obtain a diverse representation of expressed genes and to maximize the yield of immune-related gene functions. After de novo assembly, 15,250 contigs were identified, representing sequence data from a total of 8756 known avian genes (as inferred from the closely related zebra finch). The transcriptome assembly contain sequence data from nine manually annotated MHC genes, including an almost complete MHC class I coding sequence. There were 407, 303 and 68 genes overexpressed in spleen, blood and bursa, respectively. Gene ontology terms related to ribosomal function were associated with overexpression in spleen and oxygen transport functions with overexpression in blood. In addition to the transcript sequences, we provide 327 gene-linked microsatellites (SSRs) with sufficient flanking sequences for primer design, and 3177 single-nucleotide polymorphisms (SNPs) within genes, that can be used in follow-up molecular ecology studies of this ecological well-studied species.

High-utility conserved avian microsatellite markers enable parentage and population studies across a wide range of species.

BACKGROUND: Microsatellites are widely used for many genetic studies. In contrast to single nucleotide polymorphism (SNP) and genotyping-by-sequencing methods, they are readily typed in samples of low DNA quality/concentration (e.g. museum/non-invasive samples), and enable the quick, cheap identification of species, hybrids, clones and ploidy. Microsatellites also have the highest cross-species utility of all types of markers used for genotyping, but, despite this, when isolated from a single species, only a relatively small proportion will be of utility. Marker development of any type requires skill and time. The availability of sufficient "off-the-shelf" markers that are suitable for genotyping a wide range of species would not only save resources but also uniquely enable new comparisons of diversity among taxa at the same set of loci. No other marker types are capable of enabling this. We therefore developed a set of avian microsatellite markers with enhanced cross-species utility. RESULTS: We selected highly-conserved sequences with a high number of repeat units in both of two genetically distant species. Twenty-four primer sets were designed from homologous sequences that possessed at least eight repeat units in both the zebra finch (Taeniopygia guttata) and chicken (Gallus gallus). Each primer sequence was a complete match to zebra finch and, after accounting for degenerate bases, at least 86% similar to chicken. We assessed primer-set utility by genotyping individuals belonging to eight passerine and four non-passerine species. The majority of the new Conserved Avian Microsatellite (CAM) markers amplified in all 12 species tested (on average, 94% in passerines and 95% in non-passerines). This new marker set is of especially high utility in passerines, with a mean 68% of loci polymorphic per species, compared with 42% in non-passerine species. CONCLUSIONS: When combined with previously described conserved loci, this new set of conserved markers will not only reduce the necessity and expense of microsatellite isolation for a wide range of genetic studies, including avian parentage and population analyses, but will also now enable comparisons of genetic diversity among different species (and populations) at the same set of loci, with no or reduced bias. Finally, the approach used here can be applied to other taxa in which appropriate genome sequences are available.

Hybridization but no evidence for backcrossing and introgression in a sympatric population of great reed warblers and clamorous reed warblers.

Hybridization is observed frequently in birds, but often it is not known whether the hybrids are fertile and if backcrossing occurs. The breeding ranges of the great reed warbler (Acrocephalus arundinaceus) and the clamorous reed warbler (A. stentoreus) overlap in southern Kazakhstan and a previous study has documented hybridization in a sympatric population. In the present study, we first present a large set of novel microsatellite loci isolated and characterised in great reed warblers. Secondly, we evaluate whether hybridization in the sympatric breeding population has been followed by backcrossing and introgression.We isolated 181 unique microsatellite loci in great reed warblers. Of 41 loci evaluated, 40 amplified and 30 were polymorphic. Bayesian clustering analyses based on genotype data from 23 autosomal loci recognised two well-defined genetic clusters corresponding to the two species. Individuals clustered to a very high extent to either of these clusters (admixture proportions ≥ 0.984) with the exception of four previously suggested arundinaceus-stentoreus hybrid birds that showed mixed ancestry (admixture proportions 0.495-0.619). Analyses of simulated hybrids and backcrossed individuals showed that the sampled birds do not correspond to first-fourth-generation backcrosses, and that fifth or higher generation backcrosses to a high extent resemble 'pure' birds at this set of markers.We conclude that these novel microsatellite loci provide a useful molecular resource for Acrocephalus warblers. The time to reach reproductive isolation is believed to be very long in birds, approximately 5 Myrs, and with an estimated divergence time of 2 Myrs between these warblers, some backcrossing and introgression could have been expected. However, there was no evidence for backcrossing and introgression suggesting that hybrids are either infertile or their progeny inviable. Very low levels of introgression cannot be excluded, which still may be an important factor as a source of new genetic variation.

Comparison between Normalised and Unnormalised 454-Sequencing Libraries for Small-Scale RNA-Seq Studies.

Next-generation sequencing of transcriptomes (RNA-Seq) is being used increasingly in studies of nonmodel organisms. Here, we evaluate the effectiveness of normalising cDNA libraries prior to sequencing in a small-scale study of the zebra finch. We find that assemblies produced from normalised libraries had a larger number of contigs but used fewer reads compared to unnormalised libraries. Considerably more genes were also detected using the contigs produced from normalised cDNA, and microsatellite discovery was up to 73% more efficient in these. There was a positive correlation between the detected expression level of genes in normalised and unnormalised cDNA, and there was no difference in the number of genes identified as being differentially expressed between blood and spleen for the normalised and unnormalised libraries. We conclude that normalised cDNA libraries are preferable for many applications of RNA-Seq and that these can also be used in quantitative gene expression studies.