High-throughput cDNA Sequencing Research Articles

RXLR and CRN Effectors from the Sunflower Downy Mildew Pathogen Plasmopara halstedii Induce Hypersensitive-Like Responses in Resistant Sunflower Lines.

Plasmopara halstedii is an obligate biotrophic oomycete causing downy mildew disease on sunflower, Helianthus annuus, an economically important oil crop. Severe symptoms of the disease (e.g., plant dwarfism, leaf bleaching, sporulation and production of infertile flower) strongly impair seed yield. Pl resistance genes conferring resistance to specific P. halstedii pathotypes were located on sunflower genetic map but yet not cloned. They are present in cultivated lines to protect them against downy mildew disease. Among the 16 different P. halstedii pathotypes recorded in France, pathotype 710 is frequently found, and therefore continuously controlled in sunflower by different Pl genes. High-throughput sequencing of cDNA from P. halstedii led us to identify potential effectors with the characteristic RXLR or CRN motifs described in other oomycetes. Expression of six P. halstedii putative effectors, five RXLR and one CRN, was analyzed by qRT-PCR in pathogen spores and in the pathogen infecting sunflower leaves and selected for functional analyses. We developed a new method for transient expression in sunflower plant leaves and showed for the first time subcellular localization of P. halstedii effectors fused to a fluorescent protein in sunflower leaf cells. Overexpression of the CRN and of 3 RXLR effectors induced hypersensitive-like cell death reactions in some sunflower near-isogenic lines resistant to pathotype 710 and not in susceptible corresponding lines, suggesting they could be involved in Pl loci-mediated resistances.

Characterization of IgM in Norwegian cleaner fish (lumpfish and wrasses)

The use of cleaner fish in Norwegian aquaculture has to a large extent been based on wild catches, but breeding of lumpfish and ballan wrasse is currently increasing. Due to disease problems and required vaccine development, tools to study immune responses and a better understanding of the immune system in these species is demanded. The present study comprises lumpfish (Cyclopterus lumpus) and five species of wrasses: Ballan wrasse (Labrus bergylta), rock cook (Centrolabrus exoletus), cuckoo wrasse (Labrus mixtus), corkwing wrasse (Symphodus melops), and goldsinny wrasse (Ctenolabrus rupestris). We present a comparison of the IgM sequences, phylogenetic relationship to other teleosts and characteristic features of IgM in the species studied. The lumpfish IgM heavy chain sequence was assembled from high throughput cDNA sequencing whereas the wrasse sequences were determined by molecular cloning. The secreted form of the IgM heavy chain from all species consisted of four constant Ig domains. IgM was purified from lumpfish and ballan wrasse sera by gel filtration followed by anion exchange chromatography, and polyclonal sera were produced against these proteins. Antisera against ballan wrasse IgM showed cross-reactivity to all analyzed species of wrasses, some cross-reactivity to lumpfish, very low reaction to salmon, and no reaction to cod. Anti- IgM sera against lumpfish cross-reacted to the light chain of all species studied. Wrasses and lumpfish IgM showed high binding affinities for protein A. IgM concentration in adult ballan wrasse (700–800 g) was measured by single radial immunodiffusion assay and found to be 13.4 mg/ml which is about 36% of the total protein concentration. The IgM concentration in lumpfish (600–3600 g) was estimated to 1–2.6 mg/ml, which corresponds to approximately 3% of the total protein concentration.

Molecular diversity and distribution pattern of ciliates in sediments from deep-sea hydrothermal vents in the Okinawa Trough and adjacent sea areas

In comparison with the macrobenthos and prokaryotes, patterns of diversity and distribution of microbial eukaryotes in deep-sea hydrothermal vents are poorly known. The widely used high-throughput sequencing of 18S rDNA has revealed a high diversity of microeukaryotes yielded from both living organisms and buried DNA in marine sediments. More recently, cDNA surveys have been utilized to uncover the diversity of active organisms. However, both methods have never been used to evaluate the diversity of ciliates in hydrothermal vents. By using high-throughput DNA and cDNA sequencing of 18S rDNA, we evaluated the molecular diversity of ciliates, a representative group of microbial eukaryotes, from the sediments of deep-sea hydrothermal vents in the Okinawa Trough and compared it with that of an adjacent deep-sea area about 15km away and that of an offshore area of the Yellow Sea about 500km away. The results of DNA sequencing showed that Spirotrichea and Oligohymenophorea were the most diverse and abundant groups in all the three habitats. The proportion of sequences of Oligohymenophorea was the highest in the hydrothermal vents whereas Spirotrichea was the most diverse group at all three habitats. Plagiopyleans were found only in the hydrothermal vents but with low diversity and abundance. By contrast, the cDNA sequencing showed that Plagiopylea was the most diverse and most abundant group in the hydrothermal vents, followed by Spirotrichea in terms of diversity and Oligohymenophorea in terms of relative abundance. A novel group of ciliates, distinctly separate from the 12 known classes, was detected in the hydrothermal vents, indicating undescribed, possibly highly divergent ciliates may inhabit this environment. Statistical analyses showed that: (i) the three habitats differed significantly from one another in terms of diversity of both the rare and the total ciliate taxa, and; (ii) the adjacent deep sea was more similar to the offshore area than to the hydrothermal vents. In terms of the diversity of abundant taxa, however, there was no significant difference between the hydrothermal vents and the adjacent deep sea, both of which differed significantly from the offshore area. As abundant ciliate taxa can be found in several sampling sites, they are likely adapted to large environmental variations, while rare taxa are found in specific habitat and thus are potentially more sensitive to varying environmental conditions.

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.

Differential analyses for RNA-seq: transcript-level estimates improve gene-level inferences.

High-throughput sequencing of cDNA (RNA-seq) is used extensively to characterize the transcriptome of cells. Many transcriptomic studies aim at comparing either abundance levels or the transcriptome composition between given conditions, and as a first step, the sequencing reads must be used as the basis for abundance quantification of transcriptomic features of interest, such as genes or transcripts. Various quantification approaches have been proposed, ranging from simple counting of reads that overlap given genomic regions to more complex estimation of underlying transcript abundances. In this paper, we show that gene-level abundance estimates and statistical inference offer advantages over transcript-level analyses, in terms of performance and interpretability. We also illustrate that the presence of differential isoform usage can lead to inflated false discovery rates in differential gene expression analyses on simple count matrices but thatthis can be addressed by incorporating offsets derived from transcript-level abundance estimates. We also show that the problem is relatively minor in several real data sets. Finally, we provide an R package ( tximport) to help users integrate transcript-level abundance estimates from common quantification pipelines into count-based statistical inference engines.

Lower IgG somatic hypermutation rates during acute dengue virus infection is compatible with a germinal center-independent B cell response.

BackgroundThe study of human B cell response to dengue virus (DENV) infection is critical to understand serotype-specific protection and the cross-reactive sub-neutralizing response. Whereas the first is beneficial and thus represents the ultimate goal of vaccination, the latter has been implicated in the development of severe disease, which occurs in a small, albeit significant, fraction of secondary DENV infections. Both primary and secondary infections are associated with the production of poly-reactive and cross-reactive IgG antibodies.MethodsTo gain insight into the effect of DENV infection on the B cell repertoire, we used VH region high-throughput cDNA sequencing of the peripheral blood IgG B cell compartment of 19 individuals during the acute phase of infection. For 11 individuals, a second sample obtained 6 months later was analyzed for comparison. Probabilities of sequencing antibody secreting cells or memory B cells were estimated using second-order Monte Carlo simulation.ResultsWe found that in acute disease there is an increase in IgG B cell diversity and changes in the relative use of segments IGHV1-2, IGHV1-18, and IGHV1-69. Somewhat unexpectedly, an overall low proportion of somatic hypermutated antibody genes was observed during the acute phase plasmablasts, particularly in secondary infections and those cases with more severe disease.ConclusionsOur data are consistent with an innate-like antiviral recognition system mediated by B cells using defined germ-line coded B cell receptors, which could provide a rapid germinal center-independent antibody response during the early phase of infection. A model describing concurrent T-dependent and T-independent B cell responses in the context of DENV infection is proposed, which incorporates the selection of B cells using hypomutated IGHV segments and their potential role in poly/cross-reactivity. Its formal demonstration could lead to a definition of its potential implication in antibody-dependent enhancement, and may contribute to rational vaccine development efforts.Electronic supplementary materialThe online version of this article (doi:10.1186/s13073-016-0276-1) contains supplementary material, which is available to authorized users.

High-Throughput Nuclease Probing of RNA Structures Using FragSeq.

High-throughput sequencing of cDNA (RNA-Seq) can be used to generate nuclease accessibility data for many distinct transcripts in the same mixture simultaneously. Such assays accelerate RNA structure analysis and provide researchers with new technologies to tackle biological questions on a transcriptome-wide scale. FragSeq is an experimental assay for transcriptome-wide RNA structure probing using RNA-Seq, coupled with data analysis tools that allow quantitative determination of nuclease accessibility at single-base resolution. We provide a practical guide to designing and carrying out FragSeq experiments and data analysis.

Differential analyses for RNA-seq: transcript-level estimates improve gene-level inferences

High-throughput sequencing of cDNA (RNA-seq) is used extensively to characterize the transcriptome of cells. Many transcriptomic studies aim at comparing either abundance levels or the transcriptome composition between given conditions, and as a first step, the sequencing reads must be used as the basis for abundance quantification of transcriptomic features of interest, such as genes or transcripts. Several different quantification approaches have been proposed, ranging from simple counting of reads that overlap given genomic regions to more complex estimation of underlying transcript abundances. In this paper, we show that gene-level abundance estimates and statistical inference offer advantages over transcript-level analyses, in terms of performance and interpretability. We also illustrate that while the presence of differential isoform usage can lead to inflated false discovery rates in differential expression analyses on simple count matrices and transcript-level abundance estimates improve the performance in simulated data, the difference is relatively minor in several real data sets. Finally, we provide an R package (tximport) to help users integrate transcript-level abundance estimates from common quantification pipelines into count-based statistical inference engines.

Interaction of Arabidopsis Trihelix-Domain Transcription Factors VFP3 and VFP5 with Agrobacterium Virulence Protein VirF

Agrobacterium is a natural genetic engineer of plants that exports several virulence proteins into host cells in order to take advantage of the cell machinery to facilitate transformation and support bacterial growth. One of these effectors is the F-box protein VirF, which presumably uses the host ubiquitin/proteasome system (UPS) to uncoat the packaging proteins from the invading bacterial T-DNA. By analogy to several other bacterial effectors, VirF most likely has several functions in the host cell and, therefore, several interacting partners among host proteins. Here we identify one such interactor, an Arabidopsis trihelix-domain transcription factor VFP3, and further show that its very close homolog VFP5 also interacted with VirF. Interestingly, interactions of VirF with either VFP3 or VFP5 did not activate the host UPS, suggesting that VirF might play other UPS-independent roles in bacterial infection. To better understand the potential scope of VFP3 function, we used RNAi to reduce expression of the VFP3 gene. Transcriptome profiling of these VFP3-silenced plants using high-throughput cDNA sequencing (RNA-seq) revealed that VFP3 substantially affected plant gene expression; specifically, 1,118 genes representing approximately 5% of all expressed genes were significantly either up- or down-regulated in the VFP3 RNAi line compared to wild-type Col-0 plants. Among the 507 up-regulated genes were genes implicated in the regulation of transcription, protein degradation, calcium signaling, and hormone metabolism, whereas the 611 down-regulated genes included those involved in redox regulation, light reactions of photosynthesis, and metabolism of lipids, amino acids, and cell wall. Overall, this pattern of changes in gene expression is characteristic of plants under stress. Thus, VFP3 likely plays an important role in controlling plant homeostasis.

High-resolution transcriptome analysis with long-read RNA sequencing.

RNA sequencing (RNA-seq) enables characterization and quantification of individual transcriptomes as well as detection of patterns of allelic expression and alternative splicing. Current RNA-seq protocols depend on high-throughput short-read sequencing of cDNA. However, as ongoing advances are rapidly yielding increasing read lengths, a technical hurdle remains in identifying the degree to which differences in read length influence various transcriptome analyses. In this study, we generated two paired-end RNA-seq datasets of differing read lengths (2×75 bp and 2×262 bp) for lymphoblastoid cell line GM12878 and compared the effect of read length on transcriptome analyses, including read-mapping performance, gene and transcript quantification, and detection of allele-specific expression (ASE) and allele-specific alternative splicing (ASAS) patterns. Our results indicate that, while the current long-read protocol is considerably more expensive than short-read sequencing, there are important benefits that can only be achieved with longer read length, including lower mapping bias and reduced ambiguity in assigning reads to genomic elements, such as mRNA transcript. We show that these benefits ultimately lead to improved detection of cis-acting regulatory and splicing variation effects within individuals.

Identifizierung von Promotoren in Corynebacterium glutamicum

The high-throughput sequencing of cDNA derived from whole transcript — omes (RNA-Seq) has revolutionized studies of bacterial transcription. We optimized and applied RNA-Seq for the genome-wide determination of transcription starts and their associated promoters using as an example of the biotechnologically relevant microorganism Corynebacterium glutamicum. By applying a method named 5′-RNA-Seq we were able to identify thousands of promoters in this organism.

Methods to determine the transcriptomes of trypanosomes in mixtures with mammalian cells: the effects of parasite purification and selective cDNA amplification.

Patterns of gene expression in cultured Trypanosoma brucei bloodstream and procyclic forms have been extensively characterized, and some comparisons have been made with trypanosomes grown to high parasitaemias in laboratory rodents. We do not know, however, to what extent these transcriptomes resemble those in infected Tsetse flies - or in humans or cattle, where parasitaemias are substantially lower. For clinical and field samples it is difficult to characterize parasite gene expression because of the large excess of host cell RNA. We have here examined two potential solutions to this problem for bloodstream form trypanosomes, assaying transcriptomes by high throughput cDNA sequencing (RNASeq). We first purified the parasites from blood of infected rats. We found that a red blood cell lysis procedure affected the transcriptome substantially more than purification using a DEAE cellulose column, but that too introduced significant distortions and variability. As an alternative, we specifically amplified parasite sequences from a mixture containing a 1000-fold excess of human RNA. We first purified polyadenylated RNA, then made trypanosome-specific cDNA by priming with a spliced leader primer. Finally, the cDNA was amplified using nested primers. The amplification procedure was able to produce samples in which 20% of sequence reads mapped to the trypanosome genome. Synthesis of the second cDNA strand with a spliced leader primer, followed by amplification, is sufficiently reproducible to allow comparison of different samples so long as they are all treated in the same way. However, SL priming distorted the abundances of the cDNA products and definitely cannot be used, by itself, to measure absolute mRNA levels. The amplification method might be suitable for clinical samples with low parasitaemias, and could also be adapted for other Kinetoplastids and to samples from infected vectors.

Improved structural annotation of protein-coding genes in the Meloidogyne hapla genome using RNA-Seq.

As high-throughput cDNA sequencing (RNA-Seq) is increasingly applied to hypothesis-driven biological studies, the prediction of protein coding genes based on these data are usurping strictly in silico approaches. Compared with computationally derived gene predictions, structural annotation is more accurate when based on biological evidence, particularly RNA-Seq data. Here, we refine the current genome annotation for the Meloidogyne hapla genome utilizing RNA-Seq data. Published structural annotation defines 14 420 protein-coding genes in the M. hapla genome. Of these, 25% (3751) were found to exhibit some incongruence with RNA-Seq data. Manual annotation enabled these discrepancies to be resolved. Our analysis revealed 544 new gene models that were missing from the prior annotation. Additionally, 1457 transcribed regions were newly identified on the ends of as-yet-unjoined contigs. We also searched for trans-spliced leaders, and based on RNA-Seq data, identified genes that appear to be trans-spliced. Four 22-bp trans-spliced leaders were identified using our pipeline, including the known trans-spliced leader, which is the M. hapla ortholog of SL1. In silico predictions of trans-splicing were validated by comparison with earlier results derived from an independent cDNA library constructed to capture trans-spliced transcripts. The new annotation, which we term HapPep5, is publically available at www.hapla.org.

Hyb: A bioinformatics pipeline for the analysis of CLASH (crosslinking, ligation and sequencing of hybrids) data

Associations between proteins and RNA–RNA duplexes are important in post-transcriptional regulation of gene expression. The CLASH (Cross-linking, Ligation and Sequencing of Hybrids) technique captures RNA–RNA interactions by physically joining two RNA molecules associated with a protein complex into a single chimeric RNA molecule. These events are relatively rare and considerable effort is needed to detect a small number of chimeric sequences amongst millions of non-chimeric cDNA reads resulting from a CLASH experiment. We present the “hyb” bioinformatics pipeline, which we developed to analyse high-throughput cDNA sequencing data from CLASH experiments. Although primarily designed for use with AGO CLASH data, hyb can also be used for the detection and annotation of chimeric reads in other high-throughput sequencing datasets. We examined the sensitivity and specificity of chimera detection in a test dataset using the BLAST, BLAST+, BLAT, pBLAT and Bowtie2 read alignment programs. We obtained the most reliable results in the shortest time using a combination of preprocessing with Flexbar and subsequent read-mapping using Bowtie2. The “hyb” software is distributed under the GNU GPL (General Public License) and can be downloaded from https://github.com/gkudla/hyb.

The Baboon Kidney Transcriptome: Analysis of Transcript Sequence, Splice Variants, and Abundance

The baboon is an invaluable model for the study of human health and disease, including many complex diseases of the kidney. Although scientists have made great progress in developing this animal as a model for numerous areas of biomedical research, genomic resources for the baboon, such as a quality annotated genome, are still lacking. To this end, we characterized the baboon kidney transcriptome using high-throughput cDNA sequencing (RNA-Seq) to identify genes, gene variants, single nucleotide polymorphisms (SNPs), insertion-deletion polymorphisms (InDels), cellular functions, and key pathways in the baboon kidney to provide a genomic resource for the baboon. Analysis of our sequencing data revealed 45,499 high-confidence SNPs and 29,813 InDels comparing baboon cDNA sequences with the human hg18 reference assembly and identified 35,900 cDNAs in the baboon kidney, including 35,150 transcripts representing 15,369 genic genes that are novel for the baboon. Gene ontology analysis of our sequencing dataset also identified numerous biological functions and canonical pathways that were significant in the baboon kidney, including a large number of metabolic pathways that support known functions of the kidney. The results presented in this study catalogues the transcribed mRNAs, noncoding RNAs, and hypothetical proteins in the baboon kidney and establishes a genomic resource for scientists using the baboon as an experimental model.

A comparison of methods for differential expression analysis of RNA-seq data

BackgroundFinding genes that are differentially expressed between conditions is an integral part of understanding the molecular basis of phenotypic variation. In the past decades, DNA microarrays have been used extensively to quantify the abundance of mRNA corresponding to different genes, and more recently high-throughput sequencing of cDNA (RNA-seq) has emerged as a powerful competitor. As the cost of sequencing decreases, it is conceivable that the use of RNA-seq for differential expression analysis will increase rapidly. To exploit the possibilities and address the challenges posed by this relatively new type of data, a number of software packages have been developed especially for differential expression analysis of RNA-seq data.ResultsWe conducted an extensive comparison of eleven methods for differential expression analysis of RNA-seq data. All methods are freely available within the R framework and take as input a matrix of counts, i.e. the number of reads mapping to each genomic feature of interest in each of a number of samples. We evaluate the methods based on both simulated data and real RNA-seq data.ConclusionsVery small sample sizes, which are still common in RNA-seq experiments, impose problems for all evaluated methods and any results obtained under such conditions should be interpreted with caution. For larger sample sizes, the methods combining a variance-stabilizing transformation with the ‘limma’ method for differential expression analysis perform well under many different conditions, as does the nonparametric SAMseq method.

Application of qRT-PCR and RNA-Seq analysis for the identification of housekeeping genes useful for normalization of gene expression values during Striga hermonthica development

Striga is a root parasitic weed that attacks many of the staple crops in Africa, India and Southeast Asia, inflicting tremendous losses in yield and for which there are few effective control measures. Studies of parasitic plant virulence and host resistance will be greatly facilitated by the recent emergence of genomic resources that include extensive transcriptome sequence datasets spanning all life stages of S. hermonthica. Functional characterization of Striga genes will require detailed analyses of gene expression patterns. Quantitative real-time PCR is a powerful tool for quantifying gene expression, but correct normalization of expression levels requires identification of control genes that have stable expression across tissues and life stages. Since no S. hermonthica housekeeping genes have been established for this purpose, we evaluated the suitability of six candidate housekeeping genes across key life stages of S. hermonthica from seed conditioning to flower initiation using qRT-PCR and high-throughput cDNA sequencing. Based on gene expression analysis by qRT-PCR and RNA-Seq across heterogeneous Striga life stages, we determined that using the combination of three genes, UBQ1, PP2A and TUB1 provides the best normalization for gene expression throughout the parasitic life cycle. The housekeeping genes characterized here provide robust standards that will facilitate powerful descriptions of parasite gene expression patterns.

Achieving high throughput sequencing of a cDNA library utilizing an alternative protocol for the bench top next-generation sequencing system

The development of next-generation sequencing (NGS) technologies has provided novel tools for genome analysis and expression profiling. A high throughput cDNA sequencing method using a bench top next-generation sequencing system, GS Junior, is now available. Here, we used an alternative protocol to the standard method for generating the cDNA library. This protocol can decrease the number of processing steps to manipulate RNA when constructing a cDNA library from an RNA sample, and does not require mRNA isolation from total RNA. Thus it can decrease the risk of RNA degradation and the cost for preparing a cDNA library. Also, the efficiency of sequencing data obtained with this approach is comparable to the standard method as verified by sequencing characteristics and expression levels of the reference gene glyceraldehyde-3-phosphate dehydrogenase (GAPDH).

Pervasive Initiation and 3′-End Formation of Poxvirus Postreplicative RNAs

Poxviruses are large DNA viruses that replicate within the cytoplasm and encode a complete transcription system, including a multisubunit RNA polymerase, stage-specific transcription factors, capping and methylating enzymes, and a poly(A) polymerase. Expression of the more than 200 open reading frames by vaccinia virus, the prototype poxvirus, is temporally regulated: early mRNAs are synthesized immediately after infection, whereas intermediate and late mRNAs are synthesized following genome replication. The postreplicative transcripts are heterogeneous in length and overlap the entire genome, which pose obstacles for high resolution mapping. We used tag-based methods in conjunction with high throughput cDNA sequencing to determine the precise 5'-capped and 3'-polyadenylated ends of postreplicative RNAs. Polymerase slippage during initiation of intermediate and late RNA synthesis results in a 5'-poly(A) leader that allowed the unambiguous identification of true transcription start sites. Ninety RNA start sites were located just upstream of intermediate and late open reading frames, but many more appeared anomalous, occurring within coding and non-coding regions, indicating pervasive transcription initiation. We confirmed the presence of functional promoter sequences upstream of representative anomalous start sites and demonstrated that alternative start sites within open reading frames could generate truncated isoforms of proteins. In an analogous manner, poly(A) sequences allowed accurate mapping of the numerous 3'-ends of postreplicative RNAs, which were preceded by a pyrimidine-rich sequence in the DNA coding strand. The distribution of postreplicative promoter sequences throughout the genome provides enormous transcriptional complexity, and the large number of previously unmapped RNAs may have novel functions.

Molecular and developmental contributions to divergent pigment patterns in marine and freshwater sticklebacks

Pigment pattern variation across species or populations offers a tractable framework in which to investigate the evolution of development. Juvenile threespine sticklebacks (Gasterosteus aculeatus) from marine and freshwater environments exhibit divergent pigment patterns that are associated with ecological differences. Juvenile marine sticklebacks have a silvery appearance, whereas sticklebacks from freshwater environments exhibit a pattern of vertical bars. We investigated both the developmental and molecular basis of this population-level variation in pigment pattern. Time course imaging during the transition from larval to juvenile stages revealed differences between marine and freshwater fish in spatial patterns of chromatophore differentiation as well as in pigment amount and dispersal. In freshwater fish, melanophores appear primarily within dark bars whereas iridophores appear within light bars. By contrast, in marine fish, these chromatophores are interspersed across the flank. In addition to spatially segregated chromatophore differentiation, pigment amount and dispersal within melanophores varies spatially across the flank of freshwater, but not marine fish. To gain insight into the molecular pathways that underlie the differences in pigment pattern development, we evaluated differential gene expression in the flanks of developing fish using high-throughput cDNA sequencing (RNA-seq) and quantitative PCR. We identified several genes that were differentially expressed across dark and light bars of freshwater fish, and between freshwater and marine fish. Together, these experiments begin to shed light on the process of pigment pattern evolution in sticklebacks.