Ripening of tomato fruits is triggered by the plant hormone ethylene, but its effect is restricted by an unknown developmental cue to mature fruits containing viable seeds. To determine whether this cue involves epigenetic remodeling, we expose tomatoes to the methyltransferase inhibitor 5-azacytidine and find that they ripen prematurely. We performed whole-genome bisulfite sequencing on fruit in four stages of development, from immature to ripe. We identified 52,095 differentially methylated regions (representing 1% of the genome) in the 90% of the genome covered by our analysis. Furthermore, binding sites for RIN, one of the main ripening transcription factors, are frequently localized in the demethylated regions of the promoters of numerous ripening genes, and binding occurs in concert with demethylation. Our data show that the epigenome is not static during development and may have been selected to ensure the fidelity of developmental processes such as ripening. Crop-improvement strategies could benefit by taking into account not only DNA sequence variation among plant lines, but also the information encoded in the epigenome.

INTRODUCTIONConventional Illumina RNA-Seq does not have the resolution to decode the complex eukaryote transcriptome due to the lack of RNA polarity information. Strand-specific RNA sequencing (ssRNA-Seq) can overcome these limitations and as such is better suited for genome annotation, de novo transcriptome assembly, and accurate digital gene expression analysis. This protocol describes a simple and robust method to generate ssRNA-Seq libraries for the Illumina sequencing platform. It has significantly increased the throughput to 96 libraries in a two-day preparation while simultaneously lowering the reagent costs to below ten dollars per library. It is compatible with both single-read and paired-end multiplex sequencing and, most importantly, its data can also be used with existing conventional RNA-Seq data. This is a significant advantage, because it enables researchers to switch to ssRNA-Seq even if a large amount of data has already been generated by the nonstrand specific methods.

Molecular evolutionary rate variation in Gossypium (cotton) was characterized using sequence data for 48 nuclear genes from both genomes of allotetraploid cotton, models of its diploid progenitors, and an outgroup. Substitution rates varied widely among the 48 genes, with silent and replacement substitution levels varying from 0.018 to 0.162 and from 0.000 to 0.073, respectively, in comparisons between orthologous Gossypium and outgroup sequences. However, about 90% of the genes had silent substitution rates spanning a more narrow threefold range. Because there was no evidence of rate heterogeneity among lineages for any gene and because rates were highly correlated in independent tests, evolutionary rate is inferred to be a property of each gene or its genetic milieu rather than the clade to which it belongs. Evidence from approximately 200,000 nucleotides (40,000 per genome) suggests that polyploidy in Gossypium led to a modest enhancement in rates of nucleotide substitution. Phylogenetic analysis for each gene yielded the topology expected from organismal history, indicating an absence of gene conversion or recombination among homoeologs subsequent to allopolyploid formation. Using the mean synonymous substitution rate calculated across the 48 genes, allopolyploid cotton is estimated to have formed circa 1.5 million years ago (MYA), after divergence of the diploid progenitors about 6.7 MYA.

The golden apple snail (Pomacea canaliculata) is a fresh water snail listed among the top-100 worst invasive species, worldwide and a noted agricultural and quarantine pest that causes great economic losses. It is characterized by fast growth, strong stress tolerance, a high reproduction rate, and adaptation to a broad range of environments.Here, we used long-read sequencing to produce a 440-Mb high-quality chromosome-level assembly for the P. canaliculata genome. In total, 50 Mb (11.4%) repeat sequences and 21,533 gene models were identified in the genome. The major findings of this study include the recent explosion of DNA/hAT-Charlie transposable elements (TEs), the expansion of the P450 gene family and the constitution of the cellular homeostasis system, which contributes to ecological plasticity in stress adaptation. In addition, the high transcriptional levels of perivitellin genes in the ovary and albumen gland promote the function of nutrient supply and defence ability in eggs. Furthermore, the gut metagenome also contains diverse genes for food digestion and xenobiotic degradation.

Allopolyploidization has been a driving force in plant evolution. Formation of common wheat (Triticum aestivum L.) represents a classic example of successful speciation via allopolyploidy. Nevertheless, the immediate chromosomal consequences of allopolyploidization in wheat remain largely unexplored. We report here an in-depth investigation on transgenerational chromosomal variation in resynthesized allohexaploid wheats that are identical in genome constitution to common wheat. We deployed sequential FISH, genomic in situ hybridization (GISH), and homeolog-specific pyrosequencing, which enabled unequivocal identification of each of the 21 homologous chromosome pairs in each of >1,000 individual plants from 16 independent lines. We report that wholechromosome aneuploidy occurred ubiquitously in early generations (from selfed generation S 1 to >S 20 ) of wheat allohexaploidy although at highly variable frequencies (20-100%). In contrast, other types of gross structural variations were scant. Aneuploidy included an unexpected hidden type, which had a euploid chromosome number of 2n = 42 but with simultaneous loss and gain of nonhomeologous chromosomes. Of the three constituent subgenomes, B showed the most lability for aneuploidy, followed by A, but the recently added D subgenome was largely stable in most of the studied lines. Chromosome loss and gain were also unequal across the 21 homologous chromosome pairs. Pedigree analysis showed no evidence for progressive karyotype stabilization even with multigenerational selection for euploidy. Profiling of two traits directly related to reproductive fitness showed that although pollen viability was generally reduced by aneuploidy, the adverse effect of aneuploidy on seed-set is dependent on both aneuploidy type and synthetic line.chromosome dynamics | hidden aneuploidy | synthetic wheat | wheat evolution H exaploid common wheat (Triticum aestivum L.) is a major food crop with international significance, the evolution of which is characterized by two sequential allopolyploidization events: one leading to formation of allotetraploid wheat (T. turgidum L.) and the other to allohexaploid wheat (T. aestivum) (1, 2). Despite decades of research, the mechanisms by which the initial allopolyploid individuals became stabilized, established, and accumulate to successful speciation remains largely unknown in this important crop. In theory, chromosome-level perturbation should be among the first manifestations of nascent allopolyploidization. Indeed, two recent molecular cytogenetic studies, in resynthesized allotetraploid Brassica napus lines (3) and young natural allotetraploid Tragopogon miscellus populations (4), respectively, have provided unique insights into the chromosomal dynamics associated with nascent allotetraploidy. Being at the resolution of individual chromosomes, these studies have documented a surprisingly high incidence of both structural and numerical changes in nascent allotetraploid plants (3, 4). It was found that early generations of resynthesized allotetrap...

Enhancers are important regulators of gene expression in eukaryotes. Enhancers function independently of their distance and orientation to the promoters of target genes. Thus, enhancers have been difficult to identify. Only a few enhancers, especially distant intergenic enhancers, have been identified in plants. We developed an enhancer prediction system based exclusively on the DNase I hypersensitive sites (DHSs) in the Arabidopsis thaliana genome. A set of 10,044 DHSs located in intergenic regions, which are away from any gene promoters, were predicted to be putative enhancers. We examined the functions of 14 predicted enhancers using the b-glucuronidase gene reporter. Ten of the 14 (71%) candidates were validated by the reporter assay. We also designed 10 constructs using intergenic sequences that are not associated with DHSs, and none of these constructs showed enhancer activities in reporter assays. In addition, the tissue specificity of the putative enhancers can be precisely predicted based on DNase I hypersensitivity data sets developed from different plant tissues. These results suggest that the open chromatin signature-based enhancer prediction system developed in Arabidopsis may serve as a universal system for enhancer identification in plants.

Hybridization between different species plays an important role in plant genome evolution, as well as is a widely used approach for crop improvement. McClintock has predicted that plant wide hybridization constitutes a "genomic shock" whereby cryptic transposable elements may be activated. However, direct experimental evidence showing a causal relationship between plant wide hybridization and transposon mobilization has not yet been reported. The miniature-Ping (mPing) is a recently isolated active miniature inverted-repeat transposable element transposon from rice, which is mobilized by tissue culture and gamma-ray irradiation. We show herein that mPing, together with its putative transposase-encoding partner, Pong, is mobilized in three homologous recombinant inbred lines (RILs), derived from hybridization between rice (cultivar Matsumae) and wild rice (Zizania latifolia Griseb.), harboring introgressed genomic DNA from wild rice. In contrast, both elements remain immobile in two lines sharing the same parentage to the RILs but possessing no introgressed DNA. Thus, we have presented direct evidence that is consistent with McClintock's insight by demonstrating a causal link between wide hybridization and transposon mobilization in rice. In addition, we report an atypical behavior of mPing/Pong mobilization in these lines, i.e., the exclusive absence of footprints after excision.

Cytosine methylation at CG sites ( m CG) plays critical roles in development, epigenetic inheritance, and genome stability in mammals and plants. In the dicot model plant Arabidopsis thaliana, methyltransferase 1 (MET1), a principal CG methylase, functions to maintain m CG during DNA replication, with its null mutation resulting in global hypomethylation and pleiotropic developmental defects. Null mutation of a critical CG methylase has not been characterized at a whole-genome level in other higher eukaryotes, leaving the generality of the Arabidopsis findings largely speculative. Rice is a model plant of monocots, to which many of our important crops belong. Here we have characterized a null mutant of OsMet1-2, the major CG methylase in rice. We found that seeds homozygous for OsMet1-2 gene mutation (OsMET1-2 −/− ), which directly segregated from normal heterozygote plants (OsMET1-2 +/− ), were seriously maldeveloped, and all germinated seedlings underwent swift necrotic death. Compared with wild type, genome-wide loss of m CG occurred in the mutant methylome, which was accompanied by a plethora of quantitative molecular phenotypes including dysregulated expression of diverse protein-coding genes, activation and repression of transposable elements, and altered small RNA profiles. Our results have revealed conservation but also distinct functional differences in CG methylases between rice and Arabidopsis.Oryza sativa L. | monocotyledons C ytosine methylation is an evolutionarily conserved epigenetic modification across biological kingdoms (1-3). However, fundamental differences exist for this epigenetic mark between animals and plants in many aspects (1-3). For example, in mammalian somatic cells, methylated cytosines ( m Cs) occur almost exclusively in a CG sequence context, and CG methylation ( m CG) is maintained during DNA replication by a single CG methylase, DNA methyltransferase 1 (DNMT1). Cytosine methylation in plants is more complex, occurs in all sequence contexts (CG, CHG, and CHH; where H is A, T, or C), and is established and maintained by multiple enzymes that have distinct but also overlapping functions (2, 4-6). Nevertheless, m CG stands out as the most pervasive type (2). In the dicot model plant Arabidopsis thaliana, m CG is maintained by a major CG methylase, namely DNA methyltransferase 1 (MET1), the homolog of mammalian DNMT1 (2, 7-11). Null mutation of MET1 (the homozygous met1 mutant) virtually eliminates genome-wide m CG to 1.7% that of wild type (WT), that is, from 24.6% in WT to 0.42% in mutant (12, 13). Homozygous null met1 mutant can be obtained from selfed progenies of heterozygotes (MET1 +/− ) but was often at much lower frequencies than expected by Mendelian segregation (11). Homozygous Arabidopsis met1 plants show pleiotropic developmental abnormalities of variable penetrance, which aggravate in frequency and severity with inbreeding (11). Although conditional knockout DNMT1 mutants were established in human colorectal carcinoma cell lines (14) and loss-of-function DNA methyltran...