Spatial gradients in cell wall composition and transcriptional profiles along elongating maize internodes.

BACKGROUND: The elongating maize internode represents a useful system for following development of cell walls in vegetative cells in the Poaceae family. Elongating internodes can be divided into four developmental zones, namely the basal intercalary meristem, above which are found the elongation, transition and maturation zones. Cells in the basal meristem and elongation zones contain mainly primary walls, while secondary cell wall deposition accelerates in the transition zone and predominates in the maturation zone. RESULTS: The major wall components cellulose, lignin and glucuronoarabinoxylan (GAX) increased without any abrupt changes across the elongation, transition and maturation zones, although GAX appeared to increase more between the elongation and transition zones. Microarray analyses show that transcript abundance of key glycosyl transferase genes known to be involved in wall synthesis or re-modelling did not match the increases in cellulose, GAX and lignin. Rather, transcript levels of many of these genes were low in the meristematic and elongation zones, quickly increased to maximal levels in the transition zone and lower sections of the maturation zone, and generally decreased in the upper maturation zone sections. Genes with transcript profiles showing this pattern included secondary cell wall CesA genes, GT43 genes, some ß-expansins, UDP-Xylose synthase and UDP-Glucose pyrophosphorylase, some xyloglucan endotransglycosylases/hydrolases, genes involved in monolignol biosynthesis, and NAM and MYB transcription factor genes. CONCLUSIONS: The data indicated that the enzymic products of genes involved in cell wall synthesis and modification remain active right along the maturation zone of elongating maize internodes, despite the fact that corresponding transcript levels peak earlier, near or in the transition zone.

The maize methylome influences mRNA splice sites and reveals widespread paramutation-like switches guided by small RNA.

The maize genome, with its large complement of transposons and repeats, is a paradigm for the study of epigenetic mechanisms such as paramutation and imprinting. Here, we present the genome-wide map of cytosine methylation for two maize inbred lines, B73 and Mo17. CG (65%) and CHG (50%) methylation (where H = A, C, or T) is highest in transposons, while CHH (5%) methylation is likely guided by 24-nt, but not 21-nt, small interfering RNAs (siRNAs). Correlations with methylation patterns suggest that CG methylation in exons (8%) may deter insertion of Mutator transposon insertion, while CHG methylation at splice acceptor sites may inhibit RNA splicing. Using the methylation map as a guide, we used low-coverage sequencing to show that parental methylation differences are inherited by recombinant inbred lines. However, frequent methylation switches, guided by siRNA, persist for up to eight generations, suggesting that epigenetic inheritance resembling paramutation is much more common than previously supposed. The methylation map will provide an invaluable resource for epigenetic studies in maize.

Complex regulation by Apetala2 domain-containing transcription factors revealed through analysis of the stress-responsive TdCor410b promoter from durum wheat.

Expression of the wheat dehydrin gene Cor410b is induced several fold above its non-stressed levels upon exposure to stresses such as cold, drought and wounding. Deletion analysis of the TdCor410b promoter revealed a single functional C-repeat (CRT) element. Seven transcription factors (TFs) were shown to bind to this CRT element using yeast one-hybrid screens of wheat and barley cDNA libraries, of which only one belonged to the DREB class of TFs. The remaining six encoded ethylene response factors (ERFs) belong to three separate subfamilies. Analysis of binding selectivity of these TFs indicated that all seven could bind to the CRT element (GCCGAC), and that three of the six ERFs could bind both to the CRT element and the ethylene-responsive GCC-box (GCCGCC). The TaERF4 subfamily members specifically bound the CRT element, and did not bind either the GCC-box or DRE element (ACCGAC). Molecular modeling and site-directed mutagenesis identified a single residue Pro42 in the Apetala2 (AP2) domain of TaERF4-like proteins that is conserved in monocotyledonous plants and is responsible for the recognition selectivity of this subfamily. We suggest that both DREB and ERF proteins regulate expression of the Cor410b gene through a single, critical CRT element. Members of the TaERF4 subfamily are specific, positive regulators of Cor410b gene expression.

Endo-(1,4)-ß-glucanase gene families in the grasses: temporal and spatial co-transcription of orthologous genes.

BACKGROUND: Endo-(1,4)-ß-glucanase (cellulase) glycosyl hydrolase GH9 enzymes have been implicated in several aspects of cell wall metabolism in higher plants, including cellulose biosynthesis and degradation, modification of other wall polysaccharides that contain contiguous (1,4)-ß-glucosyl residues, and wall loosening during cell elongation. RESULTS: The endo-(1,4)-ß-glucanase gene families from barley (Hordeum vulgare), maize (Zea mays), sorghum (Sorghum bicolor), rice (Oryza sativa) and Brachypodium (Brachypodium distachyon) range in size from 23 to 29 members. Phylogenetic analyses show variations in clade structure between the grasses and Arabidopsis, and indicate differential gene loss and gain during evolution. Map positions and comparative studies of gene structures allow orthologous genes in the five species to be identified and synteny between the grasses is found to be high. It is also possible to differentiate between homoeologues resulting from ancient polyploidizations of the maize genome. Transcript analyses using microarray, massively parallel signature sequencing and quantitative PCR data for barley, rice and maize indicate that certain members of the endo-(1,4)-ß-glucanase gene family are transcribed across a wide range of tissues, while others are specifically transcribed in particular tissues. There are strong correlations between transcript levels of several members of the endo-(1,4)-ß-glucanase family and the data suggest that evolutionary conservation of transcription exists between orthologues across the grass family. There are also strong correlations between certain members of the endo-(1,4)-ß-glucanase family and other genes known to be involved in cell wall loosening and cell expansion, such as expansins and xyloglucan endotransglycosylases. CONCLUSIONS: The identification of these groups of genes will now allow us to test hypotheses regarding their functions and joint participation in wall synthesis, re-modelling and degradation, together with their potential role in lignocellulose conversion during biofuel production from grasses and cereal crop residues.

Cell wall modifications in maize pulvini in response to gravitational stress.

Changes in cell wall polysaccharides, transcript abundance, metabolite profiles, and hormone concentrations were monitored in the upper and lower regions of maize (Zea mays) pulvini in response to gravistimulation, during which maize plants placed in a horizontal position returned to the vertical orientation. Heteroxylan levels increased in the lower regions of the pulvini, together with lignin, but xyloglucans and heteromannan contents decreased. The degree of substitution of heteroxylan with arabinofuranosyl residues decreased in the lower pulvini, which exhibited increased mechanical strength as the plants returned to the vertical position. Few or no changes in noncellulosic wall polysaccharides could be detected on the upper side of the pulvinus, and crystalline cellulose content remained essentially constant in both the upper and lower pulvinus. Microarray analyses showed that spatial and temporal changes in transcript profiles were consistent with the changes in wall composition that were observed in the lower regions of the pulvinus. In addition, the microarray analyses indicated that metabolic pathways leading to the biosynthesis of phytohormones were differentially activated in the upper and lower regions of the pulvinus in response to gravistimulation. Metabolite profiles and measured hormone concentrations were consistent with the microarray data, insofar as auxin, physiologically active gibberellic acid, and metabolites potentially involved in lignin biosynthesis increased in the elongating cells of the lower pulvinus.

A customized gene expression microarray reveals that the brittle stem phenotype fs2 of barley is attributable to a retroelement in the HvCesA4 cellulose synthase gene.

The barley (Hordeum vulgare) brittle stem mutants, fs2, designated X054 and M245, have reduced levels of crystalline cellulose compared with their parental lines Ohichi and Shiroseto. A custom-designed microarray, based on long oligonucleotide technology and including genes involved in cell wall metabolism, revealed that transcript levels of very few genes were altered in the elongation zone of stem internodes, but these included a marked decrease in mRNA for the HvCesA4 cellulose synthase gene of both mutants. In contrast, the abundance of several hundred transcripts changed in the upper, maturation zones of stem internodes, which presumably reflected pleiotropic responses to a weakened cell wall that resulted from the primary genetic lesion. Sequencing of the HvCesA4 genes revealed the presence of a 964-bp solo long terminal repeat of a Copia-like retroelement in the first intron of the HvCesA4 genes of both mutant lines. The retroelement appears to interfere with transcription of the HvCesA4 gene or with processing of the mRNA, and this is likely to account for the lower crystalline cellulose content and lower stem strength of the mutants. The HvCesA4 gene maps to a position on chromosome 1H of barley that coincides with the previously reported position of fs2.

Artificial chromosome formation in maize (Zea mays L.).

We report on the construction of maize minichromosomes using shuttle vectors harboring native centromeric segments, origins of replication, selectable marker genes, and telomeric repeats. These vectors were introduced into scutellar cells of maize immature embryos by microprojectile bombardment. Several independent transformation events were identified containing minichromosomes in addition to the normal diploid complement of 20 maize chromosomes. Immunostaining indicated that the minichromosomes recruited centromeric protein C, which is a specific component of the centromere/kinetochore complex. Minichromosomes were estimated to be 15-30 Mb in size based on cytological measurements. Fluorescent in situ hybridization (FISH) showed that minichromosomes contain the centromeric, telomeric, and exogenous unique marker sequences interspersed with maize retrotransposons. Minichromosomes were detected for at least a year in actively dividing callus cultures, providing evidence for their stability through numerous cell cycles. Plants were regenerated and minichromosomes were detected in root tips, providing confirmation of their normal replication and transmission during mitosis and through organogenesis. Assembly of maize artificial chromosomes may provide a tool to study centromere function and a foundation for developing new high capacity vectors for plant functional genomics and breeding.

The role of regulatory genes during maize domestication: evidence from nucleotide polymorphism and gene expression.

We investigated DNA sequence variation in 72 candidate genes in maize landraces and the wild ancestor of maize, teosinte. The candidate genes were chosen because they exhibit very low sequence diversity among maize inbreds and have sequence homology to known regulatory genes. We observed signatures of selection in 17 candidate genes, indicating that they were potential targets of artificial selection during domestication. In addition, 21 candidate genes were identified as potential targets of natural selection in teosinte. A comparison of the proportion of selected genes between our regulatory genes and genes unfiltered for their potential function (but also with very low sequence diversity among maize inbreds) provided some weak evidence that regulatory genes are overrepresented among selected genes. We detected no significant association between the positions of genes identified as potential targets of selection during domestication and quantitative trait loci (QTL) responsible for maize domestication traits. However, a subset of these genes, those identified by sequence homology as kinase/phosphatase genes, significantly cluster with the domestication QTL. We also analyzed expression profiles of genes in distinct maize tissues and observed that domestication genes are expressed on average at a significantly higher level than neutral genes in reproductive organs, including kernels.

Artificial trans-acting siRNAs confer consistent and effective gene silencing.

Manipulating gene expression is critical to exploring gene function and a useful tool for altering commercial traits. Techniques such as hairpin-based RNA interference, virus-induced gene silencing, and artificial microRNAs take advantage of endogenous posttranscriptional gene silencing pathways to block translation of designated transcripts. Here we present a novel gene silencing method utilizing artificial trans-acting small interfering RNAs in Arabidopsis (Arabidopsis thaliana). Replacing the endogenous small interfering RNAs encoded in the TAS1c gene with sequences from the FAD2 gene silenced FAD2 activity to levels comparable to the fad2-1 null allele in nearly all transgenic events. Interestingly, exchanging the endogenous miR173 target sequence in TAS1c with an miR167 target sequence led to variable, inefficient silencing of FAD2, suggesting a specific requirement for the miR173 trigger for production of small interfering RNAs from the TAS1c locus.

Whole genome scan detects an allelic variant of fad2 associated with increased oleic acid levels in maize.

We used whole genome scan association mapping to identify loci with major effect on oleic acid content in maize kernels. Single nucleotide polymorphism haplotypes at 8,590 loci were tested for association with oleic acid content in 553 maize inbreds. A single locus with major effect on oleic acid was mapped between 380 and 384 cM in the IBM2 neighbors genetic map on chromosome 4 and confirmed in a biparental population. A fatty acid desaturase, fad2, identified approximately 2 kb from the associated genetic marker, is the most likely candidate gene responsible for the differences in the phenotype. The fad2 alleles with high- and low-oleic acid content were sequenced and allelic differences in fad2 RNA level in developing embryos was investigated. We propose that a non-conservative amino acid polymorphism near the active site of fad2 contributes to the effect on oleic acid content. This is the first report of the use of a high resolution whole genome scan association mapping where a putative gene responsible for a quantitative trait was identified in plants.

Conserved noncoding genomic sequences associated with a flowering-time quantitative trait locus in maize.

Flowering time is a fundamental trait of maize adaptation to different agricultural environments. Although a large body of information is available on the map position of quantitative trait loci for flowering time, little is known about the molecular basis of quantitative trait loci. Through positional cloning and association mapping, we resolved the major flowering-time quantitative trait locus, Vegetative to generative transition 1 (Vgt1), to an approximately 2-kb noncoding region positioned 70 kb upstream of an Ap2-like transcription factor that we have shown to be involved in flowering-time control. Vgt1 functions as a cis-acting regulatory element as indicated by the correlation of the Vgt1 alleles with the transcript expression levels of the downstream gene. Additionally, within Vgt1, we identified evolutionarily conserved noncoding sequences across the maize-sorghum-rice lineages. Our results support the notion that changes in distant cis-acting regulatory regions are a key component of plant genetic adaptation throughout breeding and evolution.

Identification of transposons, retroelements, and a gene family predominantly expressed in floral tissues in chromosome 3DS of the hexaploid wheat progenitor Aegilops tauschii.

A multigene family expressed during early floral development was identified on the short arm of wheat chromosome 3D in the region of the Ph2 locus, a locus controlling homoeologous chromosome pairing in allohexaploid wheat. Physical, genetic and molecular characterisation of the Wheat Meiosis 1 (WM1) gene family identified seven members that localised within a region of 173-kb. WM1 gene family members were sequenced and they encode mainly type Ia plasma membrane-anchored leucine rich repeat-like receptor proteins. In situ expression profiling suggests the gene family is predominantly expressed in floral tissue. In addition to the WM1 gene family, a number of other genes, gene fragments and pseudogenes were identified. It has been predicted that there is approximately one gene every 19-kb and that this region of the wheat genome contains 23 repetitive elements including BARE-1 and Wis2-1 like sequences. Nearly 50% of the repetitive elements identified were similar to known transposons from the CACTA superfamily. Ty1-copia, Ty3-gypsy and Athila LTR retroelements were also prevalent within the region. The WM1 gene cluster is present on 3DS and on barley 3HS but missing from the A and B genomes of hexaploid wheat. This suggests either recent generation of the cluster or specific deletion of the cluster during wheat polyploidisation. The evolutionary significance of the cluster, its possible roles in disease response or floral and early meiotic development and its location at or near the Ph2 locus are discussed.

Evolution of DNA sequence nonhomologies among maize inbreds.

Allelic chromosomal regions totaling more than 2.8 Mb and located on maize (Zea mays) chromosomes 1L, 2S, 7L, and 9S have been sequenced and compared over distances of 100 to 350 kb between the two maize inbred lines Mo17 and B73. The alleles contain extended regions of nonhomology. On average, more than 50% of the compared sequence is noncolinear, mainly because of the insertion of large numbers of long terminal repeat (LTR)-retrotransposons. Only 27 LTR-retroelements are shared between alleles, whereas 62 are allele specific. The insertion of LTR-retrotransposons into the maize genome is statistically more recent for nonshared than shared ones. Most surprisingly, more than one-third of the genes (27/72) are absent in one of the inbreds at the loci examined. Such nonshared genes usually appear to be truncated and form clusters in which they are oriented in the same direction. However, the nonshared genome segments are gene-poor, relative to regions shared by both inbreds, with up to 12-fold difference in gene density. By contrast, miniature inverted terminal repeats (MITEs) occur at a similar frequency in the shared and nonshared fractions. Many times, MITES are present in an identical position in both LTRs of a retroelement, indicating that their insertion occurred before the replication of the retroelement in question. Maize ESTs and/or maize massively parallel signature sequencing tags were identified for the majority of the nonshared genes or homologs of them. In contrast with shared genes, which are usually conserved in gene order and location relative to rice (Oryza sativa), nonshared genes violate the maize colinearity with rice. Based on this, insertion by a yet unknown mechanism, rather than deletion events, seems to be the origin of the nonshared genes. The intergenic space between conserved genes is enlarged up to sixfold in maize compared with rice. Frequently, retroelement insertions create a different sequence environment adjacent to conserved genes.

Linkage disequilibrium and sequence diversity in a 500-kbp region around the adh1 locus in elite maize germplasm.

Linkage disequilibrium (LD) at the adh locus was examined in two sets of maize inbreds. A set of 32 was chosen to represent most of the genetic diversity in the cultivated North American elite maize breeding pool. A second set of 192 inbreds was chosen to sample more deeply the two major heterotic groups in elite maize germplasm. Analysis of several loci in the vicinity of the adh gene shows that LD as measured by D' and r2 extends greater than 500 kbp in this germplasm. The presence of this exceptionally long segment of high LD may be suggestive of selection acting on one of the genes in the vicinity of adh1 or of a locally reduced rate of recombination.

Long-range patterns of diversity and linkage disequilibrium surrounding the maize Y1 gene are indicative of an asymmetric selective sweep.

Both yellow and white corn occurs among ancestral open pollinated varieties. More recently, breeders have selected yellow endosperm variants of maize over ancestral white phenotypes for their increased nutritional value resulting from the up-regulation of the Y1 phytoene synthase gene product in endosperm tissue. As a result, diversity within yellow maize lines at the Y1 gene is dramatically decreased as compared to white corn. We analyzed patterns of sequence diversity and linkage disequilibrium in nine low copy regions located at varying distances from the Y1 gene, including a homolog of the barley Mlo gene. Patterns consistent with a selective sweep, such as significant associations of informative single-nucleotide polymorphisms with endosperm color phenotype, linkage disequilibrium, and significantly reduced diversity within the yellow endosperm haplotypes, were observed up to 600 kb downstream of Y1, whereas the upstream region showed a more rapid recovery. The starch branching enzyme 1 (sbe1) gene is the first region downstream of Y1 that does not have a highly conserved haplotype in the yellow endosperm germplasm.

Anchoring 9,371 maize expressed sequence tagged unigenes to the bacterial artificial chromosome contig map by two-dimensional overgo hybridization.

Our goal is to construct a robust physical map for maize (Zea mays) comprehensively integrated with the genetic map. We have used a two-dimensional 24 x 24 overgo pooling strategy to anchor maize expressed sequence tagged (EST) unigenes to 165,888 bacterial artificial chromosomes (BACs) on high-density filters. A set of 70,716 public maize ESTs seeded derivation of 10,723 EST unigene assemblies. From these assemblies, 10,642 overgo sequences of 40 bp were applied as hybridization probes. BAC addresses were obtained for 9,371 overgo probes, representing an 88% success rate. More than 96% of the successful overgo probes identified two or more BACs, while 5% identified more than 50 BACs. The majority of BACs identified (79%) were hybridized with one or two overgos. A small number of BACs hybridized with eight or more overgos, suggesting that these BACs must be gene rich. Approximately 5,670 overgos identified BACs assembled within one contig, indicating that these probes are highly locus specific. A total of 1,795 megabases (Mb; 87%) of the total 2,050 Mb in BAC contigs were associated with one or more overgos, which are serving as sequence-tagged sites for single nucleotide polymorphism development. Overgo density ranged from less than one overgo per megabase to greater than 20 overgos per megabase. The majority of contigs (52%) hit by overgos contained three to nine overgos per megabase. Analysis of approximately 1,022 Mb of genetically anchored BAC contigs indicates that 9,003 of the total 13,900 overgo-contig sites are genetically anchored. Our results indicate overgos are a powerful approach for generating gene-specific hybridization probes that are facilitating the assembly of an integrated genetic and physical map for maize.

Genes associated with the end of dormancy in grapes.

A grape bud EST library was constructed and 4270 ESTs sequenced. The library clones were arrayed for the purpose of investigating the level of gene expression over time, particularly leading up to the buds' release from dormancy. The arrays were hybridized with P(33)-labeled probes produced from samples of buds collected at weekly intervals. These probes covered the time from 9 weeks prior to bud burst until just after the emergence of the shoots. Expression patterns from these genes have been examined. It was found that 74% of the genes in the data set were homologous to known proteins. Genes were then assigned to functional categories according to their primary BLAST match. Of these 13% were involved with photosynthesis, 13% with disease resistance and defense, 5% energy, 12% metabolism, 20% protein production and processing, 25% cell structure and plant growth and the remaining 12% were unclassified The expression pattern of a selection of "candidate" genes retrieved from literature previously reporting an association with dormancy changes was assessed. On closer examination most of these genes relate to the oxidative processes and stress responses within the cell. The results of this study show that even in the dormant state, gene expression in the buds is high.

Monitoring genome-wide changes in gene expression in response to endogenous cytokinin reveals targets in Arabidopsis thaliana.

Cytokinins have been implicated in developmental and growth processes in plants including cell division, chloroplast biogenesis, shoot meristem initiation and senescence. The regulation of these processes requires changes in cytokinin-responsive gene expression. Here, we induced the expression of a bacterial isopentenyl transferase gene, IPT, in transgenic Arabidopsis thaliana seedlings to study the regulation of genome-wide gene expression in response to endogenous cytokinin. Using MPSS (massively parallel signature sequencing) we identified 823 and 917 genes that were up- and downregulated, respectively, following 24 h of IPT induction. When comparing the response to cytokinin after 6 and 24 h, we identified different clusters of genes showing a similar course of regulation. Our study provides researchers with the opportunity to rapidly assess whether genes of interest are regulated by cytokinins.

Genome-wide gene expression profiling in Arabidopsis thaliana reveals new targets of abscisic acid and largely impaired gene regulation in the abi1-1 mutant.

The phytohormone abscisic acid (ABA) plays important regulatory roles in many plant developmental processes including seed dormancy, germination, growth, and stomatal movements. These physiological responses to ABA are in large part brought about by changes in gene expression. To study genome-wide ABA-responsive gene expression we applied massively parallel signature sequencing (MPSS) to samples from Arabidopsis thaliana wildtype (WT) and abi1-1 mutant seedlings. We identified 1354 genes that are either up- or downregulated following ABA treatment of WT seedlings. Among these ABA-responsive genes, many encode signal transduction components. In addition, we identified novel ABA-responsive gene families including those encoding ribosomal proteins and proteins involved in regulated proteolysis. In the ABA-insensitive mutant abi1-1, ABA regulation of about 84.5% and 6.9% of the identified genes was impaired or strongly diminished, respectively; however, 8.6% of the genes remained appropriately regulated. Compared to other methods of gene expression analysis, the high sensitivity and specificity of MPSS allowed us to identify a large number of ABA-responsive genes in WT Arabidopsis thaliana. The database given in our supplementary material (http://jcs.biologists.org/supplemental) provides researchers with the opportunity to rapidly assess whether genes of interest may be regulated by ABA. Regulation of the majority of the genes by ABA was impaired in the ABA-insensitive mutant abi1-1. However, a subset of genes continued to be appropriately regulated by ABA, which suggests the presence of at least two ABA signaling pathways, only one of which is blocked in abi1-1.

SNP frequency, haplotype structure and linkage disequilibrium in elite maize inbred lines.

BACKGROUND: Recent studies of ancestral maize populations indicate that linkage disequilibrium tends to dissipate rapidly, sometimes within 100 bp. We set out to examine the linkage disequilibrium and diversity in maize elite inbred lines, which have been subject to population bottlenecks and intense selection by breeders. Such population events are expected to increase the amount of linkage disequilibrium, but reduce diversity. The results of this study will inform the design of genetic association studies. RESULTS: We examined the frequency and distribution of DNA polymorphisms at 18 maize genes in 36 maize inbreds, chosen to represent most of the genetic diversity in U.S. elite maize breeding pool. The frequency of nucleotide changes is high, on average one polymorphism per 31 bp in non-coding regions and 1 polymorphism per 124 bp in coding regions. Insertions and deletions are frequent in non-coding regions (1 per 85 bp), but rare in coding regions. A small number (2-8) of distinct and highly diverse haplotypes can be distinguished at all loci examined. Within genes, SNP loci comprising the haplotypes are in linkage disequilibrium with each other. CONCLUSIONS: No decline of linkage disequilibrium within a few hundred base pairs was found in the elite maize germplasm. This finding, as well as the small number of haplotypes, relative to neutral expectation, is consistent with the effects of breeding-induced bottlenecks and selection on the elite germplasm pool. The genetic distance between haplotypes is large, indicative of an ancient gene pool and of possible interspecific hybridization events in maize ancestry.