Wen Luan Wu

An overview of the Phalaenopsis orchid genome through BAC end sequence analysis

BackgroundPhalaenopsis orchids are popular floral crops, and development of new cultivars is economically important to floricultural industries worldwide. Analysis of orchid genes could facilitate orchid improvement. Bacterial artificial chromosome (BAC) end sequences (BESs) can provide the first glimpses into the sequence composition of a novel genome and can yield molecular markers for use in genetic mapping and breeding.ResultsWe used two BAC libraries (constructed using the BamHI and HindIII restriction enzymes) of Phalaenopsis equestris to generate pair-end sequences from 2,920 BAC clones (71.4% and 28.6% from the BamHI and HindIII libraries, respectively), at a success rate of 95.7%. A total of 5,535 BESs were generated, representing 4.5 Mb, or about 0.3% of the Phalaenopsis genome. The trimmed sequences ranged from 123 to 1,397 base pairs (bp) in size, with an average edited read length of 821 bp. When these BESs were subjected to sequence homology searches, it was found that 641 (11.6%) were predicted to represent protein-encoding regions, whereas 1,272 (23.0%) contained repetitive DNA. Most of the repetitive DNA sequences were gypsy- and copia-like retrotransposons (41.9% and 12.8%, respectively), whereas only 10.8% were DNA transposons. Further, 950 potential simple sequence repeats (SSRs) were discovered. Dinucleotides were the most abundant repeat motifs; AT/TA dimer repeats were the most frequent SSRs, representing 253 (26.6%) of all identified SSRs. Microsynteny analysis revealed that more BESs mapped to the whole-genome sequences of poplar than to those of grape or Arabidopsis, and even fewer mapped to the rice genome. This work will facilitate analysis of the Phalaenopsis genome, and will help clarify similarities and differences in genome composition between orchids and other plant species.ConclusionUsing BES analysis, we obtained an overview of the Phalaenopsis genome in terms of gene abundance, the presence of repetitive DNA and SSR markers, and the extent of microsynteny with other plant species. This work provides a basis for future physical mapping of the Phalaenopsis genome and advances our knowledge thereof.

The comparative chloroplast genomic analysis of photosynthetic orchids and developing DNA markers to distinguish Phalaenopsis orchids

The chloroplast genome of Phalaenopsis equestris was determined and compared to those of Phalaenopsis aphrodite and Oncidium Gower Ramsey in Orchidaceae. The chloroplast genome of P. equestris is 148,959 bp, and a pair of inverted repeats (25,846 bp) separates the genome into large single-copy (85,967 bp) and small single-copy (11,300 bp) regions. The genome encodes 109 genes, including 4 rRNA, 30 tRNA and 75 protein-coding genes, but loses four ndh genes (ndhA, E, F and H) and seven other ndh genes are pseudogenes. The rate of inter-species variation between the two moth orchids was 0.74% (1107 sites) for single nucleotide substitution and 0.24% for insertions (161 sites; 1388 bp) and deletions (189 sites; 1393 bp). The IR regions have a lower rate of nucleotide substitution (3.5-5.8-fold) and indels (4.3-7.1-fold) than single-copy regions. The intergenic spacers are the most divergent, and based on the length variation of the three intergenic spacers, 11 native Phalaenopsis orchids could be successfully distinguished. The coding genes, IR junction and RNA editing sites are relatively more conserved between the two moth orchids than between those of Phalaenopsis and Oncidium spp.

Genetic variation and population differentiation of Michelia formosana (Magnoliaceae) based on cpDNA variation and RAPD fingerprints: relevance to post-Pleistocene recolonization

Abstractu2002We used sequence variation of the atpB-rbcL intergenic spacer of cpDNA and nested clade analysis to assess the phylogeographic pattern of Michelia formosana, a species restricted to Taiwan and the Ryukyus. In total, 31 haplotypes were identified and clustered into four major chlorotypes. Genetic composition of nearly all populations was heterogeneous and paraphyletic phylogenetically. Although the apportionment of cpDNA variation hardly revealed a geographic pattern due to the coancestry of dominant sequences, some chlorotypes were restrictedly distributed. According to the patterns of clade dispersion and displacement, a reconstructed minimum spanning network revealed that historical events of past fragmentation and range expansion, associated with glaciation, may have shaped the phylogeographic patterns of M. formosana. Four possible refugia were identified: the Iriomote and Ishigaki Islands (the southern Ryukyus), Wulai (northern Taiwan), and Nanjen (southern Taiwan), on the basis of the interior positions of their haplotypes in the network and the high level of nucleotide diversity. Given insufficient time for coalescence at the cpDNA locus since the late Pleistocene recolonization, lineage sorting led to low levels of genetic differentiation among populations. In contrast, hierarchical examination of the random amplified polymorphic DNA (RAPD) data scored from six populations across three geographical regions, using an analysis of molecular variance (AMOVA), indicated high genetic differentiation both among populations (ΦST = 0.471) and among regions (ΦCT = 0.368). An unweighted pair group method with arithmetic mean (UPGMA) tree of the RAPD fingerprints revealed that populations of two offshore islands of eastern Taiwan (M. formosana var. kotoensis) were clustered with geographically remote populations of the Ryukyus instead of those in southern Taiwan, suggesting some historical division due to geographic barriers of the central mountain range. In contrast to the paraphyly of the nearly neutral cpDNA alleles, differentiated RAPDs may have experienced diversifying selection.

Overexpression of TMAC2, a novel negative regulator of abscisic acid and salinity responses, has pleiotropic effects in Arabidopsis thaliana

Phytohormone abscisic acid (ABA) regulates many aspects of plant development and growth. To explore the molecular mechanism of ABA, we identified the novel ABA-regulated genes in Arabidopsis thaliana by searching for genes possessing two or more ABREs (ABA-responsive elements). One of these genes, two or more ABREs-containing gene 2 (TMAC2) is highly induced by ABA and NaCl. Database searches revealed that TMAC2 encodes a protein with no domains of known function. Expression of TMAC2-GFP fusion protein in Arabidopsis mesophyll protoplasts indicated that TMAC2 is targeted to the nucleus. Although the gene has a basal level of expression in various Arabidopsis organs/tissues except for adult leaves, a high expression level was detected in roots. Constitutive overexpression of TMAC2 in plants resulted in the insensitivity to ABA and NaCl, suggesting that TMAC2 plays a negative role in ABA and salt stress responses. Furthermore, TMAC2-overexpressing plants exhibited the short roots, late flowering and starch-excess phenotypes. RT-PCR analysis showed that decreased expression of two floral- and one starch degradation-related genes, SOC1/AGL20 and SEP3/AGL9, and SEX1, respectively, may lead to altered phenotypes of TMAC2-overexpressing plants. Taken together, our data reveal that TMAC2 acts in the nucleus and is an important negative regulator of ABA and salt stress responses, and could play a critical role in controlling root elongation, floral initiation and starch degradation.

Deciphering mycorrhizal fungi in cultivated Phalaenopsis microbiome with next-generation sequencing of multiple barcodes

Identifying the species composition of a microbial ecosystem is often hampered by difficulties in culturing the organisms and in the low sequencing depth of traditional DNA barcoding. Metagenomic analysis, a huge-scale nucleotide-sequence-based tool, can overcome such difficulties. In this study, Sanger sequencing of 500 nrITS clones uncovered 29 taxa of 19 fungal genera, whereas metagenomics with next-generation sequencing identified 512 operational taxonomic units (OTUs) for ITS1/2 and 364 for ITS3/4. Nevertheless, high throughput sequencing of PCR amplicons of ITS1/2, ITS3/4, nrLSU-LR, nrLSU-U, mtLSU, and mtATP6, all with at least 1,300× coverage and about 21 million reads in total, yielded a very diverse fungal composition. The fact that 74xa0% of the OTUs were exclusively uncovered with single barcodes indicated that each marker provided its own insights into the fungal flora. To deal with the high heterogeneity in the data and to integrate the information on species composition across barcodes, a rank-scoring strategy was developed. Accordingly, 205 genera among 64 orders of fungi were identified in healthy Phalaenopsis roots. Of the barcodes utilized, ITS1/2, ITS3/4, and nrLSU-U were the most competent in uncovering the fungal diversity. These barcodes, though detecting different compositions likely due to primer preference, provided complementary and comprehensive power in deciphering the microbial diversity, especially in revealing rare species.

Characterization and adaptive evolution of α-tubulin genes in the Miscanthus sinensis complex (Poaceae)

To investigate the organization of and mode of selection in the α-tubulin genes, full-length α-tubulin genes were cloned from four intraspecific taxa of Miscanthus sinensis and its close relatives M. floridulus and M. condensatus using standard polymerase chain reaction (PCR) and rapid amplification of genomic ends (RAGE)-PCR strategies. Genealogical analysis of angiosperms recovered a monophyletic group of Miscanthus α-tubulin genes, which is homologous to the tua5 locus of maize. Two clusters of nearly equal frequency revealed paraphyly within each Miscanthus taxon. Between-cluster recombination was frequent. Additional evidence for co-occurrence of two haplotypes within individuals and a large-scale crossover all suggested a likely allelic relationship between the Miscanthus clusters. Given a long between-species divergence time in Miscanthus, wide occurrence of the trans-species polymorphisms in α-tubulin genes and the approximately equal frequency of each allelic type make it extremely unlikely that α-tubulin diversity has been maintained under neutrality. Balancing selection may have contributed to such an apportioning of genetic variability as well as to high levels of genetic variation in α-tubulin and higher substitution rates at synonymous sites of exons than at intron bases of M. sinensis. In addition, certain effects of demographic oscillation may have distorted the scenario of a functional locus operating under balancing selection.

Histone Acetylation Accompanied with Promoter Sequences Displaying Differential Expression Profiles of B-Class MADS-Box Genes for Phalaenopsis Floral Morphogenesis

Five B-class MADS-box genes, including four APETALA3 (AP3)-like PeMADS2∼5 and one PISTILLATA (PI)-like PeMADS6, specify the spectacular flower morphology in orchids. The PI-like PeMADS6 ubiquitously expresses in all floral organs. The four AP3-like genes, resulted from two duplication events, express ubiquitously at floral primordia and early floral organ stages, but show distinct expression profiles at late floral organ primordia and floral bud stages. Here, we isolated the upstream sequences of PeMADS2∼6 and studied the regulatory mechanism for their distinct gene expression. Phylogenetic footprinting analysis of the 1.3-kb upstream sequences of AP3-like PeMADS2∼5 showed that their promoter regions have sufficiently diverged and contributed to their subfunctionalization. The amplified promoter sequences of PeMADS2∼6 could drive beta-glucuronidase (GUS) gene expression in all floral organs, similar to their expression at the floral primordia stage. The promoter sequence of PeMADS4, exclusively expressed in lip and column, showed a 1.6∼3-fold higher expression in lip/column than in sepal/petal. Furthermore, we noted a 4.9-fold increase in histone acetylation (H3K9K14ac) in the translation start region of PeMADS4 in lip as compared in petal. All these results suggest that the regulation via the upstream sequences and increased H3K9K14ac level may act synergistically to display distinct expression profiles of the AP3-like genes at late floral organ primordia stage for Phalaenopsis floral morphogenesis.

The mitochondrial DNA markers for distinguishing Phalaenopsis species and revealing maternal phylogeny

Moth orchids (Phalaenopsis) are among the top-traded blooming potted plants in the world. To explore mitochondrial DNA (mtDNA) markers for species identification, we located simple sequence repeats in the mtDNA of Phalaenopsis aphrodite subsp. formosana and then pre-screened them for polymorphic markers by their comparison with corresponding mtDNA regions of P. equestris. The combination of 13 selected markers located in intergenic spacers could unambiguously distinguish 15 endemic moth orchids. Five most variable markers with polymorphic information content (PIC) ≥ 0.7 could be combined to classify 18 of 19 endemic moth orchids including parental strains most commonly used in breeding programs. The sequences of four selected mtDNA regions were highly variable, and one region (MT2) could be used to completely distinguish 19 endemic moth orchids. Though mitochondrial introns were highly conserved among moth orchids, evolutionary hotspots, such as variable simple sequence repeats and minisatellite repeats, were identified as useful markers. Furthermore, a marker technology was applied to reveal the maternal inheritance mode of mtDNA in the moth orchids. Moreover, phylogenetic analysis indicates that the mtDNA was nonmonophyletic below the Phalaenopsis genus. In summary, we have revealed a set of mtDNA markers that could be used for identification and phylogenetic study of Phalaenopsis orchids.

Genbank accession #: AJ133707

GenBank accession #: AJ133707 Authors: Wu WL and Hsiao IL Address: Department of Biology, Cheng-Kung University, Tainan, Taiwan 700 Source of sequence: Sequence was determined from a PCR product, which was ligated to a T-vector, from genomic DNA ofEchinochloa crus-galli , a wild species of grass. Trivial name: rDNA ITS of Echinochloa crus-galli. Description: ITS region of ribosomal DNA lying between 16S and 25S rRNA genes was cloned from a wild species,Echinochloa crus-galli , with 207 bp in ITS1, 163 bp in 5.8S rRNA gene and 218 bp in ITS2. Nucleotide sequence of the ITS region in Echinochloa crus-galli shares 89.5%, 89.3% and 87.2% identity with Entolasia stricta, Panicum bisulcatum, andImperata cylindrica, respectively. Acknowledgements:We wish to thank Dr C.S. Kuoh for help with the collection of E. crus-galliplants.

Genbank accession #: 133708

GenBank accession #: 133708 Authors: Wu WL and Hsiao LL Address: Department of Biology, Cheng-Kung University, Tainan, Taiwan 700 Source of sequence: Sequence was determined from a PCR product, which was ligated to a T-vector, from genomic DNA ofEchinochloa colonum, a wild species of grass. Trivial name: rDNA ITS of Echinochloa colonum. Description: ITS region of ribosomal DNA lying between 16S and 25S rRNA genes was cloned from a wild species,Echinochloa colonum, with 206 bp in ITS1, 163 bp in 5.8S rRNA gene and 216 bp in ITS2. Nucleotide sequence of ITS region in Echinochloa colonumshares 88.4%, 87.2% and 86.1% identity with Entolasia stricta, Imperata cylindrica, andSorghum versicolor , respectively. Acknowledgements:We wish to thank Dr C.S. Kuoh for help with the collection of E. colonumplants.