Search

Abstract The molecular characterization of genetically modified organisms (GMOs) is essential for ensuring safety and gaining regulatory approval for commercialization. According to CODEX standards, this characterization involves evaluating the presence of introduced genes, insertion sites, copy number, and nucleotide sequence structure. Advances in technology have led to the increased use of next-generation sequencing (NGS) over traditional methods such as Southern blotting. While both methods provide high reproducibility and accuracy, Southern blotting is labor-intensive and time-consuming due to the need for repetitive probe design and analyses for each target, resulting in low throughput. Conversely, NGS facilitates rapid and comprehensive analysis by mapping whole-genome sequencing (WGS) data to plasmid sequences, accurately identifying T-DNA insertion sites and flanking regions. This advantage allows for efficient detection of T-DNA presence, copy number, and unintended gene insertions without additional probe work. This paper reviews the current status of GMO genome characterization using NGS and proposes more efficient strategies for this purpose.

Inorganic phosphate (Pi) homeostasis plays an important role in plant growth and abiotic stress tolerance. Several MYB-CC transcription factors involved in Pi homeostasis have been identified in rice (Oryza sativa). PHOSPHATE STARVATION RESPONSE-LIKE 7 (PHL7) is a class II MYC-CC protein, in which the MYC-CC domain is located at the N terminus. In this study, we established that OsPHL7 is localized to the nucleus and that the encoding gene is induced by Pi deficiency. The Pi-responsive genes and Pi transporter genes are positively regulated by OsPHL7. The overexpression of OsPHL7 enhanced the tolerance of rice plants to Pi starvation, whereas the RNA interference-based knockdown of this gene resulted in increased sensitivity to Pi deficiency. Transgenic rice plants overexpressing OsPHL7 produced more roots than wild-type plants under both Pi-sufficient and Pi-deficient conditions and accumulated more Pi in the shoots and roots. In addition, the overexpression of OsPHL7 enhanced rice tolerance to salt stress. Together, these results demonstrate that OsPHL7 is involved in the maintenance of Pi homeostasis and enhances tolerance to Pi deficiency and salt stress in rice.

Abstract Faba bean (Vicia faba L.), a globally important grain legume providing a stable source of dietary protein, was one of the earliest plant cytogenetic models. However, the lack of draft genome annotations and unclear structural information on mRNA transcripts have impeded its genetic improvement. To address this, we sequenced faba bean leaf transcriptome using the PacBio single-molecule long-read isoform sequencing platform. We identified 28,569 nonredundant unigenes, ranging from 108 to 9669 bp, with a total length of 94.5 Mb. Many unigenes (3597, 12.5%) had 2–20 isoforms, indicating a highly complex transcriptome. Approximately 96.5% of the unigenes matched sequences in public databases. The predicted proteins and transcription factors included NB-ARC, Myb_domain, C3H, bHLH, and heat shock proteins, implying that this genome has an abundance of stress resistance genes. To validate our results, we selected WCOR413-15785, DHN2-12403, DHN2-14197, DHN2-14797, COR15-14478, and HVA22-15 unigenes from the ICE-CBF-COR pathway to analyze their expression patterns in cold-treated samples via qRT-PCR. The expression of dehydrin-related genes was induced by cold stress. The assembled data provide the first insights into the deep sequencing of full-length RNA from faba bean at the single-molecule level. This study provides an important foundation to improve gene modeling and protein prediction.

In cultivated plants, shoot morphology is an important factor that influences crop economic value. However, the effects of gene expression patterns on shoot morphology are not clearly understood. In this study, the molecular mechanism behind shoot morphology (including leaf, stem, and node) was analyzed using RNA sequencing to compare weedy (creeper) and cultivar (stand) growth types obtained in F7 derived from a cross of wild and cultivated soybeans. A total of 12,513 (in leaves), 14,255 (in stems), and 11,850 (in nodes) differentially expressed genes were identified among weedy and cultivar soybeans. Comparative transcriptome and expression analyses revealed 22 phytohormone-responsive genes. We found that GIBBERELLIN 2-OXIDASE 8 (GA2ox), SPINDLY (SPY), FERONIA (FER), AUXIN RESPONSE FACTOR 8 (ARF8), CYTOKININ DEHYDROGENASE-1 (CKX1), and ARABIDOPSIS HISTIDINE KINASE-3 (AHK3), which are crucial phytohormone response genes, were mainly regulated in the shoot of weedy and cultivar types. These results indicate that interactions between phytohormone signaling genes regulate shoot morphology in weedy and cultivar growth type plants. Our study provides insights that are useful for breeding and improving crops to generate high-yield soybean varieties.

Abstract Plants have evolved sophisticated defense systems to enhance drought tolerance. These include the microRNA (miRNA) group of small noncoding RNAs that act as post‐transcriptional regulators; however, details of the mechanisms by which they confer drought tolerance are not well understood. Here, we show that osa‐MIR171f, a member of osa‐MIR171 gene family, is mainly expressed in response to drought stress and regulates the transcript levels of SCARECROW‐LIKE6‐I (SCL6‐I) and SCL6‐II in rice (Oryza sativa). The SCL6 genes are known to be involved in shoot branching and flag leaf morphology. Osa‐MIR171f‐overexpressing (osa‐MIR171f‐OE) transgenic plants showed reduced drought symptoms compared with non‐transgenic (NT) control plants under both field drought and polyethylene glycol (PEG)‐mediated dehydration stress conditions. Transcriptome analysis of osa‐MIR171f‐OE plants and osa‐mir171f‐knockout (K/O) lines generated by clustered regularly interspaced short palindromic repeats (CRISPR/Cas9) revealed that osa‐mature‐miR171a‐f (osa‐miR171) regulates the expression of flavonoid biosynthesis genes, consequently leading to drought tolerance. This upregulation in the osa‐MIR171f‐OE plants, which did not occur in NT control plants, was observed under both normal and drought conditions. Our findings indicate that osa‐miR171 plays a role in drought tolerance by regulating SCL6‐I and SCL6‐II transcript levels.

Caulophyllum robustum, commonly named Asian blue cohosh, is a perennial herb in the family Berberidaceae. It has traditionally been used for folk medicine in China. We isolated berberine from the leaves, stem, roots, and fruits of C. robustum, and this is the first report on berberine in this species. Transcriptome analysis was conducted for the characterization of berberine biosynthesis genes in C. robustum, in which, all the genes for berberine biosynthesis were identified. From 40,094 transcripts, using gene ontology (GO) analysis, 26,750 transcripts were assigned their functions in the categories of biological process, molecular function, and cellular component. In the analysis of genes expressed in different tissues, the numbers of genes in the categories of intrinsic component of membrane and transferase activity were up-regulated in leaves versus stem. The berberine synthesis genes in C. robustum were characterized by phylogenetic analysis with corresponding genes from other berberine-producing species. The co-existence of genes from different plant families in the deepest branch subclade implies that the differentiation of berberine synthesis genes occurred early in the evolution of berberine-producing plants. Furthermore, the copy number increment of the berberine synthesis genes was detected at the species level.

Berberine and berbamine are bioactive compounds of benzylisoquinoline alkaloids (BIAs) present in Berberis species. The contents of berbamine are 20 times higher than berberine in leaf tissues in three closely related species: Berberis koreana, B. thunbergii and B. amurensis. This is the first report on the quantification of berberine compared to the berbamine in the Berberis species. Comparative transcriptome analyses were carried out with mRNAs from the leaf tissues of the three-species. The comparison of the transcriptomes of B. thunbergii and B. amurensis to those of B. koreana, B. thunbergii showed a consistently higher number of differentially expressed genes than B. amurensis in KEGG and DEG analyses. All genes encoding enzymes involved in berberine synthesis were identified and their expressions were variable among the three species. There was a single copy of CYP80A/berbamunine synthase in B. koreana. Methyltransferases and cytochrome P450 mono-oxidases (CYPs) are key enzymes for BIA biosynthesis. The current report contains the copy numbers and other genomic characteristics of the methyltransferases and CYPs in Berberis species. Thus, the contents of the current research are valuable for molecular characterization for the medicinal utilization of the Berberis species.

Short crop height is the preferred breeding trait since there is a positive correlation between lodging resistance and a crop yield increase. Alternative splicing can alter transcriptome diversity and contribute to plant adaptation to environmental stress. We characterized the transcriptomes obtained from dwarf and normal soybean lines derived from a cross of Glycine max var. Peking (G. max) and G. soja var. IT182936 in an F7 RIL population to study the differences between the isoforms. Full-length mRNA derived from leaf tissues was sequenced using the PacBio RSII platform, generating 904,474 circular consensus sequence (CCS) reads. Using the Structural and Quality Annotation of Novel Transcript Isoforms (SQANTI) process, 42,582 and 44,762 high-quality isoforms, and 91 and 179 polished low-quality isoforms were obtained in dwarf and normal cells, respectively. As a result, 832 and 36,772 nonredundant transcripts were generated. Approximately 30% of the identified genes were estimated to produce two or more isoforms. We detected an average of 166,171 splice junctions (SJs), of which 93.8% were canonical SJs. We identified that novel isoforms accounted for 19% of all isoforms, among which 12% fell within coding regions. The dwarf soybean demonstrated a greater number of isoforms in most of the annotated genes, particularly in genes related to growth hormones and defense responses. Our study provides comprehensive isoform and gene information that may accelerate transcriptome research in G. max and provide a basis to further study the impact of these isoforms on plant growth.

The annual herb Euphorbia maculata L. produces anti-inflammatory and biologically active substances such as triterpenoids, tannins, and polyphenols, and it is used in traditional Chinese medicine. Of these bioactive compounds, terpenoids, also called isoprenoids, are major secondary metabolites in E. maculata. Full-length cDNA sequencing was carried out to characterize the transcripts of terpenoid biosynthesis reference genes and determine the copy numbers of their isoforms using PacBio SMRT sequencing technology. The Illumina short-read sequencing platform was also employed to identify differentially expressed genes (DEGs) in the secondary metabolite pathways from leaves, roots, and stems. PacBio generated 62 million polymerase reads, resulting in 81,433 high-quality reads. From these high-quality reads, we reconstructed a genome of 20,722 genes, in which 20,246 genes (97.8%) did not have paralogs. About 33% of the identified genes had two or more isoforms. DEG analysis revealed that the expression level differed among gene paralogs in the leaf, stem, and root. Whole sets of paralogs and isoforms were identified in the mevalonic acid (MVA), methylerythritol phosphate (MEP), and terpenoid biosynthesis pathways in the E. maculata L. The nucleotide information will be useful for identifying orthologous genes in other terpenoid-producing medicinal plants.

Maize seed pigmentation is one of the important issue to develop maize seed breeding. The differently gene expression was characterized and compared for three inbred lines, such as the pigment accumulated seed (CM22) and non-pigmented seed (CM5 and CM19) at 10 days after pollination. We obtained a total of 63,870, 82,496, and 54,555 contigs by de novo assembly to identify gene expression in the CM22, CM5, and CM19, respectably. In differentially expressed gene analysis, it was revealed that 7,044 genes were differentially expressed by at least two-fold, with 4,067 upregulated in colored maize inbred lines and 2,977 upregulated in colorless maize inbred lines. Of them,18 genes were included to the anthocyanin biosynthesis pathways, while 15 genes were upregulated in both CM22/5 and CM22/19. Additionally, 37 genes were detected in the metabolic pathway concern to the seed pigmentation by BINs analysis using MAPMAN software. Finally, these differently expressed genes may aid in the research on seed pigmentation in maize breeding programs.

In this study, we investigated changes in the isoflavone content, total phenolic content (TPC), total flavonoid content (TFC), antioxidant activities (DPPH, ABTS), and anti-inflammatory activities of small-seeded and large-seeded soybean cultivars during germination (light/dark conditions). Total isoflavone content was higher at the seed stage in large-seeded soybeans, while it increased after 7 days of germination in small-seeded soybeans, particularly in response to light conditions, under which they had high TPC, TFC, and antioxidant activities. In large-seeded soybeans, the germination environment did not significantly affect TFC or DPPH inhibition, whereas TPC and ABTS inhibition were high under dark germination conditions. Extracts of sprouts exhibited superior anti-inflammatory activities. Nitric oxide production was slightly lower in small-seeded and large-seeded soybeans germinated under light and dark conditions, respectively. Our findings indicate that germinated soybeans improved nutritionally, and that enhancement of bioactivity under different germination environments could contribute to the selection of appropriate soybean cultivars.

Abstract Background Genetically modified crops (GM crops) have been developed to improve the agricultural traits of modern crop cultivars. Safety assessments of GM crops are of paramount importance in research at developmental stages and before releasing transgenic plants into the marketplace. Sequencing technology is developing rapidly, with higher output and labor efficiencies, and will eventually replace existing methods for the molecular characterization of genetically modified organisms. Methods To detect the transgenic insertion locations in the three GM rice gnomes, Illumina sequencing reads are mapped and classified to the rice genome and plasmid sequence. The both mapped reads are classified to characterize the junction site between plant and transgene sequence by sequence alignment. Results Herein, we present a next generation sequencing (NGS)-based molecular characterization method, using transgenic rice plants SNU-Bt9–5, SNU-Bt9–30, and SNU-Bt9–109. Specifically, using bioinformatics tools, we detected the precise insertion locations and copy numbers of transfer DNA, genetic rearrangements, and the absence of backbone sequences, which were equivalent to results obtained from Southern blot analyses. Conclusion NGS methods have been suggested as an effective means of characterizing and detecting transgenic insertion locations in genomes. Our results demonstrate the use of a combination of NGS technology and bioinformatics approaches that offers cost- and time-effective methods for assessing the safety of transgenic plants.

In this study, we report the complete chloroplast (cp) genome of Sophora koreensis and its relation with other species within the Fabaceae family. The cp genome was 154,870 bp long, with a typical quadripartite structure including a pair of inverted repeat regions (25,866 bp) separated by a large (85,037 bp) and small (18,101 bp) single-copy (SC) region. The genome encodes a total of 84 protein-coding genes, 35 tRNA genes, and 8 rRNA genes. Phylogenetic analysis suggested that S. koreensis is closely related to genus Sophora alopecuroides var. alopecuroides within Fabaceae.

Woodwardia japonica is one of the diverse members of the fern group and medicinally important genus. In Korea, the natural resources of W. japonica are being exhausted by excessive exploitation and require urgent conservation. In this study, the complete chloroplast genome of W. japonica was generated, and its structure was compared with that of other members of same family. The chloroplast genome was 153224 bp long, with a typical quadripartite structure including a pair of inverted repeat regions (24591 bp) separated by a large (82480 bp) and small (21562 bp) single-copy (SC) region. The genome encodes a total of 88 protein-coding genes, 35 tRNA genes, and eight rRNA genes. Additionally we identified 87 RNA editing sites in 52 genes; most of the substitution was U to C (50 sites), while C to U conversion occurred in 37 positions. The phylogenetic analysis strongly supported the relationship of W. japonica with W. unigemmata and. A. melanocaulon (Blechnoideae).

In this study, we report the complete chloroplast (cp) genome of Thelypteris interrupta, a fern member, and comparative analysis with its related family members. The cp genome was 155,983 bp long, with a typical quadripartite structure including a pair of inverted repeat regions (25,614 bp) separated by a large (82,769 bp) and small (21,986 bp) single-copy (SC) region. The genome encodes a total of 88 protein-coding genes, 35 tRNA genes, and 8 rRNA genes. Additionally, we identified 86 RNA editing sites in 52 genes; most of the substitution was U to C (52 sites), while C to U conversion occurred in 34 positions. The phylogenetic analysis strongly supported the relationship of T. interrupta with Ampelopteris prolifera and Christella appendiculata of Thelypteridoideae family.

Garden asparagus (Asparagus officinalis L.) is a popular vegetable cultivated worldwide. The secondary metabolites in its shoot are helpful for human health. We analyzed A. officinalis transcriptomes and identified differentially expressed genes (DEGs) involved in the biosynthesis of rutin and protodioscin, which are health-promoting functional compounds, and determined their association with stem color. We sequenced the complete mRNA transcriptome using the Illumina high-throughput sequencing platform in one white, three green, and one purple asparagus cultivars. A gene set was generated by de novo assembly of the transcriptome sequences and annotated using a BLASTx search. To investigate the relationship between the contents of rutin and protodioscin and their gene expression levels, rutin and protodioscin were analyzed using high-performance liquid chromatography. A secondary metabolite analysis using high-performance liquid chromatography showed that the rutin content was higher in green asparagus, while the protodioscin content was higher in white asparagus. We studied the genes associated with the biosynthesis of the rutin and protodioscin. The transcriptomes of the five cultivars generated 336 599 498 high-quality clean reads, which were assembled into 239 873 contigs with an average length of 694 bp, using the Trinity v2.4.0 program. The green and white asparagus cultivars showed 58 932 DEGs. A comparison of rutin and protodioscin biosynthesis genes revealed that 12 of the 57 genes associated with rutin and two of the 50 genes associated with protodioscin showed more than four-fold differences in expression. These DEGs might have caused a variation in the contents of these two metabolites between green and white asparagus. The present study is possibly the first to report transcriptomic gene sets in asparagus. The DEGs putatively involved in rutin and protodioscin biosynthesis might be useful for molecular engineering in asparagus.

Berberis koreana is a medicinal plant containing berberine, which is a bioactive compound of the benzylisoquinoline alkaloid (BIA) class. BIA is widely used in the food and drug industry for its health benefits. To investigate the berberine biosynthesis pathway, gene expression analysis was performed in leaves, flowers, and fruits at different stages of growth. This was followed by full-length cDNA sequencing analysis using the PacBio sequencer platform to determine the number of isoforms of those expressed genes. We identified 23,246 full-length unigenes, among which 8479 had more than one isoform. The number of isoforms ranged between two to thirty-one among all genes. Complete isoform analysis was carried out on the unigenes encoding BIA synthesis. Thirteen of the sixteen genes encoding enzymes for berberine synthesis were present in more than one copy. This demonstrates that gene duplication and translation into isoforms may contribute to the functional specificity of the duplicated genes and isoforms in plant alkaloid synthesis. Our study also demonstrated the streamlining of berberine biosynthesis via the absence of genes for enzymes of other BIAs, but the presence of all the genes for berberine biosynthesize in B. koreana. In addition to genes encoding enzymes for the berberine biosynthesis pathway, the genes encoding enzymes for other BIAs were not present in our dataset except for those encoding corytuberine synthase (CTS) and berbamunine synthase (BS). Therefore, this explains how B. koreana produces berberine by blocking the pathways leading to other BIAs, effectively only allowing the pathway to lead to berberine synthesis.

Kenaf is a source of fiber and a bioenergy crop that is considered to be a third world crop. Recently, a new kenaf cultivar, "Jangdae," was developed by gamma irradiation. It exhibited distinguishable characteristics such as higher biomass, higher seed yield, and earlier flowering than the wild type. We sequenced and analyzed the transcriptome of apical leaf and stem using Pacific Biosciences single-molecule long-read isoform sequencing platform. De novo assembly yielded 26,822 full-length transcripts with a total length of 59 Mbp. Sequence similarity against protein sequence allowed the functional annotation of 11,370 unigenes. Among them, 10,100 unigenes were assigned gene ontology terms, the majority of which were associated with the metabolic and cellular process. The Kyoto encyclopedia of genes and genomes (KEGG) analysis mapped 8875 of the annotated unigenes to 149 metabolic pathways. We also identified the majority of putative genes involved in cellulose and lignin-biosynthesis. We further evaluated the expression pattern in eight gene families involved in lignin-biosynthesis at different growth stages. In this study, appropriate biotechnological approaches using the information obtained for these putative genes will help to modify the desirable content traits in mutants. The transcriptome data can be used as a reference dataset and provide a resource for molecular genetic studies in kenaf.

Cirsium japonicum belongs to the Asteraceae or Compositae family and is a medicinal plant in Asia that has a variety of effects, including tumour inhibition, improved immunity with flavones, and antidiabetic and hepatoprotective effects. Silymarin is synthesized by 4-coumaroyl-CoA via both the flavonoid and phenylpropanoid pathways to produce the immediate precursors taxifolin and coniferyl alcohol. Then, the oxidative radicalization of taxifolin and coniferyl alcohol produces silymarin. We identified the expression of genes related to the synthesis of silymarin in C. japonicum in three different tissues, namely, flowers, leaves, and roots, through RNA sequencing. We obtained 51,133 unigenes from transcriptome sequencing by de novo assembly using Trinity v2.1.1, TransDecoder v2.0.1, and CD-HIT v4.6 software. The differentially expressed gene analysis revealed that the expression of genes related to the flavonoid pathway was higher in the flowers, whereas the phenylpropanoid pathway was more highly expressed in the roots. In this study, we established a global transcriptome dataset for C. japonicum. The data shall not only be useful to focus more deeply on the genes related to product medicinal metabolite including flavolignan but also to study the functional genomics for genetic engineering of C. japonicum.

Lilium dauricum Ker Gawler is a wild lily species that belongs to section Sinomartagon and is one of the ancestors of the Asiatic hybrid lilies. Unique traits such as disease resistance and early flowering make L. dauricum a desirable resource for interspecific hybridization. However, in Korea, the natural resources of L. dauricum are being exhausted by excessive exploitation and require urgent conservation. In this study, the complete chloroplast genome of L. dauricum was generated using Illumina paired-end sequencing technology, and its structure was compared with that of other Lilium species. The chloroplast genome was 152,063 bp long, with a typical quadripartite structure including a pair of inverted repeat regions (26,492 bp) separated by a large (81,485 bp) and small (17,584 bp) single-copy (SC) region. The genome encodes 131 different genes, including 85 protein-coding genes, 38 tRNA genes, and 8 rRNA genes. A phylogenetic analysis strongly supported the relationship of L. dauricum with other members of Sinomartagon and Martagon lilies.

Using next-generation sequencing, we conducted a genome-wide scan of selective sweeps associated with selection toward genetic improvement in Thoroughbreds. We investigated potential phenotypic consequence of putative candidate loci by candidate gene association mapping for the finishing time in 240 Thoroughbred horses. We found a significant association with the trait for Ral GApase alpha 2 (RALGAP2) that regulates a variety of cellular processes of signal trafficking. Neighboring genes around RALGAP2 included insulinoma-associated 1 (INSM1), pallid (PLDN), and Ras and Rab interactor 2 (RIN2) genes have similar roles in signal trafficking, suggesting that a co-evolving gene cluster located on the chromosome 22 is under strong artificial selection in racehorses.

Molecular characterization technology in genetically modified organisms, in addition to how transgenic biotechnologies are developed now require full transparency to assess the risk to living modified and non-modified organisms. Next generation sequencing (NGS) methodology is suggested as an effective means in genome characterization and detection of transgenic insertion locations. In the present study, we applied NGS to insert transgenic loci, specifically the epidermal growth factor (EGF) in genetically modified rice cells. A total of 29.3 Gb (~72× coverage) was sequenced with a 2 × 150 bp paired end method by Illumina HiSeq2500, which was consecutively mapped to the rice genome and T-vector sequence. The compatible pairs of reads were successfully mapped to 10 loci on the rice chromosome and vector sequences were validated to the insertion location by polymerase chain reaction (PCR) amplification. The EGF transgenic site was confirmed only on chromosome 4 by PCR. Results of this study demonstrated the success of NGS data to characterize the rice genome. Bioinformatics analyses must be developed in association with NGS data to identify highly accurate transgenic sites.

Thoroughbred, a relatively recent horse breed, is best known for its use in horse racing. Although myostatin (MSTN) variants have been reported to be highly associated with horse racing performance, the trait is more likely to be polygenic in nature. The purpose of this study was to identify genetic variants strongly associated with racing performance by using estimated breeding value (EBV) for race time as a phenotype. We conducted a two-stage genome-wide association study to search for genetic variants associated with the EBV. In the first stage of genome-wide association study, a relatively large number of markers (~54,000 single-nucleotide polymorphisms, SNPs) were evaluated in a small number of samples (240 horses). In the second stage, a relatively small number of markers identified to have large effects (170 SNPs) were evaluated in a much larger number of samples (1,156 horses). We also validated the SNPs related to MSTN known to have large effects on racing performance and found significant associations in the stage two analysis, but not in stage one. We identified 28 significant SNPs related to 17 genes. Among these, six genes have a function related to myogenesis and five genes are involved in muscle maintenance. To our knowledge, these genes are newly reported for the genetic association with racing performance of Thoroughbreds. It complements a recent horse genome-wide association studies of racing performance that identified other SNPs and genes as the most significant variants. These results will help to expand our knowledge of the polygenic nature of racing performance in Thoroughbreds.

Chloroplast (cp) genomes of Lilium amabile, L. callosum, L. lancifolium, and L. philadelphicum were fully sequenced. Using these four novel cp genome sequences and five other previously sequenced cp genomes, features of the cp genomes were characterized in detail among species in the genus Lilium and other related genera in the order Liliales. The lengths and nucleotide composition showed little variation. No structural variation was found among the cp genomes in Liliales. Gene contents were conserved among four newly sequenced cp genome in Lilium species, the only differences being in two pseudogenes. We identified 112 genes in 13 functional categories, 18 of which carried introns that were conserved among the species in Liliales. There were 16-21 SSR loci (>12 bp, >3 repeats) in the cp genomes in Lilium and the genomic locations of these loci were highly variable among the species. Average mutations were 15 SNPs per 1kb and 5 indels per 1kb, respectively, in the cp genomes of the newly sequenced four Lilium species. Phylogenetic classifications revealed some discrepancies between trees based on the cp genomes and previous classifications based on the morphology and geographic distributions.

Recently, the technologies of DNA sequence variation and gene expression profiling have been used widely as approaches in the expertise of genome biology and genetics. The application to genome study has been particularly developed with the introduction of the next-generation DNA sequencer (NGS) Roche/454 and Illumina/Solexa systems, along with bioinformation analysis technologies of whole-genome de novo assembly, expression profiling, DNA variation discovery, and genotyping. Both massive whole-genome shotgun paired-end sequencing and mate paired-end sequencing data are important steps for constructing de novo assembly of novel genome sequencing data. It is necessary to have DNA sequence information from a multiplatform NGS with at least 2× and 30× depth sequence of genome coverage using Roche/454 and Illumina/Solexa, respectively, for effective an way of de novo assembly. Massive short-length reading data from the Illumina/Solexa system is enough to discover DNA variation, resulting in reducing the cost of DNA sequencing. Whole-genome expression profile data are useful to approach genome system biology with quantification of expressed RNAs from a whole-genome transcriptome, depending on the tissue samples. The hybrid mRNA sequences from Rohce/454 and Illumina/Solexa are more powerful to find novel genes through de novo assembly in any whole-genome sequenced species. The 20× and 50× coverage of the estimated transcriptome sequences using Roche/454 and Illumina/Solexa, respectively, is effective to create novel expressed reference sequences. However, only an average 30× coverage of a transcriptome with short read sequences of Illumina/Solexa is enough to check expression quantification, compared to the reference expressed sequence tag sequence.

Colletotrichum acutatum is a destructive fungal pathogen which causes anthracnose in a wide range of crops. Here we report the whole genome sequence and annotation of C. acutatum strain KC05, isolated from an infected pepper in Kangwon, South Korea. Genomic DNA from the KC05 strain was used for the whole genome sequencing using a PacBio sequencer and the MiSeq system. The KC05 genome was determined to be 52,190,760 bp in size with a G + C content of 51.73% in 27 scaffolds and to contain 13,559 genes with an average length of 1516 bp. Gene prediction and annotation were performed by incorporating RNA-Seq data. The genome sequence of the KC05 was deposited at DDBJ/ENA/GenBank under the accession number LUXP00000000.

MicroRNAs (miRNAs) are a class of non-coding RNAs approximately 21 nt in length which play important roles in regulating gene expression in plants. Although many miRNA studies have focused on a few model plants, miRNAs and their target genes remain largely unknown in hot pepper (Capsicum annuum), one of the most important crops cultivated worldwide. Here, we employed high-throughput sequencing technology to identify miRNAs in pepper extensively from 10 different libraries, including leaf, stem, root, flower, and six developmental stage fruits. Based on a bioinformatics pipeline, we successfully identified 29 and 35 families of conserved and novel miRNAs, respectively. Northern blot analysis was used to validate further the expression of representative miRNAs and to analyze their tissue-specific or developmental stage-specific expression patterns. Moreover, we computationally predicted miRNA targets, many of which were experimentally confirmed using 5' rapid amplification of cDNA ends analysis. One of the validated novel targets of miR-396 was a domain rearranged methyltransferase, the major de novo methylation enzyme, involved in RNA-directed DNA methylation in plants. This work provides the first reliable draft of the pepper miRNA transcriptome. It offers an expanded picture of pepper miRNAs in relation to other plants, providing a basis for understanding the functional roles of miRNAs in pepper.