Whole Genome Analysis of Gene Expression Reveals Coordinated Activation of Signaling and Metabolic Pathways during Pollen-Pistil Interactions in Arabidopsis

Plant reproduction depends on the concerted activation of many genes to ensure correct communication between pollen and pistil. Here, we queried the whole transcriptome of Arabidopsis (Arabidopsis thaliana) in order to identify genes with specific reproductive functions. We used the Affymetrix ATH1 whole genome array to profile wild-type unpollinated pistils and unfertilized ovules. By comparing the expression profile of pistils at 0.5, 3.5, and 8.0 h after pollination and applying a number of statistical and bioinformatics criteria, we found 1,373 genes differentially regulated during pollen-pistil interactions. Robust clustering analysis grouped these genes in 16 time-course clusters representing distinct patterns of regulation. Coregulation within each cluster suggests the presence of distinct genetic pathways, which might be under the control of specific transcriptional regulators. A total of 78% of the regulated genes were expressed initially in unpollinated pistil and/or ovules, 15% were initially detected in the pollen data sets as enriched or preferentially expressed, and 7% were induced upon pollination. Among those, we found a particular enrichment for unknown transcripts predicted to encode secreted proteins or representing signaling and cell wall-related proteins, which may function by remodeling the extracellular matrix or as extracellular signaling molecules. A strict regulatory control in various metabolic pathways suggests that fine-tuning of the biochemical and physiological cellular environment is crucial for reproductive success. Our study provides a unique and detailed temporal and spatial gene expression profile of in vivo pollen-pistil interactions, providing a framework to better understand the basis of the molecular mechanisms operating during the reproductive process in higher plants.

Genetic Transformation of Passionflower and Evaluation of R(1) and R(2) Generations for Resistance to Cowpea aphid borne mosaic virus

We report on the production and evaluation of passionflower transgenic lines for resistance to Cowpea aphid borne mosaic virus (CABMV). Genetic transformation was done using Agrobacterium tumefaciens and transgene integration was confirmed by Southern blot analyses, resulting in nine transgenic lines for `IAC 275` and three for `IAC 277`. Transgenic lines were clonally propagated and evaluated for resistance to CABMV After the third inoculation, under higher inoculum pressure, only propagated plants of the transgenic line T16 remained asymptomatic, indicating a high resistance to infection with CABMV. This transgenic line was self-pollinated and the RI generation was evaluated together with the RI generation of another resistant transgenic line (T2) identified previously. Plants were inoculated with CABMV by means of viruliferous Myzus nicotianae. All 524 T2R(1) plants became infected, whereas 13 of 279 T16R(1) remained asymptomatic after four successive inoculations. A TI6R(2) generation was obtained and plants were inoculated with CABMV mechanically or by aphids. After successive inoculations, 118 of 258 plants were symptomless, suggesting that the resistance to CABMV was maintained in the plant genome as the homozygous condition was achieved. Five selected resistant TI6R(2) plants which contained the capsid protein gene are being crossed for further analyses.

Transformation of Lotus japonicus using the herbicide resistance bar gene as a selectable marker

Transgenic plants of the model legume Lotus japonicus were regenerated by hypocotyl transformation using a bar gene as a selectable marker. The bar encodes for Phosphinothricin Acetyl Transferase that detoxifies phosphinothricin (PPT), the active ingredient of herbicides such as Ignite (AgrEvo) and Basta (Hoechst). Transgenic L. japonicus plants resistant to PPT were positive upon PCR by bar gene-specific primers. In 5 out of 7 independent lines tested, PPT resistance segregated as a single dominant allele indicating a single T-DNA insertion into the plant genome. All regenerated plants were fertile and void of visible somaclonal abnormalities contrary to 14% infertility when antibiotic selectable markers were used. The lack of somaclonal variation, ease of PPT application and low cost of PPT makes this protocol an attractive alternative for the regeneration of transgenic L. japonicus. The production of PPT herbicide-resistant L. japonicus plants may have significant commercial applications in crop production.

Versatile gene-specific sequence tags for Arabidopsis functional genomics: transcript profiling and reverse genetics applications.

Microarray transcript profiling and RNA interference are two new technologies crucial for large-scale gene function studies in multicellular eukaryotes. Both rely on sequence-specific hybridization between complementary nucleic acid strands, inciting us to create a collection of gene-specific sequence tags (GSTs) representing at least 21,500 Arabidopsis genes and which are compatible with both approaches. The GSTs were carefully selected to ensure that each of them shared no significant similarity with any other region in the Arabidopsis genome. They were synthesized by PCR amplification from genomic DNA. Spotted microarrays fabricated from the GSTs show good dynamic range, specificity, and sensitivity in transcript profiling experiments. The GSTs have also been transferred to bacterial plasmid vectors via recombinational cloning protocols. These cloned GSTs constitute the ideal starting point for a variety of functional approaches, including reverse genetics. We have subcloned GSTs on a large scale into vectors designed for gene silencing in plant cells. We show that in planta expression of GST hairpin RNA results in the expected phenotypes in silenced Arabidopsis lines. These versatile GST resources provide novel and powerful tools for functional genomics.

Preferential binding of the methyl-CpG binding domain protein 2 at methylated transcriptional start site regions.

Methyl-CpG Binding Domain (MBD) proteins are thought to be key molecules in the interpretation of DNA methylation signals leading to gene silencing through recruitment of chromatin remodeling complexes. In cancer, the MBD-family member, MBD2, may be primarily involved in the repression of genes exhibiting methylated CpG at their 5' end. Here we ask whether MBD2 randomly associates methylated sequences, producing chance effects on transcription, or exhibits a more specific recognition of some methylated regions. Using chromatin and DNA immunoprecipitation, we analyzed MBD2 and RNA polymerase II deposition and DNA methylation in HeLa cells on arrays representing 25,500 promoter regions. This first whole-genome mapping revealed the preferential localization of MBD2 near transcription start sites (TSSs), within the region analyzed, 7.5 kb upstream through 2.45 kb downstream of 5' transcription start sites. Probe by probe analysis correlated MBD2 deposition and DNA methylation. Motif analysis did not reveal specific sequence motifs; however, CCG and CGC sequences seem to be overrepresented. Nonrandom association (multiple correspondence analysis, p < 0.0001) between silent genes, DNA methylation and MBD2 binding was observed. The association between MBD2 binding and transcriptional repression weakened as the distance between binding site and TSS increased, suggesting that MBD2 represses transcriptional initiation. This hypothesis may represent a functional explanation for the preferential binding of MBD2 at methyl-CpG in TSS regions.

Dissection of the complex phenotype in cuticular mutants of Arabidopsis reveals a role of SERRATE as a mediator.

Mutations in LACERATA (LCR), FIDDLEHEAD (FDH), and BODYGUARD (BDG) cause a complex developmental syndrome that is consistent with an important role for these Arabidopsis genes in cuticle biogenesis. The genesis of their pleiotropic phenotypes is, however, poorly understood. We provide evidence that neither distorted depositions of cutin, nor deficiencies in the chemical composition of cuticular lipids, account for these features, instead suggesting that the mutants alleviate the functional disorder of the cuticle by reinforcing their defenses. To better understand how plants adapt to these mutations, we performed a genome-wide gene expression analysis. We found that apparent compensatory transcriptional responses in these mutants involve the induction of wax, cutin, cell wall, and defense genes. To gain greater insight into the mechanism by which cuticular mutations trigger this response in the plants, we performed an overlap meta-analysis, which is termed MASTA (MicroArray overlap Search Tool and Analysis), of differentially expressed genes. This suggested that different cell integrity pathways are recruited in cesA cellulose synthase and cuticular mutants. Using MASTA for an in silico suppressor/enhancer screen, we identified SERRATE (SE), which encodes a protein of RNA-processing multi-protein complexes, as a likely enhancer. In confirmation of this notion, the se lcr and se bdg double mutants eradicate severe leaf deformations as well as the organ fusions that are typical of lcr and bdg and other cuticular mutants. Also, lcr does not confer resistance to Botrytis cinerea in a se mutant background. We propose that there is a role for SERRATE-mediated RNA signaling in the cuticle integrity pathway.

Sequences of the coat protein gene from brazilian isolates of Papaya ringspot virus

Papaya ringspot virus (PRSV) is the causal agent of the main papaya (Carica papaya) disease in the world. Brazil is currently the world's main papaya grower, responsible for about 40% of the worldwide production. Resistance to PRSV on transgenic plants expressing the PRSV coat protein (cp) gene was shown to be dependent on the sequence homology between the cp transgene expressed in the plant genome and the cp gene from the incoming virus, in an isolate-specific fashion. Therefore, knowledge of the degree of homology among the cp genes from distinct PRSV isolates which are present in a given area is important to guide the development of transgenic papaya for the control of PRSV in that area. The objective of the present study was to assess the degree of homology among the PRSV cp genes of several Brazilian isolates of this virus. Papaya and PRSV are present in many different ecosystems within Brazil. Twelve PRSV isolates, collected in eight different states from four different geographic regions, were used in this study. The sequences of the cp gene from these isolates were compared among themselves and to the gene used to generate transgenic papaya for Brazil. An average degree of homology of 97.3% at the nucleotide sequence was found among the Brazilian isolates. When compared to 27 isolates from outside Brazil in a homology tree, the Brazilian isolates were clustered with Australian, Hawaiian, and Central and North American isolates, with an average degree of homology of 90.7% among them.

Integration of retrotransposons-based markers in a linkage map of barley

A deeper understanding of random markers is important if they are to be employed for a range of objectives. The sequence specific amplified polymorphism (S-SAP) technique is a powerful genetic analysis tool which exploits the high copy number of retrotransposon long terminal repeats (LTRs) in the plant genome. The distribution and inheritance of S-SAP bands in the barley genome was studied using the Steptoe × Morex (S × M) double haploid (DH) population. Six S-SAP primer combinations generated 98 polymorphic bands, and map positions were assigned to all but one band. Eight putative co-dominant loci were detected, representing 16 of the mapped markers. Thus at least 81 of the mapped S-SAP loci were dominant. The markers were distributed along all of the seven chromosomes and a tendency to cluster was observed. The distribution of S-SAP markers over the barley genome concurred with the knowledge of the high copy number of retrotransposons in plants. This experiment has demonstrated the potential for the S-SAP technique to be applied in a range of analyses such as genetic fingerprinting, marker assisted breeding, biodiversity assessment and phylogenetic analyses.

Identification and validation of candidate genes associated with domesticated and improved traits in soybean

Soybean, an important source of vegetable oils and proteins for humans, has undergone significant phenotypic changes during domestication and improvement. However, there is limited knowledge about genes related to these domesticated and improved traits, such as flowering time, seed development, alkaline-salt tolerance, and seed oil content (SOC). In this study, more than 106,000 single nucleotide polymorphisms (SNPs) were identified by restriction site associated DNA sequencing of 14 wild, 153 landrace, and 119 bred soybean accessions, and 198 candidate domestication regions (CDRs) were identified via multiple genetic diversity analyses. Of the 1489 candidate domestication genes (CDGs) within these CDRs, a total of 330 CDGs were related to the above four traits in the domestication, gene ontology (GO) enrichment, gene expression, and pathway analyses. Eighteen, 60, 66, and 10 of the 330 CDGs were significantly associated with the above four traits, respectively. Of 134 traitassociated CDGs, 29 overlapped with previous CDGs, 11 were consistent with candidate genes in previous trait association studies, and 66 were covered by the domesticated and improved quantitative trait loci or their adjacent regions, having six common CDGs, such as one functionally characterized gene Glyma15 g17480 (GmZTL3). Of the 68 seed size (SS) and SOC CDGs, 37 were further confirmed by gene expression analysis. In addition, eight genes were found to be related to artificial selection during modern breeding. Therefore, this study provides an integrated method for efficiently identifying CDGs and valuable information for domestication and genetic research.

Large-scale transcriptome analyses reveal new genetic marker candidates of head, neck, and thyroid cancer

A detailed genome mapping analysis of 213,636 expressed sequence tags (EST) derived from nontumor and tumor tissues of the oral cavity, larynx, pharynx, and thyroid was done. Transcripts matching known human genes were identified; potential new splice variants were flagged and subjected to manual curation, pointing to 788 putatively new alternative splicing isoforms, the majority (75%) being insertion events. A subset of 34 new splicing isoforms (5% of 788 events) was selected and 23 (68%) were confirmed by reverse transcription-PCR and DNA sequencing. Putative new genes were revealed, including six transcripts mapped to well-studied chromosomes such as 22, as well as transcripts that mapped to 253 intergenic regions. In addition, 2,251 noncoding intronic RNAs, eventually involved in transcriptional regulation, were found. A set of 250 candidate markers for loss of heterozygosis or gene amplification was selected by identifying transcripts that mapped to genomic regions previously known to be frequently amplified or deleted in head, neck, and thyroid tumors. Three of these markers were evaluated by quantitative reverse transcription-PCR in an independent set of individual samples. Along with detailed clinical data about tumor origin, the information reported here is now publicly available on a dedicated Web site as a resource for further biological investigation. This first in silico reconstruction of the head, neck, and thyroid transcriptomes points to a wealth of new candidate markers that can be used for future studies on the molecular basis of these tumors. Similar analysis is warranted for a number of other tumors for which large EST data sets are available.

Microsatellite loci isolation from river buffalo using enriched partial genomic libraries

The extensive use of buffalo in agriculture, especially in developing countries, begs for genetic resources to evaluate and improve traits important to local and regional economies. Brazil presents the largest water buffalo populations in the New World, with 1 1 million heads including swamp and river types. To design rational breeding strategies for optimum utilization and conservation of available genetic variability in the Brazilian buffalo's population, it is essential to understand their genetic architecture and relationship among various breeds. This depends, in part, on the knowledge of their genetic structure based on molecular markers like microsatellites. In the present study, we developed six enriched partial genomic libraries for river buffalo using selective hybridization methods. Genomic DNA was hybridized with six different arrays of repeat motif, 5' biotinylated - (CA)(15), (CT)(15), (AGG)(8), (GAAA)(8), (GATA)(8), (AAAAC)(8) - and bound to streptavidin coated beads. The cloning process generated a total of 1920 recombinant clones. Up to date, 487 were directly sequenced for the presence of repeats, from which 13 have been positive for presence of repeats as follows: 9 for di-nucleotide repeats, 3 for tri-nucleotide repeats and 1 for tetra-nucleotide repeat. PCR primer pairs for the isolated microsatellites are under construction to determine optimum annealing temperature. These microsatellites will be useful for studies involving phylogenetic relationships, genome mapping and genetic diversity analysis within buffalo populations worldwide.

Contribuição ao mapeamento do cromossomo 9 do búfalo de rio (Bubalus bubalis) utilizando painel de células somáticas híbridas irradiadas

Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (CAPES)

Construção de uma biblioteca BAC e avaliação de marcadores para caracterização de regiões alvo do genoma do búfalo

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

Structure, organization, and expression of the alpha prolamin multigenic family bring new insights into the evolutionary relationships among grasses

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)