Plant Genomics

Ever since the recent completion of the peach genome, the focus of genetic research in this area has turned to the identification of genes related to important traits, such as fruit aroma volatiles. Of the over 100 volatile compounds described in peach, lactones most likely have the strongest effect on fruit aroma, while esters, terpenoids, and aldehydes have minor, yet significant effects. The identification of key genes underlying the production of aroma compounds is of interest for any fruit-quality improvement strategy.

Volatile (52 compounds) and gene expression (4348 genes) levels were profiled in peach fruit from a maturity time-course series belonging to two peach genotypes that showed considerable differences in maturation characteristics and postharvest ripening. This data set was analyzed by complementary correlation-based approaches to discover the genes related to the main aroma-contributing compounds: lactones, esters, and phenolic volatiles, among others. As a case study, one of the candidate genes was cloned and expressed in yeast to show specificity as an omega-6 Oleate desaturase, which may be involved in the production of a precursor of lactones/esters.

Our approach revealed a set of genes (an alcohol acyl transferase, fatty acid desaturases, transcription factors, protein kinases, cytochromes, etc.) that are highly associated with peach fruit volatiles, and which could prove useful in breeding or for biotechnological purposes.

The late blight pathogen Phytophthora infestans can attack both potato foliage and tubers. Although interaction transcriptome dynamics between potato foliage and various pathogens have been reported, no transcriptome study has focused specifically upon how potato tubers respond to pathogen infection. When inoculated with P. infestans, tubers of nontransformed 'Russet Burbank' (WT) potato develop late blight disease while those of transgenic 'Russet Burbank' line SP2211 (+RB), which expresses the potato late blight resistance gene RB (Rpi-blb1), do not. We compared transcriptome responses to P. infestans inoculation in tubers of these two lines.

We demonstrated the practicality of RNA-seq to study tetraploid potato and present the first RNA-seq study of potato tuber diseases. A total of 483 million paired end Illumina RNA-seq reads were generated, representing the transcription of around 30,000 potato genes. Differentially expressed genes, gene groups and ontology bins that exhibited differences between the WT and +RB lines were identified. P. infestans transcripts, including those of known effectors, were also identified.

Faster and stronger activation of defense related genes, gene groups and ontology bins correlate with successful tuber resistance against P. infestans. Our results suggest that the hypersensitive response is likely a general form of resistance against the hemibiotrophic P. infestans---even in potato tubers, organs that develop below ground.

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.

Several beta-galactosidases of the Glycosyl Hydrolase 35 (GH35) family have been characterized, and many of these modify cell wall components, including pectins, xyloglucans, and arabinogalactan proteins. The phloem fibres of flax (Linum usitatissimum) have gelatinous-type cell walls that are rich in crystalline cellulose and depend on beta-galactosidase activity for their normal development. In this study, we investigate the transcript expression patterns and inferred evolutionary relationships of the complete set of flax GH35 genes, to better understand the functions of these genes in flax and other species.

Using the recently published flax genome assembly, we identified 43 beta-galactosidase-like (BGAL) genes, based on the presence of a GH35 domain. Phylogenetic analyses of their protein sequences clustered them into eight sub-families. Sub-family B, whose members in other species were known to be expressed in developing flowers and pollen, was greatly under represented in flax (p-value < 0.01). Sub-family A5, whose sole member from arabidopsis has been described as its primary xyloglucan BGAL, was greatly expanded in flax (p-value < 0.01). A number of flax BGALs were also observed to contain non-consensus GH35 active sites. Expression patterns of the flax BGALs were investigated using qRT-PCR and publicly available microarray data. All predicted flax BGALs showed evidence of expression in at least one tissue.

Flax has a large number of BGAL genes, which display a distinct distribution among the BGAL sub-families, in comparison to other closely related species with available whole genome assemblies. Almost every flax BGAL was expressed in fibres, the majority of which expressed predominately in fibres as compared to other tissues, suggesting an important role for the expansion of this gene family in the development of this species as a fibre crop. Variations displayed in the canonical GH35 active site suggest a variety of roles unique to flax, which will require further characterization.

Interspecific hybridization creates individuals harboring diverged genomes. The interaction of these genomes can generate successful evolutionary novelty or disadvantageous genomic conflict. Annual sunflowers Helianthus annuus and H. petiolaris have a rich history of hybridization in natural populations. Although first-generation hybrids generally have low fertility, hybrid swarms that include later generation and fully fertile backcross plants have been identified, as well as at least three independently-originated stable hybrid taxa. We examine patterns of transcript accumulation in the earliest stages of hybridization of these species via analyses of transcriptome sequences from laboratory-derived F1 offspring of an inbred H. annuus cultivar and a wild H. petiolaris accession.

While nearly 14% of the reference transcriptome showed significant accumulation differences between parental accessions, total F1 transcript levels showed little evidence of dominance, as midparent transcript levels were highly predictive of transcript accumulation in F1 plants. Allelic bias in F1 transcript accumulation was detected in 20% of transcripts containing sufficient polymorphism to distinguish parental alleles; however the magnitude of these biases were generally smaller than differences among parental accessions.

While analyses of allelic bias suggest that cis regulatory differences between H. annuus and H. petiolaris are common, their effect on transcript levels may be more subtle than trans-acting regulatory differences. Overall, these analyses found little evidence of regulatory incompatibility or dominance interactions between parental genomes within F1 hybrid individuals, although it is unclear whether this is a legacy or an enabler of introgression between species.

Phytophthora infestans, the cause of potato late blight, is infamous for having triggered the Irish Great Famine in the 1840s. Until the late 1970s, P. infestans diversity outside of its Mexican center of origin was low, and one scenario held that a single strain, US-1, had dominated the global population for 150 years; this was later challenged based on DNA analysis of historical herbarium specimens. We have compared the genomes of 11 herbarium and 15 modern strains. We conclude that the nineteenth century epidemic was caused by a unique genotype, HERB-1, that persisted for over 50 years. HERB-1 is distinct from all examined modern strains, but it is a close relative of US-1, which replaced it outside of Mexico in the twentieth century. We propose that HERB-1 and US-1 emerged from a metapopulation that was established in the early 1800s outside of the species' center of diversity.

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.