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

The evolutionary legacy of size-selective harvesting extends from genes to populations

Size-selective harvesting is assumed to alter life histories of exploited fish populations, thereby negatively affecting population productivity, recovery, and yield. However, demonstrating that fisheries-induced phenotypic changes in the wild are at least partly genetically determined has proved notoriously difficult. Moreover, the population-level consequences of fisheries-induced evolution are still being controversially discussed. Using an experimental approach, we found that five generations of size-selective harvesting altered the life histories and behavior, but not the metabolic rate, of wild-origin zebrafish (Danio rerio). Fish adapted to high positively size selective fishing pressure invested more in reproduction, reached a smaller adult body size, and were less explorative and bold. Phenotypic changes seemed subtle but were accompanied by genetic changes in functional loci. Thus, our results provided unambiguous evidence for rapid, harvest-induced phenotypic and evolutionary change when harvesting is intensive and size selective. According to a life-history model, the observed life-history changes elevated population growth rate in harvested conditions, but slowed population recovery under a simulated moratorium. Hence, the evolutionary legacy of size-selective harvesting includes populations that are productive under exploited conditions, but selectively disadvantaged to cope with natural selection pressures that often favor large body size.Peer reviewe

A polygenic basis for birth weight in a wild population of red deer ( Cervus elaphus )

The genetic architecture of traits under selection has important consequences for the response to selection and potentially for population viability. Early QTL mapping studies in wild populations have reported loci with large effect on trait variation. However, these results are contradicted by more recent genome-wide association analyses, which strongly support the idea that most quantitative traits have a polygenic basis. This study aims to re-evaluate the genetic architecture of a key morphological trait, birth weight, in a wild population of red deer (Cervus elaphus), using genomic approaches. A previous study using 93 microsatellite and allozyme markers and linkage mapping on a kindred of 364 deer detected a pronounced QTL on chromosome 21 explaining 29% of the variance in birth weight, suggesting that this trait is partly controlled by genes with large effects. Here, we used data for more than 2,300 calves genotyped at >39,000 SNP markers and two approaches to characterise the genetic architecture of birth weight. First, we performed a genome-wide association (GWA) analysis, using a genomic relatedness matrix to account for population structure. We found no SNPs significantly associated with birth weight. Second, we used genomic prediction to estimate the proportion of variance explained by each SNP and chromosome. This analysis confirmed that most genetic variance in birth weight was explained by loci with very small effect sizes. Third, we found that the proportion of variance explained by each chromosome was slightly positively correlated with its size. These three findings highlight a highly polygenic architecture for birth weight, which contradicts the previous QTL study. These results are probably explained by the differences in how associations are modelled between QTL mapping and GWA. Our study suggests that models of polygenic adaptation are the most appropriate to study the evolutionary trajectory of this trait

Associations of ATR and CHEK1 Single Nucleotide Polymorphisms with Breast Cancer

DNA damage and replication checkpoints mediated by the ATR-CHEK1 pathway are key to the maintenance of genome stability, and both ATR and CHEK1 have been proposed as potential breast cancer susceptibility genes. Many novel variants recently identified by the large resequencing projects have not yet been thoroughly tested in genome-wide association studies for breast cancer susceptibility. We therefore used a tagging SNP (tagSNP) approach based on recent SNP data available from the 1000 genomes projects, to investigate the roles of ATR and CHEK1 in breast cancer risk and survival. ATR and CHEK1 tagSNPs were genotyped in the Sheffield Breast Cancer Study (SBCS; 1011 cases and 1024 controls) using Illumina GoldenGate assays. Untyped SNPs were imputed using IMPUTE2, and associations between genotype and breast cancer risk and survival were evaluated using logistic and Cox proportional hazard regression models respectively on a per allele basis. Significant associations were further examined in a meta-analysis of published data or confirmed in the Utah Breast Cancer Study (UBCS). The most significant associations for breast cancer risk in SBCS came from rs6805118 in ATR (p=7.6x10-5) and rs2155388 in CHEK1 (p=3.1x10-6), but neither remained significant after meta-analysis with other studies. However, meta-analysis of published data revealed a weak association between the ATR SNP rs1802904 (minor allele frequency is 12%) and breast cancer risk, with a summary odds ratio (confidence interval) of 0.90 (0.83-0.98) [p=0.0185] for the minor allele. Further replication of this SNP in larger studies is warranted since it is located in the target region of 2 microRNAs. No evidence of any survival effects of ATR or CHEK1 SNPs were identified. We conclude that common alleles of ATR and CHEK1 are not implicated in breast cancer risk or survival, but we cannot exclude effects of rare alleles and of common alleles with very small effect sizes

DNA methylation differences between stick insect ecotypes

International audienceAbstract Epigenetic mechanisms, such as DNA methylation, can influence gene regulation and affect phenotypic variation, raising the possibility that they contribute to ecological adaptation. Beginning to address this issue requires high‐resolution sequencing studies of natural populations to pinpoint epigenetic regions of potential ecological and evolutionary significance. However, such studies are still relatively uncommon, especially in insects, and are mainly restricted to a few model organisms. Here, we characterize patterns of DNA methylation for natural populations of Timema cristinae adapted to two host plant species (i.e. ecotypes). By integrating results from sequencing of whole transcriptomes, genomes and methylomes, we investigate whether environmental, host and genetic differences of these stick insects are associated with methylation levels of cytosine nucleotides in the CpG context. We report an overall genome‐wide methylation level for T. cristinae of ~14%, with methylation being enriched in gene bodies and impoverished in repetitive elements. Genome‐wide DNA methylation variation was strongly positively correlated with genetic distance (relatedness), but also exhibited significant host‐plant effects. Using methylome‐environment association analysis, we pinpointed specific genomic regions that are differentially methylated between ecotypes, with these regions being enriched for genes with functions in membrane processes. The observed association between methylation variation and genetic relatedness, and with the ecologically important variable of host plant, suggests a potential role for epigenetic modification in T. cristinae adaptation. To substantiate such adaptive significance, future studies could test whether methylation can be transmitted across generations and the extent to which it responds to experimental manipulation in field and laboratory studies

Males and Females Contribute Unequally to Offspring Genetic Diversity in the Polygynandrous Mating System of Wild Boar

The maintenance of genetic diversity across generations depends on both the number of reproducing males and females. Variance in reproductive success, multiple paternity and litter size can all affect the relative contributions of male and female parents to genetic variation of progeny. The mating system of the wild boar (Sus scrofa) has been described as polygynous, although evidence of multiple paternity in litters has been found. Using 14 microsatellite markers, we evaluated the contribution of males and females to genetic variation in the next generation in independent wild boar populations from the Iberian Peninsula and Hungary. Genetic contributions of males and females were obtained by distinguishing the paternal and maternal genetic component inherited by the progeny. We found that the paternally inherited genetic component of progeny was more diverse than the maternally inherited component. Simulations showed that this finding might be due to a sampling bias. However, after controlling for the bias by fitting both the genetic diversity in the adult population and the number of reproductive individuals in the models, paternally inherited genotypes remained more diverse than those inherited maternally. Our results suggest new insights into how promiscuous mating systems can help maintain genetic variation

Genetic Predictions of Prion Disease Susceptibility in Carnivore Species Based on Variability of the Prion Gene Coding Region

Mammalian species vary widely in their apparent susceptibility to prion diseases. For example, several felid species developed prion disease (feline spongiform encephalopathy or FSE) during the bovine spongiform encephalopathy (BSE) epidemic in the United Kingdom, whereas no canine BSE cases were detected. Whether either of these or other groups of carnivore species can contract other prion diseases (e.g. chronic wasting disease or CWD) remains an open question. Variation in the host-encoded prion protein (PrP(C)) largely explains observed disease susceptibility patterns within ruminant species, and may explain interspecies differences in susceptibility as well. We sequenced and compared the open reading frame of the PRNP gene encoding PrP(C) protein from 609 animal samples comprising 29 species from 22 genera of the Order Carnivora; amongst these samples were 15 FSE cases. Our analysis revealed that FSE cases did not encode an identifiable disease-associated PrP polymorphism. However, all canid PrPs contained aspartic acid or glutamic acid at codon 163 which we propose provides a genetic basis for observed susceptibility differences between canids and felids. Among other carnivores studied, wolverine (Gulo gulo) and pine marten (Martes martes) were the only non-canid species to also express PrP-Asp163, which may impact on their prion diseases susceptibility. Populations of black bear (Ursus americanus) and mountain lion (Puma concolor) from Colorado showed little genetic variation in the PrP protein and no variants likely to be highly resistant to prions in general, suggesting that strain differences between BSE and CWD prions also may contribute to the limited apparent host range of the latter

Determinants of the efficacy of natural selection on coding and noncoding variability in two passerine species

Population genetic theory predicts that selection should be more effective when the effective population size (Ne) is larger, and that the efficacy of selection should correlate positively with recombination rate. Here, we analyzed the genomes of ten great tits and ten zebra finches. Nucleotide diversity at 4-fold degenerate sites indicates that zebra finches have a 2.83-fold larger Ne. We obtained clear evidence that purifying selection is more effective in zebra finches. The proportion of substitutions at 0-fold degenerate sites fixed by positive selection (α) is high in both species (great tit 48%; zebra finch 64%) and is significantly higher in zebra finches. When α was estimated on GC-conservative changes (i.e., between A and T and between G and C), the estimates reduced in both species (great tit 22%; zebra finch 53%). A theoretical model presented herein suggests that failing to control for the effects of GC-biased gene conversion (gBGC) is potentially a contributor to the overestimation of α, and that this effect cannot be alleviated by first fitting a demographic model to neutral variants. We present the first estimates in birds for α in the untranslated regions, and found evidence for substantial adaptive changes. Finally, although purifying selection is stronger in high-recombination regions, we obtained mixed evidence for α increasing with recombination rate, especially after accounting for gBGC. These results highlight that it is important to consider the potential confounding effects of gBGC when quantifying selection and that our understanding of what determines the efficacy of selection is incomplete