A computational workflow to identify allele-specific expression and epigenetic modification in maize

Allele-specific expression refers to the preferential expression of one of the two alleles in a diploid genome, which has been thought largely attributable to the associated cis-element variation and allele-specific epigenetic modification patterns. Allele-specific expression may contribute to the heterosis (or hybrid vigor) effect in hybrid plants that are produced from crosses of closely-related species, subspecies and/or inbred lines. In this study, using Illumina high-throughput sequencing of maize transcriptomics, chromatic H3K27me3 histone modification and DNA methylation data, we developed a new computational framework to identify allele-specifically expressed genes by simultaneously tracking allele-specific gene expression patterns and the epigenetic modification landscape in the seedling tissues of hybrid maize. This approach relies on detecting nucleotide polymorphisms and any genomic structural variation between two parental genomes in order to distinguish paternally or maternally derived sequencing reads. This computational pipeline also incorporates a modified Chi-square test to statistically identify allele-specific gene expression and epigenetic modification based on the Poisson distribution.

Genome-wide mapping of polyadenylation sites in fission yeast reveals widespread alternative polyadenylation

Regulatory elements in the 3′ untranslated regions (UTRs) of eukaryotic mRNAs influence mRNA localization, translation, and stability. 3′-UTR length is determined by the location at which mRNAs are cleaved and polyadenylated. The use of alternative polyadenylation sites is common, and can be regulated in different situations. I present a new method to identify cleavage and polyadenylation sites (CSs) at the genome-wide level. The approach is strand-specific, avoids RNA enzymatic modification steps that can introduce sequence-specific biases, and uses unique molecular identifiers to ensure that all identified CS originates from individual RNA molecules. I applied this method to create the first comprehensive genome-wide map of polyadenylation sites of the fission yeast Schizosaccharomyces pombe, comprising the analysis of 2,021,000 individual mRNAs that defined 8,883 CSs. CSs were identified for 90% of coding genes and 50% of ncRNAs. Alternative polyadenylation was prevalent in both groups, with 41% and 45% of all detected genes, respectively, displaying more than one CS. The specificity of the cleavage reaction was gene-specific, resulting in highly variable levels of heterogeneity in 3′-UTR lengths. Finally, I show that for both coding and non-coding genes, the most common regulatory motif associated with CSs in fission yeast is the canonical human AAUAAA sequence.

Model selection as a tool for phylogeographic inference: an example from the willow Salix melanopsis

Phylogeographic inference has typically relied on analyses of data from one or a few genes to provide estimates of demography and population histories. While much has been learned from these studies, all phylogeographic analysis is conditioned on the data, and thus, inferences derived from data that represent a small sample of the genome are unavoidably tenuous. Here, we demonstrate one approach for moving beyond classic phylogeographic research. We use sequence capture probes and Illumina sequencing to generate data from >400 loci in order to infer the phylogeographic history of Salix melanopsis, a riparian willow with a disjunct distribution in coastal and the inland Pacific Northwest. We evaluate a priori phylogeographic hypotheses using coalescent models for parameter estimation, and the results support earlier findings that identified post-Pleistocene dispersal as the cause of the disjunction in S. melanopsis. We also conduct a series of model selection exercises using IMa2, Migrate-n and ∂a∂i. The resulting ranking of models indicates that refugial dynamics were complex, with multiple regions in the inland regions serving as the source for postglacial colonization. Our results demonstrate that new sources of data and new approaches to data analysis can rejuvenate phylogeographic research by allowing for the identification of complex models that enable researchers to both identify and estimate the most relevant parameters for a given system.

Rapid construction of parallel analysis of RNA end (PARE) libraries for Illumina sequencing

MicroRNAs (miRNAs) are ∼21nt small RNAs that pair to their target mRNAs and in many cases trigger cleavage, particularly in plants. Although many computational tools can predict miRNA:mRNA interactions, it remains critical to validate cleavage events, due to miRNA function in translational repression or due to high rates of false positives (over 90%) for unvalidated target predictions. A few years ago, three laboratories described similar methods to validate cleavage of miRNA targets by the cloning en masse of 5' ends of cleaved or uncapped mRNAs. To take advantage of the recent progress in high-throughput sequencing technology, we have devised an updated protocol to (1) enable much faster library preparation, and (2) reduce the cost by pooling indexed samples together for sequencing. Here we provide a step-by-step protocol for PARE library construction, starting from total RNA. This protocol has been successfully used in our laboratory to validate miRNA targets in a variety of plant species. We also provide advice for troubleshooting on some common issues.

Prediction of the protein components of a putative Calanus finmarchicus (Crustacea, Copepoda) circadian signaling system using a de novo assembled transcriptome

Diel vertical migration and seasonal diapause are critical life history events for the copepod Calanus finmarchicus. While much is known about these behaviors phenomenologically, little is known about their molecular underpinnings. Recent studies in insects suggest that some circadian genes/proteins also contribute to the establishment of seasonal diapause. Thus, it is possible that in Calanus these distinct timing regimes share some genetic components. To begin to address this possibility, we used the well-established Drosophila melanogaster circadian system as a reference for mining clock transcripts from a 200,000+ sequence Calanus transcriptome; the proteins encoded by the identified transcripts were also deduced and characterized. Sequences encoding homologs of the Drosophila core clock proteins CLOCK, CYCLE, PERIOD and TIMELESS were identified, as was one encoding CRYPTOCHROME 2, a core clock protein in ancestral insect systems, but absent in Drosophila. Calanus transcripts encoding proteins known to modulate the Drosophila core clock were also identified and characterized, e.g. CLOCKWORK ORANGE, DOUBLETIME, SHAGGY and VRILLE. Alignment and structural analyses of the deduced Calanus proteins with their Drosophila counterparts revealed extensive sequence conservation, particularly in functional domains. Interestingly, reverse BLAST analyses of these sequences against all arthropod proteins typically revealed non-Drosophila isoforms to be most similar to the Calanus queries. This, in combination with the presence of both CRYPTOCHROME 1 (a clock input pathway protein) and CRYPTOCHROME 2 in Calanus, suggests that the organization of the copepod circadian system is an ancestral one, more similar to that of insects like Danaus plexippus than to that of Drosophila.

Tissue-specific transcriptome profiling of Plutella xylostella third instar larval midgut

The larval midgut of diamondback moth, Plutella xylostella, is a dynamic tissue that interfaces with a diverse array of physiological and toxicological processes, including nutrient digestion and allocation, xenobiotic detoxification, innate and adaptive immune response, and pathogen defense. Despite its enormous agricultural importance, the genomic resources for P. xylostella are surprisingly scarce. In this study, a Bt resistant P. xylostella strain was subjected to the in-depth transcriptome analysis to identify genes and gene networks putatively involved in various physiological and toxicological processes in the P. xylostella larval midgut. Using Illumina deep sequencing, we obtained roughly 40 million reads containing approximately 3.6 gigabases of sequence data. De novo assembly generated 63,312 ESTs with an average read length of 416 bp, and approximately half of the P. xylostella sequences (45.4%, 28,768) showed similarity to the non-redundant database in GenBank with a cut-off E-value below 10(-5). Among them, 11,092 unigenes were assigned to one or multiple GO terms and 16,732 unigenes were assigned to 226 specific pathways. In-depth analysis identified genes putatively involved in insecticide resistance, nutrient digestion, and innate immune defense. Besides conventional detoxification enzymes and insecticide targets, novel genes, including 28 chymotrypsins and 53 ABC transporters, have been uncovered in the P. xylostella larval midgut transcriptome; which are potentially linked to the Bt toxicity and resistance. Furthermore, an unexpectedly high number of ESTs, including 46 serpins and 7 lysozymes, were predicted to be involved in the immune defense.As the first tissue-specific transcriptome analysis of P. xylostella, this study sheds light on the molecular understanding of insecticide resistance, especially Bt resistance in an agriculturally important insect pest, and lays the foundation for future functional genomics research. In addition, current sequencing effort greatly enriched the existing P. xylostella EST database, and makes RNAseq a viable option in the future genomic analysis.

Genome-wide mapping of cytosine methylation revealed dynamic DNA methylation patterns associated with genes and centromeres in rice

We conducted genome-wide mapping of cytosine methylation using methylcytosine immunoprecipitation combined with Illumina sequencing. The chromosomal distribution pattern of methylated DNA is similar to the heterochromatin distribution pattern on rice chromosomes. The DNA methylation patterns of rice genes are similar to those in Arabidopsis thaliana, including distinct methylation patterns asssociated with gene bodies and promoters. The DNA sequences in the core domains of rice Cen4, Cen5 and Cen8 showed elevated methylation levels compared with sequences in the pericentromeric regions. In addition, elevated methylation levels were associated with the DNA sequences in the CENH3-binding subdomains, compared with the sequences in the flanking H3 subdomains. In contrast, the centromeric domain of Cen11, which is composed exclusively of centromeric satellite DNA, is hypomethylated compared with the pericentromeric domains. Thus, the DNA sequences associated with functional centromeres can be either hypomethylated or hypermethylated. The methylation patterns of centromeric DNA appear to be correlated with the composition of the associated DNA sequences. We propose that both hypomethylation and hypermethylation of CENH3-associated DNA sequences can serve as epigenetic marks to distinguish where CENH3 deposition will occur within the surrounding H3 chromatin.