Single-species microarrays and comparative transcriptomics

Sex-Biased Expression of Young Genes in Silurana (Xenopus) tropicalis

Sex-biased gene expression can evolve from sex-specific selection and is often associated with sex-linked genes. Gene duplication is a particularly effective mechanism for the generation of sex-biased genes, in which a new copy can help resolve intralocus sexual conflicts. This study assesses sex-biased gene expression in an amphibian with homomorphic ZW sex chromosomes, the Western clawed frog Silurana (Xenopus)tropicalis. Previous work has shown that the sex chromosomes in this species are mainly undifferentiated and pseudoautosomal. Consistent with ongoing recombination between the sex chromosomes, this study detected little evidence for the general sexualization of sex-linked regions. A subset of genes closely linked to the sex determining locus displays a tendency for male-biased expression and elevated rates of evolution relative to genes in other genomic locations. This may be a symptom of an early stage of sex chromosome differentiation driven by, for example, chromosomal degeneration or natural selection on genes in this portion of the Z chromosome. Alternatively, it could reflect variation between the sexes in allelic copy number coupled with a lack of dosage compensation. Irrespective of the genomic location, lineage-specific genes and recently duplicated genes had significantly high levels of sex-biased expression, offering insights into the early transcriptional differentiation of young genes.

Comparative transcriptomics of Atlantic Salmo salar, chum Oncorhynchus keta and pink salmon O. gorbuscha during infections with salmon lice Lepeophtheirus salmonis

Background

Salmon species vary in susceptibility to infections with the salmon louse (Lepeophtheirus salmonis). Comparing mechanisms underlying responses in susceptible and resistant species is important for estimating impacts of infections on wild salmon, selective breeding of farmed salmon, and expanding our knowledge of fish immune responses to ectoparasites. Herein we report three L. salmonis experimental infection trials of co-habited Atlantic Salmo salar, chum Oncorhynchus keta and pink salmon O. gorbuscha, profiling hematocrit, blood cortisol concentrations, and transcriptomic responses of the anterior kidney and skin to the infection.

Results

In all trials, infection densities (lice per host weight (g)) were consistently highest on chum salmon, followed by Atlantic salmon, and lowest in pink salmon. At 43 days post-exposure, all lice had developed to motile stages, and infection density was uniformly low among species. Hematocrit was reduced in infected Atlantic and chum salmon, and cortisol was elevated in infected chum salmon. Systemic transcriptomic responses were profiled in all species and large differences in response functions were identified between Atlantic and Pacific (chum and pink) salmon. Pink and chum salmon up-regulated acute phase response genes, including complement and coagulation components, and down-regulated antiviral immune genes. The pink salmon response involved the largest and most diverse iron sequestration and homeostasis mechanisms. Pattern recognition receptors were up-regulated in all species but the active components were often species-specific. C-type lectin domain family 4 member M and acidic mammalian chitinase were specifically up-regulated in the resistant pink salmon.

Conclusions

Experimental exposures consistently indicated increased susceptibility in chum and Atlantic salmon, and resistance in pink salmon, with differences in infection density occurring within the first three days of infection. Transcriptomic analysis suggested candidate resistance functions including local inflammation with cytokines, specific innate pattern recognition receptors, and iron homeostasis. Suppressed antiviral immunity in both susceptible and resistant species indicates the importance of future work investigating co-infections of viral pathogens and lice.

Electronic supplementary material

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

A comparison of host-defense peptides in skin secretions of female Xenopus laevis × Xenopus borealis and X. borealis × X. laevis F1 hybrids

Peptidomic analysis was used to compare the diversity of host-defense peptides in norepinephrine-stimulated skin secretions from laboratory-generated female F1 hybrids of Xenopus laevis and Xenopus borealis (Pipidae). Skin secretions of hybrids with maternal X. laevis (XLB) contained 12 antimicrobial peptides (AMPs), comprising 8 from X. laevis and 4 from X. borealis. Magainin-B1, XPF-B1, PGLa-B1 CPF-B2, CPF-B3 and CPF-B4 from X. borealis and XPF-1, XPF-2, and CPF-6 from X. laevis were not detected and CPF-1 and CPF-7 were present in low concentration. The secretions contained caerulein and caerulein-B1 derived from both parents but lacked X. laevis xenopsin and X. borealis caerulein-B2. Skin secretions of hybrids with maternal X. borealis (XBL) contained 14 AMPs comprising 6 from X. borealis and 8 from X. laevis. Magainin-B1, XPF-B1, PGLa-B1, CPF-B2, XPF-1, CPF-5, and CPF-7 were absent and CPF-B3, CPF-B4, CPF-1 and CPF-6 were present only in low concentration. Xenopsin and caerulein were identified in the secretions but caerulein-B2 was absent and caerulein-B1 was present in low concentration. No peptides were identified in secretions of either XLB or XBL hybrids that were not present in the parental species. The data indicate that hybridization between X. laevis and X. borealis results in increased diversity of host-defense peptides in skin secretions but point to extensive AMP gene silencing compared with previously studied female X. laevis×X. muelleri F1 hybrids and no novel peptide expression.

Expression divergence measured by transcriptome sequencing of four yeast species

Background

The evolution of gene expression is a challenging problem in evolutionary biology, for which accurate, well-calibrated measurements and methods are crucial.

Results

We quantified gene expression with whole-transcriptome sequencing in four diploid, prototrophic strains of Saccharomyces species grown under the same condition to investigate the evolution of gene expression. We found that variation in expression is gene-dependent with large variations in each gene's expression between replicates of the same species. This confounds the identification of genes differentially expressed across species. To address this, we developed a statistical approach to establish significance bounds for inter-species differential expression in RNA-Seq data based on the variance measured across biological replicates. This metric estimates the combined effects of technical and environmental variance, as well as Poisson sampling noise by isolating each component. Despite a paucity of large expression changes, we found a strong correlation between the variance of gene expression change and species divergence (R² = 0.90).

Conclusion

We provide an improved methodology for measuring gene expression changes in evolutionary diverged species using RNA Seq, where experimental artifacts can mimic evolutionary effects.

GEO Accession Number: GSE32679

Toward an unbiased evolutionary platform for unraveling Xenopus developmental gene networks

The availability of both the Xenopus tropicalis genome and the soon to be released Xenopus laevis genome provides a solid foundation for Xenopus developmental biologists. The Xenopus community has presently amassed expression data for ∼2,300 genes in the form of published images collected in the Xenbase, the principal Xenopus research database. A few of these genes have been examined in both X. tropicalis and X. laevis and the cross-species comparison has been proven invaluable for studying gene function. A recently published work has yielded developmental expression profiles for the majority of Xenopus genes across fourteen developmental stages spanning the blastula, gastrula, neurula, and the tail-bud. While this data was originally queried for global evolutionary and developmental principles, here we demonstrate its general use for gene-level analyses. In particular, we present the accessibility of this dataset through Xenbase and describe biases in the characterized genes in terms of sequence and expression conservation across the two species. We further indicate the advantage of examining coexpression for gene function discovery relating to developmental processes conserved across species. We suggest that the integration of additional large-scale datasets--comprising diverse functional data--into Xenbase promises to provide a strong foundation for researchers in elucidating biological processes including the gene regulatory programs encoding development.

Snapshot of DNA methylation changes associated with hybridization in Xenopus

Hybridization often results in dramatic genome reconfigurations including epigenetic changes that control gene expression. Here we survey methylation patterns of interspecific Xenopus F1 hybrids relative to parental species X. laevis and X. muelleri, using methyl-sensitive amplification polymorphisms (MSAPs). Out of a total of 546 MSAP markers, 364 were effective in elucidating the difference in methylation patterns between the hybrids and the parental species. Principal coordinate analysis of methylated fragments revealed four distinct clusters with the two parental species separate from hybrid males and females. On average, hybrids were characterized by a higher proportion (70.6%) of methylated fragments compared with the parental species (64.5%), and this difference was consistent with previously observed disruptions of hybrid transcriptomes. The proportion of methylated fragments did not correlate with variation in genome size, as measured with flow cytometry. The levels of methylation in sterile hybrid males (73.8%) were higher than in fertile hybrid females (68.6%), but this difference was not statistically significant. A total of 76 methylated fragments (20.9%) were hybrid-unique, presumably originating from methylation alterations in hybrid genomes.

Identification of suitable reference genes for normalization of qPCR data in comparative transcriptomics analyses in the Triticeae

Comparative transcriptomics are useful to determine the role of orthologous genes among Triticeae species. Thus they constitute an interesting tool to improve the use of wild relatives for crop breeding. Reverse transcription quantitative real-time PCR (qPCR) is the most accurate measure of gene expression but efficient normalization is required. The choice and optimal number of reference genes must be experimentally determined and the primers optimized for cross-species amplification. Our goal was to test the utility of wheat-reference genes for qPCR normalization when species carrying the following genomes (A, B, D, R, H ( v ) and H ( ch )) are compared either simultaneously or in smaller subsets of samples. Wheat/barley/rye consensus primers outperformed wheat-specific ones which indicate that consensus primers should be considered for data normalization in comparative transcriptomics. All genes tested were stable but their ranking in terms of stability differed among subsets of samples. CDC (cell division control protein, AAA-superfamily of ATPases, Ta54227) and RLI (68 kDa protein HP68 similar to Arabidopsis thaliana RNase L inhibitor protein, Ta2776) were always among the three most stable genes. The optimal number of reference genes varied between 2 and 3 depending on the subset of samples and the method used (geNorm vs. coefficient of determination between sequential normalization factors). In any case a maximum number of three reference genes would provide adequate normalization independent of the subset of samples considered. This work constitutes a substantial advance towards comparative transcriptomics using qPCR since it provides useful primers/reference genes.

Comparative functional genomic analysis of Solanum glandular trichome types

Glandular trichomes play important roles in protecting plants from biotic attack by producing defensive compounds. We investigated the metabolic profiles and transcriptomes to characterize the differences between different glandular trichome types in several domesticated and wild Solanum species: Solanum lycopersicum (glandular trichome types 1, 6, and 7), Solanum habrochaites (types 1, 4, and 6), Solanum pennellii (types 4 and 6), Solanum arcanum (type 6), and Solanum pimpinellifolium (type 6). Substantial chemical differences in and between Solanum species and glandular trichome types are likely determined by the regulation of metabolism at several levels. Comparison of S. habrochaites type 1 and 4 glandular trichomes revealed few differences in chemical content or transcript abundance, leading to the conclusion that these two glandular trichome types are the same and differ perhaps only in stalk length. The observation that all of the other species examined here contain either type 1 or 4 trichomes (not both) supports the conclusion that these two trichome types are the same. Most differences in metabolites between type 1 and 4 glands on the one hand and type 6 glands on the other hand are quantitative but not qualitative. Several glandular trichome types express genes associated with photosynthesis and carbon fixation, indicating that some carbon destined for specialized metabolism is likely fixed within the trichome secretory cells. Finally, Solanum type 7 glandular trichomes do not appear to be involved in the biosynthesis and storage of specialized metabolites and thus likely serve another unknown function, perhaps as the site of the synthesis of protease inhibitors.

Transcript profiling of two alfalfa genotypes with contrasting cell wall composition in stems using a cross-species platform: optimizing analysis by masking biased probes

BACKGROUND

The GeneChip(R) Medicago Genome Array, developed for Medicago truncatula, is a suitable platform for transcript profiling in tetraploid alfalfa [Medicago sativa (L.) subsp. sativa]. However, previous research involving cross-species hybridization (CSH) has shown that sequence variation between two species can bias transcript profiling by decreasing sensitivity (number of expressed genes detected) and the accuracy of measuring fold-differences in gene expression.

RESULTS

Transcript profiling using the Medicago GeneChip(R) was conducted with elongating stem (ES) and post-elongation stem (PES) internodes from alfalfa genotypes 252 and 1283 that differ in stem cell wall concentrations of cellulose and lignin. A protocol was developed that masked probes targeting inter-species variable (ISV) regions of alfalfa transcripts. A probe signal intensity threshold was selected that optimized both sensitivity and accuracy. After masking for both ISV regions and previously identified single-feature polymorphisms (SFPs), the number of differentially expressed genes between the two genotypes in both ES and PES internodes was approximately 2-fold greater than the number detected prior to masking. Regulatory genes, including transcription factor and receptor kinase genes that may play a role in development of secondary xylem, were significantly over-represented among genes up-regulated in 252 PES internodes compared to 1283 PES internodes. Several cell wall-related genes were also up-regulated in genotype 252 PES internodes. Real-time quantitative RT-PCR of differentially expressed regulatory and cell wall-related genes demonstrated increased sensitivity and accuracy after masking for both ISV regions and SFPs. Over 1,000 genes that were differentially expressed in ES and PES internodes of genotypes 252 and 1283 were mapped onto putative orthologous loci on M. truncatula chromosomes. Clustering simulation analysis of the differentially expressed genes suggested co-expression of some neighbouring genes on Medicago chromosomes.

CONCLUSIONS

The problems associated with transcript profiling in alfalfa stems using the Medicago GeneChip as a CSH platform were mitigated by masking probes targeting ISV regions and SFPs. Using this masking protocol resulted in the identification of numerous candidate genes that may contribute to differences in cell wall concentration and composition of stems of two alfalfa genotypes.

A simple optimization can improve the performance of single feature polymorphism detection by Affymetrix expression arrays

Background

High-density oligonucleotide arrays are effective tools for genotyping numerous loci simultaneously. In small genome species (genome size: < ~300 Mb), whole-genome DNA hybridization to expression arrays has been used for various applications. In large genome species, transcript hybridization to expression arrays has been used for genotyping. Although rice is a fully sequenced model plant of medium genome size (~400 Mb), there are a few examples of the use of rice oligonucleotide array as a genotyping tool.

Results

We compared the single feature polymorphism (SFP) detection performance of whole-genome and transcript hybridizations using the Affymetrix GeneChip^® Rice Genome Array, using the rice cultivars with full genome sequence, japonica cultivar Nipponbare and indica cultivar 93-11. Both genomes were surveyed for all probe target sequences. Only completely matched 25-mer single copy probes of the Nipponbare genome were extracted, and SFPs between them and 93-11 sequences were predicted. We investigated optimum conditions for SFP detection in both whole genome and transcript hybridization using differences between perfect match and mismatch probe intensities of non-polymorphic targets, assuming that these differences are representative of those between mismatch and perfect targets. Several statistical methods of SFP detection by whole-genome hybridization were compared under the optimized conditions. Causes of false positives and negatives in SFP detection in both types of hybridization were investigated.

Conclusions

The optimizations allowed a more than 20% increase in true SFP detection in whole-genome hybridization and a large improvement of SFP detection performance in transcript hybridization. Significance analysis of the microarray for log-transformed raw intensities of PM probes gave the best performance in whole genome hybridization, and 22,936 true SFPs were detected with 23.58% false positives by whole genome hybridization. For transcript hybridization, stable SFP detection was achieved for highly expressed genes, and about 3,500 SFPs were detected at a high sensitivity (> 50%) in both shoot and young panicle transcripts. High SFP detection performances of both genome and transcript hybridizations indicated that microarrays of a complex genome (e.g., of Oryza sativa) can be effectively utilized for whole genome genotyping to conduct mutant mapping and analysis of quantitative traits such as gene expression levels.

Transcript profiling of common bean (Phaseolus vulgaris L.) using the GeneChip Soybean Genome Array: optimizing analysis by masking biased probes

BACKGROUND

Common bean (Phaseolus vulgaris L.) and soybean (Glycine max) both belong to the Phaseoleae tribe and share significant coding sequence homology. This suggests that the GeneChip(R) Soybean Genome Array (soybean GeneChip) may be used for gene expression studies using common bean.

RESULTS

To evaluate the utility of the soybean GeneChip for transcript profiling of common bean, we hybridized cRNAs purified from nodule, leaf, and root of common bean and soybean in triplicate to the soybean GeneChip. Initial data analysis showed a decreased sensitivity and accuracy of measuring differential gene expression in common bean cross-species hybridization (CSH) GeneChip data compared to that of soybean. We employed a method that masked putative probes targeting inter-species variable (ISV) regions between common bean and soybean. A masking signal intensity threshold was selected that optimized both sensitivity and accuracy of measuring differential gene expression. After masking for ISV regions, the number of differentially-expressed genes identified in common bean was increased by 2.8-fold reflecting increased sensitivity. Quantitative RT-PCR (qRT-PCR) analysis of 20 randomly selected genes and purine-ureide pathway genes demonstrated an increased accuracy of measuring differential gene expression after masking for ISV regions. We also evaluated masked probe frequency per probe set to gain insight into the sequence divergence pattern between common bean and soybean. The sequence divergence pattern analysis suggested that the genes for basic cellular functions and metabolism were highly conserved between soybean and common bean. Additionally, our results show that some classes of genes, particularly those associated with environmental adaptation, are highly divergent.

CONCLUSIONS

The soybean GeneChip is a suitable cross-species platform for transcript profiling in common bean when used in combination with the masking protocol described. In addition to transcript profiling, CSH of the GeneChip in combination with masking probes in the ISV regions can be used for comparative ecological and/or evolutionary genomics studies.

Microarray comparative genomic hybridisation analysis incorporating genomic organisation, and application to enterobacterial plant pathogens

Microarray comparative genomic hybridisation (aCGH) provides an estimate of the relative abundance of genomic DNA (gDNA) taken from comparator and reference organisms by hybridisation to a microarray containing probes that represent sequences from the reference organism. The experimental method is used in a number of biological applications, including the detection of human chromosomal aberrations, and in comparative genomic analysis of bacterial strains, but optimisation of the analysis is desirable in each problem domain.We present a method for analysis of bacterial aCGH data that encodes spatial information from the reference genome in a hidden Markov model. This technique is the first such method to be validated in comparisons of sequenced bacteria that diverge at the strain and at the genus level: Pectobacterium atrosepticum SCRI1043 (Pba1043) and Dickeya dadantii 3937 (Dda3937); and Lactococcus lactis subsp. lactis IL1403 and L. lactis subsp. cremoris MG1363. In all cases our method is found to outperform common and widely used aCGH analysis methods that do not incorporate spatial information. This analysis is applied to comparisons between commercially important plant pathogenic soft-rotting enterobacteria (SRE) Pba1043, P. atrosepticum SCRI1039, P. carotovorum 193, and Dda3937.Our analysis indicates that it should not be assumed that hybridisation strength is a reliable proxy for sequence identity in aCGH experiments, and robustly extends the applicability of aCGH to bacterial comparisons at the genus level. Our results in the SRE further provide evidence for a dynamic, plastic 'accessory' genome, revealing major genomic islands encoding gene products that provide insight into, and may play a direct role in determining, variation amongst the SRE in terms of their environmental survival, host range and aetiology, such as phytotoxin synthesis, multidrug resistance, and nitrogen fixation.

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.