CCLS96.1, a member of a multicopy gene family, may encode a non-coding RNA preferentially transcribed in reproductive organs of Silene latifolia

NBLAST: a graphical user interface-based two-way BLAST software with a dot plot viewer

BLAST, a basic bioinformatics tool for searching local sequence similarity, has been one of the most widely used bioinformatics programs since its introduction in 1990. Users generally use the web-based NCBI-BLAST program for BLAST analysis. However, users with large sequence data are often faced with a problem of upload size limitation while using the web-based BLAST program. This proves inconvenient as scientists often want to run BLAST on their own data, such as transcriptome or whole genome sequences. To overcome this issue, we developed NBLAST, a graphical user interface-based BLAST program that employs a two-way system, allowing the use of input sequences either as "query" or "target" in the BLAST analysis. NBLAST is also equipped with a dot plot viewer, thus allowing researchers to create custom database for BLAST and run a dot plot similarity analysis within a single program. It is available to access to the NBLAST with http://nbitglobal.com/nblast.

Identification and characterization of mRNA-like noncoding RNAs in Salvia miltiorrhiza

Identification and characterization of 5,446 mlncRNAs from Salvia miltiorrhiza showed that the majority of identified mlncRNAs were stress responsive, providing a framework for elucidating mlncRNA functions in S. miltiorrhiza. mRNA-like noncoding RNAs (mlncRNAs) are transcribed by RNA polymerase II and are polyadenylated, capped and spliced. They play important roles in plant development and defense responses. However, there is no information available for mlncRNAs in Salvia miltiorrhiza Bunge, the first Chinese medicinal material entering the international market. To perform a transcriptome-wide identification of S. miltiorrhiza mlncRNAs, we assembled over 8 million RNA-seq reads from GenBank database and 5,624 ESTs from PlantGDB into 44422 unigenes. Using a computational identification pipeline, we identified 5446 S. miltiorrhiza mlncRNA candidates from the assembled unigenes. Of the 5446 mlncRNAs, 2 are primary transcripts of conserved miRNAs, and 2030 can be grouped into 470 families with at least two members in a family. Quantitative real-time PCR analysis of mlncRNAs with at least 900 nt showed that the majority were differentially expressed in roots, stems, leaves and flowers and responsive to methyl jasmonate (MeJA) treatment in S. miltiorrhiza. Analysis of published RNA-seq data showed that a total of 3,044 mlncRNAs were expressed in hairy roots of S. miltiorrhiza and the expression of 1,904 of the 3,044 mlncRNAs was altered by yeast extract and Ag(+) treatment. The results indicate that the majority of mlncRNAs are involved in plant response to stress. This study provides a framework for understanding the roles of mlncRNAs in S. miltiorrhiza.

Genome-wide identification and characterization of long intergenic non-coding RNAs in Ganoderma lucidum

Ganoderma lucidum is a white-rot fungus best-known for its medicinal activities. We have previously sequenced its genome and annotated the protein coding genes. However, long non-coding RNAs in G. lucidum genome have not been analyzed. In this study, we have identified and characterized long intergenic non-coding RNAs (lincRNA) in G. lucidum systematically. We developed a computational pipeline, which was used to analyze RNA-Seq data derived from G. lucidum samples collected from three developmental stages. A total of 402 lincRNA candidates were identified, with an average length of 609 bp. Analysis of their adjacent protein-coding genes (apcGenes) revealed that 46 apcGenes belong to the pathways of triterpenoid biosynthesis and lignin degradation, or families of cytochrome P450, mating type B genes, and carbohydrate-active enzymes. To determine if lincRNAs and these apcGenes have any interactions, the corresponding pairs of lincRNAs and apcGenes were analyzed in detail. We developed a modified 3' RACE method to analyze the transcriptional direction of a transcript. Among the 46 lincRNAs, 37 were found unidirectionally transcribed, and 9 were found bidirectionally transcribed. The expression profiles of 16 of these 37 lincRNAs were found to be highly correlated with those of the apcGenes across the three developmental stages. Among them, 11 are positively correlated (r>0.8) and 5 are negatively correlated (r<-0.8). The co-localization and co-expression of lincRNAs and those apcGenes playing important functions is consistent with the notion that lincRNAs might be important regulators for cellular processes. In summary, this represents the very first study to identify and characterize lincRNAs in the genomes of basidiomycetes. The results obtained here have laid the foundation for study of potential lincRNA-mediated expression regulation of genes in G. lucidum.

Comprehensive transcriptome analysis reveals novel genes involved in cardiac glycoside biosynthesis and mlncRNAs associated with secondary metabolism and stress response in Digitalis purpurea

Abstract

Conclusions

Through comprehensive transcriptome analysis, we not only identified 29 novel gene families potentially involved in the biosynthesis of cardiac glycosides but also characterized a large number of mlncRNAs. Our results suggest the importance of mlncRNAs in secondary metabolism and stress response in D. purpurea.

Analysis of BAC clones containing homologous sequences on the end of the Xq arm and on chromosome 7 in the dioecious plant Silene latifolia

Silene latifolia is a model dioecious plant with morphologically distinguishable XY sex chromosomes. The end of the Xq arm is quite different from that of the Yp arm, although both are located at opposite ends of their respective chromosomes relative to a pseudo-autosomal region. The Xq arm does not seem to originate from the same autosome as the Yp arm. Bacterial artificial chromosome clone #15B12 has an insert containing a 130-kb stretch in which a 313-bp satellite DNA is repeated 420 times. PCR with a single primer revealed that this 130-kb stretch consists of three reversals of the orientation of the satellite DNA. A non-long terminal repeat retroelement and two sequences that share homology with an Oryza sativa RING zinc finger and a putative Arabidopsis thaliana protein, respectively, were found in the sequences that flank the satellite DNA. Fluorescence in situ hybridization carried out using this low-copy region of #15B12 as a probe confirmed that these sequences originated from the X chromosome and that homologous sequences exist at the end of chromosome 7. The region distal to DD44X on the Xq arm is postulated to have recombined with a region containing satellite DNA on chromosome 7 during the process of sex chromosome evolution.

A SUPERMAN-like gene is exclusively expressed in female flowers of the dioecious plant Silene latifolia

To elucidate the mechanism(s) underlying dioecious flower development, the present study analyzed a SUPERMAN (SUP) homolog, SlSUP, which was identified in Silene latifolia. The sex of this plant is determined by heteromorphic X and Y sex chromosomes. It was revealed that SlSUP is a single-copy autosomal gene expressed exclusively in female flowers. Introduction of a genomic copy of SlSUP into the Arabidopsis thaliana sup (sup-2) mutant complemented the excess-stamen and infertile phenotypes of sup-2, and the overexpression of SlSUP in transgenic Arabidopsis plants resulted in reduced stamen numbers as well as the suppression of petal elongation. During the development of the female flower in S. latifolia, the expression of SlSUP is first detectable in whorls 2 and 3 when the normal expression pattern of the B-class flowering genes was already established and persisted in the stamen primordia until the ovule had matured completely. In addition, significant expression of SlSUP was detected in the ovules, suggestive of the involvement of this gene in ovule development. Furthermore, it was revealed that the de-suppression of stamen development by infection of the S. latifolia female flower with Microbotryum violaceum was accompanied by a significant reduction in SlSUP transcript levels in the induced organs. Taken together, these results demonstrate that SlSUP is a female flower-specific gene and suggest that SlSUP has a positive role in the female flower developmental pathways of S. latifolia.

Diamonds in the rough: mRNA-like non-coding RNAs

Non-coding RNAs are increasingly being identified as crucial regulators of gene expression and other cellular functions in plants. Experimental and computational methods have revealed the existence of mRNA-like non-coding RNAs (mlncRNAs), a class of non-coding RNAs that, in plants, are associated with tissue-specific expression, development and the phosphate-starvation response. Although their mechanisms of action are largely unknown, one can speculate that mlncRNAs act through secondary structures or specific sequences that bind to proteins or metabolites, or that have catalytic activity. This review summarizes the computational methods developed to identify candidate mlncRNAs, and the current experimental evidence regarding the function of several known mlncRNAs.

The Y chromosome-specific STS marker MS2 and its peripheral regions on the Y chromosome of the dioecious plant Silene latifolia

Sex determination in Silene latifolia uses the XX/XY system. The recent evolution of dioecy in S. latifolia provides a unique opportunity to study the early stages of Y chromosome evolution. However, the current Y chromosome map still contains many large gaps with no available markers. In this study, a sequence tagged site (STS) marker, MS2, was isolated and mapped to the same locus as L8 on the Y chromosome. To investigate the peripheral regions of MS2, a bacterial artificial chromosome (BAC) library was constructed from a male plant, and the BAC clone containing MS2 (MS2-9d12F) was isolated from 32 640 clones with an average insert size of 115 kb. A 109-kb insert of the BAC clone was analyzed. BLASTX analysis showed 11 sequences similar to some known proteins, most of which are retrotransposon-like elements. The ORF Finder predicted 9 ORFs within MS2-9d12F. RT-PCR analyses revealed that only 4 of the 9 predicted ORFs are expressed in both male and female plants. These 4 ORFs are candidates for genes having counterparts on both the X and Y chromosomes. Dot-matrix plot analysis and a BLASTN search revealed LTR-like sequences close to the retrotransposon-like elements and high similarity to 3 known genomic sequences of S. latifolia. These results suggest an accumulation of retrotransposons and segmental duplications in peripheral regions of MS2 during the early stage of sex chromosome evolution.

Plant sex chromosomes: molecular structure and function

Recent molecular and genomic studies carried out in a number of model dioecious plant species, including Asparagus officinalis, Carica papaya, Silene latifolia, Rumex acetosa and Marchantia polymorpha, have shed light on the molecular structure of both homomorphic and heteromorphic sex chromosomes, and also on the gene functions they have maintained since their evolution from a pair of autosomes. The molecular structure of sex chromosomes in species from different plant families represents the evolutionary pathway followed by sex chromosomes during their evolution. The degree of Y chromosome degeneration that accompanies the suppression of recombination between the Xs and Ys differs among species. The primitive Ys of A. officinalis and C. papaya have only diverged from their homomorphic Xs in a short male-specific and non-recombining region (MSY), while the heteromorphic Ys of S. latifolia, R. acetosa and M. polymorpha have diverged from their respective Xs. As in the Y chromosomes of mammals and Drosophila, the accumulation of repetitive DNA, including both transposable elements and satellite DNA, has played an important role in the divergence and size enlargement of plant Ys, and consequently in reducing gene density. Nevertheless, the degeneration process in plants does not appear to have reached the Y-linked genes. Although a low gene density has been found in the sequenced Y chromosome of M. polymorpha, most of its genes are essential and are expressed in the vegetative and reproductive organs in both male and females. Similarly, most of the Y-linked genes that have been isolated and characterized up to now in S. latifolia are housekeeping genes that have X-linked homologues, and are therefore expressed in both males and females. Only one of them seems to be degenerate with respect to its homologous region in the X. Sequence analysis of larger regions in the homomorphic X and Y chromosomes of papaya and asparagus, and also in the heteromorphic sex chromosomes of S. latifolia and R. acetosa, will reveal the degenerative changes that the Y-linked gene functions have experienced during sex chromosome evolution.

A population genomic analysis of species boundaries: neutral processes, adaptive divergence and introgression between two hybridizing plant species

Little is known about the nature of species boundaries between closely related plant species and about the extent of introgression as a consequence of hybridization upon secondary contact. To address these topics we analyzed genome-wide differentiation between two closely related Silene species, Silene latifolia and S. dioica, and assessed the strength of introgression in sympatry. More than 300 AFLP markers were genotyped in three allopatric and three sympatric populations of each species. Outlier analyses were performed separately for sympatric and allopatric populations. Both positive and negative outlier loci were found, indicating that divergent and balancing selection, respectively, have shaped patterns of divergence between the two species. Sympatric populations of the two species were found to be less differentiated genetically than allopatric populations, indicating that hybridization has led to gene introgression. We conclude that differentiation between S. latifolia and S. dioica has been shaped by a combination of introgression and selection. These results challenge the view that species differentiation is a genome-wide phenomenon, and instead support the idea that genomes can be porous and that species differentiation has a genic basis.

Floral development of an asexual and female-like mutant carrying two deletions in gynoecium-suppressing and stamen-promoting functional regions on the Y chromosome of the dioecious plant Silene latifolia

Sexual dimorphism is controlled by genes on the Y chromosome in the dioecious plant Silene latifolia. K034 is the first mutant with female flowers and asexual flowers in one individual. Its stamens are suppressed completely, and its gynoecium exhibits two suppression patterns. One gynoecium resembles a thin rod, as in wild-type males (asexual flower); the other is imperfectly suppressed, having 1-3 carpels (female-like flower). The ratio of these patterns was 9 : 1. To exclude the possibility of chimerism in K034, we crossed a female-like flower of K034 with a wild-type male. Progeny obtained from this crossing had asexual and female-like flowers in one individual. This two-flower-type phenotype was inherited without separating. To examine the identity of flower organs in K034, we analyzed the development of asexual and female-like flowers using scanning electron microscopy and in situ hybridization with SLM1 and SLM2 (orthologs of AGAMOUS and PISTILLATA, respectively) as probes. Mitotic spreads of root tip chromosomes from hairy root cultures showed that K034 had 25 chromosomes. Fluorescent in situ hybridization analysis, using a subtelomeric repetitive sequence (KpnI subfamily) as a probe, indicated that K034 possessed two X chromosomes and one Y chromosome (Y(d)), of which Y(d) had been rearranged to lose the pseudoautosomal region (PAR). PCR analysis using Y-specific sequence-tagged site (STS) markers clarified that Y(d) of K034 had two other deletions in gynoecium-suppressing and stamen-promoting regions. It is reasonable to suggest that these sex chromosomal abnormalities resulted in two abnormal sexual phenotypes: the asexual and imperfect female (female-like) flowers in K034.

Sex chromosome-linked genes in plants

Recent studies of plant sex chromosome-linked genes have revealed many interesting characteristics, although there are limited reports about heteromorphic sex chromosomes in flowering plants. Sex chromosome-linked genes in angiosperms have been characterized mainly in the dioecious plant Silene latifolia. Although all such genes were isolated from transcripts of male flower buds of S. latifolia, most seem to be housekeeping genes except for the petal- and stamen-specific MADS box gene on the Y chromosome (SlAP3Y) and the male reproductive organ-specific gene on the X chromosome (MROS3X). Recent evolutionary studies have revealed at least three evolutionary strata on the X chromosome that are related to stepwise loss of recombination between the sex chromosomes. Moreover, genetic maps showed conservation of gene organization on the X chromosome in the genus Silene and substantial pericentric inversion between the X and Y chromosomes of S. latifolia during evolution. A comparison between paralogs on the sex chromosomes revealed that introns of the Y-linked genes are longer than those of X-linked paralogs. Although analyses of sex chromosome-linked genes suggest that degeneration of the Y chromosome has occurred, the Y chromosome in flowering plants remains the largest in the male genome, unlike that of mammals. Accumulation of repetitive sequences and the entire chloroplast genome on the Y chromosome appear to have contributed to this large size. However, more detailed studies will be required to help explain the basis for the fact that heteromorphic sex chromosomes in angiosperms are large.

Expression of ASK1-like genes in arrested stamens of female Silene latifolia plants

Silene latifolia is a dioecious plant in which sex is determined by heteromorphic sex chromosomes. In female plants, stamen development is arrested before microspore mother cells are formed. In this study, we isolated four cDNAs (SlSKP1-1 to 4) encoding ASK1-like protein as expression markers to reveal when expression levels are reduced in arrested stamens of female flowers. Expression patterns of the SlSKP1 genes were analyzed by in-situ hybridization. We use the flower development classification of Grant et al. (in Plant J 6:471-480, 1994). SlSKP1 genes were highly expressed in primary parietal cells and primary sporogenous cells that develop into microspore mother cells in male flowers. Expression levels started to be reduced in the external stamens of the female flowers when stamen development was arrested at stage 7. Although microspore mother cells could not be developed in female flowers and SlSKP1 expression may be unnecessary in arrested stamens, SlSKP1 genes were still expressed in sporogenous cells of degenerated stamens at stage 8. Parietal cells stopped differentiating earlier than sporogenous cells in arrested stamens. These results suggest that not all types of cell are arrested simultaneously at a particular stage of stamen development during stamen suppression of S. latifolia.

An anther- and petal-specific gene SlMF1 is a multicopy gene with homologous sequences on sex chromosomes

A male flower-specific gene SlMF1 was isolated from male flower buds of the dioecious plant Silene latifolia. SlMF1 is expressed in all the floral meristems at the very early stage of development in both male and female flower buds. At the mature stage of development in male flower buds, SlMF1 transcripts were specifically accumulated in pollen mother cells, tapetal cells, and the developing tips of petals. Genomic Southern hybridization revealed that SlMF1 was a multicopy gene with a Y chromosome-linked homologous sequence. PCR analyses with flow-sorted chromosomes showed that SlMF1 was localized on both autosomes and the X chromosome.