Plant Genomics

Microscopy and mass spectrometry (MS) are complementary techniques: The former provides spatiotemporal information in living cells, but only for a handful of recombinant proteins at a time, whereas the latter can detect thousands of endogenous proteins simultaneously, but only in lysed samples. Here, we introduce technology that combines these strengths by offering spatially and temporally resolved proteomic maps of endogenous proteins within living cells. Our method relies on a genetically targetable peroxidase enzyme that biotinylates nearby proteins, which are subsequently purified and identified by MS. We used this approach to identify 495 proteins within the human mitochondrial matrix, including 31 not previously linked to mitochondria. The labeling was exceptionally specific and distinguished between inner membrane proteins facing the matrix versus the intermembrane space (IMS). Several proteins previously thought to reside in the IMS or outer membrane, including protoporphyrinogen oxidase, were reassigned to the matrix by our proteomic data and confirmed by electron microscopy. The specificity of peroxidase-mediated proteomic mapping in live cells, combined with its ease of use, offers biologists a powerful tool for understanding the molecular composition of living cells.

It is likely that many small ORFs (sORFs; 30–100 amino acids) are missed when genomes are annotated. To overcome this limitation, we identified ∼8,000 sORFs with high coding potential in intergenic regions of the Arabidopsis thaliana genome. However, the question remains as to whether these coding sORFs play functional roles. Using a designed array, we generated an expression atlas for 16 organs and 17 environmental conditions among 7,901 identified coding sORFs. A total of 2,099 coding sORFs were highly expressed under at least one experimental condition, and 571 were significantly conserved in other land plants. A total of 473 coding sORFs were overexpressed; ∼10% (49/473) induced visible phenotypic effects, a proportion that is approximately seven times higher than that of randomly chosen known genes. These results indicate that many coding sORFs hidden in plant genomes are associated with morphogenesis. We believe that the expression atlas will contribute to further study of the roles of sORFs in plants.

Sacred lotus is a basal eudicot with agricultural, medicinal, cultural and religious importance. It was domesticated in Asia about 7,000 years ago, and cultivated for its rhizomes and seeds as a food crop. It is particularly noted for its 1,300-year seed longevity and exceptional water repellency, known as the lotus effect. The latter property is due to the nanoscopic closely-packed protuberances on its self-cleaning leaf surface, which have been adapted for the manufacture of a self-cleaning industrial paint, Lotusan.

The genome of the China Antique variety of the sacred lotus was sequenced with Illumina and 454 technologies, at respective depths of 101x and 5.2x. The final assembly has a contig N50 of 38.8 kbp and a scaffold N50 of 3.4 Mbp, and covers 86.5% of the estimated 929 Mbp total genome size. The genome notably lacks the paleo-triplication observed in other eudicots, but reveals a lineage-specific duplication. The genome has evidence of slow evolution, with a 30% slower nucleotide mutation rate than observed in grape. Comparisons of the available sequenced genomes suggest a minimum gene set for vascular plants of 4,223 genes. Strikingly, the sacred lotus has sixteen COG2132 multi-copper oxidase family proteins with root specific expression; these are involved in root meristem phosphate starvation, reflecting adaptation to limited nutrient availability in an aquatic environment.

The slow nucleotide substitution rate makes the sacred lotus a better resource than the current standard, grape, for reconstructing the pan-eudicot genome, and should therefore accelerate comparative analysis between eudicots and monocots.

Primula species are important early spring garden plants with a centre of diversity and speciation in the East Himalaya-Hengduan Mountains in Western China. Studies on population genetics, speciation and phylogeny of Primula have been impeded by a lack of genomic resources. In the present study, we sequenced the transcriptomes of two closely related primrose species, Primula poissonii and Primula wilsonii, using short reads on the Illumina Genome Analyzer platform.

We obtained 55,284 and 55,011 contigs with N50 values of 938 and 1,085 for P. poissonii and P. wilsonii, respectively, and 6,654 pairs of putative orthologs were identified between the two species. Estimations of non-synonymous/synonymous substitution rate ratios for these orthologs indicated that 877 of the pairs may be under positive selection (Ka/Ks > 0.5), and functional enrichment analysis revealed that significant proportions of the orthologs were in the categories DNA repair, stress resistance, which may provide some hints as to how the two closely related Primula species adapted differentially to extreme environments, such as habitats characterized by aridity, high altitude and high levels of ionizing radiation. It was possible for the first time to estimate the divergence time between the radiated species pair, P. poissonii and P. wilsonii; this was found to be approximately 0.90 +/- 0.57 Mya, which falls between the Donau and Gunz glaciation in the Middle Pleistocene. Primers based on 54 pairs of orthologous SSR-containing sequences between the two Primula species were designed and verified. About half of these pairs successfully amplified for both species. Of the 959 single copy nuclear genes shared by four model plants (known as APVO genes), 111 single copy nuclear genes were verified as being present in both Primula species and exon-anchored and intron-spanned primers were designed for use.

We characterized the transcriptomes for the two Primula species, and produced an unprecedented amount of genomic resources for these important garden plants. Evolutionary analysis of these two Primula species not only revealed a more precise divergence time, but also provided some novel insights into how differential adaptations occurred in extreme habitats. Furthermore, we developed two sets of genetic markers, single copy nuclear genes and nuclear microsatellites (EST-SSR). Both these sets of markers will facilitate studies on the genetic improvement, population genetics and phylogenetics of this rapidly adapting taxon.

Sugarcane is an economically important crop contributing about 80% and 40% to the world sugar and ethanol production, respectively. The complicated genetics consequential to its complex polyploid genome, however, have impeded efforts to improve sugar yield and related important agronomic traits. Modern sugarcane cultivars are complex hybrids derived mainly from crosses among its progenitor species, S. officinarum and S. spontanuem, and to a lesser degree, S. robustom. Atypical of higher plants, sugarcane stores its photoassimilates as sucrose rather than as starch in its parenchymous stalk cells. In the sugar biosynthesis pathway, sucrose synthase (SuSy, UDP-glucose: D-fructose 2-a-D-glucosyltransferase, EC 2.4.1.13) is a key enzyme in the regulation of sucrose accumulation and partitioning by catalyzing the reversible conversion of sucrose and UDP into UDP-glucose and fructose. However, little is known about the sugarcane SuSy gene family members and hence no definitive studies have been reported regarding allelic diversity of SuSy gene families in Saccharum species.

We identified and characterized a total of five sucrose synthase genes in the three sugarcane progenitor species through gene annotation and PCR haplotype analysis by analyzing 70 to 119 PCR fragments amplified from intron-containing target region. We detected all but one (i.e. ScSuSy5) of ScSuSy transcripts in five tissue types of three Saccharum species. The average SNP frequency was one SNP per 108 bp, 81 bp, and 72 bp in S. officinarum, S. robustom, and S. spontanuem respectively. The average shared SNP is 15 between S. officinarum and S. robustom, 7 between S. officinarum and S. spontanuem , and 11 between S. robustom and S. spontanuem. We identified 27, 35, and 32 haplotypes from the five ScSuSy genes in S. officinarum, S. robustom, and S. spontanuem respectively. Also, 12, 11, and 9 protein sequences were translated from the haplotypes in S. officinarum, S. robustom, S. spontanuem, respectively. Phylogenetic analysis showed three separate clusters composed of SbSuSy1 and SbSuSy2, SbSuSy3 and SbSuSy5, and SbSuSy4.

The five members of the SuSy gene family evolved before the divergence of the genera in the tribe Andropogoneae at least 12 MYA. Each ScSuSy gene showed at least one non-synonymous substitution in SNP haplotypes. The SNP frequency is the lowest in S. officinarum, intermediate in S. robustum, and the highest in S. spontaneum, which may reflect the timing of the two rounds of whole genome duplication in these octoploids. The higher rate of shared SNP frequency between S. officinarum and S. robustum than between S. officinarum and in S. spontaneum confirmed that the speciation event separating S. officinarum and S. robustum occurred after their common ancestor diverged from S, spontaneum. The SNP and haplotype frequencies in three Saccharum species provide fundamental information for designing strategies to sequence these autopolyploid genomes.

Identification of QTL with large phenotypic effects conserved across genetic backgrounds and environments is one of the prerequisites for crop improvement using marker assisted selection (MAS). The objectives of this study were to identify meta-QTL (mQTL) for grain yield (GY) and anthesis silking interval (ASI) across 18 bi-parental maize populations evaluated in the same conditions across 2-4 managed water stressed and 3-4 well watered environments.

The meta-analyses identified 68 mQTL (9 QTL specific to ASI, 15 specific to GY, and 44 for both GY and ASI). Mean phenotypic variance explained by each mQTL varied from 1.2 to 13.1% and the overall average was 6.5%. Few QTL were detected under both environmental treatments and/or multiple (>4 populations) genetic backgrounds. The number and 95% genetic and physical confidence intervals of the mQTL were highly reduced compared to the QTL identified in the original studies. Each physical interval of the mQTL consisted of 5 to 926 candidate genes.

Meta-analyses reduced the number of QTL by 68% and narrowed the confidence intervals up to 12-fold. At least the 4 mQTL (mQTL2.2, mQTL6.1, mQTL7.5 and mQTL9.2) associated with GY under both water-stressed and well-watered environments and detected up to 6 populations may be considered for fine mapping and validation to confirm effects in different genetic backgrounds and pyramid them into new drought resistant breeding lines. This is the first extensive report on meta-analysis of data from over 3100 individuals genotyped using the same SNP platform and evaluated in the same conditions across a wide range of managed water-stressed and well-watered environments.

About 80% of today's land plants are able to establish an arbuscular mycorrhizal (AM) symbiosis with Glomeromycota fungi to improve their access to nutrients and water in the soil. On the molecular level, the development of AM symbioses is only partly understood, due to the asynchronous development of the microsymbionts in the host roots. Although many genes specifically activated during fungal colonization have been identified, genome-wide information on the exact place and time point of their activation remains limited.

In this study, we relied on a combination of laser-microdissection and the use of Medicago GeneChips to perform a genome-wide analysis of transcription patterns in defined cell-types of Medicago truncatula roots mycorrhized with Glomus intraradices. To cover major stages of AM development, we harvested cells at 5-6 and at 21 days post inoculation (dpi). Early developmental stages of the AM symbiosis were analysed by monitoring gene expression in appressorial and non-appressorial areas from roots harbouring infection units at 5-6 dpi. Here, the use of laser-microdissection for the first time enabled the targeted harvest of those sites, where fungal hyphae first penetrate the root. Circumventing contamination with developing arbuscules, we were able to specifically detect gene expression related to early infection events. To cover the late stages of AM formation, we studied arbusculated cells, cortical cells colonized by intraradical hyphae, and epidermal cells from mature mycorrhizal roots at 21 dpi. Taken together, the cell-specific expression patterns of 18014 genes were revealed, including 1392 genes whose transcription was influenced by mycorrhizal colonization at different stages, namely the pre-contact phase, the infection of roots via fungal appressoria, the subsequent colonization of the cortex by fungal hyphae, and finally the formation of arbuscules. Our cellular expression patterns identified distinct groups of AM-activated genes governing the sequential reprogramming of host roots towards an accommodation of microsymbionts, including 42 AM-activated transcription factor genes.

Our genome-wide analysis provides novel information on the cell-specific activity of AM-activated genes during both early and late stages of AM development, together revealing the road map of fine-tuned adjustments of transcript accumulation within root tissues during AM fungal colonization.

Auxin plays important roles in hormone crosstalk and the plant's stress response. The auxin-responsive Gretchen Hagen3 (GH3) gene family maintains hormonal homeostasis by conjugating excess indole-3-acetic acid (IAA), salicylic acid (SA), and jasmonic acids (JAs) to amino acids during hormone- and stress-related signaling pathways. With the sequencing of the apple (Malus x domestica) genome completed, it is possible to carry out genomic studies on GH3 genes to indentify candidates with roles in abiotic/biotic stress responses.

Malus sieversii Roem., an apple rootstock with strong drought tolerance and the ancestral species of cultivated apple species, was used as the experimental material. Following genome-wide computational and experimental identification of MdGH3 genes, we showed that MdGH3s were differentially expressed in the leaves and roots of M. sieversii and that some of these genes were significantly induced after various phytohormone and abiotic stress treatments. Given the role of GH3 in the negative feedback regulation of free IAA concentration, we examined whether phytohormones and abiotic stresses could alter the endogenous auxin level. By analyzing the GUS activity of DR5::GUS-transformed Arabidopsis seedlings, we showed that ABA, SA, salt, and cold treatments suppressed the auxin response. These findings suggest that other phytohormones and abiotic stress factors might alter endogenous auxin levels.

Previous studies showed that GH3 genes regulate hormonal homeostasis. Our study indicated that some GH3 genes were significantly induced in M. sieversii after various phytohormone and abiotic stress treatments, and that ABA, SA, salt, and cold treatments reduce the endogenous level of axuin. Taken together, this study provides evidence that GH3 genes play important roles in the crosstalk between auxin, other phytohormones, and the abiotic stress response by maintaining auxin homeostasis.

The antifungal therapy caspofungin is a semi-synthetic derivative of pneumocandin B0, a lipohexapeptide produced by the fungus Glarea lozoyensis, and was the first member of the echinocandin class approved for human therapy. The nonribosomal peptide synthetase (NRPS)-polyketide synthases (PKS) gene cluster responsible for pneumocandin biosynthesis from G. lozoyensis has not been elucidated to date. In this study, we report the elucidation of the pneumocandin biosynthetic gene cluster by whole genome sequencing of the G. lozoyensis wild-type strain ATCC 20868.

The pneumocandin biosynthetic gene cluster contains a NRPS (GLNRPS4) and a PKS (GLPKS4) arranged in tandem, two cytochrome P450 monooxygenases, seven other modifying enzymes, and genes for L-homotyrosine biosynthesis, a component of the peptide core. Thus, the pneumocandin biosynthetic gene cluster is significantly more autonomous and organized than that of the recently characterized echinocandin B gene cluster. Disruption mutants of GLNRPS4 and GLPKS4 no longer produced the pneumocandins (A0 and B0), and the Deltaglnrps4 and Deltaglpks4 mutants lost antifungal activity against the human pathogenic fungus Candida albicans. In addition to pneumocandins, the G. lozoyensis genome encodes a rich repertoire of natural product-encoding genes including 24 PKSs, six NRPSs, five PKS-NRPS hybrids, two dimethylallyl tryptophan synthases, and 14 terpene synthases.

Characterization of the gene cluster provides a blueprint for engineering new pneumocandin derivatives with improved pharmacological properties. Whole genome estimation of the secondary metabolite-encoding genes from G. lozoyensis provides yet another example of the huge potential for drug discovery from natural products from the fungal kingdom.

The sequenced genomes of cucumber, melon and watermelon have relatively few R-genes, with 70, 75 and 55 copies only, respectively. The mechanism for low copy number of R-genes in Cucurbitaceae genomes remains unknown.

Manual annotation of R-genes in the sequenced genomes of Cucurbitaceae species showed that approximately half of them are pseudogenes. Comparative analysis of R-genes showed frequent loss of R-gene loci in different Cucurbitaceae species. Phylogenetic analysis, data mining and PCR cloning using degenerate primers indicated that Cucurbitaceae has limited number of R-gene lineages (subfamilies). Comparison between R-genes from Cucurbitaceae and those from poplar and soybean suggested frequent loss of R-gene lineages in Cucurbitaceae. Furthermore, the average number of R-genes per lineage in Cucurbitaceae species is approximately 1/3 that in soybean or poplar. Therefore, both loss of lineages and deficient duplications in extant lineages accounted for the low copy number of R-genes in Cucurbitaceae. No extensive chimeras of R-genes were found in any of the sequenced Cucurbitaceae genomes. Nevertheless, one lineage of R-genes from Trichosanthes kirilowii, a wild Cucurbitaceae species, exhibits chimeric structures caused by gene conversions, and may contain a large number of distinct R-genes in natural populations.

Cucurbitaceae species have limited number of R-gene lineages and each genome harbors relatively few R-genes. The scarcity of R-genes in Cucurbitaceae species was due to frequent loss of R-gene lineages and infrequent duplications in extant lineages. The evolutionary mechanisms for large variation of copy number of R-genes in different plant species were discussed.

Here we present an overview of our existing knowledge on the function of RIN4 as a regulator of plant defense and as a guardee of multiple plant R-proteins. Domain analysis of RIN4 reveals two NOI domains. The NOI domain was originally identified in a screen for nitrate induced genes. The domain is comprised of approximately 30 amino acids and contains 2 conserved motifs (PXFGXW and Y/FTXXF). The NOI gene family contains members exclusively from the plant lineage as far back as moss. In addition to the conserved NOI domain, members within the family also contain conserved C-terminal cysteine residue(s) which are sites for acylation and membrane tethering. Other than these two characteristic features, the sequence of the family of NOI-containing proteins is diverse and, with the exception of RIN4, their functions are not known. Recently published interactome data showing interactions between RIN4 and components of the exocyst complex prompt us to raise the hypothesis that RIN4 might be involved in defense associated vesicle trafficking.

The rice blast fungus, Magnaporthe oryzae is a destructive pathogen of rice and other related crops, causing significant yield losses worldwide. Endogenous small RNAs (sRNAs), including small interfering RNAs (siRNAs) and microRNAs (miRNAs) are critical components of gene regulation in many eukaryotic organisms. Recently several new species of sRNAs have been identified in fungi. This fact along with the availability of genome sequence makes M. oryzae a compelling target for sRNA profiling. We have examined sRNA species and their biosynthetic genes in M. oryzae, and the degree to which these elements regulate fungal stress responses. To this end, we have characterized sRNAs under different physiological stress conditions, which had not yet been examined in this fungus.

The resulting libraries are composed of more than 37 million total genome matched reads mapping to intergenic regions, coding sequences, retrotransposons, inverted, tandem, and other repeated regions of the genome with more than half of the small RNAs arising from intergenic regions. The 24 nucleotide (nt) size class of sRNAs was predominant. A comparison to transcriptional data of M. oryzae undergoing the same physiological stresses indicates that sRNAs play a role in transcriptional regulation for a small subset of genes. Support for this idea comes from generation and characterization of mutants putatively involved in sRNAs biogenesis; our results indicate that the deletion of Dicer-like genes and an RNA-Dependent RNA Polymerase gene increases the transcriptional regulation of this subset of genes, including one involved in virulence.

Various physiological stressors and in planta conditions alter the small RNA profile of the rice blast fungus. Characterization of sRNA biosynthetic mutants helps to clarify the role of sRNAs in transcriptional control.

Leucine-rich repeat receptor-like kinases (LRR-RLKs) comprise the largest group within the receptor-like kinase (RLK) superfamily in plants. This gene family plays critical and diverse roles in plant growth, development and stress response. Although the LRR-RLK families in Arabidopsis and rice have been previously analyzed, no comprehensive studies have been performed on this gene family in tree species.

In this work, 379 LRR-RLK genes were retrieved from the Populus trichocarpa genome and further grouped into 14 subfamilies based on their structural and sequence similarities. Approximately 82% (312 out of 379) of the PtLRR-RLK genes are located in segmental duplication blocks indicating the role of duplication process in the expansion of this gene family. The conservation and variation in motif composition and intron/exon arrangement among PtLRR-RLK subfamilies were analyzed to provide additional support for their phylogenetic relationship and more importantly to indicate the potential divergence in their functions. Expression profiling of PtLRR-RLKs showed that they were differentially expressed in different organs and tissues and some PtLRR-RLKs were specifically expressed in meristem tissues, which indicated their potential involvement in tissue development and differentiation. For most AtLRR-RLKs with defined functions, Populus homologues exhibiting similar expression patterns could be identified, which might indicate the functional conservation during evolution. Among 12 types of environmental cues analyzed by the genome-wide microarray data, PtLRR-RLKs showed specific responses to shoot organogenesis, wounding, low ammonium feeding, hypoxia and seasonal dormancy, but not to drought, re-watering after drought, flooding, AlCl3 treatment and bacteria or fungi treatments.

This study provides the first comprehensive genomic analysis of the Populus LRR-RLK gene family. Segmental duplication contributes significantly to the expansion of this gene family. Populus and Arabidopsis LRR-RLK homologues not only share similar genetic structures but also exhibit comparable expression patterns which point to the possible functional conservation of these LRR-RLKs in two model systems. Transcriptome profiling provides the first insight into the functional divergence among PtLRR-RLK gene subfamilies and suggests that they might take important roles in growth and adaptation of tree species.