Hee-Ju Yu

Genome-wide comparative analysis of the Brassica rapa gene space reveals genome shrinkage and differential loss of duplicated genes after whole genome triplication

BackgroundBrassica rapa is one of the most economically important vegetable crops worldwide. Owing to its agronomic importance and phylogenetic position, B. rapa provides a crucial reference to understand polyploidy-related crop genome evolution. The high degree of sequence identity and remarkably conserved genome structure between Arabidopsis and Brassica genomes enables comparative tiling sequencing using Arabidopsis sequences as references to select the counterpart regions in B. rapa, which is a strong challenge of structural and comparative crop genomics.ResultsWe assembled 65.8 megabase-pairs of non-redundant euchromatic sequence of B. rapa and compared this sequence to the Arabidopsis genome to investigate chromosomal relationships, macrosynteny blocks, and microsynteny within blocks. The triplicated B. rapa genome contains only approximately twice the number of genes as in Arabidopsis because of genome shrinkage. Genome comparisons suggest that B. rapa has a distinct organization of ancestral genome blocks as a result of recent whole genome triplication followed by a unique diploidization process. A lack of the most recent whole genome duplication (3R) event in the B. rapa genome, atypical of other Brassica genomes, may account for the emergence of B. rapa from the Brassica progenitor around 8 million years ago.ConclusionsThis work demonstrates the potential of using comparative tiling sequencing for genome analysis of crop species. Based on a comparative analysis of the B. rapa sequences and the Arabidopsis genome, it appears that polyploidy and chromosomal diploidization are ongoing processes that collectively stabilize the B. rapa genome and facilitate its evolution.

Auxin response factor gene family in Brassica rapa: genomic organization, divergence, expression, and evolution

Completion of the sequencing of the Brassica rapa genome enabled us to undertake a genome-wide identification and functional study of the gene families related to the morphological diversity and agronomic traits of Brassica crops. In this study, we identified the auxin response factor (ARF) gene family, which is one of the key regulators of auxin-mediated plant growth and development in the B. rapa genome. A total of 31 ARF genes were identified in the genome. Phylogenetic and evolutionary analyses suggest that ARF genes fell into four major classes and were amplified in the B. rapa genome as a result of a recent whole genome triplication after speciation from Arabidopsis thaliana. Despite its recent hexaploid ancestry, B. rapa includes a relatively small number of ARF genes compared with the 23 members in A. thaliana, presumably due to a paralog reduction related to repetitive sequence insertion into promoter and non-coding transcribed region of the genes. Comparative genomic and mRNA sequencing analyses demonstrated that 27 of the 31 BrARF genes were transcriptionally active, and their expression was affected by either auxin treatment or floral development stage, although 4 genes were inactive, suggesting that the generation and pseudogenization of ARF members are likely to be an ongoing process. This study will provide a fundamental basis for the modification and evolution of the gene family after a polyploidy event, as well as a functional study of ARF genes in a polyploidy crop species.

Petunia actin-depolymerizing factor is mainly accumulated in vascular tissue and its gene expression is enhanced by the first intron.

Actin-depolymerizing factor (ADF) is one of the actin cytoskeleton-modulating proteins. We have characterized the accumulation pattern of petunia ADF proteins. PhADF proteins are accumulated in every petunia organ and their accumulation is differentially regulated by developmental signals. Their cellular localization is vascular tissue-preferential in vegetative organs, whereas somewhat different in reproductive organs. In reproductive organs, PhADFs are present in outer integument, endocarp of ovary wall, transmitting tissue of style, and epidermis and endothecium of young anther. From a petunia genomic library, we have isolated a genomic clone encoding PhADF1. Comparison to complementary DNA sequence revealed that the coding region of PhADF1 gene consists of three exons and two introns. Analysis of chimeric gene expression using beta-glucuronidase as a reporter gene in transgenic Arabidopsis revealed that PhADF1 was strongly expressed in every vegetative tissue except petal. In addition, expression of the gene was highly enhanced by its first intron. These results suggest that PhADF1 gene of petunia is mainly expressed in vascular tissues and its expression is regulated by intron-mediated enhancement mechanism.

Identification and profiling of novel microRNAs in the Brassica rapa genome based on small RNA deep sequencing

BackgroundMicroRNAs (miRNAs) are one of the functional non-coding small RNAs involved in the epigenetic control of the plant genome. Although plants contain both evolutionary conserved miRNAs and species-specific miRNAs within their genomes, computational methods often only identify evolutionary conserved miRNAs. The recent sequencing of the Brassica rapa genome enables us to identify miRNAs and their putative target genes. In this study, we sought to provide a more comprehensive prediction of B. rapa miRNAs based on high throughput small RNA deep sequencing.ResultsWe sequenced small RNAs from five types of tissue: seedlings, roots, petioles, leaves, and flowers. By analyzing 2.75 million unique reads that mapped to the B. rapa genome, we identified 216 novel and 196 conserved miRNAs that were predicted to target approximately 20% of the genome’s protein coding genes. Quantitative analysis of miRNAs from the five types of tissue revealed that novel miRNAs were expressed in diverse tissues but their expression levels were lower than those of the conserved miRNAs. Comparative analysis of the miRNAs between the B. rapa and Arabidopsis thaliana genomes demonstrated that redundant copies of conserved miRNAs in the B. rapa genome may have been deleted after whole genome triplication. Novel miRNA members seemed to have spontaneously arisen from the B. rapa and A. thaliana genomes, suggesting the species-specific expansion of miRNAs. We have made this data publicly available in a miRNA database of B. rapa called BraMRs. The database allows the user to retrieve miRNA sequences, their expression profiles, and a description of their target genes from the five tissue types investigated here.ConclusionsThis is the first report to identify novel miRNAs from Brassica crops using genome-wide high throughput techniques. The combination of computational methods and small RNA deep sequencing provides robust predictions of miRNAs in the genome. The finding of numerous novel miRNAs, many with few target genes and low expression levels, suggests the rapid evolution of miRNA genes. The development of a miRNA database, BraMRs, enables us to integrate miRNA identification, target prediction, and functional annotation of target genes. BraMRs will represent a valuable public resource with which to study the epigenetic control of B. rapa and other closely related Brassica species. The database is available at the following link: http://bramrs.rna.kr [1].

De novo assembly and characterization of the complete chloroplast genome of radish (Raphanus sativus L.).

Radish (Raphanus sativus L.) is an edible root vegetable crop that is cultivated worldwide and whose genome has been sequenced. Here we report the complete nucleotide sequence of the radish cultivar WK10039 chloroplast (cp) genome, along with a de novo assembly strategy using whole genome shotgun sequence reads obtained by next generation sequencing. The radish cp genome is 153,368 bp in length and has a typical quadripartite structure, composed of a pair of inverted repeat regions (26,217 bp each), a large single copy region (83,170 bp), and a small single copy region (17,764 bp). The radish cp genome contains 87 predicted protein-coding genes, 37 tRNA genes, and 8 rRNA genes. Sequence analysis revealed the presence of 91 simple sequence repeats (SSRs) in the radish cp genome. Phylogenetic analysis of 62 protein-coding gene sequences from the 17 cp genomes of the Brassicaceae family suggested that the radish cp genome is most closely related to the cp genomes of Brassica rapa and Brassicanapus. Comparisons with the B. rapa and B. napus cp genomes revealed highly divergent intergenic sequences and introns that can potentially be developed as diagnostic cp markers. Synonymous and nonsynonymous substitutions of cp genes suggested that nucleotide substitutions have occurred at similar rates in most genes. The complete sequence of the radish cp genome would serve as a valuable resource for the development of new molecular markers and the study of the phylogenetic relationships of Raphanus species in the Brassicaceae family.

The complete chloroplast genome of Phalaenopsis "Tiny Star".

Abstract We determined the complete chloroplast DNA sequence of Phalaenopsis “Tiny Star” based on Illumina sequencing. The total length of the chloroplast genome is 148,918 bp long with GC content of 36.7%. It contains 70 protein-coding genes, 30 tRNA genes, and 4 rRNA genes. Comparative analysis with the reported orchid chloroplast sequences identified unique InDel variations in the “Tiny Star” chloroplast genome that have potential as genetic markers to investigate the maternal lineage of Phalaenopsis and Doritaenopsis cultivars.

The complete chloroplast genome of Aconitum chiisanense Nakai (Ranunculaceae)

Abstract We determined the complete chloroplast DNA sequence of Aconitum chiisanense Nakai, a rare Aconitum species endemic to Korea. The chloroplast genome is 155 934 bp in length and contains 4 rRNA, 30 tRNA, and 78 protein-coding genes. Phylogenetic analysis revealed that the chloroplast genome of A. chiisanense is closely related to that of A. barbatum var. puberulum. Sequence comparison with other Ranunculaceae chloroplasts identified a unique deletion in the rps16 gene of A. chiisanense chloroplast DNA that can serve as a molecular marker for species identification.

The complete mitochondrial genome of cultivated radish WK10039 (Raphanus sativus L.)

Abstract We determined the complete nucleotide sequence of the mitochondrial genome of radish cultivar WK10039 (Raphanus sativus L.). The total length of the mtDNA sequence is 244,054 bp, with GC content of 45.3%. The radish mtDNA contains 82 protein-coding genes, 17 tRNA genes, and 3 rRNA genes. Among the protein-coding genes, 34 encode proteins with known functions. There are two 5529 bp repeats in the radish mitochondrial genome that may contribute to DNA recombination resulting in at least three different forms of mtDNA in radish.

Identification of candidate domestication regions in the radish genome based on high-depth resequencing analysis of 17 genotypes

Key messageThis study provides high-quality variation data of diverse radish genotypes. Genome-wide SNP comparison along with RNA-seq analysis identified candidate genes related to domestication that have potential as trait-related markers for genetics and breeding of radish.AbstractRadish (Raphanus sativus L.) is an annual root vegetable crop that also encompasses diverse wild species. Radish has a long history of domestication, but the origins and selective sweep of cultivated radishes remain controversial. Here, we present comprehensive whole-genome resequencing analysis of radish to explore genomic variation between the radish genotypes and to identify genetic bottlenecks due to domestication in Asian cultivars. High-depth resequencing and multi-sample genotyping analysis of ten cultivated and seven wild accessions obtained 4.0 million high-quality homozygous single-nucleotide polymorphisms (SNPs)/insertions or deletions. Variation analysis revealed that Asian cultivated radish types are closely related to wild Asian accessions, but are distinct from European/American cultivated radishes, supporting the notion that Asian cultivars were domesticated from wild Asian genotypes. SNP comparison between Asian genotypes identified 153 candidate domestication regions (CDRs) containing 512 genes. Network analysis of the genes in CDRs functioning in plant signaling pathways and biochemical processes identified group of genes related to root architecture, cell wall, sugar metabolism, and glucosinolate biosynthesis. Expression profiling of the genes during root development suggested that domestication-related selective advantages included a main taproot with few branched lateral roots, reduced cell wall rigidity and favorable taste. Overall, this study provides evolutionary insights into domestication-related genetic selection in radish as well as identification of gene candidates with the potential to act as trait-related markers for background selection of elite lines in molecular breeding.

Phenotypic analysis of parents and their reciprocal F1 hybrids in Phalaenopsis

In this study, we studied the phenotypic and breeding efficiency of Phalaenopsis reciprocal hybrids and their parents. For reciprocal hybridization, Phalaenopsis ‘KS Little Gem’ and ‘1747’ from Taiwan were used as parents. After crossing ‘KS Little Gem’ × ‘1747’ in 2010, 34 individuals from ‘KS 1059’ hybrids were developed. The crossing between ‘1747’ × ‘KS Little Gem’ resulted in the development of 63 individuals from ‘KS 1076’ hybrids. Detailed morphological characteristics of newly developed hybrids; ‘KS 1059’, ‘KS 1076’ and their parents were collected according to guidelines determined by the Korea Seed and Variety Service (KSVS). The leaf characteristics were similar in reciprocal hybrids and their parents, except in number of leaves, which as were maximal in ‘KS 1059’. In terms of flower characteristics including number of flowers, inflorescence length, pedicel length and pedicel diameter, features of the parents were combined in the hybrids. By contrast in case of flower length, both hybrids showed similarity to ‘KS Little Gem’, and a similar pattern was observed in flower width. The petal length in hybrids and parents were similar, but the petal width of hybrids was more closer to that of ‘KS Little Gem’. The length of whiskers was long in both parents, but very short whiskers were recorded in hybrids. The length and width of the apical lobe were similar in parents and hybrids, but hybrids were closer to ‘1747’. In the case of duration of of flowering, the maximum value of 123 days was recorded in ‘KS Little Gem’ similar to ‘KS 1059’, whereas ‘1747’ and ‘KS 1076’ had almost identical values. It can be concluded that to generate Phalaenopsis with a longer shelf life, ‘KS Little Gem’ might be selected as the female parent during hybridization.