Targeted analysis of orthologous phytochrome A regions of the sorghum, maize, and rice genomes using comparative gene-island sequencing

CNMS: The preferred genic markers for comparative genomic, molecular phylogenetic, functional genetic diversity and differential gene regulatory expression analyses in chickpea

The intra/inter-genomic comparative mapping-based phylogenetic footprinting identified 5 paralogous and 656 orthologous genome-wide CNMS markers in the upstream sequences of chickpea genes. These CNMS markers revealed a high-degree of gene-based syntenic relationship between chickpea and Medicago genomes while minimum between chickpea and Vitis genomes. The time of divergence and duplication estimated using CNMS markers highlight the expected phylogenetic relationships between chickpea and six dicot (legume) species as well as occurrence of ancient genome (approximately 53 Mya) with small-scale recent segmental (approximately 10 Mya) duplication events in chickpea. A wider level of functional molecular diversity (14 to 88 percent) and admixed population genetic structure was detected among desi, kabuli and wild genotypes by genic CNMS markers at a genome-wide scale suggesting their utility in large-scale genetic analysis in chickpea. The subfunctionalization at the cis-regulatory element region and TFBS (transcription factor binding site) motif levels in the upstream sequences of CNMS marker-associated orthologous genes than the paralogues was predominant. Functional constraint might have considerable effect on these CNMScontaining regulatory elements controlling consistent orthologous gene expression in dicots. A rapid subfunctionalization based on diverge differential expression of paralogous CNMS marker-associated genes particularly those that underwent recent small-scale segmental duplication events in chickpea was apparent. The differential regulation of expression and subfunctionalization potential of Ultra CNMS marker-associated genes suggest their utility in deciphering the complex gene regulatory function as well as identification and targeted mapping of potential genes/QTLs governing vital agronomic traits in chickpea. The gene-based CNMS markers with desirable inherent genetic attributes like higher degree of comparative genome mapping, functional genetic diversity and differential gene regulatory expression potential can significantly propel the genomics-assisted chickpea crop improvement.

Genome-wide conserved non-coding microsatellite (CNMS) marker-based integrative genetical genomics for quantitative dissection of seed weight in chickpea

Phylogenetic footprinting identified 666 genome-wide paralogous and orthologous CNMS (conserved non-coding microsatellite) markers from 5'-untranslated and regulatory regions (URRs) of 603 protein-coding chickpea genes. The (CT)n and (GA)n CNMS carrying CTRMCAMV35S and GAGA8BKN3 regulatory elements, respectively, are abundant in the chickpea genome. The mapped genic CNMS markers with robust amplification efficiencies (94.7%) detected higher intraspecific polymorphic potential (37.6%) among genotypes, implying their immense utility in chickpea breeding and genetic analyses. Seventeen differentially expressed CNMS marker-associated genes showing strong preferential and seed tissue/developmental stage-specific expression in contrasting genotypes were selected to narrow down the gene targets underlying seed weight quantitative trait loci (QTLs)/eQTLs (expression QTLs) through integrative genetical genomics. The integration of transcript profiling with seed weight QTL/eQTL mapping, molecular haplotyping, and association analyses identified potential molecular tags (GAGA8BKN3 and RAV1AAT regulatory elements and alleles/haplotypes) in the LOB-domain-containing protein- and KANADI protein-encoding transcription factor genes controlling the cis-regulated expression for seed weight in the chickpea. This emphasizes the potential of CNMS marker-based integrative genetical genomics for the quantitative genetic dissection of complex seed weight in chickpea.

A gene regulatory network model for floral transition of the shoot apex in maize and its dynamic modeling

The transition from the vegetative to reproductive development is a critical event in the plant life cycle. The accurate prediction of flowering time in elite germplasm is important for decisions in maize breeding programs and best agronomic practices. The understanding of the genetic control of flowering time in maize has significantly advanced in the past decade. Through comparative genomics, mutant analysis, genetic analysis and QTL cloning, and transgenic approaches, more than 30 flowering time candidate genes in maize have been revealed and the relationships among these genes have been partially uncovered. Based on the knowledge of the flowering time candidate genes, a conceptual gene regulatory network model for the genetic control of flowering time in maize is proposed. To demonstrate the potential of the proposed gene regulatory network model, a first attempt was made to develop a dynamic gene network model to predict flowering time of maize genotypes varying for specific genes. The dynamic gene network model is composed of four genes and was built on the basis of gene expression dynamics of the two late flowering id1 and dlf1 mutants, the early flowering landrace Gaspe Flint and the temperate inbred B73. The model was evaluated against the phenotypic data of the id1 dlf1 double mutant and the ZMM4 overexpressed transgenic lines. The model provides a working example that leverages knowledge from model organisms for the utilization of maize genomic information to predict a whole plant trait phenotype, flowering time, of maize genotypes.

Role of conserved non-coding regulatory elements in LMW glutenin gene expression

Transcriptional regulation of LMW glutenin genes were investigated in-silico, using publicly available gene sequences and expression data. Genes were grouped into different LMW glutenin types and their promoter profiles were determined using cis-acting regulatory elements databases and published results. The various cis-acting elements belong to some conserved non-coding regulatory regions (CREs) and might act in two different ways. There are elements, such as GCN4 motifs found in the long endosperm box that could serve as key factors in tissue-specific expression. Some other elements, such as the AACA/TA motifs or the individual prolamin box variants, might modulate the level of expression. Based on the promoter sequences and expression characteristic LMW glutenin genes might be transcribed following two different mechanisms. Most of the s- and i-type genes show a continuously increasing expression pattern. The m-type genes, however, demonstrate normal distribution in their expression profiles. Differences observed in their expression could be related to the differences found in their promoter sequences. Polymorphisms in the number and combination of cis-acting elements in their promoter regions can be of crucial importance in the diverse levels of production of single LMW glutenin gene types.

Construction of a watermelon BAC library and identification of SSRs anchored to melon or Arabidopsis genomes

A bacterial artificial chromosome (BAC) library was constructed for watermelon (Citrullus lanatus (Thunb.) Matsum. & Nakai var. lanatus) with an average insert-size of 106 kb, providing 21 haploid genome equivalents. The library was used to identify BAC clones that are anchored to probes evenly distributed on the genomes of melon or Arabidopsis. Twenty eight probes (representing 66% of the tested probes) from melon and 30 probes (65%) from Arabidopsis identified positive BAC clones. Two methods were implemented to identify SSRs from the positively hybridizing BAC clones. First, analysis of BAC end sequences revealed 37 SSRs. For the second method, pooled DNA of BACs identified by the melon probes was used to develop a shotgun library. The library was then screened with synthetic SSR oligonucleotides by hybridization. Sequence analysis of positively hybridizing shotgun clones revealed 142 different SSRs. Thirty eight SSRs were characterized using three watermelon cultivars, five plant introduction (PI) accessions of C. lanatus var lanatus and four PIs of C. lanatus var citroides. Of these, 36 (95%) were found to be polymorphic with up to six alleles per marker. Polymorphism information content values for polymorphic markers varied between 0.22 and 0.79 with an average of 0.53. The methods described herein will be valuable for the construction of a watermelon linkage map with SSRs evenly distributed on its genome that is anchored to the genomes of melon and Arabidopsis.

Comparative structural and functional characterization of sorghum and maize duplications containing orthologous myb transcription regulators of 3-deoxyflavonoid biosynthesis

Sequence characterization of the genomic region of sorghum yellow seed 1 shows the presence of two genes that are arranged in a head to tail orientation. The two duplicated gene copies, y1 and y2 are separated by a 9.084 kbp intergenic region, which is largely composed of highly repetitive sequences. The y1 is the functional copy, while the y2 may represent a pseudogene; there are several sequence indels and rearrangements within the putative coding region of y2. The y1 gene encodes a R2R3 type of Myb domain protein that regulates the expression of chalcone synthase, chalcone isomerase and dihydroflavonol reductase genes required for the biosynthesis of 3-deoxyflavonoids. Expression of y1 can be observed throughout the plant and it represents a combination of expression patterns produced by different alleles of the maize p1. Comparative sequence analysis within the coding regions and flanking sequences of y1, y2 and their maize and teosinte orthologs show local rearrangements and insertions that may have created modified regulatory regions. These micro-colinearity modifications possibly are responsible for differential patterns of expression in maize and sorghum floral and vegetative tissues. Phylogenetic analysis indicates that sorghum y1 and y2 sequences may have arisen by gene duplication mechanisms and represent an evolutionarily parallel event to the duplication of maize p2 and p1 genes.

Comprehensive molecular cytogenetic analysis of sorghum genome architecture: distribution of euchromatin, heterochromatin, genes and recombination in comparison to rice

Cytogenetic maps of sorghum chromosomes 3-7, 9, and 10 were constructed on the basis of the fluorescence in situ hybridization (FISH) of approximately 18-30 BAC probes mapped across each of these chromosomes. Distal regions of euchromatin and pericentromeric regions of heterochromatin were delimited for all 10 sorghum chromosomes and their DNA content quantified. Euchromatic DNA spans approximately 50% of the sorghum genome, ranging from approximately 60% of chromosome 1 (SBI-01) to approximately 33% of chromosome 7 (SBI-07). This portion of the sorghum genome is predicted to encode approximately 70% of the sorghum genes ( approximately 1 gene model/12.3 kbp), assuming that rice and sorghum encode a similar number of genes. Heterochromatin spans approximately 411 Mbp of the sorghum genome, a region characterized by a approximately 34-fold lower rate of recombination and approximately 3-fold lower gene density compared to euchromatic DNA. The sorghum and rice genomes exhibit a high degree of macrocolinearity; however, the sorghum genome is approximately 2-fold larger than the rice genome. The distal euchromatic regions of sorghum chromosomes 3-7 and 10 are approximately 1.8-fold larger overall and exhibit an approximately 1.5-fold lower average rate of recombination than the colinear regions of the homeologous rice chromosomes. By contrast, the pericentromeric heterochromatic regions of these chromosomes are on average approximately 3.6-fold larger in sorghum and recombination is suppressed approximately 15-fold compared to the colinear regions of rice chromosomes.

Structure and expression of maize phytochrome family homeologs

To begin the study of phytochrome signaling in maize, we have cloned and characterized the phytochrome gene family from the inbred B73. Through DNA gel blot analysis of maize genomic DNA and BAC library screens, we show that the PhyA, PhyB, and PhyC genes are each duplicated once in the genome of maize. Each gene pair was positioned to homeologous regions of the genome using recombinant inbred mapping populations. These results strongly suggest that the duplication of the phytochrome gene family in maize arose as a consequence of an ancient tetraploidization in the maize ancestral lineage. Furthermore, sequencing of Phy genes directly from BAC clones indicates that there are six functional phytochrome genes in maize. Through Northern gel blot analysis and a semiquantitative reverse transcriptase polymerase chain reaction assay, we determined that all six phytochrome genes are transcribed in several seedling tissues. However, expression from PhyA1, PhyB1, and PhyC1 predominate in all seedling tissues examined. Dark-grown seedlings express higher levels of PhyA and PhyB than do light-grown plants but PhyC genes are expressed at similar levels under light and dark growth conditions. These results are discussed in relation to phytochrome gene regulation in model eudicots and monocots and in light of current genome sequencing efforts in maize.

Phylogenetic analysis of 5'-noncoding regions from the ABA-responsive rab16/17 gene family of sorghum, maize and rice provides insight into the composition, organization and function of cis-regulatory modules

Phylogenetic analysis of sequences from gene families and homologous genes from species of varying divergence can be used to identify conserved noncoding regulatory elements. In this study, phylogenetic analysis of 5'-noncoding sequences was optimized using rab17, a well-characterized ABA-responsive gene from maize, and five additional rab16/17 homologs from sorghum and rice. Conserved 5'-noncoding sequences among the maize, sorghum, and rice rab16/17 homologs were identified with the aid of the software program FootPrinter and by screening for known transcription-factor-binding sites. Searches for 7 of 8 (7/8)bp sequence matches within aligned 5'-noncoding segments of the rab genes identified many of the cis-elements previously characterized by biochemical analysis in maize rab17 plus several additional putative regulatory elements. Differences in the composition of conserved noncoding sequences among rab16/17 genes were related to variation in rab gene mRNA levels in different tissues and to response to ABA treatment using qRT-PCR. Absence of a GRA-like element in the promoter of sorghum dhn2 relative to maize rab17 was correlated with an approximately 85-fold reduction of dhn2 RNA in sorghum shoots. Overall, we conclude that phylogenetic analysis of gene families among rice, sorghum, and maize will help identify regulatory sequences in the noncoding regions of genes and contribute to our understanding of grass gene regulatory networks.

Molecular markers from the transcribed/expressed region of the genome in higher plants

In recent years, molecular marker technology in higher plants has witnessed a shift from the so-called random DNA markers (RDMs), developed in the past arbitrarily from genomic DNA and cDNA, to the molecular markers representing the transcriptome and the other coding sequences. These markers have been described as gene targeted markers (GTMs). Another specific class of markers includes the so-called functional markers (FMs), which are supposed to have a cause and effect relationship with the traits of interest. In this review, we first describe the development of these markers representing the transcriptome or genes per se; we then discuss the uses of these markers in some detail and finally add a note on the future directions of research and the implications of the wider application of these markers in crop improvement programmes. Using suitable examples, we describe markers of different classes derived from cDNA clones, expressed sequence tags (ESTs), gene sequences and the unique (coding) sequences obtained through methyl filtration or genome normalization (high C(0) t fraction) from gDNA libraries. While we briefly describe RFLPs, SSRs, AFLPs and SNPs developed from the transcriptome (cDNA clones and EST databases), we have discussed in more detail some of the novel markers developed from the transcriptome and specific genes. These novel markers include expressed sequence tag polymorphisms (ESTPs), conserved orthologue set (COS) markers, amplified consensus genetic markers (ACGMs), gene specific tags (GSTs), resistance gene analogues (RGAs) and exon-retrotransposon amplification polymorphism (ERAP). Uses of these markers have been discussed in some detail under the following headings: development of transcript and functional maps, estimations of genetic diversity, marker-assisted selection (MAS), candidate-gene (CG) approach and map-based cloning, genetical genomics and identification of eQTLs, study of genome organization and taxonomic and phylogenetic studies. At the end, we also append a list of websites relevant to further studies on the transcriptome. For want of space, considerable information including voluminous data in the form of 12 tables, and a long list of references cited in these tables, has been placed on the Internet as electronic supplementary material (ESM), which the readers may find useful.

Maize Genome Sequencing by Methylation Filtration

Gene enrichment strategies offer an alternative to sequencing large and repetitive genomes such as that of maize. We report the generation and analysis of nearly 100,000 undermethylated (or methylation filtration) maize sequences. Comparison with the rice genome reveals that methylation filtration results in a more comprehensive representation of maize genes than those that result from expressed sequence tags or transposon insertion sites sequences. About 7% of the repetitive DNA is unmethylated and thus selected in our libraries, but potentially active transposons and unmethylated organelle genomes can be identified. Reverse transcription polymerase chain reaction can be used to finish the maize transcriptome.

Sequence-based alignment of sorghum chromosome 3 and rice chromosome 1 reveals extensive conservation of gene order and one major chromosomal rearrangement

The completed rice genome sequence will accelerate progress on the identification and functional classification of biologically important genes and serve as an invaluable resource for the comparative analysis of grass genomes. In this study, methods were developed for sequence-based alignment of sorghum and rice chromosomes and for refining the sorghum genetic/physical map based on the rice genome sequence. A framework of 135 BAC contigs spanning approximately 33 Mbp was anchored to sorghum chromosome 3. A limited number of sequences were collected from 118 of the BACs and subjected to BLASTX analysis to identify putative genes and BLASTN analysis to identify sequence matches to the rice genome. Extensive conservation of gene content and order between sorghum chromosome 3 and the homeologous rice chromosome 1 was observed. One large-scale rearrangement was detected involving the inversion of an approximately 59 cM block of the short arm of sorghum chromosome 3. Several small-scale changes in gene collinearity were detected, indicating that single genes and/or small clusters of genes have moved since the divergence of sorghum and rice. Additionally, the alignment of the sorghum physical map to the rice genome sequence allowed sequence-assisted assembly of an approximately 1.6 Mbp sorghum BAC contig. This streamlined approach to high-resolution genome alignment and map building will yield important information about the relationships between rice and sorghum genes and genomic segments and ultimately enhance our understanding of cereal genome structure and evolution.

Conserved noncoding sequences among cultivated cereal genomes identify candidate regulatory sequence elements and patterns of promoter evolution

Surveys for conserved noncoding sequences (CNS) among genes from monocot cereal species were conducted to assess the general properties of CNS in grass genomes and their correlation with known promoter regulatory elements. Initial comparisons of 11 orthologous maize-rice gene pairs found that previously defined regulatory motifs could be identified within short CNS but could not be distinguished reliably from random sequence matches. Among the different phylogenetic footprinting algorithms tested, the VISTA tool yielded the most informative alignments of noncoding sequence. VISTA was used to survey for CNS among all publicly available genomic sequences from maize, rice, wheat, barley, and sorghum, representing >300 gene comparisons. Comparisons of orthologous maize-rice and maize-sorghum gene pairs identified 20 bp as a minimal length criterion for a significant CNS among grass genes, with few such CNS found to be conserved across rice, maize, sorghum, and barley. The frequency and length of cereal CNS as well as nucleotide substitution rates within CNS were consistent with the known phylogenetic distances among the species compared. The implications of these findings for the evolution of cereal gene promoter sequences and the utility of using the nearly completed rice genome sequence to predict candidate regulatory elements in other cereal genes by phylogenetic footprinting are discussed.

Gene enrichment in plant genomic shotgun libraries

The Arabidopsis genome (about 130 Mbp) has been completely sequenced; whereas a draft sequence of the rice genome (about 430 Mbp) is now available and the sequencing of this genome will be completed in the near future. The much larger genomes of several important crop species, such as wheat (about 16,000 Mbp) or maize (about 2500 Mbp), may not be fully sequenced with current technology. Instead, sequencing-analysis strategies are being developed to obtain sequencing and mapping information selectively for the genic fraction (gene space) of complex plant genomes.