A draft genome of a field-collected Steinernema feltiae strain NW.

Advances in sequencing technologies have accelerated our understanding of the complex genetic network of organisms and genomic divergences that are linked to evolutionary processes. While many model organisms and laboratory strains have been sequenced, wild populations are underrepresented in the growing list of sequenced genomes. Here, we present a de novo assembly of Steinernema feltiae, strain NW, collected from a working agricultural field in south central Washington, USA. Leveraging Pacific Biosciences (PacBio) long reads, we sequenced strain NW to a high depth (99×). The resulting de novo assembly is significantly larger than the previous assembly generated from the laboratory strain SN, with a noticeable improvement in continuity and completeness. Comparative analysis of two assemblies revealed numerous single nucleotide polymorphisms (SNPs), breakpoints, and indels present between the two genomes. This alternative genome resource and annotation could benefit the research community to examine the genetic foundation of evolutionary processes as well as genomic variation among conspecific populations.Advances in sequencing technologies have accelerated our understanding of the complex genetic network of organisms and genomic divergences that are linked to evolutionary processes. While many model organisms and laboratory strains have been sequenced, wild populations are underrepresented in the growing list of sequenced genomes. Here, we present a de novo assembly of Steinernema feltiae, strain NW, collected from a working agricultural field in south central Washington, USA. Leveraging Pacific Biosciences (PacBio) long reads, we sequenced strain NW to a high depth (99×). The resulting de novo assembly is significantly larger than the previous assembly generated from the laboratory strain SN, with a noticeable improvement in continuity and completeness. Comparative analysis of two assemblies revealed numerous single nucleotide polymorphisms (SNPs), breakpoints, and indels present between the two genomes. This alternative genome resource and annotation could benefit the research community to examine the genetic foundation of evolutionary processes as well as genomic variation among conspecific populations.

Calcium is involved in the R Mc1 (blb)-mediated hypersensitive response against Meloidogyne chitwoodi in potato.

KEY MESSAGE: Functional characterization of the Columbia root-knot nematode resistance gene R Mc1 ( blb ) in potato revealed the R gene-mediated resistance is dependent on a hypersensitive response and involves calcium. The resistance (R) gene R Mc1(blb) confers resistance against the plant-parasitic nematode, Meloidogyne chitwoodi. Avirulent and virulent nematodes were used to functionally characterize the R Mc1(blb)-mediated resistance mechanism in potato (Solanum tuberosum). Histological observations indicated a hypersensitive response (HR) occurred during avirulent nematode infection. This was confirmed by quantifying reactive oxygen species activity in response to avirulent and virulent M. chitwoodi. To gain an insight into the signal transduction pathways mediating the R Mc1(blb)-induced HR, chemical inhibitors were utilized. Inhibiting Ca(2+) channels caused a significant reduction in electrolyte leakage, an indicator of cell death. Labeling with a Ca(2+)-sensitive dye revealed high Ca(2+) levels in the root cells surrounding avirulent nematodes. Furthermore, the calcium-dependent protein kinase (CDPK), StCDPK4 had a higher transcript level in R Mc1(blb) potato roots infected with avirulent nematodes in comparison to roots infected with virulent M. chitwoodi. The results of this study indicate Ca(2+) plays a role in the R Mc1(blb)-mediated resistance against M. chitwoodi in potato.

Broad Meloidogyne Resistance in Potato Based on RNA Interference of Effector Gene 16D10.

Root-knot nematodes (Meloidogyne spp.) are a significant problem in potato (Solanum tuberosum) production. There is no potato cultivar with Meloidogyne resistance, even though resistance genes have been identified in wild potato species and were introgressed into breeding lines. The objectives of this study were to generate stable transgenic potato lines in a cv. Russet Burbank background that carry an RNA interference (RNAi) transgene capable of silencing the 16D10 Meloidogyne effector gene, and test for resistance against some of the most important root-knot nematode species affecting potato, i.e., M. arenaria, M. chitwoodi, M. hapla, M. incognita, and M. javanica. At 35 days after inoculation (DAI), the number of egg masses per plant was significantly reduced by 65% to 97% (P < 0.05) in the RNAi line compared to wild type and empty vector controls. The largest reduction was observed in M. hapla, whereas the smallest reduction occurred in M. javanica. Likewise, the number of eggs per plant was significantly reduced by 66% to 87% in M. arenaria and M. hapla, respectively, compared to wild type and empty vector controls (P < 0.05). Plant-mediated RNAi silencing of the 16D10 effector gene resulted in significant resistance against all of the root-knot nematode species tested, whereas R Mc1(blb) , the only known Meloidogyne resistance gene in potato, did not have a broad resistance effect. Silencing of 16D10 did not interfere with the attraction of M. incognita second-stage juveniles to roots, nor did it reduce root invasion.

Nondestructive imaging of plant-parasitic nematode development and host response to nematode pathogenesis.

The secluded lifestyle of endoparasitic plant nematodes hampers progress toward a comprehensive understanding of plant-nematode interactions. A novel technique that enables nondestructive, long-term observations of a wide range of live nematodes in planta is presented here. As proof of principle, Pratylenchus penetrans, Heterodera schachtii, and Meloidogyne chitwoodi were labeled fluorescently with PKH26 and used to infect Arabidopsis thaliana grown in microscopy rhizosphere chambers. Nematode behavior, development, and morphology were observed for the full duration of each parasite's life cycle by confocal microscopy for up to 27 days after inoculation. PKH26 accumulated in intestinal lipid droplets and had no negative effect on nematode infectivity. This technique enabled visualization of Meloidogyne gall formation, nematode oogenesis, and nematode morphological features, such as the metacorpus, vulva, spicules, and cuticle. Additionally, microscopy rhizosphere chambers were used to characterize plant organelle dynamics during M. chitwoodi infection. Peroxisome abundance strongly increased in early giant cells but showed a marked decrease at later stages of feeding site development, which suggests a modulation of plant peroxisomes by root-knot nematodes during the infection process. Taken together, this technique facilitates studies aimed at deciphering plant-nematode interactions at the cellular and subcellular level and enables unprecedented insights into nematode behavior in planta.

RNA Interference of Effector Gene Mc16D10L Confers Resistance Against Meloidogyne chitwoodi in Arabidopsis and Potato.

Meloidogyne chitwoodi, a quarantine pathogen, is a significant problem in potato-producing areas worldwide. In spite of considerable genetic diversity in wild potato species, no commercial potato cultivars with resistance to M. chitwoodi are available. Nematode effector genes are essential for the molecular interactions between root-knot nematodes and their hosts. Stable transgenic lines of Arabidopsis and potato (Solanum tuberosum) with resistance against M. chitwoodi were developed. RNA interference (RNAi) construct pART27(16D10i-2) was introduced into Arabidopsis thaliana and potato to express double-stranded RNA complementary to the putative M. chitwoodi effector gene Mc16D10L. Plant-mediated RNAi led to a significant level of resistance against M. chitwoodi in Arabidopsis and potato. In transgenic Arabidopsis lines, the number of M. chitwoodi egg masses and eggs was reduced by up to 57 and 67% compared with empty vector controls, respectively. Similarly, in stable transgenic lines of potato, the number of M. chitwoodi egg masses and eggs was reduced by up to 71 and 63% compared with empty vector controls, respectively. The relative transcript level of Mc16D10L was reduced by up to 76% in M. chitwoodi eggs and infective second-stage juveniles that developed on transgenic pART27(16D10i-2) potato, suggesting that the RNAi effect is systemic and heritable in M. chitwoodi.

Nematode effector proteins: an emerging paradigm of parasitism.

Phytonematodes use a stylet and secreted effectors to modify host cells and ingest nutrients to support their growth and development. The molecular function of nematode effectors is currently the subject of intense investigation. In this review, we summarize our current understanding of nematode effectors, with a particular focus on proteinaceous stylet-secreted effectors of sedentary endoparasitic phytonematodes, for which a wealth of information has surfaced in the past 10 yr. We provide an update on the effector repertoires of several of the most economically important genera of phytonematodes and discuss current approaches to dissecting their function. Lastly, we highlight the latest breakthroughs in effector discovery that promise to shed new light on effector diversity and function across the phylum Nematoda.

Global epigenetic and transcriptional trends among two rice subspecies and their reciprocal hybrids.

The behavior of transcriptomes and epigenomes in hybrids of heterotic parents is of fundamental interest. Here, we report highly integrated maps of the epigenome, mRNA, and small RNA transcriptomes of two rice (Oryza sativa) subspecies and their reciprocal hybrids. We found that gene activity was correlated with DNA methylation and both active and repressive histone modifications in transcribed regions. Differential epigenetic modifications correlated with changes in transcript levels among hybrids and parental lines. Distinct patterns in gene expression and epigenetic modifications in reciprocal hybrids were observed. Through analyses of single nucleotide polymorphisms from our sequence data, we observed a high correlation of allelic bias of epigenetic modifications or gene expression in reciprocal hybrids with their differences in the parental lines. The abundance of distinct small RNA size classes differed between the parents, and more small RNAs were downregulated than upregulated in the reciprocal hybrids. Together, our data reveal a comprehensive overview of transcriptional and epigenetic trends in heterotic rice crosses and provide a useful resource for the rice community.

Major emerging problems with minor meloidogyne species.

Root-knot nematodes (Meloidogyne spp.) represent one of the most polyphagous genera of plant-parasitic nematodes. To date, close to 100 valid species are recognized. In contrast to the size of the genus, the majority of past research focused on a small number of species, i.e., the so-called 'major' species M. arenaria, M. hapla, M. incognita, and M. javanica. This review highlights recent work aimed at 'minor' root-knot nematodes: M. chitwoodi, M. fallax, M. minor, M. enterolobii (=M. mayaguensis), M. exigua, and M. paranaensis. Some of these species have been described only recently. After a brief profile of each species, identification methods and their application in Meloidogyne spp. are summarized. Intraspecific variation and its impact on plant resistance breeding are discussed and interactions between M. enterolobii and Fusarium solani are highlighted as an example of synergistic interactions with other plant pathogens. Future research on Meloidogyne spp. is not only shaped by recent breakthroughs such as completing the genome sequences of M. hapla and M. incognita, but is also influenced by changes in agriculture. Taken together, the aim of this review is to draw attention to previously neglected and newly described Meloidogyne spp. that are developing into major problems for agriculture in tropical and temperate climates.

Relationship Between Climatic Factors and Distribution of Pratylenchus spp. in the Dryland Wheat-Production Areas of Eastern Washington.

Field surveys were conducted by collecting soil samples to estimate nematode densities in soil from winter wheat, spring wheat, spring barley, and spring legumes (lentil, chickpea, and pea) fields during 2010 and 2011. Pratylenchus spp. were observed in 60% of sampled fields. However, nematodes were detected in nearly all of the survey fields in high numbers where crops were grown every year. To identify climatic variables associated with density of Pratylenchus spp. in soil, correlation and regression analyses were performed using climate data of survey sites from 1979 to 2010. Fifty-seven climate variables were significantly correlated with densities of Pratylenchus spp. All precipitation variables were significantly positively correlated with nematode abundance. Summer maximum air temperature was negatively correlated and winter minimum air temperature was positively correlated with nematode densities. In addition, both years' nematode densities were significantly correlated with cropping intensity. Five multivariate regression models for 2010 and seven models for 2011 nematode abundance levels were developed. The majority of the climate variables selected in the models were related to precipitation. Knowledge of root-lesion nematode distribution in the dryland region of eastern Washington and associated climate variables may be helpful to determine risk and apply management practices to minimize crop damage.

Dynamic landscapes of four histone modifications during deetiolation in Arabidopsis.

Although landscapes of several histone marks are now available for Arabidopsis thaliana and Oryza sativa, such profiles remain static and do not provide information about dynamic changes of plant epigenomes in response to developmental or environmental cues. Here, we analyzed the effects of light on four histone modifications (acetylation and trimethylation of lysines 9 and 27 on histone H3: H3K9ac, H3K9me3, H3K27ac, and H3K27me3, respectively). Our genome-wide profiling of H3K9ac and H3K27ac revealed that these modifications are nontransposable element gene-specific. By contrast, we found that H3K9me3 and H3K27me3 target nontransposable element genes, but also intergenic regions and transposable elements. Specific light conditions affected the number of modified regions as well as the overall correlation strength between the presence of specific modifications and transcription. Furthermore, we observed that acetylation marks not only ELONGATED HYPOCOTYL5 and HY5-HOMOLOG upon deetiolation, but also their downstream targets. We found that the activation of photosynthetic genes correlates with dynamic acetylation changes in response to light, while H3K27ac and H3K27me3 potentially contribute to light regulation of the gibberellin metabolism. Thus, this work provides a dynamic portrait of the variations in histone modifications in response to the plant's changing light environment and strengthens the concept that histone modifications represent an additional layer of control for light-regulated genes involved in photomorphogenesis.

Genome-wide profiling of histone H3 lysine 9 acetylation and dimethylation in Arabidopsis reveals correlation between multiple histone marks and gene expression.

Lysine residue 9 of histone H3 can either be acetylated or mono-, di-, or tri-methylated. These epigenetic states have a diverse impact on regulating gene transcriptional activity and chromatin organization. H3K9ac is invariably correlated with transcriptional activation, whereas H3K9me2 has been reported to be mainly located in constitutive heterochromatin in Arabidopsis. Here, we present epigenetic landscapes for histone H3 lysine 9 acetylation (H3K9ac) and dimethylation (H3K9me2) in Arabidopsis seedlings. The results show that H3K9ac targeted 5,206 non-transposable element (non-TE) genes and 321 transposable elements (TEs), whereas H3K9me2 targeted 2,281 TEs and 1,112 non-TE genes. H3K9ac was biased towards the 5' end of genes and peaked at the ATG position, while H3K9me2 tended to span the entire gene body. H3K9ac correlated with high gene expression, while H3K9me2 correlated with low expression. Analyses of H3K9ac and H3K9me2 with the available datasets of H3K27me3 and DNA methylation revealed a correlation between the occurrence of multiple epigenetic modifications and gene expression. Genes with H3K9ac alone were actively transcribed, while genes that were also modified by either H3K27me3 or DNA methylation showed a lower expression level, suggesting that a combination of repressive marks weakened the positive regulatory effect of H3K9ac. Furthermore, we observed a significant increase of the H3K9ac modification level of selected target genes in hda19 (histone deacetylase 19) mutant seedlings, which indicated that HDA19 plays an important role in regulating the level of H3K9ac and thereby influencing the transcriptional activity in young seedlings.

Sequence mining and transcript profiling to explore cyst nematode parasitism.

BACKGROUND: Cyst nematodes are devastating plant parasites that become sedentary within plant roots and induce the transformation of normal plant cells into elaborate feeding cells with the help of secreted effectors, the parasitism proteins. These proteins are the translation products of parasitism genes and are secreted molecular tools that allow cyst nematodes to infect plants. RESULTS: We present here the expression patterns of all previously described parasitism genes of the soybean cyst nematode, Heterodera glycines, in all major life stages except the adult male. These insights were gained by analyzing our gene expression dataset from experiments using the Affymetrix Soybean Genome Array GeneChip, which contains probeset sequences for 6,860 genes derived from preparasitic and parasitic H. glycines life stages. Targeting the identification of additional H. glycines parasitism-associated genes, we isolated 633 genes encoding secretory proteins using algorithms to predict secretory signal peptides. Furthermore, because some of the known H. glycines parasitism proteins have strongest similarity to proteins of plants and microbes, we searched for predicted protein sequences that showed their highest similarities to plant or microbial proteins and identified 156 H. glycines genes, some of which also contained a signal peptide. Analyses of the expression profiles of these genes allowed the formulation of hypotheses about potential roles in parasitism. This is the first study combining sequence analyses of a substantial EST dataset with microarray expression data of all major life stages (except adult males) for the identification and characterization of putative parasitism-associated proteins in any parasitic nematode. CONCLUSION: We have established an expression atlas for all known H. glycines parasitism genes. Furthermore, in an effort to identify additional H. glycines genes with putative functions in parasitism, we have reduced the currently known 6,860 H. glycines genes to a pool of 788 most promising candidate genes (including known parasitism genes) and documented their expression profiles. Using our approach to pre-select genes likely involved in parasitism now allows detailed functional analyses in a manner not feasible for larger numbers of genes. The generation of the candidate pool described here is an important enabling advance because it will significantly facilitate the unraveling of fascinating plant-animal interactions and deliver knowledge that can be transferred to other pathogen-host systems. Ultimately, the exploration of true parasitism genes verified from the gene pool delineated here will identify weaknesses in the nematode life cycle that can be exploited by novel anti-nematode efforts.

GmEREBP1 is a transcription factor activating defense genes in soybean and Arabidopsis.

Ethylene-responsive element-binding proteins (EREBPs) are plant-specific transcription factors, many of which have been linked to plant defense responses. Conserved EREBP domains bind to the GCC box, a promoter element found in pathogenesis-related (PR) genes. We previously identified an EREBP gene from soybean (GmEREBP1) whose transcript abundance decreased in soybean cyst-nematode-infected roots of a susceptible cultivar, whereas it increased in abundance in infected roots of a resistant cultivar. Here, we report further characterization of this gene. Transient expression analyses showed that GmEREBP1 is localized to the plant nucleus and functions as a transcriptional activator in soybean leaves. Transgenic soybean plants expressing GmEREBP1 activated the expression of the ethylene (ET)-responsive gene PR2 and the ET- and jasmonic acid (JA)-responsive gene PR3, and the salicylic acid (SA)-responsive gene PR1 but not the SA-responsive PR5. Similarly, transgenic Arabidopsis plants expressing GmEREBP1 showed elevated mRNA abundance of the ET-regulated gene PR3 and the ET- and JA-regulated defense-related gene PDF1.2 but not the ET-regulated GST2, and the SA-regulated gene PR1 but not the SA-regulated PR2 and PR5. Transgenic soybean and Arabidopsis plants inoculated with cyst nematodes did not display a significantly altered susceptibility to nematode infection. These results collectively show that GmEREBP1 functions as a transacting inducer of defense gene expression in both soybean and Arabidopsis and mediates the expression of both ET- and JA- and SA-regulated defense-related genes in these plant species.

Active uptake of cyst nematode parasitism proteins into the plant cell nucleus.

Cyst nematodes produce parasitism proteins that contain putative nuclear localisation signals (NLSs) and, therefore, are predicted to be imported into the nucleus of the host plant cell. The in planta localisation patterns of eight soybean cyst nematode (Heterodera glycines) parasitism proteins with putative NLSs were determined by producing these proteins as translational fusions with the GFP and GUS reporter proteins. Two parasitism proteins were found to be imported into the nuclei of onion epidermal cells as well as Arabidopsis protoplasts. One of these two parasitism proteins was further transported into the nucleoli. Mutations introduced into the NLS domains of these two proteins abolished nuclear import and caused a cytoplasmic accumulation. Furthermore, we observed active nuclear uptake for three additional parasitism proteins, however, only when these proteins were synthesised as truncated forms. Two of these proteins were further transported into nucleoli. We hypothesise that nuclear uptake and nucleolar localisation are important mechanisms for H. glycines to modulate the nuclear biology of parasitised cells of its host plant.

Mitochondrial genome plasticity among species of the nematode genus Meloidogyne (Nematoda: Tylenchina).

The mitochondrial (mt) genomes of the plant-parasitic root-knot nematodes Meloidogyne arenaria, Meloidogyne enterolobii and Meloidogyne javanica were sequenced and compared with those of three other root-knot nematode species in order to explore the mt genome plasticity within Meloidogyne. The mt genomes of M. arenaria, M. enterolobii and M. javanica are circular, with an estimated size of 18.8, 18.9 and 19.6 kb, respectively. Compared to other nematodes these mt genomes are larger, due to the presence of large non-coding regions. The mt genome architecture within the genus Meloidogyne varied in the position of trn genes and in the position, length and nucleotide composition of non-coding regions. These variations were observed independent of the species' natural environments or reproductive modes. M. enterolobii showed three main non-coding regions whereas Meloidogyne chitwoodi, Meloidogyne incognita, M. javanica and M. arenaria had two non-coding regions, and Meloidogyne graminicola had a unique large non-coding region interrupted by two trn genes. trn genes were positioned in different regions of the mt genomes in M. chitwoodi, M. enterolobii and M. graminicola, whereas the trn gene order was identical between M. arenaria, M. incognita and M. javanica. Importantly, M. graminicola had extra copies of trnV and trnS2. High divergence levels between the two copies of each trn might indicate duplication events followed by random loss and mutations in the anticodon. Tree-based methods based on amino acid sequences of 12 mt protein-coding genes support the monophyly for the tropical and mitotic parthenogenetic species, M. arenaria, M. enterolobii, M. incognita and M. javanica and for a clade that includes the meiotic parthenogenetic species, M. chitwoodi and M. graminicola. A comparison of the mt genome architecture in plant-parasitic nematodes and phylogenetic analyses support that Pratylenchus is the most recent ancestor of root-knot nematodes.

A Meloidogyne incognita effector is imported into the nucleus and exhibits transcriptional activation activity in planta.

Root-knot nematodes are sedentary biotrophic endoparasites that maintain a complex interaction with their host plants. Nematode effector proteins are synthesized in the oesophageal glands of nematodes and secreted into plant tissue through a needle-like stylet. Effectors characterized to date have been shown to mediate processes essential for nematode pathogenesis. To gain an insight into their site of action and putative function, the subcellular localization of 13 previously isolated Meloidogyne incognita effectors was determined. Translational fusions were created between effectors and EGFP-GUS (enhanced green fluorescent protein-ß-glucuronidase) reporter genes, which were transiently expressed in tobacco leaf cells. The majority of effectors localized to the cytoplasm, with one effector, 7H08, imported into the nuclei of plant cells. Deletion analysis revealed that the nuclear localization of 7H08 was mediated by two novel independent nuclear localization domains. As a result of the nuclear localization of the effector, 7H08 was tested for the ability to activate gene transcription. 7H08 was found to activate the expression of reporter genes in both yeast and plant systems. This is the first report of a plant-parasitic nematode effector with transcriptional activation activity.

Functional characterization of nematode effectors in plants.

Secreted effectors represent the molecular interface between the nematode and its host plant. Studies that aimed at deciphering molecular plant-nematode interactions are hampered by technical hurdles that prevent the generation of transgenic nematodes. However, RNA interference (RNAi) has proven to be a valuable tool to specifically knock-down nematode effector genes, both ex planta and in planta. Plant-mediated RNAi of nematode genes not only facilitates functional characterization of effectors but also lends itself to a novel control strategy against plant-parasitic nematodes. Here, we describe currently used methods to silence genes in plant-parasitic cyst and root-knot nematodes.

Mitochondrial genomes of Meloidogyne chitwoodi and M. incognita (Nematoda: Tylenchina): comparative analysis, gene order and phylogenetic relationships with other nematodes.

Root-knot nematodes (Meloidogyne spp.) are among the most important plant pathogens. In this study, the mitochondrial (mt) genomes of the root-knot nematodes, M. chitwoodi and M. incognita were sequenced. PCR analyses suggest that both mt genomes are circular, with an estimated size of 19.7 and 18.6-19.1kb, respectively. The mt genomes each contain a large non-coding region with tandem repeats and the control region. The mt gene arrangement of M. chitwoodi and M. incognita is unlike that of other nematodes. Sequence alignments of the two Meloidogyne mt genomes showed three translocations; two in transfer RNAs and one in cox2. Compared with other nematode mt genomes, the gene arrangement of M. chitwoodi and M. incognita was most similar to Pratylenchus vulnus. Phylogenetic analyses (Maximum Likelihood and Bayesian inference) were conducted using 78 complete mt genomes of diverse nematode species. Analyses based on nucleotides and amino acids of the 12 protein-coding mt genes showed strong support for the monophyly of class Chromadorea, but only amino acid-based analyses supported the monophyly of class Enoplea. The suborder Spirurina was not monophyletic in any of the phylogenetic analyses, contradicting the Clade III model, which groups Ascaridomorpha, Spiruromorpha and Oxyuridomorpha based on the small subunit ribosomal RNA gene. Importantly, comparisons of mt gene arrangement and tree-based methods placed Meloidogyne as sister taxa of Pratylenchus, a migratory plant endoparasitic nematode, and not with the sedentary endoparasitic Heterodera. Thus, comparative analyses of mt genomes suggest that sedentary endoparasitism in Meloidogyne and Heterodera is based on convergent evolution.

The epigenome and plant development.

The epigenomic regulation of chromatin structure and genome stability is essential for the interpretation of genetic information and ultimately the determination of phenotype. High-resolution maps of plant epigenomes have been obtained through a combination of chromatin technologies and genomic tiling microarrays and through high-throughput sequencing-based approaches. The transcriptomic activity of a plant at a certain stage of development is controlled by genome-wide combinatorial interactions of epigenetic modifications. Tissue- or environment-specific epigenomes are established during plant development. Epigenomic reprogramming triggered by the activation and movement of small RNAs is important for plant gametogenesis. Genome-wide loss of DNA methylation in the endosperm and the accompanying endosperm-specific gene expression during seed development provide a genomic insight into epigenetic regulation of gene imprinting in plants. Global changes of histone modifications during plant responses to different light environments play an important regulatory role in a sophisticated light-regulated transcriptional network. Epigenomic natural variation that developed during evolution is important for phenotypic diversity and can potentially contribute to the molecular mechanisms of complex biological phenomena such as heterosis in plants.

Next-generation sequencing reveals complex relationships between the epigenome and transcriptome in maize.

Epigenetic modifications and small RNAs play an important role in gene regulation. Here, we discuss results of our Solexa/Illumina 1G sequencing-based survey of DNA methylation, activating and repressive histone modifications, small RNAs and mRNA in the maize genome. We analyze tissue-specific epigenetic patterns, discuss antagonistic relationships between repressive epigenetic marks and highlight synergistic relationships between activating histone modifications. We discuss our observation that small RNAs show a tissue-specific distribution in maize. Whereas 24-nucleotide long small interfering RNAs (siRNAs) accumulated preferentially in shoots, 21-nucleotide long micro RNAs (miRNAs) were the most abundant group in roots, which follows the transcript level of mop1. Furthermore, we discuss the possibility that a novel class of 22-nucleotide siRNAs might originate from long double-stranded RNAs in an RNA-dependent RNA polymerase (RdRP)-independent manner. This supports the intriguing possibility that maize possesses at least two distinct pathways to generate siRNAs, one of which relies on RdRP and a second one that might be RdRP-independent.