CFM1, a member of the CRM-domain protein family, functions in chloroplast group II intron splicing in Setaria viridis.

The chloroplast RNA splicing and ribosome maturation (CRM) domain is a RNA-binding domain found in a plant-specific protein family whose characterized members play essential roles in splicing group I and group II introns in mitochondria and chloroplasts. Together, these proteins are required for splicing of the majority of the approximately 20 chloroplast introns in land plants. Here, we provide evidence from Setaria viridis and maize that an uncharacterized member of this family, CRM Family Member1 (CFM1), promotes the splicing of most of the introns that had not previously been shown to require a CRM domain protein. A Setaria mutant expressing mutated CFM1 was strongly disrupted in the splicing of three chloroplast tRNAs: trnI, trnV and trnA. Analyses by RNA gel blot and polysome association suggest that the tRNA deficiencies lead to compromised chloroplast protein synthesis and the observed whole-plant chlorotic phenotypes. Co-immunoprecipitation data demonstrate that the maize CFM1 ortholog is bound to introns whose splicing is disrupted in the cfm1 mutant. With these results, CRM domain proteins have been shown to promote the splicing of all but two of the introns found in angiosperm chloroplast genomes.

Setaria viridis chlorotic and seedling-lethal mutants define critical functions for chloroplast gene expression.

Deep insights into chloroplast biogenesis have been obtained by mutant analysis; however, in C4 plants a relevant mutant collection has only been developed and exploited for maize. Here, we report the initial characterization of an ethyl methyl sulfonate-induced mutant population for the C4 model Setaria viridis. Approximately 1000 M2 families were screened for the segregation of pale-green seedlings in the M3 generation, and a subset of these was identified to be deficient in post-transcriptional steps of chloroplast gene expression. Causative mutations were identified for three lines using deep sequencing-based bulked segregant analysis, and in one case confirmed by transgenic complementation. Using chloroplast RNA-sequencing and other molecular assays, we describe phenotypes of mutants deficient in PSRP7, a plastid-specific ribosomal protein, OTP86, an RNA editing factor, and cpPNP, the chloroplast isozyme of polynucleotide phosphorylase. The psrp mutant is globally defective in chloroplast translation, and has varying deficiencies in the accumulation of chloroplast-encoded proteins. The otp86 mutant, like its Arabidopsis counterpart, is specifically defective in editing of the rps14 mRNA; however, the conditional pale-green mutant phenotype contrasts with the normal growth of the Arabidopsis mutant. The pnp mutant exhibited multiple defects in 3' end maturation as well as other qualitative changes in the chloroplast RNA population. Overall, our collection opens the door to global analysis of photosynthesis and early seedling development in an emerging C4 model.

Genome of 'Charleston Gray', the principal American watermelon cultivar, and genetic characterization of 1,365 accessions in the U.S. National Plant Germplasm System watermelon collection.

Years of selection for desirable fruit quality traits in dessert watermelon (Citrullus lanatus) has resulted in a narrow genetic base in modern cultivars. Development of novel genomic and genetic resources offers great potential to expand genetic diversity and improve important traits in watermelon. Here, we report a high-quality genome sequence of watermelon cultivar 'Charleston Gray', a principal American dessert watermelon, to complement the existing reference genome from '97103', an East Asian cultivar. Comparative analyses between genomes of 'Charleston Gray' and '97103' revealed genomic variants that may underlie phenotypic differences between the two cultivars. We then genotyped 1365 watermelon plant introduction (PI) lines maintained at the U.S. National Plant Germplasm System using genotyping-by-sequencing (GBS). These PI lines were collected throughout the world and belong to three Citrullus species, C. lanatus, C. mucosospermus and C. amarus. Approximately 25 000 high-quality single nucleotide polymorphisms (SNPs) were derived from the GBS data using the 'Charleston Gray' genome as the reference. Population genomic analyses using these SNPs discovered a close relationship between C. lanatus and C. mucosospermus and identified four major groups in these two species correlated to their geographic locations. Citrullus amarus was found to have a distinct genetic makeup compared to C. lanatus and C. mucosospermus. The SNPs also enabled identification of genomic regions associated with important fruit quality and disease resistance traits through genome-wide association studies. The high-quality 'Charleston Gray' genome and the genotyping data of this large collection of watermelon accessions provide valuable resources for facilitating watermelon research, breeding and improvement.

Deciphering genetic factors that determine melon fruit-quality traits using RNA-Seq-based high-resolution QTL and eQTL mapping.

Combined quantitative trait loci (QTL) and expression-QTL (eQTL) mapping analysis was performed to identify genetic factors affecting melon (Cucumis melo) fruit quality, by linking genotypic, metabolic and transcriptomic data from a melon recombinant inbred line (RIL) population. RNA sequencing (RNA-Seq) of fruit from 96 RILs yielded a highly saturated collection of > 58 000 single-nucleotide polymorphisms, identifying 6636 recombination events that separated the genome into 3663 genomic bins. Bin-based QTL analysis of 79 RILs and 129 fruit-quality traits affecting taste, aroma and color resulted in the mapping of 241 QTL. Thiol acyltransferase (CmThAT1) gene was identified within the QTL interval of its product, S-methyl-thioacetate, a key component of melon fruit aroma. Metabolic activity of CmThAT1-encoded protein was validated in bacteria and in vitro. QTL analysis of flesh color intensity identified a candidate white-flesh gene (CmPPR1), one of two major loci determining fruit flesh color in melon. CmPPR1 encodes a member of the pentatricopeptide protein family, involved in processing of RNA in plastids, where carotenoid and chlorophyll pigments accumulate. Network analysis of > 12 000 eQTL mapped for > 8000 differentially expressed fruit genes supported the role of CmPPR1 in determining the expression level of plastid targeted genes. We highlight the potential of RNA-Seq-based QTL analysis of small to moderate size, advanced RIL populations for precise marker-assisted breeding and gene discovery. We provide the following resources: a RIL population genotyped with a unique set of SNP markers, confined genomic segments that harbor QTL governing 129 traits and a saturated set of melon eQTLs.

Population differentiation in allele frequencies of obesity-associated SNPs.

BACKGROUND: Obesity is emerging as a global health problem, with more than one-third of the world's adult population being overweight or obese. In this study, we investigated worldwide population differentiation in allele frequencies of obesity-associated SNPs (single nucleotide polymorphisms). RESULTS: We collected a total of 225 obesity-associated SNPs from a public database. Their population-level allele frequencies were derived based on the genotype data from 1000 Genomes Project (phase 3). We used hypergeometric model to assess whether the effect allele at a given SNP is significantly enriched or depleted in each of the 26 populations surveyed in the 1000 Genomes Project with respect to the overall pooled population. Our results indicate that 195 out of 225 SNPs (86.7%) possess effect alleles significantly enriched or depleted in at least one of the 26 populations. Populations within the same continental group exhibit similar allele enrichment/depletion patterns whereas inter-continental populations show distinct patterns. Among the 225 SNPs, 15 SNPs cluster in the first intron region of the FTO gene, which is a major gene associated with body-mass index (BMI) and fat mass. African populations exhibit much smaller blocks of LD (linkage disequilibrium) among these15 SNPs while European and Asian populations have larger blocks. To estimate the cumulative effect of all variants associated with obesity, we developed the personal composite genetic risk score for obesity. Our results indicate that the East Asian populations have the lowest averages of the composite risk scores, whereas three European populations have the highest averages. In addition, the population-level average of composite genetic risk scores is significantly correlated (R2 = 0.35, P = 0.0060) with obesity prevalence. CONCLUSIONS: We have detected substantial population differentiation in allele frequencies of obesity-associated SNPs. The results will help elucidate the genetic basis which may contribute to population disparities in obesity prevalence.

Genome re-sequencing reveals the history of apple and supports a two-stage model for fruit enlargement.

Human selection has reshaped crop genomes. Here we report an apple genome variation map generated through genome sequencing of 117 diverse accessions. A comprehensive model of apple speciation and domestication along the Silk Road is proposed based on evidence from diverse genomic analyses. Cultivated apples likely originate from Malus sieversii in Kazakhstan, followed by intensive introgressions from M. sylvestris. M. sieversii in Xinjiang of China turns out to be an "ancient" isolated ecotype not directly contributing to apple domestication. We have identified selective sweeps underlying quantitative trait loci/genes of important fruit quality traits including fruit texture and flavor, and provide evidences supporting a model of apple fruit size evolution comprising two major events with one occurring prior to domestication and the other during domestication. This study outlines the genetic basis of apple domestication and evolution, and provides valuable information for facilitating marker-assisted breeding and apple improvement.Apple is one of the most important fruit crops. Here, the authors perform deep genome resequencing of 117 diverse accessions and reveal comprehensive models of apple origin, speciation, domestication, and fruit size evolution as well as candidate genes associated with important agronomic traits.

Host-secreted antimicrobial peptide enforces symbiotic selectivity in Medicago truncatula.

Legumes engage in root nodule symbioses with nitrogen-fixing soil bacteria known as rhizobia. In nodule cells, bacteria are enclosed in membrane-bound vesicles called symbiosomes and differentiate into bacteroids that are capable of converting atmospheric nitrogen into ammonia. Bacteroid differentiation and prolonged intracellular survival are essential for development of functional nodules. However, in the Medicago truncatula-Sinorhizobium meliloti symbiosis, incompatibility between symbiotic partners frequently occurs, leading to the formation of infected nodules defective in nitrogen fixation (Fix-). Here, we report the identification and cloning of the M. truncatula NFS2 gene that regulates this type of specificity pertaining to S. meliloti strain Rm41. We demonstrate that NFS2 encodes a nodule-specific cysteine-rich (NCR) peptide that acts to promote bacterial lysis after differentiation. The negative role of NFS2 in symbiosis is contingent on host genetic background and can be counteracted by other genes encoded by the host. This work extends the paradigm of NCR function to include the negative regulation of symbiotic persistence in host-strain interactions. Our data suggest that NCR peptides are host determinants of symbiotic specificity in M. truncatula and possibly in closely related legumes that form indeterminate nodules in which bacterial symbionts undergo terminal differentiation.

Exploring key cellular processes and candidate genes regulating the primary thickening growth of Moso underground shoots.

The primary thickening growth of Moso (Phyllostachys edulis) underground shoots largely determines the culm circumference. However, its developmental mechanisms remain largely unknown. Using an integrated anatomy, mathematics and genomics approach, we systematically studied cellular and molecular mechanisms underlying the growth of Moso underground shoots. We discovered that the growth displayed a spiral pattern and pith played an important role in promoting the primary thickening process of Moso underground shoots and driving the evolution of culms with different sizes among different bamboo species. Different with model plants, the shoot apical meristem (SAM) of Moso is composed of six layers of cells. Comparative transcriptome analysis identified a large number of genes related to the vascular tissue formation that were significantly upregulated in a thick wall variant with narrow pith cavity, mildly spiral growth, and flat and enlarged SAM, including those related to plant hormones and those involved in cell wall development. These results provide a systematic perspective on the primary thickening growth of Moso underground shoots, and support a plausible mechanism resulting in the narrow pith cavity, weak spiral growth but increased vascular bundle of the thick wall Moso.

Genome-wide analysis of Dongxiang wild rice (Oryza rufipogon Griff.) to investigate lost/acquired genes during rice domestication.

BACKGROUND: It is widely accepted that cultivated rice (Oryza sativa L.) was domesticated from common wild rice (Oryza rufipogon Griff.). Compared to other studies which concentrate on rice origin, this study is to genetically elucidate the substantially phenotypic and physiological changes from wild rice to cultivated rice at the whole genome level. RESULTS: Instead of comparing two assembled genomes, this study directly compared the Dongxiang wild rice (DXWR) Illumina sequencing reads with the Nipponbare (O. sativa) complete genome without assembly of the DXWR genome. Based on the results from the comparative genomics analysis, structural variations (SVs) between DXWR and Nipponbare were determined to locate deleted genes which could have been acquired by Nipponbare during rice domestication. To overcome the limit of the SV detection, the DXWR transcriptome was also sequenced and compared with the Nipponbare transcriptome to discover the genes which could have been lost in DXWR during domestication. Both 1591 Nipponbare-acquired genes and 206 DXWR-lost transcripts were further analyzed using annotations from multiple sources. The NGS data are available in the NCBI SRA database with ID SRP070627. CONCLUSIONS: These results help better understanding the domestication from wild rice to cultivated rice at the whole genome level and provide a genomic data resource for rice genetic research or breeding. One finding confirmed transposable elements contribute greatly to the genome evolution from wild rice to cultivated rice. Another finding suggested the photophosphorylation and oxidative phosphorylation system in cultivated rice could have adapted to environmental changes simultaneously during domestication.

HAWAIIAN SKIRT regulates the quiescent center-independent meristem activity in Arabidopsis roots.

Root apical meristem (RAM) drives post-embryonic root growth by constantly supplying cells through mitosis. It is composed of stem cells and their derivatives, the transit-amplifying (TA) cells. Stem cell organization and its maintenance in the RAM are well characterized, however, their relationships with TA cells remain unclear. SHORTROOT (SHR) is critical for root development. It patterns cell types and promotes the post-embryonic root growth. Defective root growth in the shr has been ascribed to the lack of quiescent center (QC), which maintains the surrounding stem cells. However, our recent investigation indicated that SHR maintains TA cells independently of QC by modulating PHABULOSA (PHB) through miRNA165/6. PHB controls TA cell activity by modulating cytokinin levels and type B Arabidopsis Response Regulator activity, in a dosage-dependent manner. To further understand TA cell regulation, we conducted a shr suppressor screen. With an extensive mutagenesis screen followed by genome sequencing of a pooled F2 population, we discovered two suppressor alleles with mutations in HAWAIIAN SKIRT (HWS). HWS, encoding an F-box protein with kelch domain, is expressed, partly depending on SHR, in the root cap and in the pericycle of the differentiation zone. Interestingly, root growth in the shr hws was more active than the wild-type roots for the first 7 days after germination, without recovering QC. Contrary to shr phb, shr hws did not show a recovery of cytokinin signaling. These indicate that HWS affects QC-independent TA cell activities through a pathway distinctive from PHB.

Messenger RNA exchange between scions and rootstocks in grafted grapevines.

BACKGROUND: Grafting has been widely practiced for centuries in the propagation and production of many vegetable and fruit species. However, the underlying molecular and genetic mechanisms for how the graft partners interact with each other to produce a successful graft remain largely unknown. We hypothesized that genome-wide mRNA exchanges, which were recently documented in grafted model plant species, are a general phenomenon widely present in grafted plants, including those in vegetable and fruit species, and have specific genotype- and environment-dependent characteristics modulating plant performance. METHODS: Using diagnostic SNPs derived from high throughput genome sequencing, we identified and characterized the patterns of genome-wide mRNA exchanges across graft junctions in grafted grapevines grown in the in vitro and field conditions. RESULTS: We identified more than 3000 genes transporting mRNAs across graft junctions. These genes were involved in diverse biological processes and those involved in basic cellular, biosynthetic, catabolic, and metabolic activities, as well as responses to stress and signal transduction, were highly enriched. Field-grown mature grafts had much fewer genes transmitting mRNAs than the in vitro young grafts (987 vs. 2679). These mobile mRNAs could move directionally or bi-directionally between scions and rootstocks. The mRNA transmission rates of these genes were generally low, with 65% or more having transmission rates lower than 0.01. Furthermore, genotypes, graft combinations and growth environments had impact on the directions of mRNA movement as well as the numbers and species of mRNAs being exchanged. Moreover, we found evidence for the presences of both passive and selective mechanisms underlying long distance mRNA trafficking in grafted grapevines. CONCLUSIONS: We extended the studies of mRNA exchanges in model species to grapevines and demonstrated that genomic-scale mRNA exchange across graft junctions occurred in grapevines in a passive or genotype and environment-dependent manner.

Genome-Wide Mapping of Structural Variations Reveals a Copy Number Variant That Determines Reproductive Morphology in Cucumber.

Structural variations (SVs) represent a major source of genetic diversity. However, the functional impact and formation mechanisms of SVs in plant genomes remain largely unexplored. Here, we report a nucleotide-resolution SV map of cucumber (Cucumis sativas) that comprises 26,788 SVs based on deep resequencing of 115 diverse accessions. The largest proportion of cucumber SVs was formed through nonhomologous end-joining rearrangements, and the occurrence of SVs is closely associated with regions of high nucleotide diversity. These SVs affect the coding regions of 1676 genes, some of which are associated with cucumber domestication. Based on the map, we discovered a copy number variation (CNV) involving four genes that defines the Female (F) locus and gives rise to gynoecious cucumber plants, which bear only female flowers and set fruit at almost every node. The CNV arose from a recent 30.2-kb duplication at a meiotically unstable region, likely via microhomology-mediated break-induced replication. The SV set provides a snapshot of structural variations in plants and will serve as an important resource for exploring genes underlying key traits and for facilitating practical breeding in cucumber.

The genome of the stress-tolerant wild tomato species Solanum pennellii.

Solanum pennellii is a wild tomato species endemic to Andean regions in South America, where it has evolved to thrive in arid habitats. Because of its extreme stress tolerance and unusual morphology, it is an important donor of germplasm for the cultivated tomato Solanum lycopersicum. Introgression lines (ILs) in which large genomic regions of S. lycopersicum are replaced with the corresponding segments from S. pennellii can show remarkably superior agronomic performance. Here we describe a high-quality genome assembly of the parents of the IL population. By anchoring the S. pennellii genome to the genetic map, we define candidate genes for stress tolerance and provide evidence that transposable elements had a role in the evolution of these traits. Our work paves a path toward further tomato improvement and for deciphering the mechanisms underlying the myriad other agronomic traits that can be improved with S. pennellii germplasm.

A genomic variation map provides insights into the genetic basis of cucumber domestication and diversity.

Most fruits in our daily diet are the products of domestication and breeding. Here we report a map of genome variation for a major fruit that encompasses ~3.6 million variants, generated by deep resequencing of 115 cucumber lines sampled from 3,342 accessions worldwide. Comparative analysis suggests that fruit crops underwent narrower bottlenecks during domestication than grain crops. We identified 112 putative domestication sweeps; 1 of these regions contains a gene involved in the loss of bitterness in fruits, an essential domestication trait of cucumber. We also investigated the genomic basis of divergence among the cultivated populations and discovered a natural genetic variant in a ß-carotene hydroxylase gene that could be used to breed cucumbers with enhanced nutritional value. The genomic history of cucumber evolution uncovered here provides the basis for future genomics-enabled breeding.

Genomic organization, phylogenetic comparison and differential expression of the SBP-box family genes in grape.

BACKGROUND: The SBP-box gene family is specific to plants and encodes a class of zinc finger-containing transcription factors with a broad range of functions. Although SBP-box genes have been identified in numerous plants including green algae, moss, silver birch, snapdragon, Arabidopsis, rice and maize, there is little information concerning SBP-box genes, or the corresponding miR156/157, function in grapevine. METHODOLOGY/PRINCIPAL FINDINGS: Eighteen SBP-box gene family members were identified in Vitis vinifera, twelve of which bore sequences that were complementary to miRNA156/157. Phylogenetic reconstruction demonstrated that plant SBP-domain proteins could be classified into seven subgroups, with the V. vinifera SBP-domain proteins being more closely related to SBP-domain proteins from dicotyledonous angiosperms than those from monocotyledonous angiosperms. In addition, synteny analysis between grape and Arabidopsis demonstrated that homologs of several grape SBP genes were found in corresponding syntenic blocks of Arabidopsis. Expression analysis of the grape SBP-box genes in various organs and at different stages of fruit development in V. quinquangularis 'Shang-24' revealed distinct spatiotemporal patterns. While the majority of the grape SBP-box genes lacking a miR156/157 target site were expressed ubiquitously and constitutively, most genes bearing a miR156/157 target site exhibited distinct expression patterns, possibly due to the inhibitory role of the microRNA. Furthermore, microarray data mining and quantitative real-time RT-PCR analysis identified several grape SBP-box genes that are potentially involved in the defense against biotic and abiotic stresses. CONCLUSION: The results presented here provide a further understanding of SBP-box gene function in plants, and yields additional insights into the mechanism of stress management in grape, which may have important implications for the future success of this crop.

The draft genome of watermelon (Citrullus lanatus) and resequencing of 20 diverse accessions.

Watermelon, Citrullus lanatus, is an important cucurbit crop grown throughout the world. Here we report a high-quality draft genome sequence of the east Asia watermelon cultivar 97103 (2n = 2× = 22) containing 23,440 predicted protein-coding genes. Comparative genomics analysis provided an evolutionary scenario for the origin of the 11 watermelon chromosomes derived from a 7-chromosome paleohexaploid eudicot ancestor. Resequencing of 20 watermelon accessions representing three different C. lanatus subspecies produced numerous haplotypes and identified the extent of genetic diversity and population structure of watermelon germplasm. Genomic regions that were preferentially selected during domestication were identified. Many disease-resistance genes were also found to be lost during domestication. In addition, integrative genomic and transcriptomic analyses yielded important insights into aspects of phloem-based vascular signaling in common between watermelon and cucumber and identified genes crucial to valuable fruit-quality traits, including sugar accumulation and citrulline metabolism.

Genome-wide identification and analysis of grape aldehyde dehydrogenase (ALDH) gene superfamily.

BACKGROUND: The completion of the grape genome sequencing project has paved the way for novel gene discovery and functional analysis. Aldehyde dehydrogenases (ALDHs) comprise a gene superfamily encoding NAD(P)(+)-dependent enzymes that catalyze the irreversible oxidation of a wide range of endogenous and exogenous aromatic and aliphatic aldehydes. Although ALDHs have been systematically investigated in several plant species including Arabidopsis and rice, our knowledge concerning the ALDH genes, their evolutionary relationship and expression patterns in grape has been limited. METHODOLOGY/PRINCIPAL FINDINGS: A total of 23 ALDH genes were identified in the grape genome and grouped into ten families according to the unified nomenclature system developed by the ALDH Gene Nomenclature Committee (AGNC). Members within the same grape ALDH families possess nearly identical exon-intron structures. Evolutionary analysis indicates that both segmental and tandem duplication events have contributed significantly to the expansion of grape ALDH genes. Phylogenetic analysis of ALDH protein sequences from seven plant species indicates that grape ALDHs are more closely related to those of Arabidopsis. In addition, synteny analysis between grape and Arabidopsis shows that homologs of a number of grape ALDHs are found in the corresponding syntenic blocks of Arabidopsis, suggesting that these genes arose before the speciation of the grape and Arabidopsis. Microarray gene expression analysis revealed large number of grape ALDH genes responsive to drought or salt stress. Furthermore, we found a number of ALDH genes showed significantly changed expressions in responses to infection with different pathogens and during grape berry development, suggesting novel roles of ALDH genes in plant-pathogen interactions and berry development. CONCLUSION: The genome-wide identification, evolutionary and expression analysis of grape ALDH genes should facilitate research in this gene family and provide new insights regarding their evolution history and functional roles in plant stress tolerance.

Arabidopsis gene co-expression network and its functional modules.

BACKGROUND: Biological networks characterize the interactions of biomolecules at a systems-level. One important property of biological networks is the modular structure, in which nodes are densely connected with each other, but between which there are only sparse connections. In this report, we attempted to find the relationship between the network topology and formation of modular structure by comparing gene co-expression networks with random networks. The organization of gene functional modules was also investigated. RESULTS: We constructed a genome-wide Arabidopsis gene co-expression network (AGCN) by using 1094 microarrays. We then analyzed the topological properties of AGCN and partitioned the network into modules by using an efficient graph clustering algorithm. In the AGCN, 382 hub genes formed a clique, and they were densely connected only to a small subset of the network. At the module level, the network clustering results provide a systems-level understanding of the gene modules that coordinate multiple biological processes to carry out specific biological functions. For instance, the photosynthesis module in AGCN involves a very large number (> 1000) of genes which participate in various biological processes including photosynthesis, electron transport, pigment metabolism, chloroplast organization and biogenesis, cofactor metabolism, protein biosynthesis, and vitamin metabolism. The cell cycle module orchestrated the coordinated expression of hundreds of genes involved in cell cycle, DNA metabolism, and cytoskeleton organization and biogenesis. We also compared the AGCN constructed in this study with a graphical Gaussian model (GGM) based Arabidopsis gene network. The photosynthesis, protein biosynthesis, and cell cycle modules identified from the GGM network had much smaller module sizes compared with the modules found in the AGCN, respectively. CONCLUSION: This study reveals new insight into the topological properties of biological networks. The preferential hub-hub connections might be necessary for the formation of modular structure in gene co-expression networks. The study also reveals new insight into the organization of gene functional modules.

Combining comparative genomics with de novo motif discovery to identify human transcription factor DNA-binding motifs.

BACKGROUND: As more and more genomes are sequenced, comparative genomics approaches provide a methodology for identifying conserved regulatory elements that may be involved in gene regulations. RESULTS: We developed a novel method to combine comparative genomics with de novo motif discovery to identify human transcription factor binding motifs that are overrepresented and conserved in the upstream regions of a set of co-regulated genes. The method is validated by analyzing a well-characterized muscle specific gene set, and the results showed that our approach performed better than the existing programs in terms of sensitivity and prediction rate. CONCLUSION: The newly developed method can be used to extract regulatory signals in co-regulated genes, which can be derived from the microarray clustering analysis.

Combining microarray and genomic data to predict DNA binding motifs.

The ability to detect regulatory elements within genome sequences is important in understanding how gene expression is controlled in biological systems. In this work, microarray data analysis is combined with genome sequence analysis to predict DNA sequences in the photosynthetic bacterium Rhodobacter sphaeroides that bind the regulators PrrA, PpsR and FnrL. These predictions were made by using hierarchical clustering to detect genes that share similar expression patterns. The DNA sequences upstream of these genes were then searched for possible transcription factor recognition motifs that may be involved in their co-regulation. The approach used promises to be widely applicable for the prediction of cis-acting DNA binding elements. Using this method the authors were independently able to detect and extend the previously described consensus sequences that have been suggested to bind FnrL and PpsR. In addition, sequences that may be recognized by the global regulator PrrA were predicted. The results support the earlier suggestions that the DNA binding sequence of PrrA may have a variable-sized gap between its conserved block elements. Using the predicted DNA binding sequences, a whole-genome-scale analysis was performed to determine the relative importance of the interplay between the three regulators PpsR, FnrL and PrrA. Results of this analysis showed that, compared to the regulation by PpsR and FnrL, a much larger number of genes are candidates to be regulated by PrrA. The study demonstrates by example that integration of multiple data types can be a powerful approach for inferring transcriptional regulatory patterns in microbial systems, and it allowed the detection of photosynthesis-related regulatory patterns in R. sphaeroides.