CleanBSequences: an efficient curator of biological sequences in R.

This work presents a new method and tool to solve a common problem of molecular biologists and geneticists who use molecular markers in their scientific research and developments: curation of sequences. Omic studies conducted by molecular biologists and geneticists usually involve the use of molecular markers. AFLP, cDNA-AFLP, and MSAP are examples of markers that render information at the genomics, transcriptomics, and epigenomics levels, respectively. These three types of molecular markers use adaptors that are the template for PCR amplification. The sequences of the adaptors have to be eliminated for the analysis of the results. Since a large number of sequences are usually obtained in these studies, this clean-up of the data could demand long time and work. To automate this work, an R package, named CleanBSequences, was created that allows the sequences to be curated massively, quickly, without errors and can be used offline. The curating is performed by aligning the forward and/or reverse primers or ends of cloning vectors with the sequences to be removed. After the alignment, new subsequences are generated without biological fragments not desired by the user, i.e., sequences needed by the techniques. In conclusion, the CleanBSequences tool facilitates the work of researchers, reducing time, effort, and working errors. Therefore, the present tool would respond to the problems related to the curation of sequences obtained from the use of some types of molecular markers. In addition to the above, being an open source, CleanBSequences is a flexible tool that has the potential to be used in future improvements to respond to new problems.

The MAP3K-Coding QUI-GON JINN (QGJ) Gene Is Essential to the Formation of Unreduced Embryo Sacs in Paspalum.

Apomixis is a clonal mode of reproduction via seeds, which results from the failure of meiosis and fertilization in the sexual female reproductive pathway. In previous transcriptomic surveys, we identified a mitogen-activated protein kinase kinase kinase (N46) displaying differential representation in florets of sexual and apomictic Paspalum notatum genotypes. Here, we retrieved and characterized the N46 full cDNA sequence from sexual and apomictic floral transcriptomes. Phylogenetic analyses showed that N46 was a member of the YODA family, which was re-named QUI-GON JINN (QGJ). Differential expression in florets of sexual and apomictic plants was confirmed by qPCR. In situ hybridization experiments revealed expression in the nucellus of aposporous plants' ovules, which was absent in sexual plants. RNAi inhibition of QGJ expression in two apomictic genotypes resulted in significantly reduced rates of aposporous embryo sac formation, with respect to the level detected in wild type aposporous plants and transformation controls. The QGJ locus segregated independently of apospory. However, a probe derived from a related long non-coding RNA sequence (PN_LNC_QGJ) revealed RFLP bands cosegregating with the Paspalum apospory-controlling region (ACR). PN_LNC_QGJ is expressed in florets of apomictic plants only. Our results indicate that the activity of QGJ in the nucellus of apomictic plants is necessary to form non-reduced embryo sacs and that a long non-coding sequence with regulatory potential is similar to sequences located within the ACR.

Acetohydroxyacid synthase activity and transcripts profiling reveal tissue-specific regulation of ahas genes in sunflower.

Acetohydroxyacid synthase (AHAS) is the target site of several herbicides and catalyses the first step in the biosynthesis of branched chain amino acid. Three genes coding for AHAS catalytic subunit (ahas1, ahas2 and ahas3) have been reported for sunflower. The aim of this work was to study the expression pattern of ahas genes family and AHAS activity in sunflower (Helianthus annuus L.). Different organs (leaves, hypocotyls, roots, flowers and embryos) were evaluated at several developmental stages. The transcriptional profile was studied through RT-qPCR. The highest expression for ahas1 was shown in leaves, where all the induced and natural gene mutations conferring herbicide resistance were found. The maximal expression of ahas2 and ahas3 occurred in immature flowers and embryos. The highest AHAS activity was found in leaves and immature embryos. Correlation analysis among ahas gene expression and AHAS activity was discussed. Our results show that differences in ahas genes expression are tissue-specific and temporally regulated. Moreover, the conservation of multiple AHAS isoforms in sunflower seems to result from different expression requirements controlled by tissue-specific regulatory mechanisms at different developmental stages.

Characterization and expression analysis of SOMATIC EMBRYOGENESIS RECEPTOR KINASE (SERK) genes in sexual and apomictic Paspalum notatum.

The SOMATIC EMBRYOGENESIS RECEPTOR-LIKE KINASE (SERK) gene plays a fundamental role in somatic embryogenesis of angiosperms, and is associated with apomixis in Poa pratensis. The objective of this work was to isolate, characterize and analyze the expression patterns of SERK genes in apomictic and sexual genotypes of Paspalum notatum. A conserved 200-bp gene fragment was amplified from genomic DNA with heterologous primers, and used to initiate a chromosomal walking strategy for cloning the complete sequence. This procedure allowed the isolation of two members of the P. notatum SERK family; PnSERK1, which is similar to PpSERK1, and PnSERK2, which is similar to ZmSERK2 and AtSERK1. Phylogenetic analyses indicated that PnSERK1 and PnSERK2 represent paralogous sequences. Southern-blot hybridization indicated the presence of at least three copies of SERK genes in the species. qRT-PCR analyses revealed that PnSERK2 was expressed at significantly higher levels than PnSERK1 in roots, leaves, reproductive tissues and embryogenic calli. Moreover, in situ hybridization experiments revealed that PnSERK2 displayed a spatially and chronologically altered expression pattern in reproductive organs of the apomictic genotype with respect to the sexual one. PnSERK2 is expressed in nucellar cells of the apomictic genotype at meiosis, but only in the megaspore mother cell in the sexual genotype. Therefore, apomixis onset in P. notatum seems to be correlated with the expression of PnSERK2 in nucellar tissue.

Differential expression of acetohydroxyacid synthase genes in sunflower plantlets and its response to imazapyr herbicide.

Acetohydroxyacid synthase (AHAS) catalyzes the first reaction in branch chain amino acids biosynthesis. This enzyme is the target of several herbicides, including all members of the imidazolinone family. Little is known about the expression of the three acetohydroxyacid synthase genes (ahas1, ahas2 and ahas3) in sunflower. The aim of this work was to evaluate ahas gene expression and AHAS activity in different tissues of sunflower plantlets. Three genotypes differing in imidazolinone resistance were evaluated, two of which carry an herbicide resistant-endowing mutation known as Ahasl1-1 allele. In vivo and in vitro AHAS activity and transcript levels were higher in leaves than in roots. The ahas3 transcript was the less abundant in both tissues. No significant difference was observed between ahas1 and ahas2 transcript levels of the susceptible genotype but a higher ahas1 transcript level was observed in leaves of genotypes carrying Ahasl1-1 allele. Similar transcript levels were found for ahas1 and ahas2 in roots of genotypes carrying Ahasl1-1 allele whereas higher ahas2 abundance was found in the susceptible genotype. Herbicide treatment triggered tissue-specific, gene and genotype-dependent changes in ahas gene expression. AHAS activity was highly inhibited in the susceptible genotype. Differential responses were observed between in vitro and in vivo AHAS inhibition assays. These findings enhance our understanding of AHAS expression in sunflower genotypes differing for herbicide resistance and its response to herbicide treatment.

Sequence characterization, in silico mapping and cytosine methylation analysis of markers linked to apospory in Paspalum notatum.

In previous studies we reported the identification of several AFLP, RAPD and RFLP molecular markers linked to apospory in Paspalum notatum. The objective of this work was to sequence these markers, obtain their flanking regions by chromosome walking and perform an in silico mapping analysis in rice and maize. The methylation status of two apospory-related sequences was also assessed using methylation-sensitive RFLP experiments. Fourteen molecular markers were analyzed and several protein-coding sequences were identified. Copy number estimates and RFLP linkage analysis showed that the sequence PnMAI3 displayed 2-4 copies per genome and linkage to apospory. Extension of this marker by chromosome walking revealed an additional protein-coding sequence mapping in silico in the apospory-syntenic regions of rice and maize. Approximately 5 kb corresponding to different markers were characterized through the global sequencing procedure. A more refined analysis based on sequence information indicated synteny with segments of chromosomes 2 and 12 of rice and chromosomes 3 and 5 of maize. Two loci associated with apomixis locus were tested in methylation-sensitive RFLP experiments using genomic DNA extracted from leaves. Although both target sequences were methylated no methylation polymorphisms associated with the mode of reproduction were detected.

Expression of lorelei-like genes in aposporous and sexual Paspalum notatum plants.

Gametophytic apomictic plants form non-reduced embryo sacs that generate clonal embryos by parthenogenesis, in the absence of both meiosis and egg-cell fertilization. Here we report the sequence and expression analysis of a lorelei-like Paspalum notatum gene, n20gap-1, which encodes a GPI-anchored protein previously associated with apomixis in this species. Phylogeny trees showed that n20gap-1 was evolutionary related to the Arabidopsis thaliana lorelei genes At4g26466 and At5g56170. The lorelei At4g26466 disruption was shown to be detrimental to sperm cell release in arabidopsis. RFLP (Restriction Fragment Length Polymorphism) analysis revealed the occurrence of several homologous sequences in the Paspalum notatum genome, exhibiting polymorphisms genetically linked to apomixis. Real-time PCR showed that lorelei-family genes present a minor activity peak at pre-meiosis and a major one at anthesis. The apomictic genotype analyzed showed a significantly increased activity at pre-meiosis, post-meiosis and anthesis with respect to a sexual genotype. In situ hybridization assays revealed expression in integuments, nucellus and the egg-cell apparatus. Several n20gap-1 alleles differing mainly at the 3' UTR sequence were identified. Allele-specific real-time PCR experiments showed that allele 28 was significantly induced in reproductive tissues of the apomictic genotype with respect to the sexual genotype at anthesis. Our results indicate that P. notatum lorelei-like genes are differentially expressed in representative sexual (Q4188) and apomictic (Q4117) genotypes, and might play a role in the final stages of the apomixis developmental cascade. However, the association of n20gap-1 expression with the trait should be confirmed in significant number of sexual and apomictic genotypes.

Expressed sequence tag analysis and development of gene associated markers in a near-isogenic plant system of Eragrostis curvula.

Eragrostis curvula (Schrad.) Nees is a forage grass native to the semiarid regions of Southern Africa, which reproduces mainly by pseudogamous diplosporous apomixis. A collection of ESTs was generated from four cDNA libraries, three of them obtained from panicles of near-isogenic lines with different ploidy levels and reproductive modes, and one obtained from 12 days-old plant leaves. A total of 12,295 high-quality ESTs were clustered and assembled, rendering 8,864 unigenes, including 1,490 contigs and 7,394 singletons, with a genome coverage of 22%. A total of 7,029 (79.11%) unigenes were functionally categorized by BLASTX analysis against sequences deposited in public databases, but only 37.80% could be classified according to Gene Ontology. Sequence comparison against the cereals genes indexes (GI) revealed 50% significant hits. A total of 254 EST-SSRs were detected from 219 singletons and 35 from contigs. Di- and tri- motifs were similarly represented with percentages of 38.95 and 40.16%, respectively. In addition, 190 SNPs and Indels were detected in 18 contigs generated from 3 to 4 libraries. The ESTs and the molecular markers obtained in this study will provide valuable resources for a wide range of applications including gene identification, genetic mapping, cultivar identification, analysis of genetic diversity, phenotype mapping and marker assisted selection.

Genome polymorphisms and gene differential expression in a 'back-and-forth' ploidy-altered series of weeping lovegrass (Eragrostis curvula).

Molecular markers were used to analyze the genomic structure of an euploid series of Eragrostis curvula, obtained after a tetraploid dihaploidization procedure followed by chromosome re-doubling with colchicine. Considerable levels of genome polymorphisms were detected between lines. Curiously, a significant number of molecular markers showed a revertant behavior following the successive changes of ploidy, suggesting that genome alterations were specific and conferred genetic structures characteristic of a given ploidy level. Genuine reversion was confirmed by sequencing. Cluster analysis demonstrated grouping of tetraploids while the diploid was more distantly related with respect to the rest of the plants. Polymorphic revertant sequences involved mostly non-coding regions, although some of them displayed sequence homology to known genes. A revertant sequence corresponding to a P-type adenosine triphosphatase was found to be differentially represented in cDNA libraries obtained from the diploid and a colchiploid, but was not found expressed in the original tetraploid. Transcriptome profiling of inflorescence followed by real-time polymerase chain reaction validation showed 0.34% polymorphic bands between apomictic tetraploid and sexual diploid plants. Several of the polymorphic sequences corresponded to known genes. Possible correlation between the results observed here and a recently reported genome-wide non-Mendelian inheritance mechanism in Arabidopsis thaliana are discussed.

A complex gene cluster for indole-diterpene biosynthesis in the grass endophyte Neotyphodium lolii.

Lolitrems are a structurally diverse group of indole-diterpene mycotoxins synthesized by Epichloë/Neotyphodium endophytes in association with Pooid grasses. Using suppression subtractive hybridization combined with chromosome walking, two clusters of genes for lolitrem biosynthesis were isolated from Neotyphodium lolii, a mutualistic endophyte of perennial ryegrass. The first cluster contains five genes, ltmP, ltmQ, ltmF, ltmC, and ltmB, four of which appear to be orthologues of functionally characterized genes from Penicillium paxilli. The second cluster contains two genes, ltmE and ltmJ, that appear to be unique to lolitrem biosynthesis. The two clusters are separated by a 16 kb AT-rich sequence that includes two imperfect direct repeats. A previously isolated ltm cluster composed of ltmG, ltmM, and ltmK, is linked to these two new clusters by 35 kb of AT-rich retrotransposon relic sequence. All 10 genes at this complex LTM locus were highly expressed in planta but expression was very low or undetectable in mycelia. ltmM and ltmC were shown to be functional orthologues of P. paxilli paxM and paxC, respectively. This work provides a genetic foundation for elucidating the metabolic grid responsible for the diversity of indole-diterpenes synthesized by N. lolii.

Transcriptome analysis of Neotyphodium and Epichloë grass endophytes.

Large-scale gene discovery has been performed for the grass fungal endophytes Neotyphodium coenophialum, Neotyphodium lolii, and Epichloë festucae. The resulting sequences have been annotated by comparison with public DNA and protein sequence databases and using intermediate gene ontology annotation tools. Endophyte sequences have also been analysed for the presence of simple sequence repeat and single nucleotide polymorphism molecular genetic markers. Sequences and annotation are maintained within a MySQL database that may be queried using a custom web interface. Two cDNA-based microarrays have been generated from this genome resource. They permit the interrogation of 3806 Neotyphodium genes (Nchip microarray), and 4195 Neotyphodium and 920 Epichloë genes (EndoChip microarray), respectively. These microarrays provide tools for high-throughput transcriptome analysis, including genome-specific gene expression studies, profiling of novel endophyte genes, and investigation of the host grass-symbiont interaction. Comparative transcriptome analysis in Neotyphodium and Epichloë was performed.