Genome Sequence Resource for the Ramie Oomycete Pathogen Phytopythium vexans HF1.

The oomycete Phytopythium vexans is a causative agent of patch canker, damping-off, and crown, stem, and root rot in many economically important plants. P. vexans HF1 was isolated in China, where it caused brown root rot of ramie, a fiber crop broadly cultivated in Asia. The genome of HF1 was sequenced by a combination of technologies producing short (Illumina HiSeq X) and long (PacBio RS) reads. The genome is 41.73 Mbp long, assembled into 44 contigs. It has a GC content of 58.17% and contains 13,051 predicted coding genes, including 1,461 putative virulence genes and 220 putative antimicrobial resistance genes. This genome sequence provides a resource for determining the molecular mechanisms of disease development in this pathosystem.

The miRNAome of ramie (Boehmeria nivea L.): identification, expression, and potential roles of novel microRNAs in regulation of cadmium stress response.

BACKGROUND: MicroRNAs (miRNAs) regulate numerous crucial abiotic stress processes in plants. However, information is limited on their involvement in cadmium (Cd) stress response and tolerance mechanisms in plants, including ramie (Boehmeria nivea L.) that produces a number of economic valuable as an important natural fibre crop and an ideal crop for Cd pollution remediation. RESULTS: Four small RNA libraries of Cd-stressed and non-stressed leaves and roots of ramie were constructed. Using small RNA-sequencing, 73 novel miRNAs were identified. Genome-wide expression analysis revealed that a set of miRNAs was differentially regulated in response to Cd stress. In silico target prediction identified 426 potential miRNA targets that include several uptake or transport factors for heavy metal ions. The reliability of small RNA sequencing and the relationship between the expression levels of miRNAs and their target genes were confirmed by quantitative PCR (q-PCR). We showed that the expression patterns of miRNAs obtained by q-PCR were consistent with those obtained from small RNA sequencing. Moreover, we demonstrated that the expression of six randomly selected target genes was inversely related to that of their corresponding miRNAs, indicating that the miRNAs regulate Cd stress response in ramie. CONCLUSIONS: This study enriches the number of Cd-responsive miRNAs and lays a foundation for the elucidation of the miRNA-mediated regulatory mechanism in ramie during Cd stress.

Comparative transcriptomics provide insight into the morphogenesis and evolution of fistular leaves in Allium.

BACKGROUND: Fistular leaves frequently appear in Allium species, and previous developmental studies have proposed that the process of fistular leaf formation involves programmed cell death. However, molecular evidence for the role of programmed cell death in the formation of fistular leaf cavities has yet to be reported. RESULTS: In this study, we characterized the leaf transcriptomes of nine Allium species, including six fistular- and three solid-leaved species. In addition, we identified orthologous genes and estimated their Ka and Ks values, in order to ascertain their selective pattern. Phylogenetic analysis based on the transcriptomes revealed that A. tuberosum was the most ancestral among the nine species, and analysis of orthologous genes between A. tuberosum and the other eight species indicated that 149 genes were subject to positive selection; whereas >3000 had undergone purifying selection in each species. CONCLUSIONS: We found that many genes that are potentially related to programmed cell death either exhibited rapid diversification in fistular-leaved species, or were conserved in solid-leaved species in evolutionary history. These genes potentially involved in programmed cell death might play important roles in the formation of fistular leaf cavities in Allium, and the differing selection patterns in fistular- and solid-leaved species may be responsible for the evolution of fistular leaves.

Identification of Ramie Genes in Response to Pratylenchus coffeae Infection Challenge by Digital Gene Expression Analysis.

Root lesion disease, caused by Pratylenchus coffeae, seriously impairs the growth and yield of ramie, an important natural fiber crop. The ramie defense mechanism against P. coffeae infection is poorly understood, which hinders efforts to improve resistance via breeding programs. In this study, the transcriptome of the resistant ramie cultivar Qingdaye was characterized using Illumina sequence technology. About 46.3 million clean pair end (PE) reads were generated and assembled into 40,826 unigenes with a mean length of 830 bp. Digital gene expression (DGE) analysis was performed on both the control roots (CK) and P. coffeae-challenged roots (CH), and the differentially expressed genes (DEGs) were identified. Approximately 10.16 and 8.07 million cDNA reads in the CK and CH cDNA libraries were sequenced, respectively. A total of 137 genes exhibited different transcript abundances between the two libraries. Among them, the expressions of 117 and 20 DEGs were up- and down-regulated in P. coffeae-challenged ramie, respectively. The expression patterns of 15 candidate genes determined by qRT-PCR confirmed the results of DGE analysis. Time-course expression profiles of eight defense-related genes in susceptible and resistant ramie cultivars were different after P. coffeae inoculation. The differential expression of protease inhibitors, pathogenesis-related proteins (PRs), and transcription factors in resistant and susceptible ramie during P. coffeae infection indicated that cystatin likely plays an important role in nematode resistance.

De novo assembly and characterization of transcriptome using Illumina paired-end sequencing and identification of CesA gene in ramie (Boehmeria nivea L. Gaud).

BACKGROUND: Ramie fiber, extracted from vegetative organ stem bast, is one of the most important natural fibers. Understanding the molecular mechanisms of the vegetative growth of the ramie and the formation and development of bast fiber is essential for improving the yield and quality of the ramie fiber. However, only 418 expressed tag sequences (ESTs) of ramie deposited in public databases are far from sufficient to understand the molecular mechanisms. Thus, high-throughput transcriptome sequencing is essential to generate enormous ramie transcript sequences for the purpose of gene discovery, especially genes such as the cellulose synthase (CesA) gene. RESULTS: Using Illumina paired-end sequencing, about 53 million sequencing reads were generated. De novo assembly yielded 43,990 unigenes with an average length of 824 bp. By sequence similarity searching for known proteins, a total of 34,192 (77.7%) genes were annotated for their function. Out of these annotated unigenes, 16,050 and 13,042 unigenes were assigned to gene ontology and clusters of orthologous group, respectively. Searching against the Kyoto Encyclopedia of Genes and Genomes Pathway database (KEGG) indicated that 19,846 unigenes were mapped to 126 KEGG pathways, and 565 genes were assigned to http://starch and sucrose metabolic pathway which was related with cellulose biosynthesis. Additionally, 51 CesA genes involved in cellulose biosynthesis were identified. Analysis of tissue-specific expression pattern of the 51 CesA genes revealed that there were 36 genes with a relatively high expression levels in the stem bark, which suggests that they are most likely responsible for the biosynthesis of bast fiber. CONCLUSION: To the best of our knowledge, this study is the first to characterize the ramie transcriptome and the substantial amount of transcripts obtained will accelerate the understanding of the ramie vegetative growth and development mechanism. Moreover, discovery of the 36 CesA genes with relatively high expression levels in the stem bark will present an opportunity to understand the ramie bast fiber formation and development mechanisms.

Identification of drought stress-responsive transcription factors in ramie (Boehmeria nivea L. Gaud).

BACKGROUND: Ramie fiber extracted from stem bark is one of the most important natural fibers. Drought is a main environment stress which severely inhibits the stem growth of ramie and leads to a decrease of the fiber yield. The drought stress-regulatory mechanism of ramie is poorly understood. RESULT: Using Illumina sequencing, approximately 4.8 and 4.7 million (M) 21-nt cDNA tags were respectively sequenced in the cDNA libraries derived from the drought-stressed ramie (DS) and the control ramie under well water condition (CO). The tags generated from the two libraries were aligned with ramie transcriptome to annotate their function and a total of 23,912 and 22,826 ramie genes were matched by these tags of DS and CO library, respectively. Comparison of gene expression level between CO and DS ramie based on the differences of tag frequencies appearing in the two libraries revealed that there were 1516 potential drought stress-responsive genes, in which 24 genes function as transcription factor (TF). Among these 24 TFs, the unigene19721 encoding the DELLA protein which is a key negative regulator in gibberellins (GAs) signal pathway was probably markedly up-regulated under water stress for a increase of tag abundance in DS library, which is possibly responsible for the inhibition of the growth of drought-stressed ramie. In order to validate the change of expression of these potential stress-responsive TFs under water deficit condition, the unigene19721 and another eleven potential stress-responsive TFs were chosen for further expression analysis in well-watered and drought-stressed ramie by real-time quantitative PCR (qRT-PCR) and the result showed that all 12 TFs were authentically involved in the response of drought stress. CONCLUSION: In this study, twelve TFs involving in the response of drought stress were first found by Illumina tag-sequencing and qRT-PCR in ramie. The discovery of these drought stress-responsive TFs will be helpful for further understanding the drought stress-regulatory mechanism of ramie and improving the drought tolerance ability of ramie.

Development of a Loop-Mediated Isothermal Amplification Method for the Rapid Detection of Phytopythium vexans.

Brown root rot caused by Phytopythium vexans is a new destructive root disease on many plants such as Gingko, Citrus, kiwifruit, and ramie. The establishment of loop-mediated isothermal amplification (LAMP) technology for detecting P. vexans can help monitor and control brown root rot quickly, efficiently, and accurately. LAMP technology is known for its simplicity, sensitivity, and speed; and it does not require any specialized equipment - a water bath or a thermoblock is sufficient for isothermal amplifications. LAMP products can be visualized by using hydroxy naphthol blue (HNB) dye or agarose gel electrophoresis. In this study, by searching and comparing the internal transcribed spacer (ITS) sequences of P. vexans and the related species in oomycete genera Pythium, Phytopythium, and Phytophthora, we designed specific primers targeting the ITS gene region of P. vexans. Using HNB dye, we established a LAMP technique for rapid detection of P. vexans by visible color change. In addition, we optimized the protocol to enhance both sensitivity and specificity for P. vexans detection. Under the optimized condition, our protocol based on LAMP technology could detect as low as 24 copies of the P. vexans genomic DNA, which is ∼100 times more sensitive than conventional PCR. This method can successfully detect P. vexans using cell suspensions from P. vexans - infected ramie root tissues.

Selection of Reference Genes for qPCR Analyses of Gene Expression in Ramie Leaves and Roots across Eleven Abiotic/Biotic Treatments.

Quantitative real-time PCR (qPCR) is commonly used for deciphering gene functions. For effective qPCR analyses, suitable reference genes are needed for normalization. The objective of this study is to identify the appropriate reference gene(s) for qPCR analyses of the leaves and roots of ramie (Boehmeria nivea L.), an important natural fiber crop. To accomplish this goal, we investigated the expression patterns of eight common plant qPCR reference genes in ramie leaves and roots under five abiotic stresses, five hormonal treatments, and one biotic stress. The relative expression stabilities of the eight genes were evaluated using four common but different approaches: geNorm, NormFinder, BestKeeper, and RefFinder. Across the 11 tested conditions, ACT1 was the most stably expressed among the eight genes while GAPDH displayed the biggest variation. Overall, while variations in the suggested reference genes were found for different tissue x treatment combinations, our analyses revealed that together, genes ACT1, CYP2, and UBQ can provide robust references for gene expression studies of ramie leaves under most conditions, while genes EF-1α, TUB, and ACT1 can be used for similar studies of ramie roots. Our results should help future functional studies of the genes in ramie genome across tissues and environmental conditions.

Deletion of the sex-determining gene SXI1α enhances the spread of mitochondrial introns in Cryptococcus neoformans.

BACKGROUND: Homing endonuclease genes (HEGs) are widely distributed genetic elements in the mitochondrial genomes of a diversity of eukaryotes. Due to their ability to self-propagate within and between genomes, these elements can spread rapidly in populations. Whether and how such elements are controlled in genomes remains largely unknown. RESULTS: Here we report that the HEG-containing introns in the mitochondrial COX1 gene in Cryptococcus neoformans are mobile and that their spread in sexual crosses is influenced by mating type (MAT) α-specific homeodomain gene SXI1α. C. neoformans has two mating types, MATa and MATα . In typical crosses between strains of the two mating types, only a small portion (< 7%) of diploid fusants inherited the HEGs from the MATα parent. However, disruption of the SXI1α gene resulted in the majority (> 95%) of the diploid fusants inheriting the HEG-containing introns from the MATα parent, a frequency significantly higher than those of intronless mitochondrial genes. CONCLUSIONS: Our results suggest that SXI1α not only determines uniparental mitochondrial inheritance but also inhibits the spread of HEG-containing introns in the mitochondrial genome in C. neoformans.

Genome-Wide Expression Profiles of Hemp (Cannabis sativa L.) in Response to Drought Stress.

Drought is the main environmental factor impairing hemp growth and yield. In order to decipher the molecular responses of hemp to drought stress, transcriptome changes of drought-stressed hemp (DS1 and DS2), compared to well-watered control hemp (CK1 and CK2), were studied with RNA-Seq technology. RNA-Seq generated 9.83, 11.30, 11.66, and 11.31 M clean reads in the CK1, CK2, DS1, and DS2 libraries, respectively. A total of 1292 differentially expressed genes (DEGs), including 409 (31.66%) upregulated and 883 (68.34%) downregulated genes, were identified. The expression patterns of 12 selected genes were validated by qRT-PCR, and the results were accordant with Illumina analysis. Gene Ontology (GO) and KEGG analysis illuminated particular important biological processes and pathways, which enriched many candidate genes such as NAC, B3, peroxidase, expansin, and inositol oxygenase that may play important roles in hemp tolerance to drought. Eleven KEGG pathways were significantly influenced, the most influenced being the plant hormone signal transduction pathway with 15 differentially expressed genes. A similar expression pattern of genes involved in the abscisic acid (ABA) pathway under drought, and ABA induction, suggested that ABA is important in the drought stress response of hemp. These findings provide useful insights into the drought stress regulatory mechanism in hemp.

Molecular Cloning, Recombinant Expression and Antifungal Activity of BnCPI, a Cystatin in Ramie (Boehmeria nivea L.).

Phytocystatins play multiple roles in plant growth, development and resistance to pests and other environmental stresses. A ramie (Boehmeria nivea L.) phytocystatin gene, designated as BnCPI, was isolated from a ramie cDNA library and its full-length cDNA was obtained by rapid amplification of cDNA ends (RACE). The full-length cDNA sequence (691 bp) consisted of a 303 bp open reading frame (ORF) encoding a protein of 100 amino acids with deduced molecular mass of 11.06 kDa and a theoretical isoelectric point (pI) of 6.0. The alignment of genome DNA (accession no. MF153097) and cDNA sequences of BnCPI showed that an intron (~104 bp) exists in the coding region. The BnCPI protein contains most of the highly conserved blocks including Gly5-Gly6 at the N-terminal, the reactive site motif QxVxG (Q49V50V51S52G53), the L79-W80 block and the [LVI]-[AGT]-[RKE]-[FY]-[AS]-[VI]-x-[EDQV]-[HYFQ]-N (L22G23R24 F25A26V27 D28D29H30 N31) block that is common among plant cystatins. BLAST analysis indicated that BnCPI is similar to cystatins from Glycine max (77%), Glycine soja (76%), Hevea brasiliensis (75%) and Ricinus communis (75%). The BnCPI was subcloned into expression vector pSmart-I and then overexpressed in Escherichia coli BL21 (DE3) as a His-tagged recombinant protein. The purified reBnCPI has a molecular mass of 11.4 kDa determined by sodium dodecyl sulfate polyacrylamide gel electrophoresis (SDS-PAGE). Purified reBnCPI can efficiently inhibit the protease activity of papain and ficin toward BANA (Nα-benzoyl-L-arginine-2-naphthyamide), as well as the mycelium growth of some important plant pathogenic fungi. The data further contribute to our understanding of the molecular functions of BnCPI.

Transcriptome Analysis of Ramie (Boehmeria nivea L. Gaud.) in Response to Ramie Moth (Cocytodes coerulea Guenée) Infestation.

The ramie moth Cocytodes coerulea Guenée (RM) is an economically important pest that seriously impairs the yield of ramie, an important natural fiber crop. The molecular mechanisms that underlie the ramie-pest interactions are unclear up to date. Therefore, a transcriptome profiling analysis would aid in understanding the ramie defense mechanisms against RM. In this study, we first constructed two cDNA libraries derived from RM-challenged (CH) and unchallenged (CK) ramie leaves. The subsequent sequencing of the CH and CK libraries yielded 40.2 and 62.8 million reads, respectively. Furthermore, de novo assembling of these reads generated 26,759 and 29,988 unigenes, respectively. An integrated assembly of data from these two libraries resulted in 46,533 unigenes, with an average length of 845 bp per unigene. Among these genes, 24,327 (52.28%) were functionally annotated by predicted protein function. A comparative analysis of the CK and CH transcriptome profiles revealed 1,980 differentially expressed genes (DEGs), of which 750 were upregulated and 1,230 were downregulated. A quantitative real-time PCR (qRT-PCR) analysis of 13 random selected genes confirmed the gene expression patterns that were determined by Illumina sequencing. Among the DEGs, the expression patterns of transcription factors, protease inhibitors, and antioxidant enzymes were studied. Overall, these results provide useful insights into the defense mechanism of ramie against RM.

Development and characterization of 1,827 expressed sequence tag-derived simple sequence repeat markers for ramie (Boehmeria nivea L. Gaud).

Ramie (Boehmeria nivea L. Gaud) is one of the most important natural fiber crops, and improvement of fiber yield and quality is the main goal in efforts to breed superior cultivars. However, efforts aimed at enhancing the understanding of ramie genetics and developing more effective breeding strategies have been hampered by the shortage of simple sequence repeat (SSR) markers. In our previous study, we had assembled de novo 43,990 expressed sequence tags (ESTs). In the present study, we searched these previously assembled ESTs for SSRs and identified 1,685 ESTs (3.83%) containing 1,878 SSRs. Next, we designed 1,827 primer pairs complementary to regions flanking these SSRs, and these regions were designated as SSR markers. Among these markers, dinucleotide and trinucleotide repeat motifs were the most abundant types (36.4% and 36.3%, respectively), whereas tetranucleotide, pentanucleotide, and hexanucleotide motifs represented <10% of the markers. The motif AG/CT was the most abundant, accounting for 28.74% of the markers. One hundred EST-SSR markers (97 SSRs located in genes encoding transcription factors and 3 SSRs in genes encoding cellulose synthases) were amplified using polymerase chain reaction for detecting 24 ramie varieties. Of these 100 markers, 98 markers were successfully amplified and 81 markers were polymorphic, with 2-6 alleles among the 24 varieties. Analysis of the genetic diversity of all 24 varieties revealed similarity coefficients that ranged from 0.51 to 0.80. The EST-SSRs developed in this study represent the first large-scale development of SSR markers for ramie. These SSR markers could be used for development of genetic and physical maps, quantitative trait loci mapping, genetic diversity studies, association mapping, and cultivar fingerprinting.